
- •ОГЛАВЛЕНИЕ
- •ВВЕДЕНИЕ
- •РАЗДЕЛ 1. УРАВНЕНИЯ ГАЗОВОЙ ДИНАМИКИ ДЛЯ ЭЛЕМЕНТАРНОЙ СТРУЙКИ – ПОВТОРЕНИЕ ПРОЙДЕННОГО
- •1.1 Используемый математический аппарат
- •1.2 Уравнение неразрывности
- •1.3 Уравнение количества движения
- •1.4 Уравнение энергии
- •1.5 Предельная скорость движения газа. Число Маха и коэффициент скорости
- •РАЗДЕЛ 2. ОДНОМЕРНЫЕ ГАЗОВЫЕ ПОТОКИ
- •2.1 Звуковые волны. Скорость звука. Излучение звука
- •2.2 Волны конечной интенсивности. Инварианты Римана. Характеристики
- •2.3 Механизм образования скачков уплотнения
- •2.4 Прямые скачки уплотнения. Ударная адиабата
- •2.5 Скорость распространения ударной волны и спутного потока за ней
- •2.6 Элементарная теория ударной трубы
- •2.7 Косые скачки уплотнения
- •2.8 Теория мелкой воды
- •2.9 Волны горения и детонации в газах
- •2.10 Общие условия перехода от дозвукового течения к сверхзвуковому
- •РАЗДЕЛ 3. ПЛОСКИЕ ТЕЧЕНИЯ НЕВЯЗКОГО ГАЗА
- •3.1 Общие уравнения
- •3.2 Метод малых возмущений
- •3.3 Дозвуковые течения при малых возмущениях
- •3.5 Обтекание малого угла сверхзвуковым потоком
- •3.6 Обтекание тонкого профиля сверхзвуковым потоком газа
- •3.7 Волны разрежения. Центрированные волны
- •3.8 Общая задача о двумерном стационарном движении газа. Уравнение Чаплыгина
- •РАЗДЕЛ 4. ТЕЧЕНИЕ ГАЗА В СОПЛАХ И ДИФФУЗОРАХ
- •4.1 Формы сопел. Простейшая теория сопла. Истечение сверхзвуковой газовой струи из сопла на нерасчетном режиме
- •4.2 Дозвуковые и сверхзвуковые диффузоры
- •РАЗДЕЛ 5. ГАЗОВАЯ ДИНАМИКА ЭЛЕМЕНТОВ ТУРБОМАШИН
- •5.1 Обтекание конечных тел. Ламинарный и турбулентный след
- •5.2 Дозвуковое обтекание тонкого крыла. Формула Жуковского
- •5.3 Обтекание решетки профилей газа
- •5.4 Простейший расчет ступени лопаточной турбомашины
- •ЛИТЕРАТУРА
- •ПРИЛОЖЕНИЯ. ЗАДАЧИ, РЕШАЕМЫЕ НА СЕМИНАРАХ И В РАСЧЕТНОМ ЗАДАНИИ
- •П1. Гидравлическое сопротивление трубы при ее внезапном расширении
- •П2. Расчет простейшего эжектора
- •П3. Пример расчета аэродинамического нагрева
- •П4. Гидравлический удар при внезапной остановке потока жидкости
- •П5. Трубка Пито при сверхзвуковом движении газа
- •П6. Прямоточный воздушно-реактивный двигатель (ПВРД)
- •П.7. Одномерное стационарное течение газа по трубе постоянного сечения с подогревом
- •П.8. Адиабатическое течение газа с трением в трубе постоянного сечения
- •П.9. Варианты расчетного задания по курсу "Техническая газодинамика"
- •Вариант 1. Крыловой профиль в гиперзвуковом потоке
- •Вариант 2. Ракета А-4
- •Вариант 3. Снаряд

