CC BY
18
УДК 533.6.011
В. Н. Голубкин1'2, Г. Б. Сизых1
1 Московский физико-технический институт (государственный университет) 2Центральный аэрогидродинамический институт имени профессора Н. Е. Жуковского
Ои и и
некоторых свойствах газовых течении с осевой
симметрией
Исследованы осесимметричные течения идеального совершенного газа с закруткой в окружном направлении. Рассмотренные течения могут быть небаротропными и вихревыми. На основе анализа уравнений Эйлера получен ряд неизвестных ранее качественных свойств таких течений. В частности, доказано, что при наличии окружной скорости осесимметричное течение не может быть потенциальным, если хотя бы часть оси симметрии находится внутри течения (а не внутри обтекаемого тела). Также обнаружена эквивалентность двух типов течений: закрученных течений в отсутствие внешних сил и незакрученных течений в потенциальном осесимметричном поле внешних сил.
Ключевые слова: уравнения Эйлера, закрученные осесимметричные течения газа, завихренность, циркуляция, преобразование Мунка-Прима.
V. N. Golubkin1'2, G. B. Sizykh1
1 Moscow Institute of Physics and Technology (State University) 2Central Aerohydrodynamic Institute n.a. professor N. E. Zhukovskiy
Some properties of gas flows with axial symmetry
The axisymmetric flow of an ideal gas with a twist in the circumferential direction is investigated. Streams can be nonbarotropic and vortex ones. Based on the analysis of the Euler equations, a number of previously unknown qualitative properties of such flows are obtained. In particular, it is proved that (in the presence of the circumferential speed) the gas flow cannot be potential if at least a portion of the axis of symmetry is situated inside of the flow (not inside of a streamlined body). Also the equivalence of the two types of flows is discovered: swirling flows in the absence of external forces and not swirling flows in an axisymmetric potential field of external forces.
Key words: Euler's equations, axisymmetric swirling gas flow, vorticity, circulation, Munk and Prim transformation.
1. Введение
Закрученные осесимметричные течения представляют собой важный частный случай пространственных течений, описывающихся двумерными уравнениями. С одной стороны, в них в полной мере проявляется сложная связь радиально-осевого и окружного движений. С другой стороны, осевая симметрия накладывает дополнительные связи на параметры движения, которые упрощают исследование по сравнению с общим трехмерным случаем. Реальные конструкции редко бывают полностью осесимметричными. Но проверка закономерностей осесимметричных потоков может быть полезна для верификации численных методов расчета трехмерных схем, когда такие схемы используются в осесимметричном случае. Поэтому исследование закрученных осесимметричных течений представляет не только теоретический, но и практический интерес.
В дальнейшем все рассматриваемые функции считаются достаточно гладкими, насколько это требуется проводимыми математическими операциями и рассуждениями. Это относится либо целиком ко всему течению, если в нем нет скачков уплотнения, тангенциальных или слабых разрывов, либо к зонам течения, расположенным между скачками и разрывами. Для краткости будем называть такие зоны зонами гладкости.
2. Основные обозначения и уравнения движения
Рассмотрим осесимметричное течение идеального газа в потенциальном поле внешних сил. Обозначим: V - вектор скорости, р - плотность, р - давление, П - потенциал массовых сил, О = rotV - завихренность. Движение газа описывается уравнениями Эйлера, которые можно записать в форме Громеки-Ламба [1]:
дV У'р (V2
- + П х V = -- --( т + П
(V+") ,
|| + div(pV) = 0.
(1)
(2)
Введем цилиндрическую систему координат г,<р,г с началом в точке О и осью Ох, совпадающей с осью симметрии. Пусть ег,е^,ег - правая тройка единичных векторов в радиальном, окружном и осевом направлениях соответственно. Вектор скорости имеет вид суммы: V = Угег + У^е^ + Угех = Vr + + Vх, где его компоненты по осям Уг, У^ и Уг, как и другие газодинамические функции, зависят только от координат г, г и от времени ¿. Очевидно, в этом случае градиент в правой части уравнения (1) имеет нулевую ^-компоненту. Закрученное осесимметричное течение имеет отличную от нуля окружную скорость Уу. Величину 7 = 2ттгУ^ назовем окружной циркуляцией.
В силу осевой симметрии ^^ = 0^, меридиональная составляющая завихренности Огх, равная векторной сумме радиальной и осевой составляющих (Ог-г = Ог + Ох), является ротором окружной скорости У^. Аналогично, окружная составляющая завихренности О^ является ротором меридиональной скорости V,- + Vх:
Ог
^ V
О,,
Учитывая, что
О^ х Vч>
0,
Огг X V
rot (Vr + Vz).
-- (?) - (?)
получим
О х V = Огг х (V, + V) + О^ х (V, + V) - - - ег.
(?) -(?)
2
Первое слагаемое в правой части направлено вдоль е^, а последние три лежат в меридиональной плоскости, проходящей через ось симметрии. Поэтому уравнение (1) равносильно системе двух уравнений:
9 (V, + V,) дИ
+ Од
х (у, + V,) = --(£ + +П
Р
2
2
И
ОV
сЛ
+ Огг х (Vr + Vz) = 0.
Уравнение неразрывности (2) в цилиндрической системе координат имеет вид
др д , тгЧ рУг д , ^ + тт (РУг) + + 7Г №) = 0.
от ог г ох
(3)
(4)
(5)
Уравнения (3) - (5) являются основными уравнениями, описывающими закрученные осесимметричные течения сжимаемого газа. Для замыкания их нужно дополнить еще одной связью, например выражающей адиабатичность течения, то есть сохранение в движущейся частице энтропийной функции ^, где к - показатель адиабаты.
Иногда течение считают потенциальным, то есть безвихревым: О = 0. Иногда - баро-тропным: р = /(р). И так далее. Но некоторые общие свойства осесимметричных течений следуют из уравнений (3) - (5) без каких-либо дополнительных предположений. Поэтому ниже, если не оговорено иного, рассматривается самый общий случай, когда течение считается неадиабатичным, небаротропным и завихренным.
е
г
е
г
3. Закрученность и завихренность
Рассмотрим течение в зоне гладкости (см. выше). Если зона гладкости содержит (хотя бы частично) ось симметрии, то на оси, в силу непрерывности, окружная компонента скорости У^ должна обращаться в нуль. Воспользуемся этим обстоятельством.
Легко проверить, что меридиональная составляющая завихренности: = Пг + и окружная циркуляция: 7 = 2пгУ^ связаны между собой соотношением
Пгг = —[е^ х У7] 2 кг
(6)
При отсутствии завихренности равна нулю любая ее составляющая, в частности: = 0. Поэтому из (6) следует, что в безвихревом течении окружная циркуляция 7 = 2пгУ^ постоянна по пространству.
Если ось симметрии (на которой У^ = 0) расположена внутри поля течения, то и постоянная по пространству величина окружной циркуляции 7 = 2пгУ^ должна всюду обращаться в нуль. Поэтому окружная скорость У^ равна нулю не только на оси, но и во всем поле течения. А это значит, что течение является незакрученным. Таким образом, приходим к следующему выводу.
Если осесимметричное течение газа потенциально и ось симметрии находится в поле течения, то такое течение лишено закрутки.
Заметим, что при осесимметричном обтекании тел вращения конечной длины можно указать часть оси симметрии, находящуюся внутри течения. Поэтому обтекание таких тел безвихревым потоком обеспечивает отсутствие закрутки. Наличие скачков уплотнения не меняет ситуацию. Действительно, в осесимметричных течениях поверхности скачков представляют собой поверхности вращения. Любая частица газа движется в составе осесим-метричного кольца, являющегося фигурой вращения этой частицы вокруг оси симметрии Ох. Момент импульса этого кольца сохраняется при переходе через скачок уплотнения. Поэтому после прохождения частицей скачка уплотнения окружная циркуляция не изменится и останется равной нулю.
Итак, внешнее безвихревое закрученное обтекание тела вращения возможно, только если во всей рассматриваемой области течения ось симметрии находится внутри обтекаемого тела.
В качестве примера рассмотрим течения в осесимметричных каналах. На рис. 1 схематически изображено течение в осесимметричной области, ограниченной стенками канала, входным (А — А) и выходным (С — С') сечениями, а также поверхностью обтекаемого тела. Внутренность тела изображена серым цветом. Осесимметрично расположенные точки обозначены одной и той же буквой (например, точки А).
Рис. 1. Часть оси симметрии находится внутри течения
Часть оси симметрии, расположенная левее точки В, находится внутри течения, а не внутри обтекаемого тела. Поэтому при любых параметрах течения во входном и выходном сечениях и при наличии закрутки течение будет вихревым.
Подобный вывод нельзя сделать для течения в области, изображенной на рис. 2 (несмотря на то, что рис. 2 можно считать частью рис. 1). В этом случае ось симметрии целиком лежит внутри обтекаемого тела.
Рис. 2. Ось симметрии находится внутри обтекаемого тела
В изображенной на рис. 2 части осесимметричного канала ни одна точка течения не лежит на оси симметрии. Следовательно, полученный в данной статье результат не исключает возможности задать такие параметры течения во входном (А — А') и выходном (С — С') сечениях, которые могут обеспечить реализацию безвихревого обтекания с закруткой. Например, если обтекаемое тело представляет собой цилиндр радиуса г0 > 0, то одним из возможных течений является так называемый потенциальный вихрь. Это стационарное течение вне цилиндра ( г > г0), компоненты скорости которого Vr = 0, Vv = С\/г,VZ = С2, где Ci, С2 = const.
Следует подчеркнуть, что полученное выше необходимое условие существования закрутки в потенциальном течении (ось симметрии должна быть внутри тела) не является достаточным. Как и другие необходимые, но недостаточные условия, оно может эффективно использоваться только в случаях, когда это необходимое условие нарушено. То есть для утверждения о невозможности существования закрутки в потенциальном течении. Выполнение этого необходимого условия (как, например, в ситуации, изображенной на рис. 2) еще не означает, что потенциальное течение с закруткой обязательно может быть реализовано. Обсуждение возможности реализации потенциального течения с закруткой при выполнении полученного выше необходимого условия выходит за рамки данного исследования.
4. Эквивалентность стационарных тотально закрученных
и незакрученных течений в потенциальном поле массовых сил
В работе [2] без предположений о баротропности или адиабатичности показано, что для осесимметричных течений идеального газа окружная циркуляция 7 = 2nrV^ сохраняется в частице газа, то есть является первым интегралом движения:
ж + (VV)J = °
В частности, для стационарных течений это означает, что окружная циркуляция 7 = 2nrVv может быть различной на разных линиях тока, но сохраняет постоянное значение на каждой линии тока.
Тотально закрученным течением назовем течение, в котором нет точек с нулевой окружной скоростью Vy. Рассмотрим стационарное тотально закрученное течение в отсутствие внешних сил. В стационарном случае окружная циркуляция 7 постоянна вдоль линий тока (V ■ V) 7 = 0 и не обращается в нуль. Поэтому ее можно использовать в качестве функции преобразования Мунка-Прима [3], связывающего параметры течений газа с одинаковой
картиной линий тока и идентичным распределением давления. Эта связь выражается формулами
V = 7У1, р= , р = Р1. (7)
72
Непосредственной проверкой можно убедиться, что функции нового течения (с индексом 1) удовлетворяют стационарным уравнениям Эйлера:
ё1у(р^) = 0, «1 х VI = - — , «1 = го^ь
Р1 V 2 /
Обозначим —1 = У1гег + У1(е( + У1гех = —1г + —1( + . Замена (7) подобрана таким образом, что в новом течении окружная циркуляция постоянна:
У1( = и 71 = 2шУ1( = 1.
2кг
Поэтому, согласно (6)
«1гг = -77— [е( х У71] = 0, 2к г
и, следовательно, в силу стационарности, уравнение (4) вырождается. Оставшиеся уравнения (3), (5) для нового течения ( У1г, У1(, У1г, р1, р1) в отсутствие массовых сил (П = 0), с учетом равенства У1( = ^, выглядят следующим образом:
П„ х (—,,+ V. ) = -V (* + Я) + 1 (£ .
О (Р1У1Г) + ^ + О- (Р1У1,) = 0. ог г ох
Эти уравнения идентичны уравнениям незакрученного течения с параметрами
У2г = У1г, У2( = 0, У2х = У1г, р2 = р1, р2 = р1 в поле массовых сил с потенциалом П = :
«2( х (—2, + —2,) = -VI2 -V (3 + % + ) . (8)
О (Р2У2Г) + ^ + ОО- (/»У*.) = 0. (9)
о о
Таким образом, каждому решению ( У2г, У2( = 0, У2х, р2, р2) какой-либо двумерной задачи незакрученного осесимметричного обтекания с потенциалом внешних массовых сил (уравнения (8) и (9)) соответствует семейство решений ( Уг, У(,Уг, р, р) задач закрученного обтекания в отсутствие внешних массовых сил:
Уг =7У2г, Уг =7У2г, У( = ^, Р = ^, Р = Р2, (10)
2кг 72
где 7 = 0 - произвольная функция постоянная вдоль линий тока. При этом каждая конкретная функция 7 задает такое решение из указанного семейства, в котором величина 7 является окружной циркуляцией 2ттг
Только что установленная связь удобна тем, что поверхности вращения линий тока незакрученных и тотально закрученных течений, связанных отношениями (10), совпадают. Поэтому все эти течения соответствуют обтеканию одних и тех же конфигураций. Любое тотально закрученное течение ( У( = 0) внутри выбранной пространственной конфигурации может быть получено соответствующим подбором функции 7, постоянной вдоль линий тока, по формулам (10) из незакрученного решения задачи с уравнениями (8), (9) и с «незакрученными» граничными условиями.
5. Инвариант линий тока стационарных закрученных осесимметричных течений
При отсутствии радиально-осевого движения, когда все частицы вращаются вокруг оси Ог, инвариантом линии тока является любой параметр жидкой частицы. Поэтому содержательным является исследование инвариантов при наличии радиально-осевой скорости: V,- + V.; = 0.
Для стационарных течений баротропного газа справедлива теорема Эйлера-Эртеля [4], которая позволяет получить инвариант линий тока с использованием другого известного инварианта. Так, например, в работе [5] с помощью теоремы Эйлера-Эртеля был получен новый инвариант течений газа. Но в рассматриваемом течении применение теоремы Эйлера-Эртеля к энтропии, к полной энтальпии и к окружной циркуляции дает тривиальный результат. В осесимметричных течениях место инвариантов, получаемых с помощью теоремы Эйлера-Эртеля, занимает другой, нетривиальный инвариант. При этом в отличие от инвариантов Эйлера-Эртеля не предполагается баротропность или адиабатичность течения.
В стационарном случае при условии V,.+Vz = 0 из формулы (4) следует существование скалярного поля Л, такого, что
А (V, + V,) = ОГ2. Применение оператора дивергенции к этому уравнению дает
Л div (V,- + Vz) + (V + Vz) -Л = 0.
С другой стороны, из уравнения неразрывности следует, что
pdiv ^ + V) + (V + V) - р = 0.
Следовательно,
или
(V, + V) -А = 0, Р
^ + ^ >- ж+щ = а (11)
Поэтому отношение р^+у' | постоянно вдоль линий тока.
Подчеркнем, что формула (11) получена без предположений о баротропности или адиа-батичности течения. Более того, при ее выводе не использовалось уравнение импульсов (1), а только уравнение неразрывности (2). Поэтому формула (11) остается верной и для стационарных осесимметричных течений вязкого газа.
Заметим, что в винтовых течениях (в которых завихренность параллельна скорости) величина ^ | также является инвариантом вследствие инвариантности для винтовых
течений величины Л^у [4]. Действительно, поскольку в винтовых течениях завихренность параллельна скорости, отношение длины вектора О к длине вектора V равно отношению длин проекций этих векторов на меридиональную плоскость, то есть - отношению |Огг | к IV + Vх|. При этом формула (11) верна как для винтовых течений, так и для течений общего вида. Поэтому она представляет собой самостоятельный новый результат наряду с уже известными для осесимметричных течений вязкой жидкости [6, 7].
6. Заключение
Проведен анализ полных уравнений Эйлера осесимметричного закрученного течения идеального газа, которое в общем случае может быть вихревым, а энтропия может быть различной на разных линиях тока. Получен ряд новых, неизвестных ранее свойств таких течений.
Результаты могут быть применены для качественного анализа течений, а полученный новый инвариант (11) - для верификации численных расчетов при отладке расчетных комплексов.
Литература
1. Лойцянский Л.Г. Механика жидкости и газа. М.: Дрофа, 2003.
2. Голубкин В.Н., Марков В.В., Сизых Г.Б. Интегральный инвариант уравнений движения вязкого газа // ПММ. 2015. Т. 79, вып. 6. С. 808-816.
3. Munk M., Prim R. On the multiplicity of steady gas flows having the same streamline pattern // Proc. Nat. Acad. Sci. USA. 1947. V. 33. P. 137-141.
4. Truesdell C. The kinematics of vorticity. Bloomington: Indiana University Press, 1954.
5. Голубинский А.И., Голубкин В.Н. О некоторых свойствах сохранения в газовой динамике // ПММ. 1985. Т. 49, вып. 1. С. 115-119.
6. Брутян М.А., Голубкин В.Н., Крапивский П.Л. Об уравнении Бернулли для осесим-метричных течений вязкой жидкости // Ученые записки ЦАГИ. 1988. Т. 19, № 2. С. 98-100.
7. Сизых Г.Б. Эволюция завихренности в закрученных осесимметричных течениях вязкой несжимаемой жидкости // Ученые записки ЦАГИ. 2015. Т. 46, № 3. С. 14-20.
References
1. Loytsyansky L.G. Mechanics of Liquids and Gases. Moscow: Drofa, 2003. (in Russian).
2. Golubkin V.N., Markov V.V., Sizykh G.B. Integral invariant of the equations of motion of viscous gas. Journal of Applied Mathematics and Mechanics. 2015. V. 79, I. 6. P. 808-816. (in Russian).
3. Munk M., Prim R. On the multiplicity of steady gas flows having the same streamline pattern. Proc. Nat. Acad. Sci. USA. 1947. V. 33. P. 137-141.
4. Truesdell C. The kinematics of vorticity. Indiana university publications. Science series N 19. Indiana university, Bloomington, Indiana. 1954. 232 p.
5. Golubinskii A.I., Golubkin V.N. On certain conservation properties in gas dynamics.PMM. 1985. V. 49, I. 1. P. 115-119. (in Russian).
6. Brutyan M.A., Golubkin V.N., Krapivskiy P.L. About the Bernoulli's equation for axisymmetric flows of a viscous fluid. TsAGI Science Journal. 1988. V. 19, N 2. P. 98-100. (in Russian).
7. Sizykh G.B. Evolution of vorticity in swirling axisymmetric flows of a viscous incompressible fluid. TsAGI Science Journal. 2015. V. 46, I. 3. P. 14-20. (in Russian).
Поступила в редакцию 03.03.2017
