FreeFOAM The Cross-Platform CFD Toolkit
pEqn.H
Go to the documentation of this file.
1 rho = thermo.rho();
2 
4 U = rUA*UEqn.H();
5 
7 (
8  "phid",
10  *(
11  (fvc::interpolate(U) & mesh.Sf())
12  + fvc::ddtPhiCorr(rUA, rho, U, phi)
13  )
14 );
15 
16 for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
17 {
18  fvScalarMatrix pEqn
19  (
20  fvm::ddt(psi, p)
21  + fvm::div(phid, p)
22  - fvm::laplacian(rho*rUA, p)
23  );
24 
25  pEqn.solve();
26 
27  if (nonOrth == nNonOrthCorr)
28  {
29  phi = pEqn.flux();
30  }
31 }
32 
33 #include <finiteVolume/rhoEqn.H>
35 
36 U -= rUA*fvc::grad(p);
37 U.correctBoundaryConditions();