[email protected] wrote: > Is it stil right if value of nodes on two different boundary sufaces/lines > are should same value as a constraint for periodic boundary condition in > the domain? > > > I need to have a > > node number 1 10 > grid x--------------------------------x > > subject to u(1) = u(10) or possibly u(1) = u(10) - g (for inhomogenousity) >
If you use DG elements, there is no degree of freedom on the boundary. I would see the boundary face in such a case like an interior face: you have to enforce the equality of u(1) and u(10) wekly by a flux term. The difficulty then reduces to building the correct FEFaceValues objects for both ends. G _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
