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 accurate computational studies of [3]. 3) In the first step of the algorithm, the only nonzero species concentrations are those of reactants [H2], [O2]. So, the maximum reaction rate is that of the forward initiation reaction max R = Rif = kif[H2]
References
[1]
Brady, J.E., Russell, J.W. and Holum, J.R. (2000) Chemistry: Matter and Its Changes. 3rd Edition, John Wiley & Sons Inc., New York. https://www.amazon.com/Chemistry-Study-Matter-Changes-2000-01-19/dp/B01F81WBQG
[2]
Burks, T.L. and Oran, E.S. (1981) A Computational Study of the Chemical Kinetics of Hydrogen Combustion. Naval Research Laboratory, NASA, Washington DC. https://apps.dtic.mil/dtic/tr/fulltext/u2/a094348.pdf
[3]
Mott, D.R. (1999) New Quasi-Steady-State and Partial Equilibrium Methods for Integrating Chemically Reacting Systems. Ph.D. Thesis Dissertation, Aerospace Engineering and Scientific Computing, University of Michigan, Ann Arbor, MI. https://pdfs.semanticscholar.org/f26c/513c99e2fe330bd78f209d262a3502c982f2.pdf
[4]
Young, T.R. and Boris, J.P. (1977) A Numerical Technique for Solving Stiff Ordinary Differential Equations Associated with the Chemical Kinetics of Reactive-Flow Problems. The Journal of Physical Chemistry, 81, 2424-2427. https://doi.org/10.1021/j100540a018
[5]
Young Jr., T.R. (1980) CHEMEQ: A Subroutine for Solving Stiff Ordinary Differential Equations. Laboratory for Computational Physics. Naval Research Laboratory, Washington DC. https://apps.dtic.mil/dtic/tr/fulltext/u2/a083545.pdf
[6]
Radhakrishnan, K. (1984) Comparison of Numerical Techniques for Integration of Stiff Ordinary Differential Equations Arising in Combustion Chemistry. Lewis Research Center, Cleveland, OH. https://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19850001758.pdf
[7]
Mott, D.R. and Oran, E.S. (2001) CHEMEQ2 A Solvent for the Stiff Ordinary Differential Equations of Chemical Kinetics. Naval Research Laboratory, Washington DC. https://apps.dtic.mil/dtic/tr/fulltext/u2/a392490.pdf https://doi.org/10.21236/ADA385560
[8]
Sichel, M., Tonello, N., Oran, E.S. and Jones, D.A. (2002) A Two-Step Kinetics Model for Numerical Simulation of Explosions and Detonations in H2-O2 Mixtures. Proceedings of the Royal Society London A, 458, 49-82. https://doi.org/10.1098/rspa.2001.0853
[9]
Westbrook, C.K. and Dryer, F.L. (1984) Chemical Kinetics Modeling of Hydro-Carbon Combustion. Progress in Energy and Combustion Science, 10, 1-57. https://doi.org/10.1016/0360-1285(84)90118-7
[10]
Westbrook, C.K., Mizobushi, Y., Poinsot, T., Smith, P.J. and Warnatz, J. (2005) Computational Combustion. Proceedings of the Combustion Institute, 30, 125-157. https://doi.org/10.1016/j.proci.2004.08.275
[11]
Law, C.K. (2006) Combustion Physics. Cambridge University Press, Cambridge.
[12]
Wang, X. and Law, C.K. (2013) An Analysis of the Explosion Limits of Hydrogen-Oxygen Mixtures. Journal of Chemical Physics, 138, Article ID: 134305. https://doi.org/10.1063/1.4798459
[13]
Liang, W. and Law, C.K. (2018) An Analysis of the Explosion Limits of Hydrogen-Oxygen Mixtures with Non-Linear Chain Reactions. Physical Chemistry, Chemical Physics, 20, 742-751. https://doi.org/10.1039/C7CP05639G
[14]
Zhao, P. (2018) Detailed Kinetics in Combustion. Simulation: Manifestation, Model Reduction and Computational Diagnostics. In: De, S., Agarwal, A., Chaudhuri, S. and Sen, S., Eds., Modeling and Simulation of Turbulent Combustion, Springer, Berlin. https://doi.org/10.1007/978-981-10-7410-3_2
[15]
Friedrich, A., Grune, J., Sempert, K., Kuznetsov, M. and Jordan, T. (2019) Hydrogen Combustion Experiments in a Vertical Semi-Confined Channel. International Journal of Hydrogen Energy, 44, 9041-9049. https://doi.org/10.1016/j.ijhydene.2018.06.098
[16]
Kuznetsov, M. and Grune, J. (2019) Experiments on Combustion Regimes for Hydrogen/Air Mixtures in a Thin Layer Geometry. International Journal of Hydrogen Energy, 44, 8727-8742. https://doi.org/10.1016/j.ijhydene.2018.11.144
[17]
Faizal, M., Chuah, L.S., Lee, C., Hameed, A., Lee, J. and Shankar, M. (2019) Review of Hydrogen Fuel for Internal Combustion Engines. Journal of Mechanical Engineering Research and Developments, 42, 35-46. https://doi.org/10.26480/jmerd.03.2019.35.46
[18]
Wikipedia (2020) Hydrogen Fuel. https://en.wikipedia.org/wiki/Hydrogen_fuel
[19]
Enviroment and Energy Conferece. Hydrogen Days (2020) Organizer: Czech Hydrogen Technology Platform (HYTEP). To Be Held 25-27 March 2020, Prague, Czech Republic. https://www.hydrogendays.cz/2020/
