39 absorptionEmissionModel,
40 greyMeanAbsorptionEmission,
56 coeffsDict_((dict.
subDict(typeName +
"Coeffs"))),
61 fileName(coeffsDict_.lookup(
"lookUpTableFileName")),
66 EhrrCoeff_(
readScalar(coeffsDict_.lookup(
"EhrrCoeff"))),
79 const word& key = iter().keyword();
80 speciesNames_.insert(key, nFunc);
82 coeffs_[nFunc].initialise(dict);
94 if (lookUpTable_.found(iter.key()))
96 label index = lookUpTable_.findFieldIndex(iter.key());
98 Info<<
"specie: " << iter.key() <<
" found on look-up table "
99 <<
" with index: " << index <<
endl;
101 specieIndex_[iter()] = index;
111 specieIndex_[iter()] = 0;
113 Info<<
"specie: " << iter.key() <<
" is being solved" <<
endl;
119 "Foam::radiation::greyMeanAbsorptionEmission(const"
120 "dictionary& dict, const fvMesh& mesh)"
121 ) <<
"specie: " << iter.key()
122 <<
" is neither in look-up table: "
123 << lookUpTable_.tableName()
124 <<
" nor is being solved" <<
nl
132 "Foam::radiation::greyMeanAbsorptionEmission(const"
133 "dictionary& dict, const fvMesh& mesh)"
134 ) <<
"specie ft is not present " <<
nl
156 label nSpecies = speciesNames_.
size();
179 const List<scalar>& species = lookUpTable_.lookUp(ft[i]);
181 for (label n=0; n<nSpecies; n++)
185 if (specieIndex_[n] != 0)
188 Yipi = species[specieIndex_[n]]*p[i]*9.869231e-6;
201 if (coeffs_[n].invTemp())
208 ((((b[5]*Ti + b[4])*Ti + b[3])*Ti + b[2])*Ti + b[1])*Ti
251 mesh_.time().timeName(),
265 E().internalField() = EhrrCoeff_*dQ;