∂n = DÑ2 n , |
(8.11) |
¶t |
|
где D– коэффициент диффузии, n(x,y,z,t) – концентрация растворителя в газовой фазе, со следующими начально-краевыми условиями:
на поверхности капли n(0, t) = ns , t > 0, в начальный момент (t = 0)
везде вне капли концентрация определяется асимптотическим значением концентрации паров в атмосфере (для водного раствора –
относительной влажностью воздуха n(x, y, z,0) = n0 ). В сферических координатах при отсутствии зависимости от углов получаем
|
¶n |
|
¶ |
|
n |
2 ¶n |
|
|
|
|
|
|
2 |
|
|
|
|
уравнение |
¶t |
= D |
¶r |
2 + |
|
, |
решение которого дается форму- |
|
|
||||||||
|
|
|
r ¶r |
|
||||
лой
|
|
|
(ns - n0 )R |
r - R |
|
|||
n - n |
0 |
= |
|
erf |
|
|
, |
(8.12) |
|
|
|
||||||
|
|
r |
|
|
|
|
||
|
|
|
2 |
Dt |
|
|||
где
|
2 |
z |
|
erf (z) = |
∫ exp(-x2 )dx |
||
p |
|||
|
0 |
||
R – радиус капли. |
|
|
Из этих формул получаем выражение для интегрального потока,
определяемого законом Фика: |
|
|
|
|
|
|
|
|
|
dn |
|
|
|
|
R |
|
|
|
|
|
|
|
||||
Q = -4pr 2 D |
|
|
» Q |
1 |
+ |
|
, |
(8.13) |
|
|
|||||||
|
|
|
0 |
|
|
|
|
|
|
dr |
|
r =R |
|
|
pDt |
|
|
|
|
|
|
|||||
где Q0 = -4pRD(ns - n0 ) соответствует квазистационарному про-
цессу испарения.
Малость второго члена в скобках в формуле (8.13) по сравнению с единицей является условием справедливости приближения квазистационарности. Для водяных капель, испаряющихся при 20° в сухом воздухе, как показывает подстановка характерных величин, это условие выполняется с высокой точностью [42].
131
Отступление от квазистационарности может иметь место для капель с большим давлением пара (очень быстро испаряющихся жидкостей), но в этом случае нельзя пренебрегать понижением температуры поверхности капель, т. е. следует рассматривать одновременно и концентрационную и температурную релаксации системы. При этом в систему исходных уравнений должно быть включено уравнение теплопроводности.
Скорость изменения объема капли описывается формулой
dV = − Q
, (8.14)
dt nL
где nL – концентрация молекул в жидкой фазе. Отсюда для водной
капли при нормальных условиях (когда справедливо условие квазистационарности и можно пренебречь охлаждением при испарении) уменьшение радиуса со временем дается выражением
R(t) = R 2 |
− 2Dt |
ns − n0 |
, |
(8.15) |
|
||||
0 |
|
nL |
|
|
|
|
|
|
где R0 – начальный радиус капли. В координатах «площадь капли – время», как следует из (8.15), получается прямая линия, причем полное время испарения капли
t = |
R 2 n |
L |
. |
(8.16) |
|
|
|||||
2Da (nS |
− n0 ) |
||||
|
|
|
Выражение (8.16) может рассматриваться и как оценка времени испарения для несферических капель (например, капель на подложке, к описанию которых мы переходим).
8.2.2. Капля на подложке
В этом случае задача с математической точки зрения ставится и решается аналогично рассмотренному случаю сферической капли, только с краевыми и начальными условиями, отвечающими контакту с подложкой (рис. 8.1). Для капель, высота которых значительно меньше капиллярной постоянной, роль гравитационных сил мала
132
(критерий Бонда). В этом случае (речь идет о каплях воды высотой около 1 мм и меньше), используя уравнение Лапласа, можно показать, что форма капли с высокой точностью описывается шаровым сегментом с краевым углом ϕ, который определяется соотношением поверхностных энергий на межфазных поверхностях: верхней выпуклой и нижней – плоской, а также подложкой не занятой жидко-
стью [43].
Хотя при таких предположениях форма капли на подложке довольно проста, имеет место существенное осложнение, отличающее данный случай от сферической капли: необходимо учитывать физикохимию взаимодействия раствора с подложкой.
Данное взаимодействие многофакторно и, вообще говоря, нелинейно по вкладу отдельных составляющих (в нем играют роль чистота и шероховатость подложки, наличие примесей растворенных веществ, каждое из которых имеет свою специфику при высыхании и взаимодействии с подложкой, а также ионная сила растворителя, внешние условия, температурный режим подложки и т.п.); так что универсальный способ учета этих факторов для любого типа капельных систем в единой модели с практической точки зрения едва ли возможен. Но рассмотрение некоторых предельных модельных систем, где учтены только наиболее универсальные факторы, позволяет описать основные особенности процессов.
а |
б |
Рис. 8.1. а – капля воды капиллярного размера на полистирольной подложке (диаметр основания около 3 мм). б – геометрия капли на подложке: краевой угол ϕ, радиус среза R, плотность потока испарения J(r,t), профиль высоты h(r,t) (справа)
[47]
133
В последние 15 лет, ввиду осознания важности данной задачи на новом этапе развития фундаментальных наук и технологий, появилось большое число работ по гидрогазодинамике испаряющейся капли на плоской подложке [44–51].
Прежде всего, отметим работы [44–45], в которых были заложены основы понимания процессов в испаряющейся капле миллиметровых (капиллярных) размеров; была осознана важность учета поведения трехфазной границы (между раствором, подложкой и воздушной средой) – контактной линии для описания гидродинамики процессов внутри капли. В зависимости от свойств подложки и раствора, а также размера капли, определяющего относительную величину капиллярных сил, существуют два основных сценария перемещения контактной линии.
По одному из них контактная линия в процессе высыхания капли может стоять на месте или двигаться очень медленно, так что краевой угол между раствором и подложкой уменьшается, делается неравновесным. По другому сценарию движение контактной линии сохраняет постоянство краевого угла, обеспечивая равновесную форму капли.
Первый сценарий реализуется, если капиллярные силы, подтягивающие контактную линию, слабы по сравнению с силами, удерживающими ее. Это происходит в случае хорошей смачиваемости подложки раствором, например капли воды на стекле, когда высота капли много меньше ее радиуса, а капиллярные силы относительно слабы. Удержанию линии также способствуют неровности подложки, а медленному движению – быстрота испарения. Второй сценарий характерен для случая слабого смачивания подложки раствором (капля на гидрофобной поверхности) – при этом капиллярные силы сравнительно велики, а подложка не удерживает контактную линию.
Согласно Дигану, когда контактная линия удерживается подложкой, имеет место эффект удержания или прилипания (pinning) контактной линии, что приводит к установлению радиальных потоков компенсационной природы, которые перемещают растворенные в капле частицы от центра к периферии, образуя характерное кольце-
134
вое утолщение – избыток твердой фазы на периферии капли – эффект капли кофе на скатерти (рис. 8.2).
Рис. 8.2. Удержание силами сцепления контактной линии капли с подложкой приводит к возникновению радиальных потоков компенсационной природы
Объяснение этого эффекта связано с тем, что для обеспечения постоянства площади контакта при отводе жидкости за счет испарения (рис.8.2 a) необходим подвод раствора к периферии (рис. 8.2 б). Если в капле растворены коллоидные частицы, гидродинамический поток увлекает их к периферии, формируя характерную кольцеобразную структуру.
Противоположный пример, в котором радиальные потоки в чистом растворителе не возникают – полусферическая капля, которая при высыхании сохраняет свою форму. Обычно реализуется некоторый промежуточный случай, когда контактная линия хотя и смещается, но не обеспечивает сохранения контактного угла.
Кроме рассмотренных потоков компенсационной природы в высыхающей капле могут образовываться потоки конвективной природы, обусловленные перепадами температуры по поверхности (потоки Марангони) или в объеме (потоки Рэлея–Бенара ).
В первом случае потоки обусловлены, как правило, зависимостью поверхностного натяжения от температуры. Критерий Марангони имеет вид:
135