27 PtrList<surfaceScalarField>
alphafs(
fluid.phases().size());
28 PtrList<volVectorField>
HbyAs(
fluid.phases().size());
30 PtrList<volScalarField>
rAUs(
fluid.phases().size());
34 for (phaseModel& phase :
fluid.phases())
36 MRF.makeAbsolute(phase.phi().oldTime());
37 MRF.makeAbsolute(phase.phi());
63 for (phaseModel& phase :
fluid.phases())
104 MRF.makeRelative(phase.phi().oldTime());
105 MRF.makeRelative(phase.phi());
115 auto dmIter =
fluid.dragModels().cbegin();
121 dmIter.good() && dcIter.good();
125 const phaseModel *phase2Ptr =
nullptr;
127 if (&phase == &dmIter()->phase1())
129 phase2Ptr = &dmIter()->phase2();
131 else if (&phase == &dmIter()->
phase2())
133 phase2Ptr = &dmIter()->phase1();
179 for (
const phaseModel& phase :
fluid.phases())
188 surfaceScalarField::Boundary
phib(
phi.boundaryField());
191 for (
const phaseModel& phase :
fluid.phases())
195 *(
mesh.Sf().boundaryField() & phase.U().boundaryField());
202 p_rgh.boundaryFieldRef(),
205 )/(
mesh.magSf().boundaryField()*
rAUf.boundaryField())
229 if (
pimple.finalNonOrthogonalIter())
236 for (phaseModel& phase :
fluid.phases())
245 *mSfGradp/phase.rho();
261 mSfGradp = pEqnIncomp.flux()/
rAUf;
266 for (phaseModel& phase :
fluid.phases())
284 phase.U().correctBoundaryConditions();
295 #include "continuityErrs.H"
fvMatrix< scalar > fvScalarMatrix
const word zeroGradientType
A zeroGradient patch field type.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
tmp< GeometricField< typename flux< Type >::type, fvsPatchField, surfaceMesh > > ddtCorr(const GeometricField< Type, fvPatchField, volMesh > &U, const GeometricField< Type, fvsPatchField, surfaceMesh > &Uf)
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
dimensioned< vector > dimensionedVector
Dimensioned vector obtained from generic dimensioned type.
autoPtr< multiphaseSystem::dragCoeffFields > dragCoeffs(fluid.dragCoeffs())
GeometricField< vector, fvPatchField, volMesh > volVectorField
const surfaceScalarField & ghf
GeometricField< scalar, fvPatchField, volMesh > volScalarField
tmp< GeometricField< Type, faPatchField, areaMesh > > laplacian(const GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
const uniformDimensionedVectorField & g
PtrList< volVectorField > HbyAs(fluid.phases().size())
const word & name() const noexcept
Return const reference to name.
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU))
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
PtrList< volScalarField > rAUs
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
PtrList< surfaceScalarField > phiHbyAs(fluid.phases().size())
void correctBoundaryConditions()
Correct boundary field.
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
static const Foam::dimensionedScalar A("", Foam::dimPressure, 611.21)
PtrList< surfaceScalarField > rAlphaAUfs(fluid.phases().size())
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