CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   MULES in interPhaseChangeFoam (http://www.cfd-online.com/Forums/openfoam/66470-mules-interphasechangefoam.html)

isabel July 16, 2009 07:06

MULES in interPhaseChangeFoam
 
Hello everybody

In gammaEqn.H of interPhaseChangeFoam solver we have these lines:


volScalarField Sp
(
IOobject
(
"Sp",
runTime.timeName(),
mesh
),
vDotvAlphal - vDotcAlphal
);
volScalarField Su
(
IOobject
(
"Su",
runTime.timeName(),
mesh
),
divU*gamma
+ vDotcAlphal
);


MULES::implicitSolve(oneField(), gamma, phi, phiGamma, Sp, Su, 1, 0);


why isn't Sp zero, if the equation to solve (MULESTemplates.C) is this:

d(1*psi)/dt + div(phigamma) + Sp*gamma - Su = 0

Is Sp different from zero to improve convergence? Or what?


All times are GMT -4. The time now is 21:20.