applications/solvers/incompressible/pimpleFoam/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 Info<< "Reading field U\n" << endl;
  16 volVectorField U
  17 (
  18     IOobject
  19     (
  20         "U",
  21         runTime.timeName(),
  22         mesh,
  23         IOobject::MUST_READ,
  24         IOobject::AUTO_WRITE
  25     ),
  26     mesh
  27 );
  28
  29 #include <finiteVolume/createPhi.H>
  30
  31
  32 label pRefCell = 0;
  33 scalar pRefValue = 0.0;
  34 setRefCell(p, mesh.solutionDict().subDict("PIMPLE"), pRefCell, pRefValue);
  35
  36
  37 singlePhaseTransportModel laminarTransport(U, phi);
  38
  39 autoPtr<incompressible::turbulenceModel> turbulence
  40 (
  41     incompressible::turbulenceModel::New(U, phi, laminarTransport)
  42 );