1 МЕХАНИКА СПЛОШНЫХ СРЕД1 МЕХАНИКА СПЛОШНЫХ СРЕД2 Газовая динамика Автор курса лекций: Породнов Борис Трифонович, д. ф.-м. н., профессор кафедры молекулярной. - презентация

1 1 МЕХАНИКА СПЛОШНЫХ СРЕД

2 2 Газовая динамика Автор курса лекций: Породнов Борис Трифонович, д. ф.-м. н., профессор кафедры молекулярной физики УГТУ-УПИ Екатеринбург 2007 Численные решения задач механики сплошных сред

3 3 Модуль 7 Газовая динамика Численные методы…Лекция 27

4 4 Взаимодействие ударных волн и волн разрежения. Истечение газа из сужающегося насадка. Сопло Лаваля в нерасчетном режиме. Обтекание пластинки сверхзвуковым потоком под малым углом атаки. Лекция Численные методы…Лекция 27

5 5 Цели изучения: Изучить отражение ударной волны от стенки, пересечение двух ударных волн Рассмотреть взаимодействие ударной волны, и волны с тангенциальным разрывом и процесс возникновения ударной волны в пограничном. Изучить сопло Лаваля в нерасчетном режиме. Рассмотреть истечение газа из сужающегося насадка и обтекание пластинки сверхзвуковым потоком под малым углом атаки Численные методы…Лекция 27

6 6 Содержание Взаимодействие ударных волн и волн разрежения Отражение ударной волны от стенки Пересечение двух ударных волн Взаимодействие ударной волны, и волны с тангенциальным разрывом Возникновение ударной волны в пограничном слое Истечение газа из сужающегося насадка Сопло Лаваля в нерасчетном режиме Обтекание пластинки сверхзвуковым потоком под малым углом атаки Численные методы…Лекция 27

7 ВЗАИМОДЕЙСТВИЕ УДАРНЫХ ВОЛН И ВОЛН РАЗРЕЖЕНИЯ Ударные волны и волны разрежения в движущемся газе могут накладываться друг на друга, пересекаться или, падая на твердую стенку, ограничивающую движение, отражаться от нее или, встречая тангенциальный разрыв, также испытывать отражение или преломление. 7 Рис Рассмотрим качественно отражение и пересечение ударных волн и волн разрежения с тангенциальным разрывом.

8 Отражение ударной волны от стенки Пусть на твердую стенку падает косая ударная волна (рис.12.26). Можно представить, что где-то вдали от стенки в сверхзвуковой поток помещен клин, от которого отошла рассматриваемая ударная волна. Для определения движения газа за падающей ударной волной можно воспользоваться ударной полярой (рис.12.27).рис.12.26рис Рис

