applications/solvers/multiphase/multiphaseInterFoam/pEqn.H

   1 {
   2     volScalarField rUA = 1.0/UEqn.A();
   3     surfaceScalarField rUAf = fvc::interpolate(rUA);
   4
   5     U = rUA*UEqn.H();
   6
   7     surfaceScalarField phiU
   8     (
   9         "phiU",
  10         (fvc::interpolate(U) & mesh.Sf()) + fvc::ddtPhiCorr(rUA, rho, U, phi)
  11     );
  12
  13     adjustPhi(phiU, U, p);
  14
  15     phi = phiU +
  16         (
  17             mixture.surfaceTensionForce()*mesh.magSf()
  18           + fvc::interpolate(rho)*(g & mesh.Sf())
  19         )*rUAf;
  20
  21     for(int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
  22     {
  23         fvScalarMatrix pEqn
  24         (
  25             fvm::laplacian(rUAf, p) == fvc::div(phi)
  26         );
  27
  28         pEqn.setReference(pRefCell, pRefValue);
  29
  30         if (corr == nCorr-1)
  31         {
  32             pEqn.solve(mesh.solver(p.name() + "Final"));
  33         }
  34         else
  35         {
  36             pEqn.solve(mesh.solver(p.name()));
  37         }
  38
  39         if (nonOrth == nNonOrthCorr)
  40         {
  41             phi -= pEqn.flux();
  42         }
  43     }
  44
  45     U += rUA*fvc::reconstruct((phi - phiU)/rUAf);
  46     U.correctBoundaryConditions();
  47 }