Влияние термокапиллярного эффекта на скорость распространения фронта тонкого слоя испаряющейся полярной жидкости Текст научной статьи по специальности «Физика»

ВЛИЯНИЕ ТЕРМОКАПИЛЛЯРНОГО ЭФФЕКТА НА СКОРОСТЬ РАСПРОСТРАНЕНИЯ ФРОНТА ТОНКОГО СЛОЯ ИСПАРЯЮЩЕЙСЯ ПОЛЯРНОЙ ЖИДКОСТИ

В.Ю. Гордеева, А.В. Люшнин

Пермский государственный педагогический университет, 614990, Пермь, Сибирская, 24

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

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

ВВЕДЕНИЕ

Макроскопические тонкие жидкие пленки, порядка 10-100 ангстрем, играют важную роль в таких областях природы, как биофизика, физика, химия, а также находят свое применение в различных технологических приложениях [1-3].

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

© Гордеева В.Ю., Люшнин А.В., 2011

При испарении тонкого слоя полярной жидкости (например, воды) при некотором значении толщины первоначально плоская поверхность становится неустойчивой, и появляются две характерные толщины тонкого слоя жидкости и к2, (Ь1>Ь1) [5], рис. 1. В

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

Рис. 1. Зависимость химического потенциала от толщины слоя к для полярной жидкости. Пунктирными линиями определены два устойчивых значения толщины слоя

Характерной особенностью полярных жидкостей является присутствие как длиннодействующих ван-дер-ваальсовых, межмолеку-лярных сил, так и короткодействующих полярных сил, которые ассоциируются с двойным электрическим слоем. Благодаря такому взаимодействию химический потенциал жидкости ц = dg/ ёк является немонотонной функцией. В этом случае тонкий слой поляризующейся жидкости при испарении может иметь два устойчивых значения толщины слоя и . а также одно неустойчивое значение толщины слоя, расположенное между ними. Эго состояние представляет собой нестабильный фронт, благодаря динамике которого формируется пальцеобразная неустойчивость.

Когда неравномерно нагретая жидкость имеет свободную границу и находится в силовом поле, на её движение существенно влияют зависимость коэффициента поверхностного натяжения от температуры и порождаемый ею термокапиллярный эффект [6]. Эго особенно проявляется в тонком слое жидкости, где термокапиллярная сила (эффект Марангони) может оказаться главным движущим фактором.

1. ПОСТАНОВКА ЗАДАЧИ

Рассматривается тонкий горизонтальный слой полярной жидкости (например, воды) толщиной к, заключенный между твердой подложкой и паром соответствующей жидкости. Жидкость может испаряться в парообразную среду. Твердая граница имеет температуру 0О, а температуру пара обозначим через 0_, причем 0О > 0_ . Предполагается, что имеет место линейная зависимость коэффициента поверхностного натяжения от температуры <7 = <Т0 - г9(0- 0О), здесь г? - температурный коэффициент поверхностного натяжения.

Эволюция тонкого слоя испаряющейся жидкости описана в работах [2, 7] и выглядит следующим образом:

Э/г 1 Э<т

— =-------------------------® ос„ к.У

Ж 2 " ,к ,к

/г2У/г

(ка+ааИ2)

(1.1)

+У—V е [‘'«(«-аУА-ди!

3 п сИг Р Ж

где Г) - вязкость жидкости, <7 - коэффициент поверхностного натяжения границы жидкость-пар, р - плотность жидкости, Цу -

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

Свободная энергия межмолекулярного взаимодействия жидкости, как функция толщины слоя, согласно [5] имеет следующий вид:

где = -Л/(12лг) , А - постоянная Хамакера, <і0 - радиус моле-

кулярного взаимодействия и /0 - дебаевская длина. Значение параметра 8'' определяется адсорбцией ионов из подложки в жидкость. Для воды А < О и Л'7” <0 [2]. Проведем процесс обезразмеривания в следующем виде: х —> /і х, І^хрі, Н —> д Н - где в качестве единиц измерения выбираются следующие соотношения:

Тогда уравнение (1.1) в безразмерной форме будет выглядеть следующим образом:

В уравнении (1.2) используются следующие безразмерные коэффициенты:

(1.2)

