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
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
Definition: fvMeshStitchersMoving.C:102
alphaScheme
const word alphaScheme(mesh.schemes().div(divAlphaName)[1].wordToken())
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
alpha1
volScalarField & alpha1(mixture.alpha1())
mesh
fvMesh & mesh
Definition: createMeshesPostProcess.H:9
rhoPhi
rhoPhi
Definition: alphaEqn.H:25
alphaPhi
surfaceScalarField alphaPhi(phi.name()+alpha1.name(), fvc::flux(phi, alpha1, alphaScheme))
Foam::min
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
Definition: fvMeshStitchersMoving.C:108
phi
phi
Definition: correctPhi.H:3
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:58
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:37