18 max(
alpha.oldTime(), phase.residualAlpha())
28 PtrList<surfaceScalarField>
AFfs(
fluid.AFfs());
32 phaseModel& phase =
fluid.movingPhases()[movingPhasei];
39 IOobject::groupName(
"rAUf", phase.name()),
61 max(
alphafs[phase.index()], phase.residualAlpha())
81 phaseModel& phase =
fluid.movingPhases()[movingPhasei];
82 MRF.correctBoundaryFlux(phase.U(), phase.phiRef());
86 PtrList<surfaceScalarField> phigFs(
phases.size());
106 -
fluid.surfaceTension(phase)*
mesh.magSf()
122 phaseModel& phase =
fluid.movingPhases()[movingPhasei];
129 IOobject::groupName(
"phiHbyA", phase.name()),
136 - phigFs[phase.index()]
143 if (!partialElimination)
145 PtrList<surfaceScalarField> phiKdPhifs(
fluid.phiKdPhifs(
rAUfs));
149 if (phiKdPhifs.set(
phasei))
177 if (partialElimination)
179 PtrList<surfaceScalarField> phiKdPhifs(
fluid.phiKdPhifs(
rAUfs));
183 if (phiKdPhifs.set(
phasei))
214 surfaceScalarField::Boundary
phib(
phi.boundaryField());
220 alphafs[
phasei].boundaryField()*phase.phi()().boundaryField();
225 p_rgh.boundaryFieldRef(),
228 )/(
mesh.magSf().boundaryField()*
rAUf.boundaryField())
233 PtrList<fvScalarMatrix> pEqnComps(
phases.size());
234 PtrList<volScalarField> dmdts(
fluid.dmdts());
241 if (phase.compressible())
247 IOobject::groupName(
"phid", phase.name()),
275 pEqnComps[
phasei].faceFluxCorrectionPtr(
nullptr);
277 pEqnComps[
phasei].relax();
303 if (pEqnComps.set(
phasei))
319 if (pEqnComps.set(
phasei))
338 while (
pimple.correctNonOrthogonal())
352 if (pEqnComps.set(
phasei))
354 pEqn += pEqnComps[
phasei];
362 if (
pimple.finalNonOrthogonalIter())
370 phaseModel& phase =
fluid.movingPhases()[movingPhasei];
377 if (pEqnComps.set(phase.index()))
379 phase.divU(-pEqnComps[phase.index()] &
p_rgh);
383 if (partialElimination)
391 mSfGradp = pEqnIncomp.flux()/
rAUf;
395 phaseModel& phase =
fluid.movingPhases()[movingPhasei];
398 phase.URef().correctBoundaryConditions();
414 phase.thermoRef().rho() += phase.thermo().psi()*(
p_rgh -
p_rgh_0);
420 p_rgh.correctBoundaryConditions();
PtrList< surfaceScalarField > AFfs(fluid.AFfs())
tmp< fvMatrix< Type > > correction(const fvMatrix< Type > &)
Return the correction form of the given matrix by subtracting the matrix multiplied by the current fi...
fvMatrix< scalar > fvScalarMatrix
multiphaseSystem::phaseModelList & phases
Info<< "Creating field kinetic energy K\"<< endl;volScalarField K("K", 0.5 *magSqr(U));if(U.nOldTimes()){ volVectorField *Uold=&U.oldTime();volScalarField *Kold=&K.oldTime();*Kold==0.5 *magSqr(*Uold);while(Uold->nOldTimes()) { Uold=&Uold-> oldTime()
PtrList< surfaceScalarField > alphaRho0fs(phases.size())
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const dimensionedScalar & pMin
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
forAll(U.boundaryField(), patchi)
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
tmp< GeometricField< Type, faPatchField, areaMesh > > ddt(const dimensioned< Type > dt, const faMesh &mesh)
PtrList< surfaceScalarField > alpharAUfs(phases.size())
surfaceScalarField phid("phid", fvc::interpolate(psi) *(fvc::flux(HbyA)+MRF.zeroFilter(rhorAUf *fvc::ddtCorr(rho, U, phi)/fvc::interpolate(rho))))
const surfaceScalarField & ghf
GeometricField< scalar, fvPatchField, volMesh > volScalarField
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for expressions::valueTypeCode::INVALID.
tmp< GeometricField< Type, faPatchField, areaMesh > > laplacian(const GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
PtrList< surfaceScalarField > rAUfs
tmp< volScalarField > byDt(const volScalarField &vf)
const uniformDimensionedVectorField & g
const dimensionSet dimForce
const dimensionSet dimDensity
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU))
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
surfaceScalarField ghSnGradRho(ghf *fvc::snGrad(rho) *mesh.magSf())
PtrList< fvVectorMatrix > UEqns(fluid.phases().size())
PtrList< surfaceScalarField > alphafs(phases.size())
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
setSnGrad< fixedFluxPressureFvPatchScalarField >(p_rgh.boundaryFieldRef(),(phiHbyA.boundaryField() - MRF.relative(phib))/(mesh.magSf().boundaryField() *rAUf.boundaryField()))
const volScalarField & gh
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
volScalarField p_rgh_0(p_rgh)
PtrList< surfaceScalarField > phiHbyAs(fluid.phases().size())
PtrList< surfaceScalarField > phiFfs(fluid.phiFfs(rAUfs))
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
const dimensionSet dimArea(sqr(dimLength))
static constexpr const zero Zero
Global zero (0)
const dimensionSet dimVelocity