Operator splitting method for numerical solving the atmospheric pollutant dispersion problem

The paper discusses the numerical modeling of transport and diffusion of air pollutants in the atmospheric boundary layer. There was developed a mathematical model of industrial emissions spread in the atmosphere taking into account the deposition velocity of fine particles. The model is described by multidimensional partial differential equations with appropriate initial and boundary conditions. The basic laws of hydrothermodynamics were used in deriving the model. In order to obtain the numerical solution of the problem, we used one of the splitting methods according to physical processes involved (transport, diffusion and absorption), as well as a second-order implicit finite-difference scheme in time. Analysis of numerical results showed that the developed computational algorithm provides sufficient accuracy of the problem solution compared with field measurement data and it has a certain advantage over other numerical methods. In the course of computational experiments, there was determined the degree of influence of such parameters as wind speed and direction, absorption coefficient and physicomechanical properties of particles on the process of atmospheric air pollutants dispersion.