applications/solvers/multiphase/barotropicCavitatingFoam/pEqn.H

   1 {
   2     surfaceScalarField rUAf = 1.0/fvc::interpolate(UEqn.A());
   3
   4     U = UEqn.H()/UEqn.A();
   5
   6     phi = fvc::interpolate(U) & mesh.Sf();
   7
   8     // Bug fix: must change name of phi on copy to keep objectRegistry happy
   9     // HJ, 7/Nov/2010
  10     surfaceScalarField phiU
  11     (
  12         "phiU",
  13         phi
  14     );
  15
  16     surfaceScalarField phip = fvc::interpolate(psiByRho)*phi;
  17
  18     for(int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
  19     {
  20         fvScalarMatrix pEqn
  21         (
  22             fvm::laplacian(rUAf, p)
  23         );
  24
  25         solve
  26         (
  27              psiByRho*fvm::ddt(p)
  28            + fvm::div(phip, p) - fvm::Sp(fvc::div(phip), p)
  29            + fvc::div(phi) - pEqn
  30         );
  31
  32         if (nonOrth == nNonOrthCorr)
  33         {
  34             phi -= pEqn.flux();
  35         }
  36     }
  37
  38     U += fvc::reconstruct(phi - phiU);
  39     U.correctBoundaryConditions();
  40 }