Yes, you are right. I was done successfully to make a periodic boundaries for DG system in deal.II and does not need to have the constraints for discontiunous system but should do compute flux correctly between two neighboring faces of corresponding periodic cells. Thank you anyway.
On Mon, August 3, 2009 8:50 pm, Guido Kanschat wrote: > [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
