alphaEqn.H
Go to the documentation of this file.
1{
2 // Temporarily making U relative to mesh motion
3 if (mesh.moving())
4 {
5 U -= fvc::reconstruct(mesh.phi());
6 }
7
8 // Updating alpha1
9 #include "alphaSuSp.H"
10 advector.advect(Sp, Su);
11
12 // Making U absolute again after advection step
13 if (mesh.moving())
14 {
15 U += fvc::reconstruct(mesh.phi());
16 }
17
18 #include "rhofs.H"
19 rhoPhi = advector.getRhoPhi(rho1f, rho2f);
20
21 alpha2 = 1.0 - alpha1;
22 mixture.correct();
23}
24
25scalar domainFraction = 0;
27{
28 const volScalarField& porosity = tporosity.cref();
29 rhoPhi *= scalar(1)/fvc::interpolate(porosity);
30 domainFraction = alpha1.weightedAverage(mesh.Vsc()*porosity).value();
31}
32else
33{
34 domainFraction = alpha1.weightedAverage(mesh.Vsc()).value();
35}
36
37Info<< "Phase-1 volume fraction = "
39 << " Min(" << alpha1.name() << ") = " << min(alpha1).value()
40 << " Max(" << alpha1.name() << ") = " << max(alpha1).value()
41 << endl;
Y[inertIndex] max(0.0)
const volScalarField & alpha1
U
Definition: pEqn.H:72
rhoPhi
Definition: alphaEqn.H:6
alpha2
Definition: alphaEqn.H:9
surfaceScalarField rho2f(fvc::interpolate(rho2))
surfaceScalarField rho1f(fvc::interpolate(rho1))
dynamicFvMesh & mesh
const bool porosityEnabled(porosityProperties.getOrDefault< bool >("porosityEnabled", false))
tmp< volScalarField > tporosity
zeroField Su
Definition: alphaSuSp.H:1
zeroField Sp
Definition: alphaSuSp.H:2
scalar domainFraction
Definition: alphaEqn.H:25
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:372
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:33
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:39
isoAdvection advector(alpha1, phi, U)