Экспериментальное исследование термокапиллярной конвекции в растворе бромида лития при неизотермической абсорбции

Вестник МГТУ. 2016. Т. 19, № 3. С. 673-679. DOI: 10.21443/1560-9278-2016-3-673-679

УДК 536.423.4

Н. С. Буфетов, Р. А. Дехтярь, В. В. Овчинников

Экспериментальное исследование термокапиллярной конвекции в растворе бромида лития при неизотермической абсорбции

Целью данной работы являлось получение новых экспериментальных данных по развитию термокапиллярной конвекции в растворе бромида лития при неизотермической абсорбции водяного пара в ограниченном сосуде. При проведении экспериментов использовались бесконтактный метод измерения температуры межфазной поверхности, основанный на использовании тепловизионной техники, и метод визуализации течения, возникающего в слое абсорбента, путем введения меток в раствор. В качестве объекта исследования был выбран неподвижный слой раствора бромида лития с начальной концентрацией 58 % и толщиной слоя 20 мм в ограниченном сосуде диаметром 70 мм. В качестве поглощаемого газа использовался водяной пар, давление которого поддерживалось постоянным во время абсорбции и составляло примерно 2 000 Па. Процесс абсорбции начинался с подачи пара в абсорбер, начальное давление в котором соответствовало давлению насыщенных паров воды при данной температуре раствора (250-300 Па). На основании визуализации течения, возникающего вследствие термокапиллярной конвекции, измерены скорость и продолжительность конвективного движения в слое абсорбента. При развитии термоконцентрационной конвекции получены термограммы межфазной поверхности для различных режимных параметров. Определены характерные времена развития термокапиллярной конвекции. На основании анализа полученных экспериментальных данных предложены эмпирические выражения для продолжительности конвективного движения в зависимости от давления пара в абсорбере и толщины слоя абсорбента. Показано, что понижение начальной температуры раствора приводит к интенсификации термоконцентрационной конвекции. Результаты работы являются полезными для проектирования и оптимизации элементов холодильной техники абсорбционного типа.

Ключевые слова: абсорбция, тепло- и массообмен, раствор, термокапиллярная конвекция.

Введение

В абсорбционных холодильных аппаратах и тепловых насосах используется явление физической абсорбции, сопровождающееся выделением большого количества тепла при поглощении пара раствором абсорбента. В таких аппаратах процесс абсорбции осуществляется на пленках раствора, стекающих по поверхностям, через которые отводится выделившееся при абсорбции тепло. Большинство теоретических работ посвящено исследованию именно пленочной абсорбции [1-4]. В современной энергетике, химической и перерабатывающей промышленности важной задачей является интенсификация тепло- и массопереноса, поиск новых способов увеличения производительности, снижения металлоемкости именно абсорберов как основного элемента холодильной машины или теплового насоса, которые менее всего изучены и исследованы. Одним из таких способов является введение в раствор поверхностно-активного вещества (ПАВ), которое приводит к возникновению приповерхностной конвекции и к интенсификации процессов переноса. Однако в экспериментах с использованием ПАВ трудно разделить вклад в тепло- и массоперенос конвективного и вынужденного течения, поэтому удобным объектом исследования является абсорбция пара неподвижным слоем раствора, когда вынужденное течение отсутствует [5-9].

В работе [5] показано, что в начальный период абсорбции, когда тепловой и диффузионный слои достаточно тонкие, малые градиенты концентрации и температуры приводят к возникновению движения жидкости в слое абсорбента. В работе [6] показано, что при введении в раствор ПАВ в неподвижном слое жидкости возникает интенсивное конвективное перемешивание, что приводит к интенсификации переноса теплоты и массы при абсорбции по сравнению с раствором без ПАВ [7-9]. При неизотермической абсорбции в отсутствии ПАВ имеется градиент температуры и концентрации абсорбента на поверхности жидкости. Это также может привести к появлению движения жидкости в приповерхностном слое (термоконцентрационная конвекция), рассмотренная подробно в [10, § 68]. Известно из [7], что на начальном этапе процесса абсорбции на границе раздела фаз происходит резкое повышение температуры и уменьшение концентрации раствора. Неоднородность температуры и концентрации на межфазной поверхности даже в отсутствие ПАВ вызывает движение жидкости в приповерхностном слое вследствие градиента поверхностного натяжения. Следует отметить, что при неизотермической абсорбции существенное влияние на развитие конвективного течения оказывают наличие мениска вблизи стенки сосуда с раствором абсорбента, а также тепловые условия на стенке.

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

