From: David K. <dav...@ak...> - 2018-03-01 05:46:38
|
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? Thanks, David |