27 PtrList<surfaceScalarField>
alphafs(
fluid.phases().size());
28 PtrList<volVectorField>
HbyAs(
fluid.phases().size());
30 PtrList<volScalarField>
rAUs(
fluid.phases().size());
36 phaseModel& phase = iter();
38 MRF.makeAbsolute(phase.phi().oldTime());
39 MRF.makeAbsolute(phase.phi());
67 phaseModel& phase = iter();
82 zeroGradientFvPatchScalarField::typeName
86 rAUs.set(phasei, (1.0/(
UEqns[phasei].A() + dragCoeffi)).ptr());
107 MRF.makeRelative(phase.phi().oldTime());
108 MRF.makeRelative(phase.phi());
118 multiphaseSystem::dragModelTable::const_iterator dmIter =
119 fluid.dragModels().begin();
120 multiphaseSystem::dragCoeffFields::const_iterator dcIter =
129 const phaseModel *phase2Ptr =
nullptr;
131 if (&phase == &dmIter()->phase1())
133 phase2Ptr = &dmIter()->phase2();
135 else if (&phase == &dmIter()->
phase2())
137 phase2Ptr = &dmIter()->phase1();
150 (1.0/phase.rho())*
rAUs[phasei]*(*dcIter())
185 phaseModel& phase = iter();
193 surfaceScalarField::Boundary
phib(
phi.boundaryField());
198 phaseModel& phase = iter();
202 *(
mesh.Sf().boundaryField() & phase.U().boundaryField());
209 p_rgh.boundaryFieldRef(),
212 )/(
mesh.magSf().boundaryField()*rAUf.boundaryField())
236 if (
pimple.finalNonOrthogonalIter())
244 phaseModel& phase = iter();
252 *mSfGradp/phase.rho();
268 mSfGradp = pEqnIncomp.flux()/
rAUf;
275 phaseModel& phase = iter();
292 phase.U().correctBoundaryConditions();
294 U += alpha*phase.U();
303 #include "continuityErrs.H" fvMatrix< scalar > fvScalarMatrix
pimpleNoLoopControl & pimple
PtrList< volScalarField > rAUs
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
tmp< GeometricField< typename flux< Type >::type, fvsPatchField, surfaceMesh > > ddtCorr(const GeometricField< Type, fvPatchField, volMesh > &U, const GeometricField< Type, fvsPatchField, surfaceMesh > &Uf)
dimensioned< vector > dimensionedVector
Dimensioned vector obtained from generic dimensioned type.
autoPtr< multiphaseSystem::dragCoeffFields > dragCoeffs(fluid.dragCoeffs())
GeometricField< vector, fvPatchField, volMesh > volVectorField
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
forAllIter(PtrDictionary< phaseModel >, fluid.phases(), iter)
const surfaceScalarField & ghf
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.
PtrList< volVectorField > HbyAs(fluid.phases().size())
PtrList< surfaceScalarField > alphafs(phases.size())
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())
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()))
PtrList< surfaceScalarField > phiHbyAs(fluid.phases().size())
void correctBoundaryConditions()
Correct boundary field.
const volScalarField & gh
dimensioned< scalar > mag(const dimensioned< Type > &)
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const dimensionedVector & g
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
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))