applications/solvers/multiphase/compressibleInterDyMFoam/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) - divU
  28         );
  29
  30         pcorrEqn.solve();
  31
  32         if (nonOrth == nNonOrthCorr)
  33         {
  34             phi -= pcorrEqn.flux();
  35         }
  36     }
  37
  38     #include "continuityErrs.H"
  39 }