|
|
(20 intermediate revisions by the same user not shown) |
Line 1: |
Line 1: |
− | == Chemical kinetics ==
| |
− |
| |
− | This chapter reports the principles that drive the computation of combustion chemistry in most CFD softwares.
| |
− |
| |
− | === Chemkin (.scheme .therm .trans), Cantera (xml) ===
| |
− |
| |
− | ...
| |
− |
| |
− |
| |
− | === Arrhenius law ===
| |
− |
| |
− | <math>\mathcal{A}_j</math> is the pre-exponential factor, <math>\mathcal{\beta}_j</math> is the temperature exponent and <math>E_{a_j}</math> the activation energy
| |
− |
| |
− | <math>
| |
− | k_j = \mathcal{A}_j T^{\mathcal{\beta}_j} \exp \left(-\frac{E_{a_j}}{R T}\right)
| |
− | </math>
| |
− |
| |
− |
| |
− | === Three-body reactions ===
| |
− |
| |
− | In the forward direction, three-body reactions involve two species A and B as reactants and yield a single product AB. In that case, the third body M is used to stabilize the excited product AB*. On the contrary, in the reverse direction, heat provides the energy necessary to break the link between A and B.
| |
− |
| |
− | <math>......</math>
| |
− |
| |
− | The third body M can be any inert molecule.
| |
− |
| |
− |
| |
− | === Falloff reactions ===
| |
− |
| |
− | Under specific conditions, some reaction rate expressions are dependent on pressure and temperature. This is especially true for the rate associated to unimolecular/recombination fall-off reactions which increases with pressure. In such cases, if the chemical process takes place in a high or low pressure limit
| |
− |
| |
− | <code>
| |
− | <!-- reaction 0012 -->
| |
− | <reaction reversible="yes" type="falloff" id="0012">
| |
− | <equation>O + CO (+ M) [=] CO2 (+ M)</equation>
| |
− | <rateCoeff>
| |
− | <Arrhenius>
| |
− | <A>1.800000E+07</A>
| |
− | <b>0</b>
| |
− | <E units="cal/mol">2385.000000</E>
| |
− | </Arrhenius>
| |
− | <Arrhenius name="k0">
| |
− | <A>6.020000E+08</A>
| |
− | <b>0</b>
| |
− | <E units="cal/mol">3000.000000</E>
| |
− | </Arrhenius>
| |
− | <efficiencies default="1.0">AR:0.5 C2H6:3 CH4:2 CO:1.5 CO2:3.5 H2:2 H2O:6 O2:6 </efficiencies>
| |
− | <falloff type="Lindemann"/>
| |
− | </rateCoeff>
| |
− | <reactants>CO:1 O:1.0</reactants>
| |
− | <products>CO2:1.0</products>
| |
− | </reaction>
| |
− | </code>
| |
− |
| |
− | === Reaction rates ===
| |
− |
| |
− | The global rate of a reaction j (evolution in concentration per unit of time) varies depending on the proportion of the rates associated to the forward and backward directions.
| |
− |
| |
− | <math>
| |
− | \mathcal{Q}_j = \mathcal{Q}_{f,j} - \mathcal{Q}_{r,j}
| |
− | </math>
| |
− |
| |
− |
| |
− | === Production/Consumption source terms ===
| |
− |
| |
− | Species <math>Y_k</math> source terms are deduced from
| |
− |
| |
− | <math>
| |
− | \dot{\omega}_k = W_k \sum_{j=1}^{N_R} \nu_{k,j} \mathcal{Q}_j
| |
− | </math>
| |
− |
| |
| == Solver to build reference trajectories == | | == Solver to build reference trajectories == |
| | | |