КАТЕГОРИИ: Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748) |
Ol — V2X 55 страница
ТТ—pT-J+T' -77-7+Г7Г- (89'10> Отношение T2/Ti неограниченно растет вместе с р2/р\, т. е. ска- 0^ = V'1Tlp2Fl• ^V Они растут пропорционально корню из давления р2.
Наконец, приведем соотношения для ударных волн слабой интенсивности, представляющие собой первые члены разложений по степеням малого отношения г = (р2— Р\)/р\- v±lz> £i.= 1+rzl
Pi Y 2Y2 \ • / Здесь сохранены члены, дающие первую поправку к значениям, отвечающим звуковому приближению.
Задачи 1. Получить формулу
где с» — критическая скорость (L. Prandtl). Решение. Поскольку величина ш + и2/2 непрерывна на ударной волне, можно ввести критическую скорость, одинаковую для газов 1 \\ 2 согласно- YPi,», УР2 v\ у + 1 -------- + -7Г = ТТ.--------- ГГГ- + -ST" = -------------------- ГГ- с (Y-1) р, 2 <Y_l)p2 ' 2 2 (y-1) (cp. (83,7)). Определяя из этих равенств р2/р2 и р4/р, и подставляя их в-уравнение о, — v2 = — Si— 1 2 Р2»2 Pl»l (результат комбинирования (85,1) и (85,2)), получим Y+1 / с2 \ —я— (»1 - о2) I 1 --) = 0. 2у \ v,v7 / - Ввиду того что t>i Ф с2, отсюда следует искомое соотношение. 2. Определить отношение p2/pi по заданным температурам Ti, Тг для ударной волны в термодинамически идеальном газе с непостоянной теплоемкостью. Решение. Для такого газа можно лишь утверждать, что w (как и е) есть функция только от температуры и что р, V, Т связаны уравнением состояния pV = RT/\i. Решая уравнение (85,9) относительно р%1р\, получаем:
где Wi -- w(Ti), w2 — w(T2).
§ 90. Гофрировочная неустойчивость ударных волн Соблюдение условий эволюционности само по себе необходимо, но еще недостаточно для гарантирования устойчивости ударной волны. Волна может оказаться неустойчивой по отношению к возмущениям, характеризующимся периодичностью вдоль поверхности разрыва и представляющим собой как бы «рябь», или «гофрировку», на этой поверхности (такого рода возмущения рассматривались уже в § 29 для тангенциальных разрывов)1). Покажем, каким образом исследуется этот вопрос для ударных волн в произвольной среде (С. П. Дьяков, 1954). Пусть ударная волна покоится, занимая плоскость х = 0; жидкость движется сквозь нее слева направо, в положительном направлении оси х. Пусть поверхность разрыва испытывает возмущение, при котором ее точки смещаются вдоль оси х на малую величину C = C0e'(V-»')f (90,1)
где ky — волновой вектор «ряби». Эта рябь на поверхности вызывает возмущение течения позади ударной волны, в области х > 0 (течение же перед разрывом, х < 0, не испытывает возмущения в силу своей сверхзвуковой скорости). Произвольное возмущение течения складывается из энтропийно-вихревой волны и звуковой волны (см. задачу к § 82). В обоих зависимость величин от времени и координат дается множителем вида ехр [г (кг — ы()] с той же частотой со, что и в (90,1). Из соображений симметрии очевидно, что волновой вектор к лежит в плоскости ху; его «/-компонента совпадает с ky в (90,1), а х-компонента различна для возмущений двух типов.
В энтропийно-вихревой волне kv2 = со, т. е. kx = со/у2 («2 — невозмущенная скорость газа за разрывом). В этой волне возмущение давления отсутствует, возмущение удельного объема связано с возмущением энтропии, 6У(ЭНТ) = (dV/ds)p8s, а возму- щение скорости подчинено условию k 6v<9HT»=— oo«»"» + k,dv\Jm)=0. (90,2) В звуковой волне в движущемся газе связь между частотой и волновым вектором дается равенством (ш — kv)2 —с2й2 (см. (68,1)); поэтому kx в этой волне определяется уравнением («-W=cK*x+ **,)• (90,3) Возмущения давления, удельного объема и скорости связаны соотношениями: бр(зв)==_(С2/1/2)20т/(зВ)( (90>4) (ю — v2kx) 6v(3B> - V2k6p<3*\ (90,5) Возмущение в целом представляется линейной комбинацией возмущений обоих типов: So = 6о<энт> + 6о(зв>, 6V = 6У(ЭНТ> + 6V'3B>, 6р = 6р,зв>. (90,6) Оно должно удовлетворять определенным граничным условиям на возмущенной поверхности разрыва. Прежде всего, на этой поверхности должна быть непрерывна тангенциальная к ней составляющая скорости, а скачок нормальной составляющей должен выражаться через возмущенные давление и плотность равенством (85,7). Эти условия записываются как v1t = (v2 + 6v)t, v,n — (v2 + 6v) n = = [(p2-pi + bp)(vi-v2-6v)}42, где t и n — единичные векторы касательной to, = to,~^[^-^]. (90,7) Далее, возмущенные значения р2 + Ьр и У2 + 6У2 должны удовлетворять тому же уравнению адиабаты Гюгонио, что и невозмущенные р2 и V2. Отсюда получаем условие, связывающее Ър и 6V: 6Р = ЖЬУ> <90>8> где производная берется вдоль адиабаты. Наконец, еще одно соотношение возникает из связи между потоком вещества через поверхность разрыва и скачками давления и плотности на ней. Для невозмущенного разрыва это соотношение выражается формулой (85,6), а для возмущенного аналогичное соотношение есть 1 / \2 Рг — Р\ + бр —я- (V^ — Un)2 = — V2 Vl-V2-6V где u — скорость точек поверхности разрыва. В первом приближении по малым величинам имеем un = —г'ю£; разлагая написанное равенство также и по степеням бр и 6V, получим: 2tog== бр ЬУ
Равенства (90,2), (90,4—5), (90,7—9) составляют систему восьми линейных алгебраических уравнений для восьми величин £, бр, 6У<ЭНТ>, 6К(ЗВ), бо^»,. 6v^y '). Условие совместности этих уравнений (выражаемое равенством нулю определителя их коэффициентов) имеет вид: 2-f +7) " (£+*«)(ffl - (| + k) ~°' №,0) где для краткости обозначено h = ^(dVi/dpi), а / имеет обычный смысл: / = Vi/Vi = V2/V2. Величину кх в (90,10) надо понимать как функцию ky и со, определяемую равенством (90,3). Условие неустойчивости состоит в существовании возмущений, экспоненциально возрастающих со временем, причем они должны экспоненциально убывать с удалением от поверхности разрыва (т. е. при х-*со); последнее условие означает, что источником возмущения является сама ударная волна, а не какой-то внешний по отношению к ней источник. Другими словами, волна неустойчива, если уравнение (90,10) имеет решения, у которых 1ти>0, Im kx>0. (90,11)
Исследование уравнения (90,10) на предмет выяснения условий существования таких решений довольно громоздко. Мы не будем производить его здесь, ограничившись указанием окончательного результата2). Гофрировочная неустойчивость ударной
волны возникает если ;2 dp2 или
напомним, что производная берется вдоль ударной адиабаты (при заданных ри Vi) 1). Условия (90,12—13) отвечают наличию у уравнения (90,10) комплексных корней, удовлетворяющих требованиям (90,11).Но в определенных условиях это уравнение может иметь также и корни с вещественными со и kx, отвечающие «уходящим» от разрыва реальным незатухающим звуковым и энтропийным волнам, т. е. спонтанному излучению звука поверхностью разрыва. Мы будем говорить о такой ситуации как об особом виде неустойчивости ударной волны, хотя неустойчивости в буквальном смысле здесь нет, — раз созданное на поверхности разрыва возмущение (рябь) неограниченно долго продолжает излучать волны, не затухая и не усиливаясь при этом; энергия, уносимая излучаемыми волнами, черпается из всей движущейся среды2). Для определения условий возникновения этого явления, преобразуем уравнение (90,10), введя угол 0 между к и осью Jtj тогда cJkx = со0 cos 0, c0k„ = щ sin 0, <о = со0 (1 + — cos 9), «о = ct (kl + ky) («о — частота звука в системе координат, движущейся вместе с газом за ударной волной), и получаем квадратное относительно cosO уравнение: -5 г (1+^)=0. (90,13)
Скорость распространения звуковой волны в движущемся со скоростью v2 газе, по отношению к неподвижной поверхности разрыва, есть t/2-f-c2cos 0. Звуковая волна будет уходящей, если эта сумма положительна, т. е. если — v2/c2 < cos 8 < 1 (90,16) (значения cos0<0 отвечают случаям, когда вектор к направлен в сторону разрыва, но снос звуковой волны движущимся газом делает ее все же «уходящей»). Спонтанное излучение звука ударной волной возникает, если уравнение (90,15) имеет корень, лежащий в этих пределах. Простое исследование приводит к следующим неравенствам, определяющим область этой неустойчивости '): 1 — v\lc\ — v,vJc\ „ dVo о»
1 — v2/c2 + vKv2jc2 dp2 c2 (нижний и верхний пределы здесь фактически отвечают нижнему и верхнему пределам в условиях (90,16)). Область (90,17) примыкает к области неустойчивости (90,13), расширяя ее.
К происхождению неустойчивости ударных волн в области (90,17) можно подойти также и с несколько иной точки зрения, рассмотрев отражение от поверхности разрыва звука, падающего на нее со стороны сжатого газа. Поскольку ударная волна движется относительно газа впереди нее со сверхзвуковой скоростью, то в этот газ звук не проникает. В газе же позади волны будем иметь, наряду с падающей звуковой волной, еще и отраженную звуковую и энтропийно-вихревую волны (а на самой поверхности разрыва возникает рябь). Задача об определении коэффициента отражения по своей постановке близка к задаче об исследовании устойчивости. Разница состоит в том, что наряду с подлежащими определению амплитудами исходящих от разрыва (отраженных) волн в граничных условиях фигурирует еще и заданная амплитуда приходящей (падающей) звуковой волны. Вместо системы однородных алгебраических уравнений мы будем иметь теперь систему неоднородных уравнений, в которых роль неоднородности играют члены с амплитудой падающей волны. Решение этой системы дается выражениями, в знаменателях которых стоит определитель однородных уравнений,— как раз тот, приравнивание которого нулю дает дисперсионное уравнение спонтанных возмущений (90,10). Тот факт, что в области (90,17) это уравнение имеет вещественные корни для cos 0, означает, что существуют определенные значения угла отражения (и тем самым угла падения), при которых коэффициент отражения становится бесконечным. Это — другая формулировка возможности спонтанного излучения звука, т. е. излучения без падающей извне звуковой волны. То же самое относится и к коэффициенту прохождения звука, падающего на поверхность разрыва спереди, навстречу ей. В этом случае не существует отраженной волны, а позади поверхности разрыва возникают прошедшие звуковая и энтропийно-вихревая волны. В области (90,17) возможно обращение коэффициента прохождения в бесконечность1). Скажем несколько слов о некоторых возможных, в принципе, типах ударных адиабат, содержащих области рассмотренных неустойчивостей2). Условие (90,12) требует отрицательной производной dptjdV^, причем ударная адиабата должна быть наклонена (к оси абсцисс) в точке 2 менее круто, чем проведенная в нее хорда 12 (т. е. обратно тому, что имеет место в обычных случаях — рис. 53). Для этого адиабата должна перегнуться, как показано на рис. 60; условие неустойчивости (90,12) выполняется на участке ab. \Р V 1 Рис. 60 Рис. 61
Условие (90,13) требует положительности производной dpi/dVi, причем наклон адиабаты должен быть достаточно мал. На рис. 60 это условие выполняется на определенных отрезках адиабаты, непосредственно примыкающих к точкам а и Л и расширяющих, таким образом, область неустойчивости. Условие (90,13) может оказаться выполненным и на участке (cd на рис. 61) адиабаты, не содержащей участка типа аЪ. Условие (90,17) еще менее жестко, чем (90,13) и еще дополнительно расширяет область неустойчивости на адиабатах Гю-гонио с dp2/dV2 > 0. Более того, нижний предел в (90,17) может быть отрицательным, так что неустойчивость этого типа может, в принципе, иметь место и в некоторых участках адиабат обычного вида, со всюду отрицательной производной dp2/dV2. Вопрос о судьбе гофрировочно-неустойчивых ударных волн тесно связан со,следующим замечательным обстоятельством: при выполнении условий (90,12) или (90,13) решение гидродинамических уравнений оказывается неоднозначным (С. S. Gardner, 1963). Для двух состояний среды, 1 я 2, связанных друг с другом соотношениями (85,1—3), ударная волна является обычно единственным решением задачи (одномерной) о течении, переводящем среду из состояния / в 2. Оказывается, что если в состоянии 2 выполнены условия (90,12) или (90,13), то решение указанной гидродинамической задачи не однозначно: переход из состояния 1 в 2 может быть осуществлен не только в ударной волне, но и через более сложную систему волн. Это второе решение (его можно назвать распадным) состоит из ударной волны меньшей интенсивности, следующего за ней контактного разрыва и из изэнтропической нестационарной волны разрежения (см. ниже § 99), распространяющейся (относительно газа позади ударной волны) в противоположном направлении; в ударной волне энтропия увеличивается от s\ до некоторого значения s3 < s2, а дальнейшее увеличение от S3 до заданного s2 происходит скачком в контактном разрыве (эта картина относится к типу, изображенному ниже на рис. 78, б; предполагается выполненным неравенство (86,2))'). Вопрос о том, чем определяется отбор одного из двух решений в конкретных гидродинамических задачах, не ясен. Если отбирается распадное решение, то это означало бы, что неустойчивость ударной волны с самопроизвольным усилением поверхностной ряби вообще не осуществляется. По-видимому, однако, такой отбор не может быть связан именно с этой неустойчивостью, поскольку неоднозначность решения не ограничена условиями (90,12—13) 2). Задачи
1. На ударную волну падает сзади (со стороны сжатого газа) нормально к ней плоская звуковая волна. Определить коэффициент отражения звука. Решение. Рассматриваем процесс в системе координат, в которой ударная волна покоится, а газ движется через нее в положительном направлении оси х; падающая звуковая волна распространяется в отрицательном направлении оси х. При нормальном падении (а потому и отражении) в отраженной энтропийной волне скорость 6v(3HI) = 0. Возмущение давления: бр = = 6р(зв) + 6р(0), где индекс (0) относится к падающей, а индекс (зв) — к отраженной звуковым волнам. Для скорости 8vx = 8v имеем 60=-^-(бр(зв> -6р<°>) (разность вместо суммы возникает ввиду противоположных направлений распространения обоих волн). Второе из граничных условий (90,7) имеет прежний вид (но в нем теперь 6V = 6V<°> + 6К<38> + 6V<3">); с учетом (90,8) и формулы (85,6) переписываем его как 5o = _-L=i^-(6p<-> +бр<°>).
Приравняв друг другу оба выражения 6у, получим для искомого отношения амплитуд давления в отраженной и падающей звуковых волнах: 6р(зв) 1 — 2М2 — h бр(0> 1 + 2М? - h (где М2 = v2/с2). Оно обращается в бесконечность на верхней границе области (90,17) Для политропного газа h = — Mj~2. При слабой интенсивности ударной волны (р2— pi < pi) отношение (1) стремится к нулю как (р% — рО2, а в обратном случае большой интенсивности стремится к постоянному пределу 6р(зв>Уу - Уд (y - о 6р<°> ~ Уу" + V2 (y - о ' 2. На ударную волну падает спереди, нормально к ней, плоская звуковая волна. Определить коэффициент прохождения звука'). Решение. Возмущение в газе / перед ударной волной 6p,=6p(0>. 6V,=6V,0> =--^-6p,, bv^^-bpu а в газе 2 позади нее: 6р2 = бр'зв), 6У2 = 6V(3B> + 6V(3HT), 6о2 = -^-др2 (индексы (0), (зв), (энт) относятся к падающей звуковой и к прошедшим звуковой и энтропийным волнам). Возмущения бр2 и 6V2 связаны друг с другом соотношением, следующим из уравнения ударной адиабаты: если
■(^b+[-T©,t(?)J'" (индекс Н у производных указывает, что они берутся вдоль адиабаты Гю-гонио1)). Граничное условие (90,7) заменяется теперь на
2 L р2 — Pi к 1 — V2 J = ~ тг[6p? ~6р> ~jl (flfj ~6v,)1- Приравняв два выражения для ба2 — 6i»i, получим для искомого отношения амплитуд в прошедшей и падающей звуковых волнах:
6р(зв) (1 + М,)2 + <? 6р<й) 1 + 2М2 — h
где h имеет прежнее значение, а 9 = У" Для политропного газа Y-1 (М?-1)2 <7 = — „ Y + 1 М2 и коэффициент прохождения: 6р(зв)(1+М,)2 6р<0) 1 4- 2М2 + м,-2[ Y + iv м,;J При слабой интенсивности ударной волны отсюда получается
а в обратном случае большой интенсивности: «РЬв) „ 1 р2 бр(0) Y+V2v(V-l) Pi В обоих случаях амплитуда давления в прошедшей звуковой волне возрастает по сравнению с давлением в падающей волне.
§ 91. Распространение ударной волны по трубе Рассмотрим распространение ударной волны по среде, заполняющей длинную трубку с переменным сечением. Наша цель состоит при этом в выяснении влияния, оказываемого изменением площади ударной волны на ее скорость (G. В. Whitham, 1958).
Будем считать, что площадь S(x) сечения трубки лишь медленно меняется вдоль ее длины (ось я) — мало на расстояниях порядка ширины трубки. Это дает возможность применить приближение (его называют гидравлическим), которое уже было использовано в § 77: можно считать все величины в потоке постоянными вдоль каждого поперечного сечения трубки, а скорость— направленной вдоль ее оси; другими словами, течение рассматривается как квазиодномерное. Такое течение описывается уравнениями ж + °£ + т%г=°- <»и> *+•£-«'($+»£)=<>■ <9'.2>
Первое из них — уравнение Эйлера, второе — уравнение адиаба-тичности, а третье — уравнение непрерывности, представленное в виде (77,1). Для выяснения интересующего нас вопроса достаточно рассмотреть трубку, в которой изменение площади S(x) не только медленно, но и по абсолютной величине остается относительно малым на протяжении всей длины. Тогда будут малы и связанные с непостоянством сечения возмущения потока, и уравнения (91,1—3) могут быть линеаризованы. Наконец, должны быть поставлены начальные условия, исключающие появление каких-либо посторонних возмущений, которые могли бы повлиять на движение ударной волны; нас интересуют только возмущения, связанные с изменением S(x). Эта цель будет достигнута, если принять, что ударная волна первоначально движется с постоянной скоростью по трубе постоянного сечения, и площадь сечения начинает меняться только вправо от некоторой точки (которую примем за х — 0). Линеаризованные уравнения (91,1—3) имеют вид dt ' дх р дх dbp. дбр 2 (дбр,, <Э6р dt + V дх дбр. дбр, dbv. pv d6S п -дТ + ',-дГ+Р-дГ + £8--дГ = 0' где буквы без индекса обозначают постоянные значения величин в однородном потоке в однородной части трубки, а символ б обозначает изменение этих величин в трубке переменного сечения. Умножив первое и третье из этих уравнений соответственно на рс и с2 и сложив затем все три уравнения, напишем следующую их комбинацию: + + *4с-](*Р + = - ^ ^- (91,4) Общее решение этого уравнения дается суммой общего решения однородного уравнения и частного решения уравнения с правой частью. Первое есть F(x— vt — ct), где F— произвольная функция; оно описывает звуковые возмущения, приходящие слева. Но в однородной области, при х < 0, возмущений нет; поэтому надо положить F = 0. Таким образом, решение сводится к интегралу неоднородного уравнения: 6p + Pc8v = --^^-. (91,5) Ударная волна движется слева направо со скоростью v\ > С\ по неподвижной среде с заданными значениями рь р;. Движение же среды позади ударной волны (среда 2) определяется решением (91,5) во всей области трубки слева от точки, достигнутой разрывом к данному моменту времени. После прохождения волны все величины в каждом сечении трубки остаются постоянными во времени, т. е. равными тем значениям, которые qhh получили в момент прохождения разрыва: давление р2, плотность р2 и скорость vi — v2 (в соответствии с принятыми в этой главе обозначениями, v% обозначает скорость газа относительно движущейся ударной волны; скорость же его относительно стенок трубки есть тогда v\— v2). В этих обозначениях (и снова выделив переменные части этих величин) равенство (91,5) запишем в виде Ц. =, _ S P2(^i-f2)c2 Все величины 8vi, 6i>2, бр2 можно выразить через одну из них, например, 6ui. Для этого пишем варьированные соотношения (85,1—2) на разрыве (при заданных р, и pi): р^о, = v26o2 + p26v2, 2/ (во, - 6v2) = Ьр2 + v\bp2 (где / = p\V\ — р2и2 — невозмущенное значение потока); к ним надо еще присоединить соотношение
где производная берется вдоль адиабаты Гюгонио. Вычисление приводит к следующему окончательному соотношению, связывающему изменение bv\ скорости ударной волны относительно неподвижного газа перед ней, с изменением 8S площади сечения трубки: 1 6S Pi — Рг + е2 Г 1 +2t)2/c2 — h-\ 1 Г+А j' (91'7) Коэффициент при квадратной скобке в (91,7) положителен. Поэтому знак отношения 8vi/8S определяется знаком выражения в этой скобке. Для всех устойчивых ударных волн этот знак положителен, так что ovi/oS < 0. Но при выполнении какого-либо из условий (90,12—13) гофрировоч-ной неустойчивости выражение в скобках становится отрицательным, так что 8v\/oS > 0. Этот результат дает возможность наглядного истолкования происхождения неустойчивости. На рис. 62 изображена «гофрированная» поверхность ударной волны, перемещающаяся направо; стрелками схематически показано направление линий тока. При перемещении ударной волны на выдавшихся вперед участках поверхности площадь 6S растет, а на отставших участках— уменьшается. При 6i>i/6S < 0 это приводит к замедлению выступивших участков и ускорению отставших, так что поверхность стре- Рас. 62 мится сгладиться. Напротив, при 8vi/8S > 0 возмущение формы поверхности будет усиливаться: выступающие участки будут уходить все дальше, а отставшие — все более отставать').
Дата добавления: 2014-11-07; Просмотров: 408; Нарушение авторских прав?; Мы поможем в написании вашей работы! Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет |