createFields.H
Go to the documentation of this file.
1#include "createRDeltaT.H"
2
3Info<< "Reading field p_rgh\n" << endl;
5(
6 IOobject
7 (
8 "p_rgh",
9 runTime.timeName(),
10 mesh,
11 IOobject::MUST_READ,
12 IOobject::AUTO_WRITE
13 ),
14 mesh
15);
16
17Info<< "Reading field U\n" << endl;
19(
20 IOobject
21 (
22 "U",
23 runTime.timeName(),
24 mesh,
25 IOobject::MUST_READ,
26 IOobject::AUTO_WRITE
27 ),
28 mesh
29);
30
31#include "createPhi.H"
32
33Info<< "Constructing twoPhaseMixtureThermo\n" << endl;
34twoPhaseMixtureThermo mixture(U, phi);
35
38
39Info<< "Reading thermophysical properties\n" << endl;
40
41const volScalarField& rho1 = mixture.thermo1().rho();
42const volScalarField& rho2 = mixture.thermo2().rho();
43
45(
46 IOobject
47 (
48 "rho",
49 runTime.timeName(),
50 mesh,
51 IOobject::READ_IF_PRESENT,
52 IOobject::AUTO_WRITE
53 ),
55);
56
57
59(
60 "pMin",
61 dimPressure,
63);
64
65mesh.setFluxRequired(p_rgh.name());
66mesh.setFluxRequired(alpha1.name());
67
68
69#include "readGravitationalAcceleration.H"
70#include "readhRef.H"
71#include "gh.H"
72
73
74// Mass flux
75// Initialisation does not matter because rhoPhi is reset after the
76// alpha1 solution before it is used in the U equation.
78(
79 IOobject
80 (
81 "rhoPhi",
82 runTime.timeName(),
83 mesh,
84 IOobject::NO_READ,
85 IOobject::NO_WRITE
86 ),
87 fvc::interpolate(rho)*phi
88);
89
90volScalarField dgdt(alpha1*fvc::div(phi));
91
92#include "createAlphaFluxes.H"
93
94// Construct compressible turbulence model
95compressibleInterPhaseTransportModel turbulence
96(
97 rho,
98 U,
99 phi,
100 rhoPhi,
102 mixture
103);
104
105#include "createK.H"
106
107#include "createMRF.H"
108#include "createFvOptions.H"
rhoPhi
Definition: rhoEqn.H:10
volScalarField & p_rgh
surfaceScalarField & phi
const volScalarField & alpha1
volScalarField & rho2
const volScalarField & alpha2
const dimensionedScalar & pMin
volScalarField & rho1
U
Definition: pEqn.H:72
alphaPhi10
Definition: alphaEqn.H:7
dynamicFvMesh & mesh
engineTime & runTime
compressible::turbulenceModel & turbulence
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
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)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:372
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:39