Численное исследование влияния условий распыления воды на температуру в следе «Водяного снаряда»

СПИСОК ЛИТЕРАТУРЫ

1. Флеминге М. Процессы затвердевания. - М.: Мир, 1977. -423 с.

2. Рыкалин Н.Н., Углов А.А., Анищенко Л.М. Высокотемпературные технологические процессы. Теплофизические основы. - М.: Наука, 1985. - 172 с.

3. Крайнов А.В., Швалова Г.В. Основы теплоэнергетики. -Томск: Изд-во Томского политехнического университета, 2011.- 164 с.

4. Логинов В.С., Крайнов А.В., Юхнов В.Е., Феоктистов Д.В., Шабунина О.С. Примеры и задачи по тепломассообмену. 2-е изд., испр. и доп. - СПб.: Изд-во «Лань», 2011. - 256 с.

5. Рассохин Н.Г., Локтионов В.В., Мухтаров Э.С., Орлов В.Е., Молодык К.Е., Ярошенко Н.И. Особенности поведения корпуса реактора ВВЭР-1000 в условиях тяжёлой аварии // Труды II Российской нац. конф. по теплообмену: Пленарные иоб-щие проблемные доклады. Доклады на круглых столах. - М.: Изд-во МЭИ, 1998. - Т. 1. - С. 153-158.

6. Гаврилов П.М., Цыганов А.А. Тепломассообмен и внутренняя самозащищённость ядерных реакторов // Труды II Российской нац. конф. по теплообмену: Пленарные и общие проблемные доклады. Доклады на круглых столах. - М.: Изд-во МЭИ, 1998.- Т. 1. -С. 147-152.

7. Krainov A.V., Shvalova G.V. Heating plants and systems. - Tomsk: TPU Publishing House, 2012. - 115 с.

8. Пасконов В.М., Полежаев В.И., Чудов Л.А. Численное моделирование процессов тепло- и массообмена. - М.: Наука, 1984. - 288 с.

9. Кузнецов ГВ., Крайнов А.В. Исследование сопряжённого теплообмена и гидродинамики при движении вязкой несжимаемой жидкости в каверне прямоугольного типа // Прикладная механика и техническая физика. - 2001. - Т. 42. - № 5. -С. 136-142.

10. Ши Д. Численные методы для решения задач теплообмена. -М.: Наука, 1988. - 544 с.

Поступила 22.01.2013 г.

УДК 536.46

ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ВЛИЯНИЯ УСЛОВИЙ РАСПЫЛЕНИЯ ВОДЫ НА ТЕМПЕРАТУРУ В СЛЕДЕ «ВОДЯНОГО СНАРЯДА»

О.В. Высокоморная, А.О. Марков*, М.Н. Назаров**, П.А. Стрижак, С.Р. Янов

Томский политехнический университет 'ЗАО «Сибирский энергетический научно-технический центр», г. Красноярск “Сибирский федеральный университет, г. Красноярск E-mail: pavelspa@tpu.ru

С использованием нескольких моделей «водяных снарядов» малых размеров выполнено численное исследование макроскопических закономерностей испарения монодисперсной совокупности капель тонкораспыленной воды при прохождении через высокотемпературные продукты сгорания. Установлены зависимости интегральных характеристик тепломассопереноса от отношения площадей испарения и площадей, занимаемых «водяным снарядом». Определены условия распыления воды для эффективного снижения температуры в зоне горения.

Ключевые слова:

Высокотемпературные газы, продукты сгорания, испарение, капли воды, моделирование. Key words:

High temperature gases, combustion products, evaporation, water drops, simulation.

Задачи разработки флегматизаторов горения и расчета параметров их оптимального распыления в зоне горения актуальны в связи со статистическими данными [1] по возрастающим материальным потерям, сопровождающим пожары различной степени сложности (особенно лесные). В таких условиях важную роль играют результаты научных исследований в направлении совершенствования существующих технологий тушения и создания новых. Одной из наиболее перспективных можно считать технологию тушения пожаров тонкораспыленной водой [2-5]. Однако до настоящего времени физические и теоретические основы этой технологии не разработаны в полной мере. Во многом это можно объяснить отсутствием результатов фундаментальных исследований макроскопических закономерностей взаимодействия

тонкораспыленной воды с высокотемпературными продуктами сгорания [4, 5].

