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
36
39
40Info<< "Reading thermophysical properties\n" << endl;
41
42const volScalarField& rho1 = mixture.thermo1().rho();
43const volScalarField& rho2 = mixture.thermo2().rho();
44
46(
47 IOobject
48 (
49 "rho",
50 runTime.timeName(),
51 mesh,
52 IOobject::READ_IF_PRESENT,
53 IOobject::AUTO_WRITE
54 ),
56);
57
58
60(
61 "pMin",
62 dimPressure,
64);
65
66mesh.setFluxRequired(p_rgh.name());
67mesh.setFluxRequired(alpha1.name());
68
69
70#include "readGravitationalAcceleration.H"
71#include "readhRef.H"
72#include "gh.H"
73
74
75// Mass flux
76// Initialisation does not matter because rhoPhi is reset after the
77// alpha1 solution before it is used in the U equation.
79(
80 IOobject
81 (
82 "rhoPhi",
83 runTime.timeName(),
84 mesh,
85 IOobject::NO_READ,
86 IOobject::NO_WRITE
87 ),
88 fvc::interpolate(rho)*phi
89);
90
91volScalarField dgdt(alpha1*fvc::div(phi));
92
93#include "createAlphaFluxes.H"
94
96// Construct compressible turbulence model
97compressibleInterPhaseTransportModel turbulence
98(
99 rho,
100 U,
101 phi,
102 rhoPhi,
104 mixture
105);
106
107#include "createK.H"
108
109#include "createMRF.H"
110#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
An implementation of the isoAdvector geometric Volume-of-Fluid method advancing the provided volume f...
Definition: isoAdvection.H:92
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
isoAdvection advector(alpha1, phi, U)