Go to the source code of this file.
|
Info<< "Reading thermophysical properties\"<< endl;autoPtr< rhoReactionThermo > | pThermo (rhoReactionThermo::New(mesh)) |
|
thermo | validate (args.executable(), "h", "e") |
|
SLGThermo | slgThermo (mesh, thermo) |
|
const word | inertSpecie (thermo.lookup("inertSpecie")) |
|
| if (!composition.species().found(inertSpecie)) |
|
Info<< "Creating field rho\"<< endl;volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volScalarField &p=thermo.p();Info<< "\Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::turbulenceModel > | turbulence (compressible::turbulenceModel::New(rho, U, phi, thermo)) |
|
Info<< "Creating combustion model\"<< endl;autoPtr< CombustionModel< rhoReactionThermo > > | combustion (CombustionModel< rhoReactionThermo >::New(thermo, turbulence())) |
|
| forAll (Y, i) |
|
fields | add (thermo.he()) |
|
volScalarField | Qdot (IOobject("Qdot", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), mesh, dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0)) |
|
|
rhoReactionThermo & | thermo = pThermo() |
|
basicSpecieMixture & | composition = thermo.composition() |
|
PtrList< volScalarField > & | Y = composition.Y() |
|
Info<< "Creating field dpdt\"<< endl;volScalarField dpdt(IOobject("dpdt", runTime.timeName(), mesh), mesh, dimensionedScalar("dpdt", p.dimensions()/dimTime, 0));Info<< "Creating field kinetic energy K\"<< endl;volScalarField K("K", 0.5 *magSqr(U));volScalarField p_rgh(IOobject("p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);p_rgh=p - rho *gh;pressureControl pressureControl(p, rho, pimple.dict(), false);mesh.setFluxRequired(p_rgh.name());multivariateSurfaceInterpolationScheme< scalar >::fieldTable | fields |
|
◆ pThermo()
Info<< "Reading thermophysical properties\n" << endl;autoPtr<rhoReactionThermo> pThermo |
( |
rhoReactionThermo::New(mesh) |
| ) |
|
◆ validate()
thermo validate |
( |
args. |
executable(), |
|
|
"h" |
, |
|
|
"e" |
|
|
) |
| |
◆ slgThermo()
◆ inertSpecie()
const word inertSpecie |
( |
thermo. |
lookup"inertSpecie" | ) |
|
◆ if()
if |
( |
!composition. |
species).found(inertSpecie | ) |
|
Definition at line 14 of file createFields.H.
References args, combustion(), Foam::dimTime, dpdt, Foam::endl(), argList::executable(), Foam::exit(), Foam::FatalIOError, FatalIOErrorIn, gh, inertSpecie(), Foam::Info, Foam::magSqr(), mesh, Foam::New(), p, p_rgh(), phi, pimple, pressureControl, rho, runTime, thermo, turbulence, and U().
◆ turbulence()
Info<< "Creating field rho\n" << endl;volScalarField rho( IOobject ( "rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), thermo.rho());volScalarField& p = thermo.p();Info<< "\nReading field U\n" << endl;volVectorField U( IOobject ( "U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh);Info<< "Creating turbulence model\n" << endl;autoPtr<compressible::turbulenceModel> turbulence |
( |
compressible::turbulenceModel::New(rho, U, phi, thermo) |
| ) |
|
◆ combustion()
Info<< "Creating combustion model\n" << endl;autoPtr<CombustionModel<rhoReactionThermo> > combustion |
( |
CombustionModel< rhoReactionThermo > |
::Newthermo, turbulence() | ) |
|
◆ forAll()
◆ add()
◆ Qdot()
volScalarField Qdot |
( |
IOobject("Qdot", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE) |
, |
|
|
mesh |
, |
|
|
dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0) |
|
|
) |
| |
◆ thermo
rhoReactionThermo& thermo = pThermo() |
◆ composition
basicSpecieMixture& composition = thermo.composition() |
PtrList<volScalarField>& Y = composition.Y() |
◆ fields
Info<< "Creating field dpdt\n" << endl;volScalarField dpdt( IOobject ( "dpdt", runTime.timeName(), mesh ), mesh, dimensionedScalar("dpdt", p.dimensions()/dimTime, 0));Info<< "Creating field kinetic energy K\n" << endl;volScalarField K("K", 0.5*magSqr(U));volScalarField p_rgh( IOobject ( "p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh);p_rgh = p - rho*gh;pressureControl pressureControl(p, rho, pimple.dict(), false);mesh.setFluxRequired(p_rgh.name());multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields |