FreeFOAM The Cross-Platform CFD Toolkit
alphaEqn.H
Go to the documentation of this file.
1 {
2  surfaceScalarField phiAlpha
3  (
4  IOobject
5  (
6  "phiAlpha",
7  runTime.timeName(),
8  mesh
9  ),
10  phi + rhoc*(mesh.Sf() & fvc::interpolate(Vdj))
11  );
12 
13  solve
14  (
15  fvm::ddt(rho, Alpha)
16  + fvm::div(phiAlpha, Alpha)
17  - fvm::laplacian(mut, Alpha)
18  );
19 
20  Info<< "Solid phase fraction = "
21  << Alpha.weightedAverage(mesh.V()).value()
22  << " Min(Alpha) = " << min(Alpha).value()
23  << " Max(Alpha) = " << max(Alpha).value()
24  << endl;
25 
26  Alpha.min(1.0);
27  Alpha.max(0.0);
28 
29  rho == rhoc/(scalar(1) + (rhoc/rhod - 1.0)*Alpha);
30  alpha == rho*Alpha/rhod;
31 }
32 
33 // ************************ vim: set sw=4 sts=4 et: ************************ //