pEqn.H
Go to the documentation of this file.
1 volVectorField HbyA("HbyA", U);
3 {
4  HbyA = trTU()&UEqn().H();
5 }
6 else
7 {
8  HbyA = trAU()*UEqn().H();
9 }
10 
11 UEqn.clear();
13 
14 MRF.makeRelative(phiHbyA);
15 
16 adjustPhi(phiHbyA, U, p);
17 
18 while (simple.correctNonOrthogonal())
19 {
20  tmp<fvScalarMatrix> tpEqn;
21 
23  {
24  tpEqn = (fvm::laplacian(trTU(), p) == fvc::div(phiHbyA));
25  }
26  else
27  {
28  tpEqn = (fvm::laplacian(trAU(), p) == fvc::div(phiHbyA));
29  }
30 
31  tpEqn().setReference(pRefCell, pRefValue);
32 
33  tpEqn().solve();
34 
35  if (simple.finalNonOrthogonalIter())
36  {
37  phi = phiHbyA - tpEqn().flux();
38  }
39 }
40 
41 #include "continuityErrs.H"
42 
43 // Explicitly relax pressure for momentum corrector
44 p.relax();
45 
47 {
48  U = HbyA - (trTU()&fvc::grad(p));
49 }
50 else
51 {
52  U = HbyA - (trAU()*fvc::grad(p));
53 }
54 
55 U.correctBoundaryConditions();
56 fvOptions.correct(U);
phi
Definition: pEqn.H:20
tmp< volTensorField > trTU
Definition: UEqn.H:21
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcLaplacian.C:45
tmp< surfaceScalarField > interpolate(const RhoType &rho)
dynamicFvMesh & mesh
p
Definition: pEqn.H:59
fv::IOoptionList & fvOptions
const scalar pRefValue
IOMRFZoneList & MRF
fvVectorMatrix UEqn(fvm::ddt(rho, U)+fvm::div(phi, U)+turbulence->divDevRhoReff(U)==fvOptions(rho, U))
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
phiHbyA
Definition: pEqn.H:21
Switch pressureImplicitPorosity(false)
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:55
adjustPhi(phiHbyA, U, p_rgh)
const label pRefCell
const dictionary & simple
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcGrad.C:52
tmp< volScalarField > trAU
Definition: UEqn.H:20
U
Definition: pEqn.H:82
HbyA
Definition: pEqn.H:7