44 void Foam::solvers::multiphaseEuler::cellPressureCorrector()
54 PtrList<surfaceScalarField> alphafs(
phases.size());
57 const phaseModel& phase =
phases[phasei];
61 alphafs[phasei].rename(
"pEqn" + alphafs[phasei].
name());
72 PtrList<PtrList<volScalarField>> invADVs;
73 PtrList<PtrList<surfaceScalarField>> invADVfs;
84 UEqns[phase.index()].A()
87 max(phase.residualAlpha() -
alpha, scalar(0))
110 PtrList<surfaceScalarField> alphaByADfs;
111 PtrList<surfaceScalarField> FgByADfs;
113 PtrList<surfaceScalarField> Ffs(
fluid.
Fs());
140 + lalphafs[movingPhasei]
151 alphaByADfs = invADVfs & lalphafs;
152 FgByADfs = invADVfs & Fgfs;
163 PtrList<volVectorField> HbyADs;
175 UEqns[phase.index()].H()
178 max(phase.residualAlpha() -
alpha, scalar(0))
181 *phase.U()().oldTime()
184 if (HVms.set(movingPhasei))
186 Hs[movingPhasei] += HVms[movingPhasei];
190 HbyADs = invADVs & Hs;
209 alphaByADfs[movingPhasei]*mSfGradp
210 - FgByADfs[movingPhasei]
213 phase.URef().correctBoundaryConditions();
224 PtrList<volVectorField> HbyADs;
225 PtrList<surfaceScalarField> phiHbyADs;
242 UEqns[phase.index()].H()
245 max(phase.residualAlpha() -
alpha, scalar(0))
248 *phase.U()().oldTime()
251 if (HVms.set(movingPhasei))
253 Hs[movingPhasei] += HVms[movingPhasei];
259 fvc::flux(Hs[movingPhasei]) + ddtCorrs[phase.index()]
263 HbyADs = invADVs & Hs;
264 phiHbyADs = invADVfs & phiHs;
274 phiHbyADs[movingPhasei] -= FgByADfs[movingPhasei];
294 phiHbyA += alphafs[phase.index()]*phiHbyADs[movingPhasei];
317 rAf += alphafs[phase.index()]*alphaByADfs[movingPhasei];
334 alphafs[phase.index()].boundaryField()
335 *phase.phi()().boundaryField();
338 setSnGrad<fixedFluxPressureFvPatchScalarField>
348 PtrList<fvScalarMatrix> pEqnComps(compressibilityEqns(dmdts, d2mdtdps));
369 pEqn += pEqnComps[phasei];
398 phiHbyADs[movingPhasei]
399 + alphaByADfs[movingPhasei]*mSfGradp;
402 phase.divU(-pEqnComps[phase.index()] &
p_rgh);
411 mSfGradp = pEqnIncomp.flux()/rAf;
419 PtrList<volVectorField> dragCorrByADs
424 PtrList<surfaceScalarField> dragCorrByADfs
426 invADVfs & dragCorrfs
437 alphaByADfs[movingPhasei]*mSfGradp
438 - FgByADfs[movingPhasei]
439 + dragCorrByADfs[movingPhasei]
441 - dragCorrByADs[movingPhasei];
454 alphaByADfs[movingPhasei]*mSfGradp
455 - FgByADfs[movingPhasei]
464 phase.URef().correctBoundaryConditions();
493 phaseModel& phase =
phases_[phasei];
494 if (!phase.incompressible())
496 phase.rho() += phase.fluidThermo().psi()*(
p_rgh - p_rgh_0);
503 if (dmdts.set(phasei) && d2mdtdps.set(phasei))
505 dmdts[phasei] += d2mdtdps[phasei]*(
p_rgh - p_rgh_0);
#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.
void relax(const scalar alpha)
Relax field (for steady-state solution).
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.
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)
void makeRelative(volVectorField &U) const
Make the given absolute velocity relative within the MRF region.
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 > ddtCorrs() const =0
Return the flux corrections for the cell-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 void invADVs(const PtrList< volScalarField > &As, PtrList< volVectorField > &HVms, PtrList< PtrList< volScalarField >> &invADVs, PtrList< PtrList< surfaceScalarField >> &invADVfs) const =0
Return the inverse of the central + drag + virtual mass.
virtual void dragCorrs(PtrList< volVectorField > &dragCorrs, PtrList< surfaceScalarField > &dragCorrf) const =0
Set the cell and faces drag correction fields.
virtual PtrList< volScalarField > dmdts() const
Return the mass transfer rates for each phase.
virtual PtrList< surfaceScalarField > Fs() const =0
Return the force fluxes for the cell-based algorithm.
tmp< surfaceScalarField > surfaceTension(const phaseModel &phase) const
Return the surface tension force.
tmp< volScalarField > rho() const
Return the mixture density.
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
Switch predictMomentum
Momentum equation predictor switch.
surfaceScalarField & phi_
Switch dragCorrection
Cell/face drag correction for cell momentum corrector.
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)
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)
fvMatrix< scalar > fvScalarMatrix
const dimensionSet dimVolumetricFlux
SurfaceField< scalar > surfaceScalarField
tmp< volScalarField > byDt(const volScalarField &vf)
void constrainHbyA(volVectorField &HbyA, const volVectorField &U, const volScalarField &p)
const dimensionSet dimTime
const dimensionSet dimDensity
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.