Файл: Стернин Л.Е. Основы газодинамики двухфазных течений в соплах.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 27.06.2024

Просмотров: 165

Скачиваний: 4

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

В результате дифференцирования уравнения (3.35) по g в соответствии с формулой

 

Г

d f

1

0 Д Ф 1

 

/ = — Л тг-+ — £>/ —

 

I

d g

kT

dg J

определяется скорость нуклеации [158]

D

W ch W (1 — GKp)

oo

■2я ^

'■f о ^ Г

sh W

 

 

 

 

л-=1

 

 

 

 

=

d f

 

D

g-e,Кр

кр

d g

пАп exp

^ — D

tj cos (япОкр)

Рис. 3.4.

Распределение комп­

Рис. 3.5. Результаты

расчетов не­

 

лексов по времени (по данным

стационарной

нуклеации:

 

 

Куртни [143])

/ —0,05 -ІО -6 с;

2—0,10 ■10—6 с;

3—0,5Х

 

 

 

Х Ю -’ с;

4 - 1 ,0 0 .1 0 - '

с;

5—4,30 • 10—' с;

 

 

 

 

 

6—00

 

 

 

При tу оо из выражения (3.36)

следует формула для квазистационарно-

го процесса нуклеации. При практическом использовании

формулы (3.36)

fo

подбирается таким образом, чтобы в пределе tу°°

эта формула

совпадала

с формулой Фольмера (3. 16).

 

 

 

 

 

 

 

На рис.

3.5 представлены результаты расчетов по формуле (3.35). Пред­

ставлена функция распределения комплексов

по размерам

в различные

мо­

менты времени. Расчеты были проведены для паров воды при Т—300 К и р«*

«*0,1 МПа. Видно, что равновесное распределение достигается приблизительно за 5 ■10_6 с.

Кантровиц [152] пренебрег членом в правой части уравнения (3.33), обус­ ловленным энергией образования комплексов. В уравнении, выведенном Пробстином (3.34), это равносильно отбрасыванию члена, содержащего Е. В ре­ зультате было получено уравнение

 

а /

аѵ

 

dt

dg2 ’

Решение этого уравнения при тех

же начальных и граничных условиях, что

и (3.34), можно записать в следующем виде [152]:

4

3739

97


(3.37)

(В отличие от Пробстина Кантровиц приближенно считал, что

Используя формулу суммирования Пуассона и оставляя только основной член, Кантровиц это выражение записал в виде

(3.38)

Эта формула автоматически получается из формулы Пробстина (3.36) при W ’’О и g Kp = r. Из формулы (3.38) следует, что для больших скоростей нуклеации требуется время

т

D '

Влияние пренебрежения членом, обусловленным энергией образования комплекса, можно оценить, сравнивая полученное решение с решением Проб­ стина. Кантровиц подсчитал, что для достижения 0,1 величины скорости нук­ леации по Пробстину требуется на 50% времени больше. Однако ряд (3.36) намного сложнее формулы (3.32). Сравнение результатов, полученных Куртни, с результатами расчетов по приближенным формулам (3.36) и (3.38) по­ казывает, что в последнем случае промежуток времени задержки приблизи­ тельно на 10% короче, что является вполне приемлемым для оценок.

Вопросам нестационарной нуклеации посвящены также вышедшие недав­ но работы [9] и [13].

§ 3.2. Рост конденсированных частиц и обмен энергией

Формулы для скорости роста капли. Процесс конденсации состоит из процессов нуклеации и роста конденсированных час­ тиц. Разрушение метастабильного состояния системы, обуслов­ ленное возникновением зародышей, не означает, что происходит мгновенный переход в равновесное состояние. Термодинамическое равновесие достигается только после снятия пересыщения и вы­ равнивания температуры окружающей среды и капли. Это про­ исходит в процессе роста конденсированных частиц: нуклеация не снимает пересыщения, так как из-за малого размера зароды­ шей их суммарный объем весьма мал.

В зависимости от соотношения между размером капли (за­ родыша) и длиной свободного пробега молекул газа существуют три режима ее роста — свободномолекулярный (размер капли намного меньше длины свободного пробега молекул), режим со скольжением (размер капли соизмерим с длиной свободного пробега молекул) и континуальный режим роста (радиус капли

98

намного превосходит длину свободного пробега молекул) *. В об­ щем случае радиус капли может изменяться от нескольких анг­ стрем до микрона и больше, в результате чего процесс роста может протекать во всех трех режимах, что существенно услож­ няет задачу.

 

Будем считать, что капля характеризуется радиусом г, одно­

родной по всей капле температурой Ts и давлением

насыщения

Ps,

соответствующим

этим

радиусу

Т,Р

и температуре. Выведем условие баланса

массы для капли. Следуя Фуксу [124],

 

разделим область вблизи капли на две

 

области (рис. 3.6).

Предположим, что

 

в

непосредственной

близости

к капле

 

концентрация паров постоянна и процесс

 

роста определяется

свободномолекуляр­

 

ным режимом. Эта область (область I)

 

кольцевая с шириной по радиусу, рав­

 

ной А. Вне области /

области II)

 

