18 max(alpha.oldTime(), phase.residualAlpha())
29 PtrList<surfaceScalarField> AFfs(
fluid.AFfs());
33 phaseModel& phase =
fluid.movingPhases()[movingPhasei];
40 IOobject::groupName(
"rAUf", phase.name()),
62 max(
alphafs[phase.index()], phase.residualAlpha())
80 fluid.correctBoundaryFlux();
83 PtrList<surfaceScalarField> phigFs(
phases.size());
103 -
fluid.surfaceTension(phase)*
mesh.magSf()
116 PtrList<surfaceScalarField> phiHbyAs(
phases.size());
119 phaseModel& phase =
fluid.movingPhases()[movingPhasei];
137 - phigFs[phase.index()],
146 PtrList<surfaceScalarField> phiKdPhifs(
fluid.phiKdPhifs(
rAUfs));
150 if (phiKdPhifs.set(
phasei))
201 surfaceScalarField::Boundary phib
203 surfaceScalarField::Internal::null(),
212 alphafs[
phasei].boundaryField()*phase.phi()().boundaryField();
215 setSnGrad<fixedFluxPressureFvPatchScalarField>
217 p_rgh.boundaryFieldRef(),
220 )/(
mesh.magSf().boundaryField()*rAUf.boundaryField())
231 while (
pimple.correctNonOrthogonal())
249 if (
fluid.incompressible())
262 if (
pimple.finalNonOrthogonalIter())
270 phaseModel& phase =
fluid.movingPhases()[movingPhasei];
273 phiHbyAs[phase.index()]
280 if (partialElimination)
288 phaseModel& phase =
fluid.movingPhases()[movingPhasei];
290 MRF.makeRelative(phase.phiRef());
297 phaseModel& phase =
fluid.movingPhases()[movingPhasei];
304 phase.URef().correctBoundaryConditions();
316 if (p_rgh.needReference() &&
fluid.incompressible())
334 phase.thermoRef().rho() += phase.thermo().psi()*(p_rgh -
p_rgh_0);
340 p_rgh.correctBoundaryConditions();
pressureReference & pressureReference
fvMatrix< scalar > fvScalarMatrix
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
pimpleNoLoopControl & pimple
mixture MRF().makeRelative(phiHbyA)
PtrList< surfaceScalarField > alphafs(phases.size())
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
volScalarField alpha(IOobject("alpha", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
MRF makeRelative(fvc::interpolate(rho), phiHbyA)
tmp< surfaceScalarField > constrainPhiHbyA(const tmp< surfaceScalarField > &tphiHbyA, const volVectorField &U, const volScalarField &p)
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
scalar getRefCellValue(const volScalarField &field, const label refCelli)
Return the current value of field in the reference cell.
PtrList< fvScalarMatrix > pEqnComps(phases.size())
GeometricField< scalar, fvPatchField, volMesh > volScalarField
const dimensionSet dimTime
Foam::fvConstraints & fvConstraints
volScalarField p_rgh_0(p_rgh)
const dimensionSet dimFlux
const dimensionSet dimDensity
const dimensionSet dimForce
PtrList< surfaceScalarField > alpharAUfs(phases.size())
const surfaceScalarField & ghf
PtrList< surfaceScalarField > phiFfs(fluid.phiFfs(rAUfs))
const dimensionSet dimVelocity
PtrList< surfaceScalarField > rAUfs
word name(const complex &)
Return a string representation of a complex.
tmp< volScalarField > byDt(const volScalarField &vf)
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
PtrList< surfaceScalarField > alphaRho0fs(phases.size())
tmp< surfaceScalarField > absolute(const tmp< surfaceScalarField > &tphi, const volVectorField &U)
Return the given relative flux in absolute form.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
phaseSystem::phaseModelList & phases
const volScalarField & gh
dimensioned< scalar > mag(const dimensioned< Type > &)
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU))
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const dimensionedVector & g
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)