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

 

_

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