Линейные и нелинейные волны - Уизем Дж.
Скачать (прямая ссылка):
содержит двумерные волны, распространяющиеся с локальной скоростью звука
относительно локального течения, как показано на рис. 8.5.
Рис. 8.5. Цилиндрические волны, возникающие при дифракции ударной волны.
Наши приближенные уравнения некоторым образом описывают след этих
цилиндрических волн при их пересечении с ударной волной.
Волновое распространение возмущений по ударной волне проявляется при
изучении уравнений в характеристической форме.
Гл. 8. Динамика ударных волн
276
При подстановке А = А (М) в уравнения (8.59) и (8.60) они принимают вид
дв А' (М) дМ
(?Р М да
50 I 1 дМ
= 0, = 0.
(8.65)
да 1 А (М)
Характеристическая форма такова:
(ж±сж)(в±Ш=0- <8-6в)
где с - функция от М, определяемая как
с('?) = / 1F' (8'67)
Поскольку А' (М) <С 0, характеристики вещественны, и мы имеем нелинейные
волны, распространяющиеся со скоростями
*=±с
аа
в (а, (Д-пространстве. Эти волны переносят изменения формы и
интенсивности ударной волны по самой ударной волне. Инварианты Римана
даются уравнением (8.66), а именно
0 + j - const на ¦J~=c> (8.68)
0-j const на -с• (8.69)
Эти уравнения во всех смыслах аналогичны исходным уравнениям одномерной
нелинейной- газовой динамики, изученным в гл. 6, и развитые там идеи и
методы можно применить к волнам, распространяющимся вдоль ударной волны.
Зависимость А (М) выводится из уравнения (8.25), которое можно переписать
в виде
(8-70)
Отсюда
Ас
-№Г. "*•">
m-i 11/2
К(М)
и интеграл в формулах для инвариантов Римана равен м м
(8.72)
1 1
Могут оказаться полезными явные формулы для слабых ударных волн с ЛГ - 1
<1 и для сильных ударных волн с М 1.
8.5. Распространение волн по ударной волне
277
Эти формулы имеют вид
4, А. ~ 1*0=1)1
М - 1 \ 1/2
Ас ¦
при М-у 1 (8.73)
X ~ п= 5,0743 для у = 1,4, ('аГ')"'
Ac~n~l,2M, to (71/) ~ n1/2 In 71/
при 71/ -оо. (8.74)
Характеристические соотношения проще всего получаются в (а, Р)-
координатах, но в приложениях к конкретным краевым задачам иногда
предпочтительнее описание в декартовых координатах (ж, у). Выражения
(8.68) и (8.69) простыми выкладками преобразуются к нужной
форме. Заметим, что
dy __ Va + yydV/da СЯ 751
dx ха zp dpjda '
где xa = 71/ cosB, ya - 71/ sin0, x$ = -A sin0, у p - .4 cos 0.
Следовательно, характеристики dftlda = ± с записываются как
^=tg(0±m), (8.76)
где
а характеристические уравнения принимают вид
0 ± to (71/) = const на С±: ^- = tg(0±m). (8.78)
Значения т (71/) и со (71/) приведены в табл. 8.3, заимствованной из
статьи Брисона и Гросса [1].
Последние уравнения можно вывести и непосредственно иэ двумерной формы
уравнений (8.50) и (8.51). При ах = cos0/71/* ау = sin 0/71/
эквивалентная система уравнений для 0 и 71/ имеет вид
(8.79)
JL | / sin в \ д ( cos 0 \
дх \~~M~ )~1)у\ )
9_ 1 COS0 \ , д [ A J + dlT1 / sinG \
дх
Таблица 8.3
Число Маха М Площадь трубки лучей А, Ах10-/7 N
Характеристический угол т, град. Интеграл ш
1 ОО 0 0
1,000001 3,668749 +10 - 0,003
1,00001 3,668672 +8 - 0,009
1,0001 3,667902 +6 0,403 0,028
1,001 3,660213 +4 1,280 0,089
1,01 3,584696 +2 4,002 0,283
1,05 1,310728 +1 8,544 0,633
1,10 2,946288 +0 11,474 0,896
1,15 1,184152 +0 13,142 1,097
1,20 6,053638 -1 14,843 1,266
1,25 3,536658 -1 15,958 1,414
1,30 2,250720 -1 16,859 1,547
1,35 1,520662 -1 17,604 1,669
1,40 1,074028 -1 18,231 1,728
1,45 7,850741 -2 18,766 1,887
1,50 5,898186 -2 19,228 1,984
1,55 4,531934 -2 19,630 2,077
1,60 3,548150 -2 19,983 2,165
1,65 2,822580 -2 20,295 2,249
1,70 2,276434 -2 20,572 2,330
1,75 1,858064 -2 20,820 2,406
1,80 1,532637 -2 21,042 2,480
1,85 1,276079 -2 21,242 2,551
1,90 1,071389 -2 21,423 2,619
1,95 9,063299 -3 21,587 2,685
2,00 7,719471 -3 21,736 2,749
2,05 6,615861 -3 21,872 2,811
2,10 5,702352 -3 21,997 2,871
2,15 4,940726 -3 22,111 2,929
2,20 4,301517 -3 22,216 2,985
2,25 3,761766 -3 22,312 3,040
2,30 3,303423 -3 22,401 3,094
2,40 2,576553 -3 22,560 3,203
2,50 2,037086 -3 22,696 3,302
2,60 1,630023 -3 22,814 3,388
2,70 1,318343 -3 22,916 3,477
2,80 1,076566 -3 23,006 3,563
2,90 8,868121 -4 23,085 3,645
3,00 7,363072 -4 23,154 3,724
3,20 5,184216 -4 23,271 3,875
8.6. Вторичные ударные волны
279
Продолжение табл. 8.3
Число Маха М Площадь трубки лучей А, Ах10-^ N
Характеристический угол т, град. Интеграл о
3,40 3,740925 -4 23,364 4,015
3,60 2,757067 -4 23,439 4,148
3,80 2,069662 -4 23,501 4,272
4,00 1,578970 -4 23,552 4,398
4,50 8,519558 -5 23,647 4,660
5,00 4,926060 -5 23,710 4,900
6,00 1,921342 -5 23,788 5,314
7,00 8,705958 -6 23,832 5,672
8,00 4,395269 -6 23,859 5,966
9,00 2,408270 -6 23,876 6,232
10,00 1,407051 -6 23,889 6,470
15,00 1,786391 -7 23,917 7,385
20,00 4,141420 -8 23,926 8,033
100,00 1,172427 -11 23,937 11,67
оо 0 23,938 оо
Несколько более длинные непосредственные вычисления показывают, что
характеристические уравнения совпадают с приведенными выше.
8.6. Вторичные ударные волны
Функция с (М) возрастает с ростом М. Следовательно, волны, движущиеся в
положительном направлении и несущие возрастание числа Маха М и 0, должны