CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

PISO in sonicFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   June 20, 2008, 03:21
Default Hello I have a doubt in the
Join Date: Mar 2009
Location: Champaign, USA
Posts: 90
Rep Power: 10
srinath is on a distinguished road

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

fvScalarMatrix pEqn
fvm::ddt(psi, p)
+ 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?

srinath is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
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

All times are GMT -4. The time now is 12:39.