Закономерности теплоотдачи и кризиса кипения в воде, недогретой до температуры насыщения

Зар Ни Аунг

Закономерности теплоотдачи и кризиса кипения в воде, недогретой до температуры насыщения тема автореферата и диссертации по физике, 01.04.14 ВАК РФ

Зар Ни Аунг АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Москва МЕСТО ЗАЩИТЫ

2013 ГОД ЗАЩИТЫ

01.04.14 КОД ВАК РФ

Диссертация по физике на тему «Закономерности теплоотдачи и кризиса кипения в воде, недогретой до температуры насыщения»

Автореферат диссертации на тему "Закономерности теплоотдачи и кризиса кипения в воде, недогретой до температуры насыщения"

На правах рукописи

Зар Ни Аунг

ЗАКОНОМЕРНОСТИ ТЕПЛООТДАЧИ И КРИЗИСА КИПЕНИЯ В ВОДЕ, НЕДОГРЕТОЙ ДО ТЕМПЕРАТУРЫ НАСЫЩЕНИЯ

01.04.14 - «Теплофизика и теоретическая теплотехника»

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

Автор: (

21 ноя т

Москва 2013

005538575

Работа выполнена в Национальном исследовательском ядерном университете

«МИФИ»

Научный руководитель: кандидат физико-математических наук, доцент

Лаврухин Алексей Анатольевич, НИЯУ МИФИ Официальные оппоненты: доктор технических наук, профессор

Ягов Виктор Владимирович, профессор кафедры инженерной теплофизики НИУ «МЭИ»

кандидат физико-математических наук Мосунова Настасья Александровна, заведующая отделением теоретических исследований ИБРАЭ РАН

Ведущая организация: Федеральное государственное унитарное

предприятие «Государственный научный центр Российской Федерации Физико-энергетический институт имени А.И. Лейпунского»

Защита состоится 11 декабря 2013 года в 15 час. 00 мин. на заседании диссертационного совета Д 212.130.04 в НИЯУ МИФИ по адресу: 115409, Москва, Каширское шоссе, д. 31.

С диссертацией можно ознакомиться в библиотеке НИЯУ МИФИ. Автореферат разослан 7 ноября 2013 года.

Просим принять участие в работе совета или прислать отзыв в одном экземпляре, заверенном печатью организации, по адресу НИЯУ МИФИ.

Ученый секретарь диссертационного совета,

д. ф.-м.н., профессор / И.И.Чернов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

Научная новизна диссертационной работы.

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

• Разработано новое корреляционное уравнение для расчета критического теплового потока при кипении воды с недогревом на горизонтальных цилиндрических нагревателях в условиях большого объема при атмосферном давлении.

• Впервые проведен анализ температурных флуктуаций теплоотдающей поверхности. Установлено, что по характеристикам распределения отклонений перегревов нагревателя от среднего значения можно судить об устойчивости режимов теплообмена и, в частности, диагностировать начало и режимы кипения.

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

• Представлено обобщение опытных данных, описывающее влияние набрасываемой тепловой нагрузки и недогрева воды до температуры насыщения на характеристики нестационарного кризиса кипения.

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

Достоверность полученных экспериментальных данных подтверждается: отработкой методики измерений в тестовых экспериментах, воспроизводимостью опытных данных. Эмпирические уравнения и обобщающие зависимости согласуются с имеющимися в литературе данными для соответствующих диапазонов режимных параметров. Обоснованность предложенных методик подтверждается сравнением с известными представлениями о механизмах теплообмена при кипении и экспериментальными данными.

Основные положения, выносимые на защиту.

• Результаты экспериментальных исследований характеристик теплообмена в насыщенной и недогретой до температуры насыщения воде при стационарном и импульсном подводе тепла.

• Эмпирическое уравнение, описывающее влияние недогрева воды до температуры насыщения и диаметра горизонтального цилиндрического нагревателя на критический тепловой поток.

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

• Основанную на проверенных физических представлениях и экспериментальных фактах методику расчета коэффициента теплоотдачи в переходных режимах кипения жидкости.

• Новые опытные данные и их обобщение, описывающие влияние недогрева воды до температуры насыщения на характеристики кризиса кипения при быстром возрастании мощности тепловыделения.

Личный вклад автора. Постановка задачи исследования осуществлена совместно с научным руководителем. Экспериментальные исследования и анализ полученных результатов, разработка моделей и методик выполнены автором самостоятельно или в соавторстве при его непосредственном участии.

Апробация работы. Материалы диссертационной работы докладывались на следующих научных сессиях и семинаре: Научные сессии НИЯУ МИФИ 2011-2013; XVII школа-семинар по проблемам физики реакторов «Волга-2012».

Публикации. Основные результаты диссертации опубликованы в 8 работах, в том числе в 3-х статьях в реферируемых журналах из перечня ВАК.

Структура и объем диссертации. Диссертация состоит из введения, пяти глав, заключения и списка литературы из 119 наименований. Содержание диссертации изложено на 141 странице, включая 55 рисунков и 1 таблицу.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе приведен обзор российских и зарубежных работ, посвященных закономерностям кипения воды, недогретой до температуры насыщения, а также особенностям протекания процессов теплообмена при нестационарных тепловых воздействиях.

На основании анализа литературы сформулирована цель работы и поставлены следующие задачи исследования.

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

• На основе новых опытных данных описать влияние недогрева и диаметра нагревателя на первый критический тепловой поток, включая сравнение с имеющимися в литературе экспериментальными исследованиями.

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

• Разработать новую методику расчета теплоотдачи в переходном режиме от однофазной конвекции к развитому пузырьковому кипению.

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

Во второй главе описаны экспериментальная установка, тепловыделяющие рабочие участки, методика проведения экспериментов и приведен анализ погрешностей измерений.

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

В качестве рабочих участков использовались платиновая проволока марки "Экстра" диаметром <1 = ОД мм и нихромовые проволоки диаметрами 0,18; 0,4 и 0,5 мм. В случае платиновой проволоки она одновременно служила и термометром сопротивления. Измерения выполнялись с рабочим участком, расположенным горизонтально в большом объеме воды при атмосферном давлении.

Опыты со стационарным подводом тепла проводились в насыщенной воде и воде, недогретой до температуры насыщения. Тепловыделение в нагревателе задавалось источником питания постоянного тока, падение напряжения через проволоку и на образцовом сопротивлении измерялось цифровыми вольтметрами. Суммарная погрешность при определении плотности теплового потока q составила не более 5 %. Разность между средней температурой платиновой проволоки и температурой воды ЛТ= Т—Тж измерялась с точностью не менее 2 К, а расчетная погрешность в определении коэффициента теплоотдачи а = д/ДГ составляла не более 15 %.

В опытах с нестационарным тепловыделением импульс мощности в нагревателе формировался тиристорным ключом, который замыкал и размыкал цепь питания рабочего участка. Открытие и закрытие ключа производилось при помощи ТТЬ сигналов с автоматизированной системы управления. Этими же сигналами осуществлялась синхронизация системы измерений. Погрешность при измерении интервалов времени составляла менее 10"7 с. Суммарная погрешность при определении набрасываемого теплового потока <7„(т) составляла не более 5 %, перегрева АТ(х) - 10 %, отводимого в жидкость потока тепла q(z) -15%.

В третьей главе представлены результаты экспериментов, проведенных с насыщенной и недогретой до температуры насыщения водой в условиях большого объема при стационарном тепловыделении в нагревателе.

На рис. 1 представлены кривые кипения воды, полученные на горизонтально ориентированной платиновой проволоке в условиях большого объема при атмосферном давлении. Недогрев 0 = Г5 - Тж в опытах менялся в пределах от 0 до 80 К. По мере увеличения теплового потока рабочий участок охлаждался в режиме естественной конвекции (пунктирные линии на рис. 1).

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

д, Вт/м2_ д, Вт/м2_

107т

• I I I М| I I I I II М| I I I I I I 1|| • | I I I I I I I | 1

10 102 А Т, К ю А т, К

Рис. 1. Кривая кипения воды под атмосферным давлением: О-данные эксперимента; __-естественная конвекция, расчет по N11 = 1,1811а°'125; - пузырьковое кипение, расчет по формуле (1); +-усредненные данные по гистограммам (рис.5); а - насыщенная вода; б - недогрев 0 = 40 К.

Пузырьковое кипение в жидкости, недогретой до температуры насыщения, принято делить на два режима: неразвитое и развитое кипение. В области неразвитого кипения процесс теплопередачи определяется как чисто конвективным переносом теплоты, так и дополнительным переносом тепла в результате микроконвекции и фазового перехода в примыкающем к стенке перегретом слое жидкости. При увеличении плотности теплового потока число активных центров парообразования растет. То есть, постепенно подавляется влияние на теплообмен чисто конвективных факторов и увеличивается роль факторов, связанных с кипением. Начиная с некоторой граничной величины перегрева теплоотдающей поверхности относительно температуры насыщения АГ5 = Т-Г, режим кипения становится развитым. В этой области коэффициент теплоотдачи, рассчитанный относительно температуры насыщения, а, = слабо зависит от недогрева 6 и может быть рассчитан как

а 5=Д-<Л

(1)

На рис. 2 представлены данные о зависимости коэффициента теплоотдачи в режиме развитого кипения а5 от теплового потока д для различных значений недогревов 0. Можно отметить только слабую тенденцию увеличения коэффициента теплоотдачи с ростом недогрева. Для всех исследованных значений не-догрева были получены значения т = 0,8 при погрешности в определении коэффициента А до 10 %. Полученные экспериментальные данные хорошо согласуются с известными данными других авторов.

а5, Вт/(м2 К)

Недогрев воды <<

до температуры

насыщения К

о 0 + 50

105- □ 10 * 60 Л

* о 20 73>#Г

V 30

д 40

104т - а)

