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

PISO in sonicFoam

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

LinkBack Thread Tools Search this Thread Display Modes
Old   June 20, 2008, 03:21
Default Hello I have a doubt in the
Join Date: Mar 2009
Location: Champaign, USA
Posts: 91
Rep Power: 14
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 Search this Thread
Search this Thread:

Advanced Search
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 Off
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 18: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 02:57

All times are GMT -4. The time now is 23:44.