Аналитический метод исследования кинетики процесса растекания капли

Размещено на http://www.allbest.ru/

Кабардино-Балкарский государственный университет им. Х.М. Бербекова, Нальчик, Россия

Аналитический метод исследования кинетики процесса растекания капли

Лесев Вадим Николаевич

к.ф.-м.н., доцент

Введение

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

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

1. Постановка задачи

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

Выбор такого подхода подразумевает, что вблизи линии трехфазного контакта высота капли стремительно уменьшается и становится доступной действию поверхностных сил вплоть до уровня подложки, т.е. ограничимся изучением капель, геометрические характеристики профиля которых в начальный момент времени удовлетворяет условию (рис. 1). В этом случае давление в капле будет определяться равенством [10, 11]:

, (1)

где - давление газа, - изотерма расклинивающего давления.

Рис. 1. К выводу уравнения растекания капли

Отметим, что довольно часто имеет степенную зависимость [8, 12] от :

, (2)

где A - постоянная Гамакера.

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

Принимая во внимание (1) и (2), по аналогии с [9, 10], можем записать

, (3)

где - вязкость, .

Соотношение (3) представляет собой общее уравнение растекания капли по сухой подложке. Отметим, что при , т.е. в случае отсутствия расклинивающего давления, уравнение (3) переходит в известное уравнение растекания [9].

Вводя безразмерные переменные: , , , где нижним индексом “” обозначены характерные значения времени, высоты и абсциссы, запишем (3) в виде:

. (4)

, .

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

Опуская для удобства записи обозначения безразмерных переменных, перепишем (4) в виде:

. (5)

Зададим граничные условия:

, (6)

, (7)

, (8)

где .

Условия (6-8) представляют собой условия сохранения объема, симметрии профиля капли относительно вертикальной оси и обращения в ноль высоты жидкой пленки при ее растекании на линии трехфазного контакта. Здесь необходимо отметить, что при формулировке условия (6) для определения объема было принято соотношение .

2. Построение приближенного аналитического решения задачи

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

Положим . В результате уравнение (5) примет вид:

.

Интегрируя дважды последнее равенство и принимая во внимание (8), получим:

,

где - произвольная функция интегрирования.

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

.

Интегрируя последнее равенство по в границах от до , находим:

.

Пользуясь произвольностью выбора функции , перепишем последнее равенство в виде:

, (9)

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

Далее, удовлетворяя (9) условию (6), будем иметь

,

откуда находим:

. (10)

Подставляя полученное значение для радиуса растекания капли (10) в (9), получим:

. (11)

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

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

Для определения функции подставим выражение (11) в уравнение (5), предварительно определив производные от функции .

Дифференцируя (11), будем иметь:

, (12)

. (13)

Из (5), принимая во внимание соотношения (12), (13), а также выражение (11), получим

(14)

Рассмотрим отдельно выражение, стоящее в квадратных скобках в правой части равенства (14).

Имеем:



.

С учетом полученного равенства, (14) примет вид:

.

Отсюда, в результате элементарных преобразований, будем иметь

, (15)

где .

Раскрывая скобки в правой части последнего равенства, перепишем (15) в виде:

или

. (16)

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

Соотношение (16) при принимает вид:

. (17)

Разделяя переменные в уравнении (17), получим

. (18)

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

. (19)

Нелинейное дифференциальное уравнение (19) будет иметь явные решения лишь при . В этом случае (19), в результате элементарных преобразований примет вид:

. (20)

Интегрируя (20), будем иметь:

или ,

где - постоянная интегрирования.

Для определения положим в последнем равенстве . Принимая высоту безразмерной капли в апексе в начальный момент времени равной единице, с учетом (11) получим .

Таким образом, окончательно устанавливаем, что .

Следовательно, для будем иметь:

. (21)

Из (10) и (21) находим:

, (22)

где, как уже было отмечено выше, .

Для размерных переменных равенство (22) примет вид:

. (23)

Положив в (23) , получим:

. (24)

Соотношение (24) устанавливает связь между радиусом растекания капли в начальный момент времени и нормализованным горизонтальным размером.

Скорости растекания капли для размерных переменных определим из (23). В результате будем иметь:

.

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

.

Считая, что в начальный момент времени значение согласно (24) имеет приращение , а приращение по времени , из последнего соотношения заключаем:

. (25)