Материалы и методы

Экспериментальное исследование проводилось на установке (рис. 1) по методике, подробно изложенной в [7]. В качестве абсорбента использовался водный раствор бромида лития LiBr, а в качестве

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

9

Рис. 1. Схема экспериментальной установки: 1 - вакуумная камера; 2 - кювета; 3 - смотровые и инфракрасные окна; 4 - инфракрасная (ИК) камера;

5 - видеокамера; 6 - световой нож; 7 - генератор пара; 8 - клапан; 9 - персональный компьютер

Установка состояла из абсорбера 1, в который помещалась заполненная абсорбентом круглая кювета 2 диаметром 70 мм. Дно кюветы изготовлено из оптического стекла толщиной 15 мм, а стенки из тонкой (40 мкм) прозрачной пленки. Абсорбер имел дно, верхнюю крышку и боковые стенки 3, прозрачные как в видимом свете, так и в инфракрасном спектре. Это позволяло проводить измерения температуры поверхности с помощью инфракрасной камеры 4. Регистрация течения в поверхностном слое проводилась видеокамерой 5 с использованием как горизонтальной, так и вертикальной подсветки 6. Генератор водяного пара 7 соединялся с абсорбером паропроводом и клапаном 8, отсекающим их друг от друга. Экспериментальная установка была оснащена автоматизированной системой поддержания заданного давления с точностью ±5 Па, а также системой автоматического сбора и обработки информации 9.

Основным параметром, определяющим термоконцентрационное течение, является величина Л% = %0 - %е - разность между начальной концентрацией абсорбента %0 и равновесной концентрацией %е, соответствующей температуре межфазной поверхности Те при заданном давлении Р. Величина Л% варьировалась за счет изменения давления абсорбируемого пара Р и начальной температуры Т0. Начальная концентрация раствора бромида лития %0 = 58.0 ± 0.1 %, толщина слоя раствора ^ = 20 ± 0.1 мм.

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

Перед началом каждого эксперимента производилось удаление неабсорбируемых примесей и длительное термостатирование раствора. Процесс абсорбции начинался с открытия клапана 8 и подачи водяного пара из генератора 7 в абсорбер 1. При этом синхронно начиналась запись на видеокамеру 5 со скоростью 60 к/с через боковое смотровое окно 3 камеры 1 возникающих конвективных движений в слое жидкости и на инфракрасную камеру 4 через верхнее специальное окно 3 термограмм поверхности раздела фаз.

Результаты и обсуждение

Была проведена визуализация поля скоростей в растворе бромида лития при неизотермической абсорбции водяного пара в ограниченном сосуде. Для примера на рис. 2 приведены характерные кадры видеосъемки с визуализационными метками для различных моментов времени с начала абсорбции, для значений начальной температуры раствора в кювете 20.0 °С при Р = 2 000 Па, ^ = 20 мм. На рис. 2, а показана картинка при наложении нескольких кадров с интервалом между ними 0.5 с, на рис. 2, Ь - с интервалом 1 с. Судя по рисункам, возникающее при абсорбции движение меток наблюдалось в приповерхностном слое жидкости толщиной 5-7 мм. При этом поверхность раздела фаз на рис. 2 расположена в самой верхней части рисунка. Расчет скоростей движения меток производился по стандартной методике: после наложения ряда видеокадров друг на друга измерялось расстояние, которое проходили выделенные метки на кадрах, а время прохождения меток определялось по номерам видеокадров. На рис. 2 видно, что в первые 2.5 с с начала процесса абсорбции при приближении к поверхности расстояние между метками увеличивается.

