Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
ОСНОВЫ АЭРОГАЗОДИНАМИКИ.pdf
Скачиваний:
415
Добавлен:
26.03.2015
Размер:
1.25 Mб
Скачать

∂ϕ

+ 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