Numerical simulation of toxic chemical dispersion after accident at railway Текст научной статьи по специальности «Математика»

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

ЕКОЛОГ1Я НА ТРАНСПОРТ!

UDC 502.36/.504.3.054:656.2

M. M. BILIAIEV1*, L. YA. MUNTIAN2*

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 [email protected], ORCID 0000-0002-1531-7882

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 [email protected], ORCID 0000-0002-1659-7740

NUMERICAL SIMULATION OF TOXIC CHEMICAL DISPERSION AFTER ACCIDENT AT RAILWAY

Purpose. This research focuses on the development of an applied numerical model to calculate the dynamics of atmospheric pollution in the emission of dangerous chemical substances in the event of transportation by railway. Methodology. For the numerical simulation of transport process of the dangerous chemical substance in the atmosphere the equation of convection-diffusion pollutant transport is used. This equation takes into account the effect of wind, atmospheric diffusion, the power of emission source, as well as the movement of the source of emission (depressurized tank) on the process of pollutant dispersion. When carrying out computing experiment one also takes into account the profile of the speed of the wind flow. For the numerical integration of pollutant transport in the atmosphere implicit finite-difference splitting scheme is used. The numerical calculation is divided into four steps of splitting and at each step of splitting the unknown value of the concentration of hazardous substance is determined by the explicit running account scheme. On the basis of the numerical model it was created the code using the algorithmic language FORTRAN. One conducted the computational experiments to assess the level of air pollution near the railway station «Illarionovo» in the event of a possible accident during transportation of ammonia. Findings. The proposed model allows you to quickly calculate the air pollution after the emission of chemically hazardous substance, taking into account the motion of the emission source. The model makes it possible to determine the size of the land surface pollution zones and the amount of pollutants deposited on a specific area. Using the developed numerical model it was estimated the environmental damage near the railway station «Illarionovo». Originality. One can use the numerical model to calculate the size and intensity of the chemical contamination zones after accidents on transport. Practical value. The numerical model, developed by authors, can be used to estimate the size and intensity of the chemical contamination zones during emergency on transport. The developed numerical model solves the problem of assessing the impact of emergency emission of ammonia near the railway station «Illarionovo».

Keywords: railway transport accidents; air pollution; chemically dangerous substance; numerical simulation

Introduction

Many accidents at railways result in toxic chemicals emissions [1]. These emissions pose potential risk to human health and environment. To predict the damage after these emissions Government Instructions are used. These Instructions are

based on simple empirical formulae which do not take account many physical factors, such as, wind profile, diffusion, etc. For example, to predict the width of polluted area the following expression is used in these Instructions

D = 0,3 • Ln,

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

where L is the length of polluted area (the L value is determined using the special Table), n is parameter which depends on the atmosphere condition (stable, neutral, etc.).

To obtain quick information about the atmosphere pollution Gauss model is used [5-7, 13, 16] or empirical model OHfl-86 [14].

In some cases CFD models are used [2-4, 9, 10, 12, 15] but these models are comprehensive for the regular users.

For a quick evaluation of toxic chemical concentrations in atmosphere fast hazard tools are very in need. It is very important for practice to have tools which take into account important details and, on the other hand, are not time consuming.

Purpose

The purpose of this work is the development of numerical model to predict the atmosphere pollution after accidents at railways.

Methodology

Mathematical model of pollutant dispersion. To simulate the toxic chemical dispersion in atmosphere 3-D transport equation is used [2, 8, 15]

- at the entrance boundary:

dC duC dvC d(w - ws)C -+-+-+ ^-s!—

dt dx dy dz

dC

= — I M^-1 + —

dx V dx J dy

dC

ц y

-kC + cC =

_d_( dC_ dz V z dz

+X(t)5(x - x,. (t))s(y - y, (t))S(z - z, ), (1)

where C is air concentration of toxic chemical; u, v, w are the wind components in the x, y and z directions; ws is gravitational settling velocity; c is the chemical decay coefficient; k is the precipitation scavenging coefficient; ^ = , , ) are

the eddy diffusivities for the three coordinate directions; qt (t) is source terms;

ô(x - x, (t))ô(y - y, (t))S(z - z, ) are Dirac delta

function; (x,,y,, z, ) are the coordinates of the

point source.

This equation is numerically integrated using the following boundary conditions:

C|

boundary

= Св

where Ce is known concentration.

