Файл: Стернин Л.Е. Основы газодинамики двухфазных течений в соплах.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 27.06.2024
Просмотров: 163
Скачиваний: 4
(3.7) м ож но записать в виде
2 а |
(3.9) |
Ѵя,„ — ^4+ ~ «1=0. |
' кр
Используя соотношения
4A — kT \пр + у(Т ),
'Pßoo = VAo o =k T ln Ptx (T) + <f(T),
где k — постоянная Больцмана; cp (Г) — некоторая функция.
Учитывая, что k = mlP, из уравнения (3.9) получим формулу Томсона
, |
Р |
~s = |
2wi |
(3.10) |
ln ---- 7 |
------- |
|||
|
^оо(Г) |
/-крйГ |
|
Капли радиуса ткр, соответствующего формуле Томсона, называются каплями критического размера или зародышами. При заданных внешних условиях кап ли докритического размера находятся в среде ненасыщенного пара и испаря
ются, а капли сверхкритического размера |
(/•> /‘кр) — в среде |
пересыщенного |
|
пара и, следовательно, растут. |
(3.10), пар, |
будучи насыщенным |
для крупных |
Как следует из формулы |
|||
капель и пересыщенным для |
плоской поверхности, является |
ненасыщенным |
для мелких капель. Этот эффект сильно возрастает с уменьшением размеров капель. Так. например, для сравнительно крупных водяных капель — разме ром в микрон—давление насыщения лишь на 0 ,1 % больше давления над плос кой поверхностью, а при радиусе капли 6 -ІО- 4 микрона давление насыщения
для такой капли в шесть раз превосходит давление насыщения над |
плоской |
||||
поверхностью. |
|
|
|
|
|
Отношение |
s= p /p 0O(’T) называется |
пересыщением. |
Аналогичным |
образом |
|
иногда вводят |
переохлаждение Д Г = Г с о —Т, |
где Too— температура |
насыще |
||
ния над плоской поверхностью перехода при заданном давлении р. |
|
||||
Из уравнения Клапейрона—Клаузиуса, записанного в конечном виде, сле |
|||||
дует |
|
|
|
|
|
|
Т |
|
Pool |
|
|
|
= Т |
|
Pool ■ |
|
|
|
|
|
|
|
|
Используя теперь формулу Томсона |
(3.10), |
получим |
(при L = const) |
||
|
АТ |
2а |
|
|
(3.11) |
|
|
|
|
|
ТLrкр0в
Равновесное распределение Гиббса. Рассмотрим равновесное состояние системы и обозначим через g Kp число молекул в капле критического размера. Пусть z —g / g KV, тогда, учитывая уравнения (3.6), (3.8) и (3.9), выражение для работы образования капли, содержащей g молекул, можно записать в виде
ДФ |
В (3z2?3— 2г ) , |
(3.12) |
||
kT |
||||
|
|
|
||
где безразмерный множитель |
|
|
|
|
|
4яг»| |
1/3 |
|
|
В = - |
„2/3 |
(3.13) |
||
|
||||
|
•5кр • |
|||
кТ |
|
|
|
Из формулы (3.12) видно, что ДФ зависит от внешних условий (через коэффициент В) и относительного числа молекул в капле г. На рис. 3. 1 пред ставлены зависимости ДФ(г) при различных внешних условиях. Если пар
находится в термодинамически устойчивом |
состоянии, |
т. е. если |
Фа < Ф в |
(или, что то же s < l ) , зависимость ДФ(г) |
представляет |
собой |
монотонно |
возрастающую функцию. Любое скопление молекул не сможет создать устой чивую каплю, так как критического размера капли в данном случае не суще ствует. Состояние насыщения (фл = фв, s = l) также будет устойчивым, так как радиус капли критического размера равен бесконечности. Если же пар на
ходится |
в термодинамически |
неустойчивом |
состоянии |
(ф л > ф в , |
s > l ) ' то |
|||||||||||||
скопления молекул после достижения ими критического размера |
( z = 1 ), соот |
|||||||||||||||||
ветствующего максимуму ДФ, оказываются спо |
|
|
|
|
|
|||||||||||||
собными к дальнейшему росту. |
|
|
|
боль |
|
|
|
|
|
|||||||||
Таким образом, |
когда |
(рл |
становится |
|
|
|
|
|
||||||||||
шим |
фв, |
может |
начаться |
превращение |
(А-+В), |
|
|
|
|
|
||||||||
скорость |
которого |
ограничена |
необходимостью |
|
|
|
|
|
||||||||||
пройти |
через «потенциальный барьер» ДФкр= |
|
|
|
|
|
||||||||||||
= ДФ (1). |
Из последнего следует, что пар в этом |
|
|
|
|
|
||||||||||||
случае |
нельзя |
|
рассматривать |
неустойчивым |
|
|
|
|
|
|||||||||
в обычной термодинамической трактовке: он мо |
|
|
|
|
|
|||||||||||||
жет сохраняться в течение более или менее про |
|
|
|
|
|
|||||||||||||
должительного времени. Подобное его состояние |
|
|
|
|
|
|||||||||||||
называется метастабильным. Математически оно |
|
|
|
|
|
|||||||||||||
характеризуется |
максимумом |
термодинамического |
|
|
|
|
|
|||||||||||
потенциала ДФ. |
|
|
|
|
|
|
в |
однород |
|
|
|
|
|
|||||
Согласно статистической физике |
|
|
|
|
|
|||||||||||||
ной молекулярной системе происходят флуктуа |
|
|
|
|
|
|||||||||||||
ции — отклонения от нормального состояния, |
при |
Рис. |
3.1. Кривые |
зависи |
||||||||||||||
водящее ее в состояния с большим термодинами |
||||||||||||||||||
ческим потенциалом, т. е. менее вероятные. Такие |
мости ДФ от z для нена |
|||||||||||||||||
флуктуации |
приводят |
к |
появлению |
комплексов |
сыщенного |
и пересыщен |
||||||||||||
(скоплений молекул), |
которые |
при |
описанных |
|
ного |
пара |
|
|||||||||||
выше |
условиях |
могут |
превратиться |
в |
зародыши |
|
|
|
|
|
||||||||
новой |
фазы. |
Согласно |
|
распределению |
|
Гиббса |
|
|
|
|
|
|||||||
количество комплексов Ng, содержащих g |
молекул, |
при равновесном |
состоя |
|||||||||||||||
нии системы равно |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
N g ^ N ^ x p |
/ |
|
ДФ \ |
|
|
|
|
(3.14) |
||||
|
|
|
|
|
|
( — — |
J, |
|
|
|
|
|||||||
где коэффициент Ni, представляющий собой количество |
молекул |
пара, при |
||||||||||||||||
мерно равен N а , т. е. количеству молекул во всей системе. |
|
|
N g от |
|||||||||||||||
На |
рис. |
3.2 |
иллюстрируется распределение Гиббса |
(зависимость |
||||||||||||||
г = ё І ё кр)- |
Видно, |
что после достижения |
минимума |
в точке z = l |
количество |
комплексов, содержащих большее количество молекул, увеличивается. Одна ко при z > l распределение Гиббса является неверным.
Стационарная нуклеация1. Формула Фольмера. Первый важный вклад в теорию образования ядер конденсации сделал Фольмер [164], который считал, что критический размер комплексов определяет «ворота в конденсацию». Он предположил, что все капли, достигшие критического размера, удаляются из рассматриваемой системы и заменяются эквивалентным числом молекул пара. Скорость образования зародышей /„ вычислялась Фольмером как произведе ние числа капель критического размера NKр в единице объема на количество молекул Ь, конденсирующихся при столкновении с единицей поверхности в единицу времени, и на поверхность зародыша
Iv = NKpb-4nrlp.
1 Термин «нуклеация» (образование ядер конденсации или зародышей)
употребляется в ряде работ [11, 15] и других и происходит от английского слова nudeation.
89
Величина |
вычисляется по распределению Гиббса, |
а b—по кинетической |
теории (см., например, работу [150]). |
|
|
|
«кР |
(3.15) |
|
У 2nniikT |
|
|
|
|
где ак — коэффициент конденсации (отношение числа |
конденсирующихся мо |
лекул пара к числу падающих на поверхность комплекса). Пользуясь выра
жениями (3.8), (3.12) и (3.14), получим формулу для определения |
скорости |
|
образования зародышей по Фольмеру |
|
|
3 ьТ |
В ехр (—В). |
(3.16) |
= |
||
При выводе формулы (3.16) предполагалось также, что вероятность столк |
||
новения с каплей одновременно двух |
или большего числа молекул |
мала по |
сравнению с вероятностью столкновения с каплей одной молекулы. Это пред положение исключает возможность перехода капли докритического размера в каплю сверхкритического размера, минуя промежуточный этап. Формула Фольмера (3.16) обладает тем недостатком, что дает завышенное значение скорости образования зародышей, так как при ее выводе не учитывается ис парение капли и дальнейшая конденсация пара на зародышах.
Формула Беккера—Дёринга. Более точное выражение для скорости обра зования зародышей получили Беккер и Дёринг [139]. Они предположили, что образование зародышей происходит в результате реакций последовательного
присоединения к комплексам |
одной молекулы (вероятность присоединения |
|
двух |
и более молекул не учитывалась). Пусть A g представляет собой комп |
|
лекс, |
состоящий из g молекул |
(комплекс класса g). Можно считать, что кон |
денсация происходит в результате протекания следующей цепочки реакций
А\ |
Ai Ä2 |
А<і + Аі Д' A3 |
|
|
(3.17) |
Ag—1 + |
Аі Д: Ад, |
где индекс G > g KV. В процессе роста комплексов число молекул в газе долж но уменьшаться, однако, как и при выводе формулы Фольмера, предполагает ся, что комплексы, содержащие G молекул, удаляются из системы и заменя ются эквивалентным числом молекул. Если обозначить через fg 1 количество
комплексов класса g в единице объема, которые образуются в результате про текания реакций (3.17), то, очевидно, должно быть
/і = Na « Nь
(3.18)
/о = 0.
Последнее равенство говорит о том, что комплексы класса G удаляются из системы.
Количество комплексов, содержащих g молекул, образующихся в единице объема в единицу времени Ig в процессе реакции
Ag—1 + Аі П Ag, |
|
|
определяется выражением |
|
|
lg = bg—\Sg—\ f g - \ |
— a gS g f g , |
(3.19) |
где bg_ 1 — количество молекул, оседающих |
на единице поверхности |
комплек- |
:а класса g —4, в единицу времени; |
|
|
S g_ 1 — площадь поверхности этого комплекса;
1 Если реакции (3.17) протекают равновесно, то fg= N g
90
ае — количество молекул, испаряющихся с единицы |
площади комплекса |
|||||
класса g в единицу времени. (Правая |
часть соотношения |
(3.19) |
представляет |
|||
собой разность количеств комплексов |
класса g, |
образующихся |
в результате |
|||
конденсации одиночных молекул на комплексах, |
состоящих из g — 1 молекулы, |
|||||
и комплексов, убывающих из класса g |
вследствие испарения |
у каждого |
из |
|||
них по одной молекуле). |
|
|
и распадающихся |
в |
||
Из баланса |
количества комплексов, образующихся |
|||||
классах g —1, g |
и g + 1 в результате протекания двух последовательных реак |
|||||
ций, из соотношения (3.19) можно получить уравнение |
|
|
|
|||
|
02 = I g — Jgi i < |
|
|
(3.20) |
которое иногда называют уравнением Беккера—Дёринга. Будем считать про цесс нуклеации квазистационарным, т. е. d{gl d t = 0. В этом случае из уравне ния (3.20) следует, что
Ig — Ig . j = / 2= const.
Таким образом, для квазистационарного процесса, описываемого цепочкой реакций (3.17), имеют место соотношения
I = äjS i / i — a2S2f 2 = • • • = і / г —і — agS g f g =
= • • •ba - i S G - i f o - i — a G S G ^ G -
Отсюда получаем последовательно
|
|
Г |
ь, |
/ |
|
|
|
S2/2 : д 2 |
— |
= — 5 ^ ! - — |
|
|
|
|
s ifi ■ «2 |
а2 |
а2 |
|
|
|
|
SGfG |
0 - 1 |
|
|
|
|
|
J0 - l A ? - l " |
|
|
|
||
|
|
Ь0 -1-*2 |
b\b2. ..b 0-1 |
Sfli =0, |
||
а у1 |
+ Г _ І + |
+ ' aa _ r |
. . a 2 |
+ •a2a3. . . a |
о |
|
а О \ |
a G - 1 |
|
|
|
|
причем обращение в нуль этого выражения обусловлено равенством /с = 0 .
Поделим обе части этого равенства на коэффициент при SiNi и получим окончательно
bjSiN 1 |
_£г_ . |
д2 дз |
д2 дз• • • а д —х |
/ |
Ь2 + |
М з + - " |
(3.21) |
+ Ь2Ьз ...* 0_ і |
Беккер и Дёринг сделали следующие допущения: 1) количество испаряю щихся молекул с поверхности комплекса пропорционально давлению насы щения для этого комплекса; 2) для капли, находящейося в равновесии, Ькр= = «кр; 3) bg= b не зависит от g [формула (3.15)]. Таким образом,
ag |
ag |
ag |
p g |
|
1 |
_ 1 \1 |
bg |
^кр |
д кр |
Ркр |
. kT |
\ r g |
г кр / J |
|
|
= |
exp [ 2 B g - p1 ( z ~ V 3 — |
1 ) ] . |
|
|
Предполагая, |
что эта формула выполняется и для малых комплексов (g = 2), |
т. е. что малые комплексы обладают поверхностным натяжением (это являет ся недостатком жидкокапельной теории), получим выражение
91
а 2аА. . . а ё |
2 В |
|
|
|
|
------------- = |
ex p {------- |
|
|
|
|
hh-'-bg |
gKp |
л —2 |
|
|
|
|
|
|
|
||
|
g |
/ |
\_L |
|
|
|
2В |
|
^кр \ 3 |
|
|
|
exp |
|
dn ■ |
|
|
|
£кК'Р |
|
|
|
|
: exp [ß (322/3 — 2z — 3gKp2/3 |
|
|
|||
поскольку А/г=1. |
|
|
|
|
|
Подставим это соотношение при g = 2, 3,..., |
(G— 1) в уравнение (3.21) и |
||||
снова заменим суммирование интегрированием, |
получим |
||||
|
|
|
G - l |
|
|
bSjNi. ^ехр |
В [3g,кр |
■2g кр / J g Kp |
®кр |
exp [ß (З,г2/3 — 2z)] dz. |
|
^ |
|||||
I |
|
|
|
|
|
ygKP
Интеграл в правой части характеризуется тем, что подынтегральное выраже ние имеет резкий максимум при 2 = 1 . Это объясняется большим значением
коэффициента В, который для паров воды равен приблизительно 50. Данный интеграл можно вычислить, например, методом Лапласа [46], основанным на разложении экспоненты вблизи максимума в ряд Тейлора —
Зг2/3 — 2г « 1— — (г — 1)2,
О
G—1 |
|
|
|
|
|
|
^кр |
|
е |
|
|
d z |
(3.22) |
^ exp [ß (З2 2/3—2z)] d z z z e B ^ |
|
|
||||
yg,кр |
|
|
|
|
|
|
Таким образом, можно записать |
|
|
|
|
|
|
|
|
|
|
|
|
(3.23) |
После пренебрежения членами 3 g “ 2/ 3 |
и |
2g“ 1 |
по сравнению с единицей по |
|||
лучим формулу Беккера—Дёринга |
|
|
|
|
|
|
Ib- d = |
|
1f |
~~В |
|
(3.24) |
|
Skp |
- |
Зя |
|
|||
|
V |
|
|
|
||
отличающуюся от формулы Фольмера (3.16) предэкспоненциальным |
множи |
|||||
телем, или |
|
|
|
|
|
|
|
S i<j |
|
1 |
|
|
|
/B~DZ=skTgKV / Ш |
ІѴ' |
|
При обычных условиях эксперимента в камере Вильсона эта поправка (т. е. отношение/в_в//щ) имеет величину порядка ІО-2.
Формула Френкеля—Зельдовича. Зельдович [56] и Френкель [123] уточни ли формулу Беккера—Дёринга и предложили более простой ее вывод. Для
92