УДК 614.843.8
МЕХАНИКА ДВИЖЕНИЯ ДВУХФАЗНОГО ПОТОКА В ДИФФУЗОРЕ ОРОСИТЕЛЯ С ПРЕДВАРИТЕЛЬНЫМ ГАЗОНАСЫЩЕНИЕМ
Качанов И.В.*, д.т.н., профессор, Карпенчук И.В.**, к.т.н., доцент, Павлюков С.Ю.** *Белорусский национальный технический университет **Командно-инженерный институт МЧС Республики Беларусь
*e-mail: [email protected] **e-mail: [email protected]
Получены основные уравнения движения двухфазного потока, которая позволяет рассчитывать распределение осредненного давления в диффузоре по длине, а также определять потери давления в диффузоре, что дает возможность проектировать выходной участок гидравлической системы с гидродинамическим сопротивлением, обеспечивающим работу инжектора в заданном режиме.
Governing equations of motion of two-phase flow, which allows to calculate the distribution of the mean pressure in the diffuser length, as well as to determine the pressure loss in the diffuser, which makes it possible to design the output section of the hydraulic system of hydro-dynamic resistance, offers an injector in a given mode.
(Поступила в редакцию 30 июня 2012 г.)
ВВЕДЕНИЕ
Установки автоматического пожаротушения являются одним из самых эффективных средств пассивной противопожарной защиты. Применение таких систем, где огнетушащее вещество при возгорании автоматически подается в защищаемое помещение, позволяет ликвидировать пожар на ранней стадии. В качестве огнетушащего вещества в таких установках используются: вода, пена, газ, порошок, аэрозоль.
Наибольшее распространение как в Беларуси, так и за рубежом получили автоматические установки пожаротушения, в которых в качестве огнетушащего вещества используются вода или пена. В Республике Беларусь доля пенных установок в общем объеме составляет около 15 %.
Установки пенного пожаротушения применяются для защиты промышленных предприятий, складов, нефтехранилищ, транспорта и т. д. Это обусловлено такими достоинствами тушения с применением пены как: существенное сокращение расхода воды, тушение пожаров больших площадей, объемное тушение, повышенная (по сравнению с водой) смачивающая способность. Кроме того, при тушении пеной не требуется одновременное перекрытие всего зеркала горения, поскольку пена способна растекаться по поверхности горящего материала.
К числу свойств пен, которые связаны с их огнетушащей способностью, относятся: кратность, стойкость, вязкость и дисперсность.
Кратностью пены называют отношение объема пены к объему раствора, из которого она получена. Кратность показывает, сколько газа содержится в пене, полученной из определенного объема жидкости. Чем больше кратность, тем легче пена и тем больше заключено в ней газа. Пены бывают низкой (менее 10), средней (10-100) и высокой (более 100) кратности.
Средняя интенсивность разрушения пены показывает, насколько уменьшается единица объема пены за единицу времени. Для оценки качества пены используется обычно не средняя интенсивность разрушения, а величина, обратная ей, которая называется стойкостью пены.
Вязкость характеризует способность пены к растеканию по поверхности. Это свойство весьма важно, так как при тушении пламени пена в ряде случаев должна пройти весьма большое расстояние от места подачи.
Дисперсность пены дает представление о степени измельченности пузырьков пены и определяется их размером. Пена полидисперсна, т. е. состоит из пузырьков различных размеров. При разрушении пены дисперсность ее меняется. Изменение размеров пузырьков протекает неодинаково в разных местах слоя пены. На поверхности пены происходит более быстрое укрупнение пузырьков, чем в глубинных слоях.
В настоящее время для получения пены низкой кратности в автоматических установках пенного пожаротушения наиболее широко применяются (до 90 % от общего числа пенных оросителей) универсальные пенные оросители. Повышение эффективности пенных оросителей и эффективности пожаротушения может быть достигнуто предварительной аэрацией (газонасыщением) огнетушащего вещества, что необходимо для повышения кратности и дисперсности пены - основных параметров, влияющих на пожаротушение.
Проведенные в начале 80-х годов работниками ВНИИПО МВД России (И.Ф. Безродный, В.А. Меркулов, В.М. Кучер, В.Ч. Реутт) исследования по изучению основных характеристик пены, влияющих на огнетушащую эффективность, позволяют сделать вывод, что одной из важных характеристик пены, влияющей на скорость ее разрушения на поверхности горючего, является дисперсность. По результатам проведенных опытов установлено, что при малой плотности теплового потока скорость разрушения пены прямо пропорциональна размеру пузырьков в пене, а при большой - не зависит от их размера. Гидродинамические параметры пенного слоя зависят от дисперсности. Уменьшение размеров пузырьков приводит к увеличению огнетушащей эффективности только до некоторого предела, при котором еще возможно хорошее растекание до горящей поверхности. Принимая во внимание то обстоятельство, что скорость разрушения пены зависит от плотности теплового потока, можно утверждать, что огнетушащая эффективность пены при неразвившихся пожарах (пожарах с тепловым потоком малой плотности), будет напрямую зависеть от дисперсности пены. Учитывая, что автоматические установки пожаротушения предназначены для обнаружения и локализации пожара на ранней стадии, т. е. на той стадии, когда плотность теплового потока будет невелика, увеличение дисперсности пены в данных установках повысит их эффективность.
Дисперсность определяется типом пеногенерирующей аппаратуры. Необходимо подчеркнуть, что до настоящего времени задача создания отечественного оросителя с предварительной аэрацией огнетушащего вещества для автоматических установок пенного пожаротушения не решена. Это можно объяснить тем, что рабочий процесс таких оросителей детально не изучен, он также и не регламентирован действующими техническими нормативно-правовыми актами. В этой связи представляет интерес провести исследования механики движения жидкости в данных оросителях.
УРАВНЕНИЯ ДВИЖЕНИЯ ДЛЯ ДВУХФАЗНОГО ПОТОКА В ДИФФУЗОРЕ
Сформулируем основные уравнения движения двухфазной жидкости в вертикальном инжекторе. Воспользуемся представлением о движении двухфазных сред, предложенном в [1, 2]. Совместное движение газа и жидкости многообразно. Крайними формами является движение двух сплошных параллельных потоков, взаимодействующих только по одной непрерывной поверхности раздела, и движение потока пены, в которой обе фазы образуют сложную и неустойчивую структуру. Специфической особенностью рассматриваемой среды является тот факт, что даже в случае, когда обе фазы практически можно считать несжимаемыми, газожидкостная система ведет себя как сжимаемая жидкость.
Схема течения двухфазной жидкости в диффузоре инжектора представлена на рис. 1.
X
Рисунок 1 - Схема течения двухфазной жидкости в диффузоре инжектора
Как показано на рисунке, диффузор расположен вертикально. Обозначим площадь произвольного сечения диффузора инжектора через £.
Известно, что газосодержание ф определяется по формуле [1, 2]:
1 +Лt
ф(Х0>^ )= [ ф'^^ ), (1)
0
где ф'(х0, ^) - доля площади £ , занятая газообразной фракцией в момент времени 1,
Дt >> 1/пср. (пср. - средняя частота прохождения отдельных образований газообразной фракции через сечение £ ), поэтому в дальнейшем будем
рассматривать двухфазный поток в диффузоре как движение сжимаемой жидкости.
Определим силы, действующие на элементарный участок диффузора, заключенный между сечениями £х и £х + дх, а также их проекции на ось Ох.
Массовые силы:
x + Ax
G = g JJJp = g J d^ JP дф-djdz, (2)
x Sx
где g - ускорение свободного падения;
рдф. - плотность двухфазной жидкости (рдф. = рдф.(х, р );
V - объем диффузора между сечениями 8х и 8х + Дх; - координата промежуточного сечения между сечениями 8Х и 8Х + Дх.
При вертикально расположенном инжекторе массовые силы составят:
х+Дх х+А
П-122-
2
G = g J РДф^dS = g J п-tg2^Чф^, (3)
где ад - угол конусности диффузора. Силы трения:
x+Ax
T =JJ TdS = J T—J ezdl, (4)
Sr. x о ад I
x cos—— k
где S6oK. - боковая поверхность усеченного конуса между сечениями Sx и Sx + д*; т - касательные напряжения на стенке диффузора; /^ - длина окружности радиуса r = ^tg(a^2);
- вектор единичной длины, сонаправленный с Т .
Jexdl = -2 - п - e0 - cosО^ - tg- S, (5)
где е 0 - вектор единичной длины, сонаправленный с осью Ох. В проекции на ось Ох сила трения составит:
x+A a
T = - J 2-п-tg-2--T(6)
x
При определении сил реакции стенок имеем:
N = JJ adS, (7)
где о - нормальные напряжения на стенке диффузора.
В проекции на ось Ох аналогично, как и в случае определения сил трения, получим:
Sбoк
х+Ах
N = | 2• п• tg2^• £•ст^.
(8)
Определяем силы давления, которые действуют на площади сечения 8х и 8х + дх, ограничивающие участок Ах:
х+Ах
р = Р (х, г )• ^ - Р (х + Ах, г )• ^ А+Дх • [ Р (х + А г ) - Р (^ г )] + р (^ г) • ( - ^х+Ах )
(9)
р = -п •tg2 ^ •
х+Ах
(х + Ах)2 • | + Р •( х* Ах + (Ах)22
гдер - давление в перпендикулярных оси Ох сечениях.
Сделаем допущение, что проскальзывание фаз отсутствует, т. е. скорости фаз равны. Для элементарного участка диффузора, заключенного между сечениями 8х и 8х + ах (т. е. Ах << 1, при этом заменяем Ах на dx), запишем закон сохранения импульса. Для этого используем формулы (3)-(9). Учтем, что для непрерывной функции
х+Ах
| /(£Щ = /(х)ёх
(10)
(с точностью до членов 2-го порядка малости). Пренебрегая членами порядка малости больше единицы, запишем изменение импульса для массы двухфазной жидкости от сечения и 5х + Ах при пренебрежении обменом массы между фазами:
о о
g • п • tg • х •рдф •ёх - 2 • п • tg^ • х • т • dx +
+2 • п • tg2 • х •ст • ёх - п • tg2
V
х2 • — • ёх + 2 • х • Р • ёх дх
у
= ((. + ¿Чф.)) + ^)- т- • V
дф.
(11)
где тдф. - массовый расход двухфазной жидкости; V - скорость двухфазной жидкости.
По Релею [3] вследствие быстротечности процесса через поверхность пузырьков отток массы практически не происходит.
Для преобразований используем следующие соотношения:
2 а
тдф. = п - .ё 2
- х
^ С +
д, дху дф.
V - х
1 дv + дv йх V д, дх
(12)
дф = 0 (постоянство суммарного массового расхода),
Пренебрегая величинами 2 и 3-го порядка малости, из уравнения (11) после преобразования получим:
g - Рдф.
+
2 -а 2 - т др 2 - р ад дх х
х
рдф.-
+V
х-
Фдф. 1 д
дv дv
— + V—
д, дх
+
У
+^ -дх (рдф.- х 2 -—
д, х дх
(13)
К полученному уравнению добавим уравнение неразрывности, в данном случае примет следующий вид:
которое
2 ад
п - .В2 - ■
х
2 дРдф. 5
+ дх ((2 - рдф.- — )
д, дх
0.
(14)
В случае установившегося движения уравнения (13) и (14) примут следующий вид:
g - Рдф. "
2 - т + 2 - а йр 2 - р
а х йх х
х •
й— 1 й / 2 \ рдф.-—-йх +—-^-йх(рдф.'х - —);
(15)
х йх
2 а д 2
п-.ё у-Рдф.- — -х
ооиб. = т
дф.'
(16)
В результате получена незамкнутая система двух уравнений с пятью неизвестными (т, о, р, —, рдф), описывающих движение газо-жидкостной смеси в диффузоре инжектора.
Далее рассмотрим случай установившегося движения жидкости в диффузоре инжектора с малыми углами конусности (а < 15°).
Поскольку диффузор расположен вертикально и рабочие давления много больше гидростатического давления в сечении диффузора, будем считать давление с сечении диффузора одинаковым [4]. Вследствие этого нормальные напряжения на стенке можно принять равным значению давления в сечении. Тогда уравнение (15) примет следующий вид:
g • РдФ.
2-т
- tg -
dp dv
dx = vdx+v
d_
x ^^x
(рдф.- X2- V)
(17)
или
g • Рдф.
2-т
x- tg f
dp dv 2рдф.
= рдф.-v— + ф
dx
dx
x
(18)
Полученные уравнения осредненного по времени движения газожидкостной смеси в диффузоре позволяют рассчитывать распределение осредненного давления в диффузоре по его длине, а также определять давление на ее границе в конце диффузора [5].
ЗАКЛЮЧЕНИЕ
Система уравнений (16) и (18) описывает движение двухфазного потока в диффузоре инжектора. Решение этих уравнений позволит в дальнейшем определять потери давления в инжекторе, его геометрические характеристики и гидродинамические параметры. Для элементарного участка диффузора, заключенного между сечениями, записан закон сохранения импульса и уравнение неразрывности. Значение газосодержания задается в зависимости от кратности пены.
Таким образом, получены основные уравнения движения двухфазного потока, отличающиеся неравномерным характером движения газожидкостной смеси в диффузоре инжектора, которые позволят рассчитывать распределение осредненного давления в диффузоре по длине, а также определять потери давления в диффузоре, что даст возможность проектировать выходной участок гидравлической системы с гидродинамическим сопротивлением, обеспечивающим работу инжектора в заданном режиме.
2
ЛИТЕРАТУРА
1. Кутателадзе, С.С. Гидродинамика газожидкостных систем / С.С. Кутателадзе, М.М. Стырико-вич. - М. : Энергия, 1976. - 296 с.
2. Уравнения движения кавитационного двухфазного потока в диффузоре пеносмесителя ПС-5 / И.В. Карпенчук [и др.] // Чрезвычайные ситуации: предупреждение и ликвидация. - 2005. -№ 7 (17). - С. 154-160.
3. Rayleigh, O.M. On the pressure developed in the collapse of a spherical cavity. Phyl. Mag., 1917, V. 34, N 200. - Р. 94-98.
4. Рабинович, Е.З. Гидравлика / Е.З. Рабинович. - М. : Недра, 1980. - 278 с.
5. Павлюков, С.Ю. Уравнения движения двухфазного потока в диффузоре инжектора-аэратора для газонасыщения огнетушащего вещества в автоматических установках пожаротушения / С.Ю. Павлюков, С.П. Ерома, В.В. Пармон // Обеспечение безопасности жизнедеятельности: проблемы и перспективы: сб. материалов VI Междунар. науч.-практ. конф. В 2 ч. Ч. 1. - Мн. : КИИ, 2012. - С. 224-226.