createFields.H
Go to the documentation of this file.
1#include "readGravitationalAcceleration.H"
2#include "readhRef.H"
3
4Info<< "Creating twoPhaseSystem\n" << endl;
5
6twoPhaseSystem fluid(mesh, g);
7
8phaseModel& phase1 = fluid.phase1();
9phaseModel& phase2 = fluid.phase2();
10
14
15volScalarField& p = phase1.thermo().p();
16
18(
19 "pMin",
20 dimPressure,
21 fluid
22);
23
24#include "gh.H"
25
26Info<< "Reading field p_rgh\n" << endl;
28(
29 IOobject
30 (
31 "p_rgh",
32 runTime.timeName(),
33 mesh,
34 IOobject::MUST_READ,
35 IOobject::AUTO_WRITE
36 ),
37 mesh
38);
39
41(
42 IOobject
43 (
44 "U",
45 runTime.timeName(),
46 mesh,
47 IOobject::NO_READ,
48 IOobject::AUTO_WRITE
49 ),
50 fluid.U()
51);
52
53label pRefCell = 0;
54scalar pRefValue = 0.0;
56(
57 p,
58 p_rgh,
59 pimple.dict(),
62);
63mesh.setFluxRequired(p_rgh.name());
64mesh.setFluxRequired(alpha1.name());
65
66Info<< "Creating field dpdt\n" << endl;
68(
69 IOobject
70 (
71 "dpdt",
72 runTime.timeName(),
73 mesh
74 ),
75 mesh,
76 dimensionedScalar(p.dimensions()/dimTime, Zero)
77);
78
79
80Info<< "Creating field kinetic energy K\n" << endl;
81volScalarField K1(IOobject::groupName("K", phase1.name()), 0.5*magSqr(U1));
82volScalarField K2(IOobject::groupName("K", phase2.name()), 0.5*magSqr(U2));
83
84#include "createMRF.H"
85#include "createFvOptions.H"
#define K1
Definition: SHA1.C:144
#define K2
Definition: SHA1.C:145
const uniformDimensionedVectorField & g
volScalarField & p_rgh
const scalar pRefValue
const label pRefCell
pimpleControl & pimple
const volScalarField & alpha1
twoPhaseSystem & fluid
phaseModel & phase1
phaseModel & phase2
volVectorField & U1
const dimensionedScalar & pMin
volVectorField & U2
U
Definition: pEqn.H:72
volScalarField & p
dynamicFvMesh & mesh
engineTime & runTime
volScalarField & dpdt
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)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:372
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
setRefCell(p, pimple.dict(), pRefCell, pRefValue)