pEqn.H
Go to the documentation of this file.
1 volScalarField rAU(1.0/UEqn.A());
4 
5 if (pimple.nCorrPISO() <= 1)
6 {
7  tUEqn.clear();
8 }
9 
11 (
12  "phiHbyA",
14  + MRF.zeroFilter(rAUf*fvc::ddtCorr(U, phi))
15 );
16 
17 MRF.makeRelative(phiHbyA);
18 adjustPhi(phiHbyA, U, p_gh);
19 
20 // Update the pressure BCs to ensure flux consistency
22 
23 // Non-orthogonal pressure corrector loop
24 while (pimple.correctNonOrthogonal())
25 {
26  fvScalarMatrix p_ghEqn
27  (
29  );
30 
31  p_ghEqn.setReference(p_ghRefCell, p_ghRefValue);
32 
33  p_ghEqn.solve(mesh.solver(p_gh.select(pimple.finalInnerIter())));
34 
35  if (pimple.finalNonOrthogonalIter())
36  {
37  phi = phiHbyA - p_ghEqn.flux();
38  }
39 }
40 
41 #include "continuityErrs.H"
42 
43 // Explicitly relax pressure for momentum corrector
44 p_gh.relax();
45 
46 p = p_gh + (g & mesh.C());
47 
48 U = HbyA - rAU*fvc::grad(p_gh);
49 U.correctBoundaryConditions();
50 fvOptions.correct(U);
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcGrad.C:52
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
fv::options & fvOptions
pimpleNoLoopControl & pimple
p
Definition: pEqn.H:50
volScalarField rAU(1.0/UEqn.A())
IOMRFZoneList & MRF
phiHbyA
Definition: pEqn.H:20
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
tmp< GeometricField< typename flux< Type >::type, fvsPatchField, surfaceMesh > > ddtCorr(const GeometricField< Type, fvPatchField, volMesh > &U, const GeometricField< Type, fvsPatchField, surfaceMesh > &Uf)
Definition: fvcDdt.C:170
tmp< fvVectorMatrix > tUEqn(fvm::ddt(rho, U)+fvm::div(phi, U)+MRF.DDt(rho, U)+turbulence->divDevRhoReff(U)==fvOptions(rho, U))
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:55
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcLaplacian.C:45
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
dynamicFvMesh & mesh
tmp< volVectorField > constrainHbyA(const tmp< volVectorField > &tHbyA, const volVectorField &U, const volScalarField &p)
Definition: constrainHbyA.C:33
adjustPhi(phiHbyA, U, p_rgh)
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.
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU))
U
Definition: pEqn.H:72
volVectorField & HbyA
Definition: pEqn.H:13
fvVectorMatrix & UEqn
Definition: UEqn.H:13
constrainPressure(p_rgh, rho, U, phiHbyA, rhorAUf, MRF)
phi
Definition: pEqn.H:18
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
Definition: fvcFlux.C:32
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const dimensionedVector & g