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

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

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

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

Добавлен: 27.06.2024

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

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

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

Re = R/ (\ + W) получим Qea = Tea = pn= \\

 

 

 

ueи MeHУ у-g, Mg H-

 

=M,

 

J + Ю -

 

уу-gRgTg

 

V'

 

 

*'

 

Потери из-за двухфазности определяются формулой

Іе — І

д ъ = - Ѵ - .

где Іе — импульс равновесного течения.

Как показывают расчеты, выполненные для указанных выше исходных данных, при W = 3 в конце сопла (ж = 12) Д<р5~0,04, а

при W = 0,5 A<ps= 0,02.

Г л а в а VI

ВАРИАЦИОННЫЕ ЗАДАЧИ

Значительный практический интерес представляет собой за­ дача определения формы канала (сопла), обеспечивающей наи­ большую тягу при двухфазном течении с отставанием частиц. Оптимизации сверхзвуковой части осесимметричного сопла при движении в нем потока газа посвящен ряд работ [43, 135, 161, ПО, 111, 42, 71, 70], однако, решение такой вариационной за­ дачи с учетом сужающейся части сопла и двумерности течения, даже при движении газа без частиц, в настоящее время отсутст­ вует. Естественно, что и при движении газа с частицами в такой постановке вариационная задача не решена. Однако, если отка­ заться от учета двумерности течения или формы дозвуковой час­ ти сопла, то требуемые решения могут быть найдены.

Простейший случай — одномерные линеаризованные тече­ ния — был рассмотрен в работах [89 и 112]. Ф. Марбл [89] опре­ делял контур сопла, обеспечивающий максимальную скорость истечения, а Л. Е. Стернин [112] — максимальную тягу. Анало­ гичная задача для нелинеаризованного одномерного потока ре­ шена А. Н. Крайко и др. [74].

В работах {89, 112, 74] наряду со сверхзвуковой частью сопла определялась и дозвуковая, сужающаяся часть. Поскольку в ра­ ботах [89 и 112]) потери из-за рассеяния (непараллельности и не­ равномерности потока) не учитывались, то полученные экстре­ мальные контуры имели весьма большие углы конусности сопел

ввыходных сечениях. Пути уменьшения этих углов рассмотрены

вработах [112 и 74].

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

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

178


Однако в ней не были выписаны условия на линии

раздела.

В более полном виде эта же задача решена в работе

[75]. При

этом, наряду с получением условий на линии раздела,

в работе

[75] были исследованы более общие конфигурации экстремаль­ ных сопел, в частности, имеющие внутренние точки излома. Сле­ дует, однако, отметить, что в работах [39 и 75] из-за чрезмерной сложности вычислений результаты исследований до числовых значений не доведены.

Ниже приведено изложение результатов исследований, вы­ полненных в работах [112, 74, 75].

§ 6. 1. Линеаризованное одномерное течение

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

Сила, действующая на внутренний контур сопла, или тяга, создаваемая соплом в пустоте,

(6 . 1)

где давление р определяется формулой (1.46).

Будем искать контур сопла при условии, что его максимально допустимая длина ограничена, т. е. при у — у а, где X — заранее заданная величина. Поскольку в соотношение (1.46) х явно не входит, то за независимое перемен­ ное удобно принять у, а х считать искомой функцией.

Составим функционал, пропорциональный тяге Р

У (1 + 1 Г ) ( 1 - Х 2 ) ( х ~ 1 ) :

 

di,

2

укН (к) dk/dy

 

 

dy

Т2 (к) x'dy/dX

 

и в соответствий с формулой (1.46)

ру заменено

выражением

