createFields.H
Go to the documentation of this file.
1 #include "createRDeltaT.H"
2 
3 Info<< "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 
17 Info<< "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 
33 
34 Info<< "Reading transportProperties\n" << endl;
35 immiscibleIncompressibleTwoPhaseMixture mixture(U, phi);
36 
37 volScalarField& alpha1(mixture.alpha1());
38 volScalarField& alpha2(mixture.alpha2());
39 
40 const dimensionedScalar& rho1 = mixture.rho1();
41 const dimensionedScalar& rho2 = mixture.rho2();
42 
43 
44 // Need to store rho for ddt(rho, U)
46 (
47  IOobject
48  (
49  "rho",
50  runTime.timeName(),
51  mesh,
52  IOobject::READ_IF_PRESENT
53  ),
54  alpha1*rho1 + alpha2*rho2
55 );
56 rho.oldTime();
57 
58 
59 // Mass flux
61 (
62  IOobject
63  (
64  "rhoPhi",
65  runTime.timeName(),
66  mesh,
67  IOobject::NO_READ,
68  IOobject::NO_WRITE
69  ),
71 );
72 
73 
74 // Construct incompressible turbulence model
75 autoPtr<incompressible::turbulenceModel> turbulence
76 (
78 );
79 
80 
81 #include "readGravitationalAcceleration.H"
82 #include "readhRef.H"
83 #include "gh.H"
84 
85 
87 (
88  IOobject
89  (
90  "p",
91  runTime.timeName(),
92  mesh,
93  IOobject::NO_READ,
94  IOobject::AUTO_WRITE
95  ),
96  p_rgh + rho*gh
97 );
98 
99 label pRefCell = 0;
100 scalar pRefValue = 0.0;
102 (
103  p,
104  p_rgh,
105  pimple.dict(),
106  pRefCell,
107  pRefValue
108 );
109 
110 if (p_rgh.needReference())
111 {
113  (
114  "p",
115  p.dimensions(),
116  pRefValue - getRefCellValue(p, pRefCell)
117  );
118  p_rgh = p - rho*gh;
119 }
120 
121 mesh.setFluxRequired(p_rgh.name());
122 mesh.setFluxRequired(alpha1.name());
123 
124 #include "createMRF.H"
125 #include "createFvOptions.H"
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
pimpleNoLoopControl & pimple
surfaceScalarField & phi
volScalarField p_rgh(IOobject("p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh)
setRefCell(p, p_rgh, pimple.dict(), pRefCell, pRefValue)
engineTime & runTime
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:256
volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mesh, dimensionedVector("zero", dimVelocity, Zero))
volScalarField & alpha2(mixture->alpha2())
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:55
scalar getRefCellValue(const volScalarField &field, const label refCelli)
Return the current value of field in the reference cell.
Definition: findRefCell.C:136
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
rho
Definition: createFields.H:81
dynamicFvMesh & mesh
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar("pos", dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar("neg", dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
Info<< "Reading field p_rgh\"<< endl;volScalarField p_rgh(IOobject("p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Creating phaseChangeTwoPhaseMixture\"<< endl;autoPtr< phaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:33
volScalarField & p
Definition: createFields.H:23
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
label pRefCell
Definition: createFields.H:106
rhoPhi
Definition: rhoEqn.H:10
volScalarField & alpha1(mixture->alpha1())
const volScalarField & gh
messageStream Info
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
scalar pRefValue
Definition: createFields.H:107