Для решения этой проблемы в последние годы выполнен цикл численных исследований процессов тепломассопереноса и фазовых переходов при движении одиночной капли [6, 7], «тандема» [8] и группы из четырех [9] и пяти [10] капель тонкораспыленной воды через высокотемпературные продукты сгорания. Определены зависимости времен испарения (существования) капель от температуры продуктов сгорания, скорости движения, размеров и расстояния между каплями для нескольких систем «водяной снаряд-пламя» [6-10]. Проанализированы условия движения капель в струях тушащей жидкости [11]. Представляет интерес анализ условий распыления воды для эффективного снижения температуры в следе «водяного снаряда».

Цель работы - численное исследование влияния условий распыления «водяного снаряда» на температуру в следе его движения.

Результаты ранее выполненных исследований [6-11] показывают, что основные закономерности движения «водяного снаряда» через высокотемпературные газы можно рассмотреть на примерах двух движущихся последовательно капель (рис. 1, а), двух капель, движущихся параллельно на заданном удалении (рис. 1, б), и пяти капель (рис. 1, в). Установлено [8-11], что если расстояние между каплями больше некоторого предельного (Хп*), то их «совместное» влияние на температуру газов и концентрацию водяных паров в следе «снаряда» будет незначительным. Поэтому при постановке задачи (рис. 1) принималось, что капли движутся со скоростью V на некотором расстоянии друг от друга Хп<Хп*. При записи уравнений, описывающих движение «водяного снаряда», учитывалось влияние сил тяжести капель и сопротивления высокотемпературной газовой среды. Исследования выполнены для продуктов сгорания типичного жидкого (представляющего наиболее пожаровзрывоопасную

группу [1]) топлива - керосина. Начальная температура капель Т0 принималась много меньше температуры газовой среды Т. Считалось, что при движении через газы капли прогреваются за счет теплопроводности и испаряются. Пары воды вдуваются в пристенную область капли и смешиваются с продуктами сгорания. Температура газов в непосредственной близости от траектории движения капель снижается. В условиях интенсивного парообразования размеры капель уменьшаются. Через некоторое время возможно их полное испарение.

Допущения, принятые при численном моделировании в рассматриваемых системах (рис. 1), аналогичны постановкам [6-11].

Важно отметить, что решение сформулированных выше задач в пространственных постановках представляет достаточно длительные вычислительные процедуры. Результатов решения таких пространственных задач применительно к условиям реальных пожаров пока не опубликовано. Этим обусловлен переход к плоским постановкам задач тепломассопереноса, аналогичным [6-11].

Рис. 1. Схема области решения задачи для двух капель, движущихся последовательно (а), двух, движущихся параллельно (б), пяти капель (в): 1 - высокотемпературная газовая смесь, 2 - капли воды

На рис. 2 показано несколько интерпретаций поперечных сечений в рассматриваемом «водяном снаряде» (капли приведены в форме сфер и некоторых кольцевых сегментов).

При решении плоских задач тепломассопере-носа их геометрии будут соответствовать моделям, приведенным на рис. 1. Поэтому в первом приближении при описании процессов тепломассопере-носа в «водяном снаряде» можно использовать рассматриваемые осесимметричные постановки.

Математические модели, соответствующие постановкам задач тепломассопереноса (рис. 1), включают системы нестационарных дифференциальных уравнений в частных производных [6-11], сформулированные в соответствии с основными положениями математической физики [12-14]. Так, например, для модели, приведенной на рис. 1, а, система уравнений в цилиндрической системе координат имеет следующий вид:

0<r<RL, 0^<^, Z2<z<Zз, Z4<z<ZL;

R1<r<RL, Zx<z<Zъ Z3<z<Z4:

ср дТ

57? +1 ВГ\ +

дг 2 г дг

дТ12 дz 2

дС„

дт

= А

дС^ +1 дС^ дг 2 г дг

дС 2 дz 2

, С + Cw - 1;

0<r<R1, Z1<z<Z2, Z3<z<Z4:

С Р дТ2 - Л

С2"2 ~ _/12 дт

дТ22 +1 дТ2 +дТ22 г дг

_ дг2 г дг дz2 .

Здесь т - время, с; т - время полного испарения капель, с; г, z - координаты цилиндрической системы, м; RL, ZL - размеры области решения, м; C - удельная теплоемкость, ДжДкг-К); р - плотность, кг/м3; I- температура, К; Я- теплопроводность, Вт/(м-К); Су - концентрация паров воды; Б - коэффициент диффузии, м2/с; С. - концентрация продуктов сгорания; индексы 1, 2, 3 соответствуют продуктам сгорания, каплям воды, парам воды.

