CC BY
104
УДК 532.6
ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ПОВЕДЕНИЯ ЛЕЖАЩЕЙ КАПЛИ НА НАГРЕТОЙ ПЛОСКОСТИ1
А.Е. Коренченко2, А.Ж. Илимбаева*, В.П. Бескачко4
Проведено математическое моделирование поведения капли вязкой жидкости, лежащей на нагретой поверхности в атмосфере нейтрального газа. Рассмотрен динамический режим с постоянным контактным углом. Поверхность предполагается гладкой, химически инертной и изотермической, газ - идеальным и неконденсирующимся. Жидкая капля описана в приближении несжимаемой жидкости с учетом конвекции Марангони и гравитационных эффектов, в газе модель учитывает гравитационную конвекцию. Проведено исследование течений в капле и окружающей газовой среде и проанализированы движения поверхности капли.
Ключевые слова: поверхностное натяжение, лежащая капля, численное моделирование.
Введение
Исследования динамического поведения лежащей капли в настоящее время вызывают большой интерес. Это связано с важной ролью, которую этот объект играет в технологических процессах, например, при распылении чернил в струйных принтерах, капельном охлаждении поверхностей, росте биологических кристаллов, в тепло-массообменных аппаратах. Математическое описание капли осуществить сложно уже потому, что межфазная поверхность жидкость-газ может свободно двигаться, что подразумевает перестройку вычислительной области на каждом временном шаге. Задача еще более усложняется, если учесть возможность испарения с поверхности жидкости и температурную неоднородность, связанную с поглощением теплоты парообразования на границе жидкость-газ.
В работах по визуализации испаряющейся лежащей капли докладывают о существовании внутренних течений в капле [1], наличии гидротермических волн [2]. В [3] обнаружено уплощение капли в сравнении с равновесной формой при уменьшении объема жидкости из-за испарения. Эксперименты подтверждают факт, что лежащая капля находится в неравновесном состоянии даже в атмосфере своего насыщенного пара, так как из-за кривизны поверхности давление насыщенного пара на ней немного превышает таковое в объеме.
Поверхность испаряющейся капли может перемещаться не только за счет фазового перехода, но также из-за присутствия внутренних течений. Обычно в вычислительных работах неравновес-ностью формы испаряющейся капли пренебрегают. Профиль капли считается равновесным и сферическим (линзоподобным), либо рассчитывается на каждом временном шаге с помощью уравнения Лапласа. Как показывает сравнение результатов численных расчетов с экспериментальными данными, на скорость испарения эти допущения влияют мало, причем как для режима испарения с постоянным радиусом пятна смачивания, так и с постоянным контактным углом [46]. Но существует еще метрологический аспект задачи об испарении лежащей капли. Он связан с тем, что движение неравновесной из-за испарения границы раздела может повлиять на точность измерения поверхностного натяжения методом лежащей капли.
Глобальная цель проводимых авторами исследований - анализ отклонения формы испаряющейся капли от равновесной. Авторы проводят работы поэтапно. В [7] было проведено исследование свободных колебаний лежащей капли, выведенной из положения механического равновесия, и рассчитаны параметры колебаний для капель различных жидкостей. Цель настоящей работы состоит в выяснении влияния температурной неоднородности (нагретая подложка) на движение границы раздела жидкость-газ.
1 Работа поддержана грантом РФФИ № 10-03-00719-а.
2 Коренченко Анна Евгеньевна - доктор физико-математических наук, старший научный сотрудник, ИМЕТ УрО РАН.
E-mail: korenchenko@physics.susu.ac.ru
3 Илимбаева Айгуль Жепспаевна - студентка, физический факультет, Южно-Уральский государственный университет.
4 Бескачко Валерий Петрович - доктор физико-математических наук, профессор, кафедра общей и теоретической физики, Южно-
Коренченко А.Е., Илимбаева А.Ж., Численное исследование поведения лежащей капли
Бескачко В.П. на нагретой плоскости
Схема численного эксперимента приведена на рис. 1. Капля жидкости объемом V = 3РЯ3 ,
Я - радиус равновеликой сферической капли, находится на твердой подложке и окружена некон-денсирующимся газом. Границы вычислительной области являются удаленными, т.е. Я, Н£ » Я . Приняты следующие приближения: температура подложки постоянна; отсутствует
испарение; предполагается осевая симметрия формы капли, распределения давления и течений в капле и газе; граница считается бесконечно тонкой и безинерционной; жидкость является несжимаемой средой; газопаровая смесь описывается приближением идеального газа; коэффициенты переноса являются постоянными в обеих средах; краевой угол является постоянным.
В работе приняты следующие обозначения: /?? - плотность жидкости, р., - плотность газа на удаленных границах, V/, \\, - коэффициенты кинематической \2
вязкости, //?, //„ - коэффициенты динамической ВЯЗКОСТИ, А/ ,
К - коэффициенты теплопроводности, С] , с^ - удельные теплоемкости, МI, Мё - молярные массы, Д. - коэффици- Цу ент объемного расширения ЖИДКОСТИ, О) , СГ,;Л., (Г/,. - коэф-
н
фициенты поверхностного натяжения на границах жидкость-газ, газ-подложка и жидкость-подложка соответст-
I
=) и.
венно. Значения коэффициентов переноса считаются посто- Рис Л Схема эксперимента. Капля
янными. Здесь и далее индексами I, й и 5 отмечаются вели- находится на плоской поверхности
чины, относящиеся к жидкости, газу и материалу подложки. в окружении газа
Математическая модель представляет собой систему дифференциальных уравнений сохранения и сопутствующих граничных условий. Система уравнений, описывающих течения в жидкости, содержит уравнения Навье-Стокса, неразрывности и конвективной теплопроводности
+(У .V )у = -—ур+п,у¥ - йЬ (т - т,), (1)
д t ' р1
V V = 0, (2)
—+ у (Ут )=— V 2т. (3)
д t р1с1
В (1)-(3) использованы следующие обозначения: V = {у, и} - вектор скорости, Р - давление
сверх гидростатического, Т - температура, й - ускорение свободного падения. В (1) гравитационные температурные эффекты описаны в приближении Буссинеска.
Система уравнений для газовой среды запишется как
ду+(У -V )У = - - Vр+ ^- V 2У - й, д t ' ' р р
(4)
др + ? (рУ ) = 0, (5)
+у (Ут )=-^- V 2т, (6)
д ^ рйсй
Р = р-. (7)
М
Здесь буквамир и р обозначены давление и плотность газовой среды.
В модели приняты следующие граничные условия: на оси
г = 0: ди = дт = дР = др = др = 0, V = 0, (8)
дг дг дr дг дг
на твердой границе
г = 0: V = и = 0, т = \т5, Г < 3Я (9)
т - (т, - тя )/(Як - 3Я) - (г - 3Я), 3Я < г <Я,
к,
Физика
на удаленных границах в газовой среде
г = Нк ог г = Як : V = и = 0, т = т&, р = рй, р = рй . (10)
Здесь тй, Рй, рй - есть невозмущенные значения температуры, давления и плотности газа на удаленных границах, т, - температура подложки, а условие (9) для температуры выбрано для того, чтобы избежать скачка температуры при г = 0, г = Як.
На границе жидкость-газ выполнены следующие условия [8]:
1) скорости непрерывны
VI = Vg, и, = ий ; (11)
2) нормальные напряжения испытывают скачок, определяемый формулой Лапласа
ЯТй|I -р у -п -У5к )2 -^ЯТй|й -ру -п -У5к )2^ = акк, (12)
где Т - тензор напряжений, - скорость перемещения границы раздела фаз;
3) скачок касательных напряжений равен касательной составляющей градиента поверхностного натяжения
^I, -ру 1 (У, п - % ) - (Щ. -руй 1(У2'п - % )) = О ),
В (13) температурная зависимость поверхностного натяжения описывается формулой
а = °0-От (т -т0), От = ■
4) температура непрерывна
т& = т; (14)
5) поток тепла в отсутствие испарения непрерывен
(13)
к дт -к, — дп
(
рст (угЯ-у3к)- -кг|п +рсйтй\уй-п-^)
л
й
= 0. (15)
Уравнения (12), (13), (15) записаны в системе отсчета, связанной с границей раздела фаз.
Численный метод решения
Для решения системы (1)-( 15) использовалась разностная схема, построенная на основе неравномерной пространственной сетки. Уравнения (1), (3) для жидкой фазы записывались в форме законов сохранения, дискретизировались и линеаризовывались. Система полученных алгебраических разностных уравнений решалась методом исключения Г аусса. Система уравнений (4)-(6), описывающая газовую среду, решалась методом прогонки.
В численном решении задачи можно выделить несколько этапов:
• минимизацией выражения для механической энергии [9] определяется форма равновесной капли, т.е. набор радиусов {г-} и высота Н (рис. 1). Определяются равновесные значения кривизны поверхности капли К0 (г). В отсутствие испарения капля имеет постоянный объем, поэтому вычисление {К0 (г)} проводится только один раз;
• определяются давление и скорости в капле путем решения системы уравнений (1, 2, 12, 13). Форма капли, а также скорости и плотность в газе считаются неизменными на этой стадии и берутся из предыдущего временного шага. Система разностных уравнений решается методом исключения Гаусса;
• определяется новая форма капли по формулам:
г*+А- = г- + ^ -Д-, 1 = 2, т, г-+А- = г- + и|„ -А-, 1 = 2, т +1.
Радиус пятна смачивания определяется из условия, что краевой угол не зависит от того, происходит натекание или отток жидкости на подложке, и равен равновесному углу, определяемому уравнением Юнга;
• распределение скоростей в газе определяется последовательным решением уравнений (4) с граничными условиями (11);
/
Коренченко А.Е., Илимбаева А.Ж., Бескачко В.П.
• распределение плотности в газе определяется с помощью решения уравнения (5). Вблизи поверхности жидкости значения плотности газовой среды интерполируются;
• распределение температуры в жидкости определяется посредством решения уравнения (3) с граничным условием (15);
• вычисляется распределение температуры в газе с помощью решения уравнения (6) с граничными условиями (14);
• сетка перестраивается в соответствии с изменениями в форме капли.
Таким образом, использованный численный метод является одним из методов расщепления по физическим процессам [10]. Программный продукт, реализующий указанную схему, был написан на языке РоГхап90 и распараллелен с использованием библиотеки МР1. Расчеты проводились на суперкомпьютере «СКИФ-Аврора ЮУрГУ».
Вычисления проводились для следующего набора параметров: Я = 0,5 мм, жидкость - вода с
-О§5 = 0 Н-м-1, О1й = 7,3-10-2 Н-м-1, р= 1000 кг/м-3, у1= 1,0-10-6 м2/с,
кг = 0,55 Вт / (м- К), Ь= 1,27-10-3 К-1, Сг = 4200 Дж/ (кг-К), От = 1,7-10-4 Н / (м- К), Як = Нк = 5 см. Газ - воздух, температура и давление на удаленном расстоянии от подложки тй = 293 К, рй = 103 Ра, коэффициенты переноса в газе определяются по формулам, полученным из распределения Максвелла. Температура подложки т, = 308 К. В момент - = 0 температура жидкости и газа одинакова и равна т , температура подложки - т,, скорости равны нулю.
Обсуждение результатов
Разность температур между подложкой и газовой средой приводит к неоднородному распределению температур на поверхности капли. Связанный с температурой градиент поверхностного натяжения является движущей силой конвекции Марангони. Это явление существенно для небольших капель, так как в малых объемах гравитационные эффекты пренебрежимо малы. Обычно в численных расчетах конвективный поток тепла от поверхности жидкости описывается посредством числа Био (Вю1;). Однако в [11] было показано, что это приближение может существенно повлиять на результат из-за неточностей определения самого числа Био. В настоящей работе тепловой поток от поверхности описывается уравнениями (3), (6), что дает возможность получить более точное описание процесса теплообмена. Кроме того, в численных исследованиях динамики капли предполагают, что конвекция Марангони не изменяет форму капли, т.е. профиль капли остается равновесным, какие бы течения в жидкости не возникали [4-6]. Вопрос о возможных перемещениях межфазной границы, насколько нам известно, никем не поднимался, хотя это кажется важным, например, для измерений поверхностного натяжения методом большой лежащей капли.
Как показали расчеты, гравитационная конвекция в жидкости в рассмотренных условиях не развивается. Это связано с малой разностью температур между подложкой и вершиной капли - в указанных выше условиях она составляет ~0,2 К. Однако этой разницы температур достаточно для развития конвекции Марангони.
Рис. 2. Конвекция Марангони в осевом сечении капли: а) начальная стадия; б) установивший вихрь
В поведении жидкости и газа можно выделить две стадии.
1) В жидкости вблизи поверхности образуется несколько малых вихрей. Формирование первого такого вихря начинается вблизи подложки и затем охватывает всю приповерхностную область (рис. 2, а). Скорости в жидкости при этом меняются быстро как по величине, так и по направлению. В это время в газе формируется течение, направленное вниз, в горячую область с малой плотностью, из-за чего плотность в нижней части повышается и появляется обратное тече-
Физика
ние. Возникают быстро затухающие течения воздуха, в которых скорости направлены то вверх, то вниз всюду, кроме области, где они огибают каплю.
2) Течения в газе затухают, причем установившееся давление имеет распределение Больцмана вдоль оси z. Вихри в капле сливаются в один большой вихрь, жидкость в котором движется вдоль поверхности капли вверх и вниз вдоль оси в объеме (рис. 2, б). В отсутствие испарения такой вихрь существует все время, пока поддерживается постоянной разность температур между подложкой и газовой средой. Такое же поведение жидкости в капле наблюдалось в [б].
Как выяснилось из расчетов, интенсивное движение внутри капли может нарушить ее равновесную форму. На рис. 3 показана временная зависимость отклонения высоты капли от равновесной. Как видно из рисунка, эта зависимость представляет собой затухающие колебания, причем сравнение с результатами [7] указывает на возбуждение основной моды колебаний лежащей капли. Таким образом, поверхность лежащей капли при нарушении термодинамического равновесия совершает затухающие механические колебания. Амплитуда колебаний зависит от разности температур подложка-газ и может составлять несколько микронов. Колебания быстро затухают, и уже по истечении 0,1 с после начала процесса форму каплю можно считать равновесной.
Заключение
Таким образом, в работе проведено численное исследование динамики поведения лежащей капли и проанализирована эволюция формы капли при нарушении термодинамического равновесия в системе подложка-капля-газ. Неравновесное состояние таково, что температура подложки поддерживается большей температуры газа в удалении от капли. Внутри капли, лежащей на нагретой подложке, возникает вихрь конвекции Марангони, скорости в котором направлены вверх в узкой области вдоль поверхности и вниз во всем остальном объеме.
Численные расчеты показали, что на начальной стадии процесса возникают свободные колебания поверхности капли, частота которых совпадает с частотой основной моды. После затухания колебаний профиль капли остается равновесным, определяемым уравнением Лапласа. Длительность переходного интервала зависит от вязкости жидкости и для капли воды массой 0,525 мг составляет 0,1 с. Интервал установления равновесного профиля этой капли, таким образом, мал, и не может повлиять на точность измерения поверхностного натяжения методом лежащей капли.
Литература
1. Hu H. Evaporation of a sessile droplet on a substrate / H. Hu, R.G. Larson // J. Phys.Chem. B. -2002. - Vol. 10б. - P. 1334-1344.
2. Sefiane K. On hydrothermal waves observed during evaporation of sessile droplets / K. Sefiane, A. Steinchen, R. Moffat // Colloids and Surfaces A: Physicochemical and Engineering Aspects. - 2010. V. 3б5. - P. 95-10S.
3. Analysis of evaporating droplets using ellipsoidal cap geometry / H.Y. Erbil, G. McHale, S.M. Rowan, V.I. Newton // J. Adhes. Sci. Technol. - 1999. - Vol. 13. - P. 1375-13S0.
4. A mathematical model for the evaporation of a thin sessile liquid droplet: Comparison between
experiment and theory / G.J. Dunn, S.K. Wilson, B.R. Duffy et al. // Colloids and Surfaces A: Physico-
chem. Eng. Aspects. - 200S. - V. 323. - Issues 1-3. - P. 50-55.
5. Experimental and numerical investigation of the evaporation into air of a drop on a heated surface / R. Mollaret, K. Sefiane, J.R.E. Christy, D. Veyret // Chem. Eng. Res. Des. - 2004. - V. S2(A4). -P.471-4S0.
6. Evaporation and fluid dynamics of a sessile drop of capillary size / Yu. Barash, T.P. Bigioni,
V.M. Vinokur, L.N. Shchur // Phys. Rev. E. - 2009. - V. 79. - P. 04б301.
7. Коренченко А.Е. Численное исследование свободных колебаний лежащей капли / А.Е. Ко-ренченко, А.Ж. Илимбаева, В.П. Бескачко // Вестник ЮУрГУ. Серия «Математика. Механика. Физика». - 2011. - Вып. 4. - № 10(227). - С. 72-7б.
Д//(мкм)
Рис. 3. Временная зависимость отклонения высоты капли от ее значения в равновесии
Коренченко А.Е., Илимбаева А.Ж., Бескачко В.П.
8. Лабунцов, Д.А. Механика двухфазных систем / Д.А. Лабунцов, В.В. Ягов. - М.: Изд. МЭИ, 2000. - 374 с.
9. Korenchenko A.E. Numerical simulation of oscillations of a drop, bounded by two rigid planes / Korenchenko A.E., Beskachko V.P. // Meas. Sci. Technol. - 2007. - V. 18. - P. 3359-3364.
10. Белоцерковский О.М. Численное моделирование в механике сплошных сред / О.М. Белоцерковский. - М.: Физматлит, 1994. - 448 с.
11. Murisic N. Modeling evaporation of sessile drops with moving contact lines / N. Murisic, L. Kondic // Phys. Rev. E. - 2008. - V. 78. - P. 065301(R).
Поступила в редакцию 13 октября 2012 г.
NUMERICAL STUDY OF A SESSILE DROP ON A HEATED SUBSTRATE
12 3
A.E. Korenchenko , A.J. Ilimbaeva , V.P. Beskachko
The dynamic behavior of a system consisting of sessile drop, heated substrate and gas is investigated by numerical methods. The mode with constant contact angle is considered. The solid surface is assumed to be planar, chemically inert and isothermal. The gas is considered to be a non-condensing ideal gas. The thermal and solutal buoyancy effects in gas are taken into account. Buoyancy and Maran-goni convection are considered in the liquid phase. The movement of liquid-gas interface is analyzed in terms of its impact on the accuracy of sessile drop method.
Keywords: surface tension, sessile droplet, computer simulation.
References
1. Hu H., Larson R.G. Evaporation of a sessile droplet on a substrate. J. Phys.Chem. B. 2002. Vol. 106. pp. 1334-1344.
2. Sefiane K., Steinchen A., Moffat R. On hydrothermal waves observed during evaporation of sessile droplets. Colloids and Surfaces A: Physicochemical and Engineering Aspects. 2010. Vol. 365. pp. 95-108.
3. Erbil H.Y., McHale G., Rowan S.M., Newton V.I. Analysis of evaporating droplets using ellipsoidal cap geometry. J. Adhes. Sci. Technol. 1999. Vol. 13. pp. 1375-1380.
4. Dunn G.J., Wilson S.K., Duffy B.R., David S., Sefiane K. A mathematical model for the evaporation of a thin sessile liquid droplet: Comparison between experiment and theory. Colloids and Surfaces A: Physicochem. Eng. Aspects. 2008. Vol. 323. Issues 1-3. pp. 50-55. http://dx.doi.org/10.1016/-j .colsurfa.2007.09.031
5. Mollaret R., Sefiane K., Christy J.R.E., Veyret D. Experimental and numerical investigation of the evaporation into air of a drop on a heated surface. Chem. Eng. Res. Des. 2004. Vol. 82 (A4). pp. 471-480.
6. Barash Yu., Bigioni T.P., Vinokur V.M., Shchur L.N. Evaporation and fluid dynamics of a sessile drop of capillary size. Phys. Rev. E. 2009. Vol. 79. pp. 046301.
7. Korenchenko A.E., Ilimbaeva A.Zh., Beskachko V.P. Chislennoe issledovanie svobodnykh kolebanii lezhashchei kapli (Numerical study of free vibrations of a sessile drop). Vestnik YuUrGU. serya «Matematika. Mekhanika. Fizika». 2011. Issue 4. no. 10(227). pp. 72-76. (in Russ.).
8. Labuntsov D.A, Iagov V.V. Mekhanika dvukhfaznykh sistem (Mechanics of two-phase systems). Moscow: Izd. MEI, 2000. 374 p. (in Russ.).
9. Korenchenko A.E., Beskachko V.P. Numerical simulation of oscillations of a drop, bounded by two rigid planes. Meas. Sci. Technol. 2007. Vol. 18. pp. 3359-3364.
10. Belotserkovskii O.M. Chislennoe modelirovanie v mekhanike sploshnykh sred (Numerical modeling in continuum mechanics). Moscow: Fizmatlit, 1994. 448 p.
11. Murisic N., Kondic L. Modeling evaporation of sessile drops with moving contact lines. Phys. Rev. E. 2008. Vol.78. p. 065301(R).
1 Korenchenko Anna Evgenievna is Dr. Sc. (Physics and Mathematics), senior researcher, Institute of Metallurgy of the Ural Branch of Russian Academy of Sciences.
E-mail: korenchenko@physics.susu.ac.ru
2 Ilimbaeva Aigul Jepspaevna is Student, Physics Department, South Ural State University.
3 Beskachko Valeriy Petrovich is Dr. Sc. (Physics and Mathematics), Professor, General and Theoretical Physics Department, South Ural State
