applications/solvers/incompressible/channelFoam/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("PISO"), pRefCell, pRefValue);
  36
  37
  38     singlePhaseTransportModel laminarTransport(U, phi);
  39
  40     autoPtr<incompressible::LESModel> sgsModel
  41     (
  42         incompressible::LESModel::New(U, phi, laminarTransport)
  43     );