|June 29, 2009, 13:33||
gammaEqn.H in the interPhaseChangeFoam solver
Join Date: Apr 2009
Posts: 171Rep Power: 8
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?
Last edited by isabel; June 30, 2009 at 06:48.
|July 7, 2009, 05:01||
Join Date: Mar 2009
Posts: 207Rep Power: 9
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?
|Thread||Thread Starter||Forum||Replies||Last Post|
|Working directory via command line||Luiz||CFX||4||March 6, 2011 21:02|
|why the solver reject it? Anyone with experience?||bearcat||CFX||6||April 28, 2008 14:08|
|compressible two phase flow in CFX4.4||youngan||CFX||0||July 1, 2003 23:32|
|CFX 5.5||Roued||CFX||1||October 2, 2001 16:49|
|Setting a B.C using UserFortran in 4.3||tokai||CFX||10||July 17, 2001 16:25|