FreeFOAM The Cross-Platform CFD Toolkit
alphaEqn.H
Go to the documentation of this file.
1 {
2  fvScalarMatrix alpha1Eqn
3  (
4  fvm::ddt(alpha1)
5  + fvm::div(phi, alpha1)
7  (
8  Dab + alphatab*turbulence->nut(), alpha1,
9  "laplacian(Dab,alpha1)"
10  )
11  );
12 
13  alpha1Eqn.solve();
14 
15  rhoPhi = alpha1Eqn.flux()*(rho1 - rho2) + phi*rho2;
16  rho = alpha1*rho1 + (scalar(1) - alpha1)*rho2;
17 
18  Info<< "Phase 1 volume fraction = "
19  << alpha1.weightedAverage(mesh.V()).value()
20  << " Min(alpha1) = " << min(alpha1).value()
21  << " Max(alpha1) = " << max(alpha1).value()
22  << endl;
23 }