механизм роста определяется континуаль­

 

ным режимом. Область, характеризуе­

 

мую режимом со скольжением, будем ус­

 

ловно представлять себе

в

виде линии,

 

представляющей собой окружность ради­ Рис. 3.6.

К расчету роста

уса г + А. Параметрам на этой окружности

частиц

придадим индекс d. Количество молекул

пара, соударяющихся в единицу времени с единицей поверхно­ сти капли и остающихся на ней, в соответствии с выражением (3. 15), равно aKpd>y 2nm1kTd. Количество молекул, испаряющих­ ся в единицу времени с единицы поверхности капли, определяет­ ся из условия равновесия при температуре Та и давлении ps:

«к p J V ^ m xkTs.

Прирост массы капли в единицу времени

Pd

Ps

~1

МІ = 4яггт1ак[ У mlkTd

j / 2 n m lkTs

J

В этой формуле давление насыщения ps над каплей радиуса

г должно определяться по формуле Томсона

 

ps= P-{Ts)exp

Сильная зависимость экспоненциального члена от г умеет место только для очень малых капель. Поэтому в большинстве случаев будем считать ps~ p ^ ( T s) .

1 Первый и третий режимы часто называются газокинетическим (или ки­ нетическим) и диффузионным режимами.

4*

99


Формулу для прироста массы капли в единицу времени при свободномолекулярном режиме при этом условии можно запи­ сать в виде

М I

4я/"2ак

У RTs

В области II континуального режима роста перенос молекул пара по направлению к капле обусловлен диффузией. Обозна­ чим через D коэффициент диффузии и будем считать, что он не зависит от расстояния до капли. Поток массы через сферу ради­ уса 2 (рис. 3.6) с центром в точке 0 при росте капли определя­ ется формулой [124]

Мп = z 2D dQ_

d z

Интегрируя это выражение по z от г = г+ Д до бесконечности (индексы, соответствующие бесконечности, опущены) и исполь­ зуя уравнение состояния p = qRT, получим

 

(3. 39)

Условие, выражающее баланс массы, имеет вид

 

r*QBd- ^ = M, = M,h

(3. 40)

dt

 

где используется равенство потоков массы на границе областей

І и II.

Из условия Мі = Мц можно найти неизвестное давление

Р а = рЦ г + P~(TS

У TdTs ] У

f T s

 

, (г +

Д) D У 2nRTs

 

 

 

 

+!

 

RTdr2aK

 

Подставляя это выражение для pd

в формулу

(3.39)

и ис­

пользуя равенство (3.40), получим

 

 

 

 

 

 

 

 

D

т7г5

■P~{TS)

Г2

Ѵ т -т“ +

ѵ

У

Ts

 

QbR

 

 

г + Д

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.41)

Отсюда,

если принять Td~ T « Г к,

получим формулу

Фукса

[124]

 

 

 

 

 

 

 

 

 

 

 

 

° [p-P'oi'Ts)}1

 

 

\

f

-1

(3-42)

 

 

 

RT

 

 

QbRT

 

Г + A

ак

