createFields.H

Go to the documentation of this file.

Info<< "Reading field p\n" << endl;
volScalarField p
(
 IOobject
 (
  "p",
  runTime.timeName(),
  mesh,
  IOobject::MUST_READ,
  IOobject::AUTO_WRITE
 ),
 mesh
 );
 
Info<< "Reading field U\n" << endl;
volVectorField U
(
 IOobject
 (
  "U",
  runTime.timeName(),
  mesh,
  IOobject::MUST_READ,
  IOobject::AUTO_WRITE
 ),
 mesh
 );
 
 
 
 
Info<< "Reading field airAge\n" << endl;
 
volScalarField airAge
(
 IOobject
 (
  "airAge",
  runTime.timeName(),
  mesh,
  IOobject::MUST_READ,
  IOobject::AUTO_WRITE
 ),
 mesh
 );
 
 
Info<< "Reading transportProperties\n" << endl;
 
IOdictionary transportProperties
(
 IOobject
 (
  "transportProperties",
  runTime.constant(),
  mesh,
  IOobject::MUST_READ_IF_MODIFIED,
  IOobject::NO_WRITE
 )
 );
 
 
Info<< "Reading diffusivity DT\n" << endl;
 
dimensionedScalar DT
(
 "DT",
 dimArea/dimTime,
 transportProperties
 );
 
 
#include "createPhi.H"
//#include "createFvOptions.H"
 
 
label pRefCell = 0;
scalar pRefValue = 0.0;
setRefCell(p, simple.dict(), pRefCell, pRefValue);
mesh.setFluxRequired(p.name());
 
 
singlePhaseTransportModel laminarTransport(U, phi);
 
    autoPtr<incompressible::turbulenceModel> turbulence
    (
     incompressible::turbulenceModel::New(U, phi, laminarTransport)
    );