Вестник МГТУ. 2016. Т. 19, № 3. С. 673-679. Б01: 10.21443/1560-9278-2016-3-673-679

Это означает, что с удалением от поверхности скорость жидкости уменьшается. На рис. 2, Ь показаны треки частиц на интервале времени от 2.5 до 12 с в наблюдаемой области толщиной 5-6 мм. Из рисунка видно, что вследствие прорастания температурного слоя движением охватывается вся наблюдаемая область. При этом в наблюдаемой области происходит вырывание скоростей движения меток по глубине, а также дальнейшее уменьшение скорости движения меток с течением времени.

Рис. 2. Треки частиц в слое жидкости на начальной стадии абсорбции при Т0 = 20.0 °С: а - 0 < t < 2,5 с

На рис. 3 представлены результаты обработки видеоматериалов по измерению продолжительности конвективного движения т в зависимости от давления абсорбции Р для раствора бромида лития толщиной слоя 20 ± 0.1 мм. Начальные концентрация и температура составили = 58 % и Т0 = 22 ± 0.5 °С, что соответствует равновесному давлению насыщенных паров над раствором Р = 255 ± 10 Па. Из рисунка видно, что при значениях давления меньше 500 Па движение жидкости в слое абсорбента слабо выражено. Дальнейшее повышение давления приводит к развитию термоконцентрационного конвективного течения, продолжительность которого увеличивается с ростом давления и изменяется от нескольких секунд до нескольких минут. Следует отметить, что при Р > 1 500 Па имеет место большой разброс экспериментальных данных, что объясняется спонтанностью возникновения движения в слое абсорбента. Здесь же можно выделить две зоны возникновения неустойчивостей, связанных с горизонтальными и вертикальными потоками в слое жидкости: верхние точки 1 на рис. 3 свидетельствуют о продолжительном перемешивании в слое именно в вертикальном направлении, т. е. сверху от границы раздела фаз ко дну кюветы. Нижние точки 3 на графике относятся к случаю, когда перемешивание происходило только в тонком приповерхностном слое в горизонтальном направлении. Данные для вертикальной конвекции интерполируются эмпирической зависимостью т = 7 Р05, а для горизонтальной конвекции т : = 4 • 10-5 Р2. В этих зависимостях численные коэффициенты приведены для давления в [Па], а времени в [с].

Рис. 3. Зависимость продолжительности течения т от давления абсорбируемого пара:

* -5 2 * 0 5

1, 3 - опытные данные, 2 - интерполяция т= 4 • 10 Р , 4 - интерполяция т = 7 Р

В результате обработки массива видеоданных для треков частиц в слое жидкости были рассчитаны вертикальные Уу и горизонтальные Ух компоненты скорости их движения в зависимости от времени (рис. 4), при этом ось х направлена горизонтально, а ось у - вертикально. На рис. 4 показаны результаты для температуры раствора Т0 = 24.5 °С и давления Р = 2 030 Па (продолжительность движения меток составляла примерно 145 с). Видно, что в верхнем слое жидкости в первые 40 с движение в основном происходит в горизонтальном направлении. На этом интервале времени в наблюдаемой области достигаются наибольшие

значения горизонтальной скорости (до 0.8 мм/с) в направлении противоположном оси х. Затем, в течение последующих 20 с, преобладает движение в направлении оси х. В дальнейшем в исследуемой области нельзя выделить преимущественного движения меток: одновременно наблюдается движение в различных направлениях по оси х, при этом средняя скорость | V] движения снижается. Вертикальные составляющие скорости ¥х в основном имели значения на порядок меньше, чем горизонтальные ¥у. Выделить преобладающее направление вертикальных составляющих скорости в какой-либо момент времени на основании этой серии опытов не представляется возможным. Следует отметить, что со временем движение меток замедляется и при t > 150 с практически прекращается.

