1 #include "readGravitationalAcceleration.H" 4 Info<<
"Reading thermophysical properties\n" <<
endl;
18 <<
"Inert specie " <<
inertSpecie <<
" not found in available species " 38 Info<<
"\nReading field U\n" <<
endl;
54 mesh.setFluxRequired(
p.name());
59 Info<<
"Creating turbulence model\n" <<
endl;
60 autoPtr<compressible::turbulenceModel>
turbulence 71 Info<<
"Creating reaction model\n" <<
endl;
72 autoPtr<CombustionModel<rhoReactionThermo>>
reaction 77 Info<<
"Creating multi-variate interpolation scheme\n" <<
endl;
78 multivariateSurfaceInterpolationScheme<scalar>::fieldTable
fields;
100 #include "createMRF.H" Info<< "Reading thermophysical properties\"<< endl;autoPtr< psiReactionThermo > pThermo(psiReactionThermo::New(mesh))
volScalarField Qdot(IOobject("Qdot", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), mesh, dimensionedScalar(dimEnergy/dimVolume/dimTime, Zero))
const word inertSpecie(thermo.get< word >("inertSpecie"))
errorManipArg< error, int > exit(error &err, const int errNo=1)
SLGThermo slgThermo(mesh, thermo)
forAll(U.boundaryField(), patchi)
Ostream & endl(Ostream &os)
Add newline and flush stream.
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tf1, const word &name, const dimensionSet &dimensions, const bool initCopy=false)
Global function forwards to reuseTmpDimensionedField::New.
GeometricField< vector, fvPatchField, volMesh > volVectorField
multivariateSurfaceInterpolationScheme< scalar >::fieldTable fields
GeometricField< scalar, fvPatchField, volMesh > volScalarField
psiReactionThermo & thermo
const dimensionSet dimVolume(pow3(dimLength))
const word & executable() const noexcept
Name of executable without the path.
Info<< "Creating reaction model\"<< endl;autoPtr< CombustionModel< psiReactionThermo > > reaction(CombustionModel< psiReactionThermo >::New(thermo, 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("pos", dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar("neg", dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
#define FatalIOErrorIn(functionName, ios)
Report an error message using Foam::FatalIOError.
const dimensionSet dimEnergy
const dimensionSet dimDensity
volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mesh, dimensionedVector(dimVelocity, Zero))
const dictionary & simple
PtrList< volScalarField > & Y
const dimensionedScalar rhoMin("rhoMin", dimDensity, Zero, pimple.dict())
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
messageStream Info
Information stream (stdout output on master, null elsewhere)
Foam::argList args(argc, argv)
basicSpecieMixture & composition
Creates and initialises the face-flux field phi.
const dimensionedScalar rhoMax("rhoMax", dimDensity, GREAT, pimple.dict())
IOerror FatalIOError
Error stream (stdout output on all processes), with additional 'FOAM FATAL IO ERROR' header text and ...
static constexpr const zero Zero
Global zero (0)