MSC 00A72 DOI: 10.14529/mmp200401
NUMERICAL STUDY OF THE SUSUPLUME AIR POLLUTION MODEL
S.M. Elsakov1, D.A. Drozin1, A.V. Herreinstein1, T.G. Krupnova1, S.G. Nitskaya1, T.Yu. Olenchikova1, A.A. Zamyshlyaeva1
1 South Ural State University, Chelyabinsk, Russian Federation E-mails: elsakovsm@susu.ru, drozinda@susu.ru, gerenshteinav@susu.ru, krupnovatg@susu.ru, nitckaiasg@susu.ru, olenchikovati@susu.ru, zamyshliaevaaa@susu.ru
In this paper, we propose a SUSUPLUME air pollution as a modern application of the classical Gaussian plume model. The presented model takes into account meteorological conditions and parameters of the pollution sources. The classical model is supplemented by the equations of motion of the center of mass of a single emission. A numerical study has shown that in stationary weather conditions the presented model qualitatively coincides with other known models. The results of calculating the concentrations of pollutants do not contradict the obtained values based on the official methodology for calculating the maximum concentrations of pollutants approved for usege in the territory of the Russian Federation. The SUSUPLUME model contains a number of identifiable parameters and it can be adapted to real conditions. The computational model consists of two blocks: a block for recording measurement information and a block for calculating the concentrations of pollutants. The measurement information registration unit has a low labor intensity (over a million registrations per second). The pollutant concentrations calculation block is laborious (400 points of calculations per second). Concentrations are calculated independently, it allows to use parallelization of the computational process in the future.
Keywords: air pollution model; Gaussian plume model; Romberg's method.
Introduction
At present, the problem of air pollution in large cities is extremely acute. For example, in Chelyabinsk, up to one third of days per year the maximum permissible concentrations of pollutants are exceeded [1]. In this situation the issue of identifying sources of air pollution becomes particularly relevant. The existing methods for monitoring the state of atmospheric air based on empirical methods for calculating the concentrations of pollutants [2] do not meet the current needs. The deployment of an extensive network of physicochemical air condition monitoring and the use of modern methods of mathematical modelling allow taking into account the variety of pollution sources in large cities, as well as the peculiarities of meteorological conditions formed under the influence of urban development and life.
Air pollution models can be divided into: empirical, computational fluid dynamics and semi-empirical models. Models are also distinguished by the scale of the computational grid: mesoscale models (grid spacing at least 10 km), as well as microscale models (grid spacing about 2 meters). Mesoscale models are used to calculate weather changes, e.g. WRF/LSM/Urban modelling system [3], NU-WRF model [4], WRF/Chem-NCSU scale
transport [5, 6] and WRF/hem-ROMS [6] and others. Microscale models are used to calculate air pollution of a city, street, enterprise, for example, the WACFD-RANS model which was applied in Madrid to two different areas of the city [7-10]. On the basis of empirical models, the Russian official methodology for calculating MRR-2017 [11] has been developed. Such models are well applicable for calculating quasi-stationary processes, for example, when the time of movement of pollutants is insignificant compared to the time of their release. The disadvantage of such models is their "rigid" structure and a large number of simplifications, which undoubtedly leads to inaccuracies in the calculation. Moreover, the introduction of additional correction factors does not increase their accuracy [12]. Compared to other models, computational fluid dynamics models allow to take into account complex flows in an urban environment, air flow around obstacles on a micro scale, and much more. For example, the WACFD-RANS model used in Spain [7-9] or the TAPM model developed and used in Australia [13, 14]. CFD models look very promising, but they have one significant drawback. The calculation of these models requires significant computational resources. Semi-empirical models have found wide practical application in calculating air pollution in an urban environment. The most widespread is the PasquillGifford model [15], which is based on the hypothesis of the normal distribution of the pollutants concentration along spatial coordinates. This model was recommended in 1986 for the creation of national local models for IAEA member countries [16]. In Russia, national models DV-2010 [17] and NPO Typhoon [18, 19] were developed to calculate the propagation of radionuclides. Although the Pasquill-Gifford model was developed in the 50s of the last century, it has not lost its relevance up to the present time, in view of its simplicity and ability to adapt to local problems. In this paper, numerical modelling is presented within the framework of the SUSUPLUME model proposed by the authors, which is based on the classical Gauss model, supplemented by the equations of motion of the center of mass of a single emission. The motivation for the development of the model was the creation of a toolkit for an industrial enterprise, which allows calculating the concentration of pollutants within the boundaries of the sanitary zone to improve emission control. The main attention was paid to solving the key problem of constructing a tunable model for specific parameters of plant emissions, in contrast to the widely used AERMOD [20] and CALPUFF [21] models recommended by the US EPA.
1. Model of Air Pollution from a Point Source
Consider a point pollution source. Let it be located in the coordinates (xi, yi5 hi) , i = 1, N, where xi,yi are coordinates in meters in the horizontal plane. The X axis is directed from West to East, and the Y axis goes from South to North, hi is the height of the mouth of the pollution source above the earth's surface in meters, i is the number of the pollution source, N is the total number of point sources. At each moment of time let us know the wind in the horizontal plane. Let u (t) be the projection of the wind speed on the X axis and v (t) be the projection of the wind speed on the Y axis.
Consider a single emission of pollutants during rather small time period At. Let us describe the motion dynamics of the center of mass for the single emission. Let xi (to, t) ,yi (t0,t) , zi (t0,t) be coordinates of the single emission from the i-th source at time t, that left the source mouth at time t0. Then the movement of the center of mass of
the single emission in the horizontal plane is described by the relations
t t Xi (to, t) = Xi + J u (t) dt, yi (to,t) = yi + J v (t) dt.
to to
On the vertical axis, the emission is influenced by the forces of gravity, Archimedes and viscous friction. Write down the corresponding second Newton's law
m,gZi (to,t) = -mgg + mag - ¡imgZi (to,t), (1)
where mg is the mass of the single emission, ma is the mass of displaced air, ^ is the viscosity coefficient, which describes the dependence of the medium resistance to the movement of the single emission. Now write the ideal gas equation for the single emission and displaced air volume. Since the pressures of these gases coincide, then
m9 nrji _ ma j^rp
wRTg —RTa,
where Ta is the air temperature, Tg is the temperature of the single emission, ma is the mass of displaced air, Ma is the molar mass of air, Mg is the molar mass of the single emission, R = 8, 314 J ■ K-1 ■ mol-1 is the universal gas constant. Hence
матд
ma = m9TTTfT- (2)
Mg T a
Since the exhaust temperature Tg can be different from the ambient air temperature Ta, the exhaust cooling process takes place simultaneously. This can be described by the relation
T (to,t) = Y (Ta - T (to,t)) , (3)
where T (t0 ,t) is the temperature of the single emission from the i-th source, y is the coefficient characterizing the rate of emission cooling depending on external conditions. Taking into account (3), as well as the initial discharge temperature T at the moment of leaving the source, we can write
Tg = T (to ,t) = Ta + (Ti - Ta) e-Y(t-t0). (4)
Substituting (2) and (4) into (1), we get
mgZi (to,t) = —mgg + mg g (Ta + (Tj - Ta) e- fimgZi (t0,t).
Therefore
* (*>, t) = ( ^ - l) g + - /х* (to, t) •
Mg J* Mg Ta
Introduce the notation a = —¡i, b = (j^ — lj g, c = TTa'>9 ■ Taking into account that zi (t0,t) > 0, the solution to equation (1) can be written as
zi (to ,t) = max [Zi (to ,t) , 0] ,
where
, , c / b c \ fea(t-t0) - 1\
= + r + + {—ä—) +
,ce-T(i-io) b(t-to)
Y (a + Y) a
hi is the height of the i-th source mouth, and wi is the initial vertical velocity of the single emission when leaving the i-th source mouth.
Concentration ci (x, y, z, to,t) at point (x, y, z) at time t, created by the single emission with a center of mass at point (xi(to,t),yi(to,t),zi(to,t)) is described by a certain distribution characterizing the diffusion process of the single emission. In the ideal case, this process is described by the heat conduction equation or explicitly
/ (Xi(t0,t)-X)2 (yi(tQ,t)-y)2 (Zi(t0,t)-Z)2\
Ci(x,y,z,to,t) =--9--e v 2"x(t-t0) ^i(t-to) J}
(2n) ' ax(t - to)uy(t - to)az(t - to)
where ax (t — t0) ,ay (t — t0) ,az (t — t0) are functions that describe the expansion dynamics for the "cloud" of the single emission. To take into account the earth's surface, introduce a virtual source for each point source, which has coordinates (xi,yi, —h). Then the concentration distribution from a point source will be described as
Ci (x, y, z, to,t) = Ci (x, y, z, to, t) + Ci (x, y, -z, to,t).
The concentration generated by the source can be calculated as
t
Ci (x,y,z,t)= Ci (x,y,z,to,t) dto, (5)
where Ci (x,y,z,t) is the concentration created by the i-th source at the point (x,y,z) at time t. Then the pollutant concentration at a given point at time t is determined as follows
N
C (x, y, z, t) = ^^ Ci (x, y, z, t)
i=1
2. Numerical Implementation of the SUSUPLUME Air Pollution Model
It is assumed that data on meteorological conditions are received discretely with step tw. For the numerical integration algorithms to work, the values u (t) and v (t) are approximated by piecewise linear functions.
For dispersion functions ax (t — t0) ,oy (t — t0), oz (t — t0) we consider two sets of functions: the Pasquill-Briggs functions [23] and the MESOPUFF [24] functions, where the relation
x(t,to) = I \Ju(s)2 + v(s)2ds
'to
is used to estimate X.
—œ
Table
Dispersion functions for different models
Atmospheric stability category Pasquill-Briggs model MESOPUFF model
А ах (x) = cry (х) = 0, 32х (1 + 2^ö)~U'b, az(x) = 0,24x^1 + ^ Ox (x) = Oy (x) = 0,36ж°'9, az(x) = 0, 000236Ж2'1
В ох (х) = (Ту (х) = 0, 22х (1 + 2êô)~U'b, az(x) = 0,24x^1 + ^ ox (x) = Oy (x) = 0,25ж°'9, az(x) = О^бвж1'09
С <7* (f) = оу (х) = 0, 32ж (1 + ääö)"0'0, oz(x) = 0, 2х Ox (x) = Oy (x) = 0,19ж0'9, az(x) = 0,1 If0'91
D (7х (х) = Oy (х) = 0, Ш (1 + 2êô)~U'b, az(x) = 0,Ux(l + 1§rn)-°'5 Ox (¿0 = Oy (ж) = 0, 13ж0'9, az(x) = 0, 57ж0'58
Е (7х (х) = Oy (х) = 0, Ш (1 + U'b, Ox (x) = ay(x) = 0, 096ж°'9, az(x) = 0, 85ж0'47
F (Тх (х) = Oy (х) = 0, Их (1 + ' , oz(x) = 0,08ж (l + -f^)-0'5 ox (ж) = Oy (x) = 0, ОбЗж0'9, az(x) = 0, 77х0'42
In (5), the lower limit of integration is —to, but in practice, we are interested in calculating the concentration in a limited area near pollution sources. We assume that if the center of mass of the single emission is located far enough from the area of interest, then its contribution can be neglected. Suppose that the minimum wind speed for which the developed model can be applied is umin, and the size of the concentration calculation area is D, then in time tmax = the single emission will leave the concentration calculation area. Accordingly, (5) can be written as
Ci (x,y,z,t) = J Ci (x,y,z,to,t) dto. (6)
t_J2D_
t
To calculate (6) numerically, we use the modified Romberg's method [25]. The modification
is divided into several
t —
2D
t
consists in the fact that the segment of integration subsegments and the integral is calculated on each segment by the Romberg's method until the specified relative error e is reached. Let the integral be calculated on the segment [a; b]. The integration algorithm is presented below. Step 1. Put i = 0. Calculate the scores Io,o,h,o using the formula
Ii
b — a
i, 0
/ (a) + m + g /(a + fc^)
k= 1
0,1,
(7)
2
10-' 10"6
f 10"» : 10"1»
rS
N X -
4. \
л "♦v.
4000
■D
3500 §
Ol Ш
3000 Ç ш ■c
2500 £
2000 J
- 1500 I
- 1000 °
ш n
500 I z
0
5
10 15 20
Number of subsegments
25
■ Maximum Integration error
■ Average error of Integration
■ Ttie total number of calculations of the Integrand
Fig. 1. Integration error depending on the number of subsegments
100 50 0
-50
100 150 200 250
Pollutant source
Fig. 2. Pollutant concentration for a point source
-8
-100
-150
0
-200
Step 2. Put i = 1, j = 0. Calculate the scores I1;1 using the formula
I -1 1
/,, ' У"'- г = 0,1,..., j = 0,.., min (г, 5). (8)
Step 3. If the condition
\Ii,min(i ,5) Ii-1, j— 11 fi e \Ii,min(i ,5) | (9)
is satisfied, then the algorithm is completed. Otherwise, go to step 4. Step 4. Put i = i + 1, calculate the estimates Ii>0 by formula (7).
Step 5. Calculate the estimates Iij for all j = 0,..., min (i, 5) by formula (8). Go to Step 3.
Typically, the maximum permissible concentration of the main pollutants is at least 10—3 mg ■ m—3. To calculate the concentration of pollutants with an accuracy of 1%, the required accuracy for calculating the integral (6) is e = 10—8. Determine the number of segments into which it is necessary to split the segment of integration from the following computational experiment.
Let one source be 10m high, wind speed be 10m/s (Western), emission intensity be 1 g/s. Calculate the concentrations using formula (6) on a 5 * 5 m grid for x G [0; 1000] (the emission does not spread upwind), y G [0; 100] (the concentration will be symmetrical). The height of the concentration calculation is z = 2 meters. To calculate the exact integral, we use the Romberg's method on the entire segment t--with initial number of
points equal to 100, 000 and accuracy of calculation e = 10—12.
The experimental results (Fig. 1) show that it is necessary to split the segment
into 23 subsegments, which is 49,19 s per segment, to achieve the integration
accuracy 10—8. Romberg's method completed its work by calculating the objective function every 0, 3 s.
From a computational point of view, the model consists of two blocks. The first block registers the measured values such as wind direction and strength, etc. Since the concentration calculation requires data for a limited length, a circular buffer is used to
t-^t
Ю-4 i 10-S i 10"s i ю-'
I 10"s 10"» Ю-10 10-"
-200 -150 -100 -50
50 100 150 200
■ 5U5UPLUME
Fig. 3. Pollutant concentration along the X axis (y = 0)
Fig. 4. Pollutant concentration along the Y axis (x = 41,1)
accumulate information. This allows information accumulation operations to be performed in O(1) time. The second block calculates the pollutants concentration at the given points. Calculations of the pollutants concentrations at different points are independent and can be performed in parallel on different computers. The average time for calculating the concentration at one point from one source is about 10-4 s on an Intel Core 2 Duo T8300, 2,4 GHz processor for one subsegment.
3. Investigation of the SUSUPLUME Model
Let's consider a model situation. One point pollution source 10 m high is located at point (0;0). The source emits nitrogen dioxide (Mg = 46) with intensity equal to 1 g/s, temperature equal to 300°C and velocity equal to 3 m/s. The wind is considered to be constant with speed equal to 3 m/s, West. Consider atmospheric class - "D", air temperature 27°C. Let us choose the Pasquill-Briggs functions as the dispersion functions (Table). Parameter values are y = 0,1, ^ = 1, 0. Let's calculate the concentration (Fig. 2) at a height of 2 m for the area [0; 400] x [-200; 200]. The maximum concentration is
n-3
reached at the point (49,1; 0) and is equal to 1, 5 ■ 10 $ ■ m
Let us compare the obtained concentration values with the maximum calculated values calculated in accordance with [2]. Construct sections through the point (49,1; 0) along the X-axis (Fig. 3) and Y-axis (Fig. 4).
Now compare the concentration distribution calculated using the SUSUPLUME model and results of the official MRR-2017 method used in the Russian Federation. The maximum concentration according to the MRR-2017 technique (Fig. 5) is achieved at the point (41,1;0). Qualitatively, the SUSUPLUME model repeats the concentration dynamics, however, quantitatively it may differ from MRR-2017.
The calculated values of concentrations according to the model do not exceed the values specified in the normative document. In this case, the MRR-2017 model does not take into account the type of substance that is emitted. Let us consider how the concentration of pollutants for sulfur dioxide will change (Mg = 64).
о
х
X
- SUSUPLUME ---MRR
X
- SUSUPLUME (nitrogen dioxide) ■ ■■■■ SUSUPLUME (sulfur dioxide)
---MRR
Fig. 5. Comparison of the maximum concentration of pollutants in different models
Fig. 6. Pollutant concentration for nitrogen dioxide and sulfur dioxide
In this case, the maximum concentration value exceeds the value calculated by the MRR-2017 (Fig. 6). Let us consider the possibility of ensuring a decrease in the calculated concentration due to the correct selection of the 7, ^ parameters. Let's assume that 7, ^ G [0; 4, 6]. With 7 = ^ = 4, 6, the single blowout will decrease vertical speed by 99% in 1 second and decrease temperature by 99% in 1 second in relation to the environment. Consider such pairs (7; for which the inequality max (c (x, 0, 2) — cmrr(x, 0, 2)) < 0.
xe[0;200]
The allowed range of parameters is shown in Fig. 7.
0.5 1.0 1.5
2.5 3.0 3.5 4.0 4.5
0 25 50 75 100 125 150 175 200
Selected parameters
— 5U5UPLUME (nitrogen dioxide) ...... 5U5UPLUME (sulfur dioxide)
---MRR
Fig. 7. Values 7) at which the concentration value does not exceed the MRR-2017 model data
Fig. 8. Pollutant concentration for identified parameters
2.0
Let y = 0, 003714 and ^ = 1, 66310, then the calculated concentrations are shown in Fig. 8. The maximum concentration of sulfur dioxide fits into the calculation according to MPR-2017, however, with the given parameters, the concentration of nitrogen dioxide has dropped sharply.
Let us calculate the concentration of pollutants for the dispersion functions MESOPUFF (Table) for the initial conditions (Fig. 9, 10).
0.0035 0.0030
с
S 0.0025
§ 0.0020
с
о
S 0.0015
D
° 0.0010 0.0005 0.0000
5U5UPLUME (Pasqulll-Brlggs) 5U5UPLUME (MESOPUFF)
-200 -150 -100 -50
50 100 150 200
SUSUPLUME (Pasqulll-Brlggs) 5U5UPLUME (MESOPUFF)
D
X
Fig. 9. Concentration of pollutants for different dispersion functions (along the X axis)
Fig. 10. Concentration of pollutants for different dispersion functions (along the Y axis)
Using the MESOPUFF functions simulates a distribution that is much "squeezed" on the Y axis and more "smeared" on the X axis. The identification of the dispersion function , oy, az, in the presence of a large number of measurements, will improve the quality of modelling the air pollution.
When identifying model parameters on real data, it is necessary to identify the parameters y, and also take into account that the dynamics of the center of mass of the single emission is determined by the parameters of the gas-air mixture Mg, w, which are also subject to identification.
4. Advantages and Disadvantages of the Proposed Model
The constructed model makes it possible to simulate the dispersion of pollutants taking into account the real measurement of the wind over time. In addition, the model is able to take into account the type of emitted substance, the current season (through the ambient temperature), the possibility of parallelization when calculating the concentration. The disadvantage of the constructed model is the considerable time for calculating the concentration and the lack of taking into account the local movement of air masses.
Acknowledgment. This research was funded by Ministry of Science and Higher Education of the Russian Federation (grant No. FENU-2020-0022).
References
1. Krupnova T.G., Rakova O.V., Plaksina A.L., Gavrilkina S.V., Baranov E.O., Abramyan A.D. Short Communication: Effect of Urban Greening and Land Use on Air Pollution in Chelyabinsk, Russia. Biodiversitas, 2020, vol. 21, no. 6, pp. 2716-2720. DOI: 10.13057/Biodiv/D210646
2. Ob utverzhdenii metodov raschetov rasseivaniya vybrosov vrednykh (zagryaznyayushchikh) veshchestv v atmosfernom vozdukhe [On the Approval of Methods for Calculating the Dispersion of Emissions of Harmful (Polluting) Substances in the Air]. Moscow, Rostekhnadzor, 2017. (in Russian)
3. Tao Z., Santanello J.A., Chin M., Zhou S., Tan Q., Kemp E.M., Peters-Lidard C.D. Effect of Land Cover on Atmospheric Processes and Air Quality Over the Continental United States -A NASA Unified WRF (NU-WRF) Model Study. Atmospheric Chemistry and Physics, 2013, vol. 13, pp. 6207-6226. DOI: 10.5194/acpd-13-5429-2013
4. Kai Wang, Yang Zhang, Khairunnisa Yahya, Shiang-Yuh Wu, Grell G. Implementation and Initial Application of a New Chemistry-Aerosol Option in WRF/Chem for Simulation of Secondary Organic Aerosols and Aerosol Indirect Effects. Atmospheric Environment, 2015, vol. 115, pp. 716-732. DOI: 10.1016/j.atmosenv.2014.12.007
5. Khairunnisa Yahya, Glotfelty T., Kai Wang, Yang Zhang, Nenes A. Modeling Regional Air Quality and Climate: Improving Organic Aerosol and Aerosol Activation Processes in WRF/Chem Version 3.7.1. Geoscientific Model Development, 2017, vol. 10, pp. 2333-2363. DOI: 10.5194/gmd-10-2333-2017
6. Jian He, Ruoying He, Yang Zhang. Impacts of Air-Sea Interactions on Regional Air Quality Predictions Using a Coupled Atmosphere-Ocean Model in Southeastern U.S. Aerosol and Air Quality Research, 2018, vol. 18, pp. 1044-1067. DOI: 10.4209/aaqr.2016.12.0570
7. Solazzo E., Vardoulakis S., Cai X. A Novel Methodology for Interpreting Air Quality Measurements from Urban Streets Using CFD Modelling. Atmospheric Environment, 2011, vol. 45, pp. 5230-5239. DOI: 10.1016/j.atmosenv.2011.05.022
8. Gmez-Losada A., Santos F.M., Gibert K., Pires J.C.M. A Data Science Approach for Spatiotemporal Modelling of Low and Resident Air Pollution in Madrid (Spain): Implications for Epidemiological Studies. Computers, Environment and Urban Systems, 2019, vol. 75, pp. 1-11. DOI: 10.1016/J.Compenvurbsys.2018.12.005
9. Sanchez B., Santiago J.L., Martilli A., Martin F., Borge R., Quaassdorff C., De La Paz D. Modelling NOX Concentrations Through CFD-RANS in an Urban Hotspot Using High Resolution Traffic Emissionsa and Meteorology from a Mesoscale Model. Atmospheric Environment, 2017, vol. 163, pp. 155-165. DOI: 10.1016/j.atmosenv.2017.05.022
10. Kamenetsky E.S., Radionoff A.A. Aerodynamics of Mountain Valleys with Varying Cross Sections. Boundary-Layer Meteorology, 1999, vol. 91, no. 2, pp. 191-197.
11. Metodika rascheta kontsentratsiy v atmosfernom vozdukhe vrednykh veshchestv, soderzhashchikhsya v vybrosakh predpriyatiy (OND-86) [Methodology for Calculating the Concentration of Harmful Substances in the Atmospheric Air Contained in the Emissions of Enterprises (OND-86)]. Leningrad, Gidrometeoizdat, 1987. (in Russian)
12. Genikhovich E.L. Osnovnye napravleniya dorabotki normativnogo dokumenta OND-86 po raschetu rasseivaniya v atmosfere vybrosov zagryaznyayushchikh veshchestv [The Main Directions of the Revision of the Normative Document OND-86 on the Calculation of Dispersion of Pollutant Emissions in the Atmosphere]. Saint-Petersburg, NPK "Atmosphere", 2002. (in Russian)
13. Hurley P.J., Physick W.L., Luhar A.K. TAPM: A Practical Approach to Prognostic Meteorological and Air Pollution Modelling. Environmental Modelling and Software, 2005, vol. 20, no. 6, pp. 737-752.
14. Hurley P. TAPM V4. User Manual. CSIRO Marine and Atmospheric Research Internal Report, 2008, no. 5, 35 p.
15. Gifford F.A. Turbulent Diffusion Typing Schemes: a Review. Nuclear Safety, 1976, vol. 17, no. 1, pp. 68-86.
16. Atmospheric Dispersion Models for Application in Relation to Radionuclide Release. International Atomic Energy Agency, Technical Publications, vol. 379, 1986, 138 p.
17. Bespalov M.S. [Modelling the Spread of Contaminants in the Atmosphere as a Tool for Air Protection Activities]. Problems of Ecological Monitoring And Modeling of Ecosystems, 2016, vol. 27, no. 1, pp. 74-85. (in Russian)
18. Byzova N.L., Garger E.K., Ivanov V.N. Eksperimental'nye issledovaniya atmosfernoy diffuzii i raschety rasseyaniya primesi [Experimental Studies of Atmospheric Diffusion and Calculations of Impurity Scattering]. Leningrad, Gidrometeoizdat, 1991. (in Russian)
19. Byzova N.L., Ivanov V.N., Garger E.K. Turbulentnost' v pogranichnom sloe atmosfery [Turbulence in the Atmospheric Boundary Layer]. Leningrad, Goskomgidromet, 1989. (in Russian)
20. The AMS/EPA Regulatory Model - AERMOD. U.S. Environmental Protection Agency, EPA-454/B-03-001, 2004.
21. Scire J.S., Strimaitis D.G., Yamartino R.J. Model Formulation and User's Guide for the CALPUFF Dispersion Model. California Air Resources Board, 1990.
22. Brykin S.N., Markovsky V.V., Serebryakov I.S. Metodika razrabotki normativov predel'no-dopustimykh vybrosov radioaktivnykh veshchestv v atmosfernyy vozdukh [Methodology for the Development of Standards for Maximum Permissible Emissions of Radioactive Substances into the Atmospheric Air], 2010. (in Russian)
23. Berlyand M.E. Prognoz i regulirovanie zagryazneniy atmosfery [Forecast and Regulation of Atmospheric Pollution]. Leningrad, Gidrometeoizdat, 1985. (in Russian)
24. Scire J., Lurmann F., Bass A., Hanna S. User's Guide to the MESOPUFF II Model and Related Processor Programs. U.S. Environmental Protection Agency, Washington, EPA/600/8-84/013 (NTIS PB84181775), 1984.
25. Romberg W. Vereinfachte numerische integration. Det kongelige norske videnskabers selskab forhandlinger, Trondheim, 1955, vol. 28, no. 7, pp. 30-36. (in German)
Received September 5, 2020
УДК 51-74 DOI: 10.14529/mmp200401
ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ МОДЕЛИ SUSUPLUME РАСПРОСТРАНЕНИЯ ЗАГРЯЗНЯЮЩИХ ВЕЩЕСТВ В АТМОСФЕРЕ
С.М. Елсаков1, Д.А. Дрозин1, А.В. Геренштейн1, Т.Г. Крупнова1, С.Г. Ницкая1, Т.Ю. Оленчикова1, А.А. Замышляева1
1Южно-Уральский государственный университет, г. Челябинск, Российская Федерация
В настоящей работе предложена модель рассеивания загрязняющих веществ в атмосферном воздухе SUSUPLUME, представляющая собой модификацию классической
модели гауссового шлейфа (Gaussian plume model). Представленная модель описывает динамику рассеяния с учетом различных факторов: метеорологической обстановки, параметров источника загрязнения. Уравнения классической модели дополнены уравнениями движения центра масс единичного выброса. Численное исследование показало, что в стационарных метеоусловиях представленная модель качественно совпадет с другими известными моделями. Результаты расчета концентраций загрязняющих веществ не противоречат полученным значениям на основе методики ММР 2017, утвержденной к применению на территории Российской Федерации. Модель также содержит ряд идентифицируемых параметров, с помощью которых она может быть адаптирована к реальным условиям. Вычислительно модель состоит из двух блоков: блока регистрации измерительной информации и блока расчета концентраций загрязняющих веществ. Блок регистрации измерительной информации имеет низкую трудоемкость (свыше миллиона регистраций в секунду). Блок расчета концентраций загрязняющих веществ имеет высокую трудоемкость (около 400 точек вычислений концентрации в секунду), но поскольку все концентрации в представленной модели рассчитываются независимо, то это позволяет в перспективе использовать распараллеливание вычислительного процесса.
Ключевые слова: модель распространения загрязняющих веществ; гауссова модель; метод Ромберга.
Литература
1. Krupnova, T.G. Short Communication: Effect of Urban Greening and Land Use on Air Pollution in Chelyabinsk, Russia / T.G. Krupnova, O.V. Rakova, A.L. Plaksina, S.V. Gavrilkina, E.O. Baranov, A.D. Abramyan // Biodiversitas. - 2020. - V. 21, № 6. -P. 2716-2720.
2. Об утверждении методов расчетов рассеивания выбросов вредных (загрязняющих) веществ в атмосферном воздухе. - M.: Ростехнадзор, 2017.
3. Tao, Z. Effect of Land Cover on Atmospheric Processes and Air Quality Over the Continental United States - A NASA Unified WRF (NU-WRF) Model Study / Z. Tao, J.A. Santanello, M. Chin, S. Zhou, Q. Tan, E.M. Kemp, C.D. Peters-Lidard // Atmospheric Chemistry and Physics. - 2013. - V. 13. - P. 6207-6226.
4. Kai Wang. Implementation and Initial Application of a New Chemistry-Aerosol Option in WRF/Chem For Simulation of Secondary Organic Aerosols and Aerosol Indirect Effects / Kai Wang, Yang Zhang, Khairunnisa Yahya, Shiang-Yuh Wu, G. Grell // Atmospheric Environment. - 2015. - V. 115. - P. 716-732.
5. Khairunnisa Yahya. Modeling Regional Air Quality and Climate: Improving Organic Aerosol and Aerosol Activation Processes in WRF/ Chem Version 3.7.1 / Khairunnisa Yahya, Glotfelty T., Kai Wang, Yang Zhang, A. Nenes // Geoscientific Model Development. - 2017. -V. 10. - P. 2333-2363.
6. Jian He. Impacts Of Air-Sea Interactions on Regional Air Quality Predictions Using a Coupled Atmosphere-Ocean Model in Southeastern U.S. / Jian He, Ruoying He, Yang Zhang // Aerosol and Air Quality Research. - 2018. - V. 18. - P. 1044-1067.
7. Solazzo, E. A Novel Methodology for Interpreting Air Quality Measurements from Urban Streets Using CFD Modelling / E. Solazzo, S. Vardoulakis, X. Cai // Atmospheric Environment. - 2011. - V. 45. - P. 5230-5239.
8. Gmez-Losada, A. A Data Science Approach for Spatiotemporal Modelling of Low and Resident Air Pollution in Madrid (Spain): Implications for Epidemiological Studies / A. Gmez-Losada, F.M. Santos, K. Gibert, J.C.M. Pires // Computers, Environment and Urban Systems. - 2019. - V. 75. - P. 1-11.
9. Sanchez, B. Modelling NOX Concentrations Through CFD-RANS in an Urban Hotspot Using High Resolution Traffic EmissionsaAnd Meteorology from a Mesoscale Model / B. Sanchez, J.L. Santiago, A. Martilli, F. Martin, R. Borge, C. Quaassdorff, D. De La Paz // Atmospheric Environment. - 2017. - V. 163. - P. 155-165.
10. Kamenetsky, E.S. Aerodynamics of Mountain Valleys with Varying Cross Sections /
E.S. Kamenetsky, A.A. Radionoff // Boundary-Layer Meteorology. - 1999. - V. 91, № 2. -P. 191-197.
11. Методика расчета концентраций в атмосферном воздухе вредных веществ, содержащихся в выбросах предприятий ^НД-86). - Л.: Гидрометеоиздат, 1987.
12. Генихович, Е.Л. Основные направления доработки нормативного документа ОНД-86 по расчету рассеивания в атмосфере выбросов загрязняющих веществ / Е.Л. Генихович. -СПб.: НПК «Атмосфера» при ГГО им. А.И. Воейкова, 2002.
13. Hurley, P.J. TAPM: A Practical Approach to Prognostic Meteorological and Air Pollution Modelling / P.J. Hurley, W.L. Physick, A.K. Luhar // Environmental Modelling and Software. - 2005. - V. 20, № 6. - P. 737-752.
14. Hurley, P. TAPM V4. User Manual / P. Hurley // CSIRO Marine and Atmospheric Research Internal Report. - 2008. - № 5. - 35 p.
15. Gifford, F.A. Turbulent Diffusion Typing Schemes: a Review / F.A. Gifford // Nuclear Safety - 1976. - V. 17, № 1. - P. 68-86.
16. Atmospheric Dispersion Models for Application in Relation to Radionuclide Release // International Atomic Energy Agency, Technical Publications, V. 379. - 1986. - 138 p.
17. Беспалов, М.С. Моделирование распространения примеси в атмосфере как инструмент воздухоохранной деятельности / М.С. Беспалов // Проблемы экологического мониторинга и моделирования экосистем. - 2016. - Т. 27, № 1.- C. 74-85.
18. Бызова, Н.Л. Экспериментальные исследования атмосферной диффузии и расчеты рассеяния примеси / Н.Л. Бызова, Е.К. Гаргер, В.Н. Иванов. - Л.: Гидрометеоиздат, 1991.
19. Бызова, Н.Л. Турбулентность в пограничном слое атмосферы / Н.Л. Бызова, В.Н. Иванов, Е.К. Гаргер. - Л.: Госкомгидромет, 1989.
20. The AMS/EPA Regulatory Model - AERMOD. - U.S. Environmental Protection Agency, 2004. EPA-454/B-03-001.
21. Scire, J.S. Model Formulation and User's Guide for the CALPUFF Dispersion Model / J.S. Scire, D.G. Strimaitis, R.J. Yamartino // California Air Resources Board, 1990.
22. Брыкин, С.Н. Методика разработки нормативов предельно-допустимых выбросов радиоактивных веществ в атмосферный воздух. Т. 2 / С.Н. Брыкин, В.В. Марковский, И.С. Серебряков. - 2010.
23. Берлянд, М.Е. Прогноз и регулирование загрязнений атмосферы / М.Е.Берлянд. - Л.: Гидрометеоиздат, 1985.
24. Scire, J. User's Guide to the MESOPUFF II Model and Related Processor Programs / J. Scire,
F. Lurmann, A. Bass, S. Hanna. - U.S. Environmental Protection Agency, Washington, EPA/600/8-84/013 (NTIS PB84181775), 1984.
25. Romberg, W. Vereinfachte numerische Integration / W. Romberg // Det Kongelige Norske Videnskabers Selskab Forhandlinger, Trondheim, 1955. - V. 28, № 7. - P. 30-36. (in German)
Сергей Михайлович Елсаков, кандидат физико-математических наук, кафедра «Прикладная математика и программирование», Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), elsakovsm@susu.ru.
Дмитрий Александрович Дрозин, кандидат экономических наук, кафедра «Прикладная математика и программирование», Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), drozinda@susu.ru.
Аркадий Васильевич Геренштейн, кандидат физико-математических наук, доцент, кафедра «Прикладная математика и программирование», ЮжноУральский государственный университет (г. Челябинск, Российская Федерация), gerenshteinav@susu.ru.
Татьяна Георгиевна Крупнова, кандидат химических наук, доцент, кафедра «Экология и химическая технология:», Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), krupnovatg@susu.ru.
Светлана Георгивна Ницкая, кандидат технических наук, доцент, кафедра «Экология и химическая технология: , Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), nitckaiasg@susu.ru.
Татьяна Юрьевна Оленчикова, кандидат технических наук, кафедра «Прикладная математика и программирование: , Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), olenchikovati@susu.ru.
Алена Александровна Замышляева, доктор физико-математических наук, кафедра «Прикладная математика и программирование: , Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), zamyshliaevaaa@susu.ru.
Поступила в редакцию 5 сентября 2020 г.