alphaEqn.H
Go to the documentation of this file.
1 // Update alpha1
2 #include "alphaSuSp.H"
3 advector.advect(Sp,(Su + divU*min(alpha1(), scalar(1)))());
4 
5 // Update rhoPhi
6 rhoPhi = advector.getRhoPhi(rho1, rho2);
7 alphaPhi10 = advector.alphaPhi();
8 
9 alpha2 = 1.0 - alpha1;
10 
11 Info<< "Phase-1 volume fraction = "
12  << alpha1.weightedAverage(mesh.Vsc()).value()
13  << " Min(" << alpha1.name() << ") = " << min(alpha1).value()
14  << " Max(" << alpha1.name() << ") - 1 = " << max(alpha1).value() - 1
15  << endl;
zeroField Su
Definition: alphaSuSp.H:1
alpha2
Definition: alphaEqn.H:9
rhoPhi
Definition: alphaEqn.H:6
isoAdvection advector(alpha1, phi, U)
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:40
volScalarField & rho1
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:531
dynamicFvMesh & mesh
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:26
volScalarField & rho2
alphaPhi10
Definition: alphaEqn.H:7
messageStream Info
Information stream (stdout output on master, null elsewhere)
zeroField divU
Definition: alphaSuSp.H:3
zeroField Sp
Definition: alphaSuSp.H:2
const volScalarField & alpha1