Лабораторные работы по физике аэрозолей (Чукин В.В.) / clouds
.pdf® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
Таблица 3.2 − Зависимость поправки на размер капли от ее радиуса (при t=10 °С и P=800 гПа )
r, м |
cR |
1.0E09 |
0.000 |
5.0E09 |
0.001 |
1.0E08 |
0.002 |
5.0E08 |
0.009 |
1.0E07 |
0.017 |
5.0E07 |
0.081 |
1.0E06 |
0.150 |
5.0E06 |
0.469 |
1.0E05 |
0.638 |
5.0E05 |
0.898 |
1.0E04 |
0.946 |
5.0E04 |
0.989 |
1.0E03 |
0.994 |
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
1.0 |
|
|
|
|
|
|
0.9 |
|
|
|
|
|
|
0.8 |
|
|
|
|
|
|
0.7 |
|
|
|
|
|
|
0.6 |
|
|
|
|
|
|
0.5 |
|
|
|
|
|
|
0.4 |
|
|
|
|
|
|
0.3 |
|
|
|
|
|
|
0.2 |
|
|
|
|
|
|
0.1 |
|
|
|
|
|
|
0.0 |
1.0E08 |
1.0E07 |
1.0E06 |
1.0E05 |
1.0E04 |
1.0E03 |
1.0E09 |
Рисунок 3.1 − Зависимость поправки на размер капли от ее радиуса (при t=10 °С и P=800 гПа )
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
При выводе уравнения (3.20) не учитывалось падение капель относительно воздушной среды. Падение капель приводит к эффекту вентиляции, то есть к перемешиванию воздуха вблизи поверхности капли и тем самым к увеличению концентрации молекул водяного пара у поверхности капли, что способствует увеличению скорости роста капли. Этот эффект можно учесть, вводя в уравнение роста капли (3.20) коэффициент вентиляции (называемый также ветровым множителем), равный:
1.00 0.108 |
|
|
|
|
|
|
|
|
|
при Re 2.5; |
|
|
|
|
|
|
|||||||
|
Re 3 Pr* 2 , |
|
|||||||||
cF={0.78 0.308 |
|
3 |
|
, |
|
(3.22) |
|||||
Re |
Pr* |
при Re≥2.5, |
где cF |
− коэффициент вентиляции (ветровой множитель); |
|
|||||
Re |
− число Рейнольдса, равное |
2 U ' r |
; |
|
|
||
|
|
|
|
|
|
|
|
Pr* |
− число Прандтля для потока, равное |
|
|
; |
|
||
|
D |
|
|||||
|
|
|
|
|
П |
|
|
U ' |
− скорость |
частицы |
относительно |
воздушной |
среды, |
равная V −U , м/с;
r− радиус частицы, м;
− динамическая вязкость атмосферы, кг/(м∙с);
DП − коэффициент молекулярной диффузии водяного пара в воздухе, м2/с.
Для капель с радиусом r 50 мкм эффект вентиляции можно не учитывать, то есть считать cF=1 . С увеличением размеров капель эффект
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
вентиляции становится более существенным. Так, для капель с радиусом r≈100 мкм ветровой множитель cF≈1.5 .
Другой аппроксимацией зависимости ветрового множителя от радиуса капли является формула [1]:
cF=1.00 0.24 Re . (3.23)
Результаты расчетов по формулам (3.22) и (3.23) представлены в таблице 3.3 и на рисунке 3.2.
Окончательно уравнение диффузионного роста капель принимает вид:
dn |
=4 r DП N П −N ПВ cR cF . |
(3.24) |
|
d |
|||
|
|
Для перехода от изменения числа молекул в капле dn |
к изменению |
радиуса капли dr необходимо использовать определение |
массы капли |
двумя способами.
С одной стороны, масса капли равна сумме масс всех молекул, из которых состоит капля:
M =mH2 O n . (3.25)
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
Таблица 3.3 − Зависимость ветрового множителя от радиуса капли (при t=10 °С и P=800 гПа )
r, м |
|
cF |
|
(3.22) |
|
(3.23) |
|
|
|
||
1.0E06 |
1.00 |
|
1.00 |
2.0E06 |
1.00 |
|
1.00 |
3.0E06 |
1.00 |
|
1.00 |
4.0E06 |
1.00 |
|
1.01 |
5.0E06 |
1.00 |
|
1.01 |
6.0E06 |
1.00 |
|
1.01 |
7.0E06 |
1.00 |
|
1.02 |
8.0E06 |
1.00 |
|
1.02 |
9.0E06 |
1.00 |
|
1.02 |
1.0E05 |
1.00 |
|
1.03 |
2.0E05 |
1.01 |
|
1.08 |
3.0E05 |
1.03 |
|
1.14 |
4.0E05 |
1.06 |
|
1.21 |
5.0E05 |
1.11 |
|
1.29 |
6.0E05 |
1.18 |
|
1.36 |
7.0E05 |
1.26 |
|
1.44 |
8.0E05 |
1.35 |
|
1.53 |
9.0E05 |
1.44 |
|
1.61 |
1.0E04 |
1.53 |
|
1.69 |
2.0E04 |
2.39 |
|
2.48 |
3.0E04 |
3.21 |
|
3.24 |
4.0E04 |
4.00 |
|
3.97 |
5.0E04 |
4.77 |
|
4.68 |
6.0E04 |
5.52 |
|
5.37 |
7.0E04 |
6.26 |
|
6.05 |
8.0E04 |
6.98 |
|
6.72 |
9.0E04 |
7.70 |
|
7.38 |
1.0E03 |
8.41 |
|
8.03 |
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
9.0
8.0
7.0
6.0
5.0
4.0
3.0
2.0
1.0 1.0E05
1.0E04 |
1.0E03 |
1 − формула (3.22);
2 − формула (3.23);
Рисунок 3.2 − Зависимость ветрового множителя от радиуса капли (при t=10 °С и P=800 гПа )
®Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
Сдругой стороны, масса капли может быть выражена через объем капли и ее плотность:
M = |
4 |
В r3 . |
(3.26) |
|
3 |
||||
|
|
|
Следовательно, после приравнивания (3.25) и (3.26), число молекул в капле и может быть выражено через радиус этой капли:
|
4 |
|
|
В |
3 |
|
|
|
n= |
|
|
|
|
|
r |
. |
(3.27) |
3 |
|
mH |
|
|||||
|
|
|
O |
|
|
|||
|
|
|
|
2 |
|
|
|
|
Изменение числа молекул в капле при изменении радиуса капли, с учетом (3.27) равны:
dn=4 |
В |
r2 dr . |
(3.28) |
|
|
||||
|
mH |
O |
|
|
|
2 |
|
|
|
В этом случае, скорость изменения радиуса капли с учетом (3.24) и (3.28) равна:
dr |
|
DП mH2 O |
|
|
|
|
= |
|
N П |
−N ПВ cR cF . |
(3.29) |
d |
r |
||||
|
|
В |
|
|
|
®Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
Вуравнениях (3.24) и (3.29) вместо концентрации молекул водяного
пара N П и N ПВ удобно пользоваться значениями давления водяного пара:
N П = |
e |
, |
(3.30) |
|
k T |
||||
|
|
|
и
N ПВ |
= |
E ПВ |
, |
(3.31) |
|
||||
|
|
k T В |
|
где T В − температура поверхности капли, К.
С учетом (3.30) и (3.31), предположении, что T ≈T В , уравнения (3.24)
и (3.29) принимают вид:
dn |
= |
4 r DП |
e−EПВ cR cF , |
(3.32) |
d |
k T |
и
dr |
|
DП mH2 O |
|
|
|
= |
|
e−EПВ cR cF . |
(3.33) |
d |
k T r |
|||
|
|
В |
|
|
В результате интегрирования уравнения (3.33) получается формула для расчета радиуса капли r по значению начального радиуса r0 и интервала времени :
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
2 |
2 DП mH2 O |
|
|
||
r= r0 |
|
|
|
(3.34) |
|
В k T e−E |
ПВ cR cF . |
||||
|
Параболическая форма закона диффузионного роста капель (3.34) приводит к сужению распределения капель по размерам в процессе роста.
Например, если имеются две капли с начальными радиусами 1 мкм и 10 мкм и растут при одинаковых условиях, то в то время, когда маленькая капля вырастет до размера 10 мкм, большая капля вырастет только до размера 14 мкм.
Можно сделать вывод, что маленькие капли быстрее растут и испаряются, чем большие, при одних и тех же внешних условиях.
Из уравнения (3.34) можно определить время полного испарения капли
r=0 с начальным радиусом r0 :
|
|
r2 k T |
. |
|
|
исп= |
|
0 В |
|
(3.35) |
|
|
|
|
|||
2 D |
П mH2 O E |
ПВ−e cR cF |
|
||
|
|
|
Если капля падает со скоростью V в потоке воздуха со скоростью U ,
то капля пройдет путь lисп до полного испарения, равный:
lисп= V −U исп . |
(3.36) |
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
Таблица 3.4 – Путь, проходимый падающей каплей до полного испарения
(при T = 280 K и S = 0.8)
r0 , мкм |
lисп , м |
1 |
6 10−6 |
3 |
1.7 10−4 |
10 |
2.1 10−2 |
30 |
1.69 |
100 |
208 |
150 |
1050 |
Из данных таблицы видно, что путь падения капель резко растет с
увеличением радиуса |
выше |
100 мкм. Это значение |
используется для |
разделения капель на облачные и дождевые. |
|
||
Облачные капли − это маленькие капли, которые испаряются вскоре |
|||
после выпадения из облака ( r |
< 100 мкм). |
|
|
Дождевые капли − это большие капли, которые достигают поверхности |
|||
Земли не испаряясь |
( r > 100 мкм). Дождевые капли |
с радиусом около |
|
100 мкм являются каплями мороси. |
|
3.2 Диффузионный рост капель с учетом эффекта нагрева
Скорость роста капли определяется не только скоростью диффузии водяного пара к ее поверхности, но и скоростью рассеивания теплоты, выделяющейся на капле при конденсации.