gammaEqn.H in the interPhaseChangeFoam solver
In gammaEqn.H in the interPhaseChangeFoam we can see this line:
MULES::implicitSolve(oneField(), gamma, phi, phiGamma, Sp, Su, 1, 0);
I know that:
gamma: is the actual value of gamma to be solved
phi: is the normal convective flux
phiGamma: U*gamma + gamma*(1-gamma)*U
1, 0 : max and min gamma values
But I don't unserstand what Sp and Su are. I think that Sp is a source term to solve the equation:
d(gamma)/dt + div(phigamma) = Sp
Am I all right?
What is Su term?
Hi, I think, your expression about the gamma equation is right.
In fact, to this equation,
Source term = vDotvAlphal * gamma + (1 - gamma) * vDotcAlphal = (vDotvAlphal - vDotcAlphal) * gamma + vDotcAlphal.
So, Sp = (vDotvAlphal - vDotcAlphal) * gamma , and it will be solved implicitly.
Su = vDotcAlphal, and it is an explicit term in this equation.
But, why Su also includes the term divU*gamma in gammaEqu.H ? If the MULES::implicitSolver was chose, this term should be deleted, right?
Because in MULESTemplates.C, phiPsi = phiGamma = phi*gamma + gamma*(1-gamma)*phi, the term div(phi*gamma) [=divU*gamma ??] has been actually included as follows:
+ fv::gaussConvectionScheme<scalar>(mesh, phi, UDs).fvmDiv(phi, psi)
- fvm::Sp(Sp, psi)
surfaceScalarField phiBD = psiConvectionDiffusion.flux();
surfaceScalarField& phiCorr = phiPsi;
phiCorr -= phiBD;
psiConvectionDiffusion + fvc::div(lambda*phiCorr),
Am I right?
|All times are GMT -4. The time now is 12:27.|