applications/solvers/solidMechanics/icoFsiElasticNonLinULSolidFoam/solveFluid.H

   1
   2         if(mesh.moving())
   3         {
   4             // Make the fluxes relative
   5             phi -= fvc::meshPhi(U);
   6         }
   7
   8 #       include "CourantNo.H"
   9
  10         fvVectorMatrix UEqn
  11         (
  12             fvm::ddt(U)
  13           + fvm::div(phi, U)
  14           - fvm::laplacian(nu, U)
  15         );
  16
  17         solve(UEqn == -fvc::grad(p));
  18
  19         // --- PISO loop
  20         volScalarField rUA = 1.0/UEqn.A();
  21
  22         for (int corr=0; corr<nCorr; corr++)
  23         {
  24             U = rUA*UEqn.H();
  25             phi = (fvc::interpolate(U) & mesh.Sf());
  26
  27             adjustPhi(phi, U, p);
  28
  29             for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
  30             {
  31                 fvScalarMatrix pEqn
  32                 (
  33                     fvm::laplacian(rUA, p)
  34                  == fvc::div(phi)
  35                 );
  36
  37                 pEqn.setReference(pRefCell, pRefValue);
  38                 pEqn.solve();
  39
  40                 if (nonOrth == nNonOrthCorr)
  41                 {
  42                     phi -= pEqn.flux();
  43                 }
  44             }
  45
  46 #           include "continuityErrs.H"
  47
  48             U -= rUA*fvc::grad(p);
  49             U.correctBoundaryConditions();
  50         }