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