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 p_rgh\n" << endl;
 
//volScalarField p_rgh
//(
// IOobject
// (
//  "p_rgh",
//  runTime.timeName(),
//  mesh,
//  IOobject::MUST_READ,
//  IOobject::AUTO_WRITE
// ),
// mesh
// );
 
volScalarField p_rgh 
(
 IOobject
 (
  "p_rgh",
  runTime.timeName(),
  mesh,
  IOobject::READ_IF_PRESENT,
  IOobject::AUTO_WRITE
 ),
 mesh,
 dimensionedScalar("p_rgh", dimensionSet(0,2,-2,0,0,0,0), 101325)
 );
 
 
Info<< "Reading field T\n" << endl;
volScalarField T
(
 IOobject
 (
  "T",
  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
 )
 );
 
 
// kinematic turbulent thermal thermal conductivity m2/s
Info<< "Reading field alphat\n" << endl;
volScalarField alphat
(
         IOobject
        (
            "alphat",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh,
        dimensionedScalar("alphat", dimLength*dimLength/dimTime,0.1)
);
 
 
Info<< "Reading diffusivity DT\n" << endl;
 
dimensionedScalar DT
(
 "DT",
 dimArea/dimTime,
 transportProperties
 );
 
     bool solveScalarT =
     transportProperties.lookupOrDefault("solveScalarT",true);
 
     bool solveComfort =
     transportProperties.lookupOrDefault("solveComfort",true);
 
//Info<< "Solving Temperature : " << solveScalarT <<  endl;
//Info<< "Solving PMV and PPD : " << solveComfort <<  endl;
 
#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)
    );