applications/solvers/combustion/coldEngineFoam/createFields.H

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