At the exit boundary (for example, this is the plane x = Lx, where Lx is distance from the entrance plane) the boundary condition is used in the following difference form

C (, +1, j, k ) = C ( j, k ),

where (, +1, j, k ) is the last computational cell and

(,, j, k ) is the previous computational cell. At the

ground (z=0 plane) the boundary condition is

dC C

-= a C,

dn

where n is normal to this plane, a> 0 is coefficient which takes into account toxic chemical interaction with ground.

The initial condition (at time t = 0 ) can be written as C = 0 in the computational region or C = C0,, (where C0 is the known concentration at

the site where the instant emission took place and C = 0 in the other part of the computational region.

In the developed numerical model the following approximations for wind speed and diffusion coefficient are used:

u = u

lg z/z

1 lg zJz

■, Ц y = k0 • U , ^z = ki

( z >

V zi J

where u1 is wind speed at height zx = 10m ; m «1; k, k0 are parameters [Berlyand et al].

Numerical model. To solve transport equation (1) the implicit change - triangle scheme is used. The main features of this scheme are shown below.

According to the principles of this implicit scheme development the time dependent derivative is approximated as following:

dC Cjl - qk

dt

At

At the first step of development the convective derivatives are represented as follows:

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

duC du+С du С

dx dx dx

dVC dV+С dV-C

dy dy dy

dwC dw+С dw~ С \--

Ci, j,k - Ci, j-\,k

Ay2

= MyCn+\ + МУуС

n+\

where u + =-

dz dz dz u + lui

2

u =

u -\u\ + V + V

V =-—

2

V =

V - V + w + w w =■

w - w

w =■

dx Ax

du С u- cn+\ - u- cn+\ ^ "i+\, j,k i+\, j,k ui, j,k i, j,k

dx Ax

dv+С ^ Vi,j+\,kCi, j,k - Vi, j,kCi, j-\,k

dy Ay

dv-C ^ Vi, j+\,kCi, j+\,k - Vi, j,kCi, j,k

dy Ay

dw+С ^ w+j,k+\Ci,j,k - wtj,kCi,j,k-\

dz Az

dw~C ^ wi,j,k+\ci,j,k+\ - wi,j,kCi,j,k

= LZC

n+\

n+\

dz

= L+C

= L-yC"

= L+Cn+\

Az

= L-C

n +\

dx dx ) Ax2

x

Ci, j,k - Ci-\

j,k

Ax2

-in+\

=M-xxCn+i + m+C" ,

f d_c}

dy y y dy ;

/~m+\ _ n+\

Ci, j+!,k Ci, j,k

Ay2

д f dc ) cjk+\ - cjk Tz Г z К-aZ2-

s-in+\ s-in+\

ci, j,k - ci, j,k-\

Az2

= M-Cn+\ + M+Cn+\.

2 2 2 At the second step the convective derivatives are approximated as following:

du+C ^ u+1+\,;,kC"j,k - Kj,kC"-\,j,k = LXcn+1

In these expressions L+x, Lx, L+, Ly, L+, Lz,

M+, Mx are the difference operators. Using these expressions the difference scheme for the transport equation can be written as follows:

cn+1 en

i,i,k - + L+xCn+1 + L-xCn+1 + L+,Cn+1 +

At

y

The second order derivatives are approximated as following:

<1 / 'ysi N en+1 r>n+1

d I dC ) Ci+1, j,k - Ci, j,k

. T- s~in+\ . T+ Г^n+\ . T-Г^n+\ . _ s~in+\

+LyC + LZC + LzC Ci, ],k =

= (M+xCn+l + MxCn+1 + M+yCn+l +

+M-yyCn+l + M+zCn+l + M-zCn+l) + q (t )8.

Solution of this equation is split in four steps on the time step of integration dt :

- at the first step (k = -4) the difference equation is:

j - j + 1(L+cCk + L+Ck + L+Ck) + -Ck =

At 2 x y z 4 yk

= \(M+xCk + MxCk + M+yCn + MyCn +

4

M+zzcn + M-Zcn )d,,

zz ' / < . , '

1=\ 4

- at the second step ( k = n + \ , c = n + \ ): the

2 4

difference equation is

y^k л

cijk- cijk \

At 2

+ -( L-Ck + L-Ck + LzCk ) + - Ck =

~ \ x y z > . ijk

С — I

4

2

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

= -(MxC + M+xC + MyC + M+yC +

+M-C + M+ZZCC) + ь1,

l=1 4

3 1

- at the third step (k = n + —, c = n + —) the

4 2

difference equation is

s~ik л

Cijk- Cyk 1

At

+- (L+Ck + z;;ck + zxk)+-Ck =

V X y z > . i]k

1

= T (MXxCс + M+xCk + MyCk + M+yCC +

+M-zCk + M+zCC) + 8i,

l=1 4

- at the fourth step (k = n +1, c = n + —) the

difference equation is

s~ik л

Ci]k- Ci]k 1

At

+2( + LyC + ¿ZCk)C]k =

1

k-1 = С

C (x, y, z, tn ),

k = 2,3, 4

1 k 5

where C , C , C - are concentrations at each time step.

FORTRAN language was used to code the developed numerical model.

Findings

The developed generic code «EMISSION.» was used to solve the following problem. A train with toxic chemical (NH3) moves near Railway Station «Illarionovo» (Dniepropetrovsk Region, Ukraine) and at time t=0 the instant emission of NH3 takes place. This emission results in NH3 cloud formation (Fig.2). Position of this instant emission is schematically shown as «star» in Fig.1 and «arrow» indicates the direction of the train movement. The train keeps moving after the accident and long term emission of NH3 (q = 5 kg /sec) follows the instant emission. So we have scenario «instant exhaust of toxic chemical + long term emission of it».

=-(MxxC + MXxC + MyC + M+yC +

+M-zCc + M+Ck) + 5l,

l=1 4

Function is equal to zero in all cells accept the cells where the ' l ' 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»).

The initial condition at each time step is written as following [11]:

Fig. 1. Satellite image near Railway Station «Illarionovo» : 1 - Site # 1 (farm field); 2 - Site #2 (farm field); dots show the rail way track; «star» shows the initial point of emission

Results of numerical simulation of atmosphere pollution for this scenario are shown in Fig. 2-4.

C (x, y, z, tn+1 ) = C

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

Fig. 2. Spatial distribution of NH3, (initial cloud), t = 42sec (level z = 10m)

Fig. 3. Spatial distribution of NH3, t = 70sec (level z = 10m )

Fig. 4. Spatial distribution of NH3, t =169sec (level z = 10m)

It is clear from Fig. 2-4 that the contamination area is formed under the following factors: initial cloud appearance + long term emission from the moving train + atmosphere dispersion and wind. From the practical point of view it is important

to evaluate the ground pollution intensity. To solve this problem the following expression was integrated

T

G = |dt U (ws + p,a)C(x,y,z = 0,t)ds , (2)

0 S, z=0

where G is mass of the pollutant whish felt down the ground site S, t is time, T is period of integration. «Rectangular method» was used for the numerical integration of this expression. Results of soil pollution prediction after accident are shown in Tabl.1.

Table 1

Time 112 sec 169 sec 239 sec

Area 4628*103 m2 6222*103 m2 6925*103 m2

Soi! Pollution Dynamics. In the case of private owners demands to pay their losses after the accident it is important to know the toxic chemical felt down the different sites, for example, the private fields. At present to evaluate the ecological damage in Ukraine, in the case of ground pollution the following formulae is used

D = K x n , (3)

where n is polluted area [m2], K is money coefficient (grivni, dollars) declared by Administration Regulations.

