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

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

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

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

Добавлен: 27.06.2024

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

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

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

Из уравнения (4.62) следует, что характеристиками являют­ ся линии тока газа (Л = 0 ) и линии тока частиц (ß = 0 ), что бы­ ло уже получено выше, а также линии, удовлетворяющие урав­ нению

В2[(1+/ V - Л2]+e^Q-1erV(і+г/'а) л2=о. (4.63)

Условие совместности, выполняющееся вдоль характеристик

(4.63^, имеет вид

 

 

 

 

-В— (а?уг Аи) ~\-а2ѵи' Jr {aiy''‘ — A2)uv' -j- А

N^j —

-(a*y' + Au)(Nt + y 'N a)-

ег05а2,42

В (vas uvs) qs

 

 

 

 

 

QsQB2

 

 

+A+Psy — y' { f X+b sx)-}WcB

е;е«д2л2 у

p

 

 

qIq&

 

где

+ (1 +

/ а)(Мг)'- І Ѵ з ) ] ^ 0 ,

(4. 64)

 

 

 

 

 

f'7

p,+?l f

- Q1- ?v;+ f 7P-

N*

Qs / , , N a p _JLl Jf y

 

Соотношение (4.63) является уравнением четвертой степени относительно производной у' и для анализа его решений разло­ жим корни этого уравнения по степеням отношения рг/рв. (Это отношение, как отмечалось выше, является весьма малой вели­ чиной, т. е. 10~3-=-10“5) . Представим у' в виде

і /'= ^ + ^ ( е г / е в ) + ^ ( е гА>в)2+ - • •

(4.65)

Подставим (4.65) в (4.63) и приравняем коэффициенты, не содержащие рг/рв, получим уравнение четвертой степени, из ко­ торого следует

УIО 1,2

иѵ ± а Y V 2 — а2

Vs

 

 

и2а2

УО, 3,4

tlS

 

 

 

 

 

коэф-

Рассмотрим сначала два последних корня. Приравнивая

фициенты при рг/р в и определяя Yu получим

.,2'

 

 

 

 

 

 

 

 

У3.4

V s

ОгЛѴ-« (V U S — U V S) Y Q s / Q

I Qi

 

 

U s

QrusV ( V U S uvs)2 V2a2

e2

 

Введем угол

 

б между векторами V и Vs по формуле

sin 6 =

= (usuusv) (Ws)-*. Тогда последняя формула примет вид

Vs

6 га У ^

sin

5

Qs

 

 

'3,4"

 

Qbm5

Y

Q(V2sin2 P ■

 

 

Us

 

 

 

 

 

 

 

142


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

В этом случае приравнивание коэффициентов при дг/дв при­

водит к условию Yі = 0, а при

( р гД > в ) 2 — к соотношению

 

ofesiv — y'o1)2и)4

(4.66)

У1,2 У01,2 ±

 

2е20а У Ѵ 2 — а2 у'о1;2Us)2

 

Отсюда следует, что эта

пара корней действительна при

Ѵ > а и мнима при Ѵ<.а. Полученные два семейства характери­ стик аналогичны линиям Маха (характеристикам первого и вто­ рого семейств) обычной газовой динамики.

При ц-ѵа производная у'о 2->-°о и вместо последнего равен­ ства следует использовать разложение

X.1 , 2 '

 

(-^0 1,2^ — “Ѵ*

I

(4. 67)

ѵ0 1,2

- О

 

2ö2Qö У V2

а2 ( х'0 12 v s us)2

 

 

где ѵ0 1,2

иѵ Т а У V2 — а2

 

 

 

V2а2

 

 

 

причем, если раньше штрихом обозначалась производная

по х,

то в последних двух равенствах штрих соответствует производ­ ным по у. Естественно, что корни х / и х2' действительны здесь также при Ѵ>а.

Подставляя первое слагаемое разложения (4.66) для треть­ его и четвертого корней в условие совместности (4.64), получим с точностью до QtIqb включительно

+

Ѵ Ѵ 2 — а2

p' +

У

+ Аі,2

ѴѴ

±

QT(Z

 

У

У Ѵ 2 а 2

 

 

,

. . . .

6 г Ц И і ,2 Г o f Us У I

+ -

---------- (*ш + УОІ,>*Г.) +

, In2

[^ (^ 7 )

+

 

 

 

 

 

 

QbQBIz

 

Bi,2(vus uvs)

-о; + л

+ я ,у- ^ іі2 (л + р „ ) + ^ ^ - 1 = о,

, n

^ i,2 Q i

 

 

 

 

 

 

У J

(4. 68)

где

Соотношения (4.66) и (4.68) являются равенствами, которые характеризуют изменение параметров вдоль характеристик с точностью до Рг/qb включительно.

143


Граница области существования действительных линий Ма­ ха, определяемая из разложений (4.66) и (4.67), не является точной, так как при Ѵ= а эти разложения несправедливы.

Разложим корни (4.63) по степеням sin б. (Этот угол во мно­ гих интересных для практики задачах сравнительно невелик).

Прежде всего перейдем от производной у' к параметру ^ = tga, где a — положительный или отрацительный угол между скоростью V и характеристикой. Очевидно, что

V + tu

(4. 69)

У’ = а tv

Подставляя у' из равенства (4.69) в выражения для А и В и учитывая, что vsu—wsu= KKssin б, получим

л ._

W

.

а

t (W 5) - Г Г , sin о

 

 

tv u

 

tv — u

где (VV's) — скалярное произведение векторов V и Vs-

Используя равенство

(4.69)

и последние равенства, вместо

уравнения (4.63)

получим

 

 

 

[V*V2s s i n 4 - 2 V V s sin Щ Ѵ Ѵ ^+ Р ІѴ Ѵ ,?] [а%{1 +

62) - ( 2l/2J-f-

 

+ e ^

Q

-

^ r 1 +*2)= 0.

(4. 70)

Будем искать решение этого уравнения в виде разложения t по степеням sin б

sin § —1~. . .

Подставляя это разложение в (4.70) и приравнивая коэффи­ циенты при одинаковых степенях sin б, окончательно получим

t g a ^ = +

4----------------^

---- 1~0 (sin28),

(4.71)

 

У Г2 — *2

£2 (V2 _ ЬЧ) (Ws)

 

где &2= a 2 [1 +Q?esQ rV 1^* (^ V ,)-s].

 

 

Аналогичным образом может быть

разложена по степеням

sin 6 и величина Т = ctg аз,4-

Подставляя в формулу (4.70)

t= l/T

и находя величины Т0 и Т\

в разложении 7’ = 7’0+7’1sin б,

придем

к выражению

 

 

 

 

c t g a 3,4 =

у V2 — ЬЧ

(62 _ д2) узуз sin ■ -fO (sin 28).

(4.72)

 

 

ь*(WA

 

 

Итак, границей области гиперболичности является значение Ѵ=Ь. При Ѵ>Ь кроме линий тока имеется два семейства дейст­ вительных характеристик. Как видно из выражений (4.71) и (4.72), характеристики расположены не симметрично относитель­ но линий тока газа, а повернуты в сторону линий тока частиц, хотя поворот этот весьма мал.

144


§4. 6. Односкоростные двумерные течения

Взаключение рассмотрим случай, когда механическое взаимодействие между частицами и газом столь велико, что различием в скоростях можно

пренебречь, т. е. когда во всех точках потока V— Vs, а температуры

Т и Т3

могут быть различными. При этом из уравнений неразрывности (4. 8)

следует

 

 

 

 

 

V(ös^)

<?Qа

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

' ~дГ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ѵѵ

Q

д_

 

 

 

 

 

 

 

 

(4.73)

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— + -^г(— ) = °.

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

'

да /

 

 

 

 

 

 

 

 

где es=e+QS.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Первое

из уравнений

(4.73)

получается

путем

сложения

уравнений

не­

разрывности

(4.8),

а последние два — путем вычитания из первого уравнения

(4.73) соотношений (4.8).

 

 

 

 

 

 

 

 

 

 

 

 

(4.5)

си­

 

Исключая из уравнений импульсов для газа (4.12) и для частиц

лу /,

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-*•

 

оѵ

+

 

 

-*■

° ’

 

 

(4-74)

 

 

 

 

 

Q a ( ^ v )V + Qs

 

ѴР 6а/>а =

 

 

где

Qj,Pa =

qP +

QsPs-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

d

и

исключая

из

Обозначая для краткости QjQ ^ qQ+QjQs, Ѵѵ +

 

dt

уравнений энергии

(4.6)

и (4.15) q,

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_Q_

di___Q

dp_

 

_0£_

des

■<?a =

°-

 

 

 

 

 

 

 

Qa

dt

6röa

dt

 

6a

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая связь между плотностями (4.11), это

соотношение

запишем в

виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_Q_

di_

1

dp

 

Qs

dp

^

de*

 

0,

 

 

 

 

 

Qa

dt

Qs

dt

 

Q2Qb

dt

 

Q^

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как следует из последних равенств

(4.73),

величины

Q/Qa и 6s/Qa

можно вносить под знак dfdt. Обозначая

Qs / a =

+ Qses +

р,

из послед-

него равенства получим уравнение энергии смеси в виде

 

 

 

 

 

 

 

 

 

Ѵуг'а + '

di у.

 

 

 

 

 

- Q a = 0.

 

 

(4. 75)

 

 

 

 

dt

 

6«

\v ™ +

- f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение энергии частиц можно записать в форме

 

 

 

 

 

 

 

 

 

 

Vyes +

des/ât — q Qs =

0.

 

 

 

 

(4.76)

 

Уравнения

(4.73) — (4.76)

описывают движение односкоростной

двухтем­

пературной среды с плотностью

Qs

и энтальпией

7S

для общего случая,

ког­

да объемом частиц не пренебрегается.

5* 3739


Г л а в а V

ДВУМЕРНЫЕ СТАЦИОНАРНЫЕ ТЕЧЕНИЯ МОНОДИСПЕРСНОЙ СРЕДЫ БЕЗ УЧЕТА ОБЪЕМА ЧАСТИЦ. ЧИСЛЕННЫЕ МЕТОДЫ РАСЧЕТА СВЕРХЗВУКОВЫХ ТЕЧЕНИЙ

Численное исследование двумерных и трехмерных течений аэрозолей даже без учета объема частиц является весьма слож­ ной задачей математической физики. Поэтому в 'настоящее вре­ мя точные численные решения имеются лишь для монодисперсных сверхзвуковых потоков. Расчет дозвуковых и трансзвуковых течений до последнего времени производился приближенными методами: или без учета влияния частиц на газ, или при допуще­ нии о равновесности течения, или с использованием каких-либо цругих грубых упрощающих предположений, например, усло­ вия постоянства отставания в трансзвуковой области (см. рабо­ ту [66]).

Первой работой, в которой был приведен пример численного расчета течения газа с частицами в осесимметричном сопле, яв­ ляется (опубликованная в 1961 г.) работа Клигеля и Никерсона [66], в которой выписана система уравнений, включающая урав­ нения характеристик для сверхзвукового течения монодисперс­ ной среды, при условии пренебрежения объемом частиц, и пред­ ставлены результаты расчетов. Исследования сверхзвуковых течений монодисперсмых потоков в осесимметричных соплах чис­

ленными методами с использованием уравнений

характеристик

производились в последние годы рядом авторов [26, 27, 28, 53,

54, 85, 129]. Результаты этих работ, часть которых

будет приве­

дена ниже, значительно прояснили характер течения двухфазной среды в сверхзвуковом сопле.

В работе [118] представлена двумерная линеаризованная тео­ рия плоских и осесимметричных сверхзвуковых течений газа и частиц. Выведено исходное уравнение в частных производных четвертого порядка для потенциала возмущений и дается его об­ щее решение с помощью преобразования Лапласа. Эффектив­ ность теории иллюстрируется на ряде конкретных задач. В рабо­

146