КАТЕГОРИИ: Архитектура-(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 54 страницаПри заданном начальном термодинамическом состоянии газа (т. е. заданных р\, V\) ударная волна определяется всего одним каким-либо параметром: если, например, задать давление р2 за волной, то по адиабате Гюгонио определится V2, а затем по формулам (85,4) и (85,6) — плотность потока / и скорости vi и Vi. Напомним, однако, что мы рассматриваем здесь ударную волну в системе координат, в которой газ движется нормально к ее поверхности. Если же учесть возможность расположения ударной волны под косым углом к направлению потока, то понадобится еще один параметр, например, значение касательной к ее поверхности составляющей скорости. Укажем здесь на следующее удобное графическое истолкование формулы (85,6). Если соединить хордой точку ри Vi на ударной адиабате (рис. 53) с некоторой произвольной точкой р2, V2 на ней, то (р2 — Pi) / (V2— Vi)= —j2 есть не что иное, как тангенс угла наклона этой хорды к оси абсцисс (к ее положительному направлению). Таким образом, значение /, а с ним и скорости ударной волны, определяется в каждой точке ударной адиабаты углом наклона хорды, проведенной в эту точку из начальной точки. Наряду с другими термодинамическими величинами в ударной волне испытывает разрыв также и энтропия. В силу закона возрастания энтропии последняя для газа может лишь возрастать при его движении. Поэтому энтропия s2 газа, прошедшего через ударную волну, должна быть больше его начальной энтропии Sb s2>si. (85,11) Мы увидим ниже, что это условие налагает существенные ограничения на характер изменения всех величин в ударной волне. Подчеркнем здесь следующее обстоятельство. Наличие ударных волн приводит к возрастанию энтропии при таких движениях, которые можно рассматривать во всем пространстве как движение идеальной жидкости, не обладающей вязкостью и теплопроводностью. Возрастание энтропии означает необратимость движения, т. е. наличие диссипации энергии. Таким образом, разрывы представляют собой механизм, который приводит к диссипации энергии при движении идеальной жидкости. В связи с этим для движения тел в идеальной жидкости, сопровождающегося возникновением ударных волн, не имеет места парадокс Даламбера (§ 11) — при таком движении тело испытывает силу сопротивления. Разумеется, истинный механизм возрастания энтропии в ударных волнах заключен в диссипативных процессах, происходящих в тех весьма тонких слоях вещества, которые в действительности представляют собой физические ударные волны (см. § 93). Замечательно, однако, что величина этой диссипации целиком определяется одними лишь законами сохранения массы, энергии и импульса, примененными к обеим сторонам этих слоев: их ширина устанавливается как раз такой, чтобы дать требуемое этими законами сохранения увеличение энтропии. Возрастание энтропии в ударной волне оказывает еще и другое существенное влияние на движение: если движение газа впереди ударной волны потенциально, то за ней оно, вообще говоря, становится вихревым; мы вернемся к этому обстоятельству в § П4. § 86. Ударные волны слабой интенсивности Рассмотрим ударную волну, в которой все величины испытывают лишь небольшой скачок; о таких разрывах мы будем говорить как об ударных волнах слабой интенсивности. Преобразуем соотношение (85,9), производя в нем разложение по степеням малых разностей s2— Si и р2 —Рь Мы увидим, что при таком разложении в (85,9) сокращаются члены первого и второго порядков по р2— pi\ поэтому необходимо производить разложение по р2 — р\ до членов третьего порядка включительно. По разности же s2— s\ достаточно разложить до членов первого порядка. Имеем: = (|^)р (* - st) + (|^-)s (р2 - р.) +
Но согласно термодинамическому соотношению dw — Tds -f* + Vdp имеем для производных: (*).-'• (■£).-"• Поэтому Щ — Щ = 7", (s2 — si) + Vi (р2 — Pi) +
Объем У2 достаточно разложить только по р2 — рь поскольку во втором члене уравнения (85,9) уже имеется малая разность •р2 — pi и разложение по s2—si дало бы член порядка (s2 — — Si) (Р2 — Pi), не интересующий нас. Таким образом, "•-"-(£).<* "'■>+т(3).(*-''- Подставляя эти разложения в (85,9), получим следующее соотношение: s2-s, = -!— (~\ (p2-pi)3. (86,1)
Таким образом, скачок энтропии в ударной волне слабой интенсивности является малой величиной третьего порядка по сравнению со скачком давления. Адиабатическая сжимаемость вещества — (дV/dp)s практически всегда падает с увеличением давления, т. е. вторая произ- водная ') (!?■). >°- ^ Подчеркнем, однако, что это неравенство не является термодинамическим соотношением и, в принципе, возможны его нарушения2). Как мы неоднократно увидим ниже, в газодинамике знак производной (86,2) весьма существен; в дальнейшем мы будем всегда считать его положительным. Проведем через точку / (рь Vi) на р, К-диаграмме две кривые— ударную адиабату и адиабату Пуассона. Уравнение адиабаты Пуассона есть s2 — Si = 0. Из сравнения этого уравнения с уравнением (86,1) ударной адиабаты вблизи точки / видно, что обе кривые касаются в этой точке, причем имеет место касание второго порядка — совпадают не только первые, но и вторые производные. Для того чтобы выяснить взаимное расположение обеих кривых вблизи точки /, воспользуемся тем, что согласно (86,1) и (86,2) при р2> Р\ на ударной адиабате должно быть Sa > si, между тем как на адиабате Пуассона остается s2 = si-Поэтому абсцисса точки на ударной адиабате должна быть при той же ординате р2 больше абсциссы точки на адиабате Пуассона. Это следует из того, что согласно известной термодинамической формуле V ds)р ср \ дТ)р
энтропия растет с увеличением объема при постоянном давлении — для всех тел, которые расширяются при нагревании, т. е. у которых (дУ/дТ)р>0. Аналогично убеждаемся в том, что ниже точки / (т. е. при рй < pi) абсциссы точек адиабаты Пуассона должны быть больше абсцисс ударной адиабаты. Таким образом, вблизи точки своего касания обе кривые расположены указанным на рис. 55 образом (##' — ударная адиабата, а РР' —• адиабаты Пуассона)3), причем в силу (86,2) обе обращены вогнутостью вверх. При малых р2—pi и V2—Vi формулу (85,6) можно написать в первом приближении в виде г—т. (мы пишем здесь производную при постоянной энтропии, имея в виду, что касательные к адиабатам Пуассона и ударной в точке 1 совпадают). Далее, скорости vi и v2 в том же приближении одинаковы и равны -я'-л/-^(*).-л/Ш. Но это есть не что иное, как скорость звука с. Таким образом, скорость распространения ударных волн слабой интенсивности совпадает в первом приближении со скоростью звука: v = c. (86,3) Рис. 56 Из полученных свойств ударной адиаба- ты в окрестности точки 1 можно вывести ряд существенных следствий. Поскольку в ударной волне должно выполняться условие s2 > s\, то должно быть и Р2 > Pl, т. е. точки 2 (р2, V2) должны находиться выше точки /. Далее, поскольку хорда 12 идет круче касательной к адиабате в точке / (рис. 53), а тангенс угла наклона этой касательной равен производной {dpi/dV\)si, имеем: г->-Ш,- Умножая это неравенство с обеих сторон на V\, находим:
где Ci — скорость звука, соответствующая точке /. Таким образом, Vi > С\. Наконец, из того, что хорда 12 расположена менее круто, чем касательная в точке 2, аналогичным образом следует, что V2 < С2 '). ') Последняя аргументация применима только вблизи точки 1, где тангенс угла наклона касательной к ударной адиабате в точке 2 отличается от производной {dp2idVi)s лишь на величину второго порядка малости. § 87] НАПРАВЛЕНИЕ ИЗМЕНЕНИЯ ВЕЛИЧИН В УДАРНОЙ ВОЛНЕ 463
Упомянем еще, в заключение, что при (d2V/dp2)s <. О из условия S2 > Si для ударных волн слабой интенсивности следовало бы р2 < Рь а для скоростей — те же неравенства v\ > С\, о2 < с2. § 87. Направление изменения величин в ударной волне Таким образом, в предположении положительности производной (86,2) для ударных волн слабой интенсивности можно весьма просто показать, что условие возрастания энтропии с необходимостью приводит также и к неравенствам р2 > Рь (87,1) oi > си v2 < с2. (87,2) Из замечания, сделанного по поводу формулы (85,6) следует, что если р2 > Рь то V2 < Vu (87,3) а поскольку / == V\/V\ = v2/V2, то и') vx > о2. (87,4) Неравенства (87,1) и (87,3) означают, что при прохождении газа через ударную волну происходит его сжатие — его давление и плотность возрастают. Неравенство oi > С\ означает, что ударная волна движется относительно находящегося перед ней газа со сверхзвуковой скоростью; ясно поэтому, что в этот газ не могут проникнуть никакие исходящие от ударной волны возмущения. Другими словами, наличие ударной волны вовсе не сказывается на состоянии газа впереди нее. Покажем теперь, что все неравенства (87,1—4) справедливы и для ударных волн произвольной интенсивности — при том же предположении о знаке производной (d2V/dp2)s2).
Величина j2 определяет наклон хорды, проведенной из начальной точки ударной адиабаты / в произвольную точку 2 (—j2 есть тангенс угла наклона этой хорды к оси V). Покажем, прежде всего, что направление изменения этой величины при перемещении точки 2 вдоль адиабаты однозначно связано с направлением изменения энтропии s2 при том же перемещении. Продифференцируем соотношения (85,5) и (85,8) по величинам, относящимся к газу 2 при заданном состоянии газа /. Это значит, что дифференцируются р2> V2, w2 и / при заданных значениях рь V\, w\. Из (85,5) получаем: dp2 + j2 dV2 = {V\ - V2) d (j2), (87,5) а из (85,8): dw2 + j2V2 dV2 = \ (V2 - Vl) d {f) или, раскрыв дифференциал dw2: т2 ds2 + v2 {dP2 + f dv2)=j {v] - vi) d (/). Подставив сюда dp2 + j2dV2 из (87,5), получим соотношение T2ds2 = ±{Vx-V2)2d(f). (87,6)
Отсюда видно, что d(j2)/ds2>0, (87,7) т. е. j2 и s2 меняются в одинаковом направлении. Дальнейшие рассуждения имеют своей следующей целью показать, что на ударной адиабате не может быть точек, в которых бы она касалась проведенной из точки / прямой (как это имело бы место в точке О на рис. 56. В такой точке угол наклона хорды (проведенной из точки /) имеет минимум, а /2—соответственно максимум, и потому d(j2)/dp2 = 0.
_ Из соотношения (87,6) видно, что в таком V случае будет и ds2/dp2 = 0. Далее, вычислим производную d(j2)/dp2 в произвольной точке ударной адиабаты. Подставив в соотношение (87,5) дифференциал dV2 в виде
взяв для ds2 выражение (87,6) и разделив все равенство на dp2, получим dpi (V, - V2) [l
j2 (У, - Уг) (dVi \ ds2 /pj (87,8) Отсюда видно, что обращение этой производной в нуль влечет за собой также и равенство 6V, - -2 V 3Pi Js, с2 т. е. U2 = с2. Обратно, из равенства у2 = с2 следует, что производная d(j2)/dp2 = 0; последняя могла бы не обратиться в нуль, лишь если бы вместе с числителем в (87,8) обратился бы в нуль также и знаменатель; но выражения в числителе и знаменателе представляют собой две различные функции точки 2 на ударной адиабате, их одновременное обращение в нуль могло бы произойти лишь чисто случайно и потому невероятно1). Таким образом, все три равенства чрт=*> ж-0' v*=c* <87-9> являются следствиями друг друга и имели бы место одновременно в точке О на кривой рис. 56 (имея в виду последнее из этих равенств, будем условно называть такую точку звуковой). Наконец, для производной от (v2/c2)2 в этой точке имеем
dp2\4) dP.2['\др2л,J Vдр22А/ Ввиду предполагаемой везде положительности производной (d2V/dp2)s имеем, следовательно, в звуковой точке: Тр-Л<*- <87'10> Теперь уже легко доказать невозможность существования звуковой точки на ударной адиабате. В точках, лежащих вблизи начальной точки / над ней, имеем и2 < с2 (см. конец предыдущего параграфа). Поэтому равенство v2 = с2 может быть достигнуто лишь при увеличении у2/с2; другими словами, в звуковой точке должно было бы быть d(vi/c2) /dp2 > 0, между тем как согласно (87,10) мы имеем как раз обратное неравенство. Аналогичным образом можно убедиться в невозможности обращения v2/c2 в единицу и на нижней части ударной адиабаты, под точкой /. Имея в виду доказанную таким образом невозможность существования звуковых точек, можно заключить непосредственно из графика ударной адиабаты, что угол наклона хорды 12 уменьшается при передвижении точки 2 вверх по кривой, а /2 соот- ') Подчеркнем, во избежание недоразумений, что сама производная d('f)Idp2 не является еще одной независимой функцией точки 2; выражение (87,8) есть ее определение. ветственно монотонно возрастает; ввиду неравенства (87,7) отсюда следует, что монотонно возрастает и энтропия s2. Таким образом, при соблюдении необходимого условия s2 > Sj будет и р2 > Рь Легко, далее, убедиться в том, что на верхней части ударной адиабаты справедливы также и неравенства v2 < с2> v\ > С\. Первое следует прямо из того, что оно справедливо вблизи точки /, а сделаться равным единице отношение v2/c2 нигде не может. Второе следует из того, что ввиду невозможности такого перегиба адиабаты, какой изображен на рис. 56, всякая хорда из точки / в находящуюся над ней точку 2 расположена более круто, чем касательная к ударной адиабате в точке /. Таким образом, на верхней части ударной адиабаты выполняются условие s2 > s\ и все три неравенства (87,1—2). Наоборот, на нижней части адиабаты все эти условия не выполняются. Следовательно, все эти условия оказываются эквивалентными друг другу и выполнение одного из них автоматически влечет за собой также и выполнение всех остальных. Напомним лишний раз, что в изложенных рассуждениях все время предполагалось выполненным условие положительности производной (d2V/dp2)s. Если эта производная могла бы менять знак, то из необходимого термодинамического неравенства s2 > Si уже нельзя было бы сделать никаких универсальных заключений о неравенствах для остальных величин.
§ 88. Эволюционность ударных волн Вывод неравенств (87,1—4) в §§ 86, 87 был связан с определенным предположением о термодинамических свойствах среды —с положительностью производной (d2V/dp2)s. Весьма существенно, однако, что неравенства vi > си v2 < с2 (88,1) для скоростей могут быть получены также и из совершенно иных соображений, показывающих, что ударные волны с нарушенными условиями (88,1) все равно не могли бы существовать, даже если бы это не противоречило изложенным выше чисто термодинамическим соображениям ').
Именно необходимо исследовать еще вопрос об устойчивости ударных волн. Наиболее общее необходимое условие устойчивости состоит в требовании, чтобы любое бесконечно малое воз- мущение начального (в некоторый момент t — 0) состояния при- Произвольное начальное малое возмущение определяется некоторым числом независимых параметров. Дальнейшая же эволюция возмущения определяется системой линеаризованных граничных условий, которые должны удовлетворяться на поверхности разрыва. Поставленное выше необходимое условие устойчивости будет выполнено, если число этих уравнений совпадает с числом содержащихся в них неизвестных параметров — тогда граничные условия определяют дальнейшее развитие возмущения, которое при малых t > 0 останется малым. Если же число уравнений больше или меньше числа независимых параметров, то задача о малом возмущении не имеет решений вовсе или имеет их бесконечное множество. Оба случая свидетельствовали бы о неправомерности исходного предположения (малость возмущения при малых г) и, таким образом, противоречили бы поставленному требованию. Сформулированное таким образом условие называют условием эволюционности течения.
Рассмотрим возмущение ударной волны, представляющее собой ее бесконечно малое смещение в направлении, перпендикулярном ее плоскости1). Оно сопровождается бесконечно малым возмущением также и других величин — давления, скорости и т. д. газа по обеим сторонам поверхности разрыва. Эти возмущения, возникнув вблизи волны, будут затем распространяться от нее, переносясь (относительно газа) со скоростью звука; это не относится лишь к возмущению энтропии, которое будет переноситься только с самим газом. Таким образом, произвольное возмущение данного типа можно рассматривать как совокупность звуковых возмущений, распространяющихся в газах / и 2 по обе стороны ударной волны, и возмущения энтропии; последнее, перемещаясь вместе с газом, будет, очевидно, существовать лишь в газе 2 позади ударной волны. В каждом из звуковых возмущений изменения всех величин связаны друг с другом определенными соотношениями, следующими из уравнений движения (как в любой звуковой волне; § 64); поэтому каждое такое возмущение определяется всего лишь одним параметром. Подсчитаем теперь число возможных звуковых возмущений. Оно зависит от относительной величины скоростей газа v\, v2 и скоростей звука С\, с2. Выберем направление движения газа (со стороны / на сторону 2) в качестве положительного направления оси х. Скорость распространения возмущения в газе / относительно неподвижной ударной волны есть u\—V\±C\, а в газе 2 «2 = v2 ± с2. Тот факт, что эти возмущения должны распространяться по направлению от ударной волны, означает, что должно быть «] < 0, и2 > 0.
Предположим, что v\ > си о2 < с2. Тогда ясно, что оба значения щ = v\ ± с\ будут положительными, а из двух значений м2 будут положительными лишь v2 + -f с2. Это значит, что в газе / вообще не сможет быть интересующих нас звуковых возмущений, а в газе 2 — всего одно, распространяющееся относительно самого газа со скоростью -}- с2. Аналогичным образом производится подсчёт в других случаях. Результат изображен на рис. 57, указываемую стрелкой сторону. Каждое же звуковое возмущение определяется, как было выше указано, одним параметром. Кроме того, во всех четырех случаях имеется еще по два параметра: параметр, определяющий распространяющееся в газе 2 возмущение энтропии, и параметр, определяющий самое смещение ударной волны. Для каждого из четырех случаев на рис. 57 цифрой в кружке указано получающееся таким образом полное число параметров, определяющих произвольное возмущение, возникающее при смещении ударной волны. С другой стороны, число необходимых граничных условий, которым должно удовлетворять возмущение на поверхности разрыва, равно трем (условия непрерывности потоков массы, энергии и импульса). Во всех изображенных на рис. 57 случаях, за исключением лишь первого, число имеющихся независимых параметров превышает число уравнений. Мы видим, что эволю-ционны лишь ударные волны, удовлетворяющие условиям (88,1). Эти условия, таким образом, необходимы для существования ударных волн, вне зависимости от термодинамических свойств среды. Искусственно созданный разрыв, не удовлетворяющий этим условиям, немедленно распался бы на другие разрывы1). Эволюционная ударная волна устойчива по отношению к рассмотренному типу возмущений и в обычном смысле этого слова. Если искать смещение ударной волны (а с ним и возмущения всех остальных величин) в виде, пропорциональном е~ш, то заранее очевидно, что однозначно определяемое граничными условиями значение со может быть только нулем — уже хотя бы из тех соображений, что в задаче нет никаких параметров размерности обратного времени, которые могли бы определить отличное от нуля значение со. Мы вернемся к вопросу об устойчивости ударных волн в § 90.
§ 89. Ударные волны в политропном газе Применим полученные в предыдущих параграфах общие соотношения к ударным волнам в политропном газе. Тепловая функция такого газа дается простой формулой Ll-vt — ^_ (у+ 1) pi + (у — 1>ра (Y-l)Pi + (Y+»)/»» " (89,1)
По этой формуле можно определить по трем из величин ри Vi, р2, V2 четвертую. Отношение V2/V\ является монотонно убывающей функцией отношения pi/pi, стремящейся к конечному пределу (у—1)/(у+1). Кривая, изображающая зависимость между р2 и V2 при заданных р\, V\ (ударная адиабата), представлена на рис. 58. Это — равнобочная гипербола с асимптотами
Уг у — 1 £± Y—1 Pi Y+ 1 ') Во всех перечисленных на рис. 57 неэволюционных случаях возмущение недоопределено — число произвольных параметров превышает число уравнений Упомянем, что в магнитной гидродинамике ударные волны могут быть неэволюционными в силу кзк недоопределенности, так и переопределенности возмущений (см. VIII, § 73). Реальным смыслом обладает, как мы знаем, только верхняя часть кривой над точкой V2/Vi = Р2/Р1 = 1, изображенная на рис. 58 (для у = 1,4) сплошной линией. Для отношения температур с обеих сторон разрыва имеем согласно уравнению состояния термодинамически идеального газа Т2/Тх = p2V2/piVi, так что ll — H. Г <У + !>Р' + (V —0Р2 1 /oq ™ Ti р, L(Y-1)P. + (Y+1)P2J- ( '> Для потока / получаем из (85,6) и (89,1): ;2 _ (Y — 0 Pi + (у+ О Р2 (яа о\ /------------------- 21Л----------- (оУ.-э) и отсюда для скорости распространения ударной волны относительно газов впереди и позади нее: .2 °i = Т" [(V - 1) Pi + (Y + 1) Р2] -= |г [(V ~ D + (Y + 1). (89 4) „ У, [(У+1)Р! + (У-1)Р2]2 с2 Г, п. /,пр<1 pl°-2- и для разности скоростей:
[(Y-l)Pi + (Y+Dp2],/2 V В применениях полезны формулы, выражающие отношения плотностей, давлений и температур в ударной волне через число Mi = Vi/cu эти формулы без труда выводятся из полученных выше соотношений и гласят: Р2 " (Y + 1)M' (89,6) р, о2 (v_l)Mf + 2
£2_ Pl Г2 [2YM2-(v-l)][(Y-l)M2 + 2] Г, (y + 02М2 ' 1 ' ' Число же М2 = у2/с2 выражается через число Mi посредством 2 2 + (у-1)М2 М? = —-5^ —. (89,9) 2YM2-(y- 1) v Это соотношение симметрично относительно Mi и М2, как это •становится очевидным, если записать его в виде уравнения 2yMiMl - (y — 1) (М? + Ml) = 2. Выпишем предельные формулы для ударных волн очень большой интенсивности (требуется, чтобы бьмо (7 — 1)р2"> >(V— 1)pi). Имеем из (89,1—2):
Дата добавления: 2014-11-07; Просмотров: 368; Нарушение авторских прав?; Мы поможем в написании вашей работы! Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет |