Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
9
Добавлен:
10.02.2015
Размер:
397.23 Кб
Скачать

1

2. Основные модели газовой динамики

Газовая динамика — раздел механики жидкости и газа, в котором изучают движение сжимаемой сплошной среды с большой скоростью и ее взаимодействие с твердыми телами. Малая пространственная протяженность области взаимодействия газового потока с твердым телом позволяет отбросить в уравнениях газовой динамики объемные силы, а большая скорость движения требует учета сжимаемости рассматриваемой среды.

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

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

2.1. Дифференциальная форма модели газовой динамики

Состояние движущегося в области V невязкого газа в произвольный момент времени t характеризуют локальные значения его плот-

ности ρ, давления p, массовой плотности внутренней энергии u и

проекций vi (i = 1, 2, 3) вектора v скорости на оси Oxi прямоугольной декартовой системы координат. Дифференциальная форма математической модели газовой динамики включает уравнение неразрывно-

2

сти (закон сохранения массы) в дивергентной форме

∂ρ

+

∂(ρvi)

= 0,

(1)

 

 

∂t

 

∂xi

 

уравнения движения (закон сохранения количества движения сплошной среды) в дивергентной форме

 

∂(ρvi)

+

∂(ρvivj + pδij)

= 0, i, j = 1, 2, 3,

(2)

 

 

 

 

 

∂t

 

∂xj

 

в которых учтены лишь нормальные напряжения σij = −pδij

ij

символ Кронекера) и не учитывается влияние плотности объемных сил (bi ≡ 0), и уравнение закона сохранения энергии в виде

 

∂ρ(viv∂ti

+

 

 

i i

∂xj

j

+ pv

j

= 0;

(3)

 

/2 + u)

 

 

∂ ρ(v v /2 + u)v

 

 

 

 

в нем опущены слагаемые:

qj — проекции вектора плотности те-

плового потока, qV — объемная плотность мощности внутренних источников теплоты и мощности bivi объемных сил.

Условия на поверхности разрыва S , которая может возникнуть при движении газа, согласно (1)–(3), имеют вид

[ρ]Dvnivi

i

i

 

i

 

vnivi

i j

+ pδ

ij

 

j

 

(4)

 

= [ρv

]n

, [ρv

]D = [ρv v

 

]n

,

 

ρ

 

+ u Dn

= ρ

 

+ u vj + pvj

nj ,

 

h

 

i

 

h

 

 

 

 

i

 

 

 

2

 

2

 

 

 

 

где Dn = D · n > 0 — проекция вектора D скорости движения произвольной точки M S на направление вектора n вектора нормали к поверхности S , имеющего проекции ni на координатные оси Oxi; [ ·] — скачок соответствующей величины в этом направлении. Если перейти к системе координат, связанной с этой точкой и движущейся со скоростью D , то вместо (4) получим

[ρv˜i] n = 0, [ρvij + pδij] n = 0, hρ vivi + u v˜j + pvjin = 0, (5)

i j 2 j

где v˜j — проекции вектора v˜ = v − D скорости частиц газа на оси подвижной системы координат.

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

3

энтропии постоянной, положив h = 0. Поэтому u можно отождествить с массовой плотностью A свободной энергии и записать

T

 

 

u = A = Z0

cv dT,

(6)

где cv — удельная массовая теплоемкость при постоянном объеме; T — абсолютная температура газа. Тогда при cv = const третье уравнение (5) примет вид [ρ(vivi/2 + cvT )˜vj + pvj]nj = 0.

Если объединить (3) с (6),

то с учетом равенств vj

∂p

=

∂(pvj)

−p

∂vj

∂xj

 

∂xj

 

∂xj

 

dv

 

d

 

v v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и vi

 

i

=

 

 

i i

 

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

dt

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρvj

dvj

+ ρcv

dT

+ vj

 

∂p

 

+ p

∂vj

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

∂xj

∂xj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dvj

 

∂p

 

 

dT

 

∂vj

 

dT

 

∂vj

 

 

 

 

 

 

 

 

 

= ρ

 

+

 

 

vj

+ ρcv

 

+ p

 

= ρcv

 

+ p

 

= 0,

(7)

 

 

 

 

 

 

dt

∂xj

dt

∂xj

dt

∂xj

так как справедливы уравнения Эйлера при bi ≡ 0. Далее, полагая газ совершенным, уравнение состояния которого имеет вид

p = ρRgT,

(8)

где Rg = Rµg — газовая постоянная; Rµ — универсальная газовая постоянная; µg — молекулярная масса газа, и используя (1) в виде

+ ρ

∂vj

= 0 и (7) имеем

 

 

 

 

 

 

 

 

 

 

 

 

dt

∂xj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

c dT

 

1 dρ

 

= p

cv dln T

 

dln ρ

= 0.

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

RgT

dt

ρ dt

Rg

dt

dt

Отсюда при cv = const получим T/ρκ−1 = const, где κ = cp/cv — показа-

тель адиабаты, а cp = cv + Rg — удельная массовая теплоемкость газа при постоянном давлении, или изобарная теплоемкость,

соответствующая изобарному (при постоянном давлении газа) термодинамическому процессу. Тогда с учетом (8) получаем уравнение

адиабаты Пуассона

p

= const.

(9)

ρκ

 

 

Таким образом, для изоэнтропического процесса движения газа p в силу

(9) зависит лишь от ρ. В этом случае газ называют баротропным. Отметим, что всегда κ > 1. При нормальной температуре для одноатомного газа κ = 5/3, а для двухатомного газа κ = 7/5.

4

Если газ не является совершенным, т. е. его параметры не удовлетворяют (8), то в качестве уравнения состояния можно использовать

уравнение Ван–дер–Ваальса в виде

p =

ρRgT

 

− α1ρ2.

(10)

1

β ρ

 

 

1

 

 

 

Кроме того, для реального газа следует также учитывать зависимость cv от T , которую можно задать в виде

cv = c0 + c1T + c2T 2.

 

 

Тогда вместо (7) получаем

 

 

 

 

(c0 + c1T + c2T 2) dT −

RgT

− α1

dρ = 0.

(11)

ρ − β1ρ2

Если для изоэнтропического процесса представить зависимость T от ρ

в виде T (ρ) = T (0)(ρ) + α1T (1)(ρ) + β1T (2)(ρ) + c1T (3)(ρ) + c2T (4)(ρ) и при-

нять коэффициенты α1, β1, c1, c2 малыми, а в дальнейшем пренебречь слагаемыми, содержащими произведения этих коэффициентов, то после интегрирования (11) найдем

 

T

= ρ

γ1 1 + α γ

 

1 − ρ1−γ1

β

γ

 

ρ

(1

ρ) +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T0

 

 

 

1 1 1 − γ1

 

 

1

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

γ1

 

c1T0

 

 

 

 

 

 

γ1

 

c2T0

 

 

 

 

1

 

 

 

 

 

e

e

 

 

 

 

 

 

 

 

 

+ ρ

 

 

 

 

e(1

ρ

 

 

) +

 

 

