songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

compressibleTwoPhaseEulerFoam: updated PIMPLE loop

Browse Source
This commit is contained in:
Henry 2011-09-05 19:08:31 +01:00
parent 209aad1bdc
commit d47625903e
2 changed files with 1 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H
View File

@ -84,23 +84,12 @@
+ (alpha2/rho2)*pEqnComp2() + (alpha2/rho2)*pEqnComp2()
) + ) +
pEqnIncomp, pEqnIncomp,
mesh.solver mesh.solver(p.select(pimple.finalInnerIter(corr, nonOrth)))
(
p.select(pimple.finalInnerIter(corr, nonOrth))
/*
(
oCorr == pimple.nOuterCorr()-1
&& corr == nCorr-1
&& nonOrth == nNonOrthCorr
)
*/
)
); );
if (nonOrth == pimple.nNonOrthCorr()) if (nonOrth == pimple.nNonOrthCorr())
{ {
surfaceScalarField mSfGradp = pEqnIncomp.flux()/Dp; surfaceScalarField mSfGradp = pEqnIncomp.flux()/Dp;
phi1 += rAlphaAU1f*mSfGradp/fvc::interpolate(rho1); phi1 += rAlphaAU1f*mSfGradp/fvc::interpolate(rho1);
phi2 += rAlphaAU2f*mSfGradp/fvc::interpolate(rho2); phi2 += rAlphaAU2f*mSfGradp/fvc::interpolate(rho2);
phi = alpha1f*phi1 + alpha2f*phi2; phi = alpha1f*phi1 + alpha2f*phi2;

1
applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/readTwoPhaseEulerFoamControls.H
View File

@ -1,4 +1,3 @@
#include "readTimeControls.H" #include "readTimeControls.H"
int nAlphaCorr(readInt(pimple.dict().lookup("nAlphaCorr"))); int nAlphaCorr(readInt(pimple.dict().lookup("nAlphaCorr")));
Switch correctAlpha(pimple.dict().lookup("correctAlpha"));