78 phiHbyAD += rADUcf*(Fgf +
Dcf()*
phid());
For cases which do no have a pressure boundary adjust the balance of fluxes to obey continuity....
Generic GeometricField class.
void relax(const scalar alpha)
Relax field (for steady-state solution).
bool needReference() const
Does the field need a reference level for solution.
void correctBoundaryConditions()
Correct boundary field.
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
A special matrix type and solver, designed for finite volume solutions of scalar equations....
SolverPerformance< Type > solve(const dictionary &)
Solve segregated or coupled returning the solution statistics.
tmp< volScalarField > A() const
Return the central coefficient.
tmp< VolField< Type > > H() const
Return the H operation source.
tmp< SurfaceField< Type > > flux() const
Return the face-flux field from the matrix.
void setReference(const label celli, const Type &value, const bool forceReference=false)
Set reference level for solution.
const surfaceVectorField & Sf() const
Return cell face area vectors.
bool finalNonOrthogonalIter() const
Flag to indicate the last non-orthogonal iteration.
bool correctNonOrthogonal()
Non-orthogonal corrector loop.
Provides controls for the pressure reference in closed-volume simulations.
scalar refValue() const
Return the pressure reference level.
label refCell() const
Return the cell in which the reference pressure is set.
Foam::fvModels & fvModels() const
Return the fvModels that are created on demand.
pimpleNoLoopControl pimple
PIMPLE inner-loop controls.
Foam::fvConstraints & fvConstraints() const
Return the fvConstraints that are created on demand.
const fvMesh & mesh
Region mesh.
autoPtr< surfaceScalarField > phid
Effective volumetric flux of the dispersed phase.
volScalarField p_
Pressure field.
uniformDimensionedVectorField g
Acceleration due to gravity.
const volVectorField & Uc
Reference to the continuous phase velocity field.
surfaceScalarField phic_
Continuous phase flux field.
surfaceScalarField alphacf
Interpolated continuous phase-fraction.
const volScalarField & alphac
Reference continuous phase-fraction.
volVectorField Uc_
Continuous phase velocity field.
autoPtr< surfaceVectorField > Ucf
Pointer to the surface momentum field.
tmp< fvVectorMatrix > tUcEqn
Cached momentum matrix.
autoPtr< volVectorField > Fd
Dispersed phase drag force.
virtual void correctPressure()
Construct the pressure equation.
autoPtr< volScalarField > Dc
autoPtr< surfaceScalarField > Dcf
Face continuous-dispersed phase drag coefficient.
const surfaceScalarField & phic
Reference to the continuous phase volumetric-flux field.
const volScalarField & p
Reference to the pressure field.
Calculate the first temporal derivative.
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 matrix for the laplacian of the field.
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< Type > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
tmp< VolField< typename outerProduct< vector, Type >::type > > reconstruct(const SurfaceField< Type > &ssf)
void correctUf(autoPtr< surfaceVectorField > &Uf, const volVectorField &U, const surfaceScalarField &phi, const MRFType &MRF)
tmp< VolField< Type > > div(const SurfaceField< Type > &ssf)
void makeAbsolute(surfaceScalarField &phi, const volVectorField &U)
Make the given flux absolute.
void makeRelative(surfaceScalarField &phi, const volVectorField &U)
Make the given flux relative.
tmp< SurfaceField< typename Foam::flux< Type >::type > > ddtCorr(const VolField< Type > &U, const SurfaceField< Type > &Uf)
tmp< fvMatrix< Type > > laplacian(const VolField< Type > &vf, const word &name)
void constrainPressure(volScalarField &p, const RhoType &rho, const volVectorField &U, const surfaceScalarField &phiHbyA, const RAUType &rhorAU, const MRFType &MRF)
bool adjustPhi(surfaceScalarField &phi, const volVectorField &U, volScalarField &p)
Adjust the balance of fluxes to obey continuity.
void constrainHbyA(volVectorField &HbyA, const volVectorField &U, const volScalarField &p)