(1 − ρ

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c0

 

 

 

2c0

 

и, учитывая (10),

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

e

 

 

p

= ρ1+γ1

1 + α

 

 

1 − ρ1−γ1

 

β (1 + γ

(1

 

 

 

ρ) +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p0

 

 

 

1 1 − γ1

 

 

1

 

 

 

 

1

0

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

e

 

 

 

e

 

1+γ1 c1T0

 

 

 

 

 

 

γ1

 

 

e

c2T0

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

+ ρ

 

 

 

 

 

 

 

(1

ρ

 

 

)

 

 

 

(1

ρ

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c0

 

 

 

2c0

 

 

где α1

=

α1ρ0

; γ1 =

 

Rg

; ρ =

e

; p0, ρ0 и T0 e

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

e

 

 

 

 

e

ρ

 

 

 

 

 

 

 

 

 

 

 

 

параметры газа в некотором

 

 

RgT0

 

 

 

c0

ρ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

В случае малости скорости v и ее производных, считая малыми и изменения p = p − p0 и ρ = ρ − ρ0 давления p и плотности ρ по сравнению с давлением p0 и плотностью ρ0 газа в его невозмущенном состоянии, из (1) и (2) с учетом (8) и (9) получаем

 

∂ρ

1

 

∂p

 

 

 

 

 

∂vi

 

 

∂vi

 

 

1 ∂p

 

a02

 

∂ρ

(13)

 

 

 

 

 

 

 

 

≈ −ρ0

 

,

 

 

≈ −

 

 

 

≈ −

 

 

 

,

 

∂t

a02

 

∂t

∂xi

 

∂t

ρ0

∂xi

ρ0

∂xi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dp

ρ=ρ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

κp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

κRgT0 — скорость

звука

(скорость

где a0 = r

 

= q ρ00 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

возмущений

 

p и

ρ). Заменив в (13)

распространения в газе малыхp

 

 

 

 

 

 

∂t2
= a2 2(r ρ) 0 ∂r2
2(r ρ)

5

приближенные равенства точными, приходим к известным волновым уравнениям

2ρ

= a2

2ρ

 

2p

= a2

2p

(14)

 

 

,

 

 

,

∂t2

∂xi∂xi

∂t2

 

0

 

0

∂xi∂xi

 

описывающим распространение малых возмущений. Пространственную форму распространяющегося возмущения называют волной.

Если до возникновения возмущения плотности газа векторное поле скоростей обладало потенциалом Φ0(x,t), где x — радиус-вектор точ-

ки с пространственными координатами xi, то и при распространении малого возмущения это поле остается потенциальным. Действительно, из второго соотношения (13), полагая, что возмущение возникло в момент времени t = 0 и vi|t=0 = ∂Φ0(x,0)/∂xi, интегрированием по времени t получаем

 

 

a2

t

∂Φ(x,t)

 

a2

t

vi = vi

t=0

Z ρdt =

, Φ(x,t) = Φ0(x,0) −

Z ρdt —

ρ0 ∂xi

∂xi

ρ0

 

 

0

0

 

 

0

0

 

 

 

 

 

 

 

 

потенциал поля скоростей газа в процессе распространения малого воз-

мущения плотности газа, удовлетворяющий уравнению 2Φ = a2 rx Φ.

∂t2 0

При установившемся движении невозмущенного газа потенциал Φ0 не зависит от времени t. В этом случае

∂Φ

= −a02

ρ

(15)

 

 

.

∂t

ρ0

Пусть в момент времени t = 0 в неподвижном газе плотностью ρ0

впределах шаровой области радиусом r0 возникает малое возмущение ρ плотности. Это возмущение в силу центральной симметрии поро-

ждает сферическую волну ρ(r,t) = ρ0 + ρ(r,t), распространение ко-

торой описывается первым уравнением (14) в виде

(r — расстояние от центра этой области). При этом потенциал Φ(r,t) поля скоростей будет удовлетворять аналогичному уравнению

 

 

 

2(rΦ)

= a2

2(rΦ)

 

 

(16)

 

 

 

∂t2

 

∂r2

 

 

 

 

 

 

0

 

 

 

 

с начальными условиями Φ(r,0) = 0 и

∂Φ(r,0)

=

−a02f(r), где с учетом

 

 

∂t

(15)

 

 

 

 

 

 

 

 

 

 

 

f(r) = 1 +

ρ

при r < r0,

 

f(r) = 1

при r > r0.

(17)

 

 

 

ρ0

 

 

 

 

 

 

 

 

6

Общее решение (16) имеет вид Φ(r,t) =

1

 

fΦ(a0t − r) + gΦ(a0t + r) ,

r

где fΦ и gΦ — произвольные дважды

дифференцируемые функции

 

 

 

.

Поскольку при r = 0 потенциал остается конечным, вместо последнего равенства следует записать

 

Φ(r,t) =

fΦ(a0t − r) − fΦ(a0t + r)

.

 

 

 

(18)

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

Тогда с учетом (15) получим

 

 

 

 

 

 

 

 

 

 

 

 

ρ ∂Φ(r,t)

 

ρ0

f0

(a0t + r) f0 (a0t r)

 

 

 

a02

 

∂t

 

 

 

a0r

 

ρ(r,t) =

0

 

 

 

=

 

Φ

 

Φ

 

.

(19)

 

 

 

 

 

 

 

 

 

 

 

Используя начальные условия, (18) и (19) приведем к виду fΦ(−r) −

− fΦ(r) = 0 и fΦ0 (−r) − fΦ0 (r) = −a0rf(r). Отсюда находим

 

fΦ0 (−r) = −fΦ0 (r) = −

a0rf(r)

(20)

 

.

2

При t > 0 для точек за пределами возмущенной шаровой области a0t + + r > r0, т. е. аргумент функции fΦ0 (a0t + r) всегда больше r0. Поэтому для этих точек в силу (17) и (20) и fΦ0 (a0t + r) = a0(a0t + r)/2.

 

Для вычисления

 

Φ

 

r)

 

 

выделим три промежутка времени

 

f

0 (a0t

 

 

 

 

распространения возмущения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,

 

r − r0

 

;

 

r − r0

,

r + r0

;

 

r + r0

,

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

 

 

a0

 

 

 

a0

 

 

 

 

a0

 

 

 

 

В

первом промежутке

r − a0t > r0, и

 

из

(17)

и

(20)

следует, что

f0

(a0t

 

r) = f0

 

(r

a0t) =

a0(r −a0t)

; во втором — a0t

 

r < r0, что в

силу (17) и (20) дает fΦ0

(a0t

 

 

r) = a0(a0t −r)

1 +

 

ρ

; в третьем — a0t

 

Φ

 

 

− −

 

 

 

 

 

 

2

 

 

 

a0(a0t

r)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

ρ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

r > r0

и, согласно (17)

и (20), f0

(a0t

 

r) =

 

 

 

 

 

. Таким образом,

из (19) следует, что в первом и третьем промежутках ρ(r,t) = ρ0, т. е. газ остается невозмущенным, а во втором промежутке ρ(r,t) = ρ0 +

+

ρ(r −a0t)

.

 

2r

Изменение плотности газа в рассмотренных промежутках времени представлено на рис. 0.1, причем на границах этих промежутков плотность газа изменяется скачком. При фиксированном значении r изменение плотности происходит в течение промежутка времени 2r0/a0. В момент времени t > 0 передний фронт сферической волны (позиция 1 на рис. 0.2) достигает точек, расположенных на расстоянии r = a0t + r0 от центра O шаровой области начального возмущения плотности газа. При t > 2r0/a0 за пределами этой области формируется задний

7

фронт 2 сферической волны, который достигает точек, расположенных на расстоянии r = a0t − r0 от центра O, и эта волна располагается в сферическом слое толщиной 2r0, равномерно расширяющемся со скоростью a0. Сфера 3 радиусом r = a0t делит этот слой на внешнюю (ρ > ρ0) и внутреннюю (ρ < ρ0) части.

Рис. 0.1

Рис. 0.2

а

б

в

Рис. 0.3

В случае точечного источника возмущения r0 = 0 и шаровой слой переходит в сферическую поверхность, также равномерно расширяющуюся со скоростью a0. Если такой источник движется из начальной точки O прямолинейно с постоянной скоростью v0, то для неподвижного наблюдателя сферическая волна в различные моменты времени не будет представлять собой концентрические сферы. При v0 < a0 каждому моменту времени t1, 2t1, 3t1, . . . соответствует сфера, охватывающая все сферы для предшествующих моментов времени (рис. 0.3, а). В случае v0 = a0 все сферы имеют общую точку касания, соответствующую начальному положению точечного источника (рис. 0.3, б). При v0 > a0 все сферы касаются прямого кругового конуса (рис. 0.3, в), называемого конусом Маха. Вершина этого конуса совпадает с начальным положением источника, а полуугол arcsin(a0/v0) при вершине называют

углом Маха.

8

2.2. Одномерное течение невязкого газа

При установившемся одномерном течении невязкого газа в прямолинейном канале постоянного поперечного сечения из (1) следует ρv = const, где ρ — плотность газа, а v — модуль вектора скорости. Отсюда ρdv + v dρ = 0, или

+

dv

= 0.

(21)

ρ

 

v

 

 

 

Для такого течения из (2) получим ρv dv + dp = 0, где p — давление газа. Используя последнее равенство и (21), находим

dp

= v2.

(22)

 

 

При изоэнтропическом процессе движения газа dp/dρ, согласно (9),

является функцией только p или ρ. В этом случае (22) определяет скорость распространения в канале малого возмущения неизменной фор-

мы, удовлетворяющую как уравнению неразрывности, так и уравнению движения, называемую местной скоростью звука и обозначаемую через a. Для совершенного газа с учетом его уравнения состояния (8)

и уравнения адиабаты Пуассона (9) запишем

 

a = s

 

 

= r

 

= pκRgT ,

(23)

 

ρ

 

 

dp

κp

 

 

 

 

где Rg — газовая постоянная; T — абсолютная температура газа.

Из (3) с учетом (6) получим

 

 

 

 

 

v2

p

 

 

 

 

+ cvT +

 

 

= const,

(24)

 

2

ρ

где cv = const — удельная массовая теплоемкость газа при постоянном объеме, или, учитывая (8) и (23),

 

v2

+

cvp

+

p

=

v2

+

 

 

 

κp

 

=

v2

+

a2

 

= const =

 

a02

 

,

(25)

2

 

 

 

 

− 1)ρ

 

κ −

 

κ −

1

Rgρ ρ

2

 

2

 

1

 

 

 

 

 

 

 

где κ = 1 + Rg/cv — показатель адиабаты; a0 =

 

 

=

 

 

 

 

κp00

 

 

κRgT0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

заторможенном газе при давлении тор

 

скорость звука в полностью

2

 

 

 

 

 

 

 

 

p

 

p

 

 

-

можения p0 =

ρ0

p +

κ −κ

1

ρ0

v

, которому соответствуют темпера-

ρ

2

тура торможения T0 = T + v2/(2cp) (cp — удельная массовая тепло-

емкость при постоянном давлении) и плотность ρ0 = p0/(RgT0).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

Из (25) с учетом (9) имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a02

 

 

 

 

 

1

 

 

 

 

 

 

κ − 1

 

 

 

 

κ

 

 

 

 

 

 

 

 

 

 

 

p

 

1− κ

 

 

 

 

 

 

 

v

2

 

 

 

 

 

 

v2 =

 

 

 

1

 

 

,

p = p

 

 

1

 

 

κ−1 .

(26)

 

κ − 1

p0

 

 

 

 

a0

 

 

 

 

 

 

 

0

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= a0q

 

 

 

 

Скорость

v

достигает максимального значения

vmax

 

2

 

при

 

κ −1

p = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если в (25) положить v = a, т. е. приравнять скорость потока местной скорости звука, то получим значение критической скорости

звука:

 

 

 

 

 

 

 

 

 

 

 

 

a = a0r

 

 

 

 

.

 

 

 

(27)

 

κ + 1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

κ

 

 

 

 

дозвуковым, если 0 < v < a , и

 

 

κ−1

 

,

 

a < v < vmax.

Ей соответствует давление p = p0

 

 

 

. Течение газа называют

 

 

 

 

 

κ + 1

 

 

 

 

 

сверхзвуковым

 

если

 

Безразмерный параметр M = v/a, характеризующий влияние сжима-

емости газа,

называют числом Маха. Используя это число, вместо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

κ

 

ρ

 

 

 

κ

 

1

 

κ1 1

 

 

 

 

 

 

p

2

 

 

 

 

 

второго равенства

(26)

получаем

p0

=

 

1 +

κ −1

M2

κ−1

и с учетом (9)

 

 

 

ρ

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

=

1 +

 

 

M2

 

 

, а в соответствии с (8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T0

= 1 +

κ −

1

M2,

 

 

 

 

(28)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

2

 

 

 

 

 

 

 

что позволяет, согласно (23) и (27), записать

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T (κ + 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

κ + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

= a

 

 

 

2T0

= a

 

 

2 + (κ − 1)M2

.

(29)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда следует, что a < a при M < 1 и a > a при M > 1.

При изоэнтропическом процессе неустановившегося одномерного течения невязкого газа в прямолинейном канале постоянного поперечного сечения из (1), (2) и (9) имеем

∂ρ ∂(ρv)

 

∂v

 

∂v

 

1

 

∂p

 

 

 

+

 

 

= 0,

 

+ v

 

= −

 

 

 

, p = p(ρ),

(30)

∂t

∂x

∂t

∂x

ρ

∂x

где t — время; x — координата, отсчитываемая вдоль оси канала. Для малых отклонений p = p − p0 и ρ = ρ − ρ0 давления p и плотности ρ от соответственно давления p0 и плотности ρ0 газа в невозмущенном состоянии при малой скорости v и ее производных из (13) следуют одномерные волновые уравнения

2ρ

= a2

2ρ

,

2v

= a2

2v

,

(31)

∂t2

 

∂x2

∂t2

∂x2

0

 

0

 

 

10

общие решения которых имеют вид ρ(x,t) = fρ(x − a0t) + gρ(x + a0t) и v(x,t) = fv(x − a0t) + gv(x + a0t). Анализ этих решений показывает,

что если, например,

 

∂ρ(x,t)

 

 

 

= 0 в начальный момент времени t = 0

 

∂t

 

 

t=0

и ρ(x,0) = F

 

(x) = f

 

 

 

 

 

кривая

ρ

ρ

(x) + g

 

(x), то по прошествии времени t

1

 

 

 

 

 

ρ

 

 

 

распределения плотности газа

сдвигается без искажения вдоль оси Ox

вправо и влево на расстояние a0t1. Вдоль линий x −a0t = const решения fρ(x − a0t) = const и fv(x − a0t) = const, а вдоль линий x + a0t = const

решения gρ(x + a0t) = const и gv(x + a0t) = const. Эти линии называют

характеристиками волнового уравнения.

Если возмущения скорости не являются малыми, то необходимо решать систему нелинейных уравнений (30), и, как правило, численно. В частном случае зависимости ρ только от v, т. е. ρ = f(v), удовлетворяющей первым двум уравнениям (30), с учетом третьего уравнения (30)

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f0 ∂v

 

∂v

 

 

 

 

∂v

 

 

 

 

f ∂v

 

 

 

∂v

 

 

 

∂v

 

 

 

 

 

 

 

+ v

 

 

= −

 

 

 

 

,

 

 

 

+ v

 

 

= −a2

 

 

 

.

 

∂t

∂x

f0

∂x

∂t

∂x

f

∂x

 

 

 

 

 

 

f

 

 

 

f0

 

 

 

 

dv

 

 

 

 

 

 

 

 

 

Отсюда следует

 

 

 

 

= a2

 

, или

 

 

 

= ±

 

 

, что позволяет первые два

 

f0

f

 

ρ

a

уравнения (30) представить в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ρ

 

 

 

 

 

∂ρ

 

 

 

 

∂v

 

 

 

 

 

 

∂v

 

 

 

+ (v

± a)

 

= 0,

 

 

+ (v ± a)

 

= 0.

 

 

∂t

∂x

 

∂t

∂x

Непосредственной проверкой можно установить, что этим уравнениям удовлетворяют соотношения

ρ(x,t) = Fρ x − (v ± a)t , v(x,t) = Fv x − (v ± a)t ,

(32)

где Fρ и Fv — произвольные функции.

При |v/a| → 0 соотношения (32) переходят в общие решения волновых уравнений (31). Из (32) следует, что в подвижной системе координат, начало которой имеет скорость v ±a, плотность и скорость частиц газа остаются постоянными. Отсюда можно сделать вывод, что возмущения будут распространяться со скоростью v ± a, в общем случае различной в разных точках потока газа, поэтому пространственная форма возмущений будет изменяться во времени.

 

 

Примем в (32) v(x,t) = Fv x − (v + a)t

и в момент времени t =

= 0 волнообразное распределение возмущения скорости v вдоль оси Ox (рис. 0.4, а). Поскольку возмущения в области, где v > 0, распространяются быстрее, чем в области, где v < 0, то при t > 0 гребень волны будет стремиться нагнать впадину, что приведет к более резкому изменению скорости в зоне смены ее знака (рис. 0.4, б). В рамках

Соседние файлы в папке МММСС (От Кувыркина Г.Н.)