9 9 Выберем на номограмме поляру, соответствующую заданному числу М 1. Проведем из начала координат прямую под известным углом χ, равным углу полураствора клина. Тогда, как известно, пересечение этой линии в точке С 1 с ударной полярой определит отрезок ОС 1, равный модулю скорости υ 2 за падающей ударной волной, а также угол поворота потока χ. Далее необходимо через точки В 1 С 1 провести линию и на нее опустить перпендикуляр из начала координат. Угол между этим перпендикуляром и осью υ 2х определяет угол наклона косой ударной волны φ 1. Зная М 1 и φ 1 по формуле ( ) определим М 2 – число Маха за падающей ударной волной (М 2

10 Газ за падающей ударной волной должен двигаться вдоль твердой стенки. Это значит, что от точки пересечения падающей ударной волны со стенкой (точка О) должна отойти также косая ударная волна, в которой газ должен повернуть на тот же самый угол χ и двигаться вдоль стенки. Скорость за отраженной ударной волной может быть найдена, если определена точка С 2 – точка пересечения линии ОС 1 с ударной полярой М 2. Если теперь из точки В 2 провести через точку С 2 линию, а на нее из начала координат опустить перпендикуляр, то линия ОD 2 определит угол φ 2 – угол отражения падающей ударной волны. Зная М 2 и φ 2, можно определить число Маха за отраженной волной. М 3 φ 1. Пользуясь формулой ( ), можно показать аналитически. Т.о., для косой ударной волны φ пад

11 11 Таким образом, падающая на твердую стенку ударная волна отражается от нее тоже в виде косой ударной волны, причем угол отражения не равен углу падения (см. рис ). Только очень слабые (звуковые) ударные волны М 2 =М 1 будут отражаться от твердой поверхности зеркально (φ 2 =φ 1 ). отражаетсярис Если падающая ударная волна настолько интенсивная, что φ 2 становится больше φ max, то в этом случае отражение будет неправильным (рис ).рис Это значит, что от точки О 1 может отойти тангенциальный разрыв, т.е. линия, на которой все термодинамические параметры прерывны, нормальная компонента скорости равна нулю, а тангенциальная к поверхности разрыва скорость газа изменяется скачком. Между стенкой и точкой отражения возникает участок прямой ударной волны ОО 1. Скорость газа, прошедшего через отраженную волну, может быть и не равна скорости газа, прошедшего прямую ударную волну.

12 12 Рассмотрим взаимодействие в газе двух ударных волн. Снова можно представить, что они образовались от двух клинов, обтекаемых сверхзвуковым потоком газа (рис.12.29).взаимодействие в газе двух ударных волнрис Если ударные волны одинаковой интенсивности, то в результате взаимодействия они отразятся от линии симметрии потока с одинаковыми углами. Если ударные волны будут иметь достаточно большую интенсивность, то пересечение будет неправильным. Образуется участок прямой ударной волны; от точек О и О отойдут тангенциальные разрывы. Это неправильное пересечение. пересечениебудет неправильным Пересечение ударных волн. Если ударные волны одинаковой интенсивности, то линия тока, проходящая через точку пересечения (точка О), не должна прерывать каких-либо изменений параметров потока и скоростей, и ее можно рассматривать как твердую стенку. Можно воспользоваться предыдущим рассмотрением отражения падающей ударной волны от твердой стенки.

13 13 Если интенсивность пересекающихся ударных волн неодинакова, то направление поворота потока определяется так: поток, прошедший ударную волну большей интенсивности, повернет на больший угол, и линия тангенциального разрыва повернет в сторону менее интенсивной ударной волны, а, следовательно, повернет и весь газодинамический поток с меньшей скоростью (рис ).рис Рис Рис.12.31

14 Развитие следа при возрастающем числе М с взаимодействием ударных волн можно наблюдать на рис а и рис б. По мере того как число Маха невозмущенного потока возрастает, приближаясь к единице, след за круговым цилиндром переходит от периодического схождения вихрей, характерного для малых скоростей, к квазистационарному следу с шейкой, устанавливающемуся при сверхзвуковых скоростях. Взаимодействие косых ударных волн с турбулентным следом за цилиндром приводит к образованию «неправильной» косой ударной волны с участком прямой за кормой цилиндра. [Dyment, Gryson, Ducruet, 1980]рис а рис б 14

15 15 М=0.80М=0.90 Рис а

16 16 М=0.95М=0.98 Рис б

17 17 Предположим, что в движущемся газе есть тангенциальный разрыв. Под тангенциальным разрывом понимают такую поверхность, на которой все термодинамические величины непрерывны, а тангенциальные компоненты скорости терпят разрыв. Непрерывность термодинамических величин следует из уравнений сохранения на поверхности разрыва см (12.1.5) На поверхности: P 1 =P 2, ρ 1 = ρ 2, T 1 =T 2, υ п1 = υ п2 – в системе координат, связанной с поверхностью разрыва, но υ 1t не равно υ 2t Пусть на этот тангенциальный разрыв падает косая ударная волна. В газе, прошедшем через ударную волну, давление возрастает. Но в области течения за тангенциальным разрывом давление P 0 не изменяется. Следовательно, должен существовать механизм, который понижал бы давление с Р 2 до Р 3 =Р. Это может быть только волна разрежения (давление в области, ограниченной характеристиками 1 и 2, понижается). В соответствии с изменением направления скорости в ударной волне и волне разрежения, поворачивает и скорость потока и тангенциальный разрыв Взаимодействие ударной волны и волны с тангенциальным разрывом

18 18 Косая ударная волна, падающая на тангенциальный разрыв, отражается в виде волны разрежения. И наоборот, волна разрежения, падая на тангенциальный разрыв, должна отразиться в виде косой ударной волны.

19 Вспомним явление отрыва пограничного слоя при М

20 Возникновение ударных волн в пограничном слое при дозвуковом обтекании выпуклого тела показано на рис а-г. 20 Рис а Множественные ударные волны на профиле с ламинарным пограничным слоем показаны на рис.12.35а. Область местного сверхзвукового течения обычно заканчивается ударными волнами, на форму которых оказывает влияние пограничный слой. Если этот слой ламинарен, то сначала появляется последовательность так называемых ламбдообразных ударных волн.рис.12.35а На этой шлирен-фотографии, выполненной при вертикальной кромке ножа, видна картина обтекания двояковыпуклого профиля 12%-ной толщины при числе Маха невозмущенного потока, равном 0,8 (Н. W. Liepmann)

21 Ламбдообразные ударные волны над ламинарным пограничным слоем при возрастающих числах Маха представлены на рис б. Видно, что при возрастании местного числа Маха, значения которого приводятся под каждой парой шлирен-фотографий, ламбдообразные ударные волны сливаются между собой. Здесь ламинарный пограничный слой формируется на искривленной пластинке, которая создает растущую ограниченную область сверхзвукового течения, погруженную в дозвуковой поток. Тем самым имитируется сверхкритическое обтекание профиля. Символы на снимках слева указывают, что кромка ножа вертикальна для верхнего ряда, а для нижнего горизонтальна [Ackeret, Feldmann, Rott, 1946].рис б 21 Рис б

22 Рис в 22 Рис в Рис в иллюстрирует одиночную ламбдообразную ударную волну на профиле с ламинарным пограничным слоем. При увеличении числа Маха множественные ударные волны, показанные на предыдущем снимках, сливаются в одну волну. Здесь видно (кромка ножа горизонтальна), что при М = 0,9 ламинарный пограничный слой оторвался от поверхности двояковыпуклого профиля перед ударной волной и стал турбулентным (Н. W. Liepmann).

23 На рис г показан шар в свободном полете в воздухе при числе Рейнольдса , который схвачен вспышкой при теневом фотографировании почти за пределами поля зрения. Видно, что λ- скачок уплотнения вызывает турбулентный отрыв пограничного слоя перед экватором. При малых числах Маха отрыв слоя происходит за экватором вниз по потоку в зависимости от числа М (A. Stilp, 1965).рис г 23 Рис г

24 Рассмотрим качественно истечение газа из сужающегося насадка в неподвижный газ. При достижении на выходе сопла некоторого давления, равного его критическому значению, скорость газа на срезе сопла достигает скорости звука c *. Однако, вне сопла давление может оказаться и больше и меньше этого критического значения. Рассмотрим случай, когда давление вне сопла больше критического давления, то есть Р 0 >Р *. Нарисуем качественно картину движения газа вне сопла Истечение газа из сужающегося насадка Будем полагать для простоты, что сопло плоское. Край сопла можно рассматривать как источник возмущения (рис ). Газ вне сопла неподвижен, газ выходящий из сопла имеет скорость, равную скорости звука на срезе сопла.рис

25 25 Тогда от края сопла должен отойти тангенциальный разрыв. Для того чтобы давление в выходящем газе повысилось от давления Р * до давления Р 0 >Р *, необходимо, чтобы от края сопла отошла косая ударная волна. Ударные волны, отошедшие от обеих кромок сопла, пересекутся по известному закону, а затем упадут на тангенциальный разрыв, от которого они отразятся в виде волн разрежения. В свою очередь, волны разрежения после пересечения отразятся от тангенциального разрыва в виде ударных волн и так далее. Газ, прошедший ударные волны, тормозится, и давление возрастает, а газ, прошедший волны разрежения, ускоряется, и давление в нем падает. Это чередование ударных волн и волн разрежения будет происходить до тех пор, пока в силу вязкости и теплопроводности кинетическая энергия движения не диссилируется настолько, что скорость движения газа станет дозвуковой.

26 26 Если давление газа вне сопла Р 0 меньше чем Р *, то есть Р 0

27 27 Режим движения газа в сопле Лаваля в существенной мере зависит от противодавления. Рассмотрим качественно, как будет происходить движение внутри сопла Лаваля по мере уменьшения противодавления. На рис представлен график изменения давления вдоль сопла в зависимости от противодавления.рис Сопло Лаваля в нерасчетном режиме

28 28 Если Р вых =Р 0, то очевидно движение газа внутри сопла отсутствует и давление будет одинаковым в любом сечении сопла (кривая 1). По мере уменьшения выходного давления, газ начнет двигаться вдоль сопла, и распределение давления вдоль сопла представится некоторой кривой 2. Очевидно, минимальное давление всегда будет в самом узком сечении сопла. При дальнейшем уменьшении Р вых в минимальном сечении скорость газа может достигнуть скорости звука в этом критическом сечении, а давление будет равно Р *. При этом возникнет альтернатива. Или дальше за критическим сечением газ будет двигаться дозвуковым образом, а давление вниз по потоку будет возрастать от Р * до Р вых (кривая 3), или газ будет двигаться далее сверхзвуковым образом, а давление будет продолжать падать до Р вых (кривая 4). Очевидно, что эта альтернатива разрешается величиной противодавления. Давление Р вых можно рассчитать следующим образом. В была получена формула, связывающая величину поперечного сечения сверхзвуковой части сопла Лаваля с числом Маха ( ):

29 Поэтому, зная S вых и S min, мы могли бы рассчитать число Маха в сечении S вых и, в частности, число Маха в выходном сечении М вых. Далее, используя уравнение Бернулли и уравнение для адиабатического движения можно легко получить формулу, аналогичную (12.2.3) и можно вычислить расчетное значение Р вых по формуле: Если выходное давление будет равно Р вых, рассчитанному по формуле (12.3) то газ будет двигаться за критическим сечением (S мин ) сверхзвуковым образом. Если давление на выходе из сопла будет меньше расчетного Р вых, то дальнейшее ускорение потока и падение давления будет происходить вне сопла.12.3 Если же давление на выходе сопла будет больше расчетного Р вых, то движение газа внутри сопла за критическим сечением не может быть везде сверхзвуковым. В этом случае образуется прямая ударная волна, переводящая в этом сечении давление до ударной волны в давление Р вых > Р вых. По мере увеличения давления от Р вых до Р вых, ударная волна проникает все глубже и глубже внутрь сопла (кривые 5,6,7).кривые 5,6,7 29 (12.3)

30 Пусть пластинка, расположенная под углом атаки α к оси симметрии, обтекается сверхзвуковым потоком. Газ, движущийся выше АО, воспримет пластинку как тупой угол. От точки О отойдет волна разрежения, в которой поток повернет вдоль пластинки, давление и температура уменьшатся, а скорость потока возрастет. Газ, движущийся ниже АО, воспримет пластинку как клин. Поэтому от носика клина отойдет косая ударная волна, давление и температура за которой возрастут, а скорость потока уменьшится. После пластинки оба потока должны объединиться в соответствии с законом неразрывности сплошной струи и продолжить движение в прежнем направлении вдоль АО с одинаковыми давлениями и скоростями. Поэтому от точки О отойдет косая ударная волна в верхнем полупространстве, за которой давление увеличится до давления в набегающем потоке, а в нижнем полупространстве- волна разрежения, за которой давление увеличится ОБТЕКАНИЕ ПЛАСТИНКИ СВЕРХЗВУКОВЫМ ПОТОКОМ ПОД МАЛЫМ УГЛОМ АТАКИ

31 Рассчитав ударную волну и волну разрежения, можно найти Р 2 и Р 2, а также найти полную силу F и составляющие - силу сопротивления F x и подъемную силу F y.полную силу F и составляющие 31 Рис.12.39

32 32 Выводы Изучено отражение ударной волны от стенки, пересечение двух ударных волн Рассмотрено взаимодействие ударной волны, и волны с тангенциальным разрывом и процесс возникновения ударной волны в пограничном. Изучено сопло Лаваля в нерасчетном режиме. Рассмотрено истечение газа из сужающегося насадка и обтекание пластинки сверхзвуковым потоком под малым углом атаки.

33 33 Информационное обеспечение лекции Литература по теме: Ландау Л.Д., Лившиц Е.М.. Гидродинамика. М.: Наука с. Седов Л.И. Механика сплошной среды. М.: Наука Т с.; Т.2, 568с. Фабрикант Н.Я. Аэродинамика. М.: ГИТТЛ с. Лойцянский Л.Г. Механика жидкости и газа. М.: Наука с.

34 34 Справочные данные Курс лекций является частью учебно-методического комплекса «Газовая динамика». Автор: Породнов Борис Трифонович, д. ф. – м. н., профессор кафедры молекулярной физики УГТУ-УПИ. Учебно-методический комплекс подготовлен на кафедре МФ ФТФ ГОУ ВПО УГТУ-УПИ.