FreeFOAM The Cross-Platform CFD Toolkit
pEqn.H
Go to the documentation of this file.
2 
4 (
5  "(rho*(1|A(U)))",
7 );
8 
9 U = rUA*UEqn.H();
10 phi =
12  *(
13  (fvc::interpolate(U) & mesh.Sf())
15  );
16 
17 surfaceScalarField phiU("phiU", phi);
18 phi -= ghf*fvc::snGrad(rho)*rUAf*mesh.magSf();
19 
20 for(int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
21 {
22  fvScalarMatrix p_rghEqn
23  (
25  );
26 
27  p_rghEqn.setReference(pRefCell, getRefCellValue(p_rgh, pRefCell));
28  p_rghEqn.solve();
29 
30  if (nonOrth == nNonOrthCorr)
31  {
32  phi -= p_rghEqn.flux();
33  }
34 }
35 
36 p == p_rgh + rho*gh;
37 
38 if (p_rgh.needReference())
39 {
41  (
42  "p",
43  p.dimensions(),
45  );
46  p_rgh = p - rho*gh;
47 }
48 
49 #include <finiteVolume/rhoEqn.H>
51 
53 U.correctBoundaryConditions();
54 
55 // ************************ vim: set sw=4 sts=4 et: ************************ //