createFields.H
Go to the documentation of this file.
1  Info<< "Reading field U\n" << endl;
3  (
4  IOobject
5  (
6  "U",
7  runTime.timeName(),
8  mesh,
9  IOobject::MUST_READ,
10  IOobject::AUTO_WRITE
11  ),
12  mesh
13  );
14 
15  #include "createPhi.H"
16 
17  singlePhaseTransportModel laminarTransport(U, phi);
18 
19  dimensionedScalar rhoInfValue
20  (
21  "rhoInf",
22  dimDensity,
24  );
25 
26  volScalarField rhoInf
27  (
28  IOobject
29  (
30  "rho",
31  runTime.timeName(),
32  mesh,
33  IOobject::NO_READ,
34  IOobject::AUTO_WRITE
35  ),
36  mesh,
37  rhoInfValue
38  );
39 
40  autoPtr<incompressible::turbulenceModel> turbulence
41  (
43  );
44 
46  (
47  IOobject
48  (
49  "mu",
50  runTime.timeName(),
51  mesh,
52  IOobject::NO_READ,
53  IOobject::AUTO_WRITE
54  ),
55  laminarTransport.nu()*rhoInfValue
56  );
57 
58  word kinematicCloudName("kinematicCloud");
60 
61  Info<< "Constructing kinematicCloud " << kinematicCloudName << endl;
62  basicKinematicCollidingCloud kinematicCloud
63  (
65  rhoInf,
66  U,
67  mu,
68  g
69  );
70 
71  IOobject Hheader
72  (
73  "H",
74  runTime.timeName(),
75  mesh,
76  IOobject::MUST_READ,
77  IOobject::AUTO_WRITE
78  );
79 
80  autoPtr<volVectorField> HPtr;
81 
82  if (Hheader.headerOk())
83  {
84  Info<< "\nReading field H\n" << endl;
85 
86  HPtr.reset(new volVectorField (Hheader, mesh));
87  }
88 
89  IOobject HdotGradHheader
90  (
91  "HdotGradH",
92  runTime.timeName(),
93  mesh,
94  IOobject::MUST_READ,
95  IOobject::AUTO_WRITE
96  );
97 
98  autoPtr<volVectorField> HdotGradHPtr;
99 
100  if (HdotGradHheader.headerOk())
101  {
102  Info<< "Reading field HdotGradH" << endl;
103 
104  HdotGradHPtr.reset(new volVectorField(HdotGradHheader, mesh));
105  }
106 
autoPtr< volVectorField > HdotGradHPtr
Definition: createFields.H:98
volVectorField U(IOobject( "U", runTime.timeName(), runTime, IOobject::NO_READ, IOobject::NO_WRITE ), mesh, dimensionedVector("zero", dimVelocity, vector::zero), p.boundaryField().types())
surfaceScalarField & phi
Definition: createFields.H:13
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:80
CollidingCloud< KinematicCloud< Cloud< basicKinematicCollidingParcel > > > basicKinematicCollidingCloud
singlePhaseTransportModel laminarTransport(U, phi)
messageStream Info
autoPtr< compressible::turbulenceModel > turbulence
Definition: createFields.H:23
dynamicFvMesh & mesh
const dimensionSet dimDensity
autoPtr< BasicCompressibleTurbulenceModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const typename BasicCompressibleTurbulenceModel::transportModel &transport, const word &propertiesName)
IOobject HdotGradHheader("HdotGradH", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE)
const volScalarField & mu
Atomic mass unit.
Definition: createFields.H:13
word kinematicCloudName("kinematicCloud")
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
bool optionReadIfPresent(const word &opt, T &) const
Read a value from the named option if present.
Definition: argListI.H:198
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:55
Foam::argList args(argc, argv)
const dimensionedVector & g
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.