src/navalFoam/pdEqn.H

   1 {
   2     if (nOuterCorr != 1)
   3     {
   4         pd.storePrevIter();
   5     }
   6
   7     rAU = 1.0/UEqn.A();
   8     surfaceScalarField rAUf = fvc::interpolate(rAU);
   9
  10     U = rAU*UEqn.H();
  11     U.correctBoundaryConditions();
  12
  13     surfaceScalarField phiU
  14     (
  15         "phiU",
  16         (fvc::interpolate(U) & mesh.Sf())
  17         // Restart accuracy lost with ddtCorr.  Removed, HJ, 16/Nov/2010
  18 //       + fvc::ddtPhiCorr(rAU, rho, U, phi)
  19     );
  20
  21     phi = phiU +
  22         (
  23             fvc::interpolate(interface.sigmaK())*fvc::snGrad(alpha1)
  24           - ghf*fvc::snGrad(rho)
  25         )*rAUf*mesh.magSf();
  26
  27     if (!mesh.moving())
  28     {
  29         adjustPhi(phi, U, pd);
  30     }
  31
  32     for(int nonOrth = 0; nonOrth <= nNonOrthCorr; nonOrth++)
  33     {
  34         fvScalarMatrix pdEqn
  35         (
  36             fvm::laplacian(rAUf, pd) == fvc::div(phi)
  37         );
  38
  39         pdEqn.setReference(pRefCell, pRefValue);
  40
  41         // Solve pressure equation to a tight tolerance in the last corrector
  42         if
  43         (
  44             oCorr == nOuterCorr - 1
  45          && corr == nCorr - 1
  46          && nonOrth == nNonOrthCorr
  47         )
  48         {
  49             pdEqn.solve(mesh.solver("pdFinal"));
  50         }
  51         else
  52         {
  53             pdEqn.solve();
  54         }
  55
  56         if (nonOrth == nNonOrthCorr)
  57         {
  58             phi -= pdEqn.flux();
  59         }
  60     }
  61
  62     // Explicitly relax pressure except for last corrector
  63     if (oCorr != nOuterCorr - 1)
  64     {
  65         pd.relax();
  66     }
  67
  68     U += rAU*fvc::reconstruct((phi - phiU)/rAUf);
  69     U.correctBoundaryConditions();
  70 }