pEqn.H
Go to the documentation of this file.
1rho = thermo.rho();
2
3volScalarField rAU(1.0/UEqn.A());
4surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU));
5volVectorField HbyA(constrainHbyA(rAU*UEqn.H(), U, p));
6surfaceScalarField phid
7(
8 "phid",
9 fvc::interpolate(psi)
10 *(
11 fvc::flux(HbyA)
12 + MRF.zeroFilter(rhorAUf*fvc::ddtCorr(rho, U, phi)/fvc::interpolate(rho))
13 )
14);
15
16MRF.makeRelative(fvc::interpolate(psi), phid);
17
18// Non-orthogonal pressure corrector loop
19while (pimple.correctNonOrthogonal())
20{
21 fvScalarMatrix pEqn
22 (
23 fvm::ddt(psi, p)
24 + fvm::div(phid, p)
25 - fvm::laplacian(rhorAUf, p)
26 ==
27 fvOptions(psi, p, rho.name())
28 );
29
30 pEqn.solve();
31
32 if (pimple.finalNonOrthogonalIter())
33 {
34 phi = pEqn.flux();
35 }
36}
37
38#include "rhoEqn.H"
39#include "compressibleContinuityErrs.H"
40
41U = HbyA - rAU*fvc::grad(p);
42U.correctBoundaryConditions();
43fvOptions.correct(U);
44K = 0.5*magSqr(U);
CGAL::Exact_predicates_exact_constructions_kernel K
fv::options & fvOptions
surfaceScalarField & phi
IOMRFZoneList & MRF
pimpleControl & pimple
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
U
Definition: pEqn.H:72
volScalarField & p
const volScalarField & psi
surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho *rAU))
fvVectorMatrix & UEqn
Definition: UEqn.H:13
HbyA
Definition: pcEqn.H:74
surfaceScalarField phid("phid", fvc::interpolate(psi) *(fvc::flux(HbyA)+MRF.zeroFilter(rhorAUf *fvc::ddtCorr(rho, U, phi)/fvc::interpolate(rho))))
tmp< volScalarField > rAU
Definition: initCorrectPhi.H:1