createFieldsPimple.H
Go to the documentation of this file.
1 Info << "Reading field p\n"
2  << endl;
3 pPtr_.reset(
4  new volScalarField(
5  IOobject(
6  "p",
7  runTime.timeName(),
8  mesh,
9  IOobject::MUST_READ,
10  IOobject::AUTO_WRITE),
11  mesh));
12 volScalarField& p = pPtr_();
13 
14 Info << "Reading field U\n"
15  << endl;
16 UPtr_.reset(
17  new volVectorField(
18  IOobject(
19  "U",
20  runTime.timeName(),
21  mesh,
22  IOobject::MUST_READ,
23  IOobject::AUTO_WRITE),
24  mesh));
25 volVectorField& U = UPtr_();
26 
27 #include "createPhiPython.H"
28 phi.write();
29 
30 // actuator related stuff
31 Info << "Creating source term. " << endl;
32 fvSourcePtr_.reset(
33  new volVectorField(
34  IOobject(
35  "fvSource",
36  runTime.timeName(),
37  mesh,
38  IOobject::NO_READ,
39  IOobject::AUTO_WRITE),
40  mesh,
41  dimensionedVector("fvSource", dimensionSet(0, 1, -2, 0, 0, 0, 0), vector::zero),
42  zeroGradientFvPatchField<vector>::typeName));
43 volVectorField& fvSource = fvSourcePtr_();
44 
45 setRefCell(p, pimple.dict(), pRefCell_, pRefValue_);
46 mesh.setFluxRequired(p.name());
47 
48 laminarTransportPtr_.reset(
49  new singlePhaseTransportModel(U, phi));
50 singlePhaseTransportModel& laminarTransport = laminarTransportPtr_();
51 
52 turbulencePtr_.reset(
53  incompressible::turbulenceModel::New(U, phi, laminarTransport));
54 
55 if (hasTField_)
56 {
57  Info << "Reading field T\n"
58  << endl;
59  TPtr_.reset(
60  new volScalarField(
61  IOobject(
62  "T",
63  runTime.timeName(),
64  mesh,
65  IOobject::MUST_READ,
66  IOobject::AUTO_WRITE),
67  mesh));
68 
69  // Laminar Prandtl number
70  PrPtr_.reset(
71  new dimensionedScalar("Pr", dimless, laminarTransport));
72 
73  // Turbulent Prandtl number
74  PrtPtr_.reset(
75  new dimensionedScalar("Prt", dimless, laminarTransport));
76 
77  // kinematic turbulent thermal thermal conductivity m2/s
78  Info << "Reading field alphat\n"
79  << endl;
80  alphatPtr_.reset(
81  new volScalarField(
82  IOobject(
83  "alphat",
84  runTime.timeName(),
85  mesh,
86  IOobject::MUST_READ,
87  IOobject::AUTO_WRITE),
88  mesh));
89 }
U
volVectorField & U
Definition: createRefsPimpleDyM.H:7
pimple
pimpleControlDF & pimple
Definition: createRefsPimpleDyM.H:5
fvSource
volScalarField & fvSource
Definition: createRefsHeatTransfer.H:7
mesh
fvMesh & mesh
Definition: createRefsHeatTransfer.H:4
p
volScalarField & p
Definition: createRefsPimpleDyM.H:6
phi
surfaceScalarField & phi
Definition: createRefsPimpleDyM.H:8
laminarTransport
singlePhaseTransportModel & laminarTransport
Definition: createRefsPimpleDyM.H:9
createPhiPython.H
runTime
Time & runTime
Definition: createRefsHeatTransfer.H:1