Mathematical model for studying the sulphur pollution over Europe

The long-range transport of sulphur pollutants in the atmosphere is an important environmental problem. A mathematical model, consisting of a system of partial differential equations, has been developed at the Danish Agency of Environmental Protection for studying this phenomenon. The model is handled numerically. A pseudospectral (Fourier) algorithm is used in the space discretization phase. A variable stepsize variable formula method (VSVFM), based on predictor corrector schemes, may be specified during the time-discretization part. The time-integration formulae are especially constructed so that: (a) the predictor-corrector schemes have good absolute stability properties on the imaginary axis and (b) the fundamental properties (consistency, zero-stability and convergence) of the global VSVFM are not affected by the changes of the stepsize and/or the predictor-corrector scheme. Both the accuracy requirements and the stability requirements are taken into account in the stepsize selection strategy as well as in the formula selection strategy. This double check leads to a great efficiency in the numerical treatment of problems arising in atmospheric environments. This is demonstrated by many numerical results. The question of the reliability of the results obtained in runs of real physical problems is discussed. Some results obtained in the comparison of calculated results with measurements are reported. An illustration, where the model is applied in studying the air pollution in the area around Denmark due to Danish sources, is presented.

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

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

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

[4]  G. Dahlquist A special stability problem for linear multistep methods , 1963 .

[5]  A. Hindmarsh LSODE and LSODI, two new initial value ordinary differential equation solvers , 1980, SGNM.

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

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

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

[9]  Zahari Zlatev,et al.  Variable stepsize variable formula methods based on predictor-corrector schemes , 1985 .

[10]  J. Lambert Computational Methods in Ordinary Differential Equations , 1973 .

[11]  L. Shampine,et al.  Computer solution of ordinary differential equations : the initial value problem , 1975 .

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

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

[14]  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 .

[15]  G. Dahlquist Convergence and stability in the numerical integration of ordinary differential equations , 1956 .

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

[17]  G. Dahlquist Stability and error bounds in the numerical integration of ordinary differential equations , 1961 .

[18]  H. Stetter Analysis of Discretization Methods for Ordinary Differential Equations , 1973 .