createFields.H
Go to the documentation of this file.
1Info<< "Reading field p (kinematic)\n" << endl;
3(
4 IOobject
5 (
6 "p",
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
31singlePhaseTransportModel laminarTransport(U, phi);
32
33autoPtr<incompressible::turbulenceModel> turbulence
34(
35 incompressible::turbulenceModel::New(U, phi, laminarTransport)
36);
37
38#include "readGravitationalAcceleration.H"
39
40Info<< "Creating field zeta\n" << endl;
41volVectorField zeta
42(
43 IOobject
44 (
45 "zeta",
46 runTime.timeName(),
47 mesh,
48 IOobject::READ_IF_PRESENT,
49 IOobject::AUTO_WRITE
50 ),
51 mesh,
52 dimensionedVector(dimLength, Zero)
53);
54
55Info<< "Creating field p_gh\n" << endl;
56volScalarField p_gh
57(
58 IOobject
59 (
60 "p_gh",
61 runTime.timeName(),
62 mesh,
63 IOobject::MUST_READ,
64 IOobject::AUTO_WRITE
65 ),
66 mesh
67);
68
69// Force p_gh to be consistent with p
70// Height is made relative to field 'refLevel'
71p_gh = p - (g & mesh.C());
72
73
74label p_ghRefCell = 0;
75scalar p_ghRefValue = 0.0;
76setRefCell(p_gh, pimple.dict(), p_ghRefCell, p_ghRefValue);
77mesh.setFluxRequired(p_gh.name());
78
79#include "createMRF.H"
80#include "createFvOptions.H"
const uniformDimensionedVectorField & g
surfaceScalarField & phi
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
singlePhaseTransportModel laminarTransport(U, phi)
setRefCell(p, pimple.dict(), pRefCell, pRefValue)