УДК 628.511.133.622.063
С. А. Вамболь, Ю. А. Скоб, Н. В. Нечипорук, О. А. Трухмаев
МОДЕЛИРОВАНИЕ СИСТЕМЫ УПРАВЛЕНИЯ ЭКОЛОГИЧЕСКОЙ БЕЗОПАСНОСТЬЮ С ИСПОЛЬЗОВАНИЕМ МНОГОФАЗНЫХ ДИСПЕРСНЫХ СТРУКТУР ПРИ ВЗРЫВЕ МЕТАНОВОЗДУШНОЙ СМЕСИ И УГОЛЬНОЙ ПЫЛИ В ПОДЗЕМНЫХ ГОРНЫХ ВЫРАБОТКАХ
УГОЛЬНЫХ ШАХТ
Ключевые слова: экологическая безопасность, многофазные дисперсные структуры, взрыв, горня выработка.
Исследована эффективность водяных завес как средства влияния на процесс распространения ударных волн и высокотемпературной смеси продуктов сгорания при взрыве.
Keywords: ecological safety, dispersion multiphase structure, detonation, excavation.
The efficiency of water curtains as a means of influencing the propagation of shock waves and high-temperature mixture of combustion products in the explosion.
Введение
Взрывы газа и угольной пыли относятся к авариям с наиболее тяжкими последствиями в социальном и экономическом плане. Взрыв метано-воздушной смеси и угольной пыли в шахте, сопровождающийся формированием и распространением по сети горных выработок ударных волн - очень опасный вид подземной аварии. Угрожая жизни и здоровью горнорабочих, взрыв всегда приводит к разрушению выработок и повреждению расположенного в них горного оборудования. При этом социальные, материальные и финансовые потери могут достигать катастрофических размеров.
Для обеспечения безопасности условий труда рабочих на угольных шахтах необходимо эффективно производить гашение угольной пыли, так как пыль, находясь во взвешенном состоянии, в смеси с метаном, образует взрывчатую систему. Гашение угольной пыли в воздухе забоя может осуществляться: предварительным нагнетанием воды в пласт для его увлажнения, орошением забоя водой с подачей их, по возможности, непосредственно в места отбойки угля, а также дополнительно, путем создания заградительных завес. Таким образом, актуальной является техническая проблема повышение эффективности способов гашения угольной пыли путем снижения энергетических затрат на создание заградительных завес в сети горных выработок.
В настоящей работе рассмотрена математическая модель процесса взрыва облака метано-воздушной смеси и угольной пыли в подземных горных выработках угольных шахт. Проведен анализ способа защиты горных выработок от разрушительного воздействия взрывов путем создания водяных завес для снижения концентрации угольной пыли в воздухе.
Анализ последних достижений и публикаций
Взрыв взвешенной в воздухе угольной пыли может произойти от электрической искры, от раскаленных частиц взрывчатых веществ, вылетающих из шпура, и от других источников воспламенения [1].
Существенное влияние на взрыв угольной пыли оказывает присутствие метана СН4. Так, при отсутствии метана угольная пыль взрывается при содержании ее в воздухе не менее 30 - 40 г/м3, а при наличии СН4 2 % - опасная концентрация пыли понижается до 10 г/м3 и при наличии метана 3 % - эта концентрация еще понижается до 5 г/м3. Для предупреждения взрыва достаточно снизить концентрацию пыли до 5 г/м3, а с учетом запаса - до 1 г/м3 [2].
Одним из наиболее эффективных методов предотвращение чрезвычайной ситуации является создание водяных завес для снижения концентрации угольной пыли. При расходе воды до 35 л на каждую тонну добываемого угля, давлении воды у оросителей до 10 - 15 атмосфер и рациональном расположении оросителей эффективность способов гашения пыли предварительным нагнетанием и орошением достигает 95 %.
В качестве пылеподавляющих устройств используются обычные распылители. При этом необходимо учитывать, что в процессе создания водяной завесы размер капли предназначенной для захвата угольно-пылевой частицы варьируется от 50 -250 мкм, что существенно влияет на захват и осаждение угольной пыли.
В отличие от обычных распылителей, ато-майзеры способны создавать мелкодисперсную во-довоздушную смесь с витающими каплями, характерный размер которых составляет порядка 20 - 40 мкм [3 - 4]. Натурные испытания с отдельными экспериментальными образцами атомайзеров показали не только способность создавать частицы заданного размера, но и существенно снижать расход жидкости при создании водовоздушной завесы. Создание и удержание водовоздушного облака гораздо эффективней, чем при обычных распылителях.
Адекватное описание вышеназванных процессов - горение, выброс, смешение и распространение газовой примеси в атмосфере; истечение газовой примеси в атмосферу; взрыв газообразной примеси в атмосфере с учетом химического взаимодействия (в том числе в многосвязных областях со сложным рельефом местности или в закрытых объ-
емах), использующих для снижения экологического и пожарного рисков мелкокапельные технологии возможно только с использованием системы нестационарных уравнений Навье-Стокса для сжимаемого газа. В настоящее время численное моделирование турбулентных течений осуществляют путем решения осредненных по Рейнольдсу уравнений Навье-Стокса, дополненных моделью турбулентности [5, 6]. Однако большинство моделей турбулентности не описывают с одинаковой степенью адекватности различные типы течений, особенно это касается течений с интенсивными отрывами потока и/или большими градиентами давления температуры. Поэтому существует необходимость построения новых моделей и расчетных схем для численного моделирования таких течений.
Постановка задачи
Необходимо определить эффективность водяных завес как средства влияния на процесс распространения ударных волн и высокотемпературной смеси продуктов сгорания при взрыве облака мета-новоздушной смеси и угольной пыли в подземных горных выработках угольных шахт на основе математической модели рассматриваемого процесса.
Рудничная пыль представляет собой твердое вещество полезного ископаемого или породы в состоянии тонкого раздробления (измельчения). Находясь во взвешенном состоянии в воздухе, пыль образует дисперсную систему (аэрозоль). Частички пыли составляют дисперсную фазу, воздух является дисперсионной средой.
Математическая модель экологических процессов в атмосфере
В результате структурного анализа перечисленных экологических процессов в атмосфере (выброса, смешения и распространения газовой примеси в атмосфере; истечения газовой примеси в атмосферу; взрыва газообразной примеси в атмосфере с учетом химического взаимодействия (в том числе в многосвязных областях со сложным рельефом местности), с учетом термодинамических свойств газов, а также алгоритма численного моделирования этих процессов), декомпозиции полной математической модели процесса выявлено, что основное влияние на процесс оказывает интенсивность конвективного и турбулентного переноса компонент примеси. Таким образом, для описания процессов смешения многокомпонентного газа в рассматриваемом случае достаточно использовать усеченные уравнения, полученные из уравнений Навье-Стокса путем отбрасывания вязких и диффузных членов - приближение Эйлера с источнико-выми членами.
Расчетной областью ^ является параллелепипед с прямолинейными образующими (рис. 1), расположенный в правой декартовой системе координат (X, У, 2) с основанием в плоскости Х02 (ось У ориентирована в направлении, противоположном действию сил тяжести Земли).
Рис. 1 - Расчетная схема смешения реагирующих газов с воздухом, их горения (с учетом химического взаимодействия компонент смесей) и дальнейшего рассеяния смеси
Расчетная область разбивается на пространственные ячейки, причем размеры граней подбираются в соответствии с характерным размером особенностей расчетной области (шероховатости обтекаемой поверхности, размерностью обтекаемых объектов).
Полная система уравнений, описывающая нестационарное трехмерное течение трехкомпо-нентной смеси газов в данной постановке имеет вид [7 - 9]: ^ ^
За ЗЬ Зс ^ а — + — + — + — = pf, (1) З ^ ^ Зг
где а, Ь, с, с1, ^ - вектор-столбцы вида:
а = ^,ри, рр, рр, Е]т,
Ь =
с =
ри, Р + ри2, ррт/, ри, (Е + Р)]7 ру, рру, Р + ру 2,р\лл/ , (Е + Р)у^, 8 = pwv, Р + pw2, (Е + ,
? = [0, 0, -д, 0,+ е8/рГ,
где 1 - время; и,у,ш - составляющие вектора скорости ^ ; Р, р - давление и плотность; Е - полная энергия единицы объема смеси газов
(2)
(3)
(4)
(5)
(6)
Е = р£+ 2 р2 +У2 + W2) ;
(7)
е - внутренняя энергия единицы массы газа; компоненты вектора ^ - суть проекции распределенных объемных источников, д - ускорение свободного
падения, вд - интенсивность тепловыделения в
единице объема газа вследствие химической реакции.
Закон переноса компоненты смеси (горючий газ, воздух, продукты сгорания) с учетом скорости турбулентной диффузии, имеет вид [7 - 9]:
З(рО|) , З(риО|) , З(рУР|) ,
З
Зх
Зу
+ З(рл/^) = Зг
Он
+ Рс
где Qi - относительная массовая плотность 1 той компоненты смеси, i = 1...3, (отношение плотности газообразного вещества компонента к плотности смеси): 1 - горючее, 2 - воздух, 3 - продукты сгорания; pQit - интенсивность изменения плотности
компоненты смеси вследствие турбулентной диффузии в соответствии с законом Фика PQit = div(pSDgradQ¡) (коэффициент турбулентной диффузии ^ определялся по методике, предложенной М.Е.Берляндом [10]); pQ¡s - интенсивность изменения плотности компоненты смеси, вследствие химической реакции.
Система уравнений (1 - 8) является незамкнутой. Дополним ее уравнениями, определяющими теплофизические свойства смеси:
м = сР=£сфр); V—
¡=1
Ы Mi
С. = 10;(сЛ; 10; = 1; к = с ¡=1 ¡=1 ^ Для идеального политропного газа величина е связана с Р и р смеси зависимостью:
Р
е =
(- 1)р'
Интенсивность изменения плотности горючего pQ1s, вследствие химической реакции, определялась как произведение молекулярной массы горючего щ на скорость его мольного изменения . Рассматривалась химическая реакция-«брутто»: 2
^ У3Х3 ,
¡=1
где v¡ - стехиометрические коэффициенты, ш -химические вещества.
Скорость мольного изменения определялась по формуле:
¡=1
где А^ Р^ Е^ V!, V'-, - параметры, полученные на основе обобщения экспериментальных данных [11], - мольная концентрация ьтой компоненты
смеси: [^¡] = —^^. Интенсивность изменения
плотности продуктов сгорания pQ3s определялась на основе закона движущихся масс:
РС>1з = —023 = РС>3з ^ 2 ^3
Интенсивность тепловыделения в единице объема газа вследствие химической реакции es определялась по формуле:
ез = ЧН ^Рс ^ ,
где £ - коэффициент полноты сгорания, Hu-| -
низшая теплота сгорания горючего.
При постановке граничных условий полагается, что расходная составляющая скорости не превосходит скорость звука. Граничные условия на входе будем задавать на поверхностях тех граней, примыкающих к границам расчетной области, через которые в расчетную область поступает атмосферный воздух. Набегающий поток на входе определяется величинами:
- полной энтальпии
k Р u2 + v2 + w2
—; (9)
00
к-1 р
- функции энтропии
2
S = Р..
=рк.
(10)
- направлением вектора скорости потока (углами a x, а у, а z).
- относительной массовой плотностью горючего Qi (Qi < 1, если поступает газообразное вещество компоненты смеси, например, метан).
Параметры потока на входе определяются из равенств (9, 10) с учетом заданных a x, а у, а z с
привлечением соотношения для "левого" инварианта Римана [12]. На непроницаемых участках, ограничивающих расчетную область поверхностей, выполняются условия «непротекания»: qn = 0, где n - вектор нормали к рассматриваемой границе. Граничные условия на выходе будем задавать на поверхностях тех граней конечно-разностных ячеек, которые примыкают к границам расчетной области и через которые предполагается вытекание или втекание смеси. В выходных областях, кроме атмосферного давления Ра|р, задаваемого либо взятого из эксперимента, использовались соотношения для "правого" инварианта Римана [12]. Задавались «мягкие» граничные условия для искомых относительных массовых плотностей компонент на всех ограничивающих расчетную область поверхностях.
В начальный момент времени во всех «газообразных» ячейках расчетной области принимаются параметры окружающей среды. В ячейках, занимаемых облаком газообразной компоненты смеси, которое образовалось в результате мгновенного выброса, относительная массовая концентрация компоненты смеси принимается равной Qi < 1. В ячейках с испарением или истечением газа задается закон изменения расхода компонента смеси.
Предполагалось, что горение происходит в объеме расчетной области, занимаемом горючей смесью с концентрацией горючего в диапазоне между минимальным и максимальным концентрационными пределами воспламеняемости
Q1 min < Q1 < Q1 max. Величины Qimin , Q1 max
задавались на основе обобщения экспериментальных данных [13].
Учет распыления капель (частиц пыли твердых фракций)
Качество распыла, как известно, характеризуется интегральными и дифференциальными кривыми распределения объемов (числа, поверхности) капель по их диаметрам и различными понятиями среднего диаметра капель. В большинстве случаев используемая в технике распыленная жидкость состоит из капель различных размеров, т.е. имеет полидисперсный характер. Для описания кривых распределения капель по размерам предложены различные зависимости. Наибольшее распространение нашло уравнение Розина-Раммлера [14]:
Р = 1 - ехр (- 0/С1+ 7 ;.
Здесь Р - объемная доля капель, диаметр которых меньше С; С+ - характерный размер или средний
диаметр, соответствующий определенному значению Р = 0,3679; п - константа распределения, характеризующая степень неоднородности распыла (обычно 2< п <4). Тогда плотность распределения капель по диаметрам примет вид:
Р/+>=П ¿Гехр (- ¿3;. с!+ с1+
Учет влияния аэродинамического
сопротивления капель (частиц пыли твердых фракций) на параметры движения газовой смеси
Определим силу сопротивления, действующую на каплю (частицу пыли):
э г рЧ 2 ст а
ГТр = "^^с!— СТ-
N1
тсС
2
где ст| = —4--площадь поверхности, Ыс - коэффициент аэродинамического сопротивления, д -скорость потока. Если капля сохраняет сферическую форму, то для определения Ыс можно рекомендовать соотношение [15]:
Сй = 24/Ре+ 4,4/7^6 + 0,35,
где Ре = ^^ - число Рейнольдса. Тогда осреднен-
ная по диаметрам капель сила аэродинамического сопротивления для выбранного контрольного объема определяется по формуле:
а I а
Роб = ^ Ер(С|)Роб(С|) 1=1 ,
где
^О т
- количество капель в
™рнр ¿р Й У*
контрольном объеме, б^О - общий расход воды
(предполагается, что вода непрерывно подается в расчетную область), т - шаг по времени, т - количество контрольных объемов в области с источни-
ками (капели, частиц пыли твердых фракций),
3
У|| = (тс / 6)С| - объем капли.
Количество частиц твердых фракций в контрольном объеме определяется по формуле:
р +АУ N. =—-
к I '
рс^р ^ У*
¡=1
где р с + - концентрация пыли во взвешенном состоянии, ДУ - контрольный объем, р с - плотность вещества твердой фракции.
Учет влияния аэродинамического сопротивления капель (частиц пыли твердых фракций) на параметры движения газовой смеси осуществлялся посредством введения в уравнения движения осред-ненной объемной силы сопротивления: а 1 -
т = — Д\/. При этом предполагалось, что р
сумма удельной мощности сил сопротивления и удельной диссипируемой мощности равна нулю:
р?тр)+М, = 0.
Учет влияния тепловыделения вследствие
фазового перехода при кипении капель (химической реакции горения частиц пыли твердых фракций) на параметры движения газовой смеси
Определим интенсивность изменения плотности примеси вследствие фазового перехода при кипении капель воды:
= 6н2о
= тДУ,
а интенсивность изменения плотности примеси вследствие химической реакции горения:
Рс8 = ——, С8 тДУ
где бс - общее изменение в единицу времени массы частиц пыли твердых фракций вследствие химической реакции горения.
Будем исходить из того, что известен закон, согласно которому время горения одной частицы зависит от начальной массы этой частицы:
(
ь = кс2 = к
6Мо|
2 3
где к = 6,14 • 106
тср N
М01 =р с^к1. Тогда согласно
известным из математического анализа и теории обыкновенных дифференциальных уравнений преобразованиям может быть получен закон определения массы частицы в произвольный момент времени химической реакции горения:
2
2 3
М| = [М3| -1,06 • 10-7рЫ3 • ь]2
1=1
а также зависимость изменение в единицу времени массы частиц пыли твердых фракций вследствие химической реакции горения:
7 2 2 7 2 1 Gci = 1.6 • 10 PC3(M3. -1.06.10"' PC3 • t)2
где t - текущее время.
Общее изменение в единицу времени массы частиц пыли твердых фракций вследствие химической реакции горения может быть определено по формуле:
I
Gc = mNk Zp(dj)GCi i=1 .
Учет влияния фазового перехода при кипении капель воды (химической реакции горения частиц пыли твердых фракций) на параметры движения газовой смеси осуществлялся посредством введения в уравнение энергии интенсивностей тепловыделения в контрольном объеме:
eH2Os = -r(P)PH2Os , eCs = £HuCPCs ,
где r(P) - удельная теплота парообразования, £ -
коэффициент полноты сгорания, Huc - низшая
теплота сгорания.
Таким образом, предложена математическая модель смешения газов, которая описывает экологические процессы в атмосфере, что позволит прогнозировать поля температуры, давления и концентрации вещества примеси в рассматриваемом объеме.
Метод представления и численное решение задач смешения газов в атмосфере
Векторное уравнение (1) является следствием законов сохранения массы, импульса и энергии, которые могут быть представлены в интегральной форме для каждой расчетной ячейки:
д
-ш adV + # Ada = JJjpfdV, (11) dtV a V
где V - объем элементарной расчетной ячейки; с -ограничивающая поверхность данной ячейки, которая имеет внешнюю нормаль n (a = an); А - тензор
плотности потока консервативных переменных a , столбцами которого являются векторы b,C e d , соответственно.
Закон переноса компоненты смеси (8) может быть также представлен в интегральной форме для каждой расчетной ячейки:
д_ dt
ШрОСУ + #рОдСст = Ш (роь +ро3)СУ (12) V ст V
Совокупность газодинамических параметров во всех ячейках в момент времени Ьп представляет собой известное решение на временном слое с индексом п. Параметры в момент времени
Ьп+1 = Ьп +т (на слое п +1) рассчитывались посредством применения явных разностных аппроксимаций для соотношений (11, 12) в рамках интег-ро-интерполяционного метода С.К. Годунова [12]
первого порядка точности. Устойчивость конечно-разностной схемы обеспечивается за счет выбора величины шага по времени т в соответствии с условием Куранта.
Компьютерное моделирование рассматриваемых процессов осуществлялось с использованием программного комплекса FIRE® [16], реализующего предложенную математическую модель. Программный комплекс FIRE® позволяет вычислять плотность, скорость, давление, температуру смеси, концентрации компонентов смеси (топлива, воздуха, продуктов сгорания), скорость тепловыделения в пределах каждого контрольного объема смеси на каждом дискретном временном шаге. Приложение написано в интегрированной среде разработки Visual C++ 4.0 Enterprise Edition с применением библиотеки Microsoft Foundation Classes (MFC).
Результаты исследования водяных завес как средства снижения последствий взрывов метано-воздушной смеси и угольной пыли в подземных горных выработках угольных шахт
Проведено компьютерное моделирование процесса взрыва облака метано-воздушной смеси и угольной пыли в штреке угольной шахты. Расчетная область представлена на рис. 2. Вычислительный эксперимент проводился при следующих условиях: длина штрека Lz = 31,2 м, высота штрека была равна Ly = 2,2 м, скорость ветра на входе составляла q = 6 м/с, облако метано-воздушной смеси располагалось на расстоянии Z1 = 10,1 м от входа, радиус облака составлял R1 = 1,6 м (см. рис. 2а). На расстоянии Z2 = 13 м располагалась зона, заполненная мелкодисперсной фазой - частицами угольной пыли либо каплями воды (в зависимости от рассматриваемого сценария расчета). Предполагалось, что в случае наличия в этой зоне капель воды, обеспечивается полное осаждение угольной пыли. На расстоянии Z3 = 14,9 м располагалась контрольная точка P, в которой осуществлялись изменения избыточного давления, температуры газовой смеси. Форма поперечного сечения штрека представлена на рис. 2б, здесь Lx = 3,2 м, Y1 = 0,6 м, R2 = 1,6 м.
Рассматривались три возможных сценария развития событий, связанных с процессом взрыва облака метано-воздушной смеси в штреке. В начальный момент времени в результате взрыва мета-но-воздушной смеси формировалось облако продуктов сгорания с высоким давлением и температурой. Далее реализовался процесс рассеяния продуктов сгорания, сопровождающийся конвективным переносом и турбулентным рассеянием продуктов сгорания вдоль штрека. Согласно принятым расчетным схемам в первом случае (сценарий 1) рассматривалось рассеяние продуктов сгорания вдоль штрека при наличии угольной пыли в воздухе, во втором случае (сценарий 2) - без наличия дисперсных фаз, в третьем случае (сценарий 3) - при наличии капель воды (водяной завесы) в воздухе.
о
ш
N-^
22
гз
Рис. 2 - Схема расчетной области
Как следует из анализа результатов расчетов, присутствие в воздухе угольной пыли приводит к росту избыточного давления и температуры в штреке (по сравнению со случаями, когда воздух не содержит дисперсных фаз - кривые 2, при наличии водяной завесы в воздухе - кривые 3 на рис. 3), вызванных возгоранием этой фазы (см. кривые 1на рис. 3 а и рис. 3 б).
Р, кПа'
150-
100-
т, к-
3000-
2500-
2000-
1500.
1000-
500
1
/ г / \ V 1
/ 2 у / / г * * * / V \ V
у __ I / 3 /
^ V Z__________________
0,05
0,1 б
0,15
I, с
Рис. 3 - Изменение избыточного давления (а) и температуры в контрольной точке: 1 - при наличии угольной пыли в воздухе; 2 - воздух без дисперсных фаз; 3 - при наличии водяной завесы в воздухе
Наличие капель воды обеспечивает не только полное осаждение угольной пыли, но и дополнительное снижение избыточного давления и температуры за счет перехода жидкой фазы в парообразное состояние при прохождении вдоль штрека высокотемпературного облака продуктов сгорания метано-воздушной смеси (см. кривые 3 на рис. 3а и рис. 3б). Для сравнения приводятся результаты расчетов изменения избыточного давления и температуры в контрольной точке для случая, когда воздух в штреке не содержал дисперсных фаз (см. кривые 2 а рис. 3 а и рис. 3б). Как видно из анализа данных эффективность водяных завес снижается при гашении взрывов только метано-воздушной смеси (т.е. без участия во взрыве угольной пыли).
Выводы
В данной работе рассмотрен способ защиты горных выработок угольных шахт от разрушительного воздействия взрывов путем создания водяных завес для снижения концентрации угольной пыли в воздухе.
На основе численного моделирования решена задача определения необходимого расхода жидкости в водяной завесе и ее местоположения с целью осаждения угольной пыли, но и дополнительного снижения избыточного давления и температуры за счет перехода жидкой фазы в парообразное состояние при прохождении вдоль штрека высокотемпературного облака продуктов сгорания метано-воздушной смеси в подземных горных выработках угольных шахт.
Взрыво-локализующее действие водяных завес состоит в создании гасящей среды на пути распространяющегося по горной выработке фронта пламени от взрыва метано-воздушной смеси. Эффективность водяных завес снижается при гашении взрывов только метано-воздушной смеси (т.е. без участия во взрыве угольной пыли).
а
б
а
Литература
1. Шохин В. Н. Гравитационные методы обогащения / В. Н. Шохин, А. Г. Лопатин. - М.: Недра, 1993. - 313 с.
2. Абрамов А.А. Переработка, обогащение и комплексное использование твердых полезных ископаемых. Т. I. Обогатительные процессы и аппараты: учебник для вузов. - М.: Изд-во Московск. гос. горн, ун-та, 2001.- 472 с.
3. Ляшенко А.М. Физические аспекты пылеподавления и распылительные устройства для их реализации / А.М. Ляшенко,. Н.В.Нечипорук, Н.В. Кобрина, С.А. Вамболь // Открытые информационные и компьютерные интегрированные технологии. - Х: НАКУ «ХАИ», 2010. -Вып. № 48. - С. 234-239.
4. Кобрина Н.В. Моделирование процесса пылеподавле-ния при погрузке, разгрузке и транспортировке сыпучих материалов / Н.В.Кобрина , В.Е. Костюк , В.Н. Кобрин, С. А. Вамболь / Открытые информационные и компьютерные интегрированные технологии. - Х: НАКУ «ХАИ», 2010. - Вып. № 48. - С. 248-252.
5. Шенг Дж. С. Обзор численных методов решения уравнений Навье-Стокса для течений сжимаемого газа // Аэрокосмическая техника.- 1986.- №2.- С. 65-92.
6. Вычислительная гидромеханика и теплообмен: В 2-х т. / Андерсон Д., Таннехилл Дж., Плетчер Р. - М.: Мир, 1990. - 726 с.
7. Numerical Modeling of Hydrogen Release, Mixture и Dispersion in Atmosphere / E.A. Granovskiy, V.A. Lyfar, Yu.A. Skob, M.L. Ugryumov // Abstracts Book и CD-ROM Proceedings the International Conference on Hydrogen Safety. - Pisa (Italy). - 2005. - 10 p. (ICHS Paper No. 110021)
8. Computational Modeling of Pressure Effects from Hydrogen Explosions / E.A. Granovskiy, V.A. Lyfar, Yu.A. Skob, M.L. Ugryumov // Abstracts Book и CD-ROM Proceed-
ings of the 2-nd International Conference on Hydrogen Safety. - San Sebastian (Spain). - 2007. - 15 p. (ICHS Paper No. 1.3.52)
9. Numerical Modeling of Hydrogen Deflagration Dynamics in Enclosed Space / Yu.A. Skob, M.L. Ugryumov, K.P. Ko-robchynskiy, V.V. Shentsov, E.A. Granovskiy, V.A. Lyfar // Abstracts Book and CD-ROM Proceedings of the 3-nd International Conference on Hydrogen Safety. - Ajaccio-Corsica (France). - 2009. - 12 p. (ICHS Paper No. ID 182)
10. Современные проблемы атмосферной диффузии и загрязнения атмосферы / Берлянд M.E. - Ленинград: Гидрометеоиздат, 1975. - 448 с.
11. Физико-химические процессы в газовой динамике. Комьютезированный справочник в 3-х томах / Под редакцией Г.Г. Черного и С.А. Лосева. - Т. 2. - Физико-химическая кинетика и термодинамика. - М.: Научно-издательский центр механики, 2002. - 368 с.
12. Численное решение многомерных задач газовой динамики / С.К. Годунов, А.В. Забродин, М.Я. Иванов, А.Н. Крайко, Г.П. Прокопов. - М.: Наука, 1976. - 400 с.
13. Физика взрыва / Под ред. К.П.Станюковича. - М.: Наука, Главная редакция физико-математичекой литературы, 1975. - 704 с.
14. Пажи Д.Г., Галустов В.С. Распылители жидкостей. -М.: Химия, 1979. - 216 с.
15. Нигматулин Р. И. Динамика многофазных сред. Т.1. -М.: Наука, 1987. - 384 с.
16. Комп'ютерна програма „Комп'ютерна штерактивна система шженерного аналiзу та прогнозу руху хiмiчно реагуючих газоповггряних сумшей в задачах промислово! аеродинамши та екологп атмосфери «FIRE»": Свщоцтво про реестращю авторського права на твiр № 30079 / Скоб Ю.А., Угрюмов М.Л., Короб-чинський К.П. (Укра!на). - Дата реестрацп 28.08.2009.
© С. А. Вамболь - к.т.н., доц., зав. каф. прикладной механики Национального университета гражданской защиты Украины, sergvambol@gmail.com; Ю. А. Скоб - к.т.н., доц. каф. информатики Национального аэрокосмического университета им. Н.Е. Жуковского «ХАИ», k106@mail.ru; Н. В. Нечипорук - д.т.н., проф., проректор по научно-педагогической работе ХАИ; О. А. Трухмаев - директор ООО «Универсалтехсервис», г. Мариуполь, violavambol@gmail.com.