Improvement of a nonnegative preserved efficient solver for atmospheric chemical kinetic equations

  • Authors

    • Atika RADID University Hassan II Casablanca FSAC-MACS
    • Karim RHOFIR University Sultan Moulay Slimane ENSA-LISERT
    2018-07-06
    https://doi.org/10.14419/ijet.v7i4.24694
  • Atmospheric Chemical kinetic, Stiff ODE’s, Positivity preservation, Iterative-Modified-Backward-Euler, Iterated QSSA.

  • Abstract

    Generally, chemical reactions from atmospheric chemistry models are described by a strongly coupled, stiff and nonlinear system of ordinary differential equations, which requires a good numerical solver. Several articles published about the solvers of chemical equations, during the numerical simulation, indicate that one renders the concentration null when it becomes negative. In order to preserve the positivity of the exact solutions, recent works have proposed a new solver called Modified-Backward-Euler (MBE). To improve this solver, we propose in this paper an iterative numerical scheme witch is better fitted to stiff problems. This new approach, called Iterative-Modified-Backward-Euler (IMBE), is based on iterative solution of the P-L structure of the implicit nonlinear ordinary differential equations on each time step. The efficiency of the iteration process is increased by using the Gauss and Successive-Over-Relaxation (SOR). In the case of fast/slow chemical kinetic reactions, we proposed an other variant called Iterative-Quasi-Steady-State-Approximation (IQSSA). The numerical exploration of stiff test problem shows clearly that this formalism is applicable to a wide range of chemical kinetics problems and give a good approximation compared to the recent solver. The numerical procedures give reasonable accurate solutions when compared to exact solution.

    Generally, chemical reactions from atmospheric chemistry models are described by a strongly coupled, stiff and nonlinear system of ordinary differential equations, which requires a good numerical solver. Several articles published about the solvers of chemical equations, during the numerical simulation, indicate that one renders the concentration null when it becomes negative. In order to preserve the positivity of the exact solutions, recent works have proposed a new solver called Modified-Backward-Euler (MBE). To improve this solver, we propose in this paper an iterative numerical scheme witch is better fitted to stiff problems. This new approach, called Iterative-Modified-Backward-Euler (IMBE), is based on iterative solution of the P-L structure of the implicit nonlinear ordinary differential equations on each time step. The efficiency of the iteration process is increased by using the Gauss and Successive-Over-Relaxation (SOR). In the case of fast/slow chemical kinetic reactions, we proposed an other variant called Iterative-Quasi-Steady-State-Approximation (IQSSA). The numerical exploration of stiff test problem shows clearly that this formalism is applicable to a wide range of chemical kinetics problems and give a good approximation compared to the recent solver. The numerical procedures give reasonable accurate solutions when compared to exact solution.

  • References

    1. [1] C.J. Aro , CHEMSODE: a stiffODE solver for the equations of chemical kinetics, Comput. Phys. Commun. 97 (1996) 304–314 .

      [2] F. Feng , Z.F. Wang , J. Li , G.R. Carmichael ,A nonnegativity preserved efficient algorithm for atmospheric chemical kinetic equations,Appl. Math. Comput., 271 (2015), 519–531 .

      [3] Fan Feng, Xuebin Chi, Zifa Wang, Jie Li, Jinrong Jiang, Wenyi Yang, A nonnegativity preserved efficient chemical solver applied to the air pollution forecast, Applied Mathematics and Computation, 314 (2017), 44-–57.

      [4] W. Gong and H.R. Cho. A numerical scheme for the integration of the gas phase chemical rate equations in 3D atmospheric models. Atmospheric Environment, 27A (1993), 2160–2147.

      [5] E. Hesstvedt, Ø. Hov, and I.S.A. Isaksen. Quasi-steady-state approximations in air pollution modeling: Comparison of two numerical schemes for oxidant prediction. International Journal of Chemical Kinetics, 10 (1978) 971—994.

      [6] A.C. Hindmarsh , LSODE and LSODI, two new initial value ordinary differential equation solvers, ACM-SIGNUM Newslett, 15 (1980), 10-–11

      [7] M.R. Houyoux , J.M. Vukovich , Updates to the sparse matrix operator kernel emissions (SMOKE) modeling system and integration with models-3, in: Proceedings of the Emission Inventory: Regional Strategies for the Future, Air Waste Management Association, Raleigh, NC, 1999, p. 1461

      [8] L. O. Jay, A. Sandu, F. A. Potra, and G. R. Carmichael, Improved Quasi-Steady-State-Approximation Methods for Atmospheric Chemistry Integration, SIAM J. Sci. Comput., 18(1) (1997), 182–202.

      [9] D.R. Mott , E.S. Oran , B.V. Leer, A quasi-steady-state solver for the stiff ordinary differential equations of reaction kinetics, J. Comput. Phys. 164 (2000), 407–428

      [10] Sandu, A. Positive numerical integration methods for chemical kinetic systems. J. Comput. Phys. 170 (2001), 1–14

      [11] A. Sandu,J. Verwer, M. van Loon, G. Carmichael, F. Potra, D. Dabdub, J. Seinfeld, Benchmarking stiff ODE solvers for atmospheric chemistry problems.

      I: Implicit versus explicit. Atmos. Environ. 31 (1997), 3151—3166.

      [12] A. Sandu , J.G. Verwer , J.G. Blom , E.J. Spee , G.R. Carmichael , F.A. Potra , Benchmarking stiff ode solvers for atmospheric chemistry problems II:

      Rosenbrock solvers, Atmos. Environ. 31 (1997) 3459–3472 .

      [13] J.H. Seinfeld , S.N. Pandis , Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, second ed., John Wiley and Sons, Inc., 2006 .

      [14] R. D. Saylor and G. D. Ford On the comparison of numerical methods for the integration of kinetic equations in atmospheric chemistry and transport

      models, Atmos. Environ., 29 (1995), 2585- 2593.

      [15] V.N. Shulyk, O.V. Klymenko, I.B. Svir, Numerical solution of stiff ODEs describing complex homogeneous chemical processes. J. Math. Chem. 43 (1)

      (2008), 252-–264.

      [16] A; M. Zhabotinsky, Periodic kinetics of oxidation acid in solution, Biofizika 0 (1964), 306–311

      [17] J. G. Verwer, Gauss–Seidel Iteration for Stiff ODES from Chemical Kinetics, SIAM J. Sci. Comput., 15(5), 1994, 1243—1250

      [18] R. Wolke and O. Knoth. Implicit-explicit Rung-Kutta method applied to atmospheric chemistry modeling. Environnemental Modeling and Software 15

      (2000), 711–719.

  • Downloads

  • How to Cite

    RADID, A., & RHOFIR, K. (2018). Improvement of a nonnegative preserved efficient solver for atmospheric chemical kinetic equations. International Journal of Engineering & Technology, 7(4), 6657-6664. https://doi.org/10.14419/ijet.v7i4.24694

    Received date: 2018-12-23

    Accepted date: 2019-06-08

    Published date: 2018-07-06