10»

Т" 106

А( 6)

с/, Вт/м2

Рис. 2. Теплоотдача при развитом пузырьковом кипении воды с недогревом на нагревателе из платины, с1 = 0,1 мм ( - расчет по формуле (1) при т = 0,8, А = 0,63).

Результаты наших экспериментов, а также опытные данные других авторов о зависимости критического теплового потока от диаметра горизонтального цилиндрического нагревателя при кипении насыщенной воды в диапазоне ШЬ = 0,04 - 0,40 могут быть описаны уравнением

где Ь = ^/стТ^Т—р*] - постоянная Лапласа, р', р" - плотности жидкости и пара на линии насыщения; а - коэффициент поверхностного натяжения, g-ускорение свободного падения, ^Кр1*,о ~ 1,2-10б Вт/м2 - критическии тепловой поток при недогреве 8 = О К для плоской горизонтальной поверхности больших размеров, обращенной греющей стороной вверх, который можно рассчитать как

?кР1.,о = 0,14г д/р' ^ст£(р'-р"),

(3)

где г - скрытая теплота парообразования.

На рис. 3 в безразмерных координатах показано влияние недогрева воды на критический тепловой поток. Видно, что полученные нами данные по е в исследованном диапазоне ШЬ могут быть описаны линейной зависимостью от недогрева (С.С. Кутателадзе, 1979):