As it is clear from Eq. (3) to make the adequate evaluation of the ecological damage D it is necessary to obtain the correct information about n. The developed generic code can solve this problem using formulae (2) which is calculated for marked cells of the computational region. These cells indicate the area, for example, of the farm fields (Fig.2, Site #1and Site # 2). In Tabl.2 results about these sites pollution are presented.

Table 2

Mass of Toxic Chemical felt down the Different Sites

Time Site #1 (See Fig.1) Site #2 (See Fig.1)

154 sec 0.41 x 105 grams 723 grams

211 sec 0.81 x 105 grams 6 166 grams

267 sec 0.10 x 106 grams 0.22 x 105 grams

Results from Tabl.2 show that the most contaminated area is formed at Site # 1 which is influenced by initial cloud and following emission from moving train.

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

Originality and practical value

A new numerical model to predict atmosphere pollution after accidents at railways was developed. The model is based on the 3-D transport equation. The developed model takes into account wind profile, diffusion, emission rate and source movement This model allows to evaluate the ecological damage in the case of different emissions at railways. The model can be useful in the field of safety prediction of transport routes and risk assessment.

Conclusions

The article contains results of numerical simulation of air pollution near Illarionovo station after accident at railway. To simulate the process of air pollution the 3-D developed numerical model was used. The developed numerical model takes into account the main physical processes which influence the pollutant dispersion in atmosphere. The future work in this field will be connected with development of fluid dynamics model to simulate the pollutant dispersion over the complex terrain.

LIST OF REFERENCE LINKS

1. Авария и катастрофы. Предупреждение и ликвидация последствий : учеб. пособие : в 5 кн. / под ред. В. А. Котляревского, А. В. Забегаева.

- Москва : ACB, 2001. - Кн. 5. - 416 с.

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

- Днепропетровск : Акцент ПП, 2013. - 159 с.

3. Беляев, Н. Н. Прогнозирование загрязнения приземного слоя атмосферы при горении твердого ракетного топлива в хранилище / Н. Н. Беляев, А. В. Берлов // Зб. наук. пр. НГУ. - Днш-ропетровськ, 2013. - № 42. - С. 160-167.

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

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

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

7. Израэль, Ю. А. Экология и контроль состояния природной среды / Ю. А. Израэль. - Москва : Гидрометеоиздат, 1984. - 560 с.

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

- Москва : Наука, 1982. - 320 с.

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

10. Русакова, Т. И. Прогнозирование загрязнения воздушной среды от автотранспорта на улицах и в микрорайонах города / Т. И. Русакова // Наука та прогрес транспорту. - 2013. - № 6 (48). - С. 32-45. doi: 10.15802/stp2013/19675.

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

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

- Киев : Наук. думка, 1997. - 368 с.

13. Швыряев, А. А. Оценка риска воздействия загрязнения атмосферы в исследуемом регионе / А. А. Швыряев, В. В. Меньшиков. - Москва : Изд-во МГУ, 2004. - 124 с.

14. Эльтерман, В. М. Вентиляция химических производств / В. М. Эльтерман. - Москва : Химия, 1980. - 288 с.

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. - Springer : Netherlands, 2012. - P. 87-91. doi: 10.1007/978-94-007-1359-8_15.

16. Hanna, S. Air Quality Modelling over. Short Distances / S. Hanna // College on Atmospheric Boundary Layer and Air Pollution Modelling. - 16 May-3 June, 1994. - № SMR/760-2. -P. 712-743.

Наука та прогрес транспорту. Вкник Дншропетровського нащонального ушверситету залiзничного транспорту, 2016, № 2 (62)

М. М. БШЯеВ1*, Л. Я. МУНТЯН2*

1 Каф. «Пдравлжа та водопостачання», Днiпропетровський нацюнальний унгверситет залiзничного транспорту iменi академжа В. Лазаряна, вул. Лазаряна, 2, Днгпропетровськ, Укра!на, 49010, тел. +38 (056) 273 15 09, ел. пошта [email protected], ОЯСГО 0000-0002-1531-7882

2*Каф. «Гiдравлiка та водопостачання», Днiпропетровський нацюнальний унгверситет залiзничного транспорту iменi академжа В. Лазаряна, вул. Лазаряна, 2, Днгпропетровськ, Укра1на, 49010, тел. +38 (056) 273 15 09, ел. пошта [email protected], ОЯСГО 0000-0002-1659-7740

ЧИСЕЛЬНЕ МОДЕЛЮВАННЯ РОЗПОВСЮДЖЕННЯ Х1М1ЧНО НЕБЕЗПЕЧНО1 РЕЧОВИНИ ПРИ АВАРП НА ЗАЛ1ЗНИЦ1

Мета. Дане дослвдження спрямоване на розробку прикладно! чисельно! моделi для розрахунку динамiки забруднення атмосфери при емюи хiмiчно небезпечно! речовини у разi його транспортування залiзничним транспортом. Методика. Для чисельного моделювання процесу перенесення хiмiчно небезпечно! речовини в атмосферi застосовуеться рiвняння конвективно-дифузiйного переносу домшки. Дане рiвняння враховуе вплив виру, атмосферно! дифузii, потужностi джерела емюи, а також рух джерела емюи (розгерметизовано! цистерни) на процес розсшвання забруднювача. При проведеннi обчислювального експерименту також враховуеться профiль швидкостi вирового потоку. Для чисельного iнтегрування рiвняння переносу домiшки в атмосферi застосовуеться неявна рiзницева схема розщеплення. Чисельний розрахунок розбиваеться на чотири кроки розщеплення, i на кожному крощ розщеплення невiдоме значення концентрацii небезпечно! речовини визначаеться за явною схемою «бiжучого» рахунку. На основi побудовано! чисельно! моделi створений код на алгоршмчнш мовi Фортран. Проведено обчислювальний експеримент iз оцiнки рiвня забруднення атмосфери поблизу залiзничноi' станцii «1ларюнове» у разi можливо! аварii при транспортуванш амiаку. Результата. Запропонована модель дозволяе оперативно розрахувати забруднення атмосфери при викидi хiмiчно небезпечно! речовини з урахуванням руху джерела емюи. Модель дозволяе визначити розмiри зони забруднення поверхш землi, а також шлькють забруднювача, що осiв на конкретну дiлянку. За допомогою розроблено! чисельно! моделi виконана оцiнка еколопчного збитку поблизу залiзничноi станцii «1ларюнове». Наукова новизна. Створена ефективна чисельна модель, що дозволяе оперативно розраховувати еколопчний збиток при емюи небезпечних речовин на об'ектах залiзничного транспорту. Практична значимкть. Розроблена авторами чисельна модель може бути використана для оцшки розмiрiв та iнтенсивностi зон хiмiчного забруднення при надзвичайних ситуациях на транспортi. За допомогою розроблено! чисельно! моделi розв'язана задача з оцшки наслщшв аварiйноi емiсii ашаку бiля залiзничноi станцii «Iларiонове».

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

Н. Н. БЕЛЯЕВ1*, Л. Я. МУНТЯН2*

1 Каф. «Гидравлика и водоснабжение», Днепропетровский национальный университет железнодорожного транспорта имени академика В. Лазаряна, ул. Лазаряна, 2, Днепропетровск, Украина, 49010, тел. +38 (056) 273 15 09, эл. почта [email protected], ОЯСГО 0000-0002-1531-7882

2*Каф. «Гидравлика и водоснабжение», Днепропетровский национальный университет железнодорожного транспорта имени академика В. Лазаряна, ул. Лазаряна, 2, Днепропетровск, Украина, 49010, тел. +38 (056) 273 15 09, эл. почта [email protected], ОЯСГО 0000-0002-1659-7740

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ РАСПРОСТРАНЕНИЯ ХИМИЧЕСКИ ОПАСНОГО ВЕЩЕСТВА ПРИ АВАРИИ НА ЖЕЛЕЗНОЙ ДОРОГЕ

Цель. Данное исследование направлено на разработку прикладной численной модели для расчета динамики загрязнения атмосферы при эмиссии химически опасного вещества в случае его транспортировки железнодорожным транспортом. Методика. Для численного моделирования процесса переноса химически опасного вещества в атмосфере применяется уравнение конвективно-диффузионного переноса примеси.

Наука та прогрес транспорту. Вкник Дншропетровського нащонального ушверситету залiзничного транспорту, 2016, № 2 (62)

Данное уравнение учитывает влияние ветра, атмосферной диффузии, мощности источника эмиссии, а также движение источника эмиссии (разгерметизированной цистерны) на процесс рассеивания загрязнителя. При проведении вычислительного эксперимента также учитывается профиль скорости ветрового потока. Для численного интегрирования уравнения переноса примеси в атмосфере применяется неявная разностная схема расщепления. Численный расчет разбивается на четыре шага расщепления, и на каждом шаге расщепления неизвестное значение концентрации опасного вещества определяется по явной схеме «бегущего» счета. На основе построенной численной модели создан код на алгоритмическом языке Фортран. Проведен вычислительный эксперимент по оценке уровня загрязнения атмосферы вблизи железнодорожной станции «Илларионово» в случае возможной аварии при транспортировке аммиака. Результаты. Предложенная модель позволяет оперативно рассчитать загрязнение атмосферы при выбросе химически опасного вещества с учетом движения источника эмиссии. Модель позволяет определить размеры зоны загрязнения поверхности земли, а также количество загрязнителя, осевшего на конкретный участок. С помощью разработанной численной модели выполнена оценка экологического ущерба вблизи железнодорожной станции «Илларионово». Научная новизна. Создана эффективная численная модель, позволяющая оперативно рассчитывать экологический ущерб при эмиссии опасных веществ на объектах железнодорожного транспорта. Практическая значимость. Разработанная учеными численная модель может быть использована для оценки размеров и интенсивности зон химического загрязнения при чрезвычайных ситуациях на транспорте. С помощью разработанной численной модели решена задача по оценке последствий аварийной эмиссии аммиака возле железнодорожной станции «Илларионово».

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

REFERENCES

1. Kotlyarevskiy V.A., Zabegayev A.V. Avariya i katastrofy. Preduprezhdeniye i likvidatsiya posledstviy. Kn. 5. [Accident and disasters. Warning and liquidation of consequences, Vol. 5]. Mockow, ACB Publ., 2001.

2. Belyayev N.N., Gunko Ye.Yu., Mashikhina P.B. Matematicheskoye modelirovaniye v zadachakh eko-logicheskoy bezopasnosti i monitoringa chrezvychaynykh situatsiy [Mathematical design in the tasks of ecological safety and monitoring of emergencies]. Dnipropetrovsk, Aktsent PP Publ., 2013. 159 p.

3. Belyayev N.N., Berlov A.V. Prognozirovaniye zagryazneniya prizemnogo sloya atmosfery pri gorenii tver-dogo raketnogo topliva v khranilishche [Prognostication of contamination of the ground layer of atmosphere at burning of solid rocket propellant in depository]. Zbirnyk naukovykh prats natsionalnoho hirnychoho univer-sytetu [Proc. of the National Mining University], 2013, no. 42, pp. 160-167.

4. Berlov A.V. Raschet zagryazneniya prizemnogo sloya atmosfery pri gorenii tverdogo raketnogo topliva [Zbirnyk naukovykh prats natsionalnoho hirnychoho universytetu «Naukovyi visnyk budivnytstva» [Proc. of the National Mining University «Scientific Bulletin of Building»], 2014, no. 1 (75), pp. 185-189.

5. Berlyand M.Ye. Prognoz i regulirovaniye zagryazneniya atmosfery [Prognosis and adjusting of contamination of atmosphere]. Leningrad, Gidrometeoizdat Publ., 1985. 273 p.

6. Gusev N.G., Belyayev V.A. Radioaktivnyye vybrosy v biosfere [Radioactive discharges in the biosphere]. Mockow, Energoatomizdat Publ., 1991. 257 p.

7. Izrayel Yu.A. Ekologiya i kontrol sostoyaniya prirodnoy sredy [Ecology and control of the state of natural environment]. Moscow, Gidrometeoizdat Publ., 1984. 560 p.

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

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

10. Rusakova T.I. Prognozirovaniye zagryazneniya vozdushnoy sredy ot avtotransporta na ulitsakh i v mikrorayonakh goroda [Prediction of air pollution from autotransport in streets and districts of city]. Nauka ta prohres transportu - Science and Transport Progress, 2013, no. 6 (48), pp. 32-45. doi: 10.15802/stp2013/19675.

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

Наука та прогрес транспорту. Вкник Дншропетровського нащонального ушверситету залiзничного транспорту, 2016, № 2 (62)

12. Zgurovskiy M.Z., Skopetskiy V.V., Khrushch V.K., Belyayev N.N. Chislennoye modelirovaniye rasprostraneniya zagryazneniya v okruzhayushchey srede [Numerical modeling of pollution in the environment]. Kyiv, Naukova dumka Publ., 1997. 368 p.

13. Shvyryayev A.A., Menshikov V.V. Otsenka riska vozdeystviya zagryazneniya atmosfery v issleduyemom re-gione [Risk assessment of atmosphere pollution effect in the region under study]. Moscow, MGU Publ., 2004. 124 p.

14. Elterman V.M. Ventilyatsiya khimicheskikh proizvodstv [Ventilation of chemical productions]. Moscow, Khimiya Publ., 1980. 288 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.

16. Hanna S. Air Quality Modelling over. Short Distances. College on Atmospheric Boundary Layer and Air Pollution Modelling. 16 May-3 June, 1994, no. SMR/760-2, pp. 712-743.

Prof. S. A. Pichugov, D. Sc. (Phys.-Math.); (Ukraine); Prof. V. D. Petrenko, D. Sc. (Tech.) (Ukraine)

recommended this article to be published

Accessed: Nov. 16, 2015

Received: March 20, 2016