Численное исследование образования и эволюции пузырей пара в условиях падения давления жидкости тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

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

КУМЗЕРОВА Екатерина Юрьевна

ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ОБРАЗОВАНИЯ И ЭВОЛЮЦИИ ПУЗЫРЕЙ ПАРА В УСЛОВИЯХ ПАДЕНИЯ ДАВЛЕНИЯ ЖИДКОСТИ

Специальность 01.02.05 - Механика жидкости, газа и плазмы

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

Санкт-Петербург 2004

Работа выполнена в Секторе численного моделирования Физико-технического института им. А.Ф. Иоффе РАН.

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

ст. научный сотрудник Шмидт Александр Александрович.

Официальные оппоненты: доктор технических наук,

ведущий научный сотрудник Готовский Михаил Абрамович.

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

научный сотрудник Травин Андрей Константинович.

Ведущая организация: Балтийский государственный технический университет

«Военмех»

Защита состоится « » 2004 года в /¿Г часов на заседании

диссертационного совета Д • 212.229.07 в ГОУ ВПО «Санкт-Петербургский государственный политехнический университет» (195251, Санкт-Петербург, Политехническая ул., 29, корпус 1, кафедра гидроаэродинамики).

С диссертацией можно ознакомиться в фундаментальной библиотеке ГОУ ВПО «Санкт-Петербургский государственный политехнический университет»

Автореферат разослан «<?» ^ол . 2004 года.

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

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

Актуальность проблемы

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

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

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

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

рос. национальная! БИБЛИОТЕКА 1

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

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

Цели работы

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

1. Разработаны математическая модель и численный метод для описания вскипания жидкости под действием падения давления;

2. Проведен анализ влияния механизмов нуклеации на образование паровой фазы;

3. Исследовано влияние параметров несущей и дисперсной фаз, а также формы и размеров каналов на структуру вскипающих потоков.

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

В диссертационной работе впервые:

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

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

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

4. Показано, что учет распределения ядер по размерам существенен для определения параметров несущей и дисперсной фазы.

Достоверность полученных результатов

