Projecting Solutions for Axisymmetric Flow

1

Hello,

Is there a way to project a solution from half of a rectangular domain into the other half of the domain? To reduce the computational cost of solving for the full domain of Axisymmetric Navier stokes equations in 2D, I’d like to only solve for half of the domain and project the solutions to the other half since it is axisymmetric. Please let me know if this is possible in FENICS.

image
image1215×247 27.4 KB

2

Here is an example of axisymmetric computations in solid mechanics

1 Like
3

Thanks, Bleyerj. Yes, I have seen this before. One thing that I have questions about is how to apply the boundary conditions. For example, if I only use the top half of the domain, what sort of boundary condition do I give the axis at which the symmetry occurs (if at all)?

image
image1570×362 59.2 KB

4

You can apply a slip boundary condition u_y=0 on the symmetry axis

5

Thanks again. Forgive me if this is a simple implementation but I don’t have experience implementing slip bcs. For the previous full domain, I had a no-slip BCs in my code using the following:

u_ZERO = Constant((0.0,0.0))
bc_noslip = DirichletBC(W.sub(0), u_ZERO, domainBoundaries, ID_BOTTOM)

How can I make it a slip condition as you mentioned? I can set u_y = 0 but what about u_x?