A comparison of numerical techniques for solution of atmospheric kinetic equations

Abstract Numerical modeling of atmospheric chemistry is a computationally intensive problem. The equations describing the interaction among various modeled chemical species are coupled, nonlinear ordinary differential equations. Spatial dependencies in comprehensive three-dimensional air quality models require the solution of this system at thousands of spatial points. Even with increasing computer power, there is a need for efficient and accurate numerical solvers with expanded capabilities, since the next generation of air quality simulation models needs to address the increasingly complex chemistry issues emerging in new model applications. Variants of the commonly used quasi steady-state approximation and the hybrid methods currently used in several modeling systems are examined against a reference mechanism describing chemical interactions related to tropospheric oxidant and acid formation. Additional modifications to the methods are incorporated to yield more robust integration techniques. The chemistry solution methodology used in the regional acid deposition model is also incorporated in this comparison as a base methodology for representing the reference chemical mechanism. The methods are tested against the Gear integration scheme for a variety of test cases including traditional box-model calculations and detailed three-dimensional simulations, and their relative accuracies and efficiencies are investigated. Performance and implementation issues related to chemical integration schemes are examined in the context of the demands and needs of the chemistry component of future comprehensive atmospheric chemistry/transport simulation models.

[1]  I. Isaksen,et al.  Quasi‐steady‐state approximations in air pollution modeling: Comparison of two numerical schemes for oxidant prediction , 1978 .

[2]  Adrian Sandu,et al.  Benchmarking Stiff ODE Solvers for Atmospheric Chemistry Problems I: Implicit versus Explicit , 1996 .

[3]  W. Stockwell,et al.  The second generation regional acid deposition model chemical mechanism for regional air quality modeling , 1990 .

[4]  Daewon W. Byun,et al.  The next generation of integrated air quality modeling: EPA's models-3 , 1996 .

[5]  S. Winkler,et al.  Air quality predictions of the urban airshed model containing improved advection and chemistry algorithms , 1996 .

[6]  Sandra L. Winkler,et al.  An Implicit-Explicit Hybrid Solver for a System of Stiff Kinetic Equations , 1994 .

[7]  Joke Blom,et al.  A comparison of stiff ode solvers for atmospheric chemistry problems , 1995 .

[8]  D P Chock,et al.  Comparison of stiff chemistry solvers for air quality modeling. , 1994, Environmental science & technology.

[9]  M. C. Dodge,et al.  A photochemical kinetics mechanism for urban and regional scale computer modeling , 1989 .

[10]  J. Seinfeld,et al.  Extrapolation techniques used in the solution of stiff odes associated with chemical kinetics of air quality models , 1995 .

[11]  J. Boris,et al.  A numerical technique for solving stiff ordinary differential equations associated with the chemical kinetics of reactive-flow problems , 1977 .

[12]  R. J. Yamartino,et al.  The CALGRID mesoscale photochemical grid model—I. Model formulation , 1992 .

[13]  W. Stockwell.  A homogeneous gas phase mechanism for use in a regional acid deposition model , 1986 .

[14]  Leon Lapidus,et al.  Review of Numerical Integration Techniques for Stiff Ordinary Differential Equations , 1970 .

[15]  R. Turco,et al.  SMVGEAR: A sparse-matrix, vectorized gear code for atmospheric models , 1994 .

[16]  Han‐Ru Cho,et al.  A numerical scheme for the integration of the gas-phase chemical rate equations in three-dimensional atmospheric models , 1993 .

[17]  Rick D. Saylor,et al.  On the comparison of numerical methods for the integration of kinetic equations in atmospheric chemistry and transport models , 1995 .

[18]  C F Curtiss,et al.  Integration of Stiff Equations. , 1952, Proceedings of the National Academy of Sciences of the United States of America.

[19]  J. Christensen,et al.  Test of two numerical schemes for use in atmospheric transport-chemistry models , 1993 .

[20]  C. W. Gear,et al.  Numerical initial value problem~ in ordinary differential eqttations , 1971 .