pEqn.H
Go to the documentation of this file.
1{
3
4 surfaceScalarField phiHbyA
5 (
6 "phiHbyA",
7 (
8 fvc::flux(HbyA)
9 + alphacf*rAUcf*fvc::ddtCorr(Uc, phic)
10 )
11 );
12
13 if (p.needReference())
14 {
15 adjustPhi(phiHbyA, Uc, p);
16 }
17
19
20 // Update the pressure BCs to ensure flux consistency
22
23 // Non-orthogonal pressure corrector loop
24 while (pimple.correctNonOrthogonal())
25 {
26 fvScalarMatrix pEqn
27 (
28 fvm::laplacian(alphacf*rAUcf, p)
29 ==
30 fvc::ddt(alphac) + fvc::div(alphacf*phiHbyA)
31 );
32
33 pEqn.setReference(pRefCell, pRefValue);
34
35 pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
36
37 if (pimple.finalNonOrthogonalIter())
38 {
39 phic = phiHbyA - pEqn.flux()/alphacf;
40
41 p.relax();
42
43 Uc = HbyA
44 + rAUc*fvc::reconstruct((phicForces - pEqn.flux()/alphacf)/rAUcf);
45 Uc.correctBoundaryConditions();
46 }
47 }
48}
49
50#include "continuityErrs.H"
surfaceScalarField phic(mixture.cAlpha() *mag(alphaPhic/mesh.magSf()))
surfaceScalarField rAUcf("Dp", fvc::interpolate(rAUc))
volScalarField rAUc(1.0/UcEqn.A())
surfaceScalarField phicForces(fvc::flux(rAUc *cloudVolSUSu/rhoc)+rAUcf *(g &mesh.Sf()))
fvVectorMatrix UcEqn(fvm::ddt(alphac, Uc)+fvm::div(alphaPhic, Uc) - fvm::Sp(fvc::ddt(alphac)+fvc::div(alphaPhic), Uc)+continuousPhaseTurbulence->divDevRhoReff(Uc)==(1.0/rhoc) *cloudSU)
const scalar pRefValue
const label pRefCell
pimpleControl & pimple
volScalarField & p
phiHbyA
Definition: pcEqn.H:73
HbyA
Definition: pcEqn.H:74
dynamicFvMesh & mesh
adjustPhi(phiHbyA, U, p_rgh)
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:83
Calculates and prints the continuity errors.