diff --git a/src/combustionModels/diffusionMulticomponent/diffusionMulticomponent.C b/src/combustionModels/diffusionMulticomponent/diffusionMulticomponent.C index 5216838fa1..7d740ba875 100644 --- a/src/combustionModels/diffusionMulticomponent/diffusionMulticomponent.C +++ b/src/combustionModels/diffusionMulticomponent/diffusionMulticomponent.C @@ -106,10 +106,10 @@ diffusionMulticomponent::init() const List& rhs = reactions_[k].rhs(); 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 Wox = specieThermo_[oxydantIndex].W(); + const scalar Wox = specieThermo_[oxidantIndex].W(); forAll(lhs, i) { @@ -130,7 +130,7 @@ diffusionMulticomponent::init() Info << "Fuel heat of combustion : " << qFuel_[k] << endl; s_[k] = - (Wox*mag(specieStoichCoeffs[oxydantIndex])) + (Wox*mag(specieStoichCoeffs[oxidantIndex])) / (Wu*mag(specieStoichCoeffs[fuelIndex])); Info << "stoichiometric oxygen-fuel ratio : " << s_[k] << endl; @@ -169,7 +169,7 @@ diffusionMulticomponent RijPtr_(reactions_.size()), Ci_(reactions_.size(), 1.0), fuelNames_(this->coeffs().lookup("fuels")), - oxidantNames_(this->coeffs().lookup("oxydants")), + oxidantNames_(this->coeffs().lookup("oxidants")), qFuel_(reactions_.size()), stoicRatio_(reactions_.size()), s_(reactions_.size()), @@ -275,13 +275,13 @@ diffusionMulticomponent::correct() for (label k=0; k < nReactions; k++) { const label fuelIndex = species[fuelNames_[k]]; - const label oxydantIndex = species[oxidantNames_[k]]; + const label oxidantIndex = species[oxidantNames_[k]]; const volScalarField& Yfuel = this->thermo().composition().Y(fuelIndex); const volScalarField& Yox = - this->thermo().composition().Y(oxydantIndex); + this->thermo().composition().Y(oxidantIndex); const volScalarField ft ( diff --git a/src/combustionModels/diffusionMulticomponent/diffusionMulticomponent.H b/src/combustionModels/diffusionMulticomponent/diffusionMulticomponent.H index 54559b120a..e859e8e4bb 100644 --- a/src/combustionModels/diffusionMulticomponent/diffusionMulticomponent.H +++ b/src/combustionModels/diffusionMulticomponent/diffusionMulticomponent.H @@ -35,14 +35,15 @@ Description where: - D : is a model constant. - muEff : is the effective turbulent diffusivity - prob : is a Gaussian shaped distribution around the stoichiometric value - of each reaction. The distribtion has the input parameter 'sigma' - for standard deviation. + D : is a model dynamic constant defined as C*f02 where: + C is a model constant + f02 = 1 + sqr(O2/oxidantRes), oxidantRes is an user input - The variable prob is multiplied by the factor: 1 + pow2(O2/oxidantRes), - where oxidantRes is the residual oxydant specified for each reaction. + muEff : is the effective turbulent viscosity + 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: @@ -50,7 +51,7 @@ Description { Ci (1.0 1.0); // Default to 1 fuels (CH4 CO); - oxydants (O2 O2); + oxidants (O2 O2); YoxStream (0.23 0.23); // Default to 0.23 YfStream (1.0 1.0); // Default to 1.0 sigma (0.02 0.02); // Default to 0.02 @@ -59,13 +60,6 @@ Description 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 diffusionMulticomponent.C @@ -112,7 +106,7 @@ class diffusionMulticomponent //- List of fuels for each reaction wordList fuelNames_; - //- List of oxydants for each reaction + //- List of oxidants for each reaction wordList oxidantNames_; //- Heat of combustion [J/Kg]