On Mon, 13 Apr 2015, grandrabbit wrote: > I am solving Navier Sotkes problem in a rotating reference frame. > So I just wonder how to implement the angular periodic boundary condition. > Should I just derive a class from PeriodicBoundary, add a rotation_vector > and then overload get_corresponding_pos(...)?
Hmm... for angular periodic boundaries of *scalar* values, this would be sufficient. For angular periodic boundaries with *vector* values, in general it would not. You'd be enforcing continuity between u_x on one side and u_x on the other, and between u_y on one side and u_y on the other; but u_x and u_y *shouldn't* match up after rotation. The options here are either some surgery inside the compute_periodic_constraints functor (to allow PB classes to specify a linear transformation), manual creation of boundary constraint equations, or reformulating your problem in terms of u_r and u_theta. --- Roy ------------------------------------------------------------------------------ BPM Camp - Free Virtual Workshop May 6th at 10am PDT/1PM EDT Develop your own process in accordance with the BPMN 2 standard Learn Process modeling best practices with Bonita BPM through live exercises http://www.bonitasoft.com/be-part-of-it/events/bpm-camp-virtual- event?utm_ source=Sourceforge_BPM_Camp_5_6_15&utm_medium=email&utm_campaign=VA_SF _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