(Тц+ У'г), а

условие /7з= 0 — дифференциальная

связь, вводимая с переменным множи­

телем Лагранжа р(у)\ кроме того,

в настоящем параграфе величины Я, Т(к),

Я(Х), х и у соответствуют равновесному потоку,

причем для упрощения за­

писей индекс «е» у этих величин опущен. Здесь и ниже в данном

параграфе

штрих означает дифференцирование по у.

 

 

Первая вариация функционала Ф после интегрирования по частям полу­ чает вид

179


d_

ѢФ =

dy

ун

dF

+

dx'

Множитель Лагранжа вале Уи ^ У^ Уа

 

К

 

ÖF

+ Р'

 

dF3

 

dy

dx

 

b x \ d y

 

 

дх'

 

+

dF,

 

 

dF3

 

l^a

(6. 2)

bx

и ------- - b i

 

У-'dx' j

\yn

^

dV

vК '

 

у {у) можно

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

интер­

 

 

 

 

р

 

d

(

dFz \

п

 

 

(6. 3)

 

 

 

 

F-) —---I [

j , ■] = 0,

 

 

а в точке у= уа р = 0 .

 

 

 

dy

 

di'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При у = У и х = £ = 0 , и, следовательно, бх = б £ = 0 .

хн< х < х а можем

Вследствие произвольности вариации бх на интервале

написать

 

 

 

 

E i

+

II dF3\ =

0.

 

(6. 4)

 

 

 

 

d

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

dy

 

дх'

 

дх'

 

 

 

Определим вначале функцию р(г/). Из уравнения расхода, записанного в

виде

 

 

 

1

 

 

 

 

 

 

 

 

следует

'* +

1

 

1

X — 1

 

Ü L (

 

 

 

2

 

Г 1*

 

X + 1

 

У2

 

(6. 5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

__ у (I2 — 1)

 

 

 

(6. 6)

 

 

 

 

d l

 

~

21Т(1)

'

 

 

 

 

 

 

 

 

 

 

Принимая во вінимание, что

â/r3/(3^'= 1, подставляя в

выражение

(6. 3)

значение F2 и используя формулу (6.6),

получим

 

 

 

 

dy

_

 

П^Д)У^.(1

+ уХ2)

/

2

 

 

 

d l

~

2(1 +

 

1Г)(х — 1)Х2

и

+ 1/

 

 

где ро — давление торможения равновесного потока.

т. е. при у= уа, и

Учитывая краевое условие в концевой точке контура,

интегрируя последнее дифференциальное уравнение, получим

 

К- [>• («/)]

^РоУУ

 

 

 

 

Т~ — “Г +

V (>• — *а)

(6. 7>

2 (1 +

W) (* — 1) \х +

1

 

 

 

 

 

причем связь между Х н у

 

определяется равенством (6.5).

 

 

Подставим

полученное согласно

выражению

(6.7) значение р в условие

(6.4). После преобразований придем к соотношению

 

 

 

 

 

 

ІН (к) I

d l

у

 

 

 

(6. 8)

 

 

 

 

П (I)

U *

const.

 

 

 

 

 

 

/

 

 

 

 

Поскольку допустимая

вариация

бхо^ 0 ,

ра = 0, а величина (dFі/дх')а

положительна,

то, как следует из уравнения (6.2),

вариация бФ^О. Поэтому

наивыгоднейшее сопло имеет длину х—Х.

 

 

 

 

Интегрируя выражение (6.8), получим

 

 

 

 

 

х

 

 

 

 

 

 

 

 

 

 

 

х = Х ^ / ІН (I)

 

[7’ (X)]-irfX

/ Х Я (X) [7-(Х)] -Irfx

(6. 9)

 

 

 

 

 

 

 

 

Xн

 

 

 

 

180


Обе постоянные, входящие в уравнение (6.9), получены из условий: при

«/=(/„ x=0, при у —у а

х=Х.

 

 

Квадратура (6.9)

может быть выражена через эллиптические интегралы

1, 2 и 3 рода, однако,

на практике целесообразнее

выполнять численное

ин­

тегрирование.

 

сопла выполняется

по

Таким образом, расчет экстремального контура

двум достаточно простым формулам — (6.5) и (6.9). Интересно отметить, что

на координаты контура не влияет величина е (или размер частиц),

влияние

W осуществляется только через х по формуле (1.41), что характерно для ли­

неаризованной теории.

 

Необходимо отметить, что при составлении вариации Ф величина \ пола­

галась постоянной. Изменение %, как следует из соотношения (1.52),

может

иметь место при изменении (dX/dx)*, т. е. в одной точке. Таким образом, тя­ гу можно увеличить, меняя наклон контура лишь в одной точке. Это является

недостатком данной линеаризованной теории.

 

когда

квадратура

(6.9)

 

Остановимся теперь на одном частном случае,

может быть выражена через элементарные функции.

Пусть

Н(Х) =

1,

что со­

ответствует вполне реальному случаю, когда

 

 

 

 

 

 

 

 

 

 

 

 

 

З / о

Рг Сд (1 + W )

J

 

 

 

 

 

 

 

 

71 ~

Nu Cp (Cp + W cB)

2

 

 

 

 

 

 

При этом из квадратуры

(6.9) получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

x = X L ( k ) L J 1,

 

 

 

 

 

(6.10)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L ß ) .

-3/4

 

1

|п

(і + угууЛ0 0 У^н)

 

 

 

 

 

 

 

 

2 П (і - Ѵ ѵтА М і + Ѵ ѵ/ х;)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

_

4

_1

 

 

 

 

 

 

 

 

— arctg / у

/ X

+ arctg /

у /X „ J ;

 

 

 

 

 

 

 

 

 

 

 

 

La =

L (Xa).

 

 

 

 

 

 

 

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

(6.10),

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

X

X H ( X ) d X

Lg

 

 

 

 

 

 

 

 

: (X) =

2y 5

L(l),

 

 

 

( 6. 11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ТЧ (к) d x / d X

X

 

 

 

 

 

а формула (1.46) примет вид

 

 

 

 

 

 

 

 

 

 

eW

 

2

 

 

 

 

 

 

 

 

IX

 

=

1 +

 

*

 

 

 

 

X3/2 — (1 + уХ2) L (X)

 

 

 

 

 

 

 

W a l

Ре

1 + W Х + 1 (1 — Х2)Х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6. 12)

 

Аналогично этому могут быть преобразованы соотношения

(1.47)

и (1.48).

 

Выражение для \

(1.52)

можно записать в виде

 

 

 

 

 

 

 

 

 

£ =

2уі0ЛГ-і[і(1) — (% + 1)/2х].

 

 

 

(6.13)

 

Разлагая (6.12)

в ряд по степеням Д = 1 —X в окрестности

точки Я=1 и

переходя

к пределу при А

 

*0, получим после преобразований

 

 

 

 

 

-Ре)J * = > -

 

eWy.L„

 

1 + 4у — у'

 

 

 

 

1 + «7(х +

1) X

 

1 — у2

 

 

 

 

Эта формула дополняет выражение (6.12),

поскольку при Х =1

расчет по

формуле (6.12) проводить нельзя.

 

 

 

 

 

 

 

7

3739

 

 

 

 

 

 

 

 

 

 

 

 

181


Определим тягу, создаваемую экстремальным соплом. Безразмерный па­

раметр Р, пропорциональный тяге Р, создаваемой соплом, после интегрирова­ ния и преобразований может быть представлен в форме

Р

Р =

sWLa*-

X 1

Т(К)

Т (кд) 1

Т (Ха)

 

(1 + Г )Х (х + 1)

1 (1)-

х„

J

Ха La

 

 

 

 

 

(6.

14)

Первые члены этой формулы, заключенные в круглые скобки, характери­ зуют изменение импульса равновесного течения, и остальные — влияние от­

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

В работе [112] показано, что если контур сопла представить в виде

 

 

 

 

X = XLaL (X) +

8 (Хві -

X) (X - ХНі) ,

 

 

 

(6. 15)

-где

Хн .< ХНі ■< X

ХДі -< Ха ,

 

 

 

 

 

 

 

 

 

б — малая величина, по формуле

(6.15) найти dxjdX,

подставить в форму­

лу

(1.46),

а затем

определить Р, то придем

к формуле,

отличающейся от

формулы

(6.14)

лишь наличием одного слагаемого / = б27,

стоящего в фигур­

ной

скобке. При

этом J — всегда положительная

иррациональная

функция

от ЯН1 и Хаі. Поскольку перед фигурной скобкой

стоит

знак

минус,

а б —

действительное число, то это означает, что отклонение от найденного

конту­

ра в любую сторону ведет к потере тяги. Следовательно,

 

найденное решение

отвечает максимуму тяги.

 

 

 

 

 

 

 

 

 

Следует отметить, что формула (6.14) получена при условии (1.52) или

(6.13), т. е. при рассмотрении течения во всем сопле. Если

же течение

рас­

сматривается только, например, в расширяющейся части сопла,

то | вместо

условия (1.52) вычисляется по формуле (1.44) и соотношение

(6.14)

видоиз­

меняется.

 

 

 

 

 

 

 

 

 

 

 

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

тодике,

для следующих исходных данных: rj= 1 /2

и >с= 1,2.

Таблица 6.1

 

 

Параметры экстремальных сопел

 

 

 

 

X

УІѴ*

РІРе

 

Т/Те

wj w e

s WHI + W)

0,2

0,8

0,2

0,8

0,2

0,8

1,1

1,005

0,9977

0,9907

1,0034

1,0134

1,0057

1,0227

2,2

2,264

1,0133

1,0530

1,0133

1,0534

1,0001

1,0003

2,5

4,067

1,0206

1,0822

1,0194

1,0775

0,9988

0,9953

Из рассмотрения табл. 6.1 следует, что отличие первого приближения от равновесного течения невелико. Это свидетельствует о целесообразности ис­ пользования линеаризованной модели для многих интересных в практике за­ дач.

182