pEqn.H
Go to the documentation of this file.
1{
2 surfaceScalarField faceMask(localMin<scalar>(mesh).interpolate(cellMask));
3
4 volScalarField rAU(1.0/UEqn.A());
5 surfaceScalarField rAUf("rAUf", faceMask*fvc::interpolate(rAU));
6
7 volVectorField HbyA("HbyA", U);
8 HbyA = constrainHbyA(cellMask*rAU*UEqn.H(), U, p);
9
10 //mesh.interpolate(HbyA);
12 {
13 #include "interpolatedFaces.H"
14 }
15
16 tUEqn.clear();
17
18 surfaceScalarField phiHbyA("phiHbyA", fvc::flux(HbyA));
19
21
23 {
24 oversetAdjustPhi(phiHbyA, U);
25 }
26
27 // Non-orthogonal pressure corrector loop
28 while (simple.correctNonOrthogonal())
29 {
30 fvScalarMatrix pEqn
31 (
32 fvm::laplacian(rAUf, p) == fvc::div(phiHbyA)
33 );
34
35 pEqn.setReference(pRefCell, pRefValue);
36
37 pEqn.solve();
38
39 if (simple.finalNonOrthogonalIter())
40 {
41 phi = phiHbyA - pEqn.flux();
42 }
43 }
44
45 #include "continuityErrs.H"
46
47 // Explicitly relax pressure for momentum corrector
48 p.relax();
49
50 // Momentum corrector
51 volVectorField gradP(fvc::grad(p));
52 //mesh.interpolate(gradP);
53
54 U = HbyA - rAU*cellMask*gradP;
55 U.correctBoundaryConditions();
56 fvOptions.correct(U);
57}
fv::options & fvOptions
surfaceScalarField & phi
const scalar pRefValue
const label pRefCell
U
Definition: pEqn.H:72
volScalarField & p
tmp< fvVectorMatrix > tUEqn(fvm::ddt(rho, U)+fvm::div(phi, U)+MRF.DDt(rho, U)+turbulence->divDevRhoReff(U)==fvOptions(rho, U))
fvVectorMatrix & UEqn
Definition: UEqn.H:13
phiHbyA
Definition: pcEqn.H:73
HbyA
Definition: pcEqn.H:74
surfaceScalarField faceMask(localMin< scalar >(mesh).interpolate(cellMask))
volVectorField gradP(fvc::grad(p))
mesh interpolate(rAU)
dynamicFvMesh & mesh
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU))
adjustPhi(phiHbyA, U, p_rgh)
massFluxInterpolation
Definition: readControls.H:7
tmp< volScalarField > rAU
Definition: initCorrectPhi.H:1
adjustFringe
Definition: readControls.H:16
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const dictionary & simple
Calculates and prints the continuity errors.