Научная статья на тему 'Методы расчета градиента давления двухфазного потока при течении в малых каналах'

Методы расчета градиента давления двухфазного потока при течении в малых каналах Текст научной статьи по специальности «Физика»

CC BY
961
133
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
ТЕПЛООБМЕННИК / ДВУХФАЗНЫЙ ПОТОК / ХЛАДАГЕНТ / ГРАДИЕНТ ДАВЛЕНИЯ / HEAT EXCHANGER / DIPHASIC STREAM / COOLANT / PRESSURE GRADIENT

Аннотация научной статьи по физике, автор научной работы — Ховалыг Д., Бараненко Александр Владимирович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

n this article the question will be dealt with two-phase flow behavior of matter in small channels and an overview of methods for calculating the pressure gradient. A comparison of experimental data with computational methods.

Текст научной работы на тему «Методы расчета градиента давления двухфазного потока при течении в малых каналах»

УД К 621.56; 536.2

Методы расчета градиента давления двухфазного потока при течении в малых каналах

Д. ХОВАЛЫГ, д-р техн. наук К. В. БАРАНЕНКО

Санкт-Петербургский национальный исследовательский университет ИТМО Институт холода и биотехнологий 191002, Санкт-Петербург, ул. Ломоносова, 9

In this article the question will be dealt with two-phase flow behavior of matter in small channels and an overview of methods for calculating the pressure gradient. A comparison of experimental data with computational methods.

Key words: heat exchanger, diphasic stream, coolant, pressure gradient.

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

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

Серьезно изучать теплообменники с малыми каналами стали в начале 90-х годов XX в. Определение границы разделения на макро- и микротеплообменники является непростой задачей. Кандликар С. Г. в своей работе [1] в 2002 г. предложил общепринятую ныне классификацию каналов в зависимости от их гидравлического диаметра Dh:

1) традиционные каналы, > 3 мм;

2) миниканалы, 200 мк < Dh < 3 мм;

3) микроканалы, 10 мк < Dh < 200 мк.

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

Теоретическое описание поведения двухфазного потока вещества в малых каналах компактных испарителей является достаточно сложной задачей в связи с комплексным поведением кипящего потока жидкости. Вынужденное кипение в мини- и микроканалах ставит три основных вопроса перед исследователями и инженерами [1]:

1. Каким образом малый размер канала влияет на парообразование и на вынужденное течение двухфазного потока внутри его?

2. Каков характер коэффициента теплоотдачи и градиента давления?

3. Какие различия существуют в поведении потока при течении в одиночном канале и в нескольких параллельных каналах?

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

Структура двухфазного потока

Размер канала значительно влияет на характер парообразования во время кипения жидкости в потоке. На рис. 1 схематически показаны стадии роста парового пузыря внутри малого канала и характер его поведения [4]. В начальный момент при подводе теплоты через стенки канала жидкость начинает перегреваться и в отдельных точках поверхности (в центрах парообразования) образовываются пузырьки. При повышении температурного напора давление внутри пузыря увеличивается и пузырь начинает расти, занимая все сечение малого канала (см. рис. 1, б). Если движущая сила парооборазования достаточно велика для преодоления динамического давления потока и (или) давления на входе в канал, пузырь

а Г . ' о ►

Ь

" L-)

в — С 3

С Л У

д — с —► 3

— Г

Рис. 1. Схематическое изображение роста пузыря и характер его движения в малом канале

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

К настоящему времени разработано большое количество карт режимов двухфазного потока жидкости в малых каналах |4, 5]. Первое подробное изучение структуры потока в малых каналах было проведено в 1964 г. исследователями Суо М. и Гриффис П. на основе наблюдений за потоком жидкого гептана и воды, а также газообразного гелия и азота в каналах круглого сечения диаметрами 1,03 и 1,6 мм [6]. Ими была составлена карта режимов двухфазного потока (рис. 2) в зависимости от свойств вещества и расходных характеристик: динамической вязкости жидкости дж, поверхностного натяжения сг, скорости жидкости иж, объемного расхода жидкости (,}ж и газа QT.

Корнвелл К. и Кью П. А. в 1992 г. выявили четыре режима двухфазного потока хладагентов R113 и R141b в каналах прямоугольного сечения размерами 1,2x0,9 мм и 3,3x1,1 мм: отдельные пузыри, стесненные пузыри, снарядно-кольцевой режим и режим, в котором наблюдается частичное высыхание стенок.

