applications/solvers/multiphase/interDyMFoam/correctPhi.H

   1 {
   2     #include "continuityErrs.H"
   3
   4     volScalarField pcorr
   5     (
   6         IOobject
   7         (
   8             "pcorr",
   9             runTime.timeName(),
  10             mesh,
  11             IOobject::NO_READ,
  12             IOobject::NO_WRITE
  13         ),
  14         mesh,
  15         dimensionedScalar("pcorr", p.dimensions(), 0.0),
  16         pcorrTypes
  17     );
  18
  19     dimensionedScalar rAUf("(1|A(U))", dimTime/rho.dimensions(), 1.0);
  20
  21     adjustPhi(phi, U, pcorr);
  22
  23     for(int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
  24     {
  25         fvScalarMatrix pcorrEqn
  26         (
  27             fvm::laplacian(rAUf, pcorr) == fvc::div(phi)
  28         );
  29
  30         pcorrEqn.setReference(pRefCell, pRefValue);
  31         pcorrEqn.solve();
  32
  33         if (nonOrth == nNonOrthCorr)
  34         {
  35             phi -= pcorrEqn.flux();
  36         }
  37     }
  38
  39     #include "continuityErrs.H"
  40 }