createFields.H
Go to the documentation of this file.
1#include "createRDeltaT.H"
2
3Info<< "Reading thermophysical properties\n" << endl;
4
5autoPtr<psiThermo> pThermo
6(
7 psiThermo::New(mesh)
8);
9psiThermo& thermo = pThermo();
10
11volScalarField& e = thermo.he();
12
13Info<< "Reading field U\n" << endl;
14volVectorField U
15(
16 IOobject
17 (
18 "U",
19 runTime.timeName(),
20 mesh,
21 IOobject::MUST_READ,
22 IOobject::AUTO_WRITE
23 ),
24 mesh
25);
26
27volScalarField rho
28(
29 IOobject
30 (
31 "rho",
32 runTime.timeName(),
33 mesh,
34 IOobject::NO_READ,
35 IOobject::AUTO_WRITE
36 ),
37 thermo.rho()
38);
39
40volVectorField rhoU
41(
42 IOobject
43 (
44 "rhoU",
45 runTime.timeName(),
46 mesh,
47 IOobject::NO_READ,
48 IOobject::NO_WRITE
49 ),
50 rho*U
51);
52
53volScalarField rhoE
54(
55 IOobject
56 (
57 "rhoE",
58 runTime.timeName(),
59 mesh,
60 IOobject::NO_READ,
61 IOobject::NO_WRITE
62 ),
63 rho*(e + 0.5*magSqr(U))
64);
65
66surfaceScalarField pos
67(
68 IOobject
69 (
70 "pos",
71 runTime.timeName(),
72 mesh
73 ),
74 mesh,
75 dimensionedScalar("pos", dimless, 1.0)
76);
77
78surfaceScalarField neg
79(
80 IOobject
81 (
82 "neg",
83 runTime.timeName(),
84 mesh
85 ),
86 mesh,
87 dimensionedScalar("neg", dimless, -1.0)
88);
89
90surfaceScalarField phi("phi", fvc::flux(rhoU));
91
92Info<< "Creating turbulence model\n" << endl;
93autoPtr<compressible::turbulenceModel> turbulence
94(
95 compressible::turbulenceModel::New
96 (
97 rho,
98 U,
99 phi,
100 thermo
101 )
102);
surfaceScalarField & phi
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
U
Definition: pEqn.H:72
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
volScalarField & rhoE
Info<< "Reading thermophysical properties\n"<< endl;autoPtr< psiReactionThermo > pThermo(psiReactionThermo::New(mesh))
volScalarField & e
Definition: createFields.H:11