Следовательно, выражение (21) для может быть записано в виде:

. (26)

Также могут быть преобразованы выражения для радиуса растекания капли в безразмерных (22):

и размерных переменных (23):

. (27)

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

и размерных переменных:

. (28)

Из (27), (28) следует, что выбор нормализующих параметров , равносилен заданию геометрических параметров размерной капли и в начальный момент времени , определение которых для натурных экспериментов также может быть определено на основе методов математического моделирования [14, 15]. И, наоборот, при априори известных и переход к безразмерным параметрам должен осуществляться по соответствующим характерным значениям длины , высоты и времени .

При этом функция , принимая во внимание (11) и (21), окончательно, для случая безразмерных переменных, примет вид:

.

Заключение

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

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

растекание капля изотерма давление

1. Гребенник И.П. Кинетические особенности растекания расплавов по поверхности тонких пленок // Украинский физический журнал. 1980. Т. 25. № 3. С. 497-503.

2. Brochard F., De Gennes P.G. Spreading laws for liquid polymer droplets: interpretation of the "foot" // J. Physique Lett. 1984. Vol. 45. № 12. P. 597-602.

3. Joanny J.-F. Wetting of a liquid substrate // Physicochem. Hydrodyn. 1987. Vol. 9. № 1. P. 183-196.

4. Старов В.М., Чураев Н.В. Толщина смачивающих пленок на шероховатых подложках // Коллоидный журнал. 1977. Т. 39. № 6. С. 1112-1117.

5. Cazabat A.-M., Cohen Stuart M.A. Dynamics of wetting on smooth and rough surfaces // Progress in Coll. and Polymer Sci. 1987. Vol. 74. P. 69-75.

6. Cox R.G. The spreading of a liquids on a rough solid surface // J. Fluid Mech. 1983. Vol. 131. P. 1-26.

7. Воинов О.В. Динамическая теория смачивания твердых тел вязкими жидкостями под действием сил Ван-дер-Ваальса // ПМТФ. 1994. № 6. С. 69-85.

8. Дерягин Б.В., Чураев Н.В. К вопросу об определении понятия расклинивающего давления // Коллоидный журнал. 1976. Т. 38. № 3. С. 438-448.

9. Старов В.М. О растекании капель нелетучих жидкостей по плоской твердой поверхности // Коллоидный журнал. 1983. Т. 45. № 6. С. 1154-1161.

10. Калинин В.В., Старов В.М. Растекание капель жидкости с учетом действия поверхностных сил // Коллоидный журнал. 1988. Т. 50. № 1. С. 25-32.

11. Старов В.М., Чураев Н.В., Хворостянов А.Г. О форме движущегося мениска в плоских капиллярах // Коллоидный журнал. 1977. Т. 39. № 1. С. 201-205.

12. Дерягин Б.В., Чураев Н.В., Муллер В.М. Поверхностные силы. М.: Наука, 1985. 398 с.

13. Найфе А. Введение в методы возмущений. М.: Мир, 1984. 535 с.

14. Лесев В.Н., Созаев В.А. О новом методе обработки экспериментальных данных для малых капель расплавов // Известия КБГУ. 2011. Т. 1. № 1. С. 3-8.

15. Директор Л.Б., Зайченко В.М., Майков И.Л. Усовершенствованный метод лежащей капли для определения поверхностного натяжения жидкостей // ТВТ. 2010. Т. 48. № 2. С. 193-197.

16. Бородин С.А. Исследование процесса растекания капли жидкости, наносимой на поверхность подложки // Компьютерная оптика. 2005. № 28. С. 66-68.

17. Вавкушевский А.А., Арсланов В.В., Огарев В.Л. Растекание капель полимеров по гладким твердым поверхностям // Коллоидный журнал. 1984. Т. 46. № 6. С. 1076-1081.

18. Лесев В.Н., Бжеумихова О.И. Математическая модель растекания капли под действием гравитационных сил // XLIV Всероссийская конференция по проблемам математики, информатики, физики и химии: Тезисы докладов; РУДН. М., 2008. С. 33-34.

19. Пономарева М.А., Якутенок В.А. Моделирование растекания капли вязкой жидкости в плоской постановке при больших числах бонда // Вестник томского государственного университета. Математика и механика. 2007. № 1. С. 79-83.

Размещено на Allbest.ru

...