Файл: Стернин Л.Е. Основы газодинамики двухфазных течений в соплах.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)

 

АТ

 

 

(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/32z)] 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