9
Рисунок 1.3.1. Эпюра скорости вязкого потока у стенки
В случае идеальной жидкости (вязкость μ = 0, плотность ρ = const) уравнение (1.3.1) переходит в уравнение Эйлера
|
|
|
|
ρ |
D W |
|
(1.3.2) |
Dt |
=ρ F − P . |
В последующем мы достаточно часто будем использовать уравнение движения в простейшей форме для стационарного течения в канале постоянного сечения S
(P1−P2 )S=G(W 2−W 1)+F friction+F mech . |
(1.3.3) |
Здесь F friction , F mech – соответственно сила трения, действующая на боковые поверхности канала и сила, с которой газ действует вдоль
направления своего движения на машину внутри канала (лопатки компрессора и т. п.), получающую от газа техническую работу.
При отсутствии трения и совершения технической работы (1.3.3) сводится к соотношению
P2+ρ 2 W 22=P1+ρ1 W 12=const . |
(1.3.4) |
1.4Уравнение энергии
Всамом общем виде данное уравнение (закон сохранения энергии) записывается в следующей дифференциальной форме
|
|
|
|
|
10 |
|
|
|
|
2 |
|
|
|
|
|
|
D (u+ |
W |
) |
|
|
. |
(1.4.1) |
|
2 |
|
|
||||
ρ |
|
|
|
=ρ F W + (̂τ W )+ρ q |
|
|
|
D t |
|
|
|
||||
|
|
|
|
|
|
Здесь u – удельная внутренняя энергия среды; q – массовая плотность источников тепла.
Однако использовать уравнение (1.4.1) в самом общем виде в данном курсе мы будем весьма редко. Значительно чаще нам потребуется уравнение энергии в стационарном случае. Последнее записывается, обычно для элементарной струйки, рис. 1.2.1. Для входного (1 – 1') и выходного (2 – 2') сечений элементарной струйки, по которой за время dt проходит массовый расход газа dG=ρ1 W 1 S1 dt=ρ 2 W 2 S2 dt=idem ,
общий баланс энергии в стационарном |
случае может быть записан в виде |
dA+dQ=dE K +dE P+dET +dl +dl fric . |
(1.4.2) |
Здесь dA – работа, совершаемая над газом; dQ – количество тепла, переданного газу извне; dl – механическая работа, совершаемая газом; dlfric – работа сил трения; dEK – изменение кинетической энергии; dEP – изменение потенциальной энергии, если течение происходит в поле массовых сил; dET – изменение внутренней энергии.
Изменение кинетической энергии:
|
(W 2−W 2 ) |
. |
|
dEK =dG |
2 1 |
||
2 |
|||
|
|
Изменение потенциальной энергии:
dEP=dG g ( z2−z1 ) .
Изменение внутренней энергии в случае постоянной изохорной теплоемкости CV:
dET =dG(u2−u1)=dG CV (T 2−T 1) .
Работа, совершаемая над газом:
dA=A2+A1 .
Здесь A1 – работа сжатия, совершаемая над газом при входе в трубку тока. Она положительна, так как газ сжимается,

11
A1=P1 dG .
ρ 1
A2 – работа расширения, совершаемая газом при выходе из трубки тока. Она отрицательна, так как газ расталкивает окружающую среду,
A2 =−P2 dG . |
|
|
|
|
|
|
|
|||||||
|
|
ρ 2 |
|
|
|
|
|
|
|
|
||||
Таким образом |
|
|
|
|
|
|
|
|
||||||
dA=dG( |
|
P1 |
− |
P2 |
) . |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|||||||
|
|
ρ 1 |
ρ 2 |
|
|
|
|
|||||||
В таком случае (1.4.2) можно записать в виде |
|
|||||||||||||
|
|
P |
1 |
|
|
P |
2 |
|
|
(W 2−W 2) |
|
|
||
dQ+dG( |
|
− |
|
)=dl +dl fric |
+dG |
2 1 |
+dG g (z |
2−z1)+dG (u2−u1) . |
||||||
|
|
|
|
2 |
||||||||||
|
|
ρ 1 |
ρ 2 |
|
|
|
Разделив левую и правую части на dG получим соотношение, связывающее уже удельные величины (то есть отнесенные к 1 кг вещества)
|
P |
1 |
|
P |
2 |
|
|
W |
2 |
|
W |
2 |
|
Q+ |
|
− |
|
=L+L fric |
+ |
|
2 |
− |
|
1 |
+g ( z2−z1)+u2 −u1 . |
||
|
|
|
|
2 |
|
2 |
|
||||||
|
ρ1 |
ρ2 |
|
|
|
|
|
Выделение тепла Q в движущейся среде происходит вследствие двух различных процессов. Во-первых, некое количество тепла Qout выделяется благодаря действию внешних факторов (нагрев стенок канала, протекание электрического тока и другие). Во-вторых, вследствие наличия трения, которое совершает работу Lfric, часть механической энергии переходит в
тепловую Qfric = Lfric.
Учитывая, что удельная энтальпия вещества определяется как
h=u+ Pρ , а для газа с постоянной теплоемкостью (наиболее часто рассматриваемый нами случай)
RT = Pρ , u=CV T ,
12
где R – газовая постоянная, C P=CV +R (формула Майера), уравнение
энергии в стационарном случае принимает вид |
|
||||||
|
W |
2 |
|
W |
2 |
+g ( z2−z1) . |
(1.4.3) |
Qout −L=h2−h1+ |
|
2 |
− |
|
1 |
||
2 |
|
2 |
|
||||
|
|
|
|
|
|
Для энергетически изолированных течений, когда среда не совершает механической работы и к ней не подводится тепло, а изменение потенциальной энергии пренебрежимо мало (в газовых машинах это весьма частый случай), уравнение энергии принимает наиболее простой вид
|
W |
2 |
|
W |
2 |
=h0=const . |
(1.4.4) |
h2+ |
|
2 |
=h1+ |
|
1 |
||
2 |
|
2 |
|
||||
|
|
|
|
|
|
Величина h0 называется полной энтальпией или, иначе, энтальпией полного торможения, так как соответствует термодинамической энтальпии газа после его полного адиабатического торможения. Соответствующая этому температура T0 , которую принимает газ, называется температурой полного торможения.
Легко видеть, что для случая несжимаемой среды ( ρ =const ), когда при течении жидкости внутренняя энергия не меняет своего значения,
закон сохранения полной энтальпии (1.4.4) сводится к известному интегралу Бернулли
|
ρ W 2 |
|
ρ W 2 |
|
(1.4.5) |
|
P2+ |
2 |
=P1+ |
1 |
=const . |
||
2 |
2 |
|||||
|
|
|
|
1.5 Предельная скорость движения газа. Число Маха и коэффициент скорости
Как уже было получено в предыдущем параграфе, для газового потока при отсутствии энергетического обмена с окружающей средой,
h0=h+ |
W 2 |
=const . |
(1.5.1) |
|
|||
2 |
|
|
В таком случае, если газ истекает через сопло, например, в вакуум, из объема, где он покоится и имеет температуру T0 , то можно определить ту максимальную скорость, которую он способен приобрести. Так как

