Файл: Стернин Л.Е. Основы газодинамики двухфазных течений в соплах.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 — VУ У^н) |
|
|
|
|||||
|
|
|
|
|
2 П (і - Ѵ ѵтА М і + Ѵ ѵ/ х;) |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
4 |
_ |
4 |
— |
_1 |
|
|
|
|
|
|
|
|
— arctg / у |
/ X |
+ arctg / |
у /X „ J ; |
|
|
|
|
||||
|
|
|
|
|
|
|
|
La = |
L (Xa). |
|
|
|
|
|
|
|
Используя равенство |
(6.10), |
получим |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
X |
X H ( X ) d X |
2у 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)- 2х |
х„ |
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