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

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

1. Артемьева Н.В., Чернецкий М.Ю., Хрусталев Г.Н. Сравнительный анализ по математической модели аэродинамики, теплообмена и температурных полей двухвихревой топки с наличием и без заднего холостого дутья // Минеральная часть топлива, шлакование, очистка котлов, удаление и использование золы: Сборник докладов V научно-практ. конф. - Челябинск, 7-9 июня 2011. - Челябинск: ОАО «ИЦЭУ», 2011. - Т. 3. -С. 157-166.

2. Майданик М.Н., Вербовецкий Э.Х., Дектерев А.А., Чернецкий М.Ю., Гаврилов А.А., Бойков Д.В., Бердин С.В. Математическое моделирование топки и поворотного газохода котла П-50Р при совместном сжигании твердого и газообразного топлива // Теплоэнергетика. - 2011. - № 6. - С. 37-42

3. Siegel R., Howell J.R. Thermal Radiation Heat Transfer. - Washington D.C.: Hemisphere Publishing Corp., 1992. - 1072 p.

4. Magnussen B.F., Hjertager B.W. On the structure of turbulence and a generalized eddy dissipation concept for chemical reaction in turbulent flow // XIXth AIAA Aerospace Meeting. - St. Louis, USA, 1981. - P. 544-552

5. Чернецкий М.Ю., Дектерев А.А. Математическая модель процессов теплообмена и горения пылеугольного топлива при факельном сжигании // Физика горения и взрыва. - 2011. -№ 3. - С. 37-46.

6. Crow C.T., Sharma M.P., Stock D.E. The particle source in cell (PSI-CELL) model for gas-droplet flows // Journal of Fluids Engineering. - 1977. - V. 99. - P. 325-332.

7. De Soete G. Overall reaction rates of NO and N2 formation from fuel nitrogen // Proc. XVth Int. Symposium on Combust. - Pittsburgh: The Combustion Institute, 1975. - P. 1093-1102.

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

УДК 536.2:51-74

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

В.Ю. Половников, А.М. Хабибулин

Томский политехнический университет E-mail: polov@tpu.ru

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

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

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

Key words:

Loss of cold, mathematical modeling, thermal insulation, phase transitions, cryogenic liquid.

Введение

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

Одним из направлений реализации Стратегии

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

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

Тепловая изоляция ограждений холодильных сооружений эксплуатируется в тяжелых условиях

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

[3]. Накопление влаги в теплоизоляционных конструкциях низкотемпературного оборудования может привести к столь значительному увеличению потерь холода, что необходима будет их полная замена [4].

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

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

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

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

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

Постановка задачи

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

Внешний контур изоляции резервуара контактирует с влажным воздухом (рис. 1), а влага, содержащаяся в нем, конденсируется на поверхности рассматриваемого объекта. Влагоперенос происходит только в талой (увлажненной) зоне путем миграции влаги к фронту промерзания за счет пленочно-диффузионного механизма [4]. Для задачи теплопроводности на внутренней Я1 и внешней Л2 границах слоя тепловой изоляции вводятся граничные условия первого и третьего рода соответственно. Для задачи переноса влаги на границе фазового перехода £ выставляется условие идеальной гидроизоляции, а на внешнем контуре взаимодействия Л2 - условие массообмена поверхности с окружающей средой (влажный воздух). В начальный момент времени температура Т0 и относительное влагосодержание изоляции по объёму Ж0 равны постоянным величинам. Предполагается, что температура Т0 равна температуре окружающей среды Тос, а температура внутренней поверхности изоляции равна температуре холодной жидкости Тж, которая ниже температуры замерзания Тпл. В результате охлаждения изоляции образуется промерзший слой (рис. 1) переменной толщины, причем подвижная граница фазового перехода всегда имеет постоянную температуру замерзания Тпл. Предполагается, что перенос тепла в изоляции происходит только вследствие теплопроводности.

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

Математическая модель

Математическая постановка задачи базируется на апробированных подходах к моделированию физических процессов [2, 7, 8], ив предлагаемой постановке процессы переноса тепла и массы в рассматриваемой области решения (рис. 1) будут описываться следующими уравнениями:

