createFields.H
Go to the documentation of this file.
1 #include "readGravitationalAcceleration.H"
2 
3 Info<< "Reading thermophysical properties\n" << endl;
4 
5 autoPtr<fluidThermo> pThermo
6 (
7  fluidThermo::New(mesh)
8 );
9 fluidThermo& thermo = pThermo();
10 
11 volScalarField rho
12 (
13  IOobject
14  (
15  "rho",
16  runTime.timeName(),
17  mesh,
18  IOobject::NO_READ,
19  IOobject::AUTO_WRITE
20  ),
21  thermo.rho()
22 );
23 
24 Info<< "\nReading field U\n" << endl;
25 volVectorField U
26 (
27  IOobject
28  (
29  "U",
30  runTime.timeName(),
31  mesh,
32  IOobject::MUST_READ,
33  IOobject::AUTO_WRITE
34  ),
35  mesh
36 );
37 
38 #include "compressibleCreatePhi.H"
39 
40 Info<< "Creating turbulence model\n" << endl;
41 autoPtr<compressible::momentumTransportModel> turbulence
42 (
43  compressible::momentumTransportModel::New
44  (
45  rho,
46  U,
47  phi,
48  thermo
49  )
50 );
51 
52 Info<< "Constructing clouds" << endl;
53 parcelCloudList clouds(rho, U, g, thermo);
Foam::compressible::New
autoPtr< CompressibleMomentumTransportModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const viscosity &viscosity)
Definition: compressibleMomentumTransportModelTemplates.C:34
pThermo
Info<< "Reading thermophysical properties\"<< endl;autoPtr< fluidReactionThermo > pThermo(fluidReactionThermo::New(mesh))
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
Foam::volVectorField
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:59
mesh
fvMesh & mesh
Definition: createMeshesPostProcess.H:9
Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:58
rho
rho
Definition: createFields.H:81
phi
phi
Definition: correctPhi.H:3
turbulence
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::momentumTransportModel > turbulence(compressible::momentumTransportModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
U
volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mesh, dimensionedVector(dimVelocity, Zero))
Foam::Info
messageStream Info
g
const dimensionedVector & g
Definition: setRegionFluidFields.H:27
compressibleCreatePhi.H
Creates and initialises the face-flux field phi.