|December 10, 2011, 08:49||
how to apply a symmetry boundary?
Join Date: Jul 2011
Posts: 66Rep Power: 6
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
|December 12, 2011, 06:41||
Join Date: Mar 2009
Posts: 147Rep Power: 8
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.
|Thread||Thread Starter||Forum||Replies||Last Post|
|Wind turbine simulation||Saturn||CFX||34||October 16, 2014 05:27|
|CFX13 Post Periodic interface||EtaEta||CFX||7||December 8, 2011 18:15|
|URGENT!!! F_UDSI at symmetry boundary||gemini||Fluent UDF and Scheme Programming||4||March 11, 2010 02:23|
|The behavior of droplets at the symmetry boundary||SJLee||CD-adapco||0||April 17, 2007 01:50|
|Dynamic Mesh & Symmetry boundary condition||Matt3o||FLUENT||0||January 27, 2006 11:07|