|June 20, 2008, 03:21||
Hello I have a doubt in the
Join Date: Mar 2009
Location: Champaign, USA
Posts: 90Rep Power: 8
I have a doubt in the PISO loop.
In sonicFoam.C, "phid" is defined by
surfaceScalarField phid =
(fvc::interpolate(rho*U) & mesh.Sf())
+ fvc::ddtPhiCorr(rUA, rho, U, phi)
later in solving for pressure, the code says
+ fvm::div(phid, p, "div(phid,p)")
- fvm::laplacian(rho*rUA, p)
Why are we dividing by p, when defining phid and later multiplying by p, when doing fvm::div(phid, p). The corresponding term in icoFoam is treated expicitly as fvc::div(phi).
So my second question is why is this term being treated implicitly in sonicFoam, but explicitly in icoFoam?
|Thread||Thread Starter||Forum||Replies||Last Post|
|SonicFoam||stvtiegh||OpenFOAM Running, Solving & CFD||1||June 26, 2008 09:30|
|New solver from sonicFoam||dimi||OpenFOAM Running, Solving & CFD||2||December 5, 2007 19:08|
|SonicFoam divphiU||dimi||OpenFOAM Running, Solving & CFD||3||June 25, 2007 05:47|
|SonicFoam question||shuo||OpenFOAM Running, Solving & CFD||0||October 15, 2006 22:53|
|Help me with sonicFoam||marcelo||OpenFOAM Running, Solving & CFD||6||December 10, 2005 03:57|