Simple Program for Step-by-Step Time Integration in Chemical Kinetics, Applied to Simple Model for Hydrogen Combustion

• 作者: Christopher G. Koutitas Christos G. Karayannis Panagis G. Papadopoulos Panos D. Kiousis Yannis N. Dimitropoulos
• 来源期刊: 物理化学期刊(英文)
• 学科: 化学
• 关键词: Chemical Kinetics STEP-BY-STEP Time Integration ELEMENTARY REACTIONS Hydrogen COMBUSTION IGNITION DELAY Explosion
• 年，卷（期）: wlhxqkyw,（2）
• 页码范围: 99-116
• 页数: 18页
• 分类号: O64

摘要

A simple algorithm is proposed for step-by-step time integration of stiff ODEs in Chemical Kinetics. No predictor-corrector technique is used within each step of the algorithm. It is assumed that species concentrations less than 10-6 mol·L-1 do not activate any chemical reaction. So, within each step, the time steplength Δt of the algorithm is determined from the fastest reaction rate maxR by the formula Δt = 10-6mol·L-1/max R. All the reversible elementary reactions occur simultaneously;however, by a simple book-keeping technique, the updating of species concentrations, within each step of the algorithm, is performed within each elementary reaction separately. The above proposed simple algorithm for Chemical Kinetics is applied to a simple model for hydrogen combustion with only five reversible elementary reactions (Initiation, Propagation, First and Second Branching, Termination by wall destruction) with six species (H2, O2, H, O, HO, H2O). These five reversible reactions are recommended in the literature as the most significant elementary reactions of hydrogen combustion [1] [2]. Based on the proposed here simple algorithm for Chemical Kinetics, applied to the global mechanism of proposed five reversible elementary reactions for hydrogen combustion, a simple and short computer program has been developed with only about 120 Fortran instructions. By this proposed program, the following are obtained: 1) The total species concentration of hydrogen combustion, starting from the sum of initial reactants concentrations [H2] + [O2], gradually diminishes, due to termination reaction by wall destruction, and tends to the final concentration of the product [H2O], that is to the 2/3 of its initial value, in accordance to the established overall stoichiometric reaction of hydrogen combustion 2H2 + O2 → 2H2O. 2) Time-histories for concentrations of main species H2, O2, H, H2O of hydrogen combustion, in explosion and equilibrium regions, obtained by the proposed program, are compared to corresponding ones obtained by a