Файл: Стернин Л.Е. Основы газодинамики двухфазных течений в соплах.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 27.06.2024
Просмотров: 167
Скачиваний: 4
rKp до оо. Если вспомнить обозначения (3.53) и воспользоваться соотношениями
dQn |
|
|
|
|
|
|
|
|
dt |
,f r n ^ r d r - ' V V « + |
|
, |
|
|
|||
5 |
ГП17 ^ dr==~ |
^ p/ . kp^kp - |
ti ^ 'rf rn~ld r -\-\im |
[ Г - / Г ] , |
(3.57) |
|||
|
||||||||
г кр |
|
|
|
|
гкр |
|
|
|
оо |
|
|
|
|
|
|
|
|
г кр |
|
|
|
|
|
|
|
|
то вместо выражения (3.56) |
получим уравнение |
|
|
|||||
|
dQn |
tl |
J 'rrn~lf d r -f |
rn . |
|
(3. 58) |
||
|
|
dt |
|
rKP |
|
6a KP |
|
|
|
|
|
|
|
|
|
|
|
Обозначим через |
as— Qs/Qz |
массовую концентрацию конден |
||||||
сирующегося вещества, тогда из уравнений (3.54) |
и (3.58) при |
|||||||
п = 3 получим соотношение |
|
|
|
|
|
|||
|
^ = |
Т |
яг«Р' Г /+4я (?вІ |
rr* fd r ' |
|
(з-59) |
из которого ясно, что прирост массы конденсирующегося веще ства определяется как процессом образования зародышей, так
ипроцессом их роста.
Впоследующих параграфах кинетическое уравнение (3.56) будет использовано при исследовании стационарных одномерных, плоских и осесимметричных течений. Для общности результатов удобно при изучении двумерных течений вместо одной независи мой координаты ввести функцию тока ф. В переменных х, ф, г уравнение (3.56) запишется в виде [120]
« 4 ^ + д - |
( ' / ) = — 8 (г ~ гкр) |
(3. 60) |
|
дх |
дг |
Qs |
|
или, раскрывая второй член левой части,
и |
д / |
, |
• df |
/ |
Ъ(Г~ Гкр) - / |
дг |
дх |
{- г —— |
п |
д |
|||
|
1 |
дг |
е2 |
дг |
Это уравнение является линейным уравнением в частных про изводных первого порядка. Уравнения характеристик для него имеют вид
107
d x __ |
dr |
______________ df_____________ |
gi^ |
“ |
r |
I Q s ' 1 (r — г*?) — f d r \ d r |
|
Решение уравнения (3.60) эквивалентно решению системы уравнений (3.61). Если известны первые интегралы системы (3.61), то решение основного кинетического уравнения опреде ляется как произвольная функция от этих интегралов.
§ 3.4. Случаи интегрирования основного кинетического уравнения
Континуальный режим роста капель. Систему уравнений
(3.61) можно проинтегрировать в некоторых практически важ ных случаях. Ниже будут рассмотрены два из них. Метод на хождения двух первых интегралов системы (3.61) был предло жен Бахановым и Буйковым в работе [11] для случая, когда ско
рость роста капли можно представить |
в виде произведения |
|||
г= ф1(г)ф2(х). Такое представление |
возможно, например, при |
|||
использовании формулы Максвелла |
(3.44) для роста крупных |
|||
капель. Представим формулу Максвелла в виде |
||||
|
• __ |
uZ' (х) |
|
(3. 62) |
|
Г ~ |
2 (V + г ) |
’ |
|
|
|
|||
где ѵ = |
D_ |
|
|
|
|
ак |
|
|
|
|
2D |
Р — Роо (Ts) |
dx, |
|
|
б в |
u R T |
|
здесь Хі •— значение координаты х, соответствующее началу нуклеации.
Интегрирование уравнения (3. 62) возможно при двух услови ях: ѵ = const и TS = T. Первое условие выполняется при не очень
сильном изменении температуры, так как v ~ D /]/7 \Второе усло вие можно несколько ослабить, если температура капли незначи тельно отличается от температуры газа. Как показали Баханов и Буйков [12], при T t t T s формулу (3.62) можно заменить экви валентной с эффективными значениями ѵЭф и Д,ф таким образом, что правая часть перестает зависеть от Ts.
Рассмотрим первое обыкновенное дифференциальное уравне ние системы (3.61)
dx dr
После подстановки в это уравнение выражения для скорости роста капли (3.62) придем к дифференциальному уравнению с
108
разделяющимися переменными. Общий интеграл этого |
уравне |
ния имеет вид |
|
(v Jr r f - Z { x ) = C1, |
(3.63) |
где С1 — произвольная постоянная.
Как видно из этого соотношения, оно определяет рост частиц вдоль координаты х, причем радиус частиц при некотором зна чении координаты £ может быть выражен через его значение в другой точке X
г (£)= У Z (I) + [ѵ+ |
г (X)]2 —Z(x) —V . |
(3.64) |
||||||
Второе обыкновенное дифференциальное уравнение |
системы |
|||||||
(3.61) при помощи формулы (3.62) |
можно записать в виде |
|||||||
|
d f |
__ |
f Z ' |
. |
lb {г |
r Kp) |
|
|
|
dx |
2 (г + v)2 |
|
|
Qu |
|
|
|
Это линейное дифференциальное уравнение первого порядка. |
||||||||
С учетом уравнения (3.63) |
его общий интеграл запишется следу |
|||||||
ющим образом |
|
|
|
|
|
|
|
|
— |
- f |
/(£) |
8 (г |
гКр) |
<І\= Сг, |
|
||
Qs (£) и (£) |
|
г + V |
|
|||||
r + v |
J |
|
|
|
причем стоящие под знаком интеграла величины г и гкр соответ ствуют значению х = \.
Общее решение кинетического уравнения (3.60) можно запи сать в виде произвольной функции типа КДСі, С2)= 0, или С2 =
= F(Ci), т. |
е. |
|
|
|
|
|
ь (Г— Гкр) |
r + |
V |
QjMJ |
(3. 65) |
Г + V |
|||
— |
- ^ |
1 - |
|
Для определения произвольной функции необходима задать граничное условие. Пусть в сечении х = х\ функция распределе ния задана
/0 * і. r) = f 0(г),
тогда из выражения (3.65) при х = хи получим
Г+ V |
I г )2]. |
|
( 3 . 6 6 ) |
|
|
|
|
Отсюда следует, что если заменить г выражением |
|
||
у “(ѵ_|_ r f — Z {х) — V , то соотношение (3.66) |
примет вид |
|
|
F[(v + r f - Z ( x j \ |
/ о [ /( ѵ + г)2—ZQ*) —ѵ] |
(3. 67) |
|
|
|
K(V + r)2 — Z (x)
109
Заменяя в уравнении (3.65) правую часть с помощью равен ства (3.67), а значения г и гкр, стоящие под интегралом, по фор муле (3.64), получим
/ ( e ) 8 [ W |
+ |
v ) 2 - Z (;c ) + |
Z($)- |
1кр(£)] d\ -f- |
f i x , г) = ( г + ѵ ) |
|
|
|
|
е (6) и (о у z (6) + ( г |
+ ѵ)2 _ z (jf) |
|||
(Г Г ѵ) / о |
[ |
/(г + ѵ)2 — Z (Jf) — v] |
(3. 68) |
|
|
|
|
|
У (Г + v)2 — Z (X)
Преобразуем 6 — функцию, стоящую под знаком интеграла, по следующему правилу [91]
Mg' (х,г, 5)]= 8 (г - X) dg
дг
где r = x(x, 1) является корнем уравнения g(x, г, £)=0. В рас сматриваемом случае
g (X, г, І )= Ѵ ( г + v)s— Z (х) + Z (6)— V — гкр (5),
X(x, £)= Ѵ[гкр(Ю~і~ѵ]2-і~^ (x) —Z |
— |
(3. 69) |
|
Теперь окончательно можно записать |
|
|
|
f ( x , r ) = {r + v) /o ^ (r + v)2— -Z -(x)— |
M-----\-{ — |
b ( r - y ) d i . |
|
V(r + v)2 — z {X) |
І, |
|
(3.70) |
Если граничное условие таково, что при х = Хі капли отсутст вуют, то /о=0 и уравнение (3.70) упрощается. Вычисление f(x, г) по формуле (3.70) обычно проводится без использования опе рации интегрирования. Воспользуемся основным свойством функции Дирака [91]
J |
< 8(р |
(& е |
) - |
* |
) |
< я = |
< р |
( * ) |
и, считая для простоты /о= 0, |
запишем уравнение |
(3.70) |
в виде |
|||||
|
/[еос.*)] |
|
|
8 (г - 1) |
dy, |
|
||
|
ea [S(X-*)] “ [£(*.*)] |
ді |
|
|||||
|
|
|
||||||
|
|
|
|
|
|
öS |
|
|
или |
|
|
|
|
|
|
|
|
f i x, r) = |
|
/ ( S o ) |
|
f / |
d x |
\ |
|
|
e 2 |
( S o ) « |
( £ o ) |
LI ö S |
A=:e0 |
|
|
||
|
|
|
п о
где £о(/\ X) является корнем уравнения х — г>т- е-
(Г+ |
V )2 - [гкр (£ ) + ѵ]а - |
Z (X ) + |
Z |
(So) = |
О , |
||
а производная от % по | |
находится из уравнения |
(3.69) |
|||||
дх |
1 |
|
drKpг кр (£) |
_L d Z & |
|
||
öS |
Х+ ѵ |
( г кр + Ь ) |
• ~df |
2 |
di |
|
|
Таким образом функция f(x, |
г) |
зависит лишь от одного па |
раметра |о. определяемого двумя независимыми переменными х и г. Величина go представляет собой значение х, при котором произошло образование зародыша, имеющего для рассматривае мого X радиус г.
Чтобы получить формулу для массовой концентрации кон денсата, подставим выражение для функции распределения
(3.70) при /о=0 в уравнение (3.51)
00 X
|
6л |
_ |
4 |
JtQ B [ |
Г ъй Г Г |
— |
S ( r — |
x ) d |
i |
|
|
6 s |
|
3 |
i |
J |
6 |
|
|
|
|
|
|
|
|
Ч р |
|
|
|
|
|
|
|
|
|
|
X |
00 |
|
|
|
|
|
|
|
4 |
|
в f — |
d t \ |
r 38 ( r - x ) f l l r , |
|
|
||
|
|
Т3 |
|
|
|
|||||
|
= |
Я 6 |
J 6 s « |
•' |
|
|
|
|
|
|
или |
|
|
|
|
'к р |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< Ч = -^“ Я6в\ ---- (3.71) |
|
|
|
|
|||||
|
|
3 |
J |
es (5) “ (5) |
|
|
|
|
||
|
|
|
|
X i |
|
|
|
|
|
|
Выражение для скорости изменения массы |
конденсата по х |
|||||||||
молено получить подстановкой (3.62) |
и (3.70) |
при /0 = 0 в соотно |
||||||||
шение (3.59). В результате подстановки получим |
|
|
||||||||
das |
|
I |
|
|
X |
|
/(£) |
|
l2 |
|
|
-2nQBZ'(x) |
|
|
d l (3. 72) |
||||||
дх |
QSU |
6s (£)И (5) |
V+ 1 |
■ X i
В левой части этого уравнения стоит частная производная по х, потому что, как уже отмечалось, можно пользоваться всеми формулами данного параграфа для расчета как одномерных, так и двумерных течений. В первом случае частную производную следует заменить полной, так как as = as(x). Для плоских и осе симметричных течений as= as{x, ф) и соотношение (3.72), как и все остальные в этом параграфе, выполняется вдоль любой ли нии тока ф = const.
Свободномолекулярный режим роста капель. Рассмотрим вто рой практически важный случай, когда основное кинетическое уравнение (3.60) может быть проинтегрировано. Это можно сде
111