Running Air Pollution Models on Massively Parallel Machines

Abstract The concentration of air pollutants have in general been steadily increasing during the last three decades. To correctly gauge the impact of various sources of pollutants requires careful modelling of the complex physics processes associated with the chemistry and transport of air pollution. These models are computationally demanding and require today's fastest high performance computers for practical implementations. The damaging effects are normally due to the combined effects of several air pollutants. Therefore a reliable mathematical model must study simultaneously all relevant air pollutants. This requirement increases the size of the air pollution models. The discretization of models that contain many air pollutants leads to huge computational problems. After the discretization of such an air pollution model, systems of several hundred thousands (or even several millions) of equations arise and these have to be treated numerically during many time-steps (typically several thousands). Such big computational problems can successfully be treated only on big modern vector and/or parallel computers. However, access to a fast high-speed computer is not sufficient. One must also ensure that the great potential power of the computer is correctly exploited. Very often this is a rather difficult task. The efforts to solve this task in the case where the computer under consideration is a massively parallel machine will be discussed in this paper. Tests performed on several such computers will be presented.

[1]  Gene H. Golub,et al.  Numerical Linear Algebra, Digital Signal Processing and Parallel Algorithms , 1991, NATO ASI Series.

[2]  Zahari Zlatev,et al.  Comparison of numerical techniques for use in air pollution models with non-linear chemical reactions , 1989 .

[3]  D. Pepper,et al.  Modeling the dispersion of atmospheric pollution using cubic splines and Chapeau functions. [Environmental transport of chemical and radioactive gaseous wastes at Savannah River Plant] , 1977 .

[4]  O. Nevanlinna,et al.  Stability of explicit time discretizations for solving initial value problems , 1982 .

[5]  Zahari Zlatev,et al.  Studying high ozone concentrations by using the Danish Eulerian model , 1993 .

[6]  J. N. Lyness.  Computational Techniques Based on the Lanczos Representation , 1974 .

[7]  Z. Zlatev,et al.  NUMERICAL TREATMENT OF LARGE-SCALE AIR POLLUTION MODELS , 1988 .

[8]  S. Tauber,et al.  Modelling the dispersion of atmospheric pollution using cubic splines and chapeau functions , 1979 .

[9]  R. Courant,et al.  Über die partiellen Differenzengleichungen der mathematischen Physik , 1928 .

[10]  William R. Goodin,et al.  Numerical solution of the atmospheric diffusion equation for chemically reacting flows , 1982 .

[11]  Leen Hordijk,et al.  Integrated analysis of acidification in Europe. , 1985 .

[12]  Zahari Zlatev,et al.  Three-dimensional advection-diffusion modelling for regional scale , 1983 .

[13]  D. Vandeneynde.  Swell Prediction Off the Belgian Coast - An Operational Application of Supercomputing , 1992 .

[14]  Zahari Zlatev,et al.  Mathematical model for studying the sulphur pollution over Europe , 1985 .

[15]  Creon Levit.  Grid Communication on the Connection Machine: Analysis, Performance, and Improvements , 1989, Int. J. High Speed Comput..

[16]  Z. Zlatev,et al.  Studying the Sulphur and Nitrogen Pollution over Europe , 1989 .

[17]  Zahari Zlatev,et al.  Implementation of a variable stepsize variable formula method in the time-integration part of a code for treatment of long-range transport of air pollutants , 1984 .

[18]  Zahari Zlatev,et al.  Stability Restrictions on Time-Stepsize for Numerical Integration of First-Order Partial Differential Equations* , 1983 .

[19]  Steven A. Orszag,et al.  Numerical simulation of incompressible flows within simple boundaries: accuracy , 1971, Journal of Fluid Mechanics.

[20]  Zahari Zlatev.  Running large air pollution models on high speed computers , 1990 .

[21]  Z. Zlatev,et al.  A comparison of the predictions of an eulerian atmospheric transport — chemistry model with experimental measurements over the North sea , 1994 .

[22]  Zahari Zlatev,et al.  A Eulerian air pollution model for Europe with nonlinear chemistry , 1992 .

[23]  A. J. Baker,et al.  a Simple One-Dimensional Finite-Element Algorithm with Multidimensional Capabilities , 1979 .

[24]  Patrick J. Roache,et al.  A pseudo-spectral FFT technique for non-periodic problems , 1978 .

[25]  Zahari Zlatev,et al.  Treatment of some mathematical models describing long-range transport of air pollutants on vector processors , 1988, Parallel Comput..

[26]  H. Kreiss,et al.  Comparison of accurate methods for the integration of hyperbolic equations , 1972 .

[27]  Jack J. Dongarra,et al.  Performance of various computers using standard linear equations software in a FORTRAN environment , 1988, CARN.

[28]  Jack Dongarra,et al.  LINPACK Users' Guide , 1987 .

[29]  Zahari Zlatev,et al.  Testing subroutines solving advection-diffusion equations in atmospheric environments , 1983 .

[30]  Lennart Johnsson.  Supercomputers: Past and Future , 1990 .

[31]  C. Molenkamp,et al.  Accuracy of Finite-Difference Methods Applied to the Advection Equation , 1968 .

[32]  Z. Zlatev,et al.  Advance in the theory of variable stepsize variable formula methods for ordinary differential equations , 1989, Conference on Numerical Ordinary Differential Equations.

[33]  Zahari Zlatev,et al.  Modified Diagonally Implicit Runge–Kutta Methods , 1981 .

[34]  W. P. Crowley,et al.  NUMERICAL ADVECTION EXPERIMENTS1 , 1968 .

[35]  S. Zalesak.  Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[36]  Zahari Zlatev,et al.  Application of predictor-corrector schemes with several correctors in solving air pollution problems , 1984 .

[37]  B. Fornberg.  On a Fourier method for the integration of hyperbolic equations , 1975 .

[38]  Zahari Zlatev,et al.  Computations with symmetric, positive definite and band matrices on a parallel vector processor , 1988, Parallel Comput..

[39]  Zahari Zlatev,et al.  Vectorizing codes for studying long-range transport of air pollutants , 1991 .

[40]  Horst D. Simon.  Scientific Applications of the Connection Machine , 1992 .

[41]  S. Lennart Johnsson,et al.  Optimal communication channel utilization for matrix transposition and related permutations on binary cubes , 1994, Discret. Appl. Math..

[42]  Jack J. Dongarra,et al.  A proposal for an extended set of Fortran Basic Linear Algebra Subprograms , 1985, SGNM.

[43]  Z. Zlatev.  Zero-stability properties of the three-ordinate variable stepsize variable formula methods , 1981 .

[44]  Zahari Zlatev,et al.  Stability properties of variable stepsize variable formula methods , 1978 .

[45]  Charles L. Lawson,et al.  Basic Linear Algebra Subprograms for Fortran Usage , 1979, TOMS.

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