mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-11-28 03:28:01 +00:00
365f9b0006
to simplify construction of dimensionedScalar properties and avoid the duplication of the name string in the constructor call.
108 lines
2.1 KiB
C
108 lines
2.1 KiB
C
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
|
|
);
|
|
|
|
word kinematicCloudName("kinematicCloud");
|
|
args.optionReadIfPresent("cloudName", kinematicCloudName);
|
|
|
|
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.headerOk())
|
|
{
|
|
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.headerOk())
|
|
{
|
|
Info<< "Reading field HdotGradH" << endl;
|
|
|
|
HdotGradHPtr.reset(new volVectorField(HdotGradHheader, mesh));
|
|
}
|
|
|
|
#include "createNonInertialFrameFields.H"
|