how to apply a symmetry boundary?

houkensjtu · December 10, 2011, 07:49

hey guys
i am working on coding a simulation by SIMPLE method, the flow is laminar and axissymmetric, all the momentum and pressure correction equations are written in cylindrical form.
however i found it conquers error especially near the axis.the velocity derivative du/dr is not zero.

here is how i applied the symmetry boundary condition;
for all u-velocity(parallel to the symmetry axis) cells near the axis, the discrete momentum equation is:
ap*u(i,j) = ae*u(i,j+1) + aw*u(i,j-1) + an*u(i+1,j) + as*u(i-1,j) +dp/dz*r*dz*dr
then i set u(i-1,j) = u(i,j) and repeat

i found that the problem maybe caused by that the "as" coefficient is very close to zero near the axis because radius is near to zero.
any ideas about it?

or just i have any thing other wrong?

plz tell me your advice,thx

Rami · December 12, 2011, 05:41

The FVM (as well as the FEM) does not need a BC at r=0 for an axisymmetric problems as it is a natural BC. You may notice that the face area there is 0, so the relevant coefficient would be 0. Using u(i-1)=u(i) [I assume i-1 is a ghost cell] will violate the true BC. Actually you do not need ghost cells either in the FVM.

December 10, 2011, 07:49	how to apply a symmetry boundary?	#1
houkensjtu Member HouKen Join Date: Jul 2011 Posts: 67 Rep Power: 14	hey guys i am working on coding a simulation by SIMPLE method, the flow is laminar and axissymmetric, all the momentum and pressure correction equations are written in cylindrical form. however i found it conquers error especially near the axis.the velocity derivative du/dr is not zero. here is how i applied the symmetry boundary condition; for all u-velocity(parallel to the symmetry axis) cells near the axis, the discrete momentum equation is: apu(i,j) = aeu(i,j+1) + awu(i,j-1) + anu(i+1,j) + asu(i-1,j) +dp/dzrdzdr then i set u(i-1,j) = u(i,j) and repeat i found that the problem maybe caused by that the "as" coefficient is very close to zero near the axis because radius is near to zero. any ideas about it? or just i have any thing other wrong? plz tell me your advice,thx

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Wind turbine simulation	Saturn	CFX	58	July 3, 2020 01:13
CFX13 Post Periodic interface	EtaEta	CFX	7	December 8, 2011 17:15
URGENT!!! F_UDSI at symmetry boundary	gemini	Fluent UDF and Scheme Programming	4	March 11, 2010 01:23
The behavior of droplets at the symmetry boundary	SJLee	Siemens	0	April 17, 2007 01:50
Dynamic Mesh & Symmetry boundary condition	Matt3o	FLUENT	0	January 27, 2006 10:07

December 12, 2011, 05:41		#2
Rami Senior Member Rami Ben-Zvi Join Date: Mar 2009 Posts: 155 Rep Power: 17	The FVM (as well as the FEM) does not need a BC at r=0 for an axisymmetric problems as it is a natural BC. You may notice that the face area there is 0, so the relevant coefficient would be 0. Using u(i-1)=u(i) [I assume i-1 is a ghost cell] will violate the true BC. Actually you do not need ghost cells either in the FVM.