applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/createFields.H

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