[alundilong]

Hi All:

Recently, I have gone through the interPhaseChangeFoam coming with OpenFOAM-2.3.0.

I think the term associated with divU in explicitSolve() is redundant. And My derivation is associated.

99 MULES::explicitSolve
100 (
101 geometricOneField(),
102 alpha1,
103 phi,
104 tphiAlphaCorr(),
105 vDotvmcAlphal,
106 (divU*alpha1 + vDotcAlphal)(), // divU*alpha1 is redundant, right?
107 1,
108 0
109 );
110
111 tphiAlpha = tphiAlphaCorr;
112 }

Please correct me if I was wrong.
Thanks!

[rapierrz]

Hi yijin,

Do you find reason of putting alpha1*divU in Su?

Why alpha1*divU there is in explicit source term?

is this source term always positive?

Thank you.