createFields.H
Go to the documentation of this file.
1Info<< "Reading field p_rgh\n" << endl;
3(
4 IOobject
5 (
6 "p_rgh",
7 runTime.timeName(),
8 mesh,
9 IOobject::MUST_READ,
10 IOobject::AUTO_WRITE
11 ),
12 mesh
13);
14
15Info<< "Reading field U\n" << endl;
17(
18 IOobject
19 (
20 "U",
21 runTime.timeName(),
22 mesh,
23 IOobject::MUST_READ,
24 IOobject::AUTO_WRITE
25 ),
26 mesh
27);
28
29#include "createPhi.H"
30
31multiphaseMixture mixture(U, phi);
32
33// Need to store rho for ddt(rho, U)
35(
36 IOobject
37 (
38 "rho",
39 runTime.timeName(),
40 mesh,
41 IOobject::READ_IF_PRESENT
42 ),
43 mixture.rho()
44);
45rho.oldTime();
46
47// Construct incompressible turbulence model
48autoPtr<incompressible::turbulenceModel> turbulence
49(
50 incompressible::turbulenceModel::New(U, phi, mixture)
51);
52
53
54#include "readGravitationalAcceleration.H"
55#include "readhRef.H"
56#include "gh.H"
57
58
59volScalarField p
60(
61 IOobject
62 (
63 "p",
64 runTime.timeName(),
65 mesh,
66 IOobject::NO_READ,
67 IOobject::AUTO_WRITE
68 ),
69 p_rgh + rho*gh
70);
71
72label pRefCell = 0;
73scalar pRefValue = 0.0;
75(
76 p,
77 p_rgh,
78 pimple.dict(),
81);
82
83if (p_rgh.needReference())
84{
85 p += dimensionedScalar
86 (
87 "p",
88 p.dimensions(),
89 pRefValue - getRefCellValue(p, pRefCell)
90 );
91}
92
93mesh.setFluxRequired(p_rgh.name());
94
95#include "createMRF.H"
96#include "createFvOptions.H"
volScalarField & p_rgh
surfaceScalarField & phi
const scalar pRefValue
const label pRefCell
const volScalarField & gh
pimpleControl & pimple
U
Definition: pEqn.H:72
volScalarField & p
dynamicFvMesh & mesh
engineTime & runTime
compressible::turbulenceModel & turbulence
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
setRefCell(p, pimple.dict(), pRefCell, pRefValue)
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:39