Atmosphere protection in case of emergency during transportation of dangerous cargo

Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету залiзничного транспорту, 2016, № 1

UDC 502.36/.504.3.054:656.2

O. V. BERLOV1*

1 Dep. «Hydraulics and Water Supply», Dnipropetrovsk National University of Railway Transport named after Academician V. Lazaryan, Lazaryan St., 2, Dnipropetrovsk, Ukraine, 49010, tel. +38 (056) 273 15 09, e-mail berlov@bigmir.net, ORCID 0000-0002-7442-0548

ATMOSPHERE PROTECTION IN CASE OF EMERGENCY DURING TRANSPORTATION OF DANGEROUS CARGO

Purpose. The paper highlights the development of numerical models for prediction of atmospheric pollution in case of burning of the solid rocket propellant in a railway car, situated near the building on railway territory. These models can be used in predicting the effectiveness of neutralization upon the atmosphere protection for this type of accidents. Methodology. To solve this problem the numerical models based on the use of Navier-Stokes equations, to determine the velocity field of the wind flow near cars and buildings, and contaminants-transfer equations in the atmosphere were developed. For the numerical integration of pollutant transport equation was used implicit «change - triangle» difference scheme. When constructing a difference scheme physical and geometric cleavage of the transfer equation is carried out in four steps. Unknown value of pollutant concentration at each step of cleavage is determined by the explicit scheme - the method of «point-to-point computation». For the numerical integration of the Navier-Stokes equations are used implicit difference schemes. When carrying out computing experiment also takes into account: the velocity profile of wind flow; interaction between the building and the wind flow and flame jet of solid rocket propellant; the presence of a railroad car; inside which there is a source of pollution; instability of pollutant emissions. On the basis of constructed numerical models was performed the computer experiment for assessing the level of air pollution at dangerous cargo rail transportation in case of emergency at railway territory.The application calculations for the timely combustion products neutralization of solid rocket propellant were carried out. Findings. The numerical models that let promptly calculate air contamination in case of emergency during solid rocket propellant transportation, as well as calculate the rational parameters of pollutant neutralization process were developed by the researcher. These models can be used for routine calculations of various accident scenarios simulation. Originality. Numerical models were developed; they take into account significant factors, influencing the pollutant dispersion process in the atmosphere. On their base a pollutant neutralization method was offered in emergency situations on the railway transport. Practical value. Efficient numerical models, so called «diagnostic models» were considered for the rapid calculation of the air pollution level and air protection technology in emergency situations, in particular, in the case of railway transportation the solid rocket propellant.

Keywords: train accidents; air pollution; pollutant neutralization; numerical simulation

Introduction

Accidents during transportation of toxic chemicals, for example, solid rocket propellant (Fig.1) can intensively impact on the environment [2-4, 9].

In the case of this propellant burning the great amount of toxic chemicals are dispersed in atmosphere. So it is important to predict the atmosphere pollution in the case of the possible accident and develop methods to protect atmosphere from pollution (Fig. 2) [1, 11].

That is why the prediction of the atmosphere pollution in the case of the mighty emissions during transportation is of great interest.

Fig. 1. Transport of solid rocket propellant in a railway wagon

Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету залiзничного транспорту, 2016, № 1

da "dt

dua dva

dx dy

J_ Re

д 2a

д V

dx dy

d2y d2y ~öxT ~öyr

= -a.

Fig. 2. General view of settlement areas: 1 - building on the territory of station «Pavlohrad-1», 2 - train

Nowadays to predict the dimensions of the hitting area in the case of outdoor toxic chemical release the special standard model is used in Ukraine [10]. This model is formed on the basis of some empirical models and it has a lot of lacks and is, without doubt, unrealistic. The model does not take into account the influence of the wind velocity and the atmosphere diffusion on the concentration dispersion in the atmosphere. The main lack of this model is that the standard model cannot calculate the change of toxic chemical concentration in the atmosphere after the accident. In some cases Gauss model is used to predict atmosphere pollution level [5-7, 13]. It is very important to develop CFD (Computational Fluid Dynamics) models to solve the problem discussed [1, 8, 14, 15].

Purpose

The purpose of this work is the development of the numerical model to predict the atmosphere pollution in micro scale level in the case of solid pro-pellant burning which also can be used to predict the effectiveness of neutralization for atmosphere protection for this type of accident.

Methodology

