Go to the documentation of this file.00001 {
00002 fvScalarMatrix alpha1Eqn
00003 (
00004 fvm::ddt(alpha1)
00005 + fvm::div(phi, alpha1)
00006 - fvm::laplacian
00007 (
00008 Dab + alphatab*turbulence->nut(), alpha1,
00009 "laplacian(Dab,alpha1)"
00010 )
00011 );
00012
00013 alpha1Eqn.solve();
00014
00015 rhoPhi = alpha1Eqn.flux()*(rho1 - rho2) + phi*rho2;
00016 rho = alpha1*rho1 + (scalar(1) - alpha1)*rho2;
00017
00018 Info<< "Phase 1 volume fraction = "
00019 << alpha1.weightedAverage(mesh.V()).value()
00020 << " Min(alpha1) = " << min(alpha1).value()
00021 << " Max(alpha1) = " << max(alpha1).value()
00022 << endl;
00023 }
