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 Urel\n" << endl;
17(
18 IOobject
19 (
20 "Urel",
21 runTime.timeName(),
22 mesh,
23 IOobject::MUST_READ,
24 IOobject::AUTO_WRITE
25 ),
26 mesh
27);
28
29Info<< "Reading/calculating face flux field phi\n" << endl;
31(
32 IOobject
33 (
34 "phi",
35 runTime.timeName(),
36 mesh,
37 IOobject::READ_IF_PRESENT,
38 IOobject::AUTO_WRITE
39 ),
41);
42
43label pRefCell = 0;
44scalar pRefValue = 0.0;
46mesh.setFluxRequired(p.name());
47
48Info<< "Creating SRF model\n" << endl;
49autoPtr<SRF::SRFModel> SRF
50(
51 SRF::SRFModel::New(Urel)
52);
53
54// Create the absolute velocity
55volVectorField U
56(
57 IOobject
58 (
59 "U",
60 runTime.timeName(),
61 mesh,
62 IOobject::NO_READ,
63 IOobject::AUTO_WRITE
64 ),
65 Urel + SRF->U()
66);
67
68singlePhaseTransportModel laminarTransport(U, phi);
69
70autoPtr<incompressible::turbulenceModel> turbulence
71(
72 incompressible::turbulenceModel::New(U, phi, laminarTransport)
73);
74
75#include "createFvOptions.H"
surfaceScalarField & phi
const scalar pRefValue
const label pRefCell
pimpleControl & pimple
U
Definition: pEqn.H:72
volScalarField & p
dynamicFvMesh & mesh
engineTime & runTime
compressible::turbulenceModel & turbulence
Urel
Definition: pEqn.H:56
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)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:372
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > linearInterpolate(const GeometricField< Type, fvPatchField, volMesh > &vf)
Definition: linear.H:112
singlePhaseTransportModel laminarTransport(U, phi)
Info<< "Reading field p\n"<< endl;volScalarField p(IOobject("p", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field Urel\n"<< endl;volVectorField Urel(IOobject("Urel", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading/calculating face flux field phi\n"<< endl;surfaceScalarField phi(IOobject("phi", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), linearInterpolate(Urel) &mesh.Sf());label pRefCell=0;scalar pRefValue=0.0;setRefCell(p, pimple.dict(), pRefCell, pRefValue);mesh.setFluxRequired(p.name());Info<< "Creating SRF model\n"<< endl;autoPtr< SRF::SRFModel > SRF(SRF::SRFModel::New(Urel))
setRefCell(p, pimple.dict(), pRefCell, pRefValue)