19 while (
simple.correctNonOrthogonal())
21 tmp<fvScalarMatrix> tpEqn;
36 pressureReference.refCell(),
37 pressureReference.refValue()
42 if (
simple.finalNonOrthogonalIter())
62 U.correctBoundaryConditions();
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
tmp< volTensorField > trTU
tmp< fvVectorMatrix > tUEqn(fvm::div(phi, U)+MRF.DDt(rho, U)+turbulence->divDevTau(U)==fvModels.source(rho, U))
tmp< volScalarField > trAU
Calculates and prints the continuity errors.
Foam::fvConstraints & fvConstraints(Foam::fvConstraints::New(mesh))
simpleControl simple(mesh)
Switch pressureImplicitPorosity(false)
tmp< volVectorField > tHbyA
surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(rho) *fvc::flux(HbyA))
tmp< SurfaceField< typename innerProduct< vector, Type >::type > > flux(const VolField< Type > &vf)
Return the face-flux field obtained from the given volVectorField.
tmp< VolField< Type > > laplacian(const VolField< Type > &vf, const word &name)
tmp< VolField< typename outerProduct< vector, Type >::type > > grad(const SurfaceField< Type > &ssf)
tmp< VolField< Type > > div(const SurfaceField< Type > &ssf)
VolField< vector > volVectorField
fvMatrix< scalar > fvScalarMatrix
SurfaceField< scalar > surfaceScalarField
void constrainHbyA(volVectorField &HbyA, const volVectorField &U, const volScalarField &p)