applications/solvers/engine/turbDyMEngineFoam/createFields.H

   1
   2     Info<< "Reading field p\n" << endl;
   3     volScalarField p
   4     (
   5         IOobject
   6         (
   7             "p",
   8             runTime.timeName(),
   9             mesh,
  10             IOobject::MUST_READ,
  11             IOobject::AUTO_WRITE
  12         ),
  13         mesh
  14     );
  15
  16
  17     Info<< "Reading field U\n" << endl;
  18     volVectorField U
  19     (
  20         IOobject
  21         (
  22             "U",
  23             runTime.timeName(),
  24             mesh,
  25             IOobject::MUST_READ,
  26             IOobject::AUTO_WRITE
  27         ),
  28         mesh
  29     );
  30
  31 #   include "createPhi.H"
  32
  33
  34     label pRefCell = 0;
  35     scalar pRefValue = 0.0;
  36     setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);
  37
  38     scalar totalVolume = sum(mesh.V()).value();
  39
  40     singlePhaseTransportModel laminarTransport(U, phi);
  41
  42     autoPtr<incompressible::turbulenceModel> turbulence
  43     (
  44         incompressible::turbulenceModel::New(U, phi, laminarTransport)
  45     );
  46
  47     Info<< "Reading field rUA if present\n" << endl;
  48     volScalarField rUA
  49     (
  50         IOobject
  51         (
  52             "rUA",
  53             runTime.timeName(),
  54             mesh,
  55             IOobject::READ_IF_PRESENT,
  56             IOobject::AUTO_WRITE
  57         ),
  58         mesh,
  59         runTime.deltaT(),
  60         zeroGradientFvPatchScalarField::typeName
  61     );