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 }
rhoPhi
rhoPhi
Definition: alphaEqn.H:6
Foam::fvc::flux
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
Foam::max
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:40
rho1
volScalarField & rho1
Definition: setRegionFluidFields.H:25
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:531
Foam::vtk::Tools::zeroField
vtkSmartPointer< vtkFloatArray > zeroField(const word &name, const label size)
Create named field initialized to zero.
Definition: foamVtkToolsTemplates.C:320
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:6
Foam::min
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:26
rho2
volScalarField & rho2
Definition: setRegionFluidFields.H:28
Foam::Info
messageStream Info
Information stream (stdout output on master, null elsewhere)
Foam::surfaceScalarField
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Definition: surfaceFieldsFwd.H:64
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:34
alpha1
const volScalarField & alpha1
Definition: setRegionFluidFields.H:6