src/navalFoamVariant2/correctPhi.H

   1 {
   2 #   include "continuityErrs.H"
   3
   4     wordList pcorrTypes(pd.boundaryField().types());
   5
   6     for (label i=0; i < pd.boundaryField().size(); i++)
   7     {
   8         if (pd.boundaryField()[i].fixesValue())
   9         {
  10             pcorrTypes[i] = fixedValueFvPatchScalarField::typeName;
  11         }
  12     }
  13
  14     volScalarField pcorr
  15     (
  16         IOobject
  17         (
  18             "pcorr",
  19             runTime.timeName(),
  20             mesh,
  21             IOobject::NO_READ,
  22             IOobject::NO_WRITE
  23         ),
  24         mesh,
  25         dimensionedScalar("pcorr", pd.dimensions(), 0.0),
  26         pcorrTypes
  27     );
  28
  29     // Flux predictor: evaluate velocity on moving boundary and assemble new
  30     // fluxes.  Correct the fluxes to achieve zero divergence.
  31     // HJ, experimental 21/Mar/2011
  32     U.correctBoundaryConditions();
  33
  34     phi = fvc::interpolate(U) & mesh.Sf();
  35
  36     for(int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
  37     {
  38         fvScalarMatrix pcorrEqn
  39         (
  40             fvm::laplacian(rAU/rho, pcorr) == fvc::div(phi)
  41         );
  42
  43         pcorrEqn.setReference(pRefCell, pRefValue);
  44         pcorrEqn.solve(mesh.solver(pd.name()));
  45
  46         if (nonOrth == nNonOrthCorr)
  47         {
  48             phi -= pcorrEqn.flux();
  49         }
  50
  51         // Fluxes are corrected to absolute velocity and further corrected
  52         // later.  HJ, 6/Feb/2009
  53     }
  54
  55 #   include "continuityErrs.H"
  56 }