About the code “alphaPhi1 = alpha1Eqn.flux() ” in the solver compressibleTwoPhaseEule
In “alphaEqn.H” file, after solving alpha equation, it defines:
alphaPhi1 = alpha1Eqn.flux();
alphaPhi2 = phi - alphaPhi1;
Do someone know how to calculate alpha1Eqn.flux()? Is it different with “fvc::interpolate(alpha1)*phi1”?
In the UEqns.H,
+ fvm::div(alphaPhi1, U1)
- fvm::Sp(fvc::ddt(alpha1) + fvc::div(alphaPhi1), U1)
Do the codes mean “d(alpha1*U1)/dt + div(alpha1*U1* U1)” ? Why does it add the third term?
Thanks for your time.
|All times are GMT -4. The time now is 15:40.|