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