CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Main CFD Forum (http://www.cfd-online.com/Forums/main/)
-   -   implementation boundary condition for non-orthogonal mesh (http://www.cfd-online.com/Forums/main/75347-implementation-boundary-condition-non-orthogonal-mesh.html)

 m.iran April 22, 2010 10:40

implementation boundary condition for non-orthogonal mesh

hi everybody
I am writing a code that solves diffusion for non-orthogonal mesh and I don't know how could I insert boundary condition
I mean for example for 2D orthogonal system we have:

ap*phip=de*phie+dw*phiw+dn*phin+ds*phis+su

but in non-orthogonal system some terms will be added to source term, that is we will have:

ap*phip = deksi*phie+dwksi*phiw+dneta*phin+dseta*phis+deeta( phine-phise)-dweta(phinw-phisw)+dnksi(phine-phinw)-dsksi(phise-phisw)+su

and it is simiplified to :

ap*phip=deksi*phie+dwksi*phiw+dneta*phin+dseta*phi s+su

in which ;
su=su+deeta(phine-phise)-dweta(phinw-phisw)+dnksi(phine-phinw)-dsksi(phise-phisw)

NOW I DON'T KNOW HOW SHOULD I INSERT BOUNDARY CONDITION TYPE I AND TYPE II FOR THIS SYSTEM:confused:

THANKS EVERYBODY

 arjun April 22, 2010 16:54

look into cfd wiki and how diffusion term is descretized. Work the same way.

 All times are GMT -4. The time now is 15:28.