applications/solvers/compressible/steadyCompressibleSRFFoam/rhoFromP.H

   1 {
   2     // Calculate density from pressure
   3     rho.storePrevIter();
   4     rho = thermo.rho()();
   5
   6     // Bound rho
   7     volScalarField R = thermo.Cp() - thermo.Cv();
   8
   9     volScalarField rhoMin = pMin/(R*TMax);
  10     volScalarField rhoMax = pMax/(R*TMin);
  11
  12     rho = Foam::min(rho, rhoMax);
  13     rho = Foam::max(rho, rhoMin);
  14     rho.relax();
  15 }