ЗТ1 д ?

дТ2 д ?

иэф1

а 1-1

эф 2

2

1 д Т 1 дТ -------1 +----------1

ч дг2 Г дг у

2

д Т2 1 д Т2

----- +----------2

ч дг2 г дг J

дЖ2 д ?

- = Я

дг

г дг

я1 < г < £

£

£

промерзание --------->

оттаивание <---------

диффузия <--------

конденсация <---------

криожидкость

т

©

я,

влажный воздух

Рис. 1. Схематическое изображение области решения: 1 - зона промерзшей изоляции; 2 - зона увлажненной изоляции

Начальные условия:

Т = % Я1 — г — Я2, г = 0; ж = ж0, я — г — я2, г = 0.

Граничные условия:

Т = Тж, г = г > 0;

дТ

Я = а

- я

дг дЖ

дг

= в(Ж - Жос), г = я2, г> 0;

дЖ2

дг

= 0, г = £(г), г > 0;

Т =

пар

241,2 • 1п вппарч 0,6112

17,504 - 1п -

Ж =

ос 1 , т

1 + Л Лв

Влагосодержание влажного воздуха [10]:

а = 0,6112------^-----.

вл.в гвп

1 дТ я дТ2 — Р й£

Яэф1 , Яэф2 - = ЖбплР , ,

дг дг аг

Т = Т2 = ТпЛ, г = £

Объемная доля льда в изоляции, вследствие расширения воды при замерзании, увеличивалась на 9% [9]:

Жл = 1,09Ж.

Значения теплофизических свойств тепловой изоляции рассчитывались как эффективные с учетом объемных долей каждой компоненты. Например, эффективные коэффициенты теплопроводности находились из соотношений [2]:

