createFields.H
Go to the documentation of this file.
1 Info<< "Reading thermophysical properties\n" << endl;
2 
3 autoPtr<rhoThermo> pThermo(rhoThermo::New(mesh));
4 rhoThermo& thermo = pThermo();
5 thermo.validate(args.executable(), "h", "e");
6 
8 (
9  IOobject
10  (
11  "rho",
12  runTime.timeName(),
13  mesh,
14  IOobject::NO_READ,
15  IOobject::NO_WRITE
16  ),
17  thermo.rho()
18 );
19 
20 volScalarField& p = thermo.p();
21 const volScalarField& psi = thermo.psi();
22 
23 
24 Info<< "Reading field U\n" << endl;
26 (
27  IOobject
28  (
29  "U",
30  runTime.timeName(),
31  mesh,
32  IOobject::MUST_READ,
33  IOobject::AUTO_WRITE
34  ),
35  mesh
36 );
37 
38 #include "compressibleCreatePhi.H"
39 
40 
41 Info<< "Creating turbulence model\n" << endl;
42 autoPtr<compressible::turbulenceModel> turbulence
43 (
45  (
46  rho,
47  U,
48  phi,
49  thermo
50  )
51 );
52 
53 
54 #include "readGravitationalAcceleration.H"
55 #include "readhRef.H"
56 #include "gh.H"
57 
58 
59 Info<< "Reading field p_rgh\n" << endl;
61 (
62  IOobject
63  (
64  "p_rgh",
65  runTime.timeName(),
66  mesh,
67  IOobject::MUST_READ,
68  IOobject::AUTO_WRITE
69  ),
70  mesh
71 );
72 
73 // Force p_rgh to be consistent with p
74 p_rgh = p - rho*gh;
75 
76 mesh.setFluxRequired(p_rgh.name());
77 
78 Info<< "Creating field dpdt\n" << endl;
80 (
81  IOobject
82  (
83  "dpdt",
84  runTime.timeName(),
85  mesh
86  ),
87  mesh,
88  dimensionedScalar("dpdt", p.dimensions()/dimTime, 0)
89 );
90 
91 Info<< "Creating field kinetic energy K\n" << endl;
92 volScalarField K("K", 0.5*magSqr(U));
volScalarField & dpdt
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
CGAL::Exact_predicates_exact_constructions_kernel K
dimensioned< scalar > magSqr(const dimensioned< Type > &)
Creates and initialises the face-flux field phi.
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
Definition: dimensionSets.H:51
volScalarField & p_rgh
messageStream Info
autoPtr< compressible::turbulenceModel > turbulence
Definition: createFields.H:23
dynamicFvMesh & mesh
autoPtr< BasicCompressibleTurbulenceModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const typename BasicCompressibleTurbulenceModel::transportModel &transport, const word &propertiesName)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
rho
Definition: createFields.H:79
Info<< "Reading thermophysical properties\n"<< endl;autoPtr< psiThermo > pThermo(psiThermo::New(mesh))
const word & executable() const
Name of executable without the path.
Definition: argListI.H:30
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:55
Foam::argList args(argc, argv)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const volScalarField & gh