42 absorptionEmissionModel,
43 greyMeanAbsorptionEmission,
61 specieIndex_(
label(0)),
64 EhrrCoeff_(
readScalar(coeffsDict_.lookup(
"EhrrCoeff"))),
67 if (!isA<basicSpecieMixture>(thermo_))
70 <<
"Model requires a multi-component thermo package" 85 const word& key = iter().keyword();
86 speciesNames_.insert(key, nFunc);
88 coeffs_[nFunc].initialise(dict);
92 if (coeffsDict_.found(
"lookUpTableFileName"))
94 const word name = coeffsDict_.lookup(
"lookUpTableFileName");
101 fileName(coeffsDict_.lookup(
"lookUpTableFileName")),
110 <<
"specie ft is not present to use with " 111 <<
"lookUpTableFileName " <<
nl 123 if (!lookUpTablePtr_.empty())
125 if (lookUpTablePtr_().found(iter.key()))
127 label index = lookUpTablePtr_().findFieldIndex(iter.key());
129 Info<<
"specie: " << iter.key() <<
" found on look-up table " 130 <<
" with index: " << index <<
endl;
132 specieIndex_[iter()] = index;
137 specieIndex_[iter()] = 0;
139 Info<<
"specie: " << iter.key() <<
" is being solved" <<
endl;
144 <<
"specie: " << iter.key()
145 <<
" is neither in look-up table: " 146 << lookUpTablePtr_().tableName()
147 <<
" nor is being solved" <<
nl 154 specieIndex_[iter()] = 0;
160 <<
" there is not lookup table and the specie" <<
nl 162 <<
" is not found " <<
nl 193 "aCont" +
name(bandI),
201 extrapolatedCalculatedFvPatchVectorField::typeName
213 if (specieIndex_[n] != 0)
219 const List<scalar>& Ynft = lookUpTablePtr_().lookUp(ft[celli]);
221 Xipi = Ynft[specieIndex_[
n]]*
paToAtm(p[celli]);
228 invWt += mixture.
Y(
s)[celli]/mixture.
W(
s);
232 scalar Xk = mixture.
Y(index)[celli]/(mixture.
W(index)*invWt);
239 scalar Ti = T[celli];
241 if (coeffs_[n].invTemp())
248 ((((b[5]*Ti + b[4])*Ti + b[3])*Ti + b[2])*Ti + b[1])*Ti
253 ta.
ref().correctBoundaryConditions();
274 "ECont" +
name(bandI),
275 mesh_.time().timeName(),
292 E.
ref().primitiveFieldRef() = EhrrCoeff_*Qdot/mesh_.V();
296 E.
ref().primitiveFieldRef() = EhrrCoeff_*
Qdot;
303 <<
"Incompatible dimensions for Qdot field" <<
endl;
310 <<
"Qdot field not found in mesh" <<
endl;
#define forAll(list, i)
Loop across all elements in list.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
A class for handling file names.
tmp< volScalarField > ECont(const label bandI=0) const
Emission contribution for continuous phase.
errorManipArg< error, int > exit(error &err, const int errNo=1)
A list of keyword definitions, which are a keyword followed by any number of values (e...
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
T & ref() const
Return non-const reference or generate a fatal error.
Unit conversion functions.
Type & lookupObjectRef(const word &name) const
Lookup and return the object reference of the given Type.
defineTypeNameAndDebug(cloudAbsorptionEmission, 0)
Ostream & endl(Ostream &os)
Add newline and flush stream.
const speciesTable & species() const
Return the table of species.
bool foundObject(const word &name) const
Is the named Type found?
const Time & time() const
Return the top-level database.
Specialization of basicMultiComponentMixture for a mixture consisting of a number for molecular speci...
Model to supply absorption and emission coefficients for radiation modelling.
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type.
tmp< volScalarField > aCont(const label bandI=0) const
Absorption coefficient for continuous phase.
PtrList< volScalarField > & Y()
Return the mass-fraction fields.
Macros for easy insertion into run-time selection tables.
const word dictName() const
Return the local dictionary name (final part of scoped name)
const dimensionSet dimVolume(pow3(dimLength))
const dimensionSet & dimensions() const
Return dimensions.
gmvFile<< "tracers "<< particles.size()<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().x()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().y()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
const dictionary & optionalSubDict(const word &) const
Find and return a sub-dictionary if found.
const dimensionedScalar b
Wien displacement law constant: default SI units: [m.K].
A class for handling words, derived from string.
Fundamental fluid thermodynamic properties.
const word & constant() const
Return constant name.
tmp< volScalarField > eCont(const label bandI=0) const
Emission coefficient for continuous phase.
bool readScalar(const char *buf, doubleScalar &s)
Read whole of buf as a scalar. Return true if successful.
greyMeanAbsorptionEmission(const dictionary &dict, const fvMesh &mesh)
Construct from components.
forAllConstIter(PtrDictionary< phaseModel >, mixture.phases(), phase)
errorManip< error > abort(error &err)
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
word name(const complex &)
Return a string representation of a complex.
const dimensionSet dimEnergy
dimensionedScalar pow3(const dimensionedScalar &ds)
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
scalar paToAtm(const scalar pa)
Conversion from atm to Pa.
#define WarningInFunction
Report a warning using Foam::Warning.
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
Mesh data needed to do the Finite Volume discretisation.
addToRunTimeSelectionTable(absorptionEmissionModel, cloudAbsorptionEmission, dictionary)
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
const dimensionSet dimMass(1, 0, 0, 0, 0, 0, 0)
A class for managing temporary objects.
virtual scalar W(const label speciei) const =0
Molecular weight of the given specie [kg/kmol].
virtual ~greyMeanAbsorptionEmission()
Destructor.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
autoPtr< radiation::radiationModel > radiation(radiation::radiationModel::New(T))