39 absorptionEmissionModel,
40 wideBandAbsorptionEmission,
56 coeffsDict_((dict.
subDict(typeName +
"Coeffs"))),
61 fileName(coeffsDict_.lookup(
"lookUpTableFileName")),
80 dict.
lookup(
"bandLimits") >> iBands_[nBand];
81 dict.
lookup(
"EhrrCoeff") >> iEhrrCoeffs_[nBand];
82 totalWaveLength_ += iBands_[nBand][1] - iBands_[nBand][0];
88 const word& key = iter().keyword();
91 speciesNames_.insert(key, nSpec);
95 if (!speciesNames_.found(key))
99 "Foam::radiation::wideBandAbsorptionEmission(const"
100 "dictionary& dict, const fvMesh& mesh)"
101 ) <<
"specie: " << key <<
"is not in all the bands"
105 coeffs_[nSpec][nBand].initialise(specDicts.
subDict(key));
118 if (lookUpTable_.found(iter.key()))
120 label index = lookUpTable_.findFieldIndex(iter.key());
121 Info<<
"specie: " << iter.key() <<
" found in look-up table "
122 <<
" with index: " << index <<
endl;
123 specieIndex_[iter()] = index;
132 specieIndex_[iter()] = 0.0;
134 Info<<
"species: " << iter.key() <<
" is being solved" <<
endl;
140 "radiation::wideBandAbsorptionEmission(const"
141 "dictionary& dict, const fvMesh& mesh)"
142 ) <<
"specie: " << iter.key()
143 <<
" is neither in look-up table : "
144 << lookUpTable_.tableName() <<
" nor is being solved"
167 label nSpecies = speciesNames_.
size();
190 const List<scalar>& species = lookUpTable_.lookUp(ft[i]);
192 for (label n=0; n<nSpecies; n++)
196 if (specieIndex_[n] != 0)
199 Yipi = species[specieIndex_[n]]*p[i]*9.869231e-6;
211 coeffs_[n][bandI].coeffs(T[i]);
213 if (coeffs_[n][bandI].invTemp())
221 ((((b[5]*Ti + b[4])*Ti + b[3])*Ti + b[2])*Ti + b[1])*Ti
275 if (
mesh().foundObject<volScalarField>(
"hrr"))
278 E().internalField() =
281 *(iBands_[bandI][1] - iBands_[bandI][0])
295 return ILambda*(iBands_[i][1] - iBands_[i][0])/totalWaveLength_;
307 for (label j=0; j<nBands_; j++)
309 Info<<
"Calculating absorption in band: " << j <<
endl;
310 aLambda[j].internalField() = this->a(j);
311 Info<<
"Calculated absorption in band: " << j <<
endl;
313 aLambda[j].internalField()
314 *(iBands_[j][1] - iBands_[j][0])