From: John P. <jwp...@gm...> - 2018-03-01 16:43:16
|
On Wed, Feb 28, 2018 at 8:11 PM, David Knezevic <dav...@ak...> wrote: > I have only seen periodic boundary condition examples in libMesh where > boundary B is a translation of boundary A. This is what is implemented in > PeriodicBoundary::get_corresponding_pos(), where we have "return pt + > translation_vector;" > > I would like to implement a periodic boundary condition on a model with > circular symmetry, e.g. solve on sector of a disk with periodicity. To > implement this it seems like all I'd need to do is subclass > PeriodicBoundary and override get_corresponding_pos() to impose the > appropriate rotation rather than just a translation, and then add the > PeriodicBoundary subclass object to the system in the usual way. Is that > indeed all that's required, or would we need something more? Has anyone > tried this case before? > We have support for user-defined forward and inverse periodic boundary transform functions in MOOSE if you'd like to check and see how it's done there. -- John |