Numerical modeling of atmospheric pollution in the approaches of random selection of discrete dusting sites and interval distribution of dust size

The results of numerical simulation of atmospheric pollution in Apatity are presented with variations in the dusting area of discrete spatially spaced areas selected randomly and the wind flow velocity. CFD modeling in the volumetric formulation was performed using the COMSOL program. To calculate the aerodynamic characteristics, an incompressible fluid approximation was used using the standard ( k -ε) turbulence model. The process of fine dust propagation is modeled by numerical solution of the convective-diffusion impurity transfer equation taking into account the deposition rate. Numerical experiments (with a total number of more than 1,400) were carried out with a variation of the wind flow velocity from 5 to 23 m/s and a dusting area from 2 to 10 ha with a random selection of 20 discrete sites. Dynamic velocity distributions for specific areas of dusting, interval and total spatial distributions of dust pollution (dust particles with a diameter from 0 to 70 microns in increments of 10 microns) are obtained. The peculiarities of the influence on the levels of atmospheric pollution of specific areas of Apatity depending on the location of dusty areas on the surface of the beach of the tailings dump are noted. The calculated levels of atmospheric pollution in the center of Apatity averaged by the number of combinations of dusting sites are analyzed and generalized to functional dependencies. The calculated dependences of the dust concentration on the dusting area at a fixed wind speed are described by linear functions. The dependence of the dust concentration on the wind flow velocity at a fixed dusting area can be approximated by a power function. The generalized functional dependence makes it possible to predict the dust concentration in Apatity depending on the dusting area of randomly selected discrete areas on the surface of the tailings dump and the wind flow velocity. The obtained dependence permits to make a forecast of the critical dusting area at which the level of atmospheric pollution reaches the maximum permissible concentrations, depending on the speed of the wind flow.