> I have a partial differential equation in spherical coordinates of the > complex function Psi(theta,phi). So one will need a mesh which is spherical > manifold. Can I use dealii to design the periodic boundary conditions in > the phi angle such that my complex function Psi is periodic in 4*Pi not > 2*Pi ?
That's a bit tricky if you really want to solve on a spherical manifold because you have to cut it at phi=0 to implement your periodicity. Can you solve it on the rectangle [-pi,pi]x[0,4pi] instead? There are two examples that go in the directions you care about: step-45 deals with periodic boundary conditions, step-38 with solving on manifolds. > Another question: Can I have a new operator in the differential equation > which reverses the sign of one coordinate? For example: D Psi(theta , phi)= > Psi(theta , -phi) You mean for any point in the interior of your domain? That's a nonlocal operator that makes everything much more complicated... W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
