alphaEqn.H
Go to the documentation of this file.
1 {
2  word alphaScheme("div(phi,alpha)");
3 
4  surfaceScalarField alphaPhi
5  (
6  phi.name() + alpha1.name(),
7  fvc::flux
8  (
9  phi,
10  alpha1,
11  alphaScheme
12  )
13  );
14 
15  MULES::explicitSolve
16  (
17  geometricOneField(),
18  alpha1,
19  phi,
20  alphaPhi,
21  oneField(),
22  zeroField()
23  );
24 
25  rhoPhi = alphaPhi*(rho1 - rho2) + phi*rho2;
26 
27  Info<< "Phase-1 volume fraction = "
28  << alpha1.weightedAverage(mesh.Vsc()).value()
29  << " Min(" << alpha1.name() << ") = " << min(alpha1).value()
30  << " Max(" << alpha1.name() << ") = " << max(alpha1).value()
31  << endl;
32 }
rho1
const volScalarField & rho1
Definition: createFieldRefs.H:5
Foam::max
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
alpha1
volScalarField & alpha1(mixture.alpha1())
rhoPhi
rhoPhi
Definition: alphaEqn.H:25
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:18
alphaPhi
surfaceScalarField alphaPhi(phi.name()+alpha1.name(), fvc::flux(phi, alpha1, alphaScheme))
phi
phi
Definition: correctPhi.H:3
Foam::min
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
alphaScheme
const word alphaScheme(mesh.divScheme(divAlphaName)[1].wordToken())
rho2
const volScalarField & rho2
Definition: createFieldRefs.H:6
Foam::Info
messageStream Info
Foam::fvc::flux
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
Definition: fvcFlux.C:32
Foam::surfaceScalarField
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Definition: surfaceFieldsFwd.H:57
Foam::MULES::explicitSolve
void explicitSolve(const RdeltaTType &rDeltaT, const RhoType &rho, volScalarField &psi, const surfaceScalarField &phiPsi, const SpType &Sp, const SuType &Su)
Definition: MULESTemplates.C:38