applications/solvers/combustion/fireFoam/UEqn.H

   1 fvVectorMatrix UEqn
   2 (
   3     fvm::ddt(rho, U)
   4   + fvm::div(phi, U)
   5   + turbulence->divDevRhoReff(U)
   6 );
   7
   8 UEqn.relax();
   9
  10 if (oCorr == nOuterCorr - 1)
  11 {
  12     solve
  13     (
  14         UEqn
  15       ==
  16         fvc::reconstruct
  17         (
  18             fvc::interpolate(rho)*(g & mesh.Sf())
  19           - fvc::snGrad(p)*mesh.magSf()
  20         ),
  21         mesh.solver("UFinal")
  22     );
  23 }
  24 else
  25 {
  26     solve
  27     (
  28         UEqn
  29       ==
  30         fvc::reconstruct
  31         (
  32             fvc::interpolate(rho)*(g & mesh.Sf())
  33           - fvc::snGrad(p)*mesh.magSf()
  34         )
  35     );
  36 }