applications/solvers/compressible/rhoPisoFoam/createFields.H

   1     Info<< "Reading thermophysical properties\n" << endl;
   2
   3     autoPtr<basicPsiThermo> pThermo
   4     (
   5         basicPsiThermo::New(mesh)
   6     );
   7     basicPsiThermo& thermo = pThermo();
   8
   9     volScalarField& p = thermo.p();
  10     volScalarField& h = thermo.h();
  11     const volScalarField& psi = thermo.psi();
  12
  13     volScalarField rho
  14     (
  15         IOobject
  16         (
  17             "rho",
  18             runTime.timeName(),
  19             mesh,
  20             IOobject::NO_READ,
  21             IOobject::AUTO_WRITE
  22         ),
  23         thermo.rho()
  24     );
  25
  26     Info<< "\nReading field U\n" << endl;
  27     volVectorField U
  28     (
  29         IOobject
  30         (
  31             "U",
  32             runTime.timeName(),
  33             mesh,
  34             IOobject::MUST_READ,
  35             IOobject::AUTO_WRITE
  36         ),
  37         mesh
  38     );
  39
  40 #   include "compressibleCreatePhi.H"
  41
  42
  43     Info<< "Creating turbulence model\n" << endl;
  44     autoPtr<compressible::turbulenceModel> turbulence
  45     (
  46         compressible::turbulenceModel::New
  47         (
  48             rho,
  49             U,
  50             phi,
  51             thermo
  52         )
  53     );
  54
  55
  56     Info<< "Creating field DpDt\n" << endl;
  57     volScalarField DpDt =
  58         fvc::DDt(surfaceScalarField("phiU", phi/fvc::interpolate(rho)), p);