9кр1,в — <7кр1«,0

1 +А

3/4 п\

ср 9

(4)

где ср - изобарная теплоемкость жидкости на линии насыщения, А ~ 0,1 для нагревателей с большим относительным диаметром, <НЬ > 1.

Рис. 3. Влияние недогрева воды и диаметра нагревателя на критический тепловой поток: <>-¿ = 0,1 мм; □-¿ = 0,18 мм; Д - <1 - 0,4 мм; V - <1 = 0,5 мм;

- линейные аппроксимации.

Анализ результатов экспериментов с нагревателями диаметром <¿ = 0,11,0 мм позволяет предложить следующее эмпирическое соотношение для расчета коэффициента А в формуле (4):

А = 0,03-

(5)

Окончательно формула для расчета критического теплового потока при кипении воды на горизонтально ориентированных цилиндрических нагревателях в условиях большого объема при атмосферном давлении имеет вид:

9кр1,е

= 1,91 —

9кр1»,0

1 + 0,036/

И 3/4 СР°

и V

(6)

Формула (6) может использоваться в диапазонах 0,1 < ¡1 < 1,0 мм; 0 < 0 < 80 К. Сравнение результатов расчета по формуле (6) с данными наших измерений и опытными данными других авторов приведено на рис. 4.

Ч кр! , эксп , расч

Рис. 4. Сравнение результатов расчета по формуле (б) с опытными данными. Данные эксперимента: О - ^ = 0,1 мм; □ -(/ = 0,18 мм; Д-£/=0,4мм;

V - (I = 0,5 мм; данные Бобрович Г.И. и др., 1964 О-¿=0,2 мм; 3-</ = 0,5мм; данные 1пое Т. е(. а!., 1988 С - ¿= 1,0 мм; — - расчет по формуле (6).

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

Для определения статистических характеристик колебаний температуры нагревателя на различных участках кривых кипения (помечены крестиками на рис. 1) были построены гистограммы вероятности регистрации перегрева в диапазоне АТ ± 0,05 К от мгновенного перегрева рабочего участка (рис. 5). Выборка данных по одному измерению тепловой мощности >уге„ и перегрева ДТ включала 105 точек, время измерения для статистического анализа составляло 100 с.

0,24

0,16

0,08-

<ДГ) = 30,4К а = 0,3 К

0,016

28 Ш

(АТ) = 71,5 К ст = 4Д К

0,008-

0,005-

(АТ ) = 69,3 К

1 ст = 4,2 К

j 1 б)

0,04-

0,02-

(АТ ) = 73^4~К с = 1,1 К

АТ, К

75 АТ, К

г)

70 АТ, К

Рис. 5. Гистограммы отклонений перегревов рабочего участка при значениях тепловой мощности: а - 290, б - 1290, в - 1280, г - 6040 кВт/м2; — - нормальное распределение, (АТ) и а — среднее значение и дисперсия мгновенных перегревов в выборке соответственно. Недогрев воды б = 40 К.

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

Полученные в опытах данные о дисперсии в зависимости от среднего значения подводимого к жидкости теплового потока представлены на рис. 6. Сравнивая эти данные с кривыми кипения (см. рис. 1) можно сделать вывод о том, что чем больше о, тем менее устойчив режим. В точках, где дисперсия максимальна, происходит смена режимов теплообмена.

а, К

105 Ч, Вт/м2 10« </, Вт/м2

Рис. 6. Значения дисперсий ст в зависимости от среднего теплового потока д для недогретой до температуры насыщения воды: а - 6 = 0, б - 40 К.

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

Для анализа температурных флуктуаций рассматривалось уравнение теплового баланса в системе «нагреватель — окружающая жидкость» для мгновенных значений входящих в это уравнение величин:

СР1

¿АТ

~сПГ

:Чен(дг)-4(дг),

где q{AT) - отводимый в жидкость тепловой поток.

К правой части уравнения (7) добавлялся случайный шум с малой амплитудой колебаний температуры В(АТа) около равновесного состояния с АТ = ЛГ0:

М/(т,ДГ0) = ^

ср-| -В(АТ0)-<р(т).

(8)

где -1 < ф(т) < 1 — стохастическая функция. В точке равновесия ,Р(ДГо) = и'генСДГо) - д(Д7"0) = 0, и, при разложении функции Т7 в ряд Тейлора в окрестности АТ0 с сохранением только линейных слагаемых, получено безразмерное стохастическое дифференциальное уравнение в вариациях 5Г*:

ёт 4-АТ '

(9)

где т* = ^т, 5Т* = БТ/АТ* - безразмерные время и температура, 1 <1^

4 =

срП С1ДГ

1Д7" = 1К.

Для эволюции плотности вероятности со(5Г, х) использовано уравнение Фоккера-Планка:

