Термографическое исследование пленки жидкости стекающей по поверхности с локальным источником тепла

Актуальность работы.

Использование тонких пленок жидкости. Тонкие пленки жидкости широко используются в технологических процессах, так как обеспечивают высокую интенсивность тепло-массопереноса и значительную поверхность контакта фаз при малых удельных расходах жидкости, и в ряде случаев незначительные потери давления [Тананайко и Воронцов 1975, Алексеенко Накоряков Покусаев 1992]. Пленочные течения специально создаются в различных аппаратах химической технологии, пищевой, фармацевтической промышленности, в криогенной индустрии. Гравитационно стекающие пленки, например, используются в испарителях низкого давления, применяемых для концентрирования пищевых продуктов, в аппаратах для опреснения морской воды, а также в ректификационных колоннах. Перспективными являются пленочные системы охлаждения электронных компонентов. Тонкие пленки жидкости могут также возникать при движении двухфазных потоков в каналах испарительно-конденсационных систем. Снижение толщины стекающего слоя жидкости позволяет интенсифицировать процесс теплообмена, однако тонкие пленки подвержены разрывам. Свободные от жидкости участки поверхности исключаются из процесса тепло-массообмена, что резко снижает эффективность аппаратов. Чтобы избежать разрывов, в пленочных аппаратах обычно расход жидкости в несколько раз превосходит количество испаряемой жидкости.

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

В 1994 г. в Институте теплофизики СО РАН д.ф.-.м.н O.A. Кабовым было обнаружено явление формирования пространственной самоорганизующейся структуры в тонкой движущейся под действием гравитации пленке жидкости при ее локальном нагреве со стороны подложки. Структуры представляют собой вал жидкости в области верхней кромки нагревателя, из которого с определенной периодичностью стекают струи жидкости и тонкая пленка между ними. Гипотезу о термокапиллярной природе механизма образования структур было предложено проверить на основании измерений температуры поверхности пленки с помощью инфракрасного сканера.

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

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

Научная новизна работы заключается в том, что автором впервые:

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

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

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

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

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

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

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

Публикации и апробация работы. По теме диссертации опубликовано 1 печатных работ. Результаты работы докладывались автором диссертации на Международной конференции молодых ученых «Актуальные вопросы теплофизики физической гидрогазодинамики» (Новосибирск, 1997), на семинаре «Физическа гидродинамика» в Институте теплофизики под руководством чл.-корр. РАН С. 1.

Алексеенко (Новосибирск 1997, 2000), на семинаре «Прикладная гидродинамика» в Институте гидродинамики им. М. А. Лаврентьева СО РАН под руководством чл.-корр. РАН В. В. Пухначева (Новосибирск 2000), на VI Всероссийской конференции молодых ученых «Актуальные вопросы теплофизики и физической гидрогазодинамики» (Новосибирск, 2000). Так же результаты, изложенные в диссертации, докладывались соавторами работ на различных международных и российских конференциях и семинарах.

Личное участие автора. Данная работа выполнена в лаборатории «интенсификации процессов теплообмена» (1994;2000, заведующий д.ф.-м.н. O.A. Кабов) Института теплофизики СО РАН. В диссертации изложены результаты исследований, выполненных автором самостоятельно и в сотрудничестве с сотрудниками лаборатории. Лично автором выполнено экспериментальное исследование температурного поля на поверхности пленки жидкости и численные расчеты течения локально нагреваемой стекающей пленки. Анализ и обобщение экспериментальных данных, а также сопоставление с численными расчетами выполнено совместно с O.A. Кабовым. Кроме того, И. В. Марчуку принадлежит методика термографического исследования.

Выводы.

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

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

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

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

ОСНОВНЫЕ ОБОЗНАЧЕНИЯ, а — коэффициент температуропроводности, м /с В — ширина нагревателя, м В=аЫХу/ - число Био сризобарная теплоемкость жидкости, Дж/(кг К) Сконцентрация, %.

Сгкритерий деформируемости поверхности термокапиллярными силами, =-итЛТ1оо <Шч т /¿-А, -спектральная поверхностная плотность излучения абсолютно черного тела, Вт/м3 gускорение свободного падения, м/с2 ктолщина пленки, м г0 — начальная толщина пленки, м.

О 1 / ^.

К= (у / g) — масггаб вязко-гравитационного взаимодействия, м 1а = сг / pg — масштаб капиллярно-гравитационного взаимодействия, м Ьдлина нагревателя, характерный линейный размер, м Ьь — длина начального участка развития теплового пограничного слоя, м Ьудлина гладкой, безволновой зоны течения пленки, м.

ЛД (<�За1дТдТ1ду)к2.

Мп=—число Марангони ра.

Ии — число Нуссельта.

N — число узлов сетки по координате X. р — давление, Н/м2 ра — атмосферное давление, Н/м2.

Ре — число Пекле.

Рг = рср / Я — число Прандтля — плотность теплового потока, Вт/см2.

— мощность, выделяемая нагревателем, Вт.

Ке=Г/рчисло Рейнольдса пленки.

Яе^, — число Рейнольдса при ??=?/2.

Яе/, — число Рейнольдса при ??=1.

Ттемпература, °С.

То — температура натекающей пленки жидкости, °С.

Ти — средневзвешенная температура пленки, °С.

Ттах — максимальная температура на поверхности пленки, °С.

Т — температура поверхности пленки, °С и, V — компоненты скорости в направлении х, ух соответственно, м/с и&bdquo-, у0 — масштаб скорости жидкости, и0=5х (др%к ()2/р0, у0=е и0, м/с и, V — безразмерные компонент ы скорости в направлении х, у, г соответственно х, у, г — декартовы координаты, м.

X, У, 2- безразмерные декартовы координаты.

Греческие символы акоэффициент теплоотдачи, Вт/(м2 К) а (Я) — коэффициент поглощения жидкост и, м" 1 Г — удельный расход жидкост и, кг/мс.

Гу — объемный расход жидкости на единицу ширины пленки, м /с.

АЛ — диапазон длин волн, регистрируемый прибором.

Л) — коэффициент излучения жидкости, безразмерный е^(Л) — коэффициент излучения стенки, безразмерный s — соотношение линейных масштабов в пленке.

Лкоэффициент теплопроводности, Вт/мК.

X — функция Хевисайда.

0- угол наклона пластины к горизонту, градус втемпература, безразмерная.

Л — длина волны, расстояние между гребнями струй, м.

1 — динамическая вязкость жидкости, кг/мс л0 — динамическая вязкость жидкости, =ju (T0), кг/мс д — динамическая вязкость жидкости, = ц / /л0, безразмерная vкоэффициент кинематической вязкости, м /с. рплотность, Kr/MJ сткоэффициент поверхностного натяжения, Н/ м утлинейный коэффициент, определяющий зависимость поверхностного натяжения от температуры, Н/мК ткасательное напряжение, Н/м.

Нижние индексы.

Оначальное значение величины или величина для вертикальной поверхности рассчитаная по теории Нуссельта surвеличина на поверхности пленки Wстенка х, у, г — производная величины по переменной х, у и касательной к поверхности.