Go to the source code of this file.
|
Info<< "Reading thermophysical properties\"<< endl;autoPtr< psiReactionThermo > | pThermo (psiReactionThermo::New(mesh)) |
|
thermo | validate (args.executable(), "h", "e") |
|
const word | inertSpecie (thermo.lookup("inertSpecie")) |
|
| if (!composition.species().found(inertSpecie)) |
|
volScalarField | rho (IOobject("rho", runTime.timeName(), mesh), thermo.rho()) |
|
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField &p=thermo.p();pressureControl pressureControl(p, rho, pimple.dict(), false);mesh.setFluxRequired(p.name());Info<< "Creating turbulence model.\"<< nl;autoPtr< compressible::turbulenceModel > | turbulence (compressible::turbulenceModel::New(rho, U, phi, thermo)) |
|
Info<< "Creating reaction model\"<< endl;autoPtr< CombustionModel< psiReactionThermo > > | reaction (CombustionModel< psiReactionThermo >::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)) |
|
|
psiReactionThermo & | 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));multivariateSurfaceInterpolationScheme< scalar >::fieldTable | fields |
|
◆ pThermo()
Info<< "Reading thermophysical properties\n" << endl;autoPtr<psiReactionThermo> pThermo |
( |
psiReactionThermo::New(mesh) |
| ) |
|
◆ validate()
thermo validate |
( |
args. |
executable(), |
|
|
"h" |
, |
|
|
"e" |
|
|
) |
| |
◆ inertSpecie()
const word inertSpecie |
( |
thermo. |
lookup"inertSpecie" | ) |
|
◆ if()
if |
( |
!composition. |
species).found(inertSpecie | ) |
|
Definition at line 12 of file createFields.H.
References args, Foam::dimTime, dpdt, Foam::endl(), argList::executable(), Foam::exit(), Foam::FatalIOError, FatalIOErrorIn, inertSpecie(), Foam::Info, Foam::magSqr(), mesh, Foam::New(), Foam::nl, p, phi, pimple, pressureControl, reaction(), rho, runTime, thermo, turbulence, and U().
◆ rho()
volScalarField rho |
( |
IOobject("rho", runTime.timeName(), mesh) |
, |
|
|
thermo. |
rho() |
|
) |
| |
◆ turbulence()
Info<< "Reading field U\n" << endl;volVectorField U( IOobject ( "U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh);volScalarField& p = thermo.p();pressureControl pressureControl(p, rho, pimple.dict(), false);mesh.setFluxRequired(p.name());Info << "Creating turbulence model.\n" << nl;autoPtr<compressible::turbulenceModel> turbulence |
( |
compressible::turbulenceModel::New(rho, U, phi, thermo) |
| ) |
|
◆ reaction()
Info<< "Creating reaction model\n" << endl;autoPtr<CombustionModel<psiReactionThermo> > reaction |
( |
CombustionModel< psiReactionThermo > |
::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
psiReactionThermo& 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));multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields |