CC BY
18
Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету залiзничного транспорту, 2017, № 3 (69)
ЕКОЛОГ1Я НА ТРАНСПОРТ!
UDC [502.3:504.5]:519.872
L. V. AMELINA1*, M. M. BILIAIEV2
1 Dep. «Hydraulics and Water Supply», Dnipropetrovsk National University of Railway Transport named after Academician V. Lazaryan, Lazaryan St., 2, Dnipro, Ukraine, 49010, tel. +38 (056) 273 15 09, e-mail gidravlika2013@mail.ru, ORCID 0000-0002-8525-7096
2Dep. «Hydraulics and Water Supply», Dnipropetrovsk National University of Railway Transport named after Academician V. Lazaryan, Lazaryan St., 2, Dnipro, Ukraine, 49010, tel. +38 (056) 273 15 09, e-mail gidravlika2013@mail.ru, ORCID 0000-0002-1531-7882
NUMERICAL SIMULATION OF AIR POLLUTION IN CASE OF UNPLANNED AMMONIA RELEASE
Purpose. Development fast calculating model which takes into account the meteorological parameters and buildings which are situated near the source of toxic chemical emission. Methodology. The developed model is based on the equation for potential flow and equation of pollutant dispersion. Equation of potential flow is used to compute wind pattern among buildings. To solve equation for potential flow Samarskii implicit difference scheme is used. The implicit change - triangle difference scheme is used to solve equation of mass transfer. Numerical integration is carried out using the rectangular difference grid. Method of porosity technique («markers method») is used to create the form of comprehensive computational region. Emission of ammonia is modeled using Delta function for point source. Findings. Developed 2D numerical model belongs to the class of «diagnostic models». This model takes into account the main physical factors affecting the process of dispersion of pollutants in the atmosphere. The model takes into account the influence of buildings on pollutant dispersion. On the basis of the developed numerical models a computational experiment was carried out to estimate the level of toxic chemical pollution in the case of unplanned ammonia release at ammonia pump station. Originality. Developed numerical model allows to calculate the 2D wind pattern among buildings and pollutant dispersion in the case unplanned ammonia release. Model allows to perform fast calculations of the atmosphere pollution. Practical value. The model can be used when developing the PLAS (Emergency Response Plan). Keywords: air pollution; unplanned release; toxic chemical; numerical modeling
Introduction
The ammonia pipeline Toliatti-Odessa was built in the late 70s-early 80s specifically for the transportation of the main products of the Toliatti nitrogen plant for export. The end point of ammonia pipeline is Odessa Port. There are several pumping stations along the route of this pipeline (Fig. 1). These pump stations support the correct pressure in ammonia pipeline. From the point view of industrial safety these pump stations are the chemically dangerous objects [4, 6, 15]. According to the Law of Ukraine for high-risk objects, a PLAS (Emergency Response Plan) document should be developed for such industrial
object. Prediction of contaminated zones and detection of Dangerous Level of Contamination is the basis of this document. Therefore, the actual task is to estimate the level of contamination in working areas of the pump station in the case of unplanned ammonia release.
Review of literature sources
To solve the problem of chemical contamination zones formation in the case of unplanned ammonia emissions analytical models are widely used. For example, Berland model was used to predict air pollution in the case of ammonia pipe rupture [6]:
Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету залiзничного транспорту, 2017, № 3 (69)
C =
exp
(zH j
2 J z
2(2 + n - m')k1^J%k0
о m-n i
y2 U1Z1 (
,2+n-m
H
2+n-m
4 kо x
k1 (2 + n - m )x
Ij-
2 + n - m
2+n-m
2uj z'm-nHz 2 (2 + n - m )2 kjx
Fig. 1. GOOGLE's image of Ammonia Pump Station:
1 - pump station
Mathematical formulation
To simulate the pollutant dispersion in the atmosphere 2D transport model is used [5, 7]
ÔC duC ôvC — +-+--
Ôt Ôx Ôy
C Ôf ÔC oC =—I цx—
Ôx I Ôx
ÔfV y ^ ) + £ Q (t )8(x - x )8(y - y ),
Ôy
(1)
QT - emission rate; H - height of emission source; k0, ku m - empirical constants; I - Bessel function; C-concentration.
Another approach for assessing the zones of chemical contamination is the application of the Gaussian plume model [2, 3, 10-14]. The use of the analytical models or Gaussian models allow to calculate quickly zones of chemical contamination. On the other hand, these models have significant lacks because they cannot be used when we model toxic chemical dispersion among buildings. For this purpose, it is necessary to use numerical models [1, 8, 9] which are based on Fluid Dynamics equations. In Ukraine, there is a certain deficit of such models [8, 9]. Worthy of note that the application Navier-Stokes equations for this purpose demands using of very fine computational grid and much computational time.
Purpose
The purpose of this paper is to develop a numerical model for computing the chemical contamination of air on the territory of the ammonia pump station for unplanned ammonia release (accidental release or terror act).
where C is mean concentration; u, v are the wind velocity components; c is the parameter taking into account the process of pollutant chemical decay or washout; , ) are the diffusion coefficients; Q is intensity of point source emission; 5 (r - ri) are
Dirak delta function; r= (xt, yt) are the coordinates of the point source .
To simulate the wind flow in the case of the buildings at the territory of Pump Station the 2D model of potential flow is used [7]
Ô2P Ô2P
Ôx2 Ôy2
= о,
(2)
where P is the potential of velocity. The wind velocity components are calculated as follows:
u =-
ÔP Ôx '
V = -
ÔP_
Ôy
Boundary conditions for modeling equations are discussed in [5, 7].
Numerical model
The computation of wind pattern and pollutant dispersion is carried out on rectangular grid. To create the form of buildings we use porosity technique or so called «markers method» [1, 7]. Markers are used to separate the computational cells where flow takes place from the cells which correspond to buildings.
Main features of the finite difference schemes which we use for the numerical integration of modeling equations are shown below.
To solve equation (1) we use change - triangle difference scheme [1, 7]. The time dependent derivative in Eq. (1) is approximated as follows:
m
X
Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету залiзничного транспорту, 2017, № 3 (69)
C „ en+1 - en dt ~
At
At the first step convective derivatives are represented in the following way:
duC du+C du C
dx dx
dx
dvC dv+C dv С
dy dy dy
where u+ = -
u + u
u — u
V + V
■; u =■
V =■
V— V
V = -
2
At the second step the convective derivatives are approximated as follows:
du+С м u++i,j Of — u+ С—11J dx
Ax
du С ^ ui+i, j См, j — uij С"+ d x
= L+x С
n+1
Ax
= L- С
n+1
dV+e j+iej—Vijeu—i r+en+i. dy Ay
ôV—e ^ Vi,j+iei, j+i— V—eij dy
Ay
=L+e
= L—yen+1.
The second order derivatives are approximated as follows:
d dx
d
dy
c
Vx dx
dC v ¥
^x-
Cn+1 — cn+1 cn+1 — cn+1 4+\,j 4j ~ 4,j 4—1, j
"Vx-
Ax2 ™ Ax2
=M—xCn+1 + M+xCn+1,
Cn+1_си+1
ьi,j 4,i—1
Cn+1 _C^+1 C^+1_C^+1
" " 1 J i,J
-M+1 ^y
written as follows:
С n+1 с n ij_ij_
At
+ L+Сn+1 + L—Cn+1 + L+ Cn+1 +
+L—y С n+l +-Cj =
= M+ Cn+1 + M - Cn+1 + M+Cn+1 + M-Cn+1
XX XX yy yy
Solution of the transport equation in finite -difference form is split in four steps on the time step of integration dt:
- at the first step ( k = 4 ) the difference equation is:
' n+k /~<n
(~< n+k _fn ij ij
At
+ \{L+X Ck + L+ Ck) + -С
k
4 -
= - (M+C n+k + M - Cn + M+C n+k + M-C" ) (3)
4 V xx » yy yy/V/
- at the second step ( k = n + 2; c = n + 4):the
difference equation is
Ck■ - Cc 1
4 +1L Сk + Ly Сk ) + - Ск, At 2V x y /4 IJ
--4(MxCk + M+Сc + My Сk + M+yCc) (4)
- at the third step (k = n + —• c = n + —) the
4' 2
expression (4) is used;
3
- at the fourth step (k = n +1; c = n + 4) the
expression (3) is used.
At the fifth step (at this step the influence of the source of pollutant emission is taken into account) the following approximation is used:
5 n+1 5 n
Ay2
V y
Ay2
Сi, j — Сi, j = N Ql (t" ) s
At
1=1 Ax Ay
=M~ Cn+1 + M+Cn+1
yy yy '
In these expressions
L+, L-, L+, Ly, M+x, M„, M+,, My are the difference operators. Using these expressions, the difference scheme for the transport equation can be
Function is equal to zero in all cells accept the cells where source of emission is situated.
This difference scheme is implicit and absolutely steady but the unknown concentration C is calculated using the explicit formulae at each step (so called «method of running calculation»).
2
2
2
Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету з^зничного транспорту, 2017, № 3 (69)
To solve equation (2) we transform it to the «evolution type»
ÔP Ô2P Ô2P
(3)
ôrç ôx2 dy2 ' where ^ is 'fictitious' time. For the solution of equation (3) tends to
the solution of equation (2).
To solve equation (3) A. A. Samarskii's change-triangle difference scheme is used. According to this scheme the solution of equation (3) is split into two steps:
- at the first step the difference equation is
pn+l/2 _pn pn _pn _™+]/2 + pn+1/2 pn _pn
'+ ' '' '' i_lj + •,j+/ i,j +
Ay2
_ pn+V2 + pn+/2
hi_', j_/
0,5Arç
i,j _ i+j,j i,j + i,j
Ax2
Ax2
+
Ay2
- at the second step the difference equation is
nn+1 - n> i.j i.j
0,5A^
+42
nn+1 - nn+1 . rl+1,j ri. j
' Ax2
V2
- Pf2 +
Ax2
pn+1 — pn+1 -Гi•.j+1 riJ
Ay2
+
-pn+1/2 + pn+l/2 ',j_',j-1
Ay2 •
From these expressions the unknown value Ptj is determined using the explicit formulae at each step of splitting («method of running calculation»).
The calculation is completed if the condition
|pn+1 - pn \ <s
Findings
Developed numerical model and code were used to compute ammonia concentrations at the territory of ammonia pump station in the case of unplanned release (Fig.2). It was supposed that release takes place near building with ammonia pumps (Fig. 3, 4). Sketch of computational region is shown in Figure 4. Figures 5, 6 show modeling results for ammonia emission. Emission rate is Q=17 kg/s and was chosen from literature [6].
Fig. 2. Ammonia pump station
Fig. 3. Buildings with ammonia pumps
is fulfilled (where s is a small number, n is the number of iteration). The components of velocity vector are calculated on the sides of computational cell as follows
ui,j =
=
P — P
ri, j ri-1,j
Ax
P — P
ri, j ri,j-1
Ay
Calculation of velocity components on the sides of computational cell allows to develop the conservative numerical scheme for pollutant dispersion.
For coding of difference formulae, we used FORTRAN language.
Fig. 4. Sketch of computational region (ammonia pump station): 1, 2 - buildings with ammonia pumps, 3, 4, 5 - industrial buildings on the territory of station; 6 - receptor position; 7 - position of ammonia release at the territory of pump station
Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету з^зничного транспорту, 2017, № 3 (69)
0.143Е-03
□
__0.215Е+03 о 0.201Е+03 г _a. 1G7E '93 it
0; Г.ТК '1! i _a. ISÖE+ЙЗ г.
ö. 14!iE+33 :
-0.132E+03 t
0,118E+fl3 ё
-а; ¡24Е+Ю
0.902E+02
-Й.7ВЗЕ+02 у
0;ё24Е+Я2
-a;405E+02
0.347E+02
-Ö.2f!ßä+82
0;694E+flt
0.3Б2Е*01 ____
coordinate x 8.24L.S+03
Fig. 5. Computed ammonia concentration, /=250 s
С. gT^ffl3
40 50 60 70 109,9 119,9
Fig. 6. Ammonia concentration at receptor spot (position No. 6, Fig. 4)
As we can see from Fig. 5 plume of toxic chemical quickly covers all the territory of ammonia pump station. To estimate the danger of such unplanned release we computed the ammonia concentration near one building at the territory of this station. Receptor position (position of person) is shown in Fig.4. Dynamics of ammonia concentration at this spot is shown in Fig.6. It is clear that the people affected will die at the territory in the case of this unplanned release because the ammonia concentration exceeds Level of Concern which is 20 mg/m3.
Worthy of note that computational time was 5 sec. It allows to use the developed model for prediction of air pollution during PLAS development.
Originality and practical value
A 2D numerical model has been developed to compute contamination zones among buildings during the accidental emission of a hazardous substance. The presented 2D numerical model is based on the application of the fundamental equations of aerodynamics and mass transfer.
The peculiarity of the developed model is the use of standard meteorological information and quick calculation.
Conclusions
Numerical 2D numerical model for estimating the level of atmospheric air pollution during the emergency emission of hazardous substances is proposed. Proposed numerical model allows to predict level of pollution of atmospheric air among buildings. The solution of the aerodynamic problem is based on the numerical integration of the equation for the velocity potential. To predict the air pollution, the equation of mass transfer is used. The mass transfer equation takes into account the convective and diffusive transport of pollutants in atmosphere, taking into account buildings situated near the source of emission. Emission of a dangerous substance is simulated by a point source, which is modeled using Dirac's delta function.
Further improvement of the model should be carried out in the direction of creating a 3D numerical model that takes into account the formation of vortices in the air flow.
LIST OF REFERENCE LINKS
1. Беляев, Н. Н. Моделирование нестационарных процессов аварийного загрязнения атмосферы : монография / Н. Н. Беляев, А. В. Берлов, П. Б. Машихина. - Днепропетровск : Акцент ПП, 2014. - 127 с.
2. Берлянд, М. Е. Прогноз и регулирование загрязнения атмосферы / М. Е. Берлянд. - Ленинград : Гидро-метеоиздат, 1985. - 273 с.
3. Бруяцкий, Е. В. Теория атмосферной диффузии радиоактивных выбросов / Е. В. Бруяцкий. - Киев : Ин-т гидромеханики НАН Украины, 2000. - 443 с.
4. Заказнов, В. Ф. Распространение аммиака при разгерметизации аммиакопровода, емкостей / В. Ф. Заказнов, Л. А. Куршева // Исследования и разработки по созданию магистральных аммиакопро-водов и складов жидкого аммиака : тр. ГИАП. - Москва, 1985. - С. 57.
5. Марчук, Г. И. Математическое моделирование в проблеме окружающей среды / Г. И. Марчук. - Москва : Наука, 1982. - 320 с.
Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету залiзничного транспорту, 2017, № 3 (69)
6. Цыкало, А. Л. Испарение и рассеивание аммиака при его разливах и утечках. Серия: Азотная промышленность / А. Л. Цыкало, И. И. Стрижевский, А. Д. Баглет. - Москва : НИИТЭХИМ, 1982. - 48 с.
7. Численное моделирование распространения загрязнения в окружающей среде / М. З. Згуровский, В. В. Скопецкий, В. К. Хрущ, Н. Н. Беляев. - Ки!в : Наук. думка, 1997. - 368 с.
8. Biliaiev, M. M. Numerical simulation of the atmosphere pollution after accident at the «Tolliaty-Odessa» ammonia pipe / M. M. Biliaiev, L. V. Amelina, M. M. Kharytonov // NATO Science for Peace and Security. Series C: Environmental security. - 2013. - P. 391-395. doi: 10.1007/978-94-007-5577-2_66.
9. Biliaiev, M. M. The Numeric Forecast of Air Pollution Caused by a Blasting Accident in the Enterprise Responsible for Rocket Fuel Utilization in Ukraine / M. M. Biliaiev, M. M. Kharytonov // NATO Science for Peace and Security. Series C: Environmental Security. - 2012. - P. 313-327. doi: 10.1007/978-94-007-5034-0_25.
10. Daly, A. Accident reconstruction and plume modeling of an unplanned ammonia release / A. Daly, P. Zanetti, M. Jennings // Air Pollution XX. WIT Transactions on Ecology and The Environment. - 2013. - Vol. 174. -P. 3-13. doi:10.2495/AIR130011.
11. Dispersion Modeling of Hydrogen Sulfide at Cimarex Rands Butte Project Using ALOHA / Bureau of Land Management Pinedale Field Office, SWCA Environmental Consultants. - Wyoming, 2010. - 26 p.
12. Janos, T. Atmostheric spreading model for ammonia released from the poultry house [Electronic resource] / T. Janos. E. Gorliczay, J. Borbely // Ecotoxicologie, Zootehnie si Tehnologii de Industrie Alimentara. -2016. - Vol. XV/B. - P. 331-337. - Available at: http://protmed.uoradea.ro/facultate/publicatii/ecotox_zooteh_ind_alim/2016B/ipa/17%20Tamas_Janos.pdf. -Title from the screen. - Accessed : 30.05.2017.
13. Mellsen, S. B. A Fortran Program for Calculating Chemical Hazards Using the NATO Stanag 2103/ATP-45 Algorithm [Electronic resource] / S. B. Mellsen // Suffield memorandum No. 1275. - Alberta : Defence Research Establishment Suffield, 1989. - 34 p. - Available at: http://oai.dtic.mil/oai/oai?verb=getRecord&metadataPrefix=html&identifier=ADA214763. - Title from the screen. - Accessed : 30.05.2017.
14. The analysis of the use of mathematical modeling for emergency planning purposes [Electronic resource] / O. Zavila, P. Dobes, J. Dlabka, J. Bitta // Bezpecnostni vyzkum. The science for population protection. -2015. - No. 2. - P. 1-9. - Available at: http://www.population-protection.eu/prilohy/casopis/30/213.pdf. - Title from the screen. - Accessed : 30.05.2017.
15. The Pentagon Shield Field Program: Toward Critical Infrastructure Protection / T. Warner, P. Benda, S. Swerdlin [et al.] // Bulletin of the American Meteorological Society. - 2007. - Vol. 88. - Iss. 2. -P. 167-176. doi: 10.1175/BAMS-88-2-167.
Л. В. АМЕЛ1НА1*, М. М. ыдлев2
1 Каф. «Пдравлжа та водопостачання», Дншропетровський нацюнальний ушверситет з^зничного транспорту iменi академжа В. Лазаряна, вул. Лазаряна, 2, Дншро, Украша, 49010, тел. +38 (056) 273 15 09, ел. пошта gidravlika2013@mail.ru, ORCID 0000-0002-8525-7096
2*Каф. «Пдравлжа та водопостачання», Дншропетровський нацюнальний утверситет залiзничного транспорту iменi академжа В. Лазаряна, вул. Лазаряна, 2, Дншро, Украша, 49010, тел. +38 (056) 273 15 09, ел. пошта gidravlika2013@mail.ru, ORCID 0000-0002-1531-7882
ЧИСЕЛЬНЕ МОДЕЛЮВАННЯ ЗАБРУДНЕННЯ ПОВ1ТРЯ В РАЗ1 РАПТОВОГО ВИКИДУ АМ1АКУ
Мета. Дослщження спрямоване на розробку моделi швидкого обчислення, яка враховувала б метеороло-пчш параметри та будiвлi, котрi знаходяться поблизу джерела токсичного викиду. Методика. Розроблена модель заснована на рiвняннi для потенцтного потоку та рiвняннi розсшвання забруднюючих речовин. Рiвняння потенцтного потоку використовуеться для обчислення моделi виру мiж будiвлями, а для вирь шення рiвняння потенцтно! течп - неявна рiзницева схема Самарського. Неявна поперемтно-трикутна рiз-ницева схема застосовуеться для виршення рiвняння масопереносу. Чисельне ттегрування здтснюеться за допомогою прямокутно! рiзницевоi сггки. Метод маркування («метод маркерiв») вживаеться для створення форми велико! розрахунково! области. Викид амiаку моделюеться з використанням Дельта функцп для точ-кового джерела. Результата. Розроблена двомiрна чисельна модель вщноситься до класу <^агностичш
Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету залiзничного транспорту, 2017, № 3 (69)
модел». Ця модель мае на yBa3i ochobhî фiзичнi фактори, що впливають на процес розсшвання забруднюю-чих речовин в атмосфера Модель враховуе вплив бyдiвель на розсiювaння забруднюючих речовин. На осно-Bi розроблених чисельних моделей був проведений обчислювальний експеримент для оцшки рiвня токсичного хiмiчного забруднення в рaзi раптового викиду ашаку на aмiaчнiй нaсоснiй стaнцiï. Наукова новизна. Розроблена авторами чисельна модель дозволяе розрахувати двомiрнy модель виру серед бyдiвель та розсшвання забруднюючих речовин у рaзi раптового викиду aмiaкy. Модель дозволяе ви-конувати швидк1 розрахунки забруднення атмосфери. Практична значимкть. Модель можливо використо-вувати при розробщ ПЛАСа (план лшшдацп aвaрiйних ситyaцiй).
Ключовi слова: забруднення повиря; раптовий викид; токсичнi хiмiчнi речовини; чисельне моделювання
Л. В. АМЕЛИНА1*, Н. Н. БЕЛЯЕВ2
1 Каф. «Гидравлика и водоснабжение», Днепропетровский национальный университет железнодорожного транспорта имени академика В. Лазаряна, ул. Лазаряна, 2, Днипро, Украина, 49010, тел. +38 (056) 273 15 09, эл. почта gidravlika2013@mail.ru, ORCID 0000-0002-8525-7096
2Каф. «Гидравлика и водоснабжение», Днепропетровский национальный университет железнодорожного транспорта имени академика В. Лазаряна, ул. Лазаряна, 2, Днипро, Украина, 49010, тел. +38 (056) 273 15 09, эл. почта gidravlika2013@mail.ru, ORCID 0000-0002-1531-7882
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ЗАГРЯЗНЕНИЯ ВОЗДУХА В СЛУЧАЕ ВНЕЗАПНОГО ВЫБРОСА АММИАКА
Цель. Исследование направлено на разработку быстрой расчетной модели, которая учитывала б метеорологические параметры и здания, расположенные вблизи источника токсического химического выброса. Методика. Разработанная модель основана на уравнении для потенциального потока и уравнении дисперсии загрязняющих веществ. Уравнение потенциального потока используется для вычисления модели ветра между зданиями, а для решения уравнения потенциального течения - неявная разностная схема Самарского. Неявная попеременно-треугольная разностная схема применяется для решения уравнения массопереноса. Численное интегрирование осуществляется с помощью прямоугольной разностной сетки. Метод маркировки («метод маркеров») используется для создания формы большой расчетной области. Эмиссия аммиака моделируется с использованием Дельта функции для точечного источника. Результаты. Разработанная двумерная численная модель относится к классу «диагностические модели». Эта модель имеет в виду основные физические факторы, влияющие на процесс рассеивания загрязняющих веществ в атмосфере. Модель учитывает влияние зданий на дисперсию загрязняющих веществ. На основе разработанных численных моделей был проведен вычислительный эксперимент для оценки уровня токсического химического загрязнения в случае внезапного выброса аммиака на аммиачной насосной станции. Научная новизна. Разработанная авторами численная модель позволяет рассчитать двумерную модель ветра между зданиями и дисперсию загрязняющих веществ в случае внезапного выброса аммиака. Модель позволяет выполнять быстрые расчеты загрязнения атмосферы. Практическая значимость. Модель можно применять при разработке ПЛАСа (план ликвидации аварийных ситуаций).
Ключевые слова: загрязнение воздуха; внезапный выброс; токсичные химические вещества; численное моделирование
REFERENCES
1. Biliaiev, M. M., Berlov, A. V., & Mashikhina, P. B. (2014). Modelirovaniye nestatsionarnykh protsessov avariynogo zagryazneniya atmosfery [Monograph]. Dnipropetrovsk: Aktsent PP.
2. Berlyand, M. Y. (1985). Prognoz i regulirovaniye zagryazneniya atmosfery. Leningrad: Gidrometeoizdat.
3. Bruyatskiy, Y. V. (2000). Teoriya atmosfernoy diffuzii radioaktivnykh vybrosov. Kyiv: Institut gidromekhaniki NAN Ukrainy.
4. Zakaznov, V. F., & Kursheva, L. A. (1985). Rasprostraneniye ammiaka pri razgermetizatsii ammiakoprovoda, emkostey. In Issledovaniya i razrabotki po sozdaniyu magistralnykh ammiakoprovodov i skladov zhidkogo ammiaka. Moscow: The State Research and Design Institute of the Nitric Industry and Organic Synthesis Products.
Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету з^зничного транспорту, 2017, № 3 (69)
5. Marchuk, G. I. (1982). Matematicheskoye modelirovaniye v probleme okruzhayushchey sredy. Moscow: Nau-ka.
6. Tsykalo, A. L., Strizhevskiy, I. I., & Baglet, A. D. (1982). Azotnaya promyshlennost: Ispareniye i rasseivaniye ammiakapri ego razlivakh i utechkakh. Moscow: NIITEKHIM.
7. Zgurovskiy, M. Z., Skopetskiy, V. V., Khrushch, V. K., & Biliaiev M. M. (1997). Chislennoye modelirovaniye rasprostraneniya zagryazneniya v okruzhayushchey srede. Kyiv: Naukova dumka.
8. Biliaiev, M. M., Amelina, L.V., & Kharitonov, M. M. (2013). Numerical simulation of the atmosphere pollution after accident at the "Tolliaty-Odessa" ammonia pipe. NATO Science for Peace and Security Series C: Environmental Security, 391-395. doi: 10.1007/978-94-007-5577-2_66
9. Biliaiev, M. M., & Kharytonov, M. M. (2012). The Numeric Forecast of Air Pollution Caused by a Blasting Accident in the Enterprise Responsible for Rocket Fuel Utilization in Ukraine. NATO Science for Peace and Security Series C: Environmental Security, 313-327. doi: 10.1007/978-94-007-5034-0_25
10. Daly, A., Zanetti, P., & Jennings, M. (2013). Accident reconstruction and plume modeling of an unplanned ammonia release. Air Pollution XXI. WIT Transactions on Ecology and the Environment, 174, 3-13. doi: 10.2495/AIR130011
11. SWCA Environmental Consultants. (2010). Dispersion Modeling of Hydrogen Sulfide at Cimarex Rands Butte Project Using ALOHA. Wyoming.
12. Janos, T., Gorliczay, E., & Borbely, J. (2016). Atmostheric spreading model for ammonia released from the poultry house. Ecotoxicologie, Zootehnie si Tehnologii de Industrie Alimentara, XV/B, 331-337. Retrieved from
http://protmed.uoradea.ro/facultate/publicatii/ecotox_zooteh_ind_alim/2016B/ipa/17%20Tamas_Janos.pdf
13. Mellsen, S. B. (1989). A Fortran Program for Calculating Chemical Hazards Using the NATO Stanag 2103/ATP-45 Algorithm: Suffield memorandum 1275. Alberta: Defence Research Establishment Suffield. Retrieved from http://oai.dtic.mil/oai/oai?verb=getRecord&metadataPrefix=html&identifier=ADA214763
14. Zavila, O., Dobes, P., Dlabka, J., & Bitta, J. (2015). The analysis of the use of mathematical modeling for emergency planning purposes. Bezpecnostni vyzkum, 2. Retrieved from http://www.population-protection. eu/prilohy/casopis/30/213.pdf
15. Warner, T., Benda, P., Swerdlin, S., Knievel, J., Copeland, J., Crook, A., ..., & Weil, J. (2007). The Pentagon Shield Field Program: Toward Critical Infrastructure Protection. Bulletin of the American Meteorological Society, 88 (2), 167-176. doi: 10.1175/BAMS-88-2-167
Prof. S. A. Pichugov, Dr. Sc. in Phys.-and-Math. (Ukraine); Prof. S. Z. Polishchuk, D. Sc. (Tech.),
(Ukraine) recommended this article to be published
Accessed: Feb. 10, 2017
Received: May 18, 2017
