createFields.H
Go to the documentation of this file.
1#include "readGravitationalAcceleration.H"
2
3Info<< "Reading field U\n" << endl;
5(
6 IOobject
7 (
8 "U",
9 runTime.timeName(),
10 mesh,
11 IOobject::MUST_READ,
12 IOobject::AUTO_WRITE
13 ),
14 mesh
15);
16
17#include "createPhi.H"
18
19singlePhaseTransportModel laminarTransport(U, phi);
20
21dimensionedScalar rhoInfValue
22(
23 "rhoInf",
24 dimDensity,
26);
27
28volScalarField rhoInf
29(
30 IOobject
31 (
32 "rho",
33 runTime.timeName(),
34 mesh,
35 IOobject::NO_READ,
36 IOobject::AUTO_WRITE
37 ),
38 mesh,
39 rhoInfValue
40);
41
42autoPtr<incompressible::turbulenceModel> turbulence
43(
44 incompressible::turbulenceModel::New(U, phi, laminarTransport)
45);
46
47volScalarField mu
48(
49 IOobject
50 (
51 "mu",
52 runTime.timeName(),
53 mesh,
54 IOobject::NO_READ,
55 IOobject::AUTO_WRITE
56 ),
57 laminarTransport.nu()*rhoInfValue
58);
59
61(
62 args.getOrDefault<word>("cloud", "kinematicCloud")
63);
64
65Info<< "Constructing kinematicCloud " << kinematicCloudName << endl;
67(
69 rhoInf,
70 U,
71 mu,
72 g
73);
74
75IOobject Hheader
76(
77 "H",
78 runTime.timeName(),
79 mesh,
80 IOobject::MUST_READ,
81 IOobject::AUTO_WRITE
82);
83
84autoPtr<volVectorField> HPtr;
85
86if (Hheader.typeHeaderOk<volVectorField>(true))
87{
88 Info<< "\nReading field H\n" << endl;
89
90 HPtr.reset(new volVectorField (Hheader, mesh));
91}
92
94(
95 "HdotGradH",
96 runTime.timeName(),
97 mesh,
98 IOobject::MUST_READ,
99 IOobject::AUTO_WRITE
100);
101
102autoPtr<volVectorField> HdotGradHPtr;
103
104if (HdotGradHheader.typeHeaderOk<volVectorField>(true))
105{
106 Info<< "Reading field HdotGradH" << endl;
107
108 HdotGradHPtr.reset(new volVectorField(HdotGradHheader, mesh));
109}
110
const uniformDimensionedVectorField & g
surfaceScalarField & phi
T getOrDefault(const word &optName, const T &deflt) const
Get a value from the named option if present, or return default.
Definition: argListI.H:307
Cloud class to introduce kinematic colliding parcels.
U
Definition: pEqn.H:72
const volScalarField & mu
dynamicFvMesh & mesh
engineTime & runTime
compressible::turbulenceModel & turbulence
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:83
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:82
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:372
Foam::argList args(argc, argv)
singlePhaseTransportModel laminarTransport(U, phi)
autoPtr< volVectorField > HdotGradHPtr
Definition: createFields.H:102
IOobject HdotGradHheader("HdotGradH", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE)
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
Definition: createFields.H:84
const word kinematicCloudName(args.getOrDefault< word >("cloud", "kinematicCloud"))