где М > 0 и В > 0 являются безразмерными числами Марангони и Био. Параметры О > 0 , % > 0 характеризуют уровень испарения, химический потенциал пара и дебаевскую длину, соответственно.

2. РЕЗУЛЬТАТЫ ЧИСЛЕННОГО ЭКСПЕРИМЕНТА

Уравнение (1.2) решалось численно с использованием пятиточечной конечно-разностной схемы методом Кранка-Николсена, граничные условия предполагались одинаковыми И = кххх = 0 на

обоих концах отрезка, для которого проводились вычисления. В работе [7] было показано, что при малых значениях параметра О может существовать пальцеобразная неустойчивость, а с увеличением значения параметра испарения фронт движения стабилизируется. В работе рассматривается влияние термокапиллярного эффекта на динамику распространения тонкого слоя жидкости, как при малых значениях параметра испарения, где существует область неустойчивости, так и для больших значений безразмерного параметра испарения О . Решение уравнения (1.2) описывает распространение тонкого слоя жидкости в диапазоне значений двух равновесных значений толщины слоя и /?2 .

3. АНАЛИЗ РЕЗУЛЬТАТОВ

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

Численное решение уравнения (1.2), представленное на рис. 2, а и б, показывает динамику движения фронта тонкого слоя жидкости с постоянной скоростью, которая может быть вычислена непосредственно из компьютерных данных траектории движения профиля. На рис. 3 представлена зависимость скорости распространения тонкого слоя жидкости от величины испарения для различных значений безразмерного числа Марангони. С увеличением числа М возрастает скорость движения фронта. При малых значениях термокапиллярного параметра М , как и в работе [7], при О —> 0 скорость движения фронта также стремится к нулю. При больших зна-

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

а б

Рис. 2. Профили распространения фронта при /= 1.085, 0 = 0.07, $ = 0.003, В = 0.001 для М= 0.1 (а) и 25 (б). Штриховой линией обозначено начальное распределение

Рис. 3. Зависимость скорости распространения II от параметра О. для В = 0.001: дляМ = 0.1 (штриховая кривая) и25 (сплошная)

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

БИБЛИОГРАФИЧЕСКИЕ ССЫЛКИ

1. Kistler S.F. and Schweizer P. Liquid Film Coating-Scientific Principles and Their Techological Implications. Kluwer, 1997. 476 p.

2. Oron A., Davis S.H., and Bankoff S.G. Long-scale evolution of thin liquid films // Rev. Mod. Phys. 1997. Vol. 69. P. 931-980.

3. Sharma A., Khanna R. Pattern formation in unstable thin liquid films//Phys. Rev. Lett. 1998. Vol. 81. P. 3463-3466.

4. Troian S.M., Herbolzheimer E., Safran S.A., and Joanny J.F. Fingering instabilities of driven spreading films // Europhys. Lett. 1989. Vol. 10. P. 25-30.

5. Samid-Merzel N., Lipson S.G., and Tannhauser D.S. Pattern formation in drying water films // Phys. Rev E. 1998. Vol. 57. P. 2906-2913.

6. Копбосынов Б.К., Пухначев В.В. Термокапиллярное движение в тонком слое жидкости // Гидромеханика и процессы переноса в невесомости. УНЦ АН СССР. 1983. С. 116-125.

7. Lyushnin А. V., Golovin А.А., and Pismen L.M. Fingering instability of thin evaporation film // Phys. Rev. E. 2002. Vol. 65. P. 021602-021608.

INFLUENCE OF THERMOCAPILLARY EFFECT ON THE PROPAGATION VELOCITY OF THE FRONT EVAPORATION THIN FILM OF A POLAR LIDUID

V.Y. GORDEEVA, A.V. LYUSHNIN

Abstract. This paper considers the dynamics of the evaporation process of a thin film of a polar liquid with a

free surface, is located on a solid substrate. At the free border of the liquid-vapor takes place thermocapillary instability. The chemical potential of a thin liquid film is non-monotonic function and the sum of the interaction of van der Waals forces and the electric double layer. The influence of Marangoni effect on the velocity front evaporating liquid film, are studied.

Key words', the interphase surface, thin liquid film, chemical potential, thermocapillary instability.