Файл: Стернин Л.Е. Основы газодинамики двухфазных течений в соплах.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 _JL„l 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 = |
||||
= (usu—usv) (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 |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
6а |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
— + -^г(— ) = °. |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
6а |
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 = |
Q£ + 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