Достоверность работы подтверждается тем, что полученные в ходе тестирования модели и алгоритма результаты совпадают с известными экспериментальными данными. В частности, результаты расчетов были сопоставлены с экспериментальными данными по глубине провала давления в задаче о разгерметизации трубы, заполненной водой при больших давлении и температуре, (Alamgir et al..// Journal of Heat Transfer. -1980. -V.102. -N.3. -P.433-438; Alamgir et al. // Journal of Heat Transfer. -1981. -V.103. -N.1. -P.52-55) и критическим расходам двухфазных течений в соплах (Abuaf et al. / Experimental. Reports NUREG/CR-1864, BNL-NUREG-51317 -1981; Xu et al. // Experimental Thermal and Fluid Science -1997. -Vol.14, -P. 166-173).

Практическая ценность работы

Практическая значимость диссертационной работы состоит в следующем:

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

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

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

Апробация работы

Основные результаты работы докладывались на восьми российских и международных конференциях и семинарах:

4ой Международной конференции по многофазным течениям (Нью-Орлеан, США, 2001); XIII Школе-семинаре молодых ученых и специалистов под руководством акад. РАН А.И.Леонтьева «Физические основы экспериментального и математического моделирования процессов газодинамики и теплообмена в энергетических установках» (Санк-Петербург, 2001); XIX Международном семинаре по струйным, отрывным и нестационарным течениям (Санкт-Петербург, 2002); III российской национальной конференции по теплообмену (Москва, 2002), Русско-японском международном симпозиуме по актуальным проблемам вычислительной механики (Санкт-Петербург,

2002); XIV Школе-семинаре молодых ученых и специалистов под руководством акад. РАН А.И.Леонтьева «Физические основы экспериментального и математического моделирования процессов газодинамики и теплообмена в энергетических установках» (Рыбинск, 2003); Международном коллоквиуме Euromech Colloquium 447 (Таллинн,

2003); XII Международной конференции по вычислительной механике и современным прикладным программным системам (Владимир, 2003).

Публикации по теме диссертации

Основные результаты диссертации изложены в шести опубликованных работах.

Структура и объем работы

Диссертация состоит их четырех глав, заключения и списка литературы из 95 наименований. Работа изложена на 103 страницах машинописного текста, включая 47 рисунков.

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

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

• - адекватное моделирование вскипающих потоков должно учитывать как

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

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

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

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

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

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

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

предположении отсутствия проскальзывания между фазами.

5

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

где г - вектор консервативных переменных. РХ,РГ - потоки, / - источниковый член:

(2)

Здесь а - объемное содержание пара, индексы I и v относятся, соответственно,

к несущей и дисперсной фазе, р и р - давление и истинная плотность,

- удельная полная энергия несущей фазы, - численная плотность

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

определению величин, описывающих межфазные взаимодействия (обмен массой,

импульсом и энергией) и скорость зародышеобразования Я. Интенсивности

межфазного массообмена и межфазного переноса энергии определяются с

помощью соотношения Герца-Кнудсена. Определение скорости нуклеации //, то есть

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

6

Эйлеров этап замыкает уравнение состояния воды в форме Тейта:

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

(4)

Модели нуклеации, связанные с определением скорости образования зародышей Я, представлены в разделе 2.7. В работе рассматриваются три типа нуклеации — объемная гомогенная, объемная гетерогенная и пристенная. В разделе 2.7.1 описаны основные моменты теории гомогенной нуклеации, согласно которой скорость нуклеации пропорциональна вероятности возникновения критического кластера вследствие термодинамической флуктуации (Френкель, 1957):

(5)

Н = Н0 ехр

{"ад}'

где И^ - работа по образованию критического зародыша, кь - постоянная Больцмана.

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

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

Ключевым вопросом моделирования

1.6

1.2

0.8

0.4

5

1 - Иогта!

2 - Ьодпогта!

3 - ЕдшргоЬаЫе

2 Бона-О,=1.0

А\ ОавЬес! = 5.0

1 • \ М, = 1.0

1 \\3

■ ■>• • "V. ■ ........"

-1

1

г/г.

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

Рис.1. Вид функций различных распределений ПЛОЩАДИ ПОД кривой ПЛОТНОСТИ

зародышей по размерам. вероятности (от текущего критического

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

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

определяемой через эти величины скорости нуклеации на стенке

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

В третьей главе представлено описание численного метода.

Уравнения сохранения для несущей фазы совместно с уравнением нуклеации решаются с помощью метода Родионова второго порядка точности по времени и пространству в областях гладких решений, рассмотренного в разделе 3 1. Раздел 3 2 посвящен решению задачи о распаде произвольного разрыва для воды с уравнением состояния в форме Тейта, что используется для определения потоков массы и импульса на гранях вычислительной ячейки. Обыкновенные дифференциальные уравнения динамики пузырьков решаются одновременно с уравнениями для жидкости. Возможные методы решения уравнений лагранжева этапа рассматриваются в разделе 3.3. Из-за существенно различных характерных времен процессов в жидкости и в пузырях уравнения динамики пузырей являются жесткими. Показано, что наиболее эффективным в данном случае является применение метод Гира для жестких систем обыкновенных дифференциальных уравнений.

Для реализации алгоритма используется структурированная неравномерная расчетная сетка. Увеличение точности достигается путем увеличения количества ячеек в областях больших градиентов полей физических величин, в разделе 3 4 показана сходимость по сетке.

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

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

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

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

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

С целью проверки разработанного алгоритма и сделанных предположений было проведено сравнение результатов моделирования с известными экспериментальными данными, результаты которых можно видеть на Рис. 2 и 3. На Рис.2 представлено изменение во времени давления жидкости в сечении х=0.43т для моделей модифицированной гомогенной нуклеации и объемной гетерогенной с нормальным распределением зародышем по размерам и экспериментальных данных (Alamgir et э1. , 1980). Видно, что обе модели дают хорошие совпадение по минимальному давлению, достигаемому жидкостью, однако последующее увеличение давления, связанное с ростом возникших зародышей, лучше описывается второй моделью, что можно объяснить учетом распределения зародышей по размерам.

Разницу между давлением насыщения и минимальным давлением принято называть глубиной провала давления, которая является важной характеристикой для таких задач и определение которой представляет одну из важнейших задач для практических приложений. Расчетные зависимости глубины провала давления от температуры при разных скоростях снижения давления и корреляции экспериментальных данных для этой величины (Alamgir et al., 1981) представлены на Рис.3.

Рис 2 Изменение во времени давления жидкости в Рис 3 Зависимость глубины провала давления от

сечении трубы х=0 43т для модифицированной начальной температуры жидкости. Линии - расчет,

модели гомогенной нуклеации (штриховая линия) и точки - эксперимент (Alamgir et al., 1981), 1, 2, 3

модели, учитывающей нормальное распределение соответствуют скоростям падения давления 200

зародышей по размерам (сплошная линия) Точки • GPa/s, 100 GPa/s и 10 GPa/s экспериментальные данные (Alamgir et al., 1980).

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

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

Рис. 4 и 5 демонстрируют влияние выбора функции распределения. Рис 4 показывает, что давление жидкости для равновероятного распределения (модель 3), завышено по сравнению с результатами моделей, использующих нормальное (модель 1) и логнормальное (модель 2) распределения. Причины такого поведения связаны, по-видимому, со ступенчатым характером функции распределения. Из Рис 5 можно видеть, что максимальные объемные доли пара близки для разных моделей, но также видно, что для равновероятного распределения область кипения более локализована, что является прямым следствием ступенчатого характера функции плотности вероятности.

Рис 4 Распределения давления жидкостей для Рис 5 Распределения объемной доли паровой разных моделей. фазы для разных моделей

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

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

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

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

На рис.6 и 7 показано сравнение с

экспериментальными данными (Abuaf et

al., 1981; Xu et al., 1997) для течений в

соплах. Рис.6 демонстрирует влияние

механизмов пристенной и объемной

гетерогенной нуклеации на давление

жидкости. Можно видеть, что действие

только одного из указанных механизмов

приводит к недооцениванию

Рис 6. Сравнение с экспериментальными данными

(Abuaf et al, 1981) осредненных вдоль сопла минимального давления, достигаемого профилей давления для развитой нуклеации

400000

0 01 02 03 04 05

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

осредненный вдоль сопла профиль давления. Сопоставления по

максимальному расходу дало следующие результаты: 8.71 kg/s согласно

экспериментальным данным, 8.64kg/s -рассчитанное по модели (т.е.

относительная разница 0.8%), что еще раз Рис.7. Сравнение с экспериментальными

данными (Xu et al.. 1997) по критическим

подтверждает применимость приведенных расходам. модели и алгоритма.

На Рис.7 представлено сравнение критических расходов с экспериментальными данными (Xu et al., 1997). Разница между данными расчета и эксперимента не превышает 7%. Таким образом, можно сделать вывод, что сформулированная модель применима для рассматриваемых двухфазных течений, а для адекватного моделирования вскипания в жидкости необходим учет, как объемной гетерогенной, так и пристенной нуклеации.

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

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

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

Рис 8 Распределения численной концентрации пузырей и объемной доли пара по каналу для Тщ-480К иТи=500К

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

В заключении кратко сформулированы основные результаты настоящей работы

• Проанализированы физические явления, происходящие при вскипании жидкости под действием резкого понижения давления;

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

• Численное моделирование позволяет определить зависимость от времени и от пространственной координаты параметров несущей фазы, объемного содержания, температуры и давления паровой фазы, концентраций и размеров пузырей;

• Исследовано влияние различных механизмов нуклеации (модифицированной гомогенной, объемной гетерогенной и пристенной) на структуру вскипающего потока;

• Выявлено, что для описания объемной гетерогенной нуклеации учет распределения зародышей по размерам существенен для определения параметров несущей и дисперсной фаз;

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

Основные результаты диссертации опубликованы в работах:

1. Е.Ю.Кумзерова, ААШмидт. Численное исследование поведения пузырьков при быстром изменении давления жидкости // Труды конференции «Физические основы экспериментального и математического моделирования процессов газодинамики и теплообмена в энергетических установках», Санкт-Петербург, -2001. -ТОМ.2.-С.130-134.

2. Е.Ю.Кумзерова, А.А.Шмидт. Численное моделирование нуклеации и динамики пузырьков при быстром падении давления жидкости // Журнал Технической Физики. -2002. -Том 72. -Вып.7. -С.36-40.

3. Е.Ю.Кумзерова, А.А.Шмидт. Численное исследование влияния модели образования паровой фазы на структуру вскипающих потоков // Труды III российской национальной конференции по теплообмену, Москва. -2002. -Том 4. -С.131-135.

4. Е.Ю.Кумзерова, А.А.Шмидт. Механизмы нуклеации при распространении волны разрежения в жидкости // Тезисы докладов XII Международной конференции по вычислительной механики и современным прикладным программным системам, Владимир. -2003. -С.404-406.

5. Е.Ю.Кумзерова. О влиянии распределения центров нуклеации по размерам на моделирование вскипания жидкости // Труды конференции «Физические основы экспериментального и математического моделирования процессов газодинамики и теплообмена в энергетических установках», Рыбинск, -2003. -Том 1 .-С.259-262.

6. E.Yu. Kumzerova and A.A. Schmidt. Effect of bubble nucleation mechanisms on flashing flow structure (Numerical Simulation) // Computational Fluid Dynamics. -2003. -Vol.11. -No.4.-P.507-512.

Лицензия ЛР №020593 от 07.08.97.

Подписано в печать ОЦ.&ООУ. Объем в п.л. Тираж /00. Заказ

Отпечатано с готового оригинал-макета, предоставленного автором, в типографии Издательства СПбГПУ 195251, Санкт-Петербург, Политехническая ул., 29.

Отпечатано на ризографе RN-2000 ЕР Поставщик оборудования — фирма "Р-ПРИНТ" Телефон: (812) 110-65-09 Факс: (812) 315-23-04

»-93 3 7

МОДЕЛИРОВАНИЯ ВСКИПАЮЩИХ ПОТОКОВ

1.1. Изучаемые физические процессы

1.2. Исследования вскипающих потоков

1.3. Обзор литературы по моделированию нуклеации

1.4. Выводы из обзора литературы

1.5. Цели данной работы

ГЛАВА 2. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ТЕЧЕНИЯ ВСКИПАЮЩЕЙ ЖИДКОСТИ

2.1. Физические явления, сопровождающие течения жидкости в условиях резкого падения давления

2.2. Формулировка рассматриваемых задач динамики вскипающей жидкости

2.3. Основные допущения

2.4. Лагранжево-эйлеровский подход к описанию двухфазной системы.

2.5. Уравнения эйлеровского этапа

2.5.1. Описание несущей фазы

2.5.2. Уравнение возникновения и конвективного переноса пузырей

2.5.3. Математические модели процессов межфазного переноса

2.6. Уравнения лагранжева этапа - описание динамики пузырьков

2.6.1. Уравнение сохранения массы и внутренней энергии

2.6.2. Уравнение движения межфазной границы — уравнение Рэлея-Лэмба

2.6.3. Некоторые замечания, связанные с влиянием теплопроводности жидкости на эволюцию пузыря

2.7. Модели нуклеации

2.7.1. Теория гомогенной нуклеации

2.7.2. Модифицированные модели гомогенной нуклеации

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

2.7.4. Модель пристенной нуклеации

2.8. Полная система уравнений, начальные и граничные условия, замыкание модели

2.8.1. Полная система уравнений эйлеровского этапа для осесимметричного случая

2.8.2. Полная система уравнений лагранжева этапа

2.8.3. Замыкание модели '

2.8.4. Начальные и граничные условия

ГЛАВА 3. ЧИСЛЕННЫЙ МЕТОД

3.1. Решение эйлеровского этапа

3.1.1. Схема Родионова

3.1.2. Задача о распаде разрыва для воды

3.2. Решение лагранжева этапа

3.2.1. Метод Адамса

3.2.2. Метод Гира

3.2.3. Сравнение методов

3.3. Расчетная сетка

3.4. Сходимость численной схемы на расчетной сетке

3.5. Тестирование численного метода

3.5.1. Одномерная задача о распаде разрыва

3.5.2. Двумерная задача о распаде разрыва

ГЛАВА 4. РЕЗУЛЬТАТЫ

4.1. Одномерные задачи

4.1.1. «Разгерметизация трубы»

4.1.2. Задача о закрытии клапана

4.2. Двумерные задачи

4.2.1. Течение в канале переменного сечения

4.2.2. Задача о «внезапном расширении»

РИСУНКИ К ГЛАВЕ

Основные результаты диссертационной работы состоят в следующем:

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

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

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

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

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

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

ЗАКЛЮЧЕНИЕ

1. Richter S.t Fleischer S.,Aritomi M., Hampel R. Transient two-phase flow in arbitrary inclined tubes caused by depressurization of liquid with dissolved gases // Int. J. Heat and Mass Transfer. -2001. -Vol.44, -P. 1-15.

2. Deligiannis P., Cleaver J. W. The role of nucleation in the initial phases of rapid depressurization of a subcooled liquid // Int. J. Multiphase Flow. 1990. -N 16. -P.468-475.

3. Xu J.L., Chen Т.К., Chen X.J. Critical flow in convergent-divergent nozzles with cavity nucleation model // Experimental thermal and fluid science. -1997. -N 14. -P. 166-173.

4. Hahne E., Barthou G. Evaporation waves in flashing processes// Int. J. Multiphase Flow. -2000. -N 26. -P.531-547.

5. Kumzerova E.Yu., Schmidt A.A. Effect of bubble nucleation mechanisms on flashing flow structure (Numerical Simulation) // Computational Fluid Dynamics. -2003. -Vol.11. -N 4. -P.507-512.

6. Hsu Y.Y. Review of critical flow, propagation of pressure pulse, and sonic velocity / Report NASA TND-6814. 1972. -38p.

7. Saha P. Review of two phase stream-water critical flow models with emphasis on thermal nonequilibrium / Reports BUREG/CR 0417& BNL-NUREG-50907, 1978. -35p.

8. Abdollahian D., Healzer J., Janssen E., Amos C. Critical flow data review and analysis / Report EPRI NP-2192. 1982. -46p.

9. Bailey J.F. Metastable flow of saturated water // Trans. ASME 73. -1951. -P.1109-1116.

10. Brown R.A. Flashing expansion of water through a converging-diverging nozzle. M.S. Thesis, Univ. of California, Berkley CA.UKAEC Report UCRL-6665-T. 1961. -167p.

11. Schrock V.E., Starkman E.S., Brown R.A. Flashing flow of initially subcooled water in converging-diverging nozzles // Journal of Heat Transfer. -1977. -V. 99. -P.263-268.

12. Powell A.W. Flow of subcooled water through nozzles / Report WAPD-PT-(V)-90, Weistinghouse Electic Corp.1961. -51 p.

13. Abuaf N., Jones 0,C, Wu Jr. Critical flashing flow in nozzles with subcooled inlet conditions / Polyphase flow and transport technology (edited by Bajura R.A.), ASME, New York, 1980. -P.65-74.

14. Abuaf N., Jones 0,C, Wu Jr. Critical flashing flow in nozzles with subcooled inlet conditions // Journal of Heat Transfer. -1983. -V.105. -P.379-383.

15. Alamgir Md., Lienhard J. H. Correlation of pressure undershoot during hot-water depressurization // Journal of Heat Transfer. -1981. -V. 103. -N 1. -P.52-55.

16. Simpson H.C., Silver R.S. Theory of one-dimensional two-phase homogeneous nonequilibrium flow // Proc. Inst. Mech. Engrs. Symp. on Two-phase Flow. 1962. -P.45-53.

17. Plesset M.S. Zwick S.A. The growth of vapor bubbles in superheated liquids //Journal Applied Physics. -1954. -P.493-500.

18. Edwards A.R. Conduction controlled flashing of a flui and the prediction of critical flow rates in a one-dimensional system / Report AHSB (S) R-147, UKAEA AERE.-1968. -69p.

19. Moody F.J. Maximum two-phase vessel blowdown from pipes // Journal of Heat Transfer. -1966. -V. 88. -P.285-295.

20. Henry R.E. Fauske H.K. The two-phase critical flow of one component mixtures in nozzle orifices, and short tubes // Journal of Heat Transfer. -1971. -V. 93.-P.179-187.

21. Jones O.C. Saha P. Nonequilibrium aspects of water reactor safety // In Thermal Hydraulic Aspects of Nuclear Reactor Safety. Vol 1, Light Water Reactors (Edited by Jones O.C., Bankoff S.G.) ASME, New York. 1977. -18p.

22. Riebold W.L. Reocreux M., Jones O.C. Blowdown phase / In Unclear Reactor Safety Heat Ransfer. Hemisphere, New York. 1981. -P. 325-378.

23. Jones O.C. Toward a unified approach for thermal nonequilibrium in gas-liquid systems // Nuclear Engineering Design. -1982. -Vol.69. -P.57-73.

24. Malnes D. Critical two-phase flow based on nonequilibrium model / In Nonequilibrium pwo-phase flow (Edited by Lahev R.T., Wallis G.B.). ASME. New York. 1975. -P. 11-16.

25. Henry R.E. Fauske H.K. McComas S.T. Two-phase critical flow at low qualities. Part l:experimental // Nuclear Engineering Design. -1970. -Vol.41. -P.79-91.

26. Fauske H.K. The discharge of saturated water through tubes // Proc.of the 7th National Heat Transfer Conference. ASME, OH. 1964. -P.23-29.

27. Rohatgi U.S. Reshotko E. Nonequilibrium one-dimensional two-phase flow in variable area channels / In Nonequilibrium pwo-phase flow (Edited by Lahev R.T., Wallis G.B.). ASME. New York. 1975. -P. 47-54.

28. Fritz G. Riebold W. Scultze W. Studies on thermodymanic nonequilibrium in flashing flow / Presented on OECD/NEA Specialist MTG on Transient Two-Phase Flow, Toronto, Ontario. 1976.

29. Wolfert K. The simulation of blowdown processes with condensation of thermodynamic nonequilibrium phenomena / Presented at the OECD/NEA Specialists MTG on Transient Two-phase Flow. Toronto. Ontario. 1976.

30. Aguilar F. Thompson S. Nonequilibrium flashing model of rapid pressure transient / ASME preprint. Paper No. 81-HT-35. 1981. 41 p.

31. Ardron K.H. A two-fluid model for critical vapor-liquid flow // Int. J.Multiphase Flow. -1978. -Vol.4. -P.323-337.

32. Rivard W.C. Travis J.R. A nonequilibrium vapor production model for critical flow// Nuclear Engineering Design. -1980. -Vol. 74. -P. 40-48.

33. Richter H.J. Separated two-phase flow model: application to critical two-phase flow / Report EPRI NP-1800. 1981.-53p.

34. The Marviken full scale critical-flow tests / NUREG/CR-2671, MXC-301, 1982. -57p.

35. Reocreus M. Experimental study of steam-water choked flow / Presented on OECD/NEA Specialist MTG on Transient Two-Phase Flow, Toronto, Ontario. 1976.

36. Sozzi G.L. Sutherland W.A. Critical flow of saturated and subcooled water at high pressure / GE Report NEDO-13418. 1975. -35p.

37. Winters W. Merte H. Experiments and nonequilibrium analysis of pipe blowdown // Nuclear Sci. Engng. -1979. -Vol.69. -P.411-429.

38. Zimmer G.A., Wu B.J.C., Leonhard W.L., Abuaf N., Jones O.C. Pressure and void distributions in a converging-diverging nozzle with nonequilibrium water vapor generation / Report BNL-NUREG-26003. 1979.

39. Simoneau R.J. Pressure distribution in a converging-diverging nozzle during two-phase choked flow of subcooled nitrogen / In Nonequilibrium two-phase flows (Edited by Lahev R.T. and Wallis G.B.), ASME. New York. 1975. -P. 37-44.

40. Hendricks R.C., Simoneau R.J. Barrows R.F. Two-phase choked flow of subcooled oxygen and nitrogen / Report NASA-TD-D-8169. 1976, -68p.

41. Lienhard J. H., Alamgir Md., Trela M. Early response of hot water to sudden release from high pressure // Journal of Heat Transfer. -1977. -V.100. -N 3. -P.473-479.

42. Alamgir Md., Kan C.Y., Lienhard J. H. An experimental study of the rapid depressurization of hot water // Journal of Heat Transfer. -1980. -V. 102. -N 3. -P.433-438.

43. Jones O.C. Flashing inception in flowing liquids // Journal of Heat Transfer. -1980. -V.102. -N 3. -P.439-444.

44. Saha P. Abuaf N. Wu B.J.C. A nonequilibirium vapor generation model for flashing flows / ASME preprint. Paper No. 81-HT-84. 1981. -33p.

45. Levy S., Abdollahian D. Homogeneous nonequilibrium critical flows // Int. Journal of Heat and Mass Transfer. -1982. -Vol. 25. -P.759-778.

46. Albagli D., Gany A. High speed bubble nozzle flow with heat, mass, and momentum interactions // Int. J. Heat and Mass Transfer. -2003. -Vol.46. -P.1993-2003.

47. Cokliat D., Ivanov V.A., Vasquez S.A. Two-phase model for cavitating flows// Proc. of the 2th International Conference on Multiphase Flow 1998, Lion. 1998 (on CD).

48. Iben U., Wrona F., Munz C. -D., Beck M. Cavitation in hydraulic tools based on thermodynamic properties of liquid and gas // Journal of Fluids Engineering. -2002. -Vol. 124. -P.2-7.

49. Скрипов В.П. и др. Теплофизические свойства жидкостей в метастабильном состоянии. М., Атомиздат, 1980. -208с.

50. Lee S.Y. Schrock V.E. Critical two-phase flow in pipes for subcooled stagnation states with cavity flooding incipient flashing model // Journal of Heat Transfer. -1990.-Vol.112. -P.1032-1040.

51. Elias E. Chambre P.L. Flashing inception in water during rapid decompression // Journal of Heat Transfer. -1993. -Vol. 115. -P. 231-238.

52. Abuaf, N., Zimmer G.A., We B. J. C. A study of non-equilibrium flashing of water in a converging-diverging nozzle, Vol. 1 / Experimental. Reports NUREG/CR-1864, BNL-NUREG-51317. 1981.-52p.

53. Volmer V., Weber A. Keimbildung in uebersaetigen Daempfen. // Z.Phys.Chem. -1926. -N119. -P. 277-301.

54. Френкель Я.И. Кинетическая теория жидкостей. Л., Наука, 1959. -586С.

55. Каверин A.M., Байдаков В.Г., Скрипов В.П. Частота спонтанного зародышеобразования в перегретых жидких ксеноне и криптоне // Инженерно-физический журнал. -1980. -Том XXXVIII. -№4. -С.680-684.

56. Schmelzer J.W.P., Schmelzer J.Jr. Kinetics of nucleation at increasing supersaturation // J. Colloid and Interface Science. -1999. -Vol.215. -P.345-355.

57. Павлов П.А. Теплофизика метастабильных систем. Свердловск. УНЦ. 1977.-69с.

58. Павлов П.А. Динамика вскипания сильно перегретых жидкостей. Свердловск. УрО АН СССР. 1988. -74с.

59. Lienhard J. Н., Karimi A. Homogeneous nucleation and the spinodal line // Journal of Heat Transfer. -1981. -V. 103. -N 1. -P.61-64.

60. Kendoush A.A. The delay during depressurization of saturated water // Int. J. Heat and Mass Transfer. -1989. -Vol.32. -N11. -P.2149-2154.

61. Bartak J. A study of the rapid depressurization of hot water and the dynamics of vapor bubble generation in superheated water // Int. J. Multiphase Flow. -1990. -N 5. -P.789-798.

62. Elias E., Chambre P.L. Bubble transport in flashing flow // Int. J. Multiphase Flow. -2000. -N 26. -P.191-206.

63. Kumzerova E.Yu., Schmidt A.A. Numerical simulation of homogeneous nucleation and bubble dynamics in a depressurized liquid // Proc. of the 4th International Conference on Multiphase Flow 2001, New Orleans. 2001. (на CD).

64. Kwak H.-Y., Kim Y.W. Homogeneous nucleation and macroscopic growth of gas bubble in organic solutions // Int. J. Heat and Mass Transfer. -1998. -Vol. 41. -No.4-5. -P.757-767.

65. Kiselev S.B. Kinetic boundary of metastable states in superheated and stretched liquids// Physica A, -1999. -N 269. -P.252-268.

66. Wiesche S., Rembe C., Hofer E.P. Boiling of superheated liquids near.the spinodal: I general theory // Heat and mass transfer. -1999. -Vol.35. -P.25-31.

67. Kotake S., Glass I.I. Survey of flows with nucleation and condensation / UTIAS Review No.42, CN ISSN 0082-5247, Oct.1978. -39p.

68. Peng X.F et al. Cluster dynamics and fictitious boiling in microchannels // Int. J. Heat and Mass Transfer. -2000. -N 43. -P.4259-4265.

69. Bankoff S.G. Entrapment pf gas in the spreading of a liquid over a rough surface //AlChE Journal. -1951. -Vol. 4. -P.24-26.

70. Hsu Y.Y. On the size range of active nucleation cavities on a heating surface // Journal of Heat Transfer. -1962. -Vol 94. -P. 207-212.

71. Shoukri M.S.M. Judd R.L. A theoretical model for bubble frequency in nucleate pool boiling including surface effects / Presented at 6th Int. Heat Transfer Conference, Toronto, Ontario, 1978.

72. Jones O.C. Shin T.S. An active cavity model for flashing // Nuclear Engineering Design. -1986. -Vol. 95. -P. 185-189.

73. Kocamustafaogullari G. Ishii M. Interfacial arear and nucleation site density in boiling systems // Int. J. Heat and Mass Transfer. -1982. -Vol.26. -P.1377-1387.

74. Zuber N. Nucleate boilng the region of isolated bubbles: similarity with natural convection // Int.J. Heat Mass Transfer. -1963. -N 6. -P.53-78.

75. Хирс Д., Паунд Г. Испарение и конденсация. Пер. с англ. М., Металлургия. 1966. -462С.

76. Campos F.B., Lage P.L.C. Heat and mass transfer during the formation and ascension of superheated bubbles // Int. J. Heat and Mass Transfer. -2000. -Vol.43. -P.2883-2894.

77. Shin T.S., Jones O.C. Nucleation and flashing in nozzles-1. A distributed nucleation model // Int. J. Multiphase Flow. -1993. -Vol.19. -N 6. P.943-964.

78. Blinkov V.N., Jones O.C., Nigmatulin B.I. Nucleation and flashing in nozzles-2. Comparison with experiments using a five-equation model for vapor void development// Int. J. Multiphase Flow. -1993. -Vol.19. -N 6. -P.965-986.

79. Maksic S., Mewes D. Flashing due pipe break in a sealing water cooler// Proc. of the 4th International Conference on Multiphase Flow 2001, New Orleans. 2001 (на CD).

80. Kedrinskii V.K. The lordansky-Kogarko-van Wijngaarden model: shock and rarefaction wave interactions in bubbly media //Applied Scientific Research. -1998.-Vol. 58. -P. 115-130.

81. Нигматулин P.И. Основы механики гетерогенных сред. М., Наука, 1978, -335с.

82. Рахматулин Х.А., ПММ 20. №2. 1956. -38с.

83. Крайко А.Н., Стернин Л.Е. ПММ 29. №3. 1965. -79с.

84. Нигматулин Р.И. Динамика мнгофазных сред. Москва. Наука. 1987. Т.1. -359с.

85. Pan Y., Dudukovic М.Р., Chang М. Dynamic simulation of bubbly flow in bubble columns// Chemical Eng. Science, 54. -1999. -P.2481-2489.

86. Lapin A., Lubert A. Numerical simulation of the dynamics of two-phase gas-liquid flows in bubble columns // Chemical Eng. Science. -1994. -N 49. -P.3661-3674.

87. Sommerfeld, M., Bourloutski, E. & Broder, D. Euler/Lagrange calculations of bubbly flows with consideration of bubble coalescence // The Canadian Journal of Chemical Engineering. -2003. -Vol. 81. -P.508-518.

88. Sokolichin A., Eigenberger G., Lapin A., Lubbert A. Dynamic numerical simulation of gas-liquid two-phase flows: Euler/Eler versus Euler/Lagrange // Chemical Eng. Science. -1997. -Vol. 52. -P.611-626.

89. Maeno K., Kouchi Т., Sato H. Behavior of pulse-laser-induced cavitation bubble in liquid nitrogen // Proc. of the ISSW'21. -Marcel, 1995. (on CD).

90. Лабунцов Д.А., Ягов В.В. Механика двухфазных систем. М., «МЭИ». 2000. -369с.

91. Родионов А.В. Повышение порядка аппроксимации схемы С.К.Годунова //ЖВМиМФ. -1987. -Том.27. -№3. -С. 1853-1859.

92. Численное решение многомерных задач газовой динамики. Под ред. С.К. Годунова. М.,Наука. 1976. -400с.

93. Oran E.S., Boris J.P. Numerical simulation of reactive flows. Elsevier Science Publ., 1987. -655p.