Яэф1 = ЖЯ + (1 - ж я, Яэф2 = ЖвЯв + (1 - Жв )ЯИ.

Температура замерзания воды вычислялась по формуле [10]:

Тш = (0,075 +1,382 • 10-8 • р^)(1 - 0,0009869 • ртм).

Удельная теплота плавления (кристаллизации) [10]: Оё = 334,11 + 2,261- Т,л.

Удельная теплота парообразования (конденсации) [10]:

Опар = 2500,64 - 2,369 • Тпар,

Температура парообразования (конденсации) водяного пара [10]:

Интенсивность массоотдачи на поверхности изоляции определялась по формуле [11]:

= (Рвп - Р^пов) к .

Давление насыщенного водяного пара на поверхности изоляции [10]:

17 504•Т

рвпт = 0,6112 • ехр ,-------*=-.

241,2 + Т

пов

Коэффициент массоотдачи [12]: а

Р=-

Свл.в 'Рвл,

Теплоёмкость влажного воздуха [10]: с =

влв (1+а) .

Плотность влажного воздуха [10]: 0,00348

рвл.в

-(Р -0,376• Р ).

\-^вл.в ’ -Т^вп/

0,6112

Парциальное давление водяного пара в воздухе [10]: = Р

^вп.парч 100 Ри.

Давление насыщенного водяного пара в воздухе

[10]:

17 504 • Т

о = 0,6112 • ехр —,--------^.

241,2 + ТоС

Объёмная концентрация влаги в воздухе [10]:

Обозначения: Я - граница области расчета; £ -координата границы фазового перехода, м; Т -температура, К; Ж - относительное объемное вла-госодержание изоляции; а - коэффициент температуропроводности, м2/с; / - время, с; г - координата, м; Б - коэффициент диффузии, м2/с; Я - коэффициент теплопроводности, Вт/(м-К); С - теплоемкость, ДжДкг^К); р - плотность, кг/м3; а -коэффициент теплоотдачи, Вт/(м2-К); у - плотность потока массы, кг/(м2-с); 0 - теплота фазового перехода, кДж/кг; в - коэффициент массоотдачи, м/с; р - относительная влажность воздуха, %; к - сопротивление влагообмену, (кШ^м^/кг; р -давление, Па; й - влагосодержание, (кг. влаги)/(кг. сухого воздуха); л - молярная масса, кг/кмоль; р=(рв+рл)/2 - средняя плотность, кг/м3;

Ж=(Жв+Жл)/2 - среднее относительное влагосо-держание изоляции по объёму.

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

Метод решения и исходные данные

Задача (1)—(10) решена методом конечных разностей [13] с использованием метода прогонки

ос

вп.парч

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

В качестве примера рассматривался типичный объект для хранения сжиженного углеводородного газа - аппарат диаметром 2400 мм [6]. В качестве изоляционного материала выбран полистирол ПС-1 толщиной 50 мм [6]. Температура внутренней границы Я1 равна температуре криожидкости в резервуаре Тж=230 К [14]. Температура окружающей среды варьировалась в диапазоне Гос=290...300 К, относительная влажность воздуха р=60...100 %, атмосферное давление раім=101325 Па, начальная объемная влажность изоляции 1^0=1 % [3]. Коэффициент теплоотдачи от воздуха к поверхности изоляции в соответствии с [11] принимался равным а=5,8 Вт/(м2-К), сопротивление влагообмену у поверхности изоляции составляло к=96 (МПа-с-м2)/кг

[11], коэффициент диффузии влаги в полистироле ПС-1 А=2-10-6 м2/ч [4].

В табл. 1 приведены значения теплофизических характеристик [3, 5, 10] веществ и материалов, использовавшихся при проведении численного моделирования.

Таблица 1. Теплофизические свойства материалов и веществ

Характеристика X, Вт/(м-К) С, Дж/(кг-К) р, кг/м3

Вода 0,6 4186 994,04

Лед 2,4 1924 916,8

Полистирол ПС-1 0,0342 1183 100

Результаты численного моделирования

Основные результаты численного моделирования тепломассопереноса в ограждающих конструкциях резервуаров для хранения криожидкостей в условиях миграции влаги к фронту промерзания влажной теплоизоляции приведены в табл. 2 и на рис. 2. В табл. 2 приведены величины тепло-притоков к резервуару с учетом дг и без учета д2 промерзания тепловой изоляции, сравнение этих величин между собой 8, погрешность по балансу энергии 82, значения толщины слоя промерзшей изоляции 8, объемное содержание влаги в атмосферном воздухе Жос и время выхода процессов на стационарный режим при различных значениях температуры окружающей среды Тос и относительной влажности воздуха р.

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

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

Таблица 2. Результаты численного моделирования

Гс, К р, % Я, Вт/м Я, Вт/м ^Х, % 3, мм 3, % 3, % Ъ, ч

290 60 447,6 282,6 1,14 44,7 36,9 0,36 21

80 456,7 1,53 43,7 38,1 0,28 52

100 466,8 1,91 42,7 39,5 0,25 56

295 60 479,7 306,1 1,56 41,6 36,2 0,34 72

80 492,2 2,09 40,2 37,8 0,31 130

100 511,0 2,61 39,0 40,1 0,36 138

300 60 516,7 329,7 2,11 38,5 36,2 0,24 149

80 539,1 2,82 36,8 38,8 0,34 248

100 563,4 3,54 35,2 41,5 0,26 309

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

Результаты расчета потерь холода (табл. 2) свидетельствуют о том, что увлажнение и последующее промерзание слоя тепловой изоляции приводят к росту потерь холода на 8=36,2...41,5 % в рассматриваемом диапазоне изменения Тос и р. Соотношение величин и д2 позволяет утверждать, что при эксплуатации резервуаров для хранения криожидкостей даже при относительно невысоких значениях Тос и р возможно существенное увеличение теплопритоков к рассматриваемым объектам.

Следует также отметить, что в рамках рассматриваемой задачи толщина слоя промерзшей тепловой изоляции имеет максимальное значение 8=44,7 мм (табл. 2), что соответствует практически полному промерзанию теплоизоляционного покрытия рассматриваемого объекта. Анализ результатов моделирования свидетельствует о том, что толщина слоя промерзшей изоляции заметно из-

5

Г, час

Рис. 2. Изменение теплопритоков к резервуару во времени: 1 - 290; 2 - 295; 3 - 300 К.

меняется (в пределах 20 %) в зависимости от значений Тос и р. Эти обстоятельства обуславливают необходимость учета нестационарности процессов переноса и возможного промерзания слоя тепловой изоляции при проектировании и анализе тепловых режимов работы резервуаров для хранения криожидкостей.

На рис. 2 приведены величины теплопритоков к резервуару в зависимости от времени при температуре окружающей среды Тос=290...300 К, значении относительной влажности воздуха р=100 %.

Анализ нестационарности процессов переноса в холодильной изоляции из полистирола ПС-1 свидетельствует о том, что выход на стационарный режим происходит через ?от=21—309 часов (табл. 2) в зависимости от Тос и р.

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

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

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

Выводы

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

Установлены масштабы увеличения теплопритоков (до 41,5 %) в рассматриваемой системе вследствие увлажнения, промерзания и наличия фазовых переходов в структуре и на поверхности слоя тепловой изоляции.

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

Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 12-08-00201-а и гранта Президента РФ № МК-1652.2013.8.

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

1. Энергетическая стратегия России на период до 2030 года // Официальный сайт Министерства энергетики РФ. 2013. URL: http://minenergo.gov.ru/activity/energostrategy/ (дата обращения: 22.01.2013).

2. Половников В.Ю., Хабибулин А.М. Численное моделирование теплопереноса в ограждающих конструкциях резервуаров для хранения криожидкостей с учетом промерзания тепловой изоляции // Известия Томского политехнического университета. - 2012. - Т. 320. - № 4. - C. 33-36.

3. Каганер М.Г. Тепловая изоляция в технике низких температур. - М.: Машиностроение, 1966. - 275 с.

4. Петров-Денисов В.Г., Масленников Л.А. Процессы тепло-и влагообмена в промышленной изоляции. - М.: Энергоато-миздат, 1983. - 193 с.

5. СНиП 41-03-2003. Тепловая изоляция оборудования и трубопроводов. - М.: Изд-во стандартов, 2004. - 28 с.

6. Архаров А.М. Криогенные системы: основы проектирования аппаратов, установок и систем. - М.: Машиностроение, 1999. - 720 с.

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

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

9. Курылев Е.С., Оносовский В.В., Румянцев Ю.Д. Холодильные установки. - СПб.: Политехника, 2004. - 576 с.

10. Бурцев С.И., Цветков Ю.Н. Влажный воздух. Состав и свойства. - СПб.: СПбГАХПТ, 1998. - 146 с.

11. Фокин К.Ф. Строительная теплотехника ограждающих частей зданий. - М.: Стройиздат, 1973. - 287 с.

12. Цветков Ф.Ф. Тепломассообмен. - М.: Изд-во МЭИ, 2005. - 549 с.

13. Самарский А.А., Гулин А.Н. Численные методы математической физики. - М.: Научный мир, 2000. - 316 с.

14. Яковлев Е.И. Тепловые режимы хранилищ сжиженных газов. - СПб.: Недра, 1992. - 182 с.

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

УДК 536.2:532/533; 532.516

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАЗРУШЕНИЯ ПОЛОСТИ ОТКРЫТОГО ТИПА ПРИ ВЗАИМОДЕЙСТВИИ С ВЯЗКОЙ ЖИДКОСТЬЮ В УСЛОВИЯХ ОХЛАЖДЕНИЯ ВНЕШНЕГО КОНТУРА

А.В. Крайнов

Томский политехнический университет E-mail: Lux_veritatis@mail.ru

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

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

Сопряжённый теплоперенос, вязкая несжимаемая жидкость, гидродинамика, контур кратера, разрушение полости.

Key words:

Conjugate heat transfer, viscous incompressible fluid, hydrodynamics, crater contour, cavity damage.

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

пени опасности по мере развития опаснейших производств и усложнения технологических циклов в таких отраслях промышленности, как атомная, химическая и других [6, 7].

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

Изучение процесса проводилось с использованием математической модели на основе системы уравнений Навье-Стокса в переменных вихрь (ю)-функция тока (^), уравнения энергии, уравнения теплопроводности для материала полости с соответствующими краевыми условиями