[qi.yang@polimi.it]

In twophasesystem.C, the following codes are in order to solve the phase fraction, right? but why there is no convective term (alpha1*phi1)?

if (pPrimeByA_.valid())
{
fvScalarMatrix alpha1Eqn
(
fvm::ddt(alpha1) - fvc::ddt(alpha1)
- fvm::laplacian(alpha1alpha2f()*pPrimeByA_(), alpha1, "bounded")
);

alpha1Eqn.relax();
alpha1Eqn.solve();

phase1_.alphaPhi() += alpha1Eqn.flux();
}
Hi everyone,

I have a question that I need to add one term in phase fraction equation (alpha1Eqn), The attachment is my continuity equation.

Thus, could I add the term like "fvc::laplacian(phase2.turbulence().nut/mysigma, alpha1)? Thanks for your attention.
Another thing is how to use the new twophasesystem.C after I compiled it. Just include the new file into twophaseeulerfoam.C? Thanks for your attention.

[qi.yang@polimi.it]

Nobody would like to give some suggestions?