Начальные (т=0) условия (рис. 1, а): 1=10при 0<К^, Z1<z<Z2, Z3<z<Z4; 1=1,, С,=1, Су=0при 0<r<RL, 0^<^, Z2<z<Z3, Z4<z<ZL; R1<r<RL, Z1<z<Z2, Z3<z<Z4.

Граничные условия: на границах «жидкость-высокотемпературная газовая среда» (r=R1, Zl<z<Z2, Zз<z<Z4; *=Д, z=Z2, z=Zз, z=Z4, 0<К^) принимались граничные условия IV рода для уравнений энергии с учетом парообразования, для уравнения диффузии задавались граничные условия II рода с учетом вдува паров воды; на оси симметрии (г=0, 0^<Д) и внешних границах (r=RL, 0^^; ^0, z=ZL, 0<r<RL) для всех уравнений принималось условие равенства нулю градиентов соответствующих функций.

Уравнение движения капель в условиях парообразования с учетом действия сил сопротивления и тяжести имеет следующий вид [15]:

3р3

^ 4р22

где УА - скорость движения капель, м/с; сх - безразмерный коэффициент сопротивления движению капель; Уе - линейная скорость паров воды вблизи торцевых поверхностей капель, м/с; Rd - радиус капель, м; g - ускорение свободного падения, м/с2.

Безразмерный коэффициент сопротивления сх, зависящий в общем случае от конфигурации поверхности тела и его положения относительно направления движения обтекающего потока, определялся согласно [15].

Для вычисления массовой скорости испарения капель воды использовалось следующее выражение [16]:

в

1

-кев ^тПЯгТе /М ’

где в - безразмерный коэффициент конденсации (испарения); кв - безразмерный коэффициент (кв«0,4 [17]); Р - давление насыщенных паров воды, Н/м2; Р - давление паров воды вблизи границ испарения, Н/м2; М - молекулярная масса воды, кг/кмоль; Те - температура капли вблизи границы испарения, К.

Для вычисления давления паров воды использовалось уравнение Клапейрона-Менделеева [12].

Толщина испаряющегося слоя капли определя-Ш т

лась по формуле [16]: Ье - —^.

' Р2

Системы дифференциальных уравнений для рассматриваемых постановок (рис. 1) решались методом конечных разностей [18]. Разностные аналоги дифференциальных уравнений решались локальноодномерным методом [18]. Для решения одномерных разностных уравнений применялся метод прогонки с использованием неявной четырехточечной схемы [18]. Для решения разностных аналогов дифференциальных уравнений с нелинейными граничными условиями применялся метод простых итераций [19]. Для повышения точности решения системы дифференциальных уравнений использовались неравномерные шаги по времени (10-6-10-4 с) и координатной сетке (0,001...0,1 мм). Вблизи границ фазового перехода разностная сетка сгущалась.

Методика оценки достоверности результатов выполненных теоретических исследований, основанная на проверке консервативности применяемой разностной схемы, аналогична используемой в [20-23].

