FreeFOAM The Cross-Platform CFD Toolkit
pEqn.H
Go to the documentation of this file.
1 rho = thermo.rho();
2 
4 surfaceScalarField rhorUAf("(rho*(1|A(U)))", fvc::interpolate(rho*rUA));
5 U = rUA*UEqn.H();
6 
8 (
10  *(
11  (fvc::interpolate(U) & mesh.Sf())
12  + fvc::ddtPhiCorr(rUA, rho, U, phi)
13  )
14 );
15 
17 
18 for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
19 {
21 
22  fvScalarMatrix p_rghEqn
23  (
25  + fvc::div(phi)
26  - fvm::laplacian(rhorUAf, p_rgh)
27  );
28 
29  p_rghEqn.solve
30  (
31  mesh.solver(p_rgh.select(corr == nCorr-1 && nonOrth == nNonOrthCorr))
32  );
33 
34  if (nonOrth == nNonOrthCorr)
35  {
36  phi += p_rghEqn.flux();
37  }
38 }
39 
40 p = p_rgh + rho*gh;
41 
42 #include <finiteVolume/rhoEqn.H>
44 
45 U += rUA*fvc::reconstruct((phi - phiU)/rhorUAf);
46 U.correctBoundaryConditions();
47