applications/solvers/incompressible/pimpleFoam/SRFPimpleFoam/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 Urel\n" << endl;
  16 volVectorField Urel
  17 (
  18     IOobject
  19     (
  20         "Urel",
  21         runTime.timeName(),
  22         mesh,
  23         IOobject::MUST_READ,
  24         IOobject::AUTO_WRITE
  25     ),
  26     mesh
  27 );
  28
  29 Info<< "Reading/calculating face flux field phi\n" << endl;
  30 surfaceScalarField phi
  31 (
  32     IOobject
  33     (
  34         "phi",
  35         runTime.timeName(),
  36         mesh,
  37         IOobject::READ_IF_PRESENT,
  38         IOobject::AUTO_WRITE
  39     ),
  40     linearInterpolate(Urel) & mesh.Sf()
  41 );
  42
  43 label pRefCell = 0;
  44 scalar pRefValue = 0.0;
  45 setRefCell(p, mesh.solutionDict().subDict("PIMPLE"), pRefCell, pRefValue);
  46
  47 singlePhaseTransportModel laminarTransport(Urel, phi);
  48
  49 autoPtr<incompressible::turbulenceModel> turbulence
  50 (
  51     incompressible::turbulenceModel::New(Urel, phi, laminarTransport)
  52 );
  53
  54 Info<< "Creating SRF model\n" << endl;
  55 autoPtr<SRF::SRFModel> SRF
  56 (
  57     SRF::SRFModel::New(Urel)
  58 );
  59
  60 // Create the absolute velocity
  61 volVectorField U
  62 (
  63     IOobject
  64     (
  65         "U",
  66         runTime.timeName(),
  67         mesh,
  68         IOobject::NO_READ,
  69         IOobject::AUTO_WRITE
  70     ),
  71     Urel + SRF->U()
  72 );