STY: Expanding header for model description

This commit is contained in:
sergio
2016-01-05 15:12:02 -08:00
parent c6fe2f29fb
commit 16e3f5c1e4
2 changed files with 16 additions and 22 deletions

View File

@ -106,10 +106,10 @@ diffusionMulticomponent<CombThermoType, ThermoType>::init()
const List<specieCoeffs>& rhs = reactions_[k].rhs(); const List<specieCoeffs>& rhs = reactions_[k].rhs();
const label fuelIndex = species[fuelNames_[k]]; const label fuelIndex = species[fuelNames_[k]];
const label oxydantIndex = species[oxidantNames_[k]]; const label oxidantIndex = species[oxidantNames_[k]];
const scalar Wu = specieThermo_[fuelIndex].W(); const scalar Wu = specieThermo_[fuelIndex].W();
const scalar Wox = specieThermo_[oxydantIndex].W(); const scalar Wox = specieThermo_[oxidantIndex].W();
forAll(lhs, i) forAll(lhs, i)
{ {
@ -130,7 +130,7 @@ diffusionMulticomponent<CombThermoType, ThermoType>::init()
Info << "Fuel heat of combustion : " << qFuel_[k] << endl; Info << "Fuel heat of combustion : " << qFuel_[k] << endl;
s_[k] = s_[k] =
(Wox*mag(specieStoichCoeffs[oxydantIndex])) (Wox*mag(specieStoichCoeffs[oxidantIndex]))
/ (Wu*mag(specieStoichCoeffs[fuelIndex])); / (Wu*mag(specieStoichCoeffs[fuelIndex]));
Info << "stoichiometric oxygen-fuel ratio : " << s_[k] << endl; Info << "stoichiometric oxygen-fuel ratio : " << s_[k] << endl;
@ -169,7 +169,7 @@ diffusionMulticomponent
RijPtr_(reactions_.size()), RijPtr_(reactions_.size()),
Ci_(reactions_.size(), 1.0), Ci_(reactions_.size(), 1.0),
fuelNames_(this->coeffs().lookup("fuels")), fuelNames_(this->coeffs().lookup("fuels")),
oxidantNames_(this->coeffs().lookup("oxydants")), oxidantNames_(this->coeffs().lookup("oxidants")),
qFuel_(reactions_.size()), qFuel_(reactions_.size()),
stoicRatio_(reactions_.size()), stoicRatio_(reactions_.size()),
s_(reactions_.size()), s_(reactions_.size()),
@ -275,13 +275,13 @@ diffusionMulticomponent<CombThermoType, ThermoType>::correct()
for (label k=0; k < nReactions; k++) for (label k=0; k < nReactions; k++)
{ {
const label fuelIndex = species[fuelNames_[k]]; const label fuelIndex = species[fuelNames_[k]];
const label oxydantIndex = species[oxidantNames_[k]]; const label oxidantIndex = species[oxidantNames_[k]];
const volScalarField& Yfuel = const volScalarField& Yfuel =
this->thermo().composition().Y(fuelIndex); this->thermo().composition().Y(fuelIndex);
const volScalarField& Yox = const volScalarField& Yox =
this->thermo().composition().Y(oxydantIndex); this->thermo().composition().Y(oxidantIndex);
const volScalarField ft const volScalarField ft
( (

View File

@ -35,14 +35,15 @@ Description
where: where:
D : is a model constant. D : is a model dynamic constant defined as C*f02 where:
muEff : is the effective turbulent diffusivity C is a model constant
prob : is a Gaussian shaped distribution around the stoichiometric value f02 = 1 + sqr(O2/oxidantRes), oxidantRes is an user input
of each reaction. The distribtion has the input parameter 'sigma'
for standard deviation.
The variable prob is multiplied by the factor: 1 + pow2(O2/oxidantRes), muEff : is the effective turbulent viscosity
where oxidantRes is the residual oxydant specified for each reaction. prob : is a normalized Gaussian shaped distribution around the stoichiometric
value of each reaction. The distribtion is controled by 'sigma'
for standard deviation and ftCorr for correction of the stoichiometric
value.
In the combustion properties dictionary: In the combustion properties dictionary:
@ -50,7 +51,7 @@ Description
{ {
Ci (1.0 1.0); // Default to 1 Ci (1.0 1.0); // Default to 1
fuels (CH4 CO); fuels (CH4 CO);
oxydants (O2 O2); oxidants (O2 O2);
YoxStream (0.23 0.23); // Default to 0.23 YoxStream (0.23 0.23); // Default to 0.23
YfStream (1.0 1.0); // Default to 1.0 YfStream (1.0 1.0); // Default to 1.0
sigma (0.02 0.02); // Default to 0.02 sigma (0.02 0.02); // Default to 0.02
@ -59,13 +60,6 @@ Description
laminarIgn false; // Default false laminarIgn false; // Default false
} }
ftCorr is used to shift the location of the flame and corrects the
stochimetric mixture value when the mesh resolution is not enough
to resolve the combustion zone.
NOTE: Optionally the switch 'laminarIgn' can be used to ignite the mixture
using laminar combustion.
SourceFiles SourceFiles
diffusionMulticomponent.C diffusionMulticomponent.C
@ -112,7 +106,7 @@ class diffusionMulticomponent
//- List of fuels for each reaction //- List of fuels for each reaction
wordList fuelNames_; wordList fuelNames_;
//- List of oxydants for each reaction //- List of oxidants for each reaction
wordList oxidantNames_; wordList oxidantNames_;
//- Heat of combustion [J/Kg] //- Heat of combustion [J/Kg]