Файл: Стернин Л.Е. Основы газодинамики двухфазных течений в соплах.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ак[ У 2яmlkTd |
j / 2 n m lkTs |
J |
В этой формуле давление насыщения ps над каплей радиуса
г должно определяться по формуле Томсона |
|
ps= P-{Ts)exp |
• |
Сильная зависимость экспоненциального члена от г умеет место только для очень малых капель. Поэтому в большинстве случаев будем считать ps~ p ^ ( T s) .
1 Первый и третий режимы часто называются газокинетическим (или ки нетическим) и диффузионным режимами.
4* |
99 |
Формулу для прироста массы капли в единицу времени при свободномолекулярном режиме при этом условии можно запи сать в виде
М I
4я/"2ак
У 2яRTs
В области II континуального режима роста перенос молекул пара по направлению к капле обусловлен диффузией. Обозна чим через D коэффициент диффузии и будем считать, что он не зависит от расстояния до капли. Поток массы через сферу ради уса 2 (рис. 3.6) с центром в точке 0 при росте капли определя ется формулой [124]
Мп = 4яz 2D dQ_
d z
Интегрируя это выражение по z от г = г+ Д до бесконечности (индексы, соответствующие бесконечности, опущены) и исполь зуя уравнение состояния p = qRT, получим
|
(3. 39) |
Условие, выражающее баланс массы, имеет вид |
|
4я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 |
Ѵ т -т“ + |
ѵ |
У |
2я Ts |
|
|||
QbR |
|
|
г + Д |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
(3.41) |
Отсюда, |
если принять Td~ T « Г к, |
получим формулу |
Фукса |
||||||||
[124] |
|
|
|
|
|
|
|
|
|
|
|
|
° [p-P'oi'Ts)}1 |
|
— |
|
\ |
f |
2я |
-1 |
(3-42) |
||
|
|
|
RT |
|
|||||||
|
QbRT |
|
Г + A |
ак |
У |
|
|
|
Если радиус капель мал, то первым слагаемым в круглой скобке уравнения (3.42) можно пренебречь, при этом получим
100
формулу Кнудсена для скорости роста капли в свободномолеку лярном режиме
j aK P— /’«о (Ts)
(3. 43)
r ~~QB Ѵ Ш Т
В этом случае скорость роста капли не зависит от ее размера. Если же радиус капли настолько велик, что А<СК то из форму лы (3.42) получим
D |
[р— р~{?а)\ |
D |
- f |
2я |
|
|
« к |
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