Рис. 4. Зависимость скорости течения раствора в наблюдаемой области от времени: ¥х - горизонтальная составляющая скорости; Уу - вертикальная составляющая; | V] - абсолютные значения скорости

При фиксированных значениях Р и Л% уменьшение начальной толщины абсорбента ^ от 20 мм до 3 мм приводит к увеличению продолжительности конвективного течения (рис. 5). Данные для зависимости

* 2

продолжительности конвективного движения т от толщины слоя о интерполируются эмпирической зависимостью т = 95 + 5 000 • ехр[-(8/0.5-0.4)°'5]. В этой зависимости численные коэффициенты приведены для толщины слоя [8] в [мм], времени в [с].

10000

0 5 10 15 20

5, мм

Рис. 5. Зависимость продолжительности конвективного течения т* от толщины слоя: = 58 %, Т0 = 22 °С, Р = 2 000 Па: 1 - опытные данные, 2 - интерполяция т* = 95 + 5 000-ехр[-(8/0.5-0.4)°'5]

Для раствора бромида лития при P = const коэффициент поверхностного натяжения возрастает с ростом температуры на поверхности, находящейся в состоянии термодинамического равновесия [12]. Так как движение слоя жидкости вследствие термоконцентационной конвекции направлено от мест с меньшим поверхностным натяжением к местам с большим поверхностным натяжением [10], можно утверждать, что жидкость будет двигаться от более холодной зоны на поверхности раздела фаз к более горячей зоне. Непрерывные измерения температуры поверхности абсорбирующего раствора с помощью тепловизионной техники позволили подробно проследить за изменением температурного поля межфазной поверхности.

Вестник МГТУ. 2016. Т. 19, № 3. С. 673-679. DOI: 10.21443/1560-9278-2016-3-673-679

На рис. 6 приведены термограммы межфазной поверхности раствора через 0.3 с после начала абсорбции, ^о = 58 %, Р = 2 000 Па. Показаны две характерные термограммы: одна при начальной температуре раствора Т0 = 24.5 °С (рис. 6, а) и при Т0 = 20.4 °С (рис. 6, Ь). Можно видеть, что для раствора с меньшей начальной температурой наблюдаются наиболее заметные неоднородности распределения температуры по поверхности, которые вызваны тепловыделением при абсорбции над мениском около стенки кюветы. Такое температурное поле приводит к концентрационной неоднородности на поверхности, вызывающей термоконцентрационную конвекцию различной интенсивности. Таким образом, чем ниже начальная температура раствора, тем выше возникающие градиенты температур на поверхности раздела фаз и тем выше сорбционные свойства раствора.

Рис. 6. Распределение температуры на поверхности раствора через 0.3 с после начала абсорбции, = 58 %, Р = 2 000 Па: а - Т0 = 24.5 °С; Ь - Т0 = 20.4 °С

Заключение

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

Библиографический список

1. Накоряков В. Е., Григорьева Н. И. Точное решение задачи о совместном тепломассопереносе при пленочной абсорбции // Инж. физ. журнал. 1977. Т. 33 (5). С. 893-896.

2. Накоряков В. Е., Григорьева Н. И. Расчет тепломассообмена при неизотермической абсорбции на начальном участке стекающей пленки // Теоретические основы химической технологии. 1980. Т. 14 (4). С. 483-488.

3. Grossman G. Simultaneous heat and mass transfer in film absorption under laminar flow // Int. J. Heat and Mass Transfer. 1983. V. 26 (3). P. 357-371.

4. Накоряков В. Е., Григорьева Н. И. Тепломассообмен при пленочной абсорбции с изменением объема жидкой фазы // Теоретические основы химической технологии. 1995. Т. 29 (3). С. 242-248.

5. Daiguji H., Hihara E., Saito T. Mechanism of absorption enhancement by surfactant // Int. J. Heat and Mass Transfer. 1997. V. 40 (8). P. 1743-1752.

6. Nakoryakov V. E., Bufetov N. S., Dekhtyar R. A. Effect of surfactant added in small amounts on nonisothermal absorption: an experimental study // Journal of Applied Mechanics and Technical Physics. 2004. V. 45 (2). P. 276-280.

7. Heat and mass transfer at vapor absorption by an immobile layer of solution / V. E. Nakoryakov, N. S. Bufetov, N. I. Grigoryeva, R. A. Dekhtyar // J. Applied Mechanics and Technical Physics. 2003. V. 44 (2). P. 111-118.

8. Буфетов Н. С., Дехтярь Р. А. Абсорбция газов, сопровождающаяся значительным выделением тепла // Ползуновский вестник. 2004. № 1. С. 53-56.

9. Vapor absorption by immobile solution layer / V. E. Nakoryakov, N. I. Grigoryeva, N. S. Bufetov, R. A. Dekhtyar, I. V. Marchuk // Int. J. Heat and Mass Transfer. 2004. V. 47 (6-7). P. 1525-1533.

10. Левич В. Г. Физико-химическая гидродинамика. М. : Гос. изд-во физико-математ. литературы. 1959. 699 с.

11. Löver H. Thermodynamischen und physikalische Eigenschaften der wässrigen Lithiumbromid Lösung. Dissertation. Karlsruhe. 1960. 116 S.

12. Grigoryeva N. I. About the methods for determination of Marangoni numbers at investigation of absorption under the conditions of heat pump operation // Theoretical Foundations of Chemical Engineering. 2005. V. 39 (6). P. 561-565.

References

1. Nakoryakov V. E., Grigoreva N. I. Tochnoe reshenie zadachi o sovmestnom teplomassoperenose pri plenochnoy absorbtsii [The exact solution of the problem of heat and mass transfer at the joint of the film absorption] // Inzh. fiz. zhurnal. 1977. V. 33 (5). P. 893-896.

2. Nakoryakov V. E., Grigoreva N. I. Raschet teplomassoobmena pri neizotermicheskoy absorbtsii na nachalnom uchastke stekayuschey plenki [The calculation of heat and mass transfer under non-isothermal absorption in the initial part of the flowing-down film] // Teoreticheskie osnovy himicheskoy tehnologii. 1980. V. 14 (4). P. 483-488.

3. Grossman G. Simultaneous heat and mass transfer in film absorption under laminar flow // Int. J. Heat and Mass Transfer. 1983. V. 26 (3). P. 357-371.

4. Nakoryakov V. E., Grigoreva N. I. Teplomassoobmen pri plenochnoy absorbtsii s izmeneniem ob'ema zhidkoy fazy [Heat and mass transfer at film absorption with changes in the volume of the liquid phase] // Teoreticheskie osnovy himicheskoy tehnologii. 1995. V. 29 (3). P. 242-248.

5. Daiguji H., Hihara E., Saito T. Mechanism of absorption enhancement by surfactant // Int. J. Heat and Mass Transfer. 1997. V. 40 (8). P. 1743-1752.

6. Nakoryakov V. E., Bufetov N. S., Dekhtyar R. A. Effect of surfactant added in small amounts on nonisothermal absorption: an experimental study // Journal of Applied Mechanics and Technical Physics. 2004. V. 45 (2). P. 276-280.

7. Heat and mass transfer at vapor absorption by an immobile layer of solution / V. E. Nakoryakov, N. S. Bufetov, N. I. Grigoryeva, R. A. Dekhtyar // J. Applied Mechanics and Technical Physics. 2003. V. 44 (2). P. 111-118.

8. Bufetov N. S., Dehtyar R. A. Absorbtsiya gazov, soprovozhdayuschayasya znachitelnym vydeleniem tepla [Absorption of gases accompanied by a significant heat release] // Polzunovskiy vestnik. 2004. N 1. P. 53-56.

9. Vapor absorption by immobile solution layer / V. E. Nakoryakov, N. I. Grigoryeva, N. S. Bufetov, R. A. Dekhtyar, I. V. Marchuk // Int. J. Heat and Mass Transfer. 2004. V. 47 (6-7). P. 1525-1533.

10. Levich V. G. Fiziko-himicheskaya gidrodinamika [Physical and chemical hydrodynamics]. M. : Gos. izd-vo fiziko-matemat. literatury. 1959. 699 p.

11. Löver H. Thermodynamischen und physikalische Eigenschaften der wässrigen Lithiumbromid Lösung. Dissertation. Karlsruhe. 1960. 116 S.

12. Grigoryeva N. I. About the methods for determination of Marangoni numbers at investigation of absorption under the conditions of heat pump operation // Theoretical Foundations of Chemical Engineering. 2005. V. 39 (6). P. 561-565.

Сведения об авторах

Буфетов Николай Сергеевич - пр-т Академика Лаврентьева, 1, г. Новосибирск, Россия, 630090; Институт теплофизики им. С. С. Кутателадзе СО РАН, канд. техн. наук, ст. науч. сотрудник; e-mail: bufetov@itp.nsc.ru

Bufetov N. S. - 1, Academician Lavrent'ev Avenue, Novosibirsk, Russia, 630090; Institute of Thermophysics SB RAS, Cand. of Tech. Sci., Senior Researcher; e-mail: bufetov@itp.nsc.ru

Дехтярь Руслан Антальевич - пр-т Академика Лаврентьева, 1, г. Новосибирск, Россия, 630090; Институт теплофизики им. С. С. Кутателадзе СО РАН, канд. техн. наук, ст. науч. сотрудник; e-mail: dekhtyar@itp.nsc.ru

Dekhtyar R. A. - 1, Academician Lavrent'ev Avenue, Novosibirsk, Russia, 630090; Institute of Thermophysics SB RAS, Cand. of Tech. Sci., Senior Researcher; e-mail: dekhtyar@itp.nsc.ru

Овчинников Валерий Викторович - пр-т Академика Лаврентьева, 1, г. Новосибирск, Россия, 630090; Институт теплофизики им. С. С. Кутателадзе СО РАН, канд. техн. наук, ст. науч. сотрудник; e-mail: avk@itp.nsc.ru

Ovchinikov V. V. - 1, Academician Lavrent'ev Avenue, Novosibirsk, Russia, 630090;

Institute of Thermophysics SB RAS, Cand. of Tech. Sci., Senior Researcher; e-mail: avk@itp.nsc.ru

BecTHHK MrTy. 2016. T. 19, № 3. C. 673-679. DOI: 10.21443/1560-9278-2016-3-673-679

N. S. Bufetov, R. A. Dekhtyar, V. V. Ovchinnikov

An experimental investigation of thermo-capillary convection in solution of lithium bromide at non-isothermal absorption

The aim of this work is to obtain new experimental data on the development of thermo-capillary convection in a solution of lithium bromide with non-isothermal absorption of water vapor in a confined vessel. In experiments the contactless method of measuring the temperature of the interface based on using thermal imaging equipment, and the visualization method of the flow arising in the absorbent layer by introducing markers into solution have been used. The immobile layer of lithium bromide solution with an initial concentration of 58 % and a layer thickness of 20 mm in a confined vessel of 70 mm in diameter has been chosen as an object of investigation. The water vapor has been used as the absorbed gas, whose pressure was kept constant during the absorption (approximately 2 000 Pa). The absorption process has started after the steam supply to the absorber where the initial pressure has corresponded to saturated water vapor pressure at the temperature of the solution (250-300 Pa). The speed and duration of convective motion in the layer absorption on the basis of imaging arising due to thermo-capillary convection have been measured. The thermograms of the interface for various operating parameters have been obtained with the development of thermo-concentration convection. The characteristic time of the thermo-capillary convection development has been determined. Based on the experimental data analysis the empirical expressions for duration of the convective motion depending on the vapor pressure in the absorber and the thickness of the absorbent layer have been proposed. It has been shown that lowering the initial temperature of the solution leads to enhancement of thermo-concentration convection. The results are useful for the design and optimization of elements of refrigeration engineering of the absorption type.

Key words: absorption, heat and mass transfer, solution, thermocapillary convection.