ЛАБОРАТОРНАЯ РАБОТА №5 Исследование процесса кристаллизации переохлажденных капель в облаках Цель работы:

ЛАБОРАТОРНАЯ РАБОТА №5
Исследование процесса кристаллизации
переохлажденных капель в облаках
Цель работы:
Исследовать зависимость процесса замерзания переохлажденных
капель от параметров состояния атмосферы в различных типах облаков.
Основные теоретические положения
Зависимость скорости гомогенной (без участия инородных частиц)
кристаллизации капель водного раствора от температуры определяется по
формуле:
max
T   exp   LВЛ  ,
  GВЛ
( гом)
J ВЛ
 J 0( гом)  exp 

 k  T  
 k T  
где
(1)
(гом)
 скорость гомогенной кристаллизации капель, м-3·с-1;
J ВЛ
J 0( гом)  постоянная, равная 10 30 м-3·с-1;
max
T   энергия, затрачиваемая на образование кристалла льда в
GВЛ
воде при температуре T  , Дж;
k  постоянная Больцмана, равная 1.38  10 23 Дж/К;
T   температура воздуха, при которой скорость гомогенной
кристаллизации
гомогенной
водного
раствора
кристаллизации
равна
чистой
скорости
воды
при
температуре T, К;
LВЛ
 теплота фазового перехода вода-лед при температуре T  ,
определяемая
по
формуле
 1004209 .186  7766 .15  T   10.5  T   m
2
H 2O
Дж;
Энергия, затрачиваемая на образование в воде кристаллов льда с
равновесным радиусом при температуре T  равна:
G
где
max
ВЛ
3
16   ВЛ
 mH2 O
T  
,
2
2  T0 
2
3  LВЛ  ln     Л
T
2
(2)
 ВЛ  коэффициент поверхностного натяжения на границе вода-лед
при
T,
температуре
определяемый
по
формуле
 0.0397875  0.00025  T  , Дж/м2;
mH 2O  масса молекулы воды, равная 2.99152  10 26 кг;
T0 – температура плавления чистого льда, равная 273.16 К;
 Л – плотность льда при температуре T  , определяемая по
формуле 956.5756 -0.144886  T  кг/м3.
Температура, при которой скорость кристаллизации водного раствора
равна
скорости
кристаллизации
чистой
воды
при
температуре T,
определяется формулой:




T   273.16  103.6  ln 1  aw  aw*  15.613  ln 2 1  aw  aw* 


 54.118  ln 3 1  aw  aw* ,
где
a w* – равновесное
соответствии с формулой:
значение
активности
(3)
воды,
определяемое
в
3323730




210368

131
.
438

T


41729
.
1

ln
T


T
aw*  exp 
,
8
.
31441

T




где
(4)
T  температура воды, К.
Зависимость скорости гетерогенной (с участием нерастворимых
инородных частиц в воде) кристаллизации капель водного раствора от
температуры определяется по формуле:
( гет )
J ВЛ
 J 0( гет)  exp 0.575  T   T0  ,
(5)
( гет )
где J ВЛ
 скорость гетерогенного замерзания капель, м-3·с-1;
J 0( гет )  постоянная, равная 10 6 м-3·с-1.
Зависимость полной (суммарной) скорости кристаллизации капель J ВЛ
от температуры определяется формулой:
( гом)
( гет)
.
J ВЛ  J ВЛ
 J ВЛ
(6)
Изменение концентрации переохлажденных капель за единицу времени
определяется выражением:
dN В
4
    r 3  J ВЛ  N В ,
d
3
где
N В  концентрация капель, м-3;
r  радиус замерзающей капли, м;
(7)
J ВЛ  скорость кристаллизации капли, м-3·с-1;
  время, c.
Изменение концентрации кристаллов льда (замерзших капель) за
единицу времени определяется выражением:
dN Л 4
   r 3  J ВЛ  N В ,
d
3
где
(8)
 концентрация кристаллов льда, м-3.
NЛ
Численное интегрирование дифференциальных уравнений (7) и (8)
можно осуществить, например, методом прямоугольников:
dN В(i )
  ,
d
(9)
dN Л(i )

  ,
d
(10)
N В(i 1)  N В(i ) 
N
( i 1)
Л
N
(i )
Л
 (i1)   (i )   ,
где
i
–
(11)
номер итерации.
Начальный момент времени соответствует моменту понижения
температуры капли ниже температуры плавления льда.
Скорость
охлаждения
поднимающегося
скоростью w можно определить по формуле:
воздуха
в
облаках
со
dT
  ВА  w ,
d
где
(12)
dT
d
 скорость охлаждения капли, K/с;
 ВА
 влажноадиабатический градиент температуры воздуха, K/м;
w
 скорость восходящих движений воздуха, м/с.
Значение влажноадиабатического градиента можно вычислить по
формуле:
LВЛ  E ПВ
R T
 ВА   А 
,
2
LВЛ  E ПВ
P  0.622 
C P  RП  T 2
P  0.622 
где
А
(13)
 сухоадиабатический градиент температуры воздуха, равный
0.0098 K/м;
P
 атмосферное давление, Па;
LВЛ
 удельная
теплота
фазового
перехода
вода-лед,
равная
2.5  10 6 Дж/кг;
E ПВ
 давление
насыщения
водяного
пара
над
плоской
поверхностью чистой воды (см. лабораторную работу №1), Па;
R
 газовая постоянная сухого воздуха, равная 287 Дж/(кгK);
T
 температура воздуха, K;
CP
 удельная теплоемкость воздуха при постоянном давлении,
равная 1005 Дж/(кгK);
RП
 газовая постоянная водяного пара, равная 461 Дж/(кгK).
Получается, что при подъеме воздуха в облаках, каждому моменту
времени
соответствует
определенная
температура
объема
воздуха,
следовательно, можно рассчитать скорость замерзания капель в зависимости
от температуры:
dN В
dN В
 d ,
dT
dT
d
(14)
и скорость появления кристаллов льда в зависимости от температуры:
dN Л
dN Л
 d .
dT
dT
d
(15)
Численное интегрирование дифференциальных уравнений (14) и (15)
можно также осуществить методом прямоугольников:
N
( i 1)
В
N
( i 1)
Л
N
(i )
В
dN В(i )

 T ,
dT
(16)
N
(i )
Л
dN Л(i )

 T ,
dT
(17)
T (i1)  T (i )  T .
(18)
За начальную температуру принимается температура плавления льда и
далее осуществляется уменьшение T с шагом T .
Скорость замерзания капель зависит от их размеров. По этому, с
течением времени концентрации капель разных размеров убывают с разной
скоростью. Момент времени, когда концентрация кристаллов льда станет
равной концентрации переохлажденных капель, называется временем
кристаллизации  КР капель радиусом r , а соответствующая этому моменту
времени температура, называется температурой интенсивной кристаллизации
капель TИК .
При расчетах удобно использовать относительное число замерзших
капель радиусом r , которое определяется по формуле:
P
где
P
NЛ
,
NВ  NЛ
(19)
 доля замерзших капель, [0, 1].
Поскольку P зависит от отношения концентраций, то P  0 при
отсутствии
кристаллов
льда
( N Л  0)
переохлажденных капель воды ( N В  0).
и
P 1
при
отсутствии
Вариант №1
Исследование вклада гомогенного и гетерогенного механизмов в процесс
кристаллизации переохлажденных капель воды
Цель работы:
Исследовать зависимость от температуры эффективности механизма
гомогенного и гетерогенного замерзания капель переохлажденной воды.
Порядок выполнения работы
Рассчитать
зависимость
скорости
гомогенного
замерзания
переохлажденных капель воды от температуры.
Рассчитать зависимость скорости гетерогенного замерзания капель от
температуры.
Рассчитать
зависимость
суммарной
скорости
кристаллизации
переохлажденных капель от температуры.
Провести расчеты для диапазона температур от 233.15 до 273.15 К с
шагом 1 К при a w  0.95 и 1 .
( гом)
( гет)
(T ) , J ВЛ
(T ) и J ВЛ (T ) .
Построить графики зависимости J ВЛ
Выполнить анализ полученных результатов. Определить в каком
диапазоне температур преобладает гетерогенная кристаллизация, а в каком
гомогенная кристаллизация воды.
Вариант №2
Исследование размера замерзающих переохлажденных капель от
температуры воздуха
Цель работы:
Исследовать зависимость радиуса замерзающих переохлажденных
капель от температуры воздуха при различных скоростях охлаждения
воздуха.
Порядок выполнения работы
Рассчитать зависимость радиуса замерзающей капли от температуры
воздуха.
Для этого необходимо для капель с радиусом r  10 7 ; 10 6 ; 10 5 ;
10 4 ; 10 3 м провести расчеты доли замерзших капель W в диапазоне
температур от 273.15 до 233.15 К с шагом -1 К при значении скорости
охлаждения воздуха
dT
= -0.001 K/с. Из условия равенства концентраций
d
переохлажденных капель и кристаллов льда, то есть P  0.5 , определить
температуру замерзания капель. Начальную концентрацию капель считать
равной N В  108 м-3, а кристаллов льда N Л  0 м-3.
Построить графики зависимости r (T ) .
Выполнить анализ полученных результатов.
Вариант №3
Исследование фазового состояния облаков
Цель работы:
Исследовать соотношение переохлажденных капель воды и замерзших
капель (кристаллов льда) в облаках в зависимости от температуры воздуха и
типа облачности.
Порядок выполнения работы
Рассчитать относительное число замерзших капель P радиусом r .
При расчете скорости охлаждения воздуха считать изменение
температуры с высотой влажноадиабатическим.
Повести расчеты по формуле (19) с учетом формул (16) – (18) для
диапазона температур от 273.15 до 233.15 К с шагом -1 К при параметрах
слоистообразных и конвективных облаков. Начальную концентрацию
капель считать равной N В  108 м-3, а кристаллов льда N Л  0 м-3.
Для слоистообразных облаков принять следующие значения:
- атмосферное давление P  700 гПа;
- вертикальная скорость w  0.01 м/с;
- средний радиус капель r  5 мкм.
Для конвективных облаков принять следующие значения:
- атмосферное давление P  500 гПа;
- вертикальная скорость w  10 м/с;
- средний радиус капель r  25 мкм.
Построить графики зависимости PСЛ (T ) и PКОНВ (T ) .
Выполнить анализ полученных результатов.
Вариант №4
Исследование процесса замерзания капель в облаках слоистых форм и
тумане
Цель работы:
Исследовать зависимость времени полного замерзания капель от
температуры воздуха и размера капель.
Порядок выполнения работы
Рассчитать относительное число замерзших капель P радиусом r .
Провести расчеты по формулам (9) – (11), (19) от момента времени
  0 с до момента, когда концентрации переохлажденных капель и
кристаллов льда станут равны, то есть P  0.5 ,
при
трех
значениях
температуры воздуха T  243.15; 253.15; 263.15 К для капель с радиусом
r  6 мкм (средний радиус капель в слоистых облаках) и r  3 мкм (средний
радиус капель в тумане) и начальной концентрацией капель N В  108 м-3.
Построить графики зависимостей PСЛ ( ) и PТУМ ( ) .
Выполнить анализ полученных результатов.
Определить время, необходимое для полного замерзания капель при
данных температурах воздуха.
Вариант №5
Исследование зависимости времени кристаллизации капель от температуры
Цель работы:
Исследовать процесс замерзания капель чистой воды при различных
температурах.
Порядок выполнения работы
Провести расчет времени кристаллизации капель чистой воды
радиусом r  1 мкм при температурах от 233.15 до 273.15 К с шагом 5 К.
Расчеты произвести отдельно для гомогенного и гетерогенного механизма
кристаллизации капель воды. Начальную концентрацию капель считать
равной N В  108 м-3, а кристаллов льда N Л  0 м-3.
( гом)
( гет)
Построить на одном графике зависимости  КР
(T ) и  КР
(T ) , а на
( гом)
( гет )
другом графики зависимости J ВЛ
(T ) и J ВЛ
(T ) .
Выполнить анализ полученных результатов.
Вариант №6
Исследование влияния растворенных веществ на процесс кристаллизации
переохлажденных капель водных растворов
Цель работы:
Исследовать зависимость от активности воды скорости гомогенной и
гетерогенной кристаллизации капель водных растворов.
Порядок выполнения работы
Рассчитать
зависимость
скорости
гомогенной
и
гетерогенной
кристаллизации переохлажденных капель раствора от активности воды.
Провести расчеты для диапазона активности воды от aw* до 1
с шагом 0.01 при T  243.15; 253.15; 263.15 К .
( гом)
( гет )
Построить графики зависимости J ВЛ
(aw ) и J ВЛ
(aw ) .
Выполнить анализ полученных результатов.
Вариант №7
Исследование процесса замерзания капель в облаках
Цель работы:
Исследовать зависимость времени полного замерзания капель от
температуры воздуха и наличия растворенных веществ.
Порядок выполнения работы
Рассчитать относительное число замерзших капель P радиусом r .
Провести расчеты по формулам (9) – (11), (19) от момента времени
  0 с до момента, когда замерзнут все капли, то есть
трех
значениях
N В  0 м-3, при
температуры воздуха T  243.15; 253.15; 263.15 К и
значениях активности воды a w  0.9 и 1 для капель с радиусом r  5 мкм и
начальной концентрацией капель N В  108 м-3.
Построить графики зависимостей P( ) .
Выполнить анализ полученных результатов. Определить время,
необходимое для полного замерзания капель.