- •ВВЕДЕНИЕ
- •1. ИСХОДНЫЕ ПОНЯТИЯ
- •1.1. О методе размерностей
- •1.2. Некоторые сведения из кинематики жидкости
- •1.3. Дифференцирование по времени интеграла, взятого по объему, который движется вместе с жидкостью
- •1.4. О записи уравнений движения в различных системах координат
- •2. СИСТЕМА УРАВНЕНИЙ ДВИЖЕНИЯ ЖИДКОСТИ КАК СПЛОШНОЙ СРЕДЫ
- •2.1. Уравнение закона сохранения массы
- •2.2. Уравнение количества движения
- •2.3. Уравнение энергии
- •2.4. Уравнения термодинамического состояния
- •2. 5. Вязкие напряжения и теплопроводность
- •2.6. Некоторые дополнительные сведения из термодинамики
- •3. ОСНОВЫ ДИНАМИКИ ИДЕАЛЬНОЙ ЖИДКОСТИ
- •3.1. Интегралы уравнений движения идеальной жидкости
- •3.2. Скорость звука
- •3.3. Различные формы записи уравнения Бернулли
- •3.4. Одномерное установившееся течение газа в канале переменного сечения
- •3.5. Некоторые примеры применения уравнения Бернулли
- •3.6. Установившееся истечение газа из сосуда через отверстие
- •3.7. Неустановившееся истечение газа из сосуда (опорожнение сосуда)
- •4. ВОЛНЫ В СВЕРХЗВУКОВОМ ТЕЧЕНИИ
- •4.1. О качественном различии дозвуковых и сверхзвуковых течений. Течение разрежения. Ударные волны
- •4.2. Течение Прандтля - Майера
- •4.3. Основы теории ударных волн
- •4.5. Конус в сверхзвуковом потоке
- •4.6. Теория Ньютона
- •4.7. Метод характеристик
- •Библиографический список
∂ϕ |
+ gz + |
p |
+ |
v 2 |
= const , |
(3.15) |
|
∂t |
ρ |
2 |
|||||
|
|
|
|
а уравнение неразрывности для несжимаемой жидкости в терминах потенциала скоростей сведется к уравнению Лапласа:
div ( ϕ) ≡ Δϕ ≡ |
∂2ϕ |
+ |
∂2ϕ |
+ |
∂2ϕ |
= 0 . |
(3.16) |
|
∂x 2 |
∂y2 |
∂z2 |
||||||
|
|
|
|
|
В результате из (3.16) можно определить ϕ и затем из (3.15) рассчитать величину давления p в каждой точке потока.
3.2. Скорость звука
Скоростью звука называют скорость распространения малых возмущений. Понятие скорости звука является одним из важней-
ших в газовой динамике. |
|
|
|
||||
|
|
|
|||||
Для |
вычисления |
скорости |
v |
v + dv |
|||
распространения |
малых |
возму- |
|||||
щений рассмотрим равномерный |
p |
p + dp |
|||||
поток, пересекающий слева на- |
|||||||
|
|
||||||
право |
фронт звуковой |
волны |
ρ |
ρ + dρ |
|||
(рис. 3.1), на котором параметры |
|||||||
|
|
||||||
потока претерпевают слабые воз- |
|
|
|||||
мущения. Выделим на поверхно- |
|
|
|||||
сти фронта единичную |
площадку. |
|
|
||||
Рис. 3.1. К определению понятия |
|||||||
Проведем по обе его стороны два |
|||||||
|
скорости звука |
||||||
сечения. Запишем |
для |
выделен- |
|
|
ного цилиндрического объема закон сохранения массы и закон изменения количества движения (см. разд. 2):
ρ(−v ) + (ρ + dρ)(v + dv ) = 0 , |
(3.17) |
ρv (−v ) + (ρ + dρ)(v + dv )2 = p − (p + dp). |
(3.18) |
Из этих двух уравнений, пренебрегая произведениями малых величин, можно получить v 2 = dpdρ . Очевидно, что для наблюда-
теля, движущегося вместе с жидкостью, жидкость покоится, а возмущение перемещается с такой же скоростью v . Таким образом, v действительно представляет собой скорость распростра-
168
нения звукового фронта или, другими словами, скорость звука, которую далее будем обозначать a . Тогда имеем
a 2 = |
dp |
. |
(3.19) |
|
|||
|
dρ |
|
Соотношение (3.19) определяет скорость звука через возмущения давления dp и плотности dρ на звуковом фронте. Как
видно из (3.19), скорость звука в несжимаемой среде (dρ = 0 )
равна бесконечности. Конечной величина скорости звука получается только в сжимаемой среде. Опыт показывает, что распространение малых возмущений в газе является адиабатическим процессом. Для совершенного газа имеем интеграл адиабаты в
виде p =Cργ (см. (2.71)). Дифференцируя это соотношение и
подставляя результат в (3.19), получим несколько полезных формул для вычисления скорости звука:
|
2 |
|
dp |
|
d (Cρ |
γ |
) |
|
|
|
1 |
|
γ−1 |
|
p |
|
||
a |
= |
= |
|
= γCρ |
γ−1 |
= γC |
γ |
p |
γ |
|
= γ |
= γRT . (3.20) |
||||||
|
dρ |
|
dρ |
|
|
|
|
|
|
ρ |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Последняя формула получена с учетом уравнения состояния Клапейрона.
3.3. Различные формы записи уравнения Бернулли
Получим различные формы записи интеграла Бернулли, которые будут полезны в приложениях.
Несжимаемая жидкость. Пусть ρ = const . Функция P в
этом случае имеет вид P = pρ +const. Если в качестве массовых сил выступает только сила тяжести, то Φ = gz . Тогда интеграл Бернулли (3.11) можно переписать в виде
gz + |
p |
+ |
v 2 |
= const . |
(3.21) |
|
ρ |
2 |
|||||
|
|
|
|
Рассматривая соотношение (3.21) вдоль линии тока и определяя константу в правой части (3.21) по параметрам торможения
(при v = 0 ), получаем уравнение Бернулли: |
|
||||||
gz + |
p |
+ |
v 2 |
= gz0 + |
p0 |
(3.22) |
|
ρ |
2 |
ρ |
|||||
|
|
|
|
||||
|
|
|
169 |
|
|
|
– основное уравнение гидравлики несжимаемой жидкости в поле сил тяжести. Индекс 0 в (3.22) соответствует параметрам заторможенного потока. Иногда это уравнение записывают, поделив
все слагаемые на g : |
z + |
p |
+ |
v2 |
= z0 + |
p |
0 |
. Здесь первое слагае- |
|
ρg |
2g |
ρg |
|||||||
|
|
|
|
|
мое называют геометрическим напором, второе – пьезометрическим, третье – скоростным. Если в качестве несжимаемой среды рассматривать газ (газ можно считать несжимаемой средой при условии, что сжимаемость не влияет на его параметры, например, в случае движения газа с малыми скоростями, более точно – с малыми числами Маха), то, пренебрегая массовыми силами, которые в задачах аэродинамики, как правило, не играют роли, т.е. полагая Φ = 0 , уравнение Бернулли можно записать так:
v 2 |
+ |
p |
= |
p0 |
. |
(3.23) |
|
2 |
ρ |
ρ |
|||||
|
|
|
|
||||
Сжимаемый (совершенный) |
газ. |
Перейдем теперь к рас- |
смотрению сжимаемого совершенного газа. В этом случае функ-
ция P записывается в виде P = γ γ−1 ρp +const . Пренебрегая си-
лами тяжести (полагая Φ = 0 ), получим из (3.11) уравнение Бернулли:
v 2 |
+ |
γ p |
= const . |
|||
|
|
|
|
|||
2 |
γ −1 ρ |
|||||
|
|
Используя соотношения из термодинамики (см. подразд. 2.6) и выражение (3.20) для скорости звука, перепишем это уравнение
в следующих эквивалентных формах: |
|
|
|
|
|
|
|
|
|
|||||||||
v 2 |
+ |
|
γ |
|
p |
= |
|
γ |
|
|
p0 |
|
, |
(3.24) |
||||
2 |
|
γ −1 ρ |
γ −1 ρ0 |
|||||||||||||||
|
|
|
|
|
|
|||||||||||||
v 2 |
+ |
|
a 2 |
= |
|
|
a02 |
|
, |
|
|
|
|
|
(3.25) |
|||
2 |
|
γ −1 |
|
|
γ − |
1 |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
v2 |
+ |
|
a2 |
|
= |
|
|
γ +1 |
aкр2 , |
|
(3.26) |
|||||||
2 |
γ −1 |
2(γ −1) |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
170
v 2 |
+c pT =c pT0 , |
(3.27) |
|||
2 |
|||||
|
υmax2 |
|
|
||
v 2 |
+c pT = |
. |
(3.28) |
||
2 |
2 |
||||
|
|
|
Здесь, как и ранее, индекс 0 соответствует параметрам торможения, "кp" – критическим параметрам, при которых скорость потока равна местной скорости звука, "max" – теоретическому максимальному (предельно возможному) значению скорости течения. Из приведенных уравнений вытекают следующие полезные формулы:
aкр = |
|
2γ |
p0 |
, |
|
(3.29) |
|
|
γ +1 |
ρ0 |
|
||||
|
|
|
|
|
|||
vmax |
= |
2γ |
p0 |
. |
(3.30) |
||
γ − |
1 ρ0 |
||||||
|
|
|
|
Для дальнейшего практически важно и физически целесообразно ввести безразмерные параметры, определяющие не только количественные характеристики, но и качественные свойства течения сжимаемой среды:
• M =va – число Маха (отношение скорости потока к местной скорости звука);
• M* =vaкр – коэффициент скорости или безразмерная
скорость (отношение скорости потока к критической скорости);
• Cr = vvmax – число Крокко (отношение скорости потока
к теоретически максимально возможной скорости).
Эти три безразмерные величины связаны друг с другом. Например, связь числа M* с числом Маха M выражается как
M*2 = |
|
(γ +1)M2 |
. |
(3.31) |
|
+(γ −1)M2 |
|||
2 |
|
|
Наиболее распространенным из перечисленных трех параметров является число Маха. В зависимости от его значения можно провести элементарную классификацию газовых течений: при M <1 имеем дозвуковое течение, при M >1 – сверхзвуковое, при M ≈1 – околозвуковое или трансзвуковое. Каждое из этих течений обладает специфическими физическими и математическими особенностями.
171
Используя приведенные выше уравнения, можно получить соотношения, связывающие текущие параметры газового потока с текущим числом Маха:
a |
|
|
|
|
|
γ −1 |
|
|
|
|
|
− |
1 |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
M 2 |
2 |
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
= 1 + |
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
(3.32) |
|||||||||||
a0 |
|
|
|
2 |
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
T |
|
|
|
|
|
|
γ −1 |
|
|
|
2 |
−1 |
≡ τ(M), |
|
|||||||||||||||||
|
|
|
|
= 1 |
+ |
|
|
|
|
|
|
|
M |
|
|
|
|
|
|
|
(3.33) |
||||||||||
T0 |
2 |
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
ρ |
|
|
|
|
|
γ − |
1 |
|
|
|
|
|
|
− |
|
1 |
|
|
|
|
|
||||||||||
|
|
|
|
|
M 2 |
|
γ−1 |
≡ ε(M), |
|
||||||||||||||||||||||
|
|
|
|
= 1 |
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.34) |
||||||||
ρ0 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
p |
|
|
|
|
|
|
γ −1 |
|
|
|
|
− |
|
γ |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
M 2 |
γ−1 |
≡ π(M), |
|
|||||||||||||||||||||||
|
|
|
|
= 1 + |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.35) |
||||||||||
p |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
γ −1 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
M |
2 |
|
|
|
2 |
|
|
|
||||||||||||||
|
v |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
(3.36) |
||||||||||||
v |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
max |
|
|
+ |
|
γ −1 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
1 |
|
|
|
|
|
|
|
M |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Эти соотношения также можно рассматривать как различные формы записи уравнения Бернулли. Они очень удобны при прак-
тическом решении задач. Для функций τ(M) , ε(M) и π(M)
составлены подробные таблицы. Из приведенных соотношений легко получить значения так называемых критических парамет-
ров (которые соответствуют числу Маха |
M =1): |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
aкр |
|
|
|
2 |
|
|
1 |
|
|
|
|
Tкр |
|
|
2 |
|
|
|
|
ρкр |
|
|
|
|
|
2 |
|
|
1 |
|
|
|
|||||||||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
γ−1 |
|
|
||||||||||||||||||||
|
|
|
|
|
|
= |
|
|
|
|
|
|
|
; |
|
|
|
|
|
|
= |
|
|
|
; |
|
|
|
|
|
= |
|
|
|
|
|
|
; |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
a0 |
|
|
|
γ +1 |
|
|
|
|
|
|
|
T0 |
|
|
γ +1 |
|
|
ρ0 |
|
|
|
|
|
γ +1 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
pкр |
|
|
|
|
|
|
|
|
|
|
|
γ |
|
|
|
|
|
|
|
vкр |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
2 |
|
|
|
γ−1 |
|
|
|
|
|
|
|
|
|
γ −1 |
|
2 |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
= |
|
|
|
|
|
|
|
|
|
; |
|
|
|
|
|
|
|
= |
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
p0 |
|
|
|
|
|
γ +1 |
|
|
|
|
|
|
|
|
|
|
vmax |
|
|
|
γ +1 |
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
Для γ =1,4 |
|
(воздух) эти соотношения дают следующие зна- |
||||||||||||||||||||||||||||||||||||||||
чения: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
aкр |
=0,913 |
; |
|
Tкр |
=0,833; |
ρкр |
=0,634 |
|
|
; |
|
pкр |
=0,528; |
vкр |
=0,408. |
|||||||||||||||||||||||||||||
|
a |
|
|
|
|
|
|
ρ |
|
|
|
|
|
p |
|
|
|
|
v |
|
|
|
|||||||||||||||||||||||
|
0 |
|
|
|
|
|
T |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
max |
|
||||||||||||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
172 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|