УДК 532:623.8/9
О МЕХАНИЗМАХ ФАЗОВЫХ ПРЕВРАЩЕНИЙ В КАПЛЯХ ВОДОТОПЛИВНОЙ ЭМУЛЬСИИ
Исаков А.Я., - д.т.н., профессор
Камчатский государственный технический университет
Рассмотрены физические особенности возникновения конкурентной фазы при нагревании капель водотопливной эмульсии. Показано, что процесс протекания микровзрыва капель определяется объёмным содержанием воды в топливе и дисперсностью эмульсии.
The article describes physical peculiarities of competitive phase which emerges when drops of water-fuel emulsion are heated. It is shown that the process of micro explosion in drops is determined by the amount of water contained in the fuel and by dispersity of emulsion.
Ключевые слова: МЕХАНИЗМ ФАЗОВЫЕ ПРЕВРАЩЕНИЯ КАПЛЯ ВОДОТОПЛИВНАЯ ЭМУЛЬСИЯ
Использование водотопливных эмульсий для питания энергетических установок позволяет помимо улучшения их экономических показателей существенно снижать уровень вредных выбросов в атмосферу. Наличие в штатном жидком топливе микроскопических включений воды принципиально изменяет процесс воспламенения в камерах сгорания распыленных капель [1]. При определённой степени нагревания водяные включения за счёт образования конкурентной паровой фазы увеличивают сой объём, разрушая исходную каплю топлива на более мелкие фрагменты, что благоприятно сказывается на процессе горения. Таким образом, основной положительный эффект, вносимый присутствием воды, обусловлен изменениями её фазового состояния. В этой связи, исследования процессов, предшествующих воспламенению капель эмульсии представляются актуальными, как с теоретических, так и с практических позиций.
Образование конкурентной фазы при перегреве воды может протекать по двум различным схемам, гомогенной и гетерогенной. Гомогенная схема
образования конкурентной фазы применима для чистых жидкостей, потому что предполагает нарушение сплошности вследствие обеспечения условий, при которых две соседние молекулы могут оказаться на расстоянии примерно в два и более раза превышающие их эффективные диаметры. Для каждой жидкости существуют вполне определённые температуры, при которых происходит разрушение сплошности. Для воды при эффективном радиусе молекул Го @ 310 - 10 м максимальное растягивающее напряжение, исходя из молекулярных соображений, определится как
20(1) @ 2К -°,15(т-288)].Щ-3 @ 5,з.108Па, 0)
Го Го
где о0 - коэффициент поверхностного натяжения на границе раздела вода-воздух при нормальной температуре, Т - текущее значение абсолютной температуры, г0 - эффективный радиус молекулы воды. Все реальные жидкости, включая воду, демонстрируют на практике значительно меньшие прочности на разрыв, вследствие того, что не являются сплошными средами. В дистиллированной воде, подвергнутой специальной обработке высоким давлением, и многоступенчатой фильтрации, по не вполне объяснимым в настоящее время причинам, продолжают существовать ядра конкурентной фазы в виде микрополостей, заполненных парами воды и газами, преимущественно воздухом.
Рассмотрим физические процессы, предшествующие воспламенению по-лидисперсной кали водотопливной эмульсии, изображённой на рис.1. В сферической капле топлива радиусом Яр равномерно по всему её объёму распределены включения воды примерно одинакового радиуса Я^, концентрация которых в пределах капли может быть принята постоянной. Развиваемая далее термофизическая модель процессов основывается на гетерогенной схеме возникновения конкурентной фазы в перегретой воде, предполагающей наличие ядер, в виде парогазовых микрополостей, которые при определенной критической температуре Тсг теряют устойчивость и
начинают увеличивать свой объём за счёт интенсивного испарения в них паров воды. Динамика развития ядер конкурентной фазы определяется уравнением [2]
Рис. 1. Полидисперсная капля водотопливной эмульсии
р*+р*
2о
Я
= {(я),
(2)
где р8, рё, р0 - соответственно давление насыщенных паров, газа и внешнее давление, о - коэффициент поверхностного натяжения на границе пар -вода. Ядро, содержащее пар и газ окажется неустойчивым при условии
»0. (3)
<ж
Совмещая уравнение (2) и неравенство (3) получим уравнение для критического значения радиуса ядра Ясг
_з| 2[а0 - 0,15(Тсг - 288)]-10-3 р ( ) р| 32 [а0 - 0,15(Тсг - 288)]-10-9
М КГ Р.(Т-)-Р0}-27 Т[ (Т)-Р0 ] • (4)
где Т;;, Т - начальное, критическое и текущее значение температуры, а0 -коэффициент поверхностного натяжения при температуре Т = 273 К.
Оценка количества воды, переходящей из жидкого состояния в парообразное состояние может быть произведена с помощью функции распределения ядер конкурентной фазы по размерам, которая принимается в виде
[3, 4]
Н(К)- 1,ззр(тах- ятт )з - я3, (5)
-тах?
где а - объёмное содержание в воде газов в свободном состоянии, К К-тш - максимальный и минимальный радиус ядер. Постоянная А определялась нами экспериментально термоакустическим методом [5].
В основу методики измерений концентрации ядер конкурентной фазы и распределения их по размерам положен анализ импульсного акустического излучения, сопровождающего взрывообразный рост ядер в перегретой воде. Экспериментальная термоакустическая установка (Рис. 2) состояла из термостата 1 и вспомогательного сосуда-охладителя 2, заполненных синтетическим маслом 8АЕ 25'^50 (АР1 СО - 4/81), предназначенным для использования в высокооборотных турбодизелях.
Рис. 2. Структурная схема экспериментальной установки
Масло имело температуру вспенивания порядка 395 0С и среднюю
з
плотность <рр> @ 853 кг/м . В термостате располагался плоский электронагреватель, питавшийся от регулируемого источника постоянного тока 3, обеспечивавшего стабилизацию заданной температуры масла в термостате. Акустический пьезокерамический широкополосный гидрофон (пьезокерамика Р7Т-10) в виде цилиндрической пластины толщиной 1 мм и диаметром 510 - 2 м размещался у дна вспомогательного сосуда с целью предотвращения перегрева выше точки Кюри.
Перед началом измерений в механическом смесителе, работающем в режиме пузырьковой гидродинамической кавитации, приготавливалась водотопливная эмульсия заданной обводнённости и дисперсности. Порционное устройство 9 генерировало каплю эмульсии в термостат, в котором температура ступенчато с интервалом АТ @ 10 К могла изменяться от Тт,п = 387 К до Ттах = 570 К. Нагревание эмульсии сопровождалось вскипанием микрообъёмов воды, содержащих ядра конкурентной фазы с размерами равными или превышающими критический для данной температуры.
Взрывообразное изменение фазового состояния воды сопровождалось генерированием акустического импульса (Рис. 3) вида
РА М - РА(тах) ехР(- ка(ю! + у) , (6)
где Т - текущее время, кА - постоянная затухания, ю - циклическая частота, у - начальная фаза. Следует отметить, что экспоненциально затухающие колебания с периодом около @ 10 - 4 с обусловлены резонансным возбуждением термостата и вспомогательного сосуда. Короткий импульс, вызванный микровзрывом микрообъёма воды, возбуждает радиальные колебания цилиндрической
Рис. 3. Акустический импульс
поверхности термостата и вспомогательного сосуда. Акустический импульс, преобразованный пьзоэлементом 8 в электрическое напряжение, усиливался широкополосным усилителем 7 и подавался на фильтр нижних частот 6 с целью подавления низкочастотных структурных шумов и на вход запоминающего осциллографа 10. Преобразователь 4 служил триггерный преобразователь экспоненциально затухающих импульсов с заполнением в прямоугольные импульсы малой длительности и большой крутизны фронтов. Количество импульсов регистрировалось цифровым частотомером 5, работающим в режиме суммирования. Таким образом, аппаратурный комплекс и применяемая методика позволяли измерять число ядер конкурентной фазы в единичном объёме воды, терявших устойчивость при заданной степени перегрева рабочей жидкости.
Эксперименты проводились с водотопливными эмульсиями на основе дизельного топлива, авиационного керосина, флотских мазутов и топочных мазутов. На рис. 4 точками показана зависимость О = ^Т) для мазута марки М-60 со среднеарифметическим диаметром включений воды <dw> = 510 - 5 м. Теоретическая зависимость получена подстановкой решения уравнения (4) в распределение (5) при экспериментальном значении коэффициента А = 310 - 5 м - 1. Пересчёт значений О(Т) в соответствующие значения Я;;(Т) осуществлялся так же по уравнению (4). Сопоставление теоретических и экспериментальных данных позволяет сделать вывод о правомерности разработанной физической модели процессов. Следует отметить, однако, что электронная реализация термоакустического метода затрудняла регистрацию при температурах, превышающих Т @ 560 К, т.к. частота следования импульсов превышала Г > 100 ; - 1, что не намного меньше температуры максимально возможного перегрева воды Ттах = 575 К. При температурах близких к Ттах изменяется природа образования ядер конкурентной фазы. Наряду с гетерогенным процессом происходит гомогенное образование паровых ядер в результате тепловых флуктуаций
плотности воды, которые приводят к возникновению кластерных образований.
4.1
1,7
0,58
0,19 Я - 10!
м
О
ю1
10'
10
10
10*
•••
Л* • 1
А ^—
•Л •••
т
/ й • •
А*
т/ / • •
•
# 0
) /
•
# /
400
440
480
520
Т. К
Рис. 4. Зависимость числа ядер конкурентной фазы в единице объёма воды от температуры Как следует из полученных результатов, для капель эмульсии с dw < 10
- 4
м вероятность присутствия гетерогенных ядер достаточно велико. Для капли эмульсии с средним арифметическим диаметром ^р> = 10 - 4 м при объёмной концентрации воды ф = 0,2 её объём Ур @ 10 - 13 м3. Как следует из данных, приведенных на рис.4, для Т < 573 К трансформация воды в пар по гетерогенной схеме происходит только в каждом пятом включении воды, остальные четыре капли воды взрываются по гетерогенной схеме, т.е.
микровзрыву таких капель предшествует процесс образования ядра конкурентной фазы при температурах, близких к максимально возможным перегревам воды.
Последнее обстоятельство делает необходимы внести некоторые коррективы в ранее опубликованные модели микровзрыва капель эмульсии, разработанные на основе анализа экспериментальных данных полученных
в лабораторных условиях для относительно крупных капель эмульсии, по— з
рядка 10 м [6, 7], которые на наш взгляд, не могут быть безоговорочно перенесены на реальные камеры сгорания, особенно, работающие в условиях повышенных давлений (дизельные двигатели). Схемы прогрева и воспламенения капель эмульсии, рассматриваемые в известной литературе, достаточно правдоподобно описывают процессы в камерах сгорания, работающих при атмосферном давлении (топочные устройства котлов различного назначения), где в распыленном топливе присутствует большое количество капель топлива крупного размера.
Процесс начальных стадий прогрева капель полидисперсных капель эмульсии (Рис.1) целесообразно рассмотреть в предположении его сферической симметрии и охарактеризовать термодинамическими параметрами эмульсии, являющимися функциями температуры, теплопроводностью 1Е, удельной теплоёмкостью сЕ, коэффициентом поверхностного натяжения се, плотностью рЕ и коэффициентом кинематической вязкости УЕ [7]
1Е = 1 р + (1 № +1 р )фу [1 - 0,72(1 - фу)], (7)
СЕ _ фУ + (1 - фУ )ср , (8)
п е =Пр (1 + 2,36фт), (9)
РЕ _ рр - фУ (рШ - рр ) , (10)
Се =ар [1 + (2,4 - 3-10-21 )фт ], (11)
где фу = УШ/(УШ + Ур), фт = тШ/(тШ+тр) - объёмная и массовая концентрация воды в топливе, 1р - коэффициенты теплопроводности воды и
топлива, сШ и ср - соответствующие удельные теплоемкости, ср - ко-
эффициенты поверхностного натяжения, 1 - температура в градусах Цельсия. Теплообмен капли эмульсии с внешней высокотемпературной средой в такой постановке описывается уравнением
ce Pi
ЭТ _ Э_
Эх _ г2 Эг
1ЕГ2 ^
Е Эг
(12)
где х - время, г - текущее значение радиуса.
При постоянстве внешней температуры Ts = const и Т(г,0) _ Т0 на поверхности капли выполняются граничные условия третьего рода
ЭгТ(0, х)_ 0,
ЭТт^,х)_ H[Ts -t(rf,х)
v(0. х)_ 0,Т(0, х)*~,
Эг
(13)
где Н - постоянная величина, определяющая величину теплового потока на границе со средой, имеющей температуру Т8.Общее приближённое решение уравнения (12) с граничными условиями (13) имеет вид
0(г, х) Ts - Т(г, х)_£
2Jn
RF sin
Jn
Rf
_ Jn - sin Jn cos Jn
xexp(-J • Fo), (14)
n _1
где Фп - корни характеристического уравнения, ро - критерий Фурье.
На рис.5 приведены результаты численного решения уравнения (12) полученные для внешней температуры Т8 = 103 К. Температура на поверхности капли эмульсии Тр зависит от объёмного содержания в ней воды. Так, например, для капли эмульсии радиусом Яр = 210 - 4 м температура её поверхности, соответствующей началу процесса образования конкурентной паровой фазы, при объёмной концентрации воды фУ = 0,05 составляет Т1 = 587 К, а при фУ = 0,2 - Т2 = 519 К, при фУ = 0,5 - Т3 @ 586 К.
Фазовая трансформация вода-пар начинается, в соответствие с развиваемой моделью, при условии
г
г
■‘•'■уах
фУ | |N(11)сШУ =
где Ур = 4/з(яЯр) - объём капли эмульсии. Уравнение (15) предполагает, что математическое ожидание числа ядер, способных терять устойчивость равно единице.
580
560
540
520
500
480
460
440
420
400
\1
4— \ < Т - 10’ К
1- 2- Фу= 0,05, Ф„= 0,2,
V 3- Фу= 0,;
\з
[23^ Я/1 О4, м
Рис. 5. Зависимость температуры поверхности капли от её диаметра при температуре внешней среды Т = 103 К
Как видно из уравнения (14) в объёме капли вдоль её радиуса будет существовать градиент температуры с максимумом на поверхности. Присутствующие в приповерхностном слое капли эмульсии частички воды будут претерпевать фазовые превращения при условии выполнения условия (15). Будет наблюдаться поверхностное вскипание исходной капли, сопровождающееся эффектом разбрызгивания. Из объёма капли вследствие множественных микровзрывов в приповерхностном слое будут вырывать-
ся мелкие сферические фрагменты - топливный туман. На практике этот эффект наблюдался в камерах сгорания котлов в виде увеличения яркости факела при использовании водотопливных эмульсий [6, 7]. Фотографии этого эффекта в настоящее время используются в рекламных целях практически всеми производителями устройств для производства водотопливных эмульсий.
Как следует из данных рис. 5 в каплях эмульсии малого диаметра испарение воды протекает при более высоких температурах, что предполагает качественно иное протекание процесса. Такие капли будут воспламеняться до момента возникновения в их объёме конкурентной фазы. Процессы трансформации воды в пар в этом случае происходят по гомогенной схеме на конечных стадиях выгорания капель.
При использовании водотопливных эмульсий в камерах сгорания котлов, работающих на тяжёлых сортах топлива, ввиду присутствия в факеле относительно крупных капель, желательно иметь дисперсность, обеспечивающую при заданной обводненности минимальное количество включений воды. На рис.6 приведена фотография эмульсии, в которой содержится два сферических включения воды существенно разного объёма: Яр @ 610 - 4 м, @ 2,510 -4 м,
@ 2,5-10 - 5 м. Фазовый переход воды в пар начнётся при температуре не превышающей значение Тсг = 440 К. Для таких капель будет характерен микровзрыв, следствием которого станет дробление исходной капли на не-
Рис. 6. Капля крупнодисперсной эмульсии
сколько более мелких фрагментов. Процесс будет протекать по гетерогенной схеме, т.к. в этом случае К(Я) >> 1.
Полученные данные позволяют заключить, что капли водотопливной эмульсии сопровождаются при их нагревании в камерах сгорания явлением микровзрыва, параметры которого определяются объёмным содержанием воды в топливе и дисперсности эмульсии. Возможны три качественно разных процесса фазовых трансформаций воды, каждый из которых сопровождается отличными друг от друга схемами дробления исходных капель топлива.
Литература
1. Лебедев О.Н. Некоторые особенности горения капель водотопливной эмульсии в дизелях // Физика горения и взрыва. -1978. -№2. - С.142-145.
2. Кнэпп Р., Дэйли Дж., Хэммит Ф. Кавитация: Пер. с англ. - М.: Мир, 1974. - 687 с.
3. Исаков А.Я. Некоторые особенности микровзрыва капли водотопливной эмульсии // Физика горения и взрыва..- 1986. -№1 -С.125 - 126.
4. Исаков А.Я., Дёминов В.И. Физическая модель процессов, предшествующих воспламенению капель водотопливной эмульсии // Физика горения и взрыва. 1986. - №6 - С. 15 - 20.
5. Исаков А.Я. Исследование условий образования конкурентной фазы в диспергированных жидкостях термоакустическим методом // Труды IV Дальневост. акуст. конф. «Акустические методы и средства исследования океана». -Владивосток: ДВНЦ АН СССР, Минвуз РСФСР, ДВПИ, 1986. -С.143 - 147.
6. Иванов В.М. Топливные эмульсии. - М.: АН СССР, 1962. - 407 с.
7. Лебедев О.Н., Сомов В.А., Сисин В.Д. Водотопливные эмульсии в судовых дизелях. Л.: Судостроение, 1988. - 105 с.