13
давление и плотность газа в вакууме равны нулю, то газ, истекая в вакуум, адиабатически расширяется, а его температура падает до абсолютного нуля. Таким образом в вакууме для газа P = 0, T = 0, ρ = 0, а значит, и энтальпия h = 0. В соответствии с (1.5.1) получаем максимально возможную скорость истечения газа в вакуум
W max=√ |
2 h0 |
. |
(1.5.2) |
Для воздуха, так как у него CP = 1000 Дж/(кг К) = const, формула упрощается,
W max=44.8 √T 0 .
Для случая T0 = 300 К получаем Wmax = 776 м/с.
Для произвольного сечения сопла, в котором термодинамическая энтальпия газа равна h,
h0−h |
= |
W 2 |
. |
(1.5.3) |
|
h |
2 h |
||||
|
|
|
В случае идеального газа (постоянная теплоемкость), воспользовавшись
формулой Майера |
C P=CV +R , вводя понятия показателя |
адиабаты |
|||||||
|
C P |
|
|
|
|
M = |
W |
|
|
k= |
, скорости звука as=√k RT , числа Маха |
(об этих |
|||||||
CV |
as |
||||||||
|
|
|
|
|
|
|
величинах подробно будет говориться в следующей главе), вместо (1.5.3) получаем
T 0−T =(k−1) M 2 .
T 2
Отсюда получаем весьма полезную формулу, позволяющую найти температуру полного торможения газа
T 0=T (1+ |
(k−1) |
M 2) . |
(1.5.4) |
|
2 |
||||
|
|
|
Число Маха M выступает в данном случае в виде безразмерной скорости газового потока. В качестве масштаба скорости as используется скорость звука (скорость распространения слабых колебаний давления) в газе при текущей его температуре.

14
В случае M < 1 течение газа называется дозвуковым, M > 1 – сверхзвуковым, а при M ~ 1 течение именуется трансзвуковым. Следует отметить, что выбор в качестве масштаба скорости звука в газе при текущей его температуре далеко не всегда удачен. Например, в задаче об истечении газа в вакуум температура газа на выходе из сопла близка к нулю, T → 0, а значит, что as → 0, в то время как M → ∞.
В такого рода задачах удобнее использовать масштаб скорости, который не зависит от конкретного рассматриваемого сечения. В качестве такового используют скорость звука в газе в критическом сечении сопла. Критическим называется сечение, в котором скорость газа в точности равна скорости звука в нем, в этом сечении M = 1. Температура в этом сечении, в соответствии с (1.5.4), равна
T crit = |
2 |
T 0 . |
(1.5.5) |
|
(k +1) |
||||
|
|
|
Скорость звука в критическом сечении acrit , определяемая через скорость звука при параметрах торможения a0=√k R T 0 ,
acrit |
=√ |
|
T crit |
|
|
2 |
|
|
|
|
|
=√ |
|
. |
|||||
a0 |
T 0 |
(k +1) |
Формулу (1.5.2) можно теперь записать в виде
W max=a0 √k 2 1 .
−
Обезразмеренная на величину acrit скорость газового потока называется
коэффициентом скорости λ,
λ = |
W |
. |
(1.5.6) |
|
|||
|
acrit |
|
Для максимальной скорости истечения коэффициент скорости, в отличие от числа Маха, конечен,
λ max= |
W max |
= |
|
|
k+1 |
. |
|
acrit |
√k−1 |
||||||
|
|
|
Связь между числом Маха и коэффициентом скорости взаимно однозначна