I wanted to know how I can impose homogeneous Neumann boundary conditions for Lagrangre multiplier, using the HDG method for the following problem:
@U/@t = Laplace(U); in Omegax(0,T]
U(x,y,0) = U_0; in Omega
U.n = 0 ; on @Omega