Э (а Э(ш-5Г*)_ Э . .^ - к — £*_.

Э т* Э 5 Г* Э (5 ТУ

(10)

где е = -

4-дг*'

Для не зависящей от времени меры вероятности уравнение (10) дает решение в виде распределения Гаусса

сз =

2 л ВАТ

т-ехр

2 ВАТ

(П)

где дисперсия о вычисляется по формуле

[ват7

°=гг

Отношение В1\\\ по физическому смыслу является амплитудой пульсаций температуры и может быть выбрано в качестве критерия, определяющего смену режимов теплообмена. На рис. 7 представлена ее качественная зависимость от перегрева нагревателя.

Рис. 7. Зависимость В!\%\ от перегрева ДТ нагревателя.

Область I - режим естественной конвекции, распределение плотности вероятности температурных флуктуаций носит характер нормального закона.

Область II - закипание, в! возрастает, в точке максимума происходит кардинальная перестройка структуры пограничного слоя.

Область Ш - режим частично развитого кипения с интенсивной активацией центров парообразования.

Область IV - полностью развитый режим кипения.

Количество центров парообразования в областях II и Ш характеризуется значительной нерегулярностью. Добавление температурных флуктуаций приводит к тому, что система переходит от одного состояния со своей плотностью центров парообразования, а значит параметрами дисперсии и среднего значения перегрева, к другому. Гистограммы будут накладываться друг на друга, что приводит к возникновению многомодального распределения.

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

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

Разработанная методика основана на следующих физических предположениях и экспериментальных фактах:

• на теплоотдающей поверхности существуют участки развитого пузырькового кипения с коэффициентом теплоотдачи акип и участки однофазной конвекции с аК0Нв соответственно;

• по мере роста температуры поверхности площадь участков, занятых пузырьковым кипением растет пропорционально количеству центров парообразования;

• количество центров парообразования пропорционально количеству микронеровностей на теплоотдающей поверхности;

• плотность вероятности распределения микронеровностей в зависимости от их глубины подчинено закону нормального распределения Гаусса.

Уравнение теплового баланса для единицы площади поверхности записывается как суперпозиция двух составляющих полного теплового потока:

<7 = «ю„в' (1 -Бт)• АГ + акип ■ • ДГ, (13)

где 5КИГ| — доля площади поверхности, занятая пузырьковым кипением.

Размеры впадин и их количество однозначно связаны с параметрами шероховатости поверхности. Все элементы шероховатости рассматривались как конические впадины с одинаковым углом полураствора (3 и глубиной /¡, изменяющейся хаотически по закону нормального распределения. Наиболее важным геометрическим параметром впадин является радиус устья г0. Если радиус

Го превышает критический радиус кривизны Лкр эф паровой фазы, то впадина будет активной (жидкость кипит на данном центре). Обычно выражение для Лкр эф задается в виде

r -в 2а(.1 + р'/р)

КРЗФ р{ехр[г(Г-Г5)/(/?г Г Г5)]-!}'

(14)

где р - давление, Лг - газовая постоянная для пара, а для воды под атмосферным давлением в среднем 5 = 7. Плотность центров парообразования пропорциональна общему количеству микронеровностей и определяется вероятностью вступления в действие впадин глубиной не меньше величины А™. = Ът^ё Р = Лкр.^ Р тогда

"Г 1

^ V 2л а,

-ехр

2 а

d$=f 1-erf

f^min _j

Ч/гсР у

(15)

2 г 2

где erf (z) = —t= \е'х dx - функция ошибок Гаусса, /гср-средняя высота мик-V;r 0j

рошероховатостей, 0А - дисперсия, Е = hcp/ah - параметр шероховатости поверхности, величина щ = 2пср, где пср - плотность центров парообразования в случае, когда /гт]П = /гср.

Зависимость количества действующих центров парообразования от температуры поверхности была получена заменой величины h в соотношении (15) на ее выражение из формулы (14):

n = ncp-i 1-erf

-^(/(Г,Гср,Г5)-1)

(16)

где f(T,T Ts) =

_ехр[гАГср/(/;гГсрГ5)]-1 ехр[гДГ5 /(/?r7Ts)]-l

, а Гср и ДГср = Гср - Ts - температура те-

плоотдающеи поверхности и ее перегрев относительно температуры насыщения в условиях, когда начинают активироваться впадины глубиной Аср и менее.

Если в среднем площадь поверхности, приходящаяся на один центр парообразования постоянна и равна 1/лср, доля площади поверхности нагрева, занятая пузырьковым кипением определяется следующим выражением:

5кип=1-егГ -Л= (/(Г,Гср,Г5) — 1) .

(17)

После подстановки выражения (17) в соотношение (13) для расчета коэффициента теплоотдачи получена формула

которую можно использовать при 0 < ДГ5 < АГср, а ДГср имеет смысл перегрева относительно температуры насыщения для полностью развитого кипения.

Сравнение предлагаемого расчетного соотношения (18) с экспериментальными данными проводилось при параметре шероховатости е = 5. Коэффициент теплоотдачи при кипении рассчитывался как акип = а5 ДГ5 / АТ. За исключением насыщенной воды все опытные данные были описаны при значении ДГср = 31 ± 3 К. Пример расчета приведен на рис. 8.

а =а,

конв

+ (аиш-авя.)]1-^ -и/(Г,Гср,7;)-1) , (18)

л/2

</, Вт/м2

107-

106-

Рис. 8. Режимы теплоотдачи от проволочного нагревателя к воде под атмосферным давлением, недогретой до температуры насыщения (0 = 40 К): О - данные эксперимента; 1 - естественная конвекция; 2 - расчет по формуле (18) при ДГср = 31 К, е = 5; 3-развитое кипение, расчет по формуле (1).

:е дгср+е

дг,к

Разработанная методика была применена для расчета температуры оболочек твэлов реакторной установки В 491. Исходные данные принимались для отклоненных параметров в горячем подканале: давление в реакторе 15,9 МПа; температура воды на входе 573,4 К; диаметр твэлов 9,1 мм; длина тепловыделяющей части 3,73 м; количество твэлов в активной зоне 50856; проходное сечение теплоносителя в активной зоне 4,14 м2; средний гидравлический диаметр пучка твэлов 10,5 мм; минимальный проектный расход теплоносителя через реактор (на входе) 83100 м3/ч; тепловая мощность на один твэл в горячем подканале 136,6 кВт. Распределение плотности теплового потока по высоте активной зоны приведено на рис. 9.

1,б.ю< 9'Вт/м2

1,2106

Рис. 9. Распределение плотности теплового потока по высоте активной ^ 8-Ю6 зоны.

0,4-10« 0

0 1 2 з г, м

При отсутствии опытных данных по значению Д7"ср для материала оболочки твэлов для расчетов принималось, что эта величина соответствует первому критическому температурному напору, для насыщенной воды при указанных выше параметрах равному примерно 9 К. Коэффициент теплоотдачи аконв составил 42 000 Вт/(м2-К) (Б.С. Петухов и др. 1958 г.), коэффициент теплоотдачи при кипении насыщенной жидкости рассчитывался по методике Н.Г. Рассохина и др. (1970 г.). Изменение температуры теплоносителя и оболочки твэл по высоте активной зоны представлено на рис. 10. Расчеты по предлагаемой методике показали, что можно выделить две характерные области. От начала кипения до интенсивной генерации пара основным механизмом теплоотдачи остается

д, Вт/м2

конвекция. Далее вследствие активации большого количества центров парообразования теплоотдача значительно улучшается. Г, К

640

620------

600-

580

560

Рис. 10. Изменение температуры воды и оболочек твэлов по длине кассеты реакторной установки В491: 1-температура воды; 2 - температура оболочек твэлов согласно методике Петухова Б.С. и др. (1958 г.) и Рассохина Н.Г. и др. (1970 г.); 3-то же по методике, основанной на формуле (18).

z, М

В пятой главе представлены и обобщены новые экспериментальные данные, расширяющие границы применимости методики расчета динамических характеристик процесса теплообмена в условиях быстрого изменения тепловой нагрузки, разработанной ранее в МИФИ (Deev V.l. et al., 2007).

При высоких набрасываемых тепловых нагрузках q„, когда кризис теплоотдачи происходит до момента отрыва пузырей, перегрев нагревателя монотонно увеличивается, а отводимый в жидкость тепловой поток достигает максимума (рис. 11). Параметры нестационарного кризиса теплообмена qyр, АТ^,

AT,К

</, Вт/м2

-25-Ю6

Рис. 11. Зависимость перегрева нагревателя и плотности теплового потока от времени. Данные эксперимента: О - А Г;

□ - 9„; Д -q\__- расчет по уравне-

—10-Ю6 нием нестационарной теплопроводности

до момента закипания;_- расчет по

модели нестационарного кризиса теплообмена Deev V.l. et al., 2007.

"1 I T 0 ОД 0,2 0,3 0,4 Т, мс

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

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

лг.Ф,* = Агпр,5 -(А7"пР,з -АТкР],5)-ехР

-к-

Чкр *?кр1

(19)

АГ„.р,К

Рис. 12. Зависимость критической температуры ДГц, от нестационарного критического теплового потока <7,р. Данные эксперимента в, К: 100 О - 0 (?»Р|.о"= 0,8 10б Вт/м2, ДГ>Р1,, = 20 К), □ -5 (94,1.5= 1,6 106Вт/м2, ДГ«р1,, = 40К), Д-10 (9ч,ио = 2,0 106 Вт/м2, ДГ,Р15 = 40 К), • - данные вЫ 8. й. а1. (2002 г.). Расчет по (19) при к = 30 и р = 0,9 для 0, К:

Вт/м2

В формуле АТщ, 5 и ДГкр1-!! - перегревы воды относительно температуры насыщения при нестационарном и первом критическом тепловом потоке; АТПр? 5 -предельный перегрев гомогенного зародышеобразования, рассчитанный отно-

сительно температуры насыщения; д.

X' г р" ДГп2р, аГ

, X' - коэффициент тепло-

проводности насыщенной воды. Коэффициент к= 30 и показатель степени

р = 0,9 получены аппроксимацией экспериментальных данных, представленных на рис. 12. Вид формулы (19) в дальнейшем, очевидно, должен уточняться по мере накопления новых экспериментальных результатов.

На рис. 11 приведен пример расчета перегрева нагревателя и отводимого в жидкость теплового потока от времени по ранее предложенной методике с учетом новой интерполяционной зависимости. Видно удовлетворительное согласование с экспериментами. Сопоставление результатов расчета критического интервала времени для насыщенной и недогретой воды при высоких тепловых нагрузках с экспериментальными данными приведено на рис. 13.

9кп5 Вт/м2

107-;

106т

100 тьр'мс

Рис. 13. Влияние критического теплового потока <7„р на интервал времени ткр перехода к пленочному кипению воды: для данных эксперимента обозначения те же, что на рис. 12; расчет по модели для 9: - 0 К, ___5 К,— . —-10 К.

0,1 1

Оказалось, что при высоких тепловых нагрузках методика, разработанная для насыщенной воды, применима и для недогретой жидкости. Это, возможно, объясняется тем, что при больших числах Якоба пузыри при своем росте приподнимает обволакивающий его перегретый пограничный слой.

ОСНОВНЫЕ ВЫВОДЫ

1. Предложена формула для расчета критического теплового потока при кипении воды с недогревом на горизонтальных цилиндрических нагревателях в условиях большого объема при атмосферном давлении. Результаты расчета по этой формуле хорошо согласуются с данными проведенных измерений и опытными данными других работ в диапазоне диаметров нагревателей от 0,1 мм до 1,0 мм и недогревах до 80 К.

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

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

3. Предложена методика расчета теплоотдачи в переходных режимах теплообмена от однофазной конвекции к развитому пузырьковому кипению. Расчеты по методике свидетельствуют о том, что можно выделить две характерные области режимов с разными механизмами теплоотдачи. Первая область характеризуется тем, что основным механизмом переноса тепла является однофазная конвекция. Во второй области происходит интенсивная активация центров парообразования, которая сопровождается значительным улучшением теплоотдачи.

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

Основные публикации по теме диссертации

1. В.И. Деев, Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин, Ю.А. Маслов, В.Н. Федосеев, B.C. Харитонов. Методика расчета теплообмена при кипения жидкости на поверхности нагрева // Ядерная физика и инжиниринг, 2011, т. 2, №5, с. 387-394.

2. В.И. Деев, Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин, В.Н. Федосеев. Статистический анализ температурных флуктуаций как метод диагности-

ки режимов теплообмена при кипении // Тепловые процессы в технике, 2013, т. 5, №4, с. 163-169.

3. В.И. Деев, Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин. Динамические характеристики процесса теплообмена в воде при импульсном подводе тепла// Естественные и технические науки, 2013, № 3 (65), с. 35-41.

4. В.И. Деев, К.В. Куценко, A.A. Лаврухин, Зар Ни Аунг. Новый метод расчета теплоотдачи в режиме кипения воды с недогревом // Научная сессия НИЯУ МИФИ-2011. Аннотации докладов, т. 1. Инновационные ядерные технологии. М.: НИЯУ МИФИ, 2010, с. 46.

5. В.И. Деев, Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин. Методика расчета теплоотдачи к недогретой до температуры насыщения воде в предкризисной области И Научная сессия НИЯУ МИФИ-2012. Аннотации докладов. т. 1. Инновационные ядерные технологии. Высокие технологии в медицине. М.: НИЯУ МИФИ, 2012, с. 55.

6. В.И. Деев, Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин. Методика расчета теплоотдачи при кипении теплоносителя в каналах ЯЭУ // XVII школа-семинар по проблемам физики реакторов «Волга-2012», 3-8 сентября 2012 г., заседание 7, сборник презентаций докладов в электронном виде.

7. Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин. Режимы теплообмена с поверхности текстурированных Ni-W лент при их различной ориентации в поле силы тяжести // Научная сессия НИЯУ МИФИ-2013. Аннотации докладов, т. 1. Инновационные ядерные технологии. Высокие технологии в медицине. М.: НИЯУ МИФИ, 2013, с. 66.

8. В.И. Деев, Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин. Методика расчета нестационарного критического теплового потока при импульсном тепло-подводе в жидкому азоту // Научная сессия НИЯУ МИФИ-2013. Аннотации докладов, т. 1. Инновационные ядерные технологии. Высокие технологии в медицине. М.: НИЯУ МИФИ, 2013, с. 69.

Подписано в печать:

01.11.2013

Заказ № 9013 Тираж - 100 экз. Печать трафаретная. Типография «11-й ФОРМАТ» ИНН 7726330900 115230, Москва, Варшавское ш., 36 (499) 788-78-56 \vvvw. autoreferat.ru

Текст научной работы диссертации и автореферата по физике, кандидата физико-математических наук, Зар Ни Аунг, Москва

Национальный исследовательский ядерный университет

«МИФИ»

На правах рукописи

04201451544

Зар Ни Аунг

ЗАКОНОМЕРНОСТИ ТЕПЛООТДАЧИ И КРИЗИСА КИПЕНИЯ В ВОДЕ, НЕДОГРЕТОЙ ДО ТЕМПЕРАТУРЫ НАСЫЩЕНИЯ

специальность 01.04.14 - «Теплофизика и теоретическая теплотехника»

Автор:

Диссертация на соискание ученой степени кандидата физико-математических наук

Научный руководитель:

кандидат физико-математических наук,

доцент Лаврухин А.А.

Москва 2013

СОДЕРЖАНИЕ

ОСНОВНЫЕ ОБОЗНАЧЕНИЯ И СОКРАЩЕНИЯ............................4

ВВЕДЕНИЕ......................................................................................................8

ГЛАВА 1 ОБЩЕЕ СОСТОЯНИЕ ПРОБЛЕМЫ ИЗУЧЕНИЯ ТЕПЛООБМЕНА В ВОДЕ НЕДОГРЕТОЙ ДО ТЕМПЕРАТУРЫ НАСЫЩЕНИЯ..............................................................................................14

1.1 Стационарный теплообмен при кипении воды с недогревом........14

1.1.1 Теплообмен при пузырьковом кипении............................................14

1.1.2 Кризис теплоотдачи при пузырьковом кипении..............................36

1.2 Нестационарный кризис кипения воды............................................46

Выводы и постановка задачи исследования...............................................58

ГЛАВА 2 МЕТОДИКА ЭКСПЕРИМЕНТАЛЬНОГО ИССЛЕДОВАНИЯ........................................................................................60

2.1 Экспериментальная установка...........................................................60

2.2 Рабочие участки..................................................................................60

2.3 Методика проведения экспериментов

при стационарном теплоподводе.................................................................64

2.4 Методика проведения экспериментов при импульсном

увеличении мощности тепловыделения......................................................68

ГЛАВА 3 ТЕПЛООБМЕН В ВОДЕ ПРИ СТАЦИОНАРНОМ ПОДВОДЕ ТЕПЛА.......................................................................................73

3.1 Кривые кипения насыщенной и недогретой до температуры насыщения воды под атмосферным давлением.........................................73

3.2 Влияние недогрева воды до температуры насыщения и диаметра нагревателя на первый критический тепловой поток................................80

3.3 Статистический анализ флуктуаций температуры

теплоотдающей поверхности.......................................................................86

3.3.1 Результаты экспериментального исследования...............................86

3.3.2 Обоснование возможности использования плотности распределения температурных шумов для диагностики смены

режимов теплообмена...................................................................................95

ГЛАВА 4 МЕТОДИКА РАСЧЕТА ТЕПЛООБМЕНА В ПЕРЕХОДНЫХ РЕЖИМАХ ПРИ КИПЕНИИ ЖИДКОСТИ НА ПОВЕРХНОСТИ НАГРЕВА......................................................................100

4.1 Модель теплообмена и методика расчета

коэффициента теплоотдачи........................................................................100

4.2 Пример расчета температуры оболочек твэлов

ядерного реактора........................................................................................110

ГЛАВА 5 ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ПРОЦЕССА ТЕПЛООБМЕНА В ВОДЕ ПРИ ИМПУЛЬСНОМ ПОДВОДЕ ТЕПЛА.....................................................................................113

5.1 Влияние величины набрасываемой тепловой нагрузки

и недогрева воды до температуры насыщения.........................................113

5.2 Расчет параметров нестационарного кризиса кипения воды........120

ОСНОВНЫЕ ВЫВОДЫ..............................................................................127

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

ОСНОВНЫЕ ОБОЗНАЧЕНИЯ И СОКРАЩЕНИЯ

ЛАТИНСКИЕ СИМВОЛЫ А коэффициерт;

а коэффициент температуропроводности, м2/с;

Ь = ТоТ^р^р^) постоянная Лапласа, м;

ср удельная теплоемкость при постоянном давлении

теплоемкость, Дж/(кг-К); Д,хр отрывной диаметр, м;

с1 диаметр, м;

Е энергия, Дж;

/ частота отрыва пузырей, с"1;

g ускорение свободного падения, м/с ;

к высота, глубина, шероховатость, м;

/ сила тока, А;

К критерий устойчивости;

I длина, м;

М* молярная масса, г/моль;

т, п показатели степени;

р давление, Па;

(2о количество теплоты, выделенное источником на

единицу поперечного сечения, Дж/м2; д плотность теплового потока, Вт/м ;

дн набрасываемая тепловая нагрузка, Вт/м2;

удельное объемное тепловыделение, Вт/м ; Я радиус парового пузыря, м;

универсальная газовая постоянная, Дж/(моль К); Дэ, /?обр электрическое, образцовое сопротивление, Ом; г скрытая теплота парообразования, Дж/кг; г0 радиус устья впадины, м; £ ---------- доля поверхности нагрева; -

Г температура, К;

АТ температурный напор, К;

и падение напряжения, В;

V объем; м ;

>уген генерируемая в нагревателе тепловая мощность, Вт/м ;

х относительное балансное паросодержание.

ГРЕЧЕСКИЕ СИМВОЛЫ а коэффициент теплоотдачи,

Вт/(м -К);

Р коэффициент объемного расширения, К"1;

у угол полураствора;

5 ширина или толщина, м;

0 = Т5 - Тж недогрев жидкости относительно температуры насыщения, К;

Л критическая длина волны, м;

X коэффициент теплопроводности, Вт/(м-К);

ц динамический коэффициент вязкости, Па с;

V коэффициент кинематической вязкости, м /с;

р плотность, кг/м ;

рэ удельное электрическое сопротивление, Ом-м;

а коэффициент поверхностного натяжения, Н/м;

т время, с.

СОКРАЩЕНИЯ БОВ начало развитого кипения;

ОЫВ начало пузырькового кипения;

ОБУ начало интенсивной генерации пара.

ИНДЕКСЫ

' жидкость на линии насыщения;

" пар на линии насыщения;

ё диаметр;

Б стенка, поверхность;

б насыщение;

выс высокий;

ген генерируемый;

гид гидравлический;

ж жидкость;

з закипание;

зап запасенный;

исп испаренный;

кип кипение;

конв конвекция;

кр критический;

кр1 первый кризис кипения;

кр2 второй кризис кипения;

макс максимальный;

мин минимальный;

н нагреватель;

нест нестационарный;

низ низкий;

обр образцовое;

отр отрывной;

п пленочный;

пар пар;

пр предельный;

расч расчетный;

с л слияние;

ср

эксп эф

средний;

экспериментальный; эффективный.

БЕЗРАЗМЕРНЫЕ КРИТЕРИИ ПОДОБИЯ

1а. = ср р' Д77 (г р") число Якоба;

N11 = а с1 / X число Нуссельта;

Рг = V / а число Прандтля;

11а = |3 § яР Д77 V а число Рэлея;

Яа = g о?3 (р' - р") / р" V а число Рэлея для пленочного кипения;

Ке = м; / / V число Рейнольдса.

Остальные условные обозначения поясняются в тексте.

ВВЕДЕНИЕ

Актуальность темы диссертации

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

• кипение является очень эффективным режимом отвода тепла и поэтому широко используется в теплообменниках, энергетическом оборудовании, для охлаждения энергонапряженных элементов различных устройств;

• кипение является чрезвычайно сложным процессом, на который оказывают влияние столь большое количество факторов, что создание строгой теории невозможно.

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

С учетом изложенного, диссертационная работа, сориентированная на получение новых методик и зависимостей, описывающих процессы

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

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

Научная новизна работы

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

Автор защищает

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

• Новые опытные данные и их обобщение, описывающее влияние недогрева воды до температуры насыщения на характеристики кризиса кипения при быстром возрастании мощности тепловыделения.

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

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

Полученные в работе результаты внедрены в учебный процесс НИЯУ МИФИ на кафедре теплофизики, а также вошли в научно-технические отчеты по следующим проектам Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России»: «Поиск предельных значений температуры и паросодержания теплоносителя ВВЭР-1200 для повышения КПД реакторной установки с учетом эксплуатационных характеристик материалов активной зоны» (ГК № 02.740.11.0046); «Ядерные реакторы нового поколения: физическое и математическое моделирование в обоснование модернизации ВВЭР и оптимизации характеристик быстрых реакторов для этапа перехода к замкнутому топливному циклу» (ГК № 14.740.11.0157).

Обоснованность и достоверность

Достоверность полученных экспериментальных данных подтверждается: отработкой методики измерений в тестовых экспериментах, воспроизводимостью опытных данных, оценками величин погрешности измерений. Эмпирические уравнения и обобщающие зависимости согласуются с имеющимися в литературе данными для соответствующих диапазонов режимных параметров. Обоснованность предложенных методик подтверждается сравнением с известными представлениями о механизмах теплообмена при кипении и экспериментальными данными.

Апробация

Материалы диссертационной работы, докладывались на следующих научных сессиях семинарах: Научные сессии НИЯУ МИФИ 2011 - 2013; XVII школа-семинар по проблемам физики реакторов «Волга-2012».

Публикации по теме диссертации Статьи в реферируемых журналах из перечня ВАК:

1. В.И. Деев, Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин, Ю.А. Маслов, В.Н. Федосеев, B.C. Харитонов Методика расчета теплообмена при кипения жидкости на поверхности нагрева // Ядерная физика и инжиниринг, 2011, том 2, № 5, с. 387-394.

2. В.И. Деев, Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин, В.Н. Федосеев Статистический анализ температурных флуктуаций как метод диагностики режимов теплообмена при кипении // Тепловые процессы в технике, 2013, Т. 5, № 4, с. 163-169.

3. В.И. Деев, Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин Динамические характеристики процесса теплообмена в воде при импульсном подводе тепла // Естественные и технические науки, 2013, №3 (65), с. 35-41.

Работы, опубликованные в других изданиях:

1. В.И. Деев, К.В. Куценко, A.A. Лаврухин, Зар Ни Аунг Новый метод расчета теплоотдачи в режиме кипения воды с недогревом //Научная сессия НИЯУ МИФИ-2011. Аннотации докладов. Т. 1. Инновационные ядерные технологии. М.: НИЯУ МИФИ, 2010.-с. 46.

2. В.И. Деев, Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин Методика расчета теплоотдачи к недогретой до температуры насыщения воде в предкризисной области //Научная сессия НИЯУ МИФИ-2012. Аннотации докладов. Т. 1. Инновационные ядерные технологии. Высокие технологии в медицине. М.: НИЯУ МИФИ, 2012,- с. 55.

3. В.И. Деев, Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин Методика расчета теплоотдачи при кипении теплоносителя в каналах ЯЭУ // XVII школа-семинар по проблемам физики реакторов, «Волга-

2012», 3-8 сентября 2012 г., заседание 7, сборник презентаций докладов в электронном виде.

4. Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин Режимы теплообмена с поверхности текстурированных Ni-W лент при их различной ориентации в поле силы тяжести //Научная сессия НИЯУ МИФИ-2013. Аннотации докладов. Т. 1. Инновационные ядерные технологии. Высокие технологии в медицине. М.: НИЯУ МИФИ, 2013,- с. 66.

5. В.И. Деев, Зар Ни Аунг, К.В. Куценко, A.A. Лаврухин Методика расчета нестационарного критического теплового потока при импульсном теплоподводе в жидкому азоту // Научная сессия НИЯУ МИФИ-2013. Аннотации докладов. Т. 1. Инновационные ядерные технологии. Высокие технологии в медицине. М.: НИЯУ МИФИ, 2013.- с. 69.

Структура и объем диссертации Диссертация состоит из введения, 5 глав, заключения и списка литературы, включающего 119 наименований. Содержание диссертации изложено на 141 странице машинописного текста, включая 55 рисунков и 1 таблицу к основному тексту.

Личное участие автора Данная работа выполнялась с 2010 г. по 2013 г. на кафедре теплофизики НИЯУ МИФИ. Постановка задачи исследования осуществлена совместно с научным руководителем. Экспериментальные исследования, обобщение и анализ полученных результатов, разработка моделей и методик выполнены автором самостоятельно или в соавторстве при его непосредственном участии.

ГЛАВА 1 ОБЩЕЕ СОСТОЯНИЕ ПРОБЛЕМЫ ИЗУЧЕНИЯ ТЕПЛООБМЕНА В ВОДЕ, НЕДОГРЕТОЙ ДО ТЕМПЕРАТУРЫ НАСЫЩЕНИЯ

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

В данной главе представлен анализ российских и зарубежных работ посвященных закономерностям кипения воды, недогретой до температуры насыщения, а также особенностям протекания процессов теплообмена при нестационарных тепловых воздействиях.

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

1.1 Стационарный теплообмен при кипении воды с недогревом 1.1.1 Теплообмен при пузырьковом кипении

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

Рис. 1.1 Схема режимов теплообмена в воде при вынужденном течении [5].

называемых центрах парообразования. Для активации этих центров температура поверхности, Т, должна превышать температуру насыщения жидкости при заданном давлении, Т5. Если температура жидкости, Тж, при этом ниже температуры насыщения, то процесс называется кипением с недогревом [4].

На рис. 1.1 схематично показаны основные области теплообмена к воде при ее течении в трубе с заданным тепловым потоком на стенке д [5]. Недогрета