24 lines
582 B
C
24 lines
582 B
C
{
|
|
word alphaScheme("div(phi,alpha)");
|
|
|
|
surfaceScalarField alphaPhi
|
|
(
|
|
phi.name() + alpha1.name(),
|
|
fvc::flux
|
|
(
|
|
phi,
|
|
alpha1,
|
|
alphaScheme
|
|
)
|
|
);
|
|
|
|
MULES::explicitSolve(alpha1, phi, alphaPhi, 1, 0);
|
|
|
|
rhoPhi = alphaPhi*(rho1 - rho2) + phi*rho2;
|
|
|
|
Info<< "Phase-1 volume fraction = "
|
|
<< alpha1.weightedAverage(mesh.Vsc()).value()
|
|
<< " Min(" << alpha1.name() << ") = " << min(alpha1).value()
|
|
<< " Max(" << alpha1.name() << ") = " << max(alpha1).value()
|
|
<< endl;
|
|
}
|