1 #include "readGravitationalAcceleration.H" 3 Info<<
"Reading thermophysical properties\n" <<
endl;
17 <<
"Inert specie " <<
inertSpecie <<
" not found in available species " 31 IOobject::READ_IF_PRESENT,
37 Info<<
"\nReading field U\n" <<
endl;
53 mesh.setFluxRequired(
p.name());
58 Info<<
"Creating turbulence model\n" <<
endl;
59 autoPtr<compressible::turbulenceModel>
turbulence 70 Info<<
"Creating combustion model\n" <<
endl;
71 autoPtr<CombustionModel<psiReactionThermo>>
combustion 76 multivariateSurfaceInterpolationScheme<scalar>::fieldTable
fields;
102 #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<< "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))
Info<< "Creating combustion model\"<< endl;autoPtr< CombustionModel< psiReactionThermo > > combustion(CombustionModel< psiReactionThermo >::New(thermo, turbulence()))
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)