createFields.H

Go to the documentation of this file.

#include "readGravitationalAcceleration.H"
 
Info<< "Reading field U\n" << endl;
volVectorField U
(
    IOobject
    (
        "U",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);
 
#include "createPhi.H"
 
singlePhaseTransportModel laminarTransport(U, phi);
 
dimensionedScalar rhoInfValue
(
    "rhoInf",
    dimDensity,
    laminarTransport
);
 
volScalarField rhoInf
(
    IOobject
    (
        "rho",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    rhoInfValue
);
 
autoPtr<incompressible::turbulenceModel> turbulence
(
    incompressible::turbulenceModel::New(U, phi, laminarTransport)
);
 
volScalarField mu
(
    IOobject
    (
        "mu",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    laminarTransport.nu()*rhoInfValue
);
 
const word kinematicCloudName
(
    args.getOrDefault<word>("cloud", "kinematicCloud")
);
 
Info<< "Constructing kinematicCloud " << kinematicCloudName << endl;
basicKinematicCollidingCloud kinematicCloud
(
    kinematicCloudName,
    rhoInf,
    U,
    mu,
    g
);
 
IOobject Hheader
(
    "H",
    runTime.timeName(),
    mesh,
    IOobject::MUST_READ,
    IOobject::AUTO_WRITE
);
 
autoPtr<volVectorField> HPtr;
 
if (Hheader.typeHeaderOk<volVectorField>(true))
{
    Info<< "\nReading field H\n" << endl;
 
    HPtr.reset(new volVectorField (Hheader, mesh));
}
 
IOobject HdotGradHheader
(
    "HdotGradH",
    runTime.timeName(),
    mesh,
    IOobject::MUST_READ,
    IOobject::AUTO_WRITE
);
 
autoPtr<volVectorField> HdotGradHPtr;
 
if (HdotGradHheader.typeHeaderOk<volVectorField>(true))
{
    Info<< "Reading field HdotGradH" << endl;
 
    HdotGradHPtr.reset(new volVectorField(HdotGradHheader, mesh));
}
 
#include "createNonInertialFrameFields.H"