Автомодельность задачи тепловой конвекции, осредненной по тонкому слою Текст научной статьи по специальности «Математика»

МЕХАНИКА MECHANICS

УДК 004.942

DOI 10.12737/22161

Автомодельность задачи тепловой конвекции, осредненной по тонкому слою* Л. В. Сахарова*

Ростовский государственный экономический университет, г. Ростов-на-Дону, Российская Федерация Self-similarity problem of thermal convection averaged over a thin layer *** L. V. Sakharova*

Rostov State University of Economics, Rostov-on-Don, Russian Federation

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

Three types of self-simulated replacements for the problem of thermal convection averaged over a thin layer of the vaporizing liquid are presented. It is a model of the drying non-viscous extended droplet specified by the non-thermal diffusivity. For the construction of self-simulated solutions, a transition to the Riemann invariants is performed. Self-simulated solutions are functions of time and position determining the drop height, the mass-transfer rate and the heat flow averaged over the drop thickness. The found self-simulated solutions are classified on the basis of the behavior of the function that describes the drop height under the evaporation-condensation. The domains of applicability of various self-simulated solutions to the simulation of different situations of drying drops and films are identified.

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

Keywords: mathematical model, self-simulated solutions, drop, evaporation- condensation.

Введение. Математическое моделирование процессов тепло- и массопереноса при испарении-конденсации капель жидкостей применимо в настоящее время во множестве современных технологий: медицинской диагностике [1], фармакологических исследованиях [2], кристаллографии белка [3], для растягивания ДНК и РНК [4], полиграфии [5], создании структурированных поверхностей [6], производстве наноструктур [7], микроматриц, в том числе однокристальных лабораторий (labs-on-a-chip) [8]. В связи с этим в последние годы разработано большое количество моделей, описывающих высыхание капель, либо пленок жидкостей. Для математического моделирования процессов в испаряющейся капле применяются различные подходы: в первую очередь, те или иные следствия системы уравнений Навье-Стокса, а также метод осреднения трехмерных задач по высоте пленки, либо капли испаряющейся жидкости.

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

* Работа выполнена в рамках инициативной НИР.

E-mail: sakharova@mail.ru

The research is done within the frame of the independent R&D.

«

S X а

X (L)

В работе [13] разработан универсальный математический аппарат для исследования процессов тепло- и массопереноса в испаряющихся каплях и жидкостях. Введен в рассмотрение оператор осреднения, на основании свойств которого осуществлено осреднение и замыкание приближения Обербека-Буссинеска, наиболее часто используемого для описания гравитационной конвекции в жидкости. Построен ряд математически корректных моделей высокого уровня общности, среди которых наиболее существенной является задача тепловой конвекции, осредненная по тонкому слою испаряющейся жидкости:

И + СИУ(И 2 л) = -К0Ф; (1)

2 4

л*--НЯ&УЯ + — ЪяЧя = 0; (2)

* 3 3

2 4

Ф* - — йфСтта + — клУф = 0; (3)

с{ + клУс = Б Ас. (4)

Здесь: И — переменная во времени высота капли, зависящая от координат; л, ф и с — осредненные по толщине капли скорость массопереноса, тепловой поток и концентрация твердой примеси в капле. Величина У0 — безразмерный параметр. Система (1)-(4), описывающая процесс испарения в общем случае, должна быть дополнена начально-краевыми условиями, определяемыми из физических допущений для каждого конкретного случая. Начальные условия могут определять начальную форму либо высоту капли, а также осредненные значения скорости массопереноса, теплового потока и концентрации примеси в начальный момент времени. Выбор краевых условий существенно зависит от физических предположений, положенных в основу конкретной модели. В зависимости от степени гладкости основания, может иметь место пиннинг границы трехфазного контакта (ее фиксация на шероховатом основании), либо депиннинг (ее отрыв на гладком основании — например, слюде, тефлоне). Угол смачивания поверхности будет острым для смачиваемой (гидрофильной) поверхности и тупым — для несмачиваемой (гидрофобной); кроме того, во втором случае изначально будет иметь место «опрокидывание» профиля капли, что существенно усложняет решение задачи. В случае, когда капля примыкает к жесткой стенке, на ней также необходимо задать угол смачивания или (и) условия теплоизолированности и непроницаемости для потока жидкости.

Как отмечается в работе [13], одним из вариантов упрощения процедуры получения решений является рассмотрение пространственно-одномерных моделей — для вращательно-симметричной капли или «протяженной капли», бесконечно вытянутой в одном из горизонтальных направлений. Оба случая соответствуют общепринятым упрощениям при моделировании капель и тонких пленок. Вращательно-симметричная модель капли отвечает случаю небольших по размеру капель (то есть радиус капли соразмерим с ее тощиной), принимающих под влиянием сил поверхностного натяжения осесимметричную форму. В этом случае в уравнениях модели отсутствует зависимость функций от полярного угла и размерность задачи понижается. Это позволяет построить аналитические решения

^ задачи для осевого сечения капли. Модель «протяженной» капли соответствует случаю тонкой пленки, обладающей относительно большой протяженностью и однородностью в одном из направлений (например, в направлении оси

и

о Оу ). Предполагается, что все физические и геометрические характеристики капли зависят только от переменной х, а

тз

г;^ решения рассматриваются в сечении капли, перпендикулярном оси Оу . Следует отметить, что все краевые эффекты

1л

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

£

Л В работе осуществлено решение задачи (1)-(4) для вращательно-симметричной капли методом характеристик

в сочетании с численными методами. Установлен ряд интересных эффектов, в том числе зависимость от времени формы профиля поверхности капли в окрестности точки трехфазного контакта. Установлено, что в начальный момент

времени форма профиля близка к линейной; однако с течением времени наблюдается ее искривление, а в некоторый фиксированный момент времени происходит «опрокидывание» профиля капли. В работе [13] показано также, что сингулярности решений, возникающих при И ^ 0 для осредненных решений, являются мнимыми и устраняются предельным переходом.

Соответственно, актуальной проблемой является получение аналитических решений для задачи (1)-(4), в первую очередь, ее автомодельных решений [14-15]. В соответствии с [14], явление рассматривается как автомодельное, если распределения его характеристик в различные моменты времени получаются друг из друга преобразованием подобия. Обнаружение автомодельности упрощает вычисления и представление характеристик явления, развивающегося во времени.

Целью настоящей работы является построение автомодельных решений задачи тепловой конвекции и их исследования с точки зрения применимости к моделированию различных ситуаций испарения-конденсации жидкости. Для построения автомодельных решений в работе выполнен переход к инвариантам Римана. Получены определяющие соотношения, позволившие построить три типа автомодельных замен, понижающих размерность задачи и сводящих ее к системе обыкновенных дифференциальных уравнений (ОДУ). Для каждой из систем ОДУ найдены различные аналитические решения, определяющие высоту капли И , осреденные по толщине капли скорость массопереноса 5 и тепловой поток ф как функции времени / и координаты х. Классификация найденных автомодельных решений

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

Задача для протяженной волны в инвариантах Римана. Для протяженной капли, то есть бесконечно вытянутой в одном из горизонтальных направлений, задача тепловой конвекции (1)-(4), принимает одномерный вид:

И( + 2И^Их + И2 5х = -У0ф; (5)

2

5/ + 3 =0; (6)

2 4

ф/ - 3 ИФ5х + 3 И5Фх = 0 (7)

Осуществим переход к инвариантам Римана Я , Р , а , связанным со старыми переменными соотношениями:

И = ЯР3/4; 5 = Р; ф = QP. (8)

В результате система (5)-(7) принимает следующий вид:

Я/ + 2ЯР1'4 Ях = -Уо PQ; (9)

2

р +3 ЯР1/4Рх =0; (10)

а, + 4 ЯР1^ =0. (11)

3 и

— — К

Для построения автомодельных замен выполним переход к функциям Я , Р , а : ев

_____ к

Я = аЯ, Р = РР, а = уа, / = 8/, х = е х, (12) ^

где а, Р, у, 8, е — параметры. При подстановке соотношений (12) в систему (9)-(11) она приводится

кследующему виду: 19

а- а2 В--1/4- --

-я* + 2я • р • Ях = -К0рур • б; (13)

8 е

В- 2 аВ5/4--1/4 -

^Р* +--— Я • Р • Рх =0; (14)

8 3 е

у — 4 аВ1/4у--1/4 —

I + 4 а^ Я • Р • ^ =0. (15)

8 3 е

Инвариантность системы (13)—(15) по отношению к (9)—(11) обеспечивается выполнением равенств:

а = ОГ= Ву; В = аВ5/4; у = а^у, (16)

8 е 8 е 8 е

На основании (16) получим два определяющих соотношения:

е

8 = —174, а = Ру8. (17)

ар

Соотношения (17) позволяют получить автомодельные замены, понижающие размерность исходной задачи (5)-(7) и приводящие ее к системам ОДУ. Рассмотрим основные три типа автомодельных замен и соответствующих им автомодельных решений задачи тепловой конвекции.

Автомодельная замена 1. Полагая а = р = 1, из уравнений (17) получим:

8 = е, у = 1. (18)

8

Соотношения (18) указывают на существование автомодельных замен задачи:

х

- = г, б* = и, (19)

*

где и и * — новые переменные. Действительно, применение (19) к (9)—(11) сводит систему дифференциальных уравнений в частных производных к следующей системе ОДУ:

Яг (2ЯР1/4 - 2) = -У0Ри; (20)

Рг |2ЯР1/4 - г | = 0; (21)

- и + и г | — г + 4ЯР1/4 | = 0; (22)

Система (20)-(22) имеет два решения, определяемых уравнением (21). Решение 1.1. Пусть

3

ЯР1/4=- г. (23)

2

Р

2 Подстановка (23) в уравнение (22) приводит к решению вида:

С и = С, г, (24)

О 1

тз

^ С1 — константа интегрирования, а подстановка (24) в (20) позволяет получить решение системы в инвариантах 'й

^ Римана:

>

- ( 81 У/5

Л Я = |С2 - — С1Ког5 , (25)

Р = 81 г4(С2 -81 С^г5Т^, (26)

20 16 I 2 32 1 0 ) К '

о = с,

(27)

Вернемся от инвариантов к обычным переменным, и на основании (25)-(27) получим решение исходной системы (5)-(7):

8 ( _ 5 81 _ _ 5

И т|С2^5 - — с,г0х5

27х \ 16

5 = 81 х4 ( Сг15 - 81 С1¥0х5

16

16

4/5

-4/5

81 „ х5 („ 5 81 54 4/5

Ф = — с, — I С2Г--с,у0х

16 1 ¿2 I 2 16 1 0

(28)

(29)

(30)

Функция И , определяющая тощину капли, ведет себя по-разному в зависимости от знаков констант С1, С2. Пусть С1; С2 > 0 (¥0 > 0). Тогда из (28) видно, что функция И обращается в ноль при

хк

2 3С2

= — 5

(31)

3 \

Точка (31) соответствует краю капли, причем это значение увеличивается с возрастанием /. Из формулы (28) следует, что при С1; С2 >0 , если 0 < х < хк, то И'х < 0 , к[ >0 то есть имеет место убывание функции по х и возрастание / (конденсация). Если х > хк , то Их > 0, к[ <0, что означает возрастание функции по х и убывание / (испарение). Таким образом, решение (28)-(30) может быть использовано непосредственно для описания испаряющейся капли, если в качестве начального момента времени положить = 0 > 0 . Предполагается, что левый край капли (стенка, к которой примыкает капля) находится в точке х = х0 > хк. Тогда константы С1, С2 можно определить, например, задавая значения И(х0) = И0, х0) = при / = /0. Если в качестве левого края капли взять точку х = 0, то формулы (28)-(30) определяют конденсацию капли. Однако их нельзя использовать целиком из-за сингулярности. Используем сшивку полученного решения с константой (рис.1):

Рис.1. Поверхность капли, описываемая Решением 1.1, С1, С2 >0 (¥0 > 0 ). хк — точка трехфазного контакта. Левая область — конденсация капли; х = 0 (начало кооординат) соответствует жесткой стенке, к которой примыкает капля. Правая область — испарение; х = х0 — жесткая стенка.

Кривые (1) и (2) соответствуют различным моментам времени < ;

ев И

к

X ев

X

К

h =

5 =

8 27^ \ СЛ5 - 81 С|У° х5 ^4/5

8 I C2t5 -Е ^ х5; 4/5

27 х

81 4 16 хс 1С^5 -11 ^ хс;

81 х 4 Г 16 1 С^5 - § с1у0 х5)" -4/5

0 < х < х„

Хс ^ Х < хк ■

0 < х < х„

Хс ^ Х < Хк ■

(32)

(33)

81^ х5_ I 5 16 1 (2

-Сх"2№5 -16^У)Х

81 с1 х51С2Г5 - 81СУх5

16 1 12 I 2 16 1 0

-4/5

-4/5

0 < х < х„

Хс ^ х < хк ■

(34)

Таким образом, при 0 < х < хс функции к , 5, ф представляют собой константы, зависящие от времени. ^отношения между константами, входящими в (32)-(34) можно получить, задавая к(0,t0) = к0, 5(0,t0) = 50, ф(0, t°) = ф0. Следовательно, Решение 1.1 при С1, С2 >0 , определяемое формулами (32)-(34), в различных областях задания может быть использовано для описания как конденсирующейся, так и испаряющейся капли.

Если же С1, С2 < 0 , то, области при х < хк и х > хк меняются ролями: левая часть отвечает за испарение, а правая — за конденсацию. Если С1 > 0, С2 < 0 (либо С1 > 0, С2 < 0), то подкоренная функция в ноль не

2 9С

обращается. Однако у функции к имеется минимум в точке х<хй = — 5-— t. Можно предположить, что функция

3 V 2У0С1

к в этом случае описывает испарение (конденсацию) пленки жидкости, расположенной неравномерным слоем между двумя твердыми стенками.

Решение 1.2. Рассмотрим второе решение уравнения (21), соответствующее Р7 = 0 :

->4

Р = А

где А — константа интегрирования. Тогда оставшиеся уравнения системы преобразуются к виду:

Я2 (2ЯА - 7) = -УйА4И;

- и + ы21 - г + —ЯБ1 | = 0.

(35)

(36)

(37)

Система (36)-(37) допускает построение асимптотического решения в виде рядов. Возвращение к исходным переменным приводит к квадратичной асимптотике автомодельного решение системы (5)-(7):

( Л Л. ГЛ П Л2

к -

Я0 +

2ДУ

= -2У0 А?щ Я )|х + АЯ

Ях =5Я0 + 2У0 В14и0.

5 = А4

ф

= А4 t

^ + _3^ х + 9я^АГ2

4Я0 А t 8

и0 +

У0 А и0

3ЯП

4.. А „2 ^

, t2,

(38)

(39)

(40)

й о тз

"¡3

и

и

О, С -Й

Исходя из требования неотрицательности функции толщины капли к в начальный момент времени, будем считать, что АЯ0 >0 . Предположим для определенности: А >0, Я0 > 0. Парабола (38) имеет вершину в точке

в(-tD1Я0/Я1;Я0 + 2А14У0/Я1), и пересечение с осью ОХ в точке хк = ^Я1

- АЯ ±

Я0 + 2У0 А4

У0 АЧ

Таким образом, с течением времени высота капли, описываемой (38)-(40) возрастает, а площадь ее увеличивается. Следовательно, автомодельное Решение 1.2 соответствует конденсации капли (Рис.2). Неизвестные константы и0 и Я0 могут быть найдены, например, из задания углов смачивания в начальный момент времени

t = t0 > 0 точках х = 0 и х = хк : к'х (х = 0) = tga|, к'х (х = хк) = tga2.

Ф

Рис.2. Поверхность капли, описываемая Решением 1.2, конденсация, ^ > 0, Я0 > 0; х = 0 — жесткая стенка, к которой примыкает капля, х = хк — точка трехфазного контакта.

Автомодельная замена 2. Пусть р = у = 1. Тогда из (17) следует, что

8 = а, е = 82.

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

х

приводящие систему (9)—(11) к виду:

2 = —, Я = V, г

V + 2 V, (- 2 + VP1/4 ) = -У0 РО. Р, | - 22 + 2 VP1/4 | = 0;

О, I- 22 + 4 VP1/4 | = 0.

Из уравнения (44) следует, что

или

а из уравнения (45) — что либо

Р =0,

^Р1/4 = 3,,

О, =0,

VPш = -z.

2

Рассмотрим решения, соответствующие различным комбинациям (46) - (49). Решение 2.1. Пусть выполняются (46) и (48), следовательно,

Р = Р/,

о=а,

где Р1, О1 — константы. Тогда решение системы (43)-(45) имеет вид:

(

х+г2

V Р/а+

Р1

А

(

V

=

1

(

\\

нрх3 + г VoР14О1 +

V

1

8 = Р4

О = ЙР4,

где Р1, О1, У1 — константы. Кривая (52) есть парабола с вершиной в точке

В

'- г2 РК;-*

К = Vo РГО. +1

(41)

(42)

(43)

(44)

(45)

(46)

(47)

(48)

(49)

(50)

(51)

(52)

(53)

(54)

(55)

и пересечением с осью ОХ в точке С (г2 Р16КР001;0). Если, например, Р1К <0 , то вершина расположена в первой координатной четверти, а ветви параболы повернуты влево. Таким образом, в соответствии с (55), в начальный

ев И К X ев X

К

2

1

момент времени г = 0 получаем параболу с вершиной в начале координат: х = Р^^й2. При увеличении времени вершина парболы смещается вправо и вверх. Таким образом, наблюдается опрокидывание профиля капли, а затем высота капли увеличивается. Следовательно, Решение 2.1, определяемое формулами (52)-(54), соответствует конденсации капли на вертикальной стенке, совпадающей с осью ординат (рис. 3).

Рис. 3. Поверхность капли, описываемая Решением 2.1, конденсация капли на вертикальной стенке.

Линии соответствуют различным моментам времени г1 < г2

3,

Решение 2.2. Пусть теперь выполняются (46) и (49). Следовательно, Р = , V =-, где Р1 — константа.

2Р

В этом случае получим второе решение системы (43)-(45):

И =

3х

2Р147'

8 = Р4,

ф:

3х

V Рхг

2 '

(56)

Решение 2.2., определяемое уравнениями (56) соответствует случаю испарения, когда поверхность капли линейна, жидкость уменьшается со временем. В данном решении по смыслу задачи > 0; в противном случае формула (56) описывает «перевернутую» каплю.

81, 4

Решение 2.3. Пусть теперь выполняются (47) и (48), О = 01; Р = ——

приобретает линейную форму:

-+V=

ск 4, 0

3

Решив (57), получим третье решение системы (43)-(45):

1/5

Я = г,-1/4,, , = 5^4 - 27 v0Qíz21/4 Р = "4,

о=а.

После возвращения к исходным переменным (8) в (58)-(59), имеем автомодельное решение исходной задачи:

г 4/5 1/5 / 27

И = г-7- 54/51 м21/2 - 27 VQíx21/4

8 = -

81х4 5 г-4/5{ л .21/2 27 т^ ..21/4

г9/5 5-4'5|л,г21/2-уVoQ,х2

4/5

4/5

й о тз

"¡3

и

и

О, £ Л

^-П 81х с-4/5{ , .21/2 27 21/4

ф = г4/5 5 I л1г VoQ\x

4/5

В этом случае уравнение (43)

(57)

(58)

(59)

(60)

(61) (62) (63)

Функция И , определяемая формулой (61) не имеет сингулярностей ни по х, ни по г. При фиксированном г

функция имеет пересечения с осью ОХ в точках хп =0, и хк =

7 А

ч 4/21

V 27^а у

2

и максимум в точке

28 Л1

675^21 у

ч 4/21

г , (0 < хт < хк).

хт

h(x)

Рис.4. Поверхность капли, описываемая Решением 2.3. хк — точка трехфазного контакта. Левая область — конденсация капли; х = 0 — фиксированный левый край капли; х = хк — движущийся правый край. Правая область — испарение; х = х0 — жесткая стенка, к которой примыкает капля. Кривые (1) и (2) соответствуют различным моментам времени < t2

Таким образом, решение (61)—(63) может быть использовано в области х0 < х < хк для описания конденсирующейся капли в условиях пиннинга ее левого края, а в области хк < х < х0 — испаряющейся капли (х0 — жесткая стенка, к которой примыкает капля) (рис.7). Соотношения между константами, входящими в (61) -(63) можно получить, задавая к(хт,t0) = к0, хт,t0) = 50. Следовательно, Решение 2.3 в различных областях задания может быть использовано для описания как конденсирующейся, так и испаряющейся капли.

Автомодельная замена 3. Пусть а = у = 1. Тогда из (17) следует, что

ß = 1, 8

е = 8

3/4

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

Х P = W/t,

t

3/4

где г и W — новые неизвестные переменные. Применение формул (65) к системе (13)-(15) приводит ее виду:

3

Rz| --j-z + 2RW1/4 | = -V0WQ

-W + Wz| -3Z + 2RW1/4 1 = 0;

Qz|-4z + 3RW1/4 | = 0.

Из (68) следует, что

или

Qz =0,

RW1/4 = - z. 16

Система (66)-(68) имеет два решения, определяемых уравнениями (69) и (70).

Решение 3.1. Рассмотрим вначале случай (70). Система сводится к виду:

3

Я-2 = ;

Решение уравнения (72) имеет вид:

Wz | - 3 z | = W.

W = B1 exp1 - 4 z2

В1 — константа интегрирования. Тогда из уравнений (70), (71) и (73) получим:

(64)

(65)

(66)

(67)

(68)

(69)

(70)

(71)

(72)

сЗ И

к

X а

X к

Я = ехр| 1 г2

16 В

Q = -

27 г

144 В 5/4Уп

1+3г 2) ехр( 3г 2

Вернемся в (74)-(75) от инвариантов к обычным переменным (8), получим:

9х

к =-ехр

16В,

В1

5 = ~у~ехр

(

( л 2 А 4х 2

4х2 "

ф:

27 7

144 В^4^

1 +

2х 2

¡— ехр 3л113 ) I ^t3

' х2 ^

(74)

(75)

(76)

(77)

(78)

Из формул (76)-(78) следует, что функция к возрастает по х, убывает по t; функция 5 убывает по х, убывает по t; функция ф возрастает по х, убывает по t (рис.5). Константа В1 может быть найдена из задания, например, угла

смачивания при х = 0 .

Рис. 5. Поверхность капли, описываемая Решением 3.1, конденсация капли в условиях пиннинга. Линии соответствуют различным моментам времени ^ < /2

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

Решение 3.2. Наконец, на основании (69) получим частное решение задачи:

к = _1 х11513/5,

5 = В

4 х

4/5

В4 t8/5 ' ' 5У0 В t8

где В есть произвольная постоянная. Решение 3.2, определяемое (79), соответствует конденсации капли в условиях пиннинга (рис.5).

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

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

1. Гольбрайх, Е. О. формировании узора трещины в свободно высыхающей пленке водного раствора белка / Е. Гольбрайх, Е. Г. Рапис, С. С. Моисеев // Журнал технической физики, 2003. — Т. 73, вып. 10. — С. 116-121.

2. Способ оценки общетоксического действия лекарственных средств на организм : патент 2232387 Рос. Федерация : G01N33/15, G01N33/49 / А. А. Ющенко, А. Д. Даудова, А. К. Аюпова, Н. Г. Урляпова , С. Н. Шатохина. — № 2002129685/15 ; заявл. 04.11.2002 ; опубл. 10.07.2004. — 7 с.

3. Рапис, Е. Белок и жизнь (самосборка и симметрия наноструктур белка) / Е. Рапис. — Иерусалим; Москва : ЗЛ. Милта-ПКП ГИТ, 2002. — 257 с.

Ф:

63 х

4/5

8/5

(79)

й о

ТЗ

'й

и

и

О, £ -Й

4. Abramchuk S.S., Khokhlov A.R., Iwataki T., Oana H., Yoshikawa K. Direct observation of DNA molecules in a convection flow of a drying droplet // Europhys. Lett. 2001. - Vol. 55. P. 294-300.

5. Harris D. J., Hu H., Conrad J. C., Lewis J. A. Patterning Colloidal Films via Evaporative Lithography // Physical Review Letters. 2007. - Apr. Vol. 98, no. 14. P. 148301.

6. Xu J., Xia J., Hong S. W. et al. Self-Assembly of Gradient Concentric Rings via Solvent Evaporation from a Capillary Bridge // Physical Review Letters. 2006.-Feb. Vol. 96, no. 6. P. 066104.

7. Helseth L. E., Fischer T. M. Particle interactions near the contact line in liquid drops // Physical Review E. 2003.-Oct. Vol. 68, no. 4. P. 042601.

8. Rieger B., van den Doel L. R., van Vliet L. J. Ring formation in nanoliter cups: Quantitative measurements of flow in micromachined wells // Physical Review E. 2003.-Sep. Vol. 68, no. 3. P. 036312.

9. Deegan R.D., Bakajin O., Dupont T.F., Huber G., Nagel S.R., Witten T.A. Contact line deposits in an evaporating drop // Physical Review E. - 2000. - vol. 62. - P. 756-765.

10. Maki K. L., Kumar S. Fast Evaporation of Spreading Droplets of Colloidal Suspensions // Langmuir. 2011. Vol. 27, no. 18. P. 11347- 11363.

11. Widjaja E., Harris M. Particle deposition study during sessile drop evaporation // AIChE J. 2008.-September. Vol. 54, no. 9. P. 2250- 2260.

12. Tarasevich Y. Y., Vodolazskaya I.V., Sakharova L.V. Mathematical modeling of pattern formation caused by drying of colloidal film under a mask // Eur. Phys. J. E. - 2016. - Vol. 39, no. 2.

13. Жуков, М. Ю. Моделирование испарения капли жидкости / М. Ю. Жуков, Е. В. Ширяева, Н. М. Полякова. — Ростов-на-Дону : Издательство Южного федерального университета, 2015. — 208 с.

14. Баренблатт, Г. И. Подобие, автомодельность, промежуточная асимптотика. Теория и приложения к геофизической гидродинамике / Г. И. Баренблатт. — Ленинград : Гидрометиоиздат. — 1982. — 257 с.

15. Рождественский, Б. Л. Системы квазилинейных уравнений и их приложения к газовой динамике / Б. Л. Рождественский, Н. Н. Яненко. — Москва : Наука — 1978.— 687 с.

References

1. Golbraykh, Е.О., Rapis, E.G., Moiseev, S.S. O formirovanii uzora treshchiny v svobodno vysykhayushchey plenke vodnogo rastvora belka. [On the formation of the crack pattern in a free drying film of aqueous protein solution. ] Technical Physics, 2003, vol. 73, iss. 10, pp. 116-121 (in Russian).

2. Yushchenko, A.A., Daudova, A.D., Ayupova, A.K., Urlyapova, N.G., Shatokhina, S.N. Sposob otsenki obshchetoksicheskogo deystviya lekarstvennykh sredstv na organizm: patent 2232387 Ros. Federatsiya : G01N33/15, G01N33/49. [A method of estimating general toxical action of drugs on the body.] Patent RF, no. 2232387, 2004 (in Russian).

3. Rapis, E. Belok i zhizn' (samosborka i simmetriya nanostruktur belka). [Protein and life (self-assembly and symmetry of protein nanostructures.)] Jerusalem; Moscow: ZL. Milta-PKP GIT, 2002, 257 p. (in Russian).

4. Abramchuk, S.S., Khokhlov, A.R., Iwataki, T., Oana, H., Yoshikawa, K. Direct observation of DNA molecules in a convection flow of a drying droplet. Europhys. Lett., 2001, vol. 55, pp. 294-300.

5. Harris, D. J., Hu, H., Conrad, J. C., Lewis, J. A. Patterning Colloidal Films via Evaporative Lithography. Physical Review Letters, 2007, vol. 98, no. 14, p. 148301.

6. Xu J., Xia J., Hong S.W., et al. Self-Assembly of Gradient Concentric Rings via Solvent Evaporation from a Capillary Bridge. Physical Review Letters, 2006, vol. 96, no. 6, p. 066104.

7. Helseth, L. E., Fischer, T. M. Particle interactions near the contact line in liquid drops. Physical Review E, 2003, vol. 68, no. 4, p. 042601.

8. Rieger, B., van den Doel, L. R., van Vliet, L. J. Ring formation in nanoliter cups: Quantitative measurements of flow in micromachined wells. Physical Review E, 2003, vol. 68, no. 3, p. 036312.

9. Deegan, R.D., Bakajin, O., Dupont, T.F., Huber, G., Nagel, S.R., Witten, T.A. Contact line deposits in an й

w

evaporating drop. Physical Review E, 2000, vol. 62,pp. 756-765.

Рн

10. Maki, K. L., Kumar, S. Fast Evaporation of Spreading Droplets of Colloidal Suspensions. Langmuir, 2011, vol. ^ 27, no. 18, pp. 11347- 11363. Ц

11. Widjaja, E., Harris, M. Particle deposition study during sessile drop evaporation. AIChE J., 2008, vol. 54, no. 9, pp. 2250- 2260.

12. Tarasevich, Y. Y., Vodolazskaya, I.V., Sakharova, L.V. Mathematical modeling of pattern formation caused by 27

drying of colloidal film under a mask. Eur. Phys. J. E, 2016, vol. 39, no. 2.

13. Zhukov, М.Y., Shiryaeva, E.V., Polyakova, N.M. Modelirovanie ispareniya kapli zhidkosti. [Simulation of the liquid drop evaporation.] Rostov-on-Don: SFU Press, 2015, 208 p. (in Russian).

14. Barenblatt, G.I. Podobie, avtomodel'nost', promezhutochnaya asimptotika. Teoriya i prilozheniya k geofizicheskoy gidrodinamike. [Similarity, self-similarity, intermediate asymptotics. Theory and applications to geophysical hydrodynamics .] Leningrad: Gidrometioizdat, 1982, 257 a (in Russian).

15. Rozhdenstvenskiy, B.L., Yanenko, N.N. Sistemy kvazilineynykh uravneniy i ikh prilozheniya k gazovoy dinamike. [Systems of quasilinear equations and their applications to gas dynamics.] Moscow: Nauka, 1978, 687 p. (in Russian).

Поступила в редакцию 15.06.2016 Сдана в редакцию 16.06.2016 Запланирована в номер 30.09.2016

Received 15.06.2016 Submitted 16.06.2016 Scheduled in the issue 30.09.2016

й о тз

"¡3

и

(U >

Л С Л