Численные исследования выполнены при следующих типичных для практики значениях параметров [15, 17, 24-27]: начальная температура капель воды 70=300 К; температура продуктов сгорания Тг= 1170 К; тепловой эффект испарения воды (2е=2,26 МДж/кг; размеры капель Д,=0,25 мм, 2==1 мм; расстояние между каплями Хп=0,01—5 мм; размеры области решения Л=0,1 м, 2=2 м; молекулярная масса воды М=18 кг/кмоль; коэффициент сопротивления движению ^=0,98; безразмерный коэффициент конденсации (испарения) в=0,1. Теплофизические характеристики взаимодействующих веществ (капель воды, продуктов сгорания жидкого топлива, паров воды) приведены в [24-27].

Для анализа влияния распределения капель жидкости в «водяном снаряде» на температуру в следе его движения рассмотрены изотермы, полученные ранее [6-11] для систем, приведенных

на рис. 1. В качестве характеристики процесса выбрана температура [6-11] в следе «водяного снаряда» на расстоянии от последней капли, равном пяти ее продольным размерам (Х*=52д). Так, например, на рис. 2 приведены изотермы для системы с двумя движущимися параллельно (рис. 1, б) на расстоянии Ьп=0,3 мм друг от друга каплями и выделено сечение Ь*. Характерный вид изотерм в окрестности оси симметрии при некотором увеличении (рис. 3) иллюстрирует выполнение сформулированных выше граничных условий (равенство нулю градиентов соответствующих функций).

Выполнен анализ распределения температуры газопаровой смеси в следе «снаряда» в сечении Ь* (рис. 3). Следует отметить, что аналогичный анализ можно выполнить для изолиний концентраций водяных паров С¥ и продуктов сгорания С, приведенных в [6-11]. Однако больший интерес представляет исследование изотерм. Это обусловлено тем, что согласно представлениям классической теории диффузии и теплопередачи в химической кинетике [14] температура в зоне реакции играет определяющую роль. Зависимость скорости реакции окисления от температуры - экспоненциальная, а от концентрации - степенная [14]. В таких условиях при тушении пожаров и возгораний целесообразно, в первую очередь, снизить температуру в зоне реакции до предельной, при которой рассматриваемая реакция затухает.

Для оценки эффективности снижения температуры в следе «снаряда» сформулированы два безразмерных параметра - 71* и 51*. Параметр 71* представляет отношение температуры в следе «водяного снаряда» (на расстоянии Ь* от последней капли) 7 к начальной температуре продуктов сгорания Т( (71*= 71/7). Параметр 5\* - это отношение площади испарения £е к площади зоны (в форме цилиндра с высотой Ь* и радиусом Л*) в следе «водяного снаряда» Бг (5*=^/^), в которой температу-

ра смеси газов и паров воды не превышает Т2. Так, например, на рис. 2 в следе капель показано характерное расстояние Ь* от последней капли. Расстояния от оси симметрии (/=0) до точек пересечения линии Ь* с изотермами представляют радиусы Я* зон в следе капель, в которых температуры не превышают значений, соответствующих пересеченным изотермам. Например, на рис. 2 в точке пересечения изотермы Т=800 К и линии Ь* определен радиус Я*=0,72 мм. Для этого радиуса рассчитывается соответствующее значение площади £2 (зоны в форме цилиндра с высотой Ь*=52й и радиусом Я*=0,72 мм).

Для вычисления площадей £е и £2 использовались соотношения:

п о

5е = £ (2пЯйЬ + 2пя|);., 5; = 2пЯ * Ь * +2пЯ *2,

I =1

где г - номер капли; п - количество капель; Я* и Ь* - характерные размеры области в следе капель с температурой Т<Т2.

На рис. 5 приведены зависимости Т1*=/(51*) для систем с двумя каплями, движущимися последовательно на расстоянии Ьп=5 мм (рис. 1, а), двумя каплями, движущимися параллельно на расстоянии Ьп=0,3 мм (рис. 1, б), и пятью каплями (рис. 1, в) при Ьп=2 мм. Значения Ьп выбраны исходя из результатов ранее [6-11] выполненных исследований, как наиболее типичные. При таких Ьп «совместное» влияние капель на время существования «водяного снаряда» значительно [6-11]. Можно отметить достаточно существенно отличающиеся диапазоны изменения 5\* (рис. 4) для рассматриваемых систем (рис. 1) при близких значениях Ьп и равных размерах капель.

Для зависимостей Т1*=/(51*), представленных на рис. 4, сформулирована группа аппроксима-ционных выражений (0,1<£1*<0,4):

• Т1*=3,698-14,282£1*+16,223£1*2 - для двух движущихся последовательно капель (рис.1, а);

• Т1*=1,069-3,257£1*+5,087£1*2 - для двух движущихся параллельно капель (рис.1, б);

• Т1*=9,289-48,793£1*+68,146£1*2 - для пяти капель (рис.1, в).

Выполненные численные исследования позволили установить (рис. 5), что наименьшие площади испарения $ для обеспечения максимальных площадей £2 с температурой Т1* в следе «водяного снаряда» соответствуют двум движущимся параллельно каплям (рис. 1, б). Максимальные значения ^ соответствуют системе с двумя последовательно перемещающимися каплями (рис. 1, а). Для пяти капель (рис. 1, в) характерны средние значения (рис. 5).

Установленные особенности можно объяснить тем, что при использовании пяти капель характерный радиус Я* больше, чем для двух капель. Однако вместе с тем и суммарная площадь испарения для пяти капель ^ существенно больше, чем для двух. Как следствие, отношения для системы с пятью каплями больше, чем для двух параллельно движущихся капель. В тоже время значения £*для пяти капель меньше, чем для двух последовательно движущихся капель. Это можно объяснить тем, что для двух последовательно движущихся капель характерный размер (Я*) области с заданной температурой Т1* существенно меньше, чем для систем, приведенных на рис. 1, би в.

Для определения расстояний между каплями (рис. 1), оптимальных с точки зрения эффективного снижения температуры продуктов сгорания, выполнены численные исследования влияния параметров Ьп на температуру газов в следе капель. Для системы с двумя параллельно движущимися каплями (рис. 1, б) параметр Ьп варьировался в диапазоне от 0,01 до 0,6 мм. В случае двух последовательно движущихся (рис. 1, а) и пяти (рис. 1, в) капель диапазон варьирования этого параметра составил Ьп=0,01...5 мм. На рис. 6 приведены зависимости Т2*=/(^2*) для рассматриваемых систем (рис. 1).

Параметр Т2* представляет отношение температуры в следе капли на расстоянии Ь* от последней капли (на оси симметрии - /=0) к температуре Тг. Параметр £,* используется для иллюстрации отно-

шения площади испарения «водяного снаряда» £е к занимаемой им площади Sw (характеризует плотность распределения капель в «водяном снаряде»). Для приведенных на рис. 1 постановок задач те-пломассопереноса площади Sw можно определить по формулам: 8^=2жК1(2Ьі+Ьп)+2жЯ12 - система с двумя последовательно движущимися каплями (рис.1, а); 5¥=2пЛ2Хй+2пЛ22- система с двумя параллельно движущимися каплями (рис. 1, б); 5,¥=2пЛ2(3Хй+2Хп)+2л:Л22 - система с пятью каплями (рис. 1, в).

Для зависимостей Г2*=/(52*), приведенных нарис. 6, сформулированы следующие аппрокси-мационные выражения (0,3<52*<0,9):

• Т2*=0,644-0,50752*+0,29152*2 - для двух движущихся последовательно капель (рис. 1, а);

• Т2*=1,444-2,78152*+1,89452*2 - для двух движущихся параллельно капель (рис. 1, б);

• Т2*=0,648-0,42852*+0,25552*2 - для пяти капель (рис. 1, в).

Установленные зависимости Т2*=/(£2*) показывают, как при снижении плотности распределения капель воды (уменьшении 52*) температура в следе «водяного снаряда» возрастает. При этом в случае достаточно больших значений характерных удалений Ьп (для двух параллельно движущихся капель при Ьп>2Яй) температура в следе несущественно снижается относительно Т£ (рис. 6).

При повышении S2* вне зависимости от геометрии расположения капель относительно друг друга температура в следе «водяного снаряда» существенно снижается относительно Т{ (рис. 6). Это особенно хорошо видно на примере системы с двумя параллельно движущимися каплями (кривая 2). При превышении параметром 52* значения

0.5.отклонения Т2* для кривых 2 и 3 не превышают 3 %. Дальнейший рост 52* не приводит к заметному изменению Т2* (рис. 6). Установленная закономерность позволяет сделать вывод о том, что при распылении «водяного снаряда» с обеспечением значения параметра 52* не меньше 0,5 могут быть достигнуты условия, достаточные для снижения температуры в зоне горения без дальнейшего повышения S2*.

В тоже время важно отметить, что для систем с двумя последовательно движущимися каплями (рис. 1, а) и пятью каплями (рис. 1, в) изменения значений Т2* с ростом 52* не так существенны, как для двух параллельно движущихся капель (рис. 1, б). Этот результат показывает, что на линиях, соответствующих траекториям движения капель, температура в следе снижается не так значительно при уменьшении Ьп, как в окрестности этих линий. Это достаточно хорошо показано в работах [6-11] и видно, в частности, на рис. 2.

Установленные зависимости Т1*=/(51*) и Т2*=/(52*) показывают, что для снижения температуры в зоне горения (не только на оси симметрии /=0 или внешнем контуре тела вращения с радиусом /=Я*) целесообразно, скорее всего, фрагментарное распыление флегматизатора. На рис. 5 показано, что наименьшие температуры Т1* в следе «снаряда» при малых значениях 51* соответствуют системе с двумя параллельно движущимися капля-

СПИСОК ЛИТЕРАТУРЫ

1. Акинин Н.И., Булхов Н.Н., Гериш В.А. Статистический анализ причин аварий и травматизма на опасных производственных объектах // Пожаровзрывобезопасность. - 2010. - № 10. -С. 53-55.

2. Карпышев А.В., Душкин А.Л., Рязанцев Н.Н. Разработка высокоэффективного универсального огнетушителя на основе генерации струй тонкораспыленных огнетушащих веществ // Пожаровзрывобезопасность. - 2007. - № 2. - С. 69-73.

3. Гаев Д.В., Ершов А.В., Прохоров В.П., Карпышев А.В., Душкин А.Л. Система противопожарной защиты салона вагона метрополитена на базе высоких технологий // Проблемы безопасности и чрезвычайных ситуаций. - 2009. - № 3. - С. 67-72.

4. Душкин А.Л., Ловчинский С.Е. Взаимодействие пламени горючей жидкости с тонкораспыленной водой // Пожаровзры-вобезопасность. - 2011. - № 11. - С. 53-55.

5. Душкин А.Л., Карпышев А.В., Ловчинский С.Е. Особенности распространения жидкостной струи в атмосфере // Пожаров-зрывобезопасность. - 2011. - № 12. - С. 45-48.

6. Волков Р.С., Кузнецов Г.В., Стрижак П.А. Численная оценка оптимальных размеров капель воды в условиях ее распыления средствами пожаротушения в помещениях // Пожаровзрыво-безопасность. - 2012. - № 5. - С. 74-78.

7. Глушков Д.О., Кузнецов ГВ., Стрижак П.А. Численное исследование тепломассопереноса при движении «тандема» капель

ми. На рис. 6 видно, что минимальные значения температуры Т2* в следе «снаряда» при варьировании 52* в достаточно широких диапазонах зафиксированы для двух последовательно движущихся капель. Можно сделать вывод о том, что для увеличения размеров зоны в следе «водяного снаряда», в которой температура может быть снижена относительно Т{ до некоторой заданной, необходимо применение «снарядов» с последовательно-параллельной (фрагментарной) системой распыления (например, как показано на рис. 1, в).

Выводы

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

Работа выполнена при финансовой поддержке Министерства образования и науки РФ (ФЦП«Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы, мероприятие 1.4, соглашение № 14.B37.21.2071).

воды в высокотемпературной газовой среде // Тепловые процессы в технике. - 2012. - № 12. - С. 531-538.

8. Стрижак П.А. Численное исследование условий испарения совокупности капель воды при движении в высокотемпературной газовой среде // Пожаровзрывобезопасность. - 2012. -№ 8. - С. 26-31

9. Волков Р.С., Высокоморная О.В., Стрижак П.А. Численное исследование условий взаимодействия диспергированного флег-матизатора горения с высокотемпературными продуктами сгорания // Безопасность труда в промышленности. - 2012. -№ 10. - С. 74-79.

10. Андреев Г.Г., Глушков Д.О., Панин В.Ф., Стрижак П.А. Тепло-массоперенос при взаимодействии диспергированного флег-матизатора горения с высокотемпературными продуктами сгорания // Бутлеровские сообщения. - 2012. - № 8. - С. 86-94.

11. Стрижак П.А. Численный анализ процесса испарения капли, движущейся в струе воды через высокоскоростные продукты сгорания // Пожаровзрывобезопасность. - 2012. - № 9. -С. 17-23.

12. Кутателадзе С.С. Основы теории теплообмена. - М.: Атомиз-дат, 1979.

13. Пасконов В.М., Полежаев В.И., Чудов Л.А. Численное моделирование процессов тепло- и массообмена. - М.: Наука, 1984.

14. Франк-Каменецкий Д.А. Диффузия и теплопередача в химической кинетике. - М.: Наука, 1987. - 490 с.

15. Лойцанский Л.Г Механика жидкости и газа. - М.: Наука, 1970.- 904 с.

16. Полежаев Ю.В., Юрьевич Ф.Б. Тепловая защита. - М.: Энергия, 1976. - 391 с.

17. Муратова Т.М., Лабунцов Д.А. Кинетический анализ процессов испарения и конденсации // Теплофизика высоких температур. - 1969. - № 5. - С. 959-967.

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

19. Коздоба Л.А. Методы решения нелинейных задач теплопроводности. - М.: Наука, 1975. - 227 с.

20. Кузнецов Г.В., Стрижак П.А. Зажигание накаленной одиночной частицей жидких углеводородных топлив // Известия Томского политехнического университета. - 2008. - Т. 312. -№ 4. - С. 5-9.

21. Кузнецов Г.В., Стрижак П.А. Численное решение задачи воспламенения жидкого пожароопасного вещества одиночной «горячей» частицей // Физика горения и взрыва. - 2009. -№ 5. - С. 42-50.

22. Kuznetsov G.V., Strizhak P.A. The influence of heat transfer conditions at the hot particle-liquid fuel interface on the ignition characteristics // Journal of Engineering Thermophysics. - 2009. -№ 2. - P. 162-167.

23. Kuznetsov G.V., Strizhak P.A. Transient heat and mass transfer at the ignition of vapor and gas mixture by a moving hot particle // International Journal of Heat and Mass Transfer. - 2010. - V. 53. -Iss. 5-6. - P. 923-930.

24. Варгафтик Н.Б. Справочник по теплофизическим свойствам газов и жидкостей. - М.: ООО «Старс», 2006. - 720 с.

25. Теплотехнический справочник: в 2 т. Т. 1 / под ред. В.Н. Юре-нева, П.Д. Лебедева. - М.: Энергия, 1975. - 743 с.

26. Теплотехнический справочник: в 2 т. Т. 2 / под ред. В.Н. Юре-нева, П.Д. Лебедева. - М.: Энергия, 1975. - 896 с.

27. Корольченко А.Я. Пожаровзрывоопасность веществ и материалов и средства их тушения: справочник. Ч. 1. - М.: Пожнаука, 2004. - 713 с.

Поступила 19.12.2012 г.

УДК 621.18

ИССЛЕДОВАНИЕ СТРУКТУРЫ ТЕЧЕНИЯ В НОВОЙ КОНСТРУКЦИИ ВИХРЕВОЙ ТОПКИ МЕТОДОМ ЦИФРОВОЙ ТРАССЕРНОЙ ВИЗУАЛИЗАЦИИ

И.С. Ануфриев, Г.В. Кузнецов*, Е.Ю. Шадрин, О.В. Шарыпов, Ю.А. Аникин, С.В. Двойнишников, Д.В. Куликов, В.В. Рахманов

Институт теплофизики СО РАН, г. Новосибирск *Томский политехнический университет E-mail: anufriev@itp.nsc.ru

Исследована структура изотермического закрученного течения в модели вихревой топки новой конструкции с горизонтальной осью вращения и распределенным тангенциальным вводом топливовоздушных струй с помощью метода цифровой трассерной визуализации. Измерены осредненные характеристики поля скорости в ряде сечений. Проанализированы особенности внутренней аэродинамики новой конструкции вихревой топки парогенератора при различных режимных параметрах.

Ключевые слова:

Турбулентные закрученные потоки, вихревая топка, тангенциальный ввод горелочных струй, метод цифровой трассерной визуализации, перспективное оборудование теплоэнергетики.

Key words:

Turbulent swirl flows, vortex furnace, tangential Input of burner streams, particle image velocimetry, promising energy equipment.

При проектировании и разработке новых аппаратов и технологий теплоэнергетики ключевым этапом является физическое моделирование, которое позволяет при низких затратах изучить основные физические процессы в элементах энергетического оборудования с целью оптимизации конструктивных и режимных параметров разрабатываемых устройств [1, 2].

В данной работе исследуется вихревое топочное устройство новой конструкции с горизонтальной осью вращения и рассредоточенным по периметру тангенциальным вводом топливовоздушных струй (патент РФ № 2042084). Основными отличительными особенностями новой конструкции вихревой топки (по сравнению с известной топкой Н.В. Голованова [3]) является: дополнительный тангенциальный подвод топлива, расположенный в ни-

жней части камеры сгорания, и увеличенная ширина горловины диффузора. Наличие в исследуемой конструкции вихревой топки рассредоточенных по периметру (условной окружности камеры горения) горелочных струй, ориентированных в противоположном направлении, обеспечивает гибкость управления структурой течения и режимными параметрами, а горизонтальная ось вращения потока увеличивает полноту выгорания топлива.

Цель данной работы - экспериментальные исследования структуры закрученного течения в новой конструкции вихревой топки. Измерение аэродинамических характеристик пространственного турбулентного течения проводилось в изотермической лабораторной модели, изготовленной из оргстекла толщиной 10 мм (рис. 1, а). Основными эл-