44 void Foam::solvers::multiphaseEuler::facePressureCorrector()
49 PtrList<surfaceScalarField> alphafs(
phases.size());
52 const phaseModel& phase =
phases[phasei];
56 alphafs[phasei].rename(
"pEqn" + alphafs[phasei].
name());
67 PtrList<PtrList<surfaceScalarField>> invADVfs;
80 max(
alpha.oldTime(), phase.residualAlpha())
81 *phase.rho().oldTime()
83 +
UEqns[phase.index()].A()
116 PtrList<surfaceScalarField> alphaByADfs;
117 PtrList<surfaceScalarField> FgByADfs;
120 PtrList<surfaceScalarField> Ffs(
fluid.
Ffs());
138 max(alphafs[phase.index()], phase.residualAlpha())
145 + lalphafs[movingPhasei]
155 alphaByADfs = invADVfs & lalphafs;
156 FgByADfs = invADVfs & Fgfs;
171 PtrList<surfaceScalarField> phiHbyADs;
186 max(
alpha.oldTime(), phase.residualAlpha())
187 *phase.rho().oldTime()
199 if (HVmfs.set(movingPhasei))
201 phiHs[movingPhasei] += HVmfs[movingPhasei];
205 phiHbyADs = invADVfs & phiHs;
214 phiHbyADs[movingPhasei] -= FgByADfs[movingPhasei];
236 phiHbyA += alphafs[phase.index()]*phiHbyADs[movingPhasei];
259 rAf += alphafs[phase.index()]*alphaByADfs[movingPhasei];
276 alphafs[phase.index()].boundaryField()
277 *phase.phi()().boundaryField();
280 setSnGrad<fixedFluxPressureFvPatchScalarField>
290 PtrList<fvScalarMatrix> pEqnComps(compressibilityEqns(dmdts, d2mdtdps));
311 pEqn += pEqnComps[phasei];
338 const label phasei = phase.index();
341 phiHbyADs[movingPhasei]
342 + alphaByADfs[movingPhasei]*mSfGradp;
345 phase.divU(-pEqnComps[phasei] &
p_rgh);
360 phase.URef().correctBoundaryConditions();
389 phaseModel& phase =
phases_[phasei];
390 if (!phase.incompressible())
392 phase.rho() += phase.fluidThermo().psi()*(
p_rgh - p_rgh_0);
399 if (dmdts.set(phasei) && d2mdtdps.set(phasei))
401 dmdts[phasei] += d2mdtdps[phasei]*(
p_rgh - p_rgh_0);
static const Foam::dimensionedScalar A("A", Foam::dimPressure, 611.21)
#define forAll(list, i)
Loop across all elements in list.
const dimensionSet & dimensions() const
Return dimensions.
static const DimensionedField< Type, GeoMesh > & null()
Return a null DimensionedField.
GeometricBoundaryField< Type, PatchField, GeoMesh > Boundary
Type of the boundary field.
const Boundary & boundaryField() const
Return const-reference to the boundary field.
bool needReference() const
Does the field need a reference level for solution.
void correctBoundaryConditions()
Correct boundary field.
Boundary & boundaryFieldRef()
Return a reference to the boundary field.
IOobject(const word &name, const fileName &instance, const objectRegistry ®istry, readOption r=NO_READ, writeOption w=NO_WRITE, bool registerObject=true)
Construct from name, instance, registry, io options.
const word & name() const
Return name.
static word groupName(Name name, const word &group)
tmp< surfaceScalarField > absolute(const tmp< surfaceScalarField > &phi) const
Return the given relative flux absolute within the MRF region.
void makeRelative(volVectorField &U) const
Make the given absolute velocity relative within the MRF region.
const FieldType & oldTime() const
Return the old time field.
bool set(const label) const
Is element set.
label size() const
Return the number of elements in the UPtrList.
const word & name() const
Return const reference to name.
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
const surfaceVectorField & Sf() const
Return cell face area vectors.
const surfaceScalarField & magSf() const
Return cell face area magnitudes.
bool finalNonOrthogonalIter() const
Flag to indicate the last non-orthogonal iteration.
virtual PtrList< surfaceScalarField > Ffs() const =0
Return the force fluxes for the face-based algorithm.
void correctBoundaryFlux()
Correct fixed-flux BCs to be consistent with the velocity BCs.
virtual bool implicitPhasePressure(const phaseModel &phase) const
Returns true if the phase pressure is treated implicitly.
virtual PtrList< volScalarField > dmdts() const
Return the mass transfer rates for each phase.
tmp< surfaceScalarField > surfaceTension(const phaseModel &phase) const
Return the surface tension force.
tmp< volScalarField > rho() const
Return the mixture density.
virtual PtrList< PtrList< surfaceScalarField > > invADVfs(const PtrList< surfaceScalarField > &Afs, PtrList< surfaceScalarField > &HVmfs) const =0
Return the inverse of the face central + drag + virtual mass.
virtual PtrList< volScalarField > d2mdtdps() const
Return the mass transfer pressure implicit coefficients.
bool incompressible() const
Return incompressibility.
bool correct()
Piso loop within outer loop.
bool correctNonOrthogonal()
Non-orthogonal corrector loop.
scalar refValue() const
Return the pressure reference level.
label refCell() const
Return the cell in which the reference pressure is set.
pimpleNoLoopControl pimple
PIMPLE inner-loop controls.
Foam::fvConstraints & fvConstraints() const
Return the fvConstraints that are created on demand.
const Time & runTime
Time.
const fvMesh & mesh
Region mesh.
uniformDimensionedVectorField g
Gravitational acceleration.
volScalarField gh
(g & h) - ghRef
uniformDimensionedScalarField pRef
Reference pressure.
surfaceScalarField ghf
(g & hf) - ghRef
volScalarField & p_rgh
Reference to the buoyant pressure for buoyant cases.
const surfaceScalarField & phi
Reference to the mass-flux field.
solvers::buoyancy buoyancy
Buoyancy force.
phaseSystem::phaseModelList & phases_
PtrList< fvVectorMatrix > UEqns
Temporary phase momentum matrices.
const phaseSystem::phaseModelPartialList & movingPhases
Reference to the moving phases.
const IOMRFZoneList & MRF
surfaceScalarField & phi_
Foam::pressureReference pressureReference
Pressure reference.
PtrList< volScalarField > rAs
Temporary storage for the reciprocal momentum equation diagonal.
const phaseSystem::phaseModelList & phases
Reference to the phases.
phaseSystem::phaseModelPartialList & movingPhases_
const volScalarField & p
Reference to the pressure field.
const phaseSystem & fluid
Reference to the multiphase fluid.
Find the reference cell nearest (in index) to the given cell but which is not on a cyclic,...
Calculate the first temporal derivative.
Calculate the divergence of the given field.
Calculate the face-flux of the given field.
Calculate the mesh motion flux and convert fluxes from absolute to relative and back.
Reconstruct volField from a face flux field.
Calculate the snGrad of the given volField.
Calculate the field for explicit evaluation of implicit and explicit sources.
Calculate the matrix for the first temporal derivative.
Calculate the matrix for the divergence of the given field and flux.
Calculate the matrix for the laplacian of the field.
Calculate the matrix for implicit and explicit sources.
surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(rho) *fvc::flux(HbyA))
volScalarField alpha(IOobject("alpha", runTime.name(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
tmp< SurfaceField< typename innerProduct< vector, Type >::type > > flux(const VolField< Type > &vf)
Return the face-flux field obtained from the given volVectorField.
static tmp< SurfaceField< Type > > interpolate(const VolField< Type > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
tmp< VolField< typename outerProduct< vector, Type >::type > > reconstruct(const SurfaceField< Type > &ssf)
tmp< VolField< Type > > div(const SurfaceField< Type > &ssf)
tmp< surfaceScalarField > absolute(const tmp< surfaceScalarField > &tphi, const volVectorField &U)
Return the given relative flux in absolute form.
void makeRelative(surfaceScalarField &phi, const volVectorField &U)
Make the given flux relative.
tmp< SurfaceField< Type > > snGrad(const VolField< Type > &vf, const word &name)
tmp< fvMatrix< Type > > laplacian(const VolField< Type > &vf, const word &name)
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
fvMatrix< scalar > fvScalarMatrix
const dimensionSet dimVolumetricFlux
SurfaceField< scalar > surfaceScalarField
tmp< volScalarField > byDt(const volScalarField &vf)
VolField< scalar > volScalarField
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
scalar getRefCellValue(const volScalarField &field, const label refCelli)
Return the current value of field in the reference cell.
void constrainPhiHbyA(surfaceScalarField &phiHbyA, const volVectorField &U, const volScalarField &p)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.