3 surfaceScalarField
rAUf(
"rAUf", fvc::interpolate(
rAU));
9 + fvc::interpolate(
rho*
rAU)*fvc::ddtCorr(
U,
phi)
30 const rhoThermo&
thermo = phase().thermo();
31 const volScalarField&
rho =
thermo.rho()();
48 while (
pimple.correctNonOrthogonal())
50 fvScalarMatrix p_rghEqnIncomp
56 tmp<fvScalarMatrix> p_rghEqnComp;
61 tmp<fvScalarMatrix> hmm
73 p_rghEqnComp.ref() += hmm;
86 if (
pimple.finalNonOrthogonalIter())
89 for (phaseModel& phase :
mixture.phases())
101 +
rAU*fvc::reconstruct((
phig + p_rghEqnIncomp.flux())/
rAUf);
102 U.correctBoundaryConditions();
118 Info<<
"max(U) " <<
max(mag(
U)).value() << endl;
119 Info<<
"min(p_rgh) " << min(
p_rgh).value() << endl;
CGAL::Exact_predicates_exact_constructions_kernel K
const surfaceScalarField & ghf
const volScalarField & gh
const dimensionedScalar & pMin
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
surfaceScalarField phig("phig", -rhorAUf *ghf *fvc::snGrad(rho) *mesh.magSf())
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU))
tmp< volScalarField > rAU
volScalarField p_rgh_0(p_rgh)
PtrList< fvScalarMatrix > p_rghEqnComps(mixture.phases().size())
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
#define forAllConstIters(container, iter)
Iterate across all elements of the container object with const access.