createFields.H
Go to the documentation of this file.
1#include "readGravitationalAcceleration.H"
2
3Info<< "Reading thermophysical properties\n" << endl;
4
5autoPtr<psiThermo> pThermo
6(
7 psiThermo::New(mesh)
8);
9psiThermo& thermo = pThermo();
10
11volScalarField rho
12(
13 IOobject
14 (
15 "rho",
16 runTime.timeName(),
17 mesh,
18 IOobject::NO_READ,
19 IOobject::AUTO_WRITE
20 ),
21 thermo.rho()
22);
23
24Info<< "\nReading field U\n" << endl;
25volVectorField U
26(
27 IOobject
28 (
29 "U",
30 runTime.timeName(),
31 mesh,
32 IOobject::MUST_READ,
33 IOobject::AUTO_WRITE
34 ),
35 mesh
36);
37
39
40Info<< "Creating turbulence model\n" << endl;
41autoPtr<compressible::turbulenceModel> turbulence
42(
43 compressible::turbulenceModel::New
44 (
45 rho,
46 U,
47 phi,
48 thermo
49 )
50);
51
53(
54 args.getOrDefault<word>("cloud", "kinematicCloud")
55);
56
57Info<< "Constructing kinematicCloud " << kinematicCloudName << endl;
58basicKinematicCloud kinematicCloud
59(
61 rho,
62 U,
63 thermo.mu(),
64 g
65);
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 parcels.
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
U
Definition: pEqn.H:72
Creates and initialises the face-flux field phi.
dynamicFvMesh & mesh
engineTime & runTime
compressible::turbulenceModel & turbulence
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)
Info<< "Reading thermophysical properties\n"<< endl;autoPtr< psiReactionThermo > pThermo(psiReactionThermo::New(mesh))
const word kinematicCloudName(args.getOrDefault< word >("cloud", "kinematicCloud"))