Введение
В природе и технике известно много процессов, сопровождающихся фазовыми переходами первого рода – жидкость-пар. Современная техника, особенно транспортная и энергетическая отрасль, построена на преобразовании потенциальной энергии топлива в кинетическую энергию, где основным технологическим процессом является переход рабочего тела из жидкого в газообразное состояние. Возрастание тепловой мощности и интенсификация процессов в тепловых аппаратах приводит к тому, что фазовому переходу жидкость-пар предшествует метастабильное (перегретое) состояние жидкости. Вскипание перегретой жидкости сопровождается значительным выделением запасенной в ней энергии. Этот факт приходится учитывать при анализе аварийных ситуаций, связанных с разгерметизацией сосудов и аппаратов заполненных жидкостью при высоких давлениях и температурах. Поэтому активно исследуются проблемы под общим названием «паровой взрыв» [1]. Анализ процессов интенсивного фазового перехода из перегретого жидкого состояния в газообразное требует детальной информации об условиях зарождения паровой фазы в перегретой жидкости и ее дальнейшего роста.
Детальное рассмотрение кинетики процессов испарения и конденсации при зародешеобразовании (теория Фольмера – Деринга – Зельдовича – Кагана) [2, 3, 7–9] показывает, что частота спонтанного зародышеобразования (число зародышей, образующихся в единице объёма в единицу времени) определяется выражением:
(1)
где k – постоянная Больцмана; N0 – число молекул в единице объёма, шт; В – коэффициент, близкий к 1010 1/с. Произведение N0В слабо изменяется с температурой. Частота флуктуационного зародышеобразования в зависимости от температуры устанавливается практически мгновенно (10-9 с), то есть при реализации температурных условий в жидкости нестационарностью можно пренебречь.
Методика экспериментов
Максимальный перегрев жидкостей ограничен температурой гомогенной нуклеации, т.е. спонтанным образованием паровых зародышей, размеры которых превосходят критический. Па таких зародышах жидкость способна чрезвычайно быстро вскипать. Частота спонтанного образования зародышей определяется формулой (1). Зависимость частоты от температуры выражена очень сильно. В случае перегретого этилового эфира, например, при температуре 143 °С частота нуклеации равна1. Выше этой температуры всего на 2 градуса вероятность появления зародыша увеличивается в 105 раз, а на 3 градуса 1010 раз. Это показывает, что, с одной стороны, жидкость чрезвычайна устойчива по отношению к появлению спонтанных зародышей, а, с другой стороны, в узком температурном диапазоне, который можно характеризовать определённой температурой, называемой температурой гомогенной нуклеации, с вероятностью, практически равной единице, она вскипает.
Практическая реализация вскипания на спонтанных центрах парообразования может быть достигнута посредством быстрого нестационарного нагрева. При стационарном нагреве осуществление подобной задачи связано с рядом сложностей из-за наличия в жидкости искусственных центров парообразования, расположенных на поверхностях дисперсных частичек, загрязняющих жидкость.
Осуществление режима вскипания на спонтанных центрах парообразования может быть достигнуто посредством быстрого нестационарного нагрева. Кратковременный нагрев погружённых в жидкость платиновых проволочек был использован Павловым и Скриповым для изучения гомогенной нуклеации жидкостей [4]. Идея эксперимента состоит в быстром нагреве граничащего с проволочкой слоя жидкости. При этом тепловыделение значительно превышает интенсивность стоков тепла на искусственных центрах парообразования. Роль этих центров, как регуляторов температуры, становится в этом случае малоэффективной и жидкость можно перегреть до температуры гомогенной нуклеации. Кипение, развившееся на флуктуационных центрах, нарушит при этом режим нагрева проволоки. А так как развитие этого процесса происходит в узком температурном диапазоне и носит ярко выраженный «сосредоточенный» характер, то это проявится в виде резкой, различимой «особенности» в температуре проволочки. Такой режим нагрева жидкости авторы [4] называют ударным. Теория метода и, в частности, критерии реализации ударного нагрева подробно рассматривались в [5]. Отметим лишь, что достижение температур гомогенной нуклеации возможно при значительных скоростях нагрева жидкости. Для органической жидкости скорости нагрева превосходят 2·105 К/с при нагреве от комнатной температуры. При меньших скоростях возможно вскипание жидкости на искусственных центрах парообразования.
Рассмотрим, следуя [6], эксперименты по определению температур гомогенной нуклеации жидкостей. Для регистрации начала нуклеации использовался резкий локальный подъём давления при распаде метастабильного слоя жидкости, граничащего с проволокой. Этот слой является источником достаточно сильного акустического излучения, который может быть воспринят приемником.
Динамика распада тонкого слоя жидкости, граничащего с проволочкой, представляется следующим образом. При достижении гомогенной нуклеации Т' в слое жидкости на случайных местах возникают флуктуационные зародыши, которые быстро вырастают до видимых размеров, сливаются между собой и образуют вокруг проволочки паровую пленку. Для количественных оценок, характеризующих образование паровой пленки на поверхности проволочки, воспользуемся [5], где описывается термодинамический кризис кипения. При импульсном кратковременном нагреве источника, погружённого в жидкость, реализуется именно этот режим кипения.
Оценка времени формирования паровой пленки для рассмотренного частного случая может быть применена для случая кратковременного нагрева проволочки импульсом тока. Физической предпосылкой использования названной оценки является то, что время образования паровой пленки после достижения температуры нуклеации существенно меньше времени нагрева поверхности от исходной температуры. Поэтому процесс распада можно отнести к постоянной температуре.
Соответствующие расчеты показывают, что время распада гомогенной фазы составляет величину порядка 10-6с. Данные оценки согласуются с экспериментами. В частности, по результатам экспериментального исследования [4], в котором проводилось фотографирование отдельных стадий процесса, от момента появления первых, видимых под микроскопом, пузырьков газа в гексане (начало кипения) до формирования парового чулка (наступления термодинамического кризиса кипения) проходит 2–3 мкс.
Регистрация времени нуклеации жидкости определяет температуру проволочки, которая отождествляется с температурой нуклеации жидкости. Так как момент нуклеации распознается по акустическому сигналу, то следует учитывать некоторое смещение значений температуры во времени. В частности, начало нуклеации относится к появлению первых зародышей флуктуационной природы, начало же интенсивного акустического излучения, регистрируемого в экспериментах, относится к появлению паровой пленки. Именно в этот момент коллективный эффект единичных паровых пузырьков, растущих на поверхности нагревателя и являющихся источниками акустического излучения, достигает наибольшей сосредоточенности.
На рисунке 1 показана принципиальная схема экспериментальной установки, применённой в исследованиях. Установка содержит: генератор прямоугольных импульсов напряжения (Г), блок задержки (БЗ), мостовую схему, в которую включается измерительный элемент. Измерительный элемент представляет собой сосуд, заполняемый исследуемой жидкостью, в которой находится проволочка, а в непосредственной близости от нее – поляризованная керамика на основе титаната бария.
Рисунок 1. Схема экспериментальной установки и характерный сигнал начала нуклеации
Рассмотрим процедуру измерений. На проволочку, находящуюся в жидкости, подаются прямоугольные импульсы напряжения малой длительности. Момент подачи импульса на проволочку синхронизирован с запуском осциллографа, на входы которого подключены выходы с датчика давления и диагонали моста. Плавно увеличивают напряжение в импульсе до тех пор, пока с датчика давления не будет зарегистрирован резко выраженный сигнал, соответствующий началу нуклеации жидкости (см. рисунок 1). Затем посредством блока задержки вводится время задержки и, соответственно, уменьшается время развёртки с тем расчётом, чтобы более подробно рассмотреть начало сигнала, обусловленного приходом акустической волны. Это позволит более точно определить момент прихода сигнала τ.
Одновременно, изменением сопротивления R напряжения в импульсе добиваются, чтобы моменту прихода сигнала τ соответствовал мгновенный баланс мостовой схемы. Из условия мгновенного баланса моста R1/Ri = R2/R3 определяется сопротивление проволочки Ri(τ)и, следовательно, ее температура. Эту температуру при определенных ограничениях на длительность импульса можно отождествить с температурой нуклеации Т'. Действительно, время прихода акустического сигнала τ несколько превышает время нагрева проволоки τ/ до температуры нуклеации. Время превышения складывается из рассмотренного выше времени формирования акустического сигнала τ0 и времени запаздывания τ3, определяемого длительностью распространения звуковой волны от проволочки до датчика
давления. Нагрев проволоки за суммарное время τ0+τ3 определит погрешность T/. Можно ввести коррекцию на время запаздывания сигнала: она тем точнее, чем сильнее неравенство τ3«τ. Коррекция же на время τ0 в значительной степени остается неопределенной. Вследствие этого время нагрева жидкости до температуры её нуклеации должно быть ограничено снизу неравенством τ»τ0. Оценочное значение τ0 определяется по формуле (12).
В экспериментах использовались вольфрамовые проволочки диаметром 20 мкм и длиной 20 мм. Длительность нагрева до температуры нуклеации составляла 400 мкс. После достижения температуры нуклеации импульс обрывался. Собственная частота колебаний пластинки титаната бария составляла 1 мГц.
Результаты исследования
В таблице1 представлены экспериментальные исследования температуры гомогенной нуклеации Т/ ряда органических диэлектрических жидкостей при атмосферном давлении. Литературные данные взяты из [5].
Таблитца 1. Температуры гомогенной нуклеации органических жидкостей
№ |
Жидклость |
Хим. формула |
Т экспер |
Т литер |
1 |
Толуол |
С7Н8 |
540,5 |
|
2 |
Н-октан |
С8Н18 |
506,8 |
|
3 |
Н-додекан |
С12Н26 |
593 |
|
4 |
Н-тридекан |
С13Н28 |
598,6 |
|
5 |
Н-тетрадекан |
С14Н30 |
618,5 |
|
6 |
Н-пентадекан |
С15Н32 |
639 |
|
7 |
Н-гексадекан |
С16Н34 |
660 |
|
8 |
Н-октадекан |
С18Н38 |
682,5 |
|
9 |
Н-тетрокозан |
С24Н50 |
747,6 |
|
10 |
1-гексен |
С6Н12 |
453,6 |
|
11 |
1-октен |
С8Н16 |
510 |
|
12 |
1-нонен |
С9Н18 |
527 |
|
13 |
Бензол |
С6Н6 |
511,4 |
|
14 |
Н-ксилол |
С8Н10 |
556 |
|
15 |
Этилбензол |
С8Н10 |
536,4 |
|
16 |
Дифенил (ВОТ) |
С12Н10 |
726 |
|
17 |
Дифениловый эфир |
С12Н10О |
681,7 |
|
18 |
Триэтиламин |
С6Н15N |
481,5 |
|
19 |
Маслянная кислота |
С4Н8O2 |
571,5 |
|
20 |
Дециловый спирт |
С10Н22O |
637 |
|
21 |
Н-гексан |
C6H14 |
457 |
406 |
22 |
Метиловый спирт |
СН4О |
456,5 |
396 |
Результирующая погрешность измерений температуры гомогенной нуклеации оценена нами в 3 градуса.
В таблице 1 приведены рассчитанные по теории Фольмера – Деринга, Зельдовича –Френкеля температуры гомогенной нуклеации исследуемых жидкостей. Рассмотрим методику расчёта. Логарифмируя выражение (1), имеем:
(2)
Или учитывая, что для жидкостей =71,
(3)
При выборе величины J1 воспользуемся следующим соображением. В реальной системе объёма V, существующей время τ, частота зародышеобразования равна J=J1Vτ. Положив J=1 (распад жидкости на одном зародыше уже даёт принципиальную возможность регистрации процесса), имеем J1=(Vτ)-1. При расчёте предполагалось, что время τ равно длительности нагрева жидкости, а объём жидкости, охваченный измерением, равен V=lπ(lЭФ)2, где l – длина проволоки, lЭФ=(аτ)0,5 , а – температуропроводность жидкости. С учётом сильной зависимости частоты зародышеобразования жидкостей от температуры (3–4 порядка на градус) такие приближения оценки величины J1 представляются приемлемыми. Для приведённых выше условий эксперимента величина J1 равна 1015 см-3с-1.
Выводы
Как видно из таблицы 1, экспериментальные и расчётные значения температур гомогенной нуклеации находятся в хорошем соотношении между собой. В таблице, по сравнениям, приведены экспериментальные данные по Т/, полученные методом импульсного нагрева с регистрацией теплового сигнала с проволочного зонда [5].
В рамках настоящей работы проведено исследование акустического излучения, сопровождающего фазовый распад при значительных скоростях нагрева проволочки (до 106 град/с). Характерное время, определяющее прохождение единичного температурного интервала, достигало при этом значений порядка 10-8 с, то есть становилось соизмеримым с временем релаксации при выделении паровой фазы.
Рецензенты:
Спирин Геннадий Георгиевич, д.т.н., профессор кафедры «Физика» ФГБОУ ВПО «Московский авиационный институт (национальный исследовательский университет)», г. Москва.
Камзолов Сергей Константинович, д.т.н., профессор кафедры «Физика» Московского государственного технического университета гражданской авиации (МГТУ ГА), г. Москва.