CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   About the code “alphaPhi1 = alpha1Eqn.flux() ” in the solver compressibleTwoPhaseEule (http://www.cfd-online.com/Forums/openfoam-programming-development/114817-about-code-alphaphi1-alpha1eqn-flux-solver-compressibletwophaseeule.html)

aegean333 March 18, 2013 10:37

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::ddt(alpha1, U1)
+ 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 07:31.