Вычислительные технологии
Том 13, Специальный выпуск 2, 2008
Импульсная система пожаротушения
*
на основе твердотопливного газогенератора*
А. Д. Рычков
Институт вычислительных технологий СО РАН, Новосибирск, Россия
e-mail: [email protected]
A scheme of pulse system of fire fighting (PSFF) for intensive heat sources extinguishing with the help of a solid-propellant gas generator is suggested. Dispersion of the finely dispersed powder of the flame retardant substance and its transport into the combustion zone in the form of an aerosol cloud is made possible due to the high pressure created by the burning products of the gas generator. Mathematical model of processes in PSFF is proposed. Results of numerical modelling suggest that this devise provides a rather stable structure of the particle cloud and which is able to support an effective transport of aerosol particles into the combustion zone.
Введение
Задача эффективного гашения интенсивных очагов пламени, внезапно возникающих при аварийных ситуациях и техногенных катастрофах, является чрезвычайно актуальной. Решение этой задачи обусловливает выбор экологически безопасного, нетоксичного и приемлемого с экономической точки зрения пламегасящего вещества (ПГВ), а также эффективных средств его доставки к очагу пламени. Наиболее перспективны для этих целей фосфорорганические соединения (ФОС), металлоорганические соединения или МОС (органические соли и комплексные соединения К, Ха. Са и др.), а также некоторые галогенсодержащие соединения (СН31, СР31). Суть механизма гашения состоит в том, что в пламени эти вещества разлагаются с образованием активных соединений, которые вступают в химические реакции с радикалами (частицами, ответственными за горение и распространение пламени), переводя их в неактивные частицы, тем самым обрывая цепные разветвленные реакции горения, в результате чего и происходит гашение пламени. Для эффективного гашения пламени достаточно объемной концентрации этих продуктов разложения порядка 1-2 %. Таким образом, эти классы соединений являются перспективными кандидатами для замены хладонов, широко используемых в настоящее время в средствах пожаротушения, но запрещенных к производству Монреальской конвенцией 1987 г.
* Работа выполнена при финансовой поддержке Программы междисциплинарных интеграционных исследований СО РАН (проекты № 28 и № 113), Программы фундаментальных исследований ОЭММ-ПУ РАН (грант № 4.1.1), Программы Государственной поддержки научных исследований, проводимых ведущими научными школами Российской Федерации (грант НШ-931.2008.9).
© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2008.
Главное препятствие для применения большинства как уже известных, так и новых эффективных ПГВ это потери при их доставке к очагу пламени. Применение мелкодисперсных аэрозолей пламегасителей с размером частиц порядка 20...40 мкм и более мелких уменьшает эти потери. Однако получение аэрозолей с такими характеристиками в средствах пожаротушения на данный момент является технически трудной задачей. Как правило, для этого используются унитарные аэрозольгенерирующие составы. Такие составы используются в системах пожаротушения на основе твердотопливных газогенераторов, в продуктах сгорания которых содержатся мелкодисперсные частицы пламегасителя в достаточно высокой концентрации, способные эффективно подавлять очаги возгорания за относительно короткое время. Однако их основной недостаток заключается в серьезных технологических трудностях, связанных с приготовлением твердотопливных составов, обладающих необходимыми свойствами, а также в относительно малой газопроизводительности, что затрудняет их использование при тушении очагов возгорания на больших площадях (например, при лесных пожарах).
Как альтернативный и более эффективный метод пожаротушения в данной работе предлагается использование устройства раздельного снаряжения, в котором распыление пламегасящего вещества и транспортировка его в зону горения осуществляются продуктами сгорания обычного твердотопливного газогенератора кратковременного (импульсного) действия, продукты сгорания которого имеют относительно низкую температуру и не содержат окислитель. Важной проблемой здесь является разработка конструкций таких устройств и составов твердотопливных зарядов для газогенераторов, которые, с одной стороны, были бы сравнительно безопасны, а с другой обеспечивали бы нужные параметры получаемого с их помощью облака аэрозолей пламегасителей, необходимые для эффективного тушения пламени.
Одна из возможных схем такой импульсной аэрозольной системы пожаротушения (ИАСП) приведена на рис. 1. Устройство представляет собой цилиндр, в центральной части которого расположен заряд твердого топлива в виде гранул. Периферийная часть цилиндра содержит распыляемый мелкодисперсный порошок пламегасителя. Между центральной и периферийной частями имеется перегородка, не оказывающая заметного гидродинамического сопротивления движению газа. Верхняя и нижняя крышки цилиндра непроницаемы и неподвижны, так что движение частиц в нем под воздействием газообразных высокотемпературных продуктов сгорания высокого давления происходит только в радиальном направлении. Тем самым в течение короткого времени (50... 100 мс) осуществляются генерация и транспорт облака ПГВ аэрозоля высокой плотности в зону горения, что не удается сделать с помощью традиционных средств пожаротушения.
При работе газогенератора происходит компактирование частиц, а после их выхода
Газогенератор
Мелкодисперсный \пламегаситель
\ Крышка
Рис. 1. Схема устройства
из цилиндра под воздействием истекающего газа начинается процесс формирования облака частиц. Возникающая при этом баллистическая волна низкой интенсивности сбивает пламя, а также распыляет и диспергирует пламегаситель, который в виде аэрозольного облака доставляется в зону горения. При испарении частиц пламегасителя из-за поглощения ими тепла температура газовой фазы в зоне облака частиц понижается и, кроме того, происходит вытеснение кислорода воздуха от очага горения, так как продукты сгорания твердотопливного газогенератора его не содержат.
Сочетание указанных выше факторов (действие баллистической волны, охлаждение газа в зоне облаке частиц, вытеснение из этой зоны кислорода воздуха, наличие паров ПГВ) позволяет тушить пожары различных типов, включая лесные. Предлагаемый метод тушения пожаров отличается взрывом, когда возникающая ударная волна вызывает довольно сильные разрушения, но эти методы не применимы для тушения пожаров в производственных помещениях. Кроме того, он имеет еще одно важное преимущество - для тушения лесных пожаров устройства подобного типа можно сбрасывать с вертолетов или метать любыми способами как на кромку пожара, так и внутрь очага, где происходит их срабатывание. При этом тушение пламени происходит не только на передней кромке пожара (что характерно для традиционных способов тушения с помощью огнетушителей или брандспойтов), но также и внутри очага пламени.
По-видимому, наиболее эффективным будет кассетное применение ИАСП, поскольку при этом охватывается большая площадь горящих материалов, в результате чего пожар можно быстрее локализовать и потушить. Также следует отметить постоянную готовность такого устройства к действию, малое время срабатывания, простоту использования и доставки. Средства пожаротушения на основе ИАСП могут найти применение не только при тушении лесных пожаров, но и на транспорте, для зашиты производственных и нежилых помещений, а также для защиты от пожаров различных важных объектов (включая стратегические).
Заметим, что такого рода устройства могут применяться кроме ИАСП также для быстрого распыления различных сорбентов и порошков при химической обработке больших площадей и во многих других практических приложениях.
1. Описание модели и основные уравнения
Продукты сгорания топливных гранул в момент выхода "пробки" насыпного материала из цилиндра имеют давление 10...15 МПа, их истечение в окружающую среду происходит в сверхзвуковом режиме, и при описании их движения необходимо учитывать турбулентный характер течения. Поэтому в данной работе для моделирования процессов формирования и движения облака мелкодисперсных частиц ПГВ использовалась осредненная по Фавру система уравнений Навье-Стокса, замыкаемая д - ^-моделью турбулентности [1]. Описание движения полидисперсной второй фазы проводилось в рамках стохастического подхода [2], учитывающего влияние турбулентной природы поля течения несущего газа на движение частиц.
При описании математической модели рассматриваемых процессов были приняты следующие допущения:
— течение является турбулентным, двухфазным, осесимметричным, нестационарным и многокомпонентным. Истечение происходит в неподвижную атмосферу стандартного состава с высокой температурой;
— продукты сгорания газогенератора представляют собой инертный газ;
— частицы ПГВ являются полидисперсными, имеют сферическую форму и состоят из N фракций, внутри каждой из которых частицы одинакового размера. Столкновениями между частицами на данном этапе моделирования пренебрегается;
— распределение температуры по объему частицы полагается однородным;
— при нагреве частицы до определенной температуры начинается ее испарение, линейная скорость которого гр предполагается заданной;
— газоприход от гранул твердого топлива и от частиц при их испарении моделируется источник!)вымп членами в уравнениях сохранения массы и энергии;
— химические реакции в среде не рассматриваются, учитывается только смешение компонентов.
Система уравнений, описывающая такое течение в цилиндрических координатах, имеет следующий вид:
др д 1 д ^^ / ¿шг \ , ,
т + Тх{ри) + гЪ~г{гру) = м^
г=1 х '
дрУо2 д 1 д . , .
+ ТГУРиУо2 + Ях,о2) + -т-г(ру¥о2 + дг,о2) = 0; (2)
дЬ дх 2 ' 2 г дг
дрУр д 1 д ^^ / ¿шг\ , ,
~дГ + Тх{ри¥р + *°>м) + гд-гГ{ру¥р + = ®
г=1 х '
дрУм д 1 д
~~д1--^ дх^ри¥м + + гдгГ^р1>¥м + = М%еп] ^
дри д , 2 ,1 д , , ^ / ¿иЛ ,,
г=1 ^ '
дру д 1 д 2 ^^ / ¿уг \ , ,
г=1 ^ '
дрЕ д..^. .1 д./гп. .
~дГ + дх + Р) + иТхх + УТхг + V*) + ~ + Р) + + + Ы =
N
=1
и { + и \ ) + (^ДМи ~
+ ^^п^^^е щ (7)
дрд д 1 д
+ + + г^РЩ + Тга) = ^ ^
дри д . .1 д . ч о
~дГ + ¿^1+ + ~гЪ~гг^риш + Тгш) = ^
+у,=1. <«»
Уравнения движения г-й частицы ПГВ вдоль ее траектории записывались виде
^ = -Рр^шу, (11)
'„.Ко .
(и + и' - щ); (12)
dui 3 ßCDiRep dt 4 ppd2
о
(v + v'~ Vi); (13)
dvi 3 ßCDiKep
dt 4 ppd2
dTj _ 6 ANu
dt Ci Pp d2
(T - Ti); (14)
dxi dyi , ,
-Ж = "" ¿r =(15)
Здесь q = v'7- и со = e/k — "турбулентные" переменные, связанные с кинетической энергией турбулентности & и со скоростью ее диссипации е; E = CvT + (u2 + v2)/2 — полная удельная энергия; ß и ßt — молекулярная и турбулентная вязкости соответственно; p и T — статическое давление и температура; R0 — газовая постоянная; величины Mgen и Qgen есть источники массы и энергии, моделирующие приход высокотемпературных продуктов сгорания от газогенератора; Nu и Ее — числа Нуссельта и Рейнольдса соответственно.
Sq = Cqi(CßS/u2 - 2D/(3u) - 1) puq, = [Cwl(CßS/u2 - ОзD/u) - O2W-Константы, входящие в описание модели турбулентности, имели такие значения: Oi = 0.55, C.2 = 0.833, Оз = 0.666, Cß = 0.09, Cqi = 0.5, ßt = Cß pq2/u,
fldru dv\ dT dT
Trx TXr ße I ö т; I , ^Px Ae — j (fir Ae — ,
\r dr dx ) dx dr
ße = ß + ßb Ae = A + A t, A t =
Pit
du 1 drv , , 2
k,k
. ßt\dYk ( ßt\dYk
Qx,k = ~ pDk + —— -7—, qr>k = - pDk + —
SctJ dx ' \ Set) dr '
k — O2, N2, M; P — пламягасящее вещество; члены в угловых скобках означают осреднение по объему ячейки разностной сетки; n — концентрация частиц г-й фракции в ячейке в данный момент времени;
24 / Re2/3\ 3 , ,
С,., < —J' 6СЛИ ** " 10 ' Re* = (16)
0.44, если ReP;i > 103,
Здесь и' я у' — компоненты случайного вектора возмущения скорости газа V =
р2к V
(V, г/); V = у—ег/с_1(А^)-—-; V — вектор осредненной скорости; Ыг — случайное
V 3 | V |
число из диапазона [-1,+1]; ет/е—1( ) — обратная функция ошибок; С — удельная теплоемкость г-й частицы, xi и у — ее координаты.
Для расчета течения несущего газа применялась противопотоковая 1.1' разностная схема второго порядка точности, обладающая ТУБ-свойствами, близкая к схеме из работы [3]. Расчет движения частиц проводился с помощью А-устойчивой разностной схемы также второго порядка точности [4]. Область решения системы уравнений (1)-(8) снизу ограничивалась непроницаемой поверхностью (подстилающая поверхность, ось ОУ направлена горизонтально, ось ОХ — вертикально), слева — осью симметрии течения, справа и сверху не ограничивалась ничем (безграничное пространство). Граничные
ОХ
условия симметрии течения, на поверхностях крышек цилиндра и на подстилающей поверхности — условия прилипания, на остальных границах — неотражающие граничные условия. Для системы (9)—(14) на подстилающей поверхности задавались условия неупругого отражения
5
= -(0.993 - 1.760 + 1.5602 - 0.49е>+„, лг~т = -лг+
7
Т 1
— + — +
где vi п, vJn, vi Т, vJТ — нормальная и тангенциальная компоненты вектора скорости г-й частицы; сим вол "+" относится к падающей, а символ "—" — к отраженной частицам; © — угол падения.
Величина массоприхода Мёеп продуктов сгорания газогенератора в (1) определялась из решения уравнения
^ — —РдгйГьд1Удд,
где рдг, Б и тьд — плотность материала топливной гранулы, ее текущая поверхность и скорость горения; М0д — число гранул в единице объема заряда твердого топлива, которое остается постоянным в процессе их горения и определяется из начальных условий.
2. Некоторые результаты расчетов
Расчеты проводились при следующих исходных условиях: диаметр цилиндра ПЛИС равен 0.33 м; высота цилиндра 0.2 м; масса частиц ПГВ составляла 10 кг; скорость испарения частиц гр = 10 мм/с; число фракций равнялось десяти; их среднемассовый размер ^43 = 30 мкм; температура испарения частиц ПГВ Тисп = 400 К. Заряд гранулированного твердого топлива газогенератора имел массу 0.1 кг; теплотворность топлива Я^п = 1780 кДж/кг. Температура окружаю щей среды Т0 = 800 К. Расстояние от нижней крышки цилиндра ИАСП до подстилающей поверхности в первом варианте расчета составляло К = 1 м и К = 0.1 м — во втором.
На рис. 2 показаны структура поля течения и положение частиц ПГВ в момент времени Ь = 10 мс. Распределение температуры свидетельствует о ее понижении в области, занятой частицами ПГВ, из-за поглощения тепла частицами, начальная температура которых 300 К, а температура испарения ниже температуры окружающей среды.
Вихревая структура за облаком частиц обеспечивает устойчивую компактную форму облака аэрозоля в форме тора, генерирующую достаточно высокую концентрацию паров ПГВ (рис. 3).
Рис. 2. Поле течения при Н = 1 м (Ъ = 10 мк)
ЯШрйр - 10.0 Т таи - 2416. Рпвк й 0.1? 2.00 1.83 1.69 Т.55: 1.42 1.29 1.17 1.04 0.92 0.80 0.68 0.56 0.45 0.33 0.21 0.01
Изолинии концентрации паров ПГВ
Л V V V V V V — -
* V'' V V — -:
А —_ ш ш. У. Ы: Щ
/Л V V, V Ш ¿г — N
А V V — —- ы
4 щ У —^ —'^ —-- V
1 * -у- Р- — — — — —
№ у — —- — — — -л-
—- V-- ■ЧСЯ, - > —- —= '■¡Г-' Щ
лч Ъг- _= & Хг- V- —= иг- V
г- & ^ '■¿к* V
£ 1 о* — >> х- Ч- 1®= ь N к к •чг- у. V
^^ о. л ■К у: у-- и
-|— •V 1 ч А* > /ч; Л Л Л уА
'л \ = Н V-- -ч- -Ч- 1%: Л, 1 Л А й V- Щ
X • ^ ^ — — Гг- ^ N А А А-
--- V —= —^ —=: Л -А -Г =-
/V И —- - -—г А:
V А т — — — .14 А.
V —— — —~ —= —~ -ЧГ А. ¥ ■А - —=
ух вк- —=- = А-. л
^ — — я® — Л:
к Л ^ к — — --- Л Ш
р Лу ** ■•А-
V ■ш — — -Л ш
■ А- -А — — — — — — — ш ^ £ -
Л. ^ — •Г -к & А А Щ
N — — —= —' — — и
Гч — —- —- — —- —- — — — — --- — р- 'Щ-
-= — — —- -= —= --- — ¿Г" и =р
------- — —= — — — — = — —~ —=■ —- я?
— —- ■у ш —
Рис. 3. Концентрации паров ПГВ (Ъ = 10 мс)
t í пе( ивI - 10.0 TjJibx - 2434. Рпвх - 0.13 2388. 2215. 2065. 1922. 1784. 1650. 1518. 1389. 1261. 1135. 1011. 887.
644.
В
Изотермы в поле течения и положения частиц ПГВ
V V V V У У У У V У у V У и У и V v- У У У y y £ í> V IS y P lé le У ¡e te J¡e te fi- К Л Лч Л = А >1 et e} ef ef ,ef. ef e?
V У у У У У и У Ы ■V У У le У ^ le lé- f ie ¡s ^ ef- ef ef ef
V У V V V У и 1/ У V У V ^ W У is V V У У ef ef ef. ef
V У V V у у V ie и Ш1 V V le y У И У У У У V ef ef- ef
V У V V У У У y к le le le Г # ^ И v У У И Ж ^ ef -ef
V V У V у у y у te y ^ te íe У У ^ V* :У ie te íe ej ef ef
V У V V У у и У ie k. le le te te Ir" У У y V У le te íe
V У V V SS у у le te y te le te te te le У y f У ¡e le Jé p
• у V V У у le у le ts У tr- le te te te te У У У У ie té íe Л- e?
V V V У у у le te y is- te te te te te te У У \У te íe íe
V У У У V le le le fe te ¿e íe ie íe ie le te У У te ie Je
V У V у у у у le íe és- ¿e íe ^ ¿e Je te У У te te J- ¿e
V V V У V к' к' le fe ¡e te íe ¿e ¿e ie V íe ¿e
у V У V ts y íe ie y ^ íe У У y íe J^
V V V У у IfS f ie í'e y ¿e te V У íe
— N| V У V у У ie íe fe JÍ- — — íe у У
Л. --- V V ie y Se y — — V ie
у У Se y Je -i V te
4 4 У? * •i- V y ¿- —= te
É ^ & ¡^ jtr —^ У
г Ъг. t¿- fp Je - —= V
1 у JJ- ы ¿r- jé- -- — -- V ■yi
№ № i % ш --Я» - ~^— — -
\ —=
Рис. 4. Поле течения при h = 0.1 м (t = 10 мс)
tiiiB(irel = 10.0 Т пвх = 2434. Рпвх = 0.13 2.00 1.83 1.69 1.55 1.42 1.29 1.17 1.04 0.92 0.80 0.68 0.56 0.45 0.33 0.21 0.01
Изолинии концентрации паров ПГВ
' V V V' Щ V V V y У У V -y y y V W № Щ1 У y 4": y y te i' f ^ - /4 л у? -г? s e?*
• V V w У У У v V" & У У р1 У /4 л А А -■yí éf-■е? et.. ■ef. e?
' V y w y У У У' У # y Ш y y ir- У y ts y У м 'épr éf ef. éf
' V V y У У У У y. y y u- Щ: к w y У Ш У •Й ef' Щ ef
' V V, V У V y pv y к У к И к k y У: y V ■У ef" ef.
' V v: V y v- v- .У у У Щ Í!^ y У У у k У У y у Ш le ef. sr
' V V y- У У У: У í>: (e У y У ^ к У У у \fe ef, ef
• V V 'Mi У V?'. y:: t« 'tí- te У: í¿r У É'i к y M у У У Ш
' V У V y У y té У y te У У У У У y V У У У У ■í> rf
' V У У' V У У У í^ y •Je y У у • fy- í^ V У у У £f\ У ef-
' V • V :V V :V- V y v^ У У У í^ 'Ш-. '¿Г te te У У ■jfe: y У ^ y íe У y. M ¡e у у у ■-■¡fe
' V ' V y У V •У: У \e y pí f У k' Í-" У ■fe íe У ie y У je- ■Vír te te y У у У: :У щ У ф;
r V < V V V •У К У te íe. íe '¿e íe ^ i ¿e ¿e y V ш У % -
к Ni -íf ¥ * ft- V ¿f V ^ Ф ^ -У? V ^ ^ íe íe 1y W te ¿e J¿» рГ '¿é ¿e & ¿r '¿té. íé ¿e iiS- - — y V у У -У. У-У у ■у ¡gr
■ ' t ! • íc- íe y ■V 1 l С ¿e e*-4 — ^ —"Л — — £ i - i — - >1 ___ __ В
■ 1 V- w 1 — —"
Рис. 5. Концентрация паров ПГВ (t = 10 мс)
Аналогичная картина приведена на рис. 4 и 5 для случая, когда ИАСП расположено вблизи подстилающей поверхности. И в этом случае видно, что предлагаемое устройство способно осуществлять транспортировку достаточно большой массы аэрозольного облака ПГВ в зону горения без заметного разрушения его структуры продуктами сгорания газогенератора.
Список литературы
fl] Coakley T.J. Turbulence modeling for high speed flows // AIAA Paper. 1992. 92-0436.
[2] crow c.T. Review - numerical models for dilute gas-particles flows // Trans, of ASME J. of Fluid Eng. 1982. Vol. 104. P. 297-303.
[3] Yoon S., Jameson A. An LU-SSOR scheme for the Euler and Navier-Stokes equations // AIAA Paper. 1987. 87-600.
[4] Рычков А.Д. Математическое моделирование газодинамических процессов в каналах и соплах. Новосибирск: Наука, 1988. 222 с.
Поступила в редакцию 14 марта 2008 г.