createFields.H
Go to the documentation of this file.
1Info<< "Reading field p\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
31
32label pRefCell = 0;
33scalar pRefValue = 0.0;
35mesh.setFluxRequired(p.name());
36
37
38Info<< "Reading field pa\n" << endl;
40(
41 IOobject
42 (
43 "pa",
44 runTime.timeName(),
45 mesh,
46 IOobject::MUST_READ,
47 IOobject::AUTO_WRITE
48 ),
49 mesh
50);
51
52Info<< "Reading field Ua\n" << endl;
54(
55 IOobject
56 (
57 "Ua",
58 runTime.timeName(),
59 mesh,
60 IOobject::MUST_READ,
61 IOobject::AUTO_WRITE
62 ),
63 mesh
64);
65
66#include "createPhia.H"
67
68
69label paRefCell = 0;
70scalar paRefValue = 0.0;
72(
73 pa,
74 simple.dict(),
75 paRefCell,
76 paRefValue
77);
78mesh.setFluxRequired(pa.name());
79
80
81singlePhaseTransportModel laminarTransport(U, phi);
82
83autoPtr<incompressible::turbulenceModel> turbulence
84(
85 incompressible::turbulenceModel::New(U, phi, laminarTransport)
86);
87
88
89dimensionedScalar zeroSensitivity(dimVelocity*dimVelocity, Zero);
90dimensionedScalar zeroAlpha(dimless/dimTime, Zero);
91
92dimensionedScalar lambda
93(
94 "lambda",
95 dimTime/sqr(dimLength),
97);
98
99dimensionedScalar alphaMax
100(
101 "alphaMax",
102 dimless/dimTime,
104);
105
106const labelList& inletCells = mesh.boundary()["inlet"].faceCells();
107//const labelList& outletCells = mesh.boundary()["outlet"].faceCells();
108
109volScalarField alpha
110(
111 IOobject
112 (
113 "alpha",
114 runTime.timeName(),
115 mesh,
116 IOobject::READ_IF_PRESENT,
117 IOobject::AUTO_WRITE
118 ),
120);
122//zeroCells(alpha, outletCells);
123
124#include "createFvOptions.H"
Y[inertIndex] max(0.0)
surfaceScalarField & phi
const scalar pRefValue
const label pRefCell
U
Definition: pEqn.H:72
volScalarField & p
dynamicFvMesh & mesh
engineTime & runTime
Creates and initialises the face-flux field phia.
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
const dictionary & simple
volScalarField & alpha
singlePhaseTransportModel laminarTransport(U, phi)
dimensionedScalar zeroAlpha(dimless/dimTime, Zero)
const labelList & inletCells
Definition: createFields.H:106
dimensionedScalar lambda("lambda", dimTime/sqr(dimLength), laminarTransport)
zeroCells(alpha, inletCells)
dimensionedScalar zeroSensitivity(dimVelocity *dimVelocity, Zero)
dimensionedScalar alphaMax("alphaMax", dimless/dimTime, laminarTransport)
setRefCell(p, pimple.dict(), pRefCell, pRefValue)