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& T = const_cast<volScalarField&>(thermo.T());
  10     volScalarField& p = thermo.p();
  11     volScalarField& e = thermo.e();
  12
  13     volScalarField psis("psi", thermo.psi()/thermo.Cp()*thermo.Cv());
  14     psis.oldTime();
  15
  16     volScalarField rho
  17     (
  18         IOobject
  19         (
  20             "rho",
  21             runTime.timeName(),
  22             mesh
  23         ),
  24         thermo.rho()
  25     );
  26
  27     Info<< "Reading 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     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     );