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
);
 
Info<< "Reading field U\n" << endl;
volVectorField U
(
    IOobject
    (
        "U",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);
 
 
Info<< "Reading field T\n" << endl;
volScalarField T
(
    IOobject
    (
        "T",
        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
    (
        transportProperties.lookup("DT")
    );
 
Info<< "Reading turbulent Schmidt number\n" << endl;
 
    scalar Sct
    (
        readScalar(transportProperties.lookup("SchmidtTurbulent"))
    );
 
    bool solveMomentum =
         transportProperties.lookupOrDefault("SolveMomentum",true);   
 
Info<< "Solving momentum equation: " << solveMomentum <<  endl;
 
 
    bool solveScalar =
         transportProperties.lookupOrDefault("SolveScalar",true);   
 
Info<< "Solving turbulent scalar transport: " << solveScalar <<  endl;
 
 
#include "createPhi.H"
 
 
label pRefCell = 0;
scalar pRefValue = 0.0;
setRefCell(p, piso.dict(), pRefCell, pRefValue);
mesh.setFluxRequired(p.name());
 
 
singlePhaseTransportModel laminarTransport(U, phi);
 
autoPtr<incompressible::turbulenceModel> turbulence
(
    incompressible::turbulenceModel::New(U, phi, laminarTransport)
);