CC BY
56
УДК 536.24
ОСОБЕННОСТИ ВЛИЯНИЯ ДАВЛЕНИЯ НА ИНТЕНСИВНОСТЬ ТЕПЛООТДАЧИ ПРИ КИПЕНИИ ПРОПАНА НА ВЫСОКОТЕПЛОПРОВОДНОЙ ПОРИСТОЙ ПОВЕРХНОСТИ
Л.Л. ВАСИЛЬЕВ, А.С. ЖУРАВЛЕВ
АНК Институт тепло- и массообмена им. А. В. Лыкова НАН Беларуси, Республика Беларусь
А.В. ОВСЯННИК, М.Н. НОВИКОВ
Учреждение образования «Гомельский государственный технический университет имени П. О. Сухого»,
Республика Беларусь
Проведено экспериментальное исследование влияния давления на теплообмен при кипении пропана на образцах с медной спеченной капиллярно-пористой структурой в диапазоне тепловых потоков ^=100^64000 Вт/м и давлений насыщения ^н=0,48^1,08 МПа (Тн=0^30°С). Выполнен анализ влияния давления на теплообмен при кипении пропана в условиях большого объема в зависимости от плотности теплового потока и характеристик пористой структуры.
Введение
В настоящее время активно ведутся поиски экологически безопасных и нетоксичных теплоносителей способных обеспечивать эффективный теплообмен. Поэтому несомненный интерес для использования представляют углеводороды (пропан, пропилен, бутан, пропан-бутановые смеси), обладающие рядом уникальных свойств и являющиеся альтернативой озоноразрушающим фреонам.
Пропан является диэлектрической жидкостью, химически не активной для большинства конструкционных материалов. Благодаря этим свойствам в сочетании с высокими теплофизическими характеристиками пропан может широко применяться в испарительных теплообменниках погружного типа для охлаждения элементов микроэлектроники, а также в качестве рабочей жидкости в тепловых трубах. Определение особенностей влияния давления на интенсивность теплоотдачи при кипении пропана в большом объеме на пористых покрытиях дает возможность выбора оптимальных эксплуатационных режимов работы теплообменного оборудования.
1. Экспериментальное оборудование
Для исследований теплообмена при кипении сжиженного пропана в АНК ИТМО НАНБ был создан специальный экспериментальный стенд, схема которого и методика проведения экспериментов подробно описаны в [1]. Экспериментальные образцы представляли собой медные горизонтальные трубы 020x1.75 мм и длиной 100 мм, на наружную поверхность которых напекалось пористое покрытие. Пористость образцов была постоянной и составляла 50^55 %. Характеристики исследованных образцов представлены в таблице 1.
Таблица 1
Характеристики исследованных образцов
№ Толщина структуры Зкс, мм Диаметр частиц ёч, мкм Средний гидравлический размер пор й0, мкм
1 0,2 3 6 0 4 18,37
2 0,3
3 0,4
4 0,5
5 0,3 63...100 24,48
6 0,4
7 0,5
8 0,6
9 0,8
10 0,4 100...+160 35,4
11 0,5
2. Результаты исследования и их обсуждение
Экспериментальное исследование давления насыщения на теплообмен при кипении пропана на образцах со спеченной капиллярно-пористой структурой из медного порошка марки ПМС-Н было проведено в диапазоне тепловых потоков д=100^64000 Вт/м2 и давлений насыщения р„=0,48^1,08 МПа (Т„= 0^30°С).
Известно, что величина коэффициента теплоотдачи при кипении жидкости зависит от давления, причем степень влияния давления по-разному проявляется в зависимости от вида теплоотдающей поверхности и теплофизических свойств жидкости. Имеющиеся в литературе опытные данные по кипению различных жидкостей на поверхностях с капиллярно-пористой структурой не дают возможности однозначно определить характер степени влияния давления на теплообмен: в одних работах [24] отмечено увеличение степени влияния давления на коэффициент теплоотдачи, результаты других [5-9] говорят о примерном равенстве или снижении этого влияния по сравнению с гладкой поверхностью. Такое качественное различие степени влияния давления можно объяснить: различием областей рабочих давлений; отличием физических свойств теплоносителей; изменением в значительном диапазоне теплопроводности скелета структуры; многообразием используемых типов покрытий и их геометрическими характеристиками. Поэтому, на наш взгляд, более корректным будет представление опытных данных в виде зависимости а=/(р/ркр). На основании проведенных экспериментальных исследований был определен показатель степени в соотношении а=С (р/ркр)к и выявлено влияние на него различных факторов, определяющих интенсивность теплоотдачи при кипении.
Одним из параметров, оказывающих существенное влияние на значения к, является плотность теплового потока. На рис. 1 показана зависимость к = _Д^) для различных фракций. Как видно из рисунка, с повышением теплового потока степень влияния давления снижается. По мере возрастания плотности теплового потока (или перегрева теплоотдающей поверхности) кипение начинается во все более мелких порах, ранее заполненных жидкостью. Очевидно, что для каждого давления существует граничный перегрев поверхности ЛТгр, начиная с которого подавляющее количество пор капиллярно-пористой структуры становятся активными. При ЛТ>ЛТгр ак-
тивизируются поры очень малых размеров, количество которых невелико. Из-за сравнительно небольшого прироста количества центров парообразования увеличение интенсивности теплоотдачи незначительно. Поскольку с увеличением давления значения ЛТгр снижаются (уменьшается критический радиус паровых зародышей), то при малых тепловых потоках влияние давления будет проявляться сильнее по сравнению с большими значениями ц. При ЛТ>ЛТгр достигается примерное равенство активных и возможных центров парообразования, что приводит к слабому или вообще отсутствию влияния давления на теплообмен. Действительно, для большинства исследованных образцов в диапазоне тепловых потоков ц>15^30 кВт/м2 (в зависимости от геометрических характеристик покрытия) наблюдалась автомодельность коэффициентов теплоотдачи от давления (к^0), что соответствует II режиму развитого кипения [1].
а) б)
в)
Рис. 1. Зависимость показателя степени к в соотношении а=С-(р/ркр)к от плотности теплового потока q при кипении пропана на поверхности со спеченной капиллярнопористой структурой для различных фракций: 1- гладкая поверхность; а) ёч=40+ 63 мкм, 2^5 - образцы № 1-ь4, соответственно (табл. 1); б) й?ч=63^100 мкм, 2^6 - образцы № 5-ь9, соответственно (табл. 1); в) й?ч=100^160 мкм, 2^3 - образцы № 10^11, соответственно (табл. 1)
Средний гидравлический размер пор й0, учитывающий влияние таких характеристик пористого покрытия, как пористость и средний размер частиц, составляющих
скелет структуры, оказывал незначительное влияние на значения к (рис. 2). Увеличение размера частиц приводит к большему размаху дифференциальной функции распределения пор по размерам и, как следствие, к увеличению значений ё0. Большим значениям ё0 соответствуют большие значения АТгр при одинаковой толщине покрытия (при ^о=18,37 мкм АТгр(^=0,5 мм)=0,7 К, а при ё0=48,8 мкм АТгр(^=0,5 мм)=1,8 К; Т„=20 °С). В то же время, число пор, приходящееся на единицу поверхности (оцененное по формуле 5&4П/(ж- ё0)), для образца с ё0=48,8 мкм почти на порядок ниже (1,9-109 м-2 и 2,67-108 м-2, соответственно). Разница в количестве активных центров парообразования при изменении давления для образца с ё0=18,37 мкм, следовательно, должна быть больше. Однако влияние этих факторов на показатель степени к взаимно компенсируется. Для образцов, состоящих из частиц мелких фракций, значительное количество пор активно уже при малых тепловых потоках («взрывное вскипание»), и темп активации центров парообразования становится соизмеримым с приростом активных центров на образцах из более крупных фракций. При увеличении плотности теплового потока степень этого влияния с изменением ё0 становится еще менее заметной, хотя следует отметить тенденцию к некоторому снижению к с ростом ё0 при небольших тепловых потоках.
0.7 - к
0.6 - & _ ~ ~ о
0.5 - о- . ^ —о .
0.4 - О ' О
0.3 - "О -
0.2 - А д "
0.1 -
0.0 1
15 20 25
О
о
Л
>
V
дск = 0,3 мм 8ск = 0,4 мм
8ск = 0,5 мм
.2
- 5=3 кВт/м
- 5=6 кВт/м2
- ^=10 кВт/м2
- ^=20 кВт/м2
- 5=64 кВт/м2
-I — \"-------
30 _ 35
ё0, мкм
—г-
40
45
50
Рис. 2. Зависимость показателя степени к в соотношении а=С- (р/ркр) от среднего гидравлического размера пор ё0 для фиксированных толщин пористого покрытия при различных тепловых потоках
Существенное влияние на показатель степени к в соотношении а=С-{р/ркр)к оказывает толщина пористого покрытия 5сК. Полученная экспериментально зависимость к=/(дск) представлена на рис. 3. Как видно из рисунка, увеличение толщины пористого покрытия при небольших тепловых потоках приводит к монотонному снижению показателя степени к. Это связано с тем, что количество парогенерирующих пор не изменяется, а гидравлическое сопротивление с ростом толщины структуры выходящему пару возрастает. Известно, что с ростом теплового потока толщина перегретого пограничного слоя снижается и в определенном диапазоне его значений (до кри-
зиса кипения) не зависит от давления. Поэтому для достаточно больших тепловых потоков (режим II) можно предположить, что увеличение толщины структуры не приводит к дополнительному приросту количества центров парообразования. Вследствие этого наблюдается автомодельность значений к от толщины капиллярнопористой структуры.
— - <1ч= 63-100 мкм - - «ч= 100-160 мкм О - <3=3 кВт/м2 О - д=6 кВт/м2 Л - д=10 кВт/м2 > - д=20 кВт/м'
< - 3=40 кВт/м2 V - з=64 кВт/м2
0.5
8,
0.6
мм
к
8
Рис. 3. Зависимость показателя степени к в соотношении а=С-(р/ркр)к от толщины пористого покрытия 5ск при различных плотностях теплового потока для различных фракций: а) й?ч=40-63 мкм; б) й?ч=63-100 мкм и й?ч=100-160 мкм
Как отмечено выше, степень влияния давления определяется темпом активации пор капиллярной структуры. Для структур, имеющих широкий диапазон размеров пор, это влияние проявляется сильнее, поскольку с изменением давления насыщения изменяется и количество активных центров парообразования. Этим можно объяснить более сильное влияние давления по сравнению с гладкой поверхностью на металловолокнистых и напыленных покрытиях. Для высокотеплопроводных медных структур, изготовленных спеканием из частиц одной фракции, количество пор, приходящихся на максимум дифференциальной функции распределения, составляет приблизительно 80 % и лежит в довольно узком интервале значений [1, 7]. При приближении критического диаметра пузырей к эффективному размеру пор наблюдается эффект «взрывного вскипания». В этом случае большинство пор капиллярной структуры являются активными. При таких условиях, в зависимости от фракции, в определенном интервале давлений насыщения его влияние должно ослабевать, что и подтверждают экспериментальные данные нашего исследования.
Выводы
В результате экспериментального исследования были выявлены особенности влияния давления на коэффициент теплоотдачи в зависимости от плотности теплового потока и параметров пористого покрытия для медных спеченных капиллярнопористых структур:
1. При увеличении плотности теплового потока значение показателя степени к в соотношении а=С-(р/ркр)к уменьшается и достигает практически нулевого значения при переходе ко второму режиму развитого кипения.
2. С ростом толщины пористого покрытия (в исследованном диапазоне значений 5ск) степень влияния давления также снижается.
3. Изменение среднего гидравлического размера пор d0 в исследованном диапазоне оказывает незначительное влияние на показатель степени к в зависимости a=C-(p/pKp)k.
4. Степень влияния давления при кипении на пористой поверхности определяется темпом активации центров парообразования, зависящим от типа покрытия и его параметров.
Литература
1. Новиков М.Н. Теплообмен при кипении пропана на горизонтальных трубах со спеченным пористым покрытием //Вестник ГГТУ. - 2000. № 2. - С. 58-66.
2. Орлов В.К., Савельев В.Н. Интенсификация теплообмена при кипении криогенных жидкостей при давлениях ниже атмосферного //Теплоэнергетика. - 1980. -№ 4. - С. 62-64.
3. Кузма-Кичта Ю.А., Комендантов А. С., Бурдунин М.Н. и др. Исследование интенсификации теплосъема в парогенерирующих каналах с пористым покрытием //Теплоэнергетика. - 1991. - № 5. - С. 42-47.
4. Hongji Z., Aihua J. Experimentation of boiling heat transfer from thin powder porous layers at atmospheric and super- atmospheric pressures. //8-th International Heat Pipe Conference. - Bejing, China, 1992. P. A-11, P. 1-6.
5. Зарипов В.К., Семена М.Г., Шаповал А.А., Левретов А.И. Интенсивность теплообмена при кипении на поверхности с пористыми покрытиями в условиях свободного движения //ИФЖ. - 1989. - Т. 57. - № 2. - С. 181-186.
6. Вершинин С.В., Майданик Ю.Ф. Исследование теплообмена при кипении ацетона на мелкопористой капиллярной структуре //Свердловск, 1987. Деп. в ВИНИТИ, № 666-В87, 20 с.
7. Васильев Л.Л., Конев С.В., Штульц П., Хорват Л. Экспериментальное исследование теплообмена при кипении жидкости в высокотеплопроводных капиллярных структурах //ИФЖ. - 1982. - Т. 42. - № 6. - С. 893-898.
8. Позняк В.Е., Савельев В.Н. Опыт применения капиллярно-пористых покрытий в криогенных системах и установках //Теплоэнергетика. - 1990. - № 12. - С. 9-12.
9. Шаповал А.А., Зарипов В.К., Семена М.Г. К расчетам интенсивности теплообмена при кипении на поверхности с пористыми покрытиями //Изв. АН СССР. Энергетика и транспорт. - 1989. - № 3. - С. 63-68.
Получено 22.05.2001 г.