When the solid propellant burns different toxic chemicals are emitted in the atmosphere. To compute the toxic chemical (product of propellant burning) dispersion the Navier-Stokes equations and equation of mass transfer are used [14]:

where y - is flow function; ra = — - — - is vor-

dx dy

ticity; u = —, v = —.

dy dx

This model is used to predict the air pollution near railway car with solid rocket propellant (Fig.2).

To simulate the process of toxic chemical dispersion in the atmosphere (for the distance about 300m) the 2-D transport equation is used [8, 14]:

C

dt

duC dvC

dx dy

- cC = div (цgradC )-

h^ôi

where C - is the concentration of toxic chemical; u, v - are the velocity components in x, y direction respectively; ^ = x, ^ y ) - are the coefficients of

turbulent diffusion in x, y direction respectively; c - is the parameter taking into account the process of toxic gas decay or rain wash out; xi, yi - are the

coordinates of point source of emission; Qi (t) - is the intensity of pollutant emission; S(x - xi )ô(y - yi) - is Dirac delta-function.

The boundary conditions for the governing equations are discussed in [14].

In the developed numerical model, the following profile of velocity component u and coefficient of diffusion is used [5, 6]:

u = u.

( Y

z.

yi

ц y = k1

r y > yi

, Цx = k0u

where u1 - is the velocity at height y1 ; k1 = 0,2; k0 = 0,1; p = 0,16; m « 1.

Numerical integration. To solve the governing equations the implicit difference schemes were used [14, 15]. These are so called «change - triangle» difference schemes.

i=1

Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету з^зничного транспорту, 2016, № 1

Findings

To minimize the extent of contamination and the risk of toxic defeat people in case of emergency during transportation of solid rocket propellant it is proposed technology, aimed at localization of the source of emission of products of combustion solid rocket propellant. It is proposed to include in the train which transports solid rocket propellant, the car with the pump unit and tank with a stock of neutralizing solution. It is also proposed to include a buffer car between the rail car where solid propellant is transported and the car with the pump unit. This buffer truck will serve as a «screen» for the protection of the car with the pump unit from heat radiation of the combustion solid rocket propellant. With this technology, was organized rapid feed of neutralizing solution into the plume of combustion products and thereby minimizing time the free flow of products of combustion solid rocket propellant in the atmosphere. In addition, unlike the traditional scheme of feeding catalyst in the form of streams offered feed through the sprinkler, which is posted on pull-rod (Fig.3). This is because the temperature of the combustion solid rocket propellant in the atmosphere is high and if you use a traditional presentation in the form of a jet, there will be rapid evaporation of droplets of reagent on the path of its movement from the installation to the jet plume.

Fig. 3. Sketch of toxic chemical neutralization:

1 - vessel with the reagent, 2 - pipeline, 3 - pump, 4 - pop-arrow, 5 - buffer car, 6 - sprinkler

Problem setting- to evaluate the effectiveness of the process of neutralization of various intensities of reagent feed and placement of the sprinkler to the source of emission of hazardous substances.

The calculation is performed for the following tasks:

Scenario 1 - no neutralization;

Scenario 2 - supply of reagent is performed at a height of 13,5m and at a distance of 20 m from the source of emission. This distance is constant in all scenarios, catalyst feed rate of 1 kg/s;

Scenario 3 - the reagent feed is done at the height of 13,5m, catalyst feed rate of 8 kg/s.

l.4e7Et0Z о I.456E+0Z r I.4Z4E+0Z d I.393E+0Z i 1.361E+02 n I.330E+0Z a I.Z99E.0Z t I.Z67E+0Z e I.Z36E.0Z I.Z04E+0Z 1.173E+0Z у 1.141E+0Z 1.110E+0Z 1.7еВЕ»01 I.471E+01 1.157E+01

0.500E+00

0.150E+01

_LL

coordinate к

Fig. 4. Concentration contours of HCl (no reagent supply, t = 120 sec): 1 - car with solid rocket propellant, 2 - building near railway station

On the base of the developed numerical models the special codes were developed. Some results of their application are shown in Fig.4-6.

The subsequent figures illustrate the area of atmospheric pollution in the case of the options considered. The first figure (Fig. 4) corresponds to the situation when there is no flow of reagent and compares the efficiency of neutralization for each option, comparing to the area of contamination with this picture.

Fig. 5. Concentration contours of HCl (reagent supply, t = 120 sec, Scenario 2: the height of the feed catalyst 13,5 m, flow rate 1 kg/s): 1 - sprinkler; 2 - building

Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету з^зничного транспорту, 2016, № 1

Fig. 6. Concentration contours of HCl (reagent supply, t = 120 sec, Scenario 3: the height of the feed catalyst 13,5 m, flow rate 8 kg/s): 1 - sprinkler; 2 - building

Results of neutralization effectiveness are shown in Tabl.1, 2.

Table 1

The efficiency of neutralization

t 45 сек 75 сек 120 сек

The efficiency of neutralization The efficiency 6,8% of neutra 8,2% ilization

t 45 сек 75 сек 120 сек

The efficiency of neutralization 50% 62% 69%

Results from Tabl.1, 2 show that the rate of reagent strongly influences the effectiveness of protection. Also it is clear from Fig.5, 6 that the height of reagent supply is too high because the region near the railway car and building is still in the plume of toxic chemical.

Originality and practical value

A new numerical model to predict atmosphere pollution and neutralization of toxic chemical after accidents at railways was developed. The model is based on the 2-D transport equation. The developed model takes into account wind profile, diffusion, emission rate. This model allows evaluating

the level of atmosphere pollution and effectiveness of neutralization. The model can be useful in the field of transport routes safety prediction.

Conclusions

The article contains results of numerical simulation of air pollution near «Pavlohrad-1» station after the accident at railway. To simulate the process of air pollution the 2-D developed numerical model was used. The developed numerical model takes into account the main physical processes which influences the pollutant dispersion in atmosphere. The developed models allow simulating the process of neutralization of toxic chemical in atmosphere. The future work in this field will be connected with development of fluid dynamics models to simulate the pollutant dispersion over the complex terrain.

LIST OF REFERENCE LINKS

1. Беляев, Н. Н. Математическое моделирование в задачах экологической безопасности и мониторинга чрезвычайных ситуаций : монография / Н. Н. Беляев, Е. Ю. Гунько, П. Б. Ма-шихина. - Днепропетровск : Акцент ПП, 2013.

- 159 с.

2. Беляев, Н. Н. Прогнозирование загрязнения приземного слоя атмосферы при горении твердого ракетного топлива в хранилище / Н. Н. Беляев, А. В. Берлов // Зб. наук. пр. НГУ.

- Дншропетровськ, 2013. - № 42. - С. 160-167.

3. Беляев, Н. Н. Моделирование аварийного загрязнения атмосферы при чрезвычайной ситуации в хранилище твердого ракетного топлива / Н. Н. Беляев, А. В. Берлов, А. В. Шевченко // Наука та прогрес транспорту. - 2014. -№ 5 (53). - С. 29-38. doi: 10.15802/stp-2014/29973.

4. Берлов, А. В. Расчет загрязнения приземного слоя атмосферы при горении твердого ракетного топлива / А. В. Берлов // Наук. Вкн. буд-ва : зб. наук. пр. / Харк. нац. ун-т буд-ва та ар-хгг. - Харшв, 2014. - №1 (75). - С. 185-189.

5. Берлянд, М. Е. Прогноз и регулирование загрязнения атмосферы / М. Е. Берлянд. - Ленинград : Гидрометеоиздат, 1985. - 273 с.

6. Бруяцкий, Е. В. Теория атмосферной диффузии радиоактивных выбросов / Е. В. Бруяцкий. -Киев : Институт гидромеханики НАН Украины, 2000. - 443 с.

7. Гусев, Н. Г. Радиоактивные выбросы в биосфере / Н. Г. Гусев, В. А. Беляев. - Москва : Энер-гоатомиздат, 1991. - 257 с.

Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету залiзничного транспорту, 2016, № 1

8. Марчук, Г. И. Математическое моделирование в проблеме окружающей среды / Г. И. Марчук. - Москва : Наука, 1982. - 320 с.

9. Машихина, П. Б. Моделирование распространения примеси в атмосфере с учетом рельефа местности / П. Б. Машихина // Вюн. Дшпро-петр. нац. ун-ту залiзн. трансп. iм. акад. В. Ла-заряна. - Днiпропетровськ, 2009. - Вип. 27. -С. 138-142.

10. Методика прогнозування наслщшв виливу (ви-киду) небезпечних хiмiчних речовин при ава-рiях на промислових об'ектах i транспортi. -Кшв, 2001. - 33 с.

11. Моделирование аварийных ситуаций на промышленных объектах и безопасность жизнедеятельности / Л. И. Антошкина, Н. Н. Беляев, Л. Ф. Долина, Е. Д. Коренюк. - Дшпропет-ровськ : Нова вдеолопя, 2011. - 123 с.

А. В. БЕРЛОВ1*

1 Каф. «Гидравлика и водоснабжение», Днепропетровский национальный университет железнодорожного транспорта имени академика В. Лазаряна, ул. Лазаряна, 2, Днепропетровск, Украина, 49010, тел. +38 (056) 273 15 09, эл. почта berlov@bigmir.net, ORCID 0000-0002-7442-0548

ЗАЩИТА АТМОСФЕРЫ В СЛУЧАЕ АВАРИИ ПРИ ТРАНСПОРТИРОВКЕ ОПАСНЫХ ГРУЗОВ

Цель. Работа посвящена разработке численных моделей для прогноза загрязнения атмосферы при горении твердого ракетного топлива в железнодорожном вагоне, расположенном вблизи здания на примагистральной территории. Данные модели могут быть использованы при прогнозировании эффективности нейтрализации по защите атмосферы для этого типа аварий. Методика. Для решения поставленной задачи разработаны численные модели, основанные на применении уравнений Навье-Стокса, для определения поля скорости ветрового потока вблизи вагонов и зданий, и уравнении переноса примеси в атмосфере. Для численного интегрирования уравнения транспорта загрязнителя использовалась неявная попеременно-треугольная разностная схема. При построении разностной схемы осуществляется физическое и геометрическое расщепление уравнения переноса на четыре шага. Неизвестное значение концентрации загрязнителя на каждом шаге расщепления определяется по явной схеме - методу «бегущего счета». Для численного интегрирования уравнений Навье-Стокса применяются неявные разностные схемы. При проведении вычислительного эксперимента также учитываются: профиль скорости ветрового потока; взаимодействие ветрового потока со зданием и струей продуктов горения твердого ракетного топлива; наличие железнодорожного вагона, внутри которого находится источник загрязнения; нестационарность эмиссии загрязнителя. На основе построенных численных моделей проведен вычислительный эксперимент по оценке уровня загрязнения атмосферы при транспортировке железнодорожным транспортом опасного груза в случае аварии на примагистральной территории. Проведены расчеты по применению своевременной нейтрализации продуктов горения твердого ракетного топлива. Результаты. Исследователем разработаны численные модели, которые позволяют оперативно рассчитать загрязнение атмосферы в случае аварии при транспортировке твердого ракетного топлива, а также рациональные параметры процесса нейтрализации загрязнителя. Эти модели могут использоваться для серийных расчетов по моделированию различных сценариев аварийных ситуаций. Научная новизна. Разработаны численные модели, позволяющие учесть существенные факторы, влияющие на процесс рассеивания загрязняющих веществ в атмосфере и предложена на их базе методика нейтрализации загрязнителя при чрезвычайных ситуациях на железнодорожном транспорте. Практическая значимость. Рассмотрены эффективные численные модели «diagnostic models» для экспресс-расчета уровня загрязнения атмосферы и разработана технология защиты

12. Самарский, А. А. Теория разностных схем / А. А. Самарский. - Москва : Наука, 1983. -616 с.

13. Уорк, К. Загрязнение воздуха. Источники и контроль / К. Уорк, С. Уорнер - Москва : Мир, 1980. - 539 с.

14. Численное моделирование распространения загрязнения в окружающей среде / М. З. Згуров-ский, В. В. Скопецкий, В. К. Хрущ, Н. Н. Беляев. - Кшв : Наук. думка, 1997. - 368 с.

15. Biliaiev, M. Numerical Simulation of Indoor Air Pollution and Atmosphere Pollution for Regions Having Complex Topography / M. Biliaiev // NATO Science for Peace and Security Series C: Environmental Security. - 2012. - P. 87-91. doi: 10.1007/978-94-007-1359-8 15.

Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету залiзничного транспорту, 2016, № 1

атмосферы при аварийных ситуациях, в частности, в случае перевозки железнодорожным транспортом твердого ракетного топлива.

Ключевые слова: аварии на железнодорожном транспорте; загрязнение атмосферы; нейтрализация загрязнителя; численное моделирование

0. В. БЕРЛОВ1*

1 Каф. «Пдравлжа та водопостачання», Дншропетровський нацюнальний ушверситет зал1зничного транспорту 1меш академжа В. Лазаряна, вул. Лазаряна 2, Дтпропетровськ, Украша, 49010, тел. +38 (056) 273 15 09, ел. пошта berlov@bigmir.net, ORCID 0000-0002-7442-0548

ЗАХИСТ АТМОСФЕРИ У ВИПАДКУ АВАРП

ПРИ ТРАНСПОРТУВАНН1 НЕБЕЗПЕЧНИХ ВАНТАЖ1В

Мета. Робота присвячена розробщ чисельних моделей для прогнозу забруднення атмосфери при горiннi твердого ракетного палива в залiзничному вагонi, розташованому поблизу будiвлi на примагiстральнiй територи. Данi моделi можуть бути використаш при прогнозуваннi ефективностi нейтралiзацil по захисту атмосфери для цього типу аварш. Методика. Для розв'язання поставлено! задачi розробленi чисельнi моделi, засноват на застосуваннi рiвнянь Нав'е-Стокса, для визначення поля швидкостi вiтрового потоку поблизу вагошв i будiвель, та рiвняннi переносу домiшки в атмосфера Для чисельного iнтегрування рiвняння транспорту забруднювача використовувалася неявна поперемiнно-трикутна рiзницева схема. При побудовi рiзницевоl схеми здшснюеться фiзичне i геометричне розщеплення рiвняння переносу на чотири кроки. Неввдоме значения концентрацй' забруднювача на кожному крощ розщеплення визначаеться за явною схемою - методу <^жучого рахунку». Для чисельного iнтегрування рiвнянь Нав'е-Стокса застосовуються неявнi рiзницевi схеми. При проведены обчислювального експерименту також враховуються: профiль швидкостi вiтрового потоку; взаемодiя вiтрового потоку з будiвлею i струменем продуктiв горiння твердого ракетного палива; наявнють залiзничного вагона, усередиш якого знаходиться джерело забруднення; нестацюнаршсть емюп забруднювача. На основi побудованих чисельних моделей проведено обчислювальний експеримент iз оцiнки рiвня забруднення атмосфери при транспортуванш залiзничним транспортом небезпечного вантажу при аварп на примапстральнш територи. Проведено розрахунки по застосуванню своечасно! нейтралiзацil продуктiв горшня твердого ракетного палива. Результати. Дослщником розроблено чисельнi моделi, яш дозволяють оперативно розрахувати забруднення атмосфери в разi аварй' при транспортуваннi твердого ракетного палива, а також рацюнальш параметри нейтралiзацil забруднювача. Цi моделi можуть використовуватись для серiйних розрахуншв рiзноманiтних сценарпв аварiйних ситуацiй. Наукова новизна. Розроблено чисельш моделi, що дозволяють врахувати iстотнi фактори, що впливають на процес розсiювання забруднюючих речовин у атмосферi та запропонована на 1'х базi методика нейтралiзацil забруднювача при надзвичайних ситуациях на залiзничному транспортi. Практична значимiсть. Розглянуто ефективнi чисельнi моделi «diagnostic models» для експрес-розрахунку рiвня забруднення атмосфери та розроблено технолопю захисту атмосфери при аваршних ситуациях, зокрема, у разi перевезення залiзничним транспортом твердого ракетного палива.

Ключовi слова: аварй на залiзничному транспорту забруднення атмосфери; нейтралiзацiя забруднювача; чисельне моделювання

REFERENCES

1. Belyayev N.N., Gunko Ye.Yu., Mashikhina P.B. Matematicheskoye modelirovaniye v zadachakh eko-logicheskoy bezopasnosti i monitoringa chrezvychaynykh situatsiy [Mathematical modeling in problems of environmental safety and monitoring emergencies]. Dnepropetrovsk, Aktsent PP Publ., 2013, 159 p.

2. Belyayev N.N., Berlov A.V. Prognozirovaniye zagryazneniya prizemnogo sloya atmosfery pri gorenii tver-dogo raketnogo topliva v khranilishche [Prediction of pollution in the atmospheric surface layer during combustion of solid rocket fuel in storage]. Zbirnyk naukovykh prats NHU [Proc of NMU]. Dnipropetrovsk, Natsionalnyi Hirnychyi Universytet Publ., 2013, no. 42, pp. 160-167

3. Belyayev N.N., Berlov A.V., Shevchenko A.V. Modelirovaniye avariynogo zagryazneniya atmosfery pri chrezvychaynoy situatsii v khranilishche tverdogo raketnogo topliva [Simulation of the accidental pollution af-

HayKa Ta nporpec TpaHcnopTy. BicHHK ^mnponeTpoBctKoro Ha^oH&ntHoro ymBepcureTy 3&ni3HHHHoro TpaHcnopTy, 2016, № 1

ter the emergency in to the storage of rocket solid propellant]. Nauka ta prohres transport - Science and Transport Progress, 2014, no. 5 (53), pp. 29-38. doi: 10.15802/stp2014/29973.

4. Berlov A.V. Raschet zagryazneniya prizemnogo sloya atmosfery pri gorenii tverdogo raketnogo topliva [Calculation of the pollution of the ground atmospheric layer during combustion of solid rocket fuel]. Naukovyi visnyk budivnytstva: zbirnyk naukovykh prats [Proc. of «Scientific Bulletin of construction»], 2014, no. 1 (75), pp.185-189.

5. Berlyand M.Ye. Prognoz i regulirovaniye zagryazneniya atmosfery [Prediction and regulation of air pollution]. Leningrad, Gidrometeoizdat Publ., 1985. 273 p.

6. Bruyatskiy Ye.V. Teoriya atmosfernoy diffuzii radioaktivnykh vybrosov [The theory of atmospheric diffusion of radioactive emissions]. Kiev, Institut gidromekhaniki NAN Ukrainy Publ., 2000. 443 p.

7. Gusev N.G., Belyaev V.A. Radioaktivnyye vybrosy v biosfere [Radioactive emissions in the biosphere]. Moscow, Energoatomizdat Publ., 1991. 257 p.

8. Marchuk G.I. Matematicheskoye modelirovaniye v probleme okruzhayushchey sredy [Mathematical modeling in the environmental problem]. Moscow, Nauka Publ., 1982. 320 p..

9. Mashikhina P.B. Modelirovaniye rasprostraneniya primesi v atmosfere s uchetom relefa mestnosti [The distribution modeling of impurities in the atmosphere with taking into account of terrain]. Visnyk Dnipropet-rovskoho natsionalnoho universytetu zaliznychnoho transportu imeni akademika V. Lazariana [Bulletin of Dnipropetrovsk National University of Railway Transport named after Academician V. Lazaryan], 2009, issue 27, pp. 138-142.

10. Metodyka prohnozuvannia naslidkiv vylyvu (vykydu) nebezpechnykh khimichnykh rechovyn pry avariiakh na promyslovykh obiektakh i transporti [Methods of prediction the consequences of the spout (emission) of hazardous chemicals in accidents at industrial objects and transport]. Kyiv, 2001. 33 p.

11. Antoshkina L.I., Belyaev N.N., Dolina L.F., Korenyuk Ye.D. Modelirovaniye avariynykh situatsiy na pro-myshlennykh obektakh i bezopasnost zhiznedeyatelnosti [Simulation of emergency situations at industrial facilities and life safety]. Dnipropetrovsk, Nova ideolohiia Publ., 2011. 123 p.

12. Samarskiy A.A. Teoriya raznostnykh skhem [The theory of difference schemes]. Moscow, Nauka Publ., 1983. 616 p.

13. Uork, K., Uorner S. Zagryazneniye vozdukha. Istochniki i kontrol [Air pollution. Sources and control]. Moscow, Mir Publ., 1980. 539 p.

14. Zgurovskiy M.Z., Skopetskiy V.V., Khrushch V.K., Belyaev N.N. Chislennoye modelirovaniye rasprostra-neniya zagryazneniya v okruzhayushchey srede [Numerical modelling of pollution spreading in the environment]. Kyiv, Naukova dumka Publ., 1997. 368 p.

15. Biliaiev M. Numerical Simulation of Indoor Air Pollution and Atmosphere Pollution for Regions Having Complex Topography. NATO Science for Peace and Security Series C: Environmental Security, 2012, pp. 87-91. doi: 10.1007/978-94-007-1359-8_15.

Prof. S. A. Pichugov, Dr. Sc.

recommended this article to be published

Accessed: Nov. 16, 2015

Received: Jan. 11, 2016