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 }
phi
surfaceScalarField & phi
Definition: setRegionFluidFields.H:28
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:256
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:18
rhoPhi
rhoPhi
Definition: alphaEqn.H:115
alpha1
const volScalarField & alpha1
Definition: createFieldRefs.H:4
Foam::min
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
Foam::Info
messageStream Info
rho2
const dimensionedScalar & rho2
Definition: createFields.H:44
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:52
rho1
const dimensionedScalar & rho1
Definition: createFields.H:43
alphaPhi
surfaceScalarField alphaPhi(phi.name()+alpha1.name(), fvc::flux(phi, alpha1, alphaScheme))
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