applications/solvers/multiphase/interFoam/interDyMFoam/correctPhi.H

   1 {
   2     #include "continuityErrs.H"
   3
   4     wordList pcorrTypes
   5     (
   6         p_rgh.boundaryField().size(),
   7         zeroGradientFvPatchScalarField::typeName
   8     );
   9
  10     forAll (p_rgh.boundaryField(), i)
  11     {
  12         if (p_rgh.boundaryField()[i].fixesValue())
  13         {
  14             pcorrTypes[i] = fixedValueFvPatchScalarField::typeName;
  15         }
  16     }
  17
  18     volScalarField pcorr
  19     (
  20         IOobject
  21         (
  22             "pcorr",
  23             runTime.timeName(),
  24             mesh,
  25             IOobject::NO_READ,
  26             IOobject::NO_WRITE
  27         ),
  28         mesh,
  29         dimensionedScalar("pcorr", p_rgh.dimensions(), 0.0),
  30         pcorrTypes
  31     );
  32
  33     dimensionedScalar rAUf("(1|A(U))", dimTime/rho.dimensions(), 1.0);
  34
  35     adjustPhi(phi, U, pcorr);
  36
  37     fvc::makeAbsolute(phi, U);
  38
  39     for (int nonOrth=0; nonOrth<=pimple.nNonOrthCorr(); nonOrth++)
  40     {
  41         fvScalarMatrix pcorrEqn
  42         (
  43             fvm::laplacian(rAUf, pcorr) == fvc::div(phi)
  44         );
  45
  46         pcorrEqn.setReference(pRefCell, pRefValue);
  47         pcorrEqn.solve();
  48
  49         if (nonOrth == pimple.nNonOrthCorr())
  50         {
  51             phi -= pcorrEqn.flux();
  52             phiAbs = phi;
  53             fvc::makeRelative(phi, U);
  54         }
  55     }
  56
  57     #include "continuityErrs.H"
  58 }