На рис. 3 изображена разработанная Янгом и Шиих [7] карта режимов хладагента R134a и смеси вода—воздух в трубке диаметром 2 мм. Они определили шесть режимов: пузырьковый (bubble), снарядный (slug), пробковый (plug), расслоенный волновой (wavy stratified), дисперсный (dispersed) и кольцевой (annular). Необходимо отме-

&

Рис. 2. Карта режимов двухфазного потока в малых каналах, разработанная Суо М. и Гриффис П.:

Чж/^г = 50; Яеж№еж = 0,75- 1(Г5

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

Существуют и более детальные разбивки режимов двухфазного потока. Например, Колеман и Гаримеллавсерии работ (1999, 2000, 2003 гг.) выделили четыре основных режима двухфазного потока хладагента Я 134а: дисперсный, прерывный, волновой, кольцевой — и разделили их суммарно на 16 режимов [5].

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

На рис. 4 показаны наблюдаемые в работе [8] режимы двухфазного потока хладагента Я245Га в канале диаметром 0,509 мм при температуре насыщения 35 °С, температуре переохлаждения 6 °С и расходе С = 500 кг/(м2 с).

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

^пр. Г» м/с

Рис. 3. Карта режимов двухфазного потока, разработанная Янгом и Шиих: иПр. ж — приведенная скорость жидкости; ипр. г — приведенная скорость газа

Пузырьковый режим

б Переходный режим

от пузырькового к снарядному

С ГТГ

! л

в Снарядный режим

ь от )

г Переходный режим

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

д Полу кольце вой режим

е Кольцевой режим

Паросодержание х = 0,038

00

1

х— 0,04

х= 0,043

х= 0,076

Рис. 4. Схематическое изображение и фотографии изменения характера двухфазного потока Я245/а в канале диаметром 0,509 мм

а б в

¿/Кан = 2,0 мм ¿кан = 0,8 мм ¿кан = 0,5 мм

Рис. 5. Иллюстрация влияния подъемной силы в малых каналах [9]

Таким образом, структура потока практически не зависит от ориентации канала, разность истинных скоростей фаз (скорость скольжения потока) уменьшается, это сдерживает появление касательного напряжения на границе жидкость—газ и подавляет образование волн на границе раздела фаз [8].

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

ралелльных каналах (рис. 6). Наряду с этим существуют значительные флуктуации давления в связи с локальным реверсивным движением жидкости.

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

d P d 2

d P\ dl)Tp +

d P d г

+

d P d г

(1)

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

(2)

где Ах — изменение массового расходного паросодержа-ния;

рж — плотность жидкости; рг — плотность газа.

В общем случае уравнение (1) имеет следующий вид [4]:

Рис. 6. Структура двухфазного потока хладагента Я 134а в шести параллельных каналах:

О/, = 1,54 мм; Т= 10 °С;х = 17%;й = 500кг/(м2 с) [10]

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

Теоретические модели описания двухфазного потока в малых каналах приведены в работах [11 — 13].

Методы описания градиента давления

Существует три фактора, которые влияют на градиент давления при движении двухфазного потока в канале: сила вязкого трения, статический напор и инерция потока [4]:

Q — угол наклона канала;

ос — истинное объемное паросодержание;

х — массовое расходное паросодержание.

Составляющая градиента давления за счет вязкого fdP\

трения —— определяется на основе различных

V d г / тр

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

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

Гомогенная модель

Падение давления, обусловленное трением, определяется из расчета потерь давления однофазного потока через коэффициент гидравлического сопротивления двухфазного потока ftp, рассчитываемый на основе свойств парожидкостной смеси. Для адиабатного потока с постоянным паросодержанием (ж = const) формула расчета падения давления определяется отношением [5]

d Р d z

тр

2 ftPLG2 Dptp

(4)

где ftp — коэффициент гидравлического сопротивления двухфазного потока (фактор Фаннинга);

L — длина канала;

G — массовый расход вещества;

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Ptp — плотность потока;

D — внутренний диаметр канала.

Для диабатного потока с переменным паросодержанием (х = var) падение давления рассчитывают, учитывая изменения паросодержания:

dP d z

тр

2 ftpLG2 D ptp

1 +

Ax

2 p*

Pv

(5)

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

- ( = ("¡¡т) + [(1-<*)/»*+ <*Рг] 5 sin fi+

d

+ dI

G2x2 G2( l-xW

+

_ржа рж(1-а) где g — ускорение свободного падения;

(3)

ftp

WRet-n,

(6)

где N = 0,079 и п = 1/4 для турбулентного течения; N = 16 и п = 1 для ламинарного течения. Средняя плотность потока определяется как

х 1-а;

Ptp = I-------1-------

\Рг Рж

Определение средней динамической вязкости потока основано на том, что вязкость паровой фазы определяется при паросодержании потока, близком к единице, а вязкость жидкой фазы определяется при паросодержании потока, близком к нулю. Соотношения для расчета вязкости в литературе [14] представлены: МакАдамсом, Вудсом и Брайаном (1942 г.); Девидсоном и коллегами (1943 г.); Сичитти и коллегами (1960 г.); Овенсом (1961 г.); Даклером, Виксом и Клевелендом (1964 г.); Битти и Уаллей (1982 г.); Гарсиа и коллегами (2003 г.); Авад и Музычка (2008 г.).

Далее приведены наиболее часто используемые соотношения:

МакАдамс, Вудс и Брайан

, X 1-х

М<р = I '—I-----------

V Мг Дж

Сичитти и др.

Ц'Ьр хрг + (1 х)рж, Даклер, Викс и Клевеленд

I Р1г , /1 \ Мж

Мір = Рір ( X------1- (1-х) —

Р г Рж

(8)

(9)

(10)

В общем виде расчет градиента давления в случае принятия гомогенной модели представления двухфазного потока имеет следующий вид [4]:

сі Р Л г

Ф2 +

ЖТ '

2/жТ

Ш2уж 2 ¿х де іпГ2

-------------+ +-----------------------

О рір аг иж + хуж

1 + Є2х

сі гіг

Тр

(П)

Ф =

ЖТ

ft.ll

1+ [ —-1 )х Рт

(12)

Модель разделенных потоков

Модель разделенных потоков предполагает, что каждая фаза движется со своей средней скоростью [15]. Модель непрерывно дополняется с момента ее разработки Локхартом и Мартинелли в 1949 г. Основной целью данной модели является нахождение эмпирической взаимосвязи множителя трения двухфазного потока Ф2 с независимыми переменными потока. Множитель Ф2 и соотношение падения давления для однофазного потока применяются для определения потерь двухфазного потока [4, 5]:

Ф2 = 1 + +

1

X X2’

(13)

где С — константа, зависящая от режима каждой фазы (табл. 1); X — параметр Мартинелли,

АР

АР

дг

1/2

Падение давления, обусловленное силами вязкого трения, определяется из соотношения

Э_-(ё' »•

(14)

где

СІР СІ 2

определяется из предположения, что весь

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

йР с! г

2/ж ьс2

Орж

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

ж — индекс, обозначающий свойства вещества, если бы только жидкость текла в канале с соответствующим расходом С(1 — х)5;

г — индекс, обозначающий свойства вещества, если бы в канале существовал только поток газа с расходом СЗх;

жг — индекс, обозначающий разность свойств газа и жидкости.

Множитель трения двухфазного потока Ф2КТ определяется из соотношения

В общем случае расчет градиента давления при использовании метода разделенных потоков определяется соотношением [4]

6Р ' її!

Ф2

2/жтС2

В Рьр

+ [(1 - а)рж + арг]ззіпП+

+с-

СІ х

2хуг 2(1— х)у»

+ -

с! а

1

(1 - х)2уж УГХ2

+

(1-а)2

(15)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

где

Л = 1 + в1

X

а

с! уг сПР

+

сіа /(1 — х)2уж угхі

йР \ (1-а)2

Таблица 1

Значение константы С в зависимости от режимов жидкости и газа [5]

Режим течения жидкости Режим течения газа Константа С

Турбулентный Турбулентный 20

Ламинарный Турбулентный 12

Турбулентный Ламинарный 10

Ламинарный Ламинарный 5

Наиболее часто используемые методы расчета градиента давления для традиционных больших каналов приведены в работе [16], в том числе расчетное соотношение Мюллера—Стейнхагена и Хека:

/ тр

I

йР

йг

где

ар

СІ 2

Р = — + 2х

6Р\

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

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

Таблица 2

Соотношения для расчета падения давления в малых каналах [16, 17]

Авторы (год) Расчетные соотношения Примечания

Мишима и Хибики (1996 г.) На основе метода Локхарта—Мартинелли константа С рассчитывается из соотношения С = 21(1 — е-0,319°). Пароводяная смесь в вертикальной трубке круглого сечения диаметром 1—4 мм, адиабатные условия

Тран и другие (2000 г.) Ш -(£) \ / тр \ / жт где Ф^т = 1 + (4,ЗУ2 - 1)[Сож°’875(1 - ж)0'875 + ж1’75], с ( а У’5 г2 (ар] '[6Р] Ш34а, Я12 и Ш 13 в каналах круглого сечения диаметрами 2,46 и 2,92 мм и в канале прямоугольного сечения 4,06х 1,7 мм

Vт>2д{рж — рг)/ ’

Ли и Л и (2001 г.) На основе метода Локхарта—Мартинелли константа С рассчитывается из соотношения С = А\чфгК&т, где Л = ; ф = ^ ржа и а (значения чисел А, д, г, 5 приведены в табл. 3). Двухфазная смесь воды с воздухом в горизонтальных каналах сечением (0,4...4) 20 мм

Жанг и Вебб (2001 г.) т -т \ / тр \ / жт где Фжт = (1 - ж)2 + 2,87ж2 ( - Р Л + \Крит / / р N-1,64 + 1,68ж0,8(1 ж)0'25 ( р ) \ -* крит / И. 134а, Я22 и Я404а в каналах круглого сечения диаметром 2,13; 3,25 и 6,25 мм, адиабатные условия

Ли и Мудавар (2005 г.) На основе метода Локхарта-Мартинелли константа С рассчитывается из следующих соотношений: — для ламинарных потоков жидкости и газа Слл = 2,1(Же^!47\¥е^; — для ламинарного потока жидкости и турбулентного газа Слт = 1,45Кеж25\Уеж23, где \¥ежт = РжСГ Ш34а в микроканалах прямоугольного сечения 231 х713 мк

Жанг и Мишима (2006 г.) На основе метода Локхарта-Мартинелли константа С рассчитывается из соотношения: С - 21(1 е-°'358/и1), где Ьа - ~ Р^°'Ъ

Сан и Мишима (2009 г.) На основе метода Локхарта-Мартинелли с ( Не, 1 -ж\ ф2 1 . \Кеж’ X ) 1 Соотношение разработано на основе статистических данных [17]

Хп X2’ /п \ 0,4 / 1 Ч 0,5 ГДеС = 1’79(кеж) ( ж ) ;П = 1Д9

Таблица 3

Значения параметров для расчетного соотношения Ли и Ли [18]

Режим потока жидкость/газ А Я Я 5 X Кежт

Ламинарное/ламинарное 6,833-Ю-8 -1,317 0,719 0,557 0,776-14,176 175-1480

Ламинарное/турбулентное 6,185-10 2 0 0 0,726 0,303-1,426 293-1506

Турбулентное/ламинарное 3,627 0 0 0,174 3,276-79,415 2606-17642

Турбулентное/турбулентное 0,408 0 0 0,451 1,309-14,781 2675-17757

Состоятельность расчетных методик при сравнении с экспериментальными данными

В работе [19] выполнен анализ определения потерь давления за счет сил вязкого трения для каналов различных диаметров, включая малые, по 12 расчетным методикам на основании более 900 экспериментальных данных для восьми веществ, полученных различными авторами. Рассматривались двухфазная смесь воды с воздухом и такие хладагенты, как Я12, Ш 13, Л 134а, Я410а и другие; массовый расход варьировался от 23 до 6000 кг/(м2 с); значение паросодержания достигало единицы; эксперименты проводились как при адиабатных, так и при диабатных условиях. Критериями сравнения являлись средняя абсолютная погрешность и попадание экспериментальных данных в ±30 %-й диапазон от расчетных значений. Наилучшее совпадение с экспериментом показала расчетная методика Мюллера-Стейнхагена и Хека [16]. Гомогенная модель с расчетом вязкости по соотношению Си-читти (9) [5, 141 и расчетный метод Мишима—Хибики (см. табл. 2) на основе модели разделенных потоков Локхарта—Мартинелли показали второй результат.

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

В работе [14] при сравнении 3908 экспериментальных данных по падению давления двухфазного потока в каналах диаметрами от 0,517 до 31,7 мм с расчетными величинами по 24 теоретико-эмпирическим корреляциям выявлено, что наилучшее предсказание потерь давления при кольцевом режиме в малых каналах показывают расчетные соотношения Ломбарди [20], при которых 87,5 % экспериментальных данных находятся в диапазоне ±30 % от расчетных; Мюллера—Стейнха-гена и Хека — 67,2 % [16], а также гомогенная модель с расчетом динамической вязкости по Сичитти — 71,8 % [5, 141.

В работе 117| были обобщены 2092 экспериментальных данных по потерям давления при двухфазном течении таких веществ, как К123, Я134а, Я22, К236еа, Я245Га, Я404а, Я407С, Я410а, Я507, С02, вода и воздух, из 18 публикаций. Гидравлические диаметры каналов варьируются от 0,506 до 12 мм; числа Рейнольдса, рассчитанные для жидкости, варьируются от 10 до 37000, для газа — от 3 до 4- 105. Сравнение экспериментальных данных с результатами расчета потерь давления за счет сил вязкого трения при двухфазном течении веществ по 12 теоретическим методам показало, что во всем диапазоне практических данных расчетное соотношение Сана и Мишима (см. табл. 2) определяет наилучшую сходимость. Это объясняется представлением множителя трения двухфазного потока Ф2 зависимым от режимов потока (от чисел Рейнольдса, рассчитанных отдельно по свойствам жидкости и газа). В области турбулентного течения расчетное соотношение Мюллера—Стейнхагена и Хека оказывается наилучшим, наряду с корреляцией Сан—Мишима. Однако для ламинарного режима (Яеж и Яег меньше 2000), характерного для малых каналов, расчетные соотношения Мишима—Хибики и Жанг—Мишима показывают наибольшую сходимость с экспериментальными данными (см. табл. 2).

В работе [21] приведены результаты исследования по определению потерь давления двухфазного потока ИЛ34а в прямоугольных каналах гидравлическим диаметром от 69,5 до 304,7 мк при массовом расходе 101,7—793,8 кг/(м2 с), паросодержании х = 0,02 ч- 0,98 и температурах насыщения 23,2—28,9 °С. Установлено, что из 266 значений наилучшее совпадение с экспериментом имеет гомогенная модель с расчетом свойств хладагента поДаклеру[5, 14]; среднее отклонение составляет 16,5%.

В работе [22] проводилось сравнение экспериментальных данных с расчетными значениями по падению давления при течении диоксида углерода в круглой трубке диаметром 0,529 мм при массовом расходе 200-1400 кг/(м2-с), паросодержании до х = 1 и температурах насыщения —10, —5, 0, 5 °С. Сравнения показали, что для всего диапазона температур наилучшее предсказание дают данные, полученные по расчетному соотношению Мюллера-Стейнхагена и Хека — 75 % экспериментальных данных находятся в диапазоне ±20 % от расчетных [16]. Наряду с этим сделано наблюдение, что при температуре —10 °С двухфазный поток лучше описывается моделью разделенных потоков, а при температуре 5 °С — гомогенной моделью. Таким образом, подтверждается очевидный вывод, что при более высоких температурах структура двухфазного потока более однородна из-за уменьшения отличий в свойствах жидкой и газовой фаз вещества (в частности, плотности), а при низких температурах двухфазный поток больше предрасположен к расслоению.

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

В общем случае сравнение экспериментальных данных с данными, полученными расчетными методами, показывает следующее. Для традиционных больших каналов расчетные методы в среднем дают завышенные результаты, с увеличением температуры отклонение увеличивается. Гомогенная модель, наряду с расчетными методами, для малых каналов дает в среднем заниженные данные с максимальным отклонением при низких температурах [22]. На данный момент ни один из существующих расчетных методов по определению градиента давления при двухфазном течении вещества в малых каналах не является универсальным. Наилучшую сходимость расчетных и экспериментальных величин показывают следующие методы: расчетное соотношение Мюллера—Стейнхагена и Хека для традиционных больших труб, корреляция для малых каналов Мишима—Хибики метода разделенных потоков Локхарта—Мартинелли (см. табл. 2) и гомогенная модель с расчетом динамической вязкости потока по соотношению Сичитти. Названные методы

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

Список литературы

1. KandlikarS. G. Fundamental issues related to flow boiling in minichannels and microchannels // Experimental thermal and fluid science. 2002. Vol. 26.

2. Mehendal S. S., Jacobi A. М., Shah R. K. Fluid flow and heat transfer at micro- and meso-scales with application to heat exchanger design // Appl. mech. rev. 2000. Vol. 53.

3. Kew P. A., Cornwell K. Correlations for the prediction of boiling heat transfer in small-diameter channels // Applied thermal engineering. 1997. Vol. 17.

4. Carey P. Van. Liquid vapor phase change phenomena: An introduction to the thermophysics of vaporization and condensation processes in heat transfer equipment // 2nd ed. Taylor & francis group, LLC. 2008.

5. Thome R. J. Engineering data book 111// Wolverine tube, Inc. 2004-2010.

6. Suo М., Griffith P. Two-phase flow in capillary tubes // ASMEJ. Basic eng. 1964. Vol. 86.

7. Yang C.-Y., Shieh C.-C. Flow pattern of air—water and two-phase R134a in small circular tubes // International Journal of Multiphase flow. 2001. Vol. 27.

8. Revelline R. Experimental two-phase fluid flow in microchannels // Ph. D. thesis N 3437. Ecole polytechnique federale de lausanne. 2005.

9. Characterization of diabatic two-phase flows in microchannels: Flow pattern results for R134a in a 0,5 мм channel / R. Revellin, V. Dupont, T. Ursenbacher, R. J. Thome, I. Zun // International journal of multiphase flow. 2006. Vol. 32.

10. Jassim E. W., Newell T. A. Prediction of two-phase pressure drop and void fraction in microchannels using probabilistic flow regime mapping // International journal of heat and mass transfer. 2006. Vol. 49.

11. Garimella S., KillionJ. D., Coleman J. W. Experimentally validated model fortwo-phase pressure drop in the intermittent flow regime for circular microchannels // Journal of fluids engineering, transactions of the ASME. 2002. Vol. 124 (1).

12. Chung P. M.-Y., Kawaji M. The effect of channel diameter on adiabatic two-phase flow characteristics in

microchannels // International journal of multiphase flow. 2004. Vol. 30.

13. Thome J. R., Dupont V., Jacobi A. M. Heat transfer model for evaporation in microchannels. Part I: Presentation of the model 11 International journal of heat and mass transfer. 2004. Vol. 47.

14. Cioncolini A., Thome R. J., Lombardi C. Unified macro-to-microscale method to predict two-phase frictional pressure drops of annular flows // International journal of multiphase flow. 2009. Vol. 35.

15. Lockhart R. W., Martinelli R. C. Proposed correlation of data for isothermal two-phase, two-component flow in pipes 11 Chem. eng. prog. 1949. Vol. 45.

16. Revellin R., Thome J. R. Adiabatic two-phase frictional pressure drops in microchannels // Experimental thermal and fluid science. 2007. Vol. 31.

17. Sun L., Mishima K. A. Evaluation analysis of prediction methods fortwo-phase flow pressure drop in mini-channels// International journal of multiphase flow. 2009. Vol. 35.

18. Lee H. J., Lee S. Y. Pressure drop correlations for two-phase flow within horizontal rectangular channels with small heights // International journal of multiphase flow. 2001. Vol. 27.

19. Ribatski G., Wojtan L., Thome J. R. An analysis of experimental data and prediction methods for two-phase frictional pressure drop and flow boiling heat transfer in microscale channels // Experimental thermal and fluid science. 2006. Vol. 31 (1).

20. Lombardi C., Carsana C. G. A dimensionless pressure drop correlation for two-phase mixtures flowing up flow in vertical ducts covering wide parameter range // Heat and technology. 1992. Vol. 10(1-2).

21. Tu X., Hrnjak P. S. Pressure drop and visualization of R 134a two-phase flow in a rectangular microchannel. Part I // ASH RAE Transactions. 2003. Vol. 109.

22. Carbon dioxide flow boiling in a single microchannel. Part I: Pressure drops / M. Ducoulombier, S. Colasson, J. Bonjour, P. Haberschill // Experimental thermal and fluid science. 2011. Vol. 35 (4).

23. Flow-regime-based model for pressure drop predictions in microchannels / G. V. Nino, W. E. Jassim, P. S. Hrnjak, A. T. Newell // HVAC& R Research. 2006. Vol. 12(1).

i Надоели баннеры? Вы всегда можете отключить рекламу.