Go to the source code of this file.
|
| Info<< "Creating combustion model\n"<< endl;autoPtr< combustionModels::psiCombustionModel > | combustion (combustionModels::psiCombustionModel::New(mesh)) |
| |
| thermo | validate (args.executable(),"h","e") |
| |
| SLGThermo | slgThermo (mesh, thermo) |
| |
| const word | inertSpecie (thermo.lookup("inertSpecie")) |
| |
| | if (!composition.contains(inertSpecie)) |
| |
| volScalarField | rho (IOobject("rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), thermo.rho()) |
| |
| Info<< "\nReading field U\n"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);mesh.setFluxRequired(p.name());dimensionedScalar rhoMax(dimensionedScalar::lookupOrDefault("rhoMax", pimple.dict(), dimDensity, GREAT));dimensionedScalar rhoMin(dimensionedScalar::lookupOrDefault("rhoMin", pimple.dict(), dimDensity, 0));Info<< "Creating turbulence model\n"<< endl;autoPtr< compressible::turbulenceModel > | turbulence (compressible::turbulenceModel::New(rho, U, phi, thermo)) |
| |
| combustion | setTurbulence (turbulence()) |
| |
| | forAll (Y, i) |
| |
| fields | add (thermo.he()) |
| |
| volScalarField | dQ (IOobject("dQ", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)) |
| |
|
| psiReactionThermo & | thermo = combustion->thermo() |
| |
| basicSpecieMixture & | composition = thermo.composition() |
| |
| PtrList< volScalarField > & | Y = composition.Y() |
| |
| volScalarField & | p = thermo.p() |
| |
| 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 |
| |
| Info<< "Creating combustion model\n" << endl;autoPtr<combustionModels::psiCombustionModel> combustion |
( |
combustionModels::psiCombustionModel::New(mesh) |
| ) |
|
| thermo validate |
( |
args. |
executable(), |
|
|
"h" |
, |
|
|
"e" |
|
|
) |
| |
| const word inertSpecie |
( |
thermo. |
lookup"inertSpecie" | ) |
|
| if |
( |
!composition. |
containsinertSpecie | ) |
|
| volScalarField rho |
( |
IOobject("rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE) |
, |
|
|
thermo. |
rho() |
|
) |
| |
| Info<< "\nReading field U\n" << endl;volVectorField U( IOobject ( "U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh);mesh.setFluxRequired(p.name());dimensionedScalar rhoMax( dimensionedScalar::lookupOrDefault ( "rhoMax", pimple.dict(), dimDensity, GREAT ));dimensionedScalar rhoMin( dimensionedScalar::lookupOrDefault ( "rhoMin", pimple.dict(), dimDensity, 0 ));Info<< "Creating turbulence model\n" << endl;autoPtr<compressible::turbulenceModel> turbulence |
( |
compressible::turbulenceModel::New(rho, U, phi, thermo) |
| ) |
|
| volScalarField dQ |
( |
IOobject("dQ", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE) |
, |
|
|
mesh |
, |
|
|
dimensionedScalar("dQ", dimEnergy/dimTime, 0.0) |
|
|
) |
| |
| basicSpecieMixture& composition = thermo.composition() |
| PtrList<volScalarField>& Y = composition.Y() |
| volScalarField& p = thermo.p() |
| 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 |