applications/solvers/coupled/MRFPorousFoam/pEqn.H

   1 // Pressure parts of the continuity equation
   2 surfaceScalarField rUAf
   3 (
   4     "rUAf",
   5     fvc::interpolate(1.0/UEqn().A())
   6 );
   7
   8 UEqn.clear();
   9
  10 surfaceScalarField presSource
  11 (
  12     "presSource",
  13     rUAf*(fvc::interpolate(fvc::grad(p, "grad(pSource)")) & mesh.Sf())
  14 );
  15
  16 fvScalarMatrix pEqn
  17 (
  18   - fvm::laplacian(rUAf, p)
  19  ==
  20   - fvc::div(presSource)
  21 );
  22
  23 pEqn.setReference(pRefCell, pRefValue);
  24
  25 UpEqn.insertEquation(3, pEqn);