У

 

 

 

Если радиус капель мал, то первым слагаемым в круглой скобке уравнения (3.42) можно пренебречь, при этом получим

100


формулу Кнудсена для скорости роста капли в свободномолеку­ лярном режиме

j aK P— /’«о (Ts)

(3. 43)

r ~~QB Ѵ Ш Т

В этом случае скорость роста капли не зависит от ее размера. Если же радиус капли настолько велик, что А<СК то из форму­ лы (3.42) получим

D

[р— р~{?а)\

D

- f

 

« к

V

(3. 44)

евRT

 

RT

При ѵ = — \/

2п

<£г приходим

к формуле Максвелла для

« к V

RT

 

 

 

 

континуального режима

гQaRTrD [p — P - (Ts)\.

Вотличие от предыдущего случая здесь скорость роста об­ ратно пропорциональна радиусу капли.

Формулы для температуры капли. Рассмотрим теперь баланс энергии между каплей и окружающим пространством. Как и в случае массообмена рассмотрим область / вблизи капли, где определяющим является свободномолекулярный режим (толщи­ на этой области А может не совпадать с толщиной, которая вво­ дилась при изучении массообмена, однако, другого обозначение вводить не будем) и область II, где определяющим механизмом передачи тепла является теплопроводность (континуальный ре­ жим). Излучением газа и капли пренебрежем. Если обозначить через / число степеней свободы молекулы, то при соударении с поверхностью она имеет среднюю энергию (j+\)kTd/2l [61]. Энергия всех падающих молекул при соударении с единицей по­ верхности капли определяется выражением

___ Pd___

] 4~ 1

у'р

V 2птргТа

2

d'

Энергия отраженных молекул равна

(1 ' «к) Pd j 1

Y2nmikTd 2 г’

где Tr —-температура, с которой молекулы отражаются от кап­ ли.

1 Для калорически совершенного газа j = 2 /(к—1), где и — показатель адиабаты.

101


Испаряющиеся молекулы уносят энергию, равную энергии конденсирующихся молекул пара, находящегося в равновесии с каплей,

____ФкPs

j ~t~1

 

У 2nm\kTs

2

5

Если учесть присутствующий инертный газ, то количество энер­ гии, передаваемое инертному газу, равно разности энергий па­ дающих и отраженных молекул. Обозначим параметры инертно­ го газа штрихом, тогда выражение для этой энергии запишется как

р

У + 1kr

_

р'

У + 1 uTd.

У~2п m[kTd

2

т

У 2nm[kTd

2

Суммируя эти составляющие, получим поток энергии в сво­ бодномолекулярном режиме на единицу поверхности капли

ак + (1 — ак) Р Pä’- ^ K T ä - T , ) -

 

У р '

J + 1k(Td— Ts)-f-

V 2 n m 1kTd

 

y'2nm'1kTd

2

 

Pd

 

Ps

I + l kT<

(3. 45)

V 2nmikTd

Y2nmikTs

 

 

Энергия, уносимая отраженными молекулами, характеризует­ ся коэффициентами аккомодации, которые определяют долю энергии, уносимую от капли или отдаваемую капле. Эти коэффи­ циенты определяются следующим образом

T r - T d

g,_ K - T ä

P s - T d

P

Ts - T d

Если Tr=Td, то ß= 0 и энергия молекулы после отражения от поверхности капли не изменяется. Если Tr= Ts, то ß = l и моле­ кула отдает или уносит (в зависимости от соотношения Td^ ,T s) максимальное количество энергии. Таким образом, O ^ ß ^ l ; аналогично изменяется и ß'.'

Поток тепла (см. рис. 3.6) при континуальном режиме тече­

ния определяется уравнением Фурье Q= •—Xy^nz2dT/dz, которое после интегрирования в области II при z— >-оо можно записать в виде

Q= 4^(r + A)XT(7rf- n ,

(3.46)

где Q — количество отводимого тепла, а — коэффициент теп­ лопроводности пара.

Подводимая к капле энергия затрачивается на увеличение ее удельной внутренней энергии е (Ts) и поверхностной энергии.

102