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(SRF::SRFModel::New(Urel));
50
51// Construct the absolute velocity
52volVectorField U
53(
54 IOobject
55 (
56 "U",
57 runTime.timeName(),
58 mesh,
59 IOobject::NO_READ,
60 IOobject::AUTO_WRITE
61 ),
62 Urel + SRF->U()
63);
64
65singlePhaseTransportModel laminarTransport(U, phi);
66
67autoPtr<incompressible::turbulenceModel> turbulence
68(
69 incompressible::turbulenceModel::New(U, phi, laminarTransport)
70);
71
72#include "createFvOptions.H"
surfaceScalarField & phi
const scalar pRefValue
const label pRefCell
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
const dictionary & simple
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)