Лабораторные работы по физике аэрозолей (Чукин В.В.) / clouds
.pdf® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
Высвобождение скрытой теплоты конденсации LПВ приводит к повышению температуры капли относительно окружающего воздуха: начинается процесс диффузии тепла в сторону от капли, который определяется уравнением, аналогичным уравнению для диффузии водяного
пара:
где А
T В cT
dQ |
=4 r А T |
В−T cT , |
(3.37) |
d |
|
|
|
−теплопроводность воздуха, Дж/(м∙с∙К);
−температура поверхности капли, К;
−поправочный коэффициент на размер капли.
Коэффициент cT к скорости диффузии тепла от капли аналогичен коэффициенту cR в уравнении диффузионного роста капель и равен:
|
|
|
1 |
СВ |
|
||||
cT |
= |
r |
|
|
|
, |
(3.38) |
||
|
|
|
|
|
|||||
|
4 A r |
|
|
||||||
|
|
1 |
СВ |
|
|
||||
|
|
r |
2 V |
|
|||||
|
|
|
|
||||||
|
|
|
T |
|
|
|
|
|
где T |
− коэффициент аккомодации тепла, равный |
T r−T |
; |
||
T |
В−T |
||||
T r |
− температура на расстоянии СВ от |
|
|||
|
поверхности капли |
||||
|
(температура молекул, отскочивших от поверхности капли), К. |
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
Поскольку T r≈T В , то коэффициент T мало отличается от 1.
Уравнение сохранения энергии для капли имеет вид:
c |
|
n |
|
dT |
В |
=L |
|
dn |
− |
dQ |
, |
(3.39) |
P H |
|
|
|
|
|
|
||||||
|
2 O |
d |
ПВ |
d |
|
d |
|
|
||||
|
|
|
|
|
|
|
|
|
где cP H |
O |
− |
теплоемкость воды, Дж/К; |
2 |
|
|
|
|
n |
− |
число молекул образующих каплю. |
Для установившегося процесса роста капли, то есть при постоянстве
температуры поверхности капли dTd В =0 выполняется условие сохранения
энергии:
dQ =L |
|
dn |
, |
(3.39) |
|
||||
d |
ПВ d |
|
то есть отток теплоты за единицу времени от капли равен притоку теплоты в результате изменения числа молекул в капле. Символом LПВ обозначена теплота, выделяющаяся при конденсации одной молекулы водяного пара.
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
Температура поверхности капли определяется из уравнения сохранения энергии (3.39) с учетом уравнения скорости изменения количества теплоты (3.37) и скорости увеличения числа молекул в капле (3.32):
T |
В=T |
L |
ПВ DП cR cF |
|
|
e |
− |
E ПВ |
|
. |
(3.40) |
|
|
A k cT |
|
|
T |
В |
|||||||
|
|
|
|
T |
|
|
|
В этом уравнении внешние условия задаются параметрами e и T , а
отношение LПВ DП cR cF слабо зависит от температуры и давления.
Неизвестными параметрами являются температура поверхности
капли T В и |
давление насыщения водяного пара на |
расстоянии СВ от |
|
поверхности |
капли при температуре поверхности |
капли |
EПВ T В ,r , |
определяемая формулой (1.31). |
|
|
Решая совместно уравнения (3.40) и (1.31) можно вычислить температуру поверхности капли и парциальное давление над каплей, а затем по формуле (3.32) рассчитать скорость конденсационного роста капли с учетом эффекта нагрева.
Температура поверхности капли определяется из совместного решения уравнений (3.40) и (1.31) итерационным методом. Задаются параметры окружающего воздуха T и e, размер капли r . Затем задается начальное значение T В (например, равное температуре окружающего воздуха T ) и
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
по (1.31) рассчитывается значение EПВ . Полученное значение подставляется в уравнение (3.40) и определяется новое значение T В . Это значение опять подставляется в (1.31) и процесс повторяется. В результате, каждый раз получаются новые значения T В , которые стремятся к значению,
одновременно удовлетворяющему уравнениям (3.40) и (1.31). Итерационный процесс (процесс последовательных приближений) прекращается, когда разница между рассчитанным значением T В и предыдущим значением T В
окажется меньше некоторого значения :
T Вi −T Вi−1 ≤ . |
(3.41) |
Рассчитанное значение T В соответствует температуре поверхности капли с точностью ± . По определенному значению T В и формуле (1.31)
рассчитывается значение EПВ . Скорость роста капли определяется по формуле (3.32).
Возможно также приближенное аналитическое решение уравнений (3.16), (3.23), (3.24). Результат такого решения определяется уравнением:
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
|
dr |
|
[e−EПВ T ,r ] |
|
|
|
|||
|
|
конд= |
|
|
|
|
|
. |
|
d |
|
L2ПВ В EПВ T |
|
В RП T |
|
(3.26) |
|||
|
|
|
|
|
|
|
r |
|
|
|
|
|
RП T 2 |
DП |
|
Первое слагаемое в знаменателе описывает влияние теплопроводности, а второе диффузии водяного пара. Скорость роста капель, описываемая уравнением (3.26), находится в хорошем согласии со скоростью роста, получаемой их совместного решения уравнений (3.16), (3.23), (3.24).
Концентрация облачных капель
Так, при скорости подъема U =0.15 м/с и ядрах конденсации с массой до m=10−12 кг получено следующее распределение относительной влажности с высотой над основанием облака.
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
Рисунок − Вертикальный профиль относительной влажности воздуха в слоистообразном облаке
Влажность достигает максимума, равного примерно 1.005 (то есть пересыщение равно 0.5 %) при подъеме на 40 м над основанием облака. На этой высоте наблюдается наибольшая скорость роста капель.
Количество образующихся облачных капель зависит как от скорости подъема, так и от концентрации ядер.
При подъеме объема воздуха в облаке, когда влажность становится примерно равна единице S≈1 крупные нерастворимые ядра конденсации увлажняются и превращаются в капли, растворимые ядра превращаются в капли растворов.
При влажности |
S 1 наиболее |
крупные из обводненных ядер |
конденсации r rкр |
продолжают расти. |
Таким образом, число растущих |
капель определяется числом ядер конденсации с радиусами более Поскольку критический размер ядра зависит от влажности воздуха S :
rкр= |
2 ПВ mH2 O |
, |
(3.29) |
|
|||
|
k T В ln S |
|
то при повышении относительной влажности активизируются ядра все меньшего размера и, соответственно, образуется большее число капель. Число
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
образующихся капель в 1 м3 облака N В оказывается равно числу частиц аэрозоля в 1 м3, с радиусом больше rкр :
∞ |
|
N В=∫nА r dr . |
(3.30) |
rкр |
|
Можно сделать вывод, что влажность достигает своего максимального значения в пределах 100 м над основанием облака, а выше этого уровня влажность убывает и приближается к примерно постоянному значению. Так как от влажности зависит количество активизирующихся ядер конденсации, то концентрация облачных капель определяется самым нижним слоем облака.
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
4 Коагуляционный рост капель
В одном кубическом метре облака находится порядка 108 капель, которые часто сталкиваются между собой.
При столкновении капли могут сливаться в одну каплю или отскакивать друг от друга. Образовавшиеся очень крупные капли могут дробиться на более мелкие. В воздухе, насыщенном водяным паром, столкновение капель с радиусом r 100 мкм сопровождается их слиянием. При столкновении более крупных капель с r 100 мкм возможен отскок или дробление.
Процесс укрупнения капель за счет столкновений и последующего слияния друг с другом называется коагуляционным ростом капель.
Существует несколько механизмов коагуляции капель:
гравитационная коагуляция;
броуновская коагуляция;
турбулентная коагуляция;
электростатическая коагуляция;
электродинамическая коагуляция.
4.1Гравитационная коагуляция
Основным механизмом коагуляции является гравитационная коагуляция капель, возникающая в результате разности скоростей падения капель разных размеров под действием силы тяжести.
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
Скорость падения крупной сферической капли радиусом r приближенно описывается формулой:
В |
[ |
] |
|
|
V |
=9.95 |
1−exp −1200 r |
, |
(4.1) |
где V В |
− |
установившаяся скорость падения капли, м/с; |
r |
− |
радиус капли, м. |
Из формулы видно, что большие капли падают быстрее, чем маленькие. Получим уравнение коагуляционного роста капли, которое описывает скорость изменения размера капли за счет слияния с более мелкими каплями. Предположим, что капля с радиусом r падает с установившейся
скоростью через облако более мелких капель с радиусом r1 .
За одну секунду капля с радиусом r пройдет путь V В r 1 с и столкнется со всеми каплями радиусом r1 , находящимися внутри объема:
r r1 2 [V В r −V В r1 ]. |
(4.2) |
В одном метре кубическом облака содержится nВ r1 dr1 капель |
с |
радиусом от r1 до r1 dr1 . Следовательно, число столкновений капли |
с |
радиусом r с каплями радиусом r1 за единицу времени (одну секунду) равно
® Чукин В. В. Физика облаков. Конспект лекций. 2005. (1е неофициальное издание)
произведению объема, проходимого каплей за секунду на концентрацию капель маленьких капель:
dn |
= r r |
1 2 [V |
В r −V В r1 ] nВ r1 dr1 . |
(4.3) |
|
||||
d |
|
|
|
Однако столкновение произойдет не со всеми каплями в этом объеме. Поскольку, по мере приближения маленькой капли к поверхности большой, траектория маленькой капли начинает искривляться и огибать большую каплю по сторонам.
Отношение числа реально захваченных капель по отношению к числу захваченных капель в ситуации, когда траектории маленьких капель не искривляются под влиянием большой капли, называется коэффициентом захвата и обозначается Э r ,r1 . Значения Э r ,r1 , как правило, находятся в пределах от 0 до 1. С учетом коэффициента захвата формула (4.3) примет вид:
dn |
= r r |
1 2 Э r ,r |
1 [V |
В r −V В r1 ] nВ r1 dr1 . |
(4.4) |
|
|||||
d |
|
|
|
|
Скорость роста массы капли равна:
dM |
∞ |
4 |
2 В r13 r r1 2 Э r ,r1 [V В r −V В r1 ] nВ r1 dr1 . |
|
|
=∫ |
(4.5) |
||||
d |
3 |
||||
r |
|
|