CC BY
12
УДК 533.6.011
В. Н. Голубкин1'2, Г. Б. Сизых1
1 Московский физико-технический институт (государственный университет) 2Центральный аэрогидродинамический институт имени профессора Н. Е. Жуковского
Точки минимума давления в осесимметричных течениях жидкости
Показано, что при отсутствии точек торможения в ограниченной области осесим-метричного вихревого течения идеальной несжимаемой жидкости давление достигает минимального значения на границе области. Нарушение этого свойства давления свидетельствует о наличии внутренней точки торможения. Полученный вывод может быть использован для верификации численных расчетов течения при обтекании тел вращения.
Ключевые слова: уравнения Эйлера, вихревые течения, принцип минимума давления, осесимметричные течения, несжимаемая жидкость.
V. N. Golubkin1'2, G. B. Sizykh1
1 Moscow Institute of Physics and Technology (State University) 2 Central Aerohydrodynamic Institute n.a. professor N. E. Zhukovskiy
Points of minimum pressure in the axisymmetric fluid
flows
It is shown that in the absence of braking points in the bounded region of the axisymmetric vortex flow of an ideal incompressible fluid pressure reaches a minimum value in the boundary region. The violation of this property of pressure indicates the presence of braking points. The result can be used to verify numerical calculations of the flow for streamlined bodies.
Key words: Euler equations, vortex flow, minimum principle for the pressure, axisymmetric flow, incompressible fluid.
1. Введение
Давление в стационарных безвихревых течениях идеальной несжимаемой жидкости не может достигать минимума во внутренней точке течения. Этот факт следует из теории гармонических функций. Для вихревых течений идеальной жидкости известен принцип Роуланда [1], предлагающий различные утверждения о минимуме или о максимуме давления, в зависимости от знака некоторой функции Q, для определения которой нужно знать первые производные компонент скорости (см. ниже). Позже этот принцип был обобщен Гамелем [2] на случай вязкой несжимаемой жидкости. В свою очередь, принцип Гамеля был обобщен Трусделлом [3] (см. также [4, §28]) на случай сжимаемой жидкости. Однако оба эти обобщения не дают новых утверждений для идеальной несжимаемой жидкости, поскольку в этом случае они сводятся к изначальному принципу Роуланда, который практически непригоден для априорных выводов, поскольку знак функции Q трудно указать заранее.
В работе [5] обнаружено свойство экстремальных значений давления в плоском вихревом течении идеальной жидкости: если в ограниченной области течения жидкости отсутствуют точки торможения и давление р не постоянно всюду, то и максимум и минимум р достигаются только на границе области.
В данной статье обоснован принцип минимума давления в осесимметричных вихревых течениях идеальной жидкости, условия которого совпадают с условиями принципа [5] и поэтому имеют более обозримый и идентифицируемый вид, чем условия Роуланда.
Полученный принцип минимума можно применить к практически важным осесиммет-ричным течениям.
2. Принцип минимума давления в осесимметричном течении идеальной жидкости
Рассмотрим пространственное стационарное течение идеальной несжимаемой жидкости в замыкании G некоторой ограниченной области G. Такое течение описывается уравнениями Эйлера:
О х V = -1 vp - 1W2, (1)
р 2
div V = 0, (2)
где p - давление, V - вектор скорости, р = const - плотность жидкости, О = rot V -завихренность, V = |V|, V - оператор Гамильтона.
Будем считать, что в рассматриваемой области отсутствуют точки торможения, то есть
V > 0. (3)
При этом допускается наличие точек торможения (V = 0) на границе области.
Также будем считать, что все гидродинамические функции являются дважды непрерывно дифференцируемыми функциями координат в области G и непрерывными на ее замыкании G.
Покажем, что для осесимметричных течений при выполнении указанных условий справедлив принцип минимума давления: если давление не постоянно в области G, то оно достигает минимального значения только на границе этой области.
Для доказательства введем трехмерную прямоугольную декартову систему координат Oxyz и обозначим через u, v, w соответствующие компоненты вектора скорости. Вычисляя дивергенцию обеих частей векторного уравнения (1), получим
1Др = Q2 - V-rotO - 1AV2, (4)
р 2 ' w
где Д - оператор Лапласа.
Учитывая, что 2AV2 = V-AV + (Vu)2 + (Vv )2 + (Vw)2 и что для несжимаемой жидкости V-rot О = -V-AV, из (4) после перестановки слагаемых получим
-Ap = Q2 - (Vu)2 - (Vv)2 - (Vw)2. (5)
р
Замечание. Уравнение (5) известно давно [1]. Функция Q, о которой говорилось во введении, есть половина правой части (5). Если во всей области Q ^ 0 или Q ^ 0, то уравнение (5) становится супер- или субгармоническим, и давление (если не постоянно всюду) не может достигать во внутренней точке течения соответственно минимума или максимума. В этом состоит принцип Роуланда. Однако до сих пор знак параметра Q можно было априорно указать только для безвихревых течений, в которых Q ^ 0. В данной статье предложен способ исследования экстремальных свойств давления для незакручен-ных осесимметричных течений без выяснения знака параметра Q. Для этого уравнение (5) преобразовано к другому виду, допускающему применение теоремы Хопфа [6, 7].
Уравнение (5) верно для общего SD-случая. Далее рассмотрим незакрученные осе-симметричные течения и воспользуемся присущим им свойством перпендикулярности
векторов скорости и завихренности. Тогда из (1) следует, что (Q V)2 = Qvp + 2VV2 j ,
или, учитывая принятое выше предположение (3) об отсутствии точек торможения (V = 0), п 2 = (vp)2 + PV2 vp •vV2 + V? ( 1vV2)2 =
= (pé?vp + pV?vV2) •vp + V?(2vV2)2.
Подставляя это выражение для квадрата завихренности в (5), получим
p ap - (ppV? vp + pV? vV2) •vp =
= V*( 1 vV2)2 - (vu)2 - (vv )2 - (vw)2. (6)
Левая часть этого уравнения инвариантна относительно выбора системы координат Oxyz. Следовательно, правая часть
I = 2 vV2)2 - (vu)2 - (vv )2 - (vw)2
также инвариантна, то есть ее значение не зависит от выбора системы координат Оxyz. Возьмем произвольную точку области G и расположим систему координат о xyz так, чтобы положительное направление оси Оx совпало с направлением скорости V. Тогда в выбранной точке соответствующие компоненты скорости будут z = V > 0, z = 0, wz = 0. Следовательно, в этой точке
I = â2(zvu + zvz + wvw)2 - (vu)2 - (vz)2 - (vw)2 = -(vz)2 - (vw)2 < 0. u2
Поскольку точка была выбрана произвольно, неравенство I ^ 0 выполнено во всей области G. Поэтому во всех точках области G уравнение (6) можно записать в виде
Ap -{Vv + hvV2} -VP = f
где правая часть неположительна: f = pI ^ 0.
Вернемся к выбранной с самого начала системе координат Oxyz. Если записать левую часть последнего уравнения в координатной форме, то второе слагаемое примет вид: b1 Ш + b2 щ + b3 dp, где коэффициенты b\, b2 и b3 суть компоненты вектора, записанного в фигурных скобках, и зависят не только от гидродинамических функций u, v, w, p, но и от их первых производных. Таким образом, уравнение (6), которое является следствием уравнений Эйлера, может быть представлено в виде
д 2p д 2p d2p d2p д 2p d2p
+ 2а\2~Л 7]--+ a22T\ 2 + 2a13^7]--+ 2a23^~7,--+ +
dx2 rJxrjy dy2 rJxrJz rjyrjz dz2
+bl dx+b2 dp+h dp = f, (7)
dx dy dz
где коэффициенты a11 = a22 = = 1, a12 = a13 = a23 = 0 и функция f ^ 0.
Для анализа свойств решения уравнения (7) (в предположении, что оно существует) воспользуемся одним из следствий теоремы Хопфа [6, 7], приведенным в работе [8]. Это следствие специально сформулировано для использования в аэрогидродинамике. Сформулируем его с использованием обозначений данной статьи.
Пусть во всех точках ограниченной области G коэффициенты a11, a12, a22, a13, a23 и a33 уравнения (7) являются коэффициентами положительно определенной квадратичной формы A. Пусть, далее, для любой точки M(x,y,z) g G существуют 'числа ш = w(x,y,z) > 0 и a = a(x,y,z) > 0 такие, что замкнутый круг U(M; ш) целиком
лежит в области G ив нем все коэффициенты уравнения (7) ограничены и выполняется неравенство det A > а. И пусть во всех точках области G выполняется неравенство f < 0.
Тогда если решение p £ C2(G) уравнения (7) непрерывно в замкнутой области G, то p ^ min p во всей области G. При этом, если p = const в G, равенство p = min p возможно
Г dG Г , г l ao
только на dG.
Если учесть, что согласно (3) скорость не обращается в нуль во внутренних точках области G, несложно проверить, что уравнение (7) удовлетворяет всем условиям приведенного следствия теоремы Хопфа. Отсюда следует, что при отсутствии точек торможения в ограниченной области осесимметричного вихревого течения идеальной несжимаемой жидкости давление достигает минимального значения только на границе области.
3. Случай наличия точек торможения
Проиллюстрируем существенность требования отсутствия точек торможения для полученного принципа минимума. Приведем пример, в котором нарушение этого требования приводит к «появлению» внутренней точки минимума давления.
Введем цилиндрическую систему координат r,p,z с началом в точке O и будем рассматривать симметричные относительно оси Oz течения. Пусть er, e^, ez - правая тройка единичных векторов в радиальном, окружном и осевом направлениях соответственно. Вектор скорости незакрученного течения имеет вид V = Vr er + Vz ez, где компоненты скорости Vr и Vz зависят только от переменных r, z.
В данном разделе все величины r, z, Vr, Vz и давление p, отнесенное к плотности, будем считать безразмерными. (Размерные величины координат V, V, компонент скорости vr, Vz и давления p можно выразить через безразмерные переменные по формулам: V = rL, V = zL, Vr = Vr U, Vz = Vz U и p = ppU2, где L и U - характерные значения длины и скорости в рассматриваемом течении.)
Будем искать точное решение уравнений (1) и (2) в виде
Vr = -(g(r) sin z)/r, Vz = — (g'(r) cos z) ¡r, (8)
p = 2 — (g(r) cos z)2 — 0.5 (Vr2 + Vz2) , (9)
где g = g(r) - некоторая функция одной координаты r. Если записать уравнения (1) и (2) в цилиндрической системе координат r,p,z и подставить в них выражения (8) и (9), то эти уравнения окажутся выполненными, если функция g удовлетворяет обыкновенному дифференциальному уравнению (здесь и далее штрих означает дифференцирование по аргументу r):
g" — r — g + 2r2g = 0. (10)
Рассмотрим задачу Коши для уравнения (10) на отрезке 0.75 ^ r ^ 1.25 с условиями, заданными в середине отрезка:
g(1) = 1,g'(1) = 0. (11)
Эта задача состоит из двух «стандартных» задач Коши. Одна из них - для уравнения (10) на отрезке 0.75 ^ r ^ 1 с условиями (11) на правом конце отрезка. Другая - для уравнения (10) на отрезке 1 ^ r ^ 1.25 с условиями (11) на левом конце отрезка.
Коэффициенты уравнения (10) ограничены на отрезке 0.75 ^ r ^ 1.25. Поэтому решение задачи Коши (10), (11) на отрезке 0.75 ^ r ^ 1.25 существует и единственно. Формулы (8) и (9), в которых g(r) является решением этой задачи, дают точное решение уравнений (1) и (2) в полосе 0.75 < r < 1.25, —те < z < Рассмотрим прямоугольную область 0.75 < r < 1.25, —0.25 < z < 0.25. В выражения (8) и (9) входят функция g(r) и ее производная g'(r). Их графики, полученные в результате численного решения задачи (10), (11), представлены на рис. 1.
Рис. 1
Из (8) и (11) следует, что точка 5 с координатами г = 1, г = 0 является точкой торможения (Уг = Уг = 0). Линии тока этого течения показаны на рис. 2.
0.25 0.00 0.25 2
Рис. 2. Линии тока
Жидкость вращается вокруг точки 5, и кривизна линий тока показывает, что давление должно уменьшаться по мере приближения к точке 5. Поэтому в точке 5 будет минимум давления. Результаты расчетов в виде картины изобар с оттенками серого цвета приведены на рис. 3.
■0.25 0.0 О 0.25 2
Рис. 3. Картина изобар
Расчеты показывают, что в рассматриваемой области точкой минимума давления является внутренняя точка 5. Таким образом, при наличии внутренней точки торможения, когда условия установленного принципа минимума нарушаются, минимум давления может располагаться внутри области, а не на границе.
Численные расчеты проведены с высокой точностью. Тем не менее предложим строгое доказательство того, что точка S является точкой локального минимума давления. Из (8) и (9) следует, что g"(1) = —1. Это позволяет найти значения первых и вторых производных p (r, z) в точке z = 0, r = 1. В итоге в окрестности точки S имеем разложение
p (1 + Ar, Az) = p (1,0) + (Ar)2 + (Az)2 + o ((Ar)2 + (Az)2) , Ar ^ 0, Az ^ 0,
показывающее, что точка S действительно является точкой локального минимума давления.
4. Точки зарождения кавитации
В монографиях [9, 10] объясняется, почему кавитация при медленном увеличении скорости набегающего потока зарождается именно на поверхности обтекаемых тел. Объяснение основано на том, что точки минимума давления в безвихревых течениях находятся на границе. Подчеркнем, что эти выводы обоснованы в [9, 10] только для безвихревых (потенциальных) течений. Принцип работы [5] и принцип минимума давления, полученный в данной статье, позволяют сделать аналогичные выводы для плоских и осесимметричных вихревых течений. Кавитация в таких вихревых течениях должна зарождаться в точках минимума давления, которые расположены именно на поверхности обтекаемых тел.
5. Заключение
Обоснован принцип минимума давления для осесимметричных вихревых течений несжимаемой жидкости. Его условия имеют более обозримый вид, чем это было приведено в [1], и поэтому данный принцип можно использовать для априорных оценок, для качественного анализа течений и для верификации численных методов расчета наряду с полученными ранее принципами экстремальных значений [5, 11].
Литература
1. Rowland H. On the motion of a perfect incompressible fluid when no bodies are present // American Journal of Mathematics. 1880. V. 3. P. 226-268.
2. Hamel G. Ein allgemeiner Satz uber den Druck bei der Bewegung volumbestandiger Flussigkeiten // Monatshefte Math. Phys. 1936. V. 43. P. 345-363.
3. Truesdell C. Two Measures of Vorticity // J Rational Mech. Anal. 1953. V. 2. P. 173-217.
4. Серрин Дж. Математические основы классической механики жидкости. М.: Изд-во ИЛ, 1963.
5. Голубкин В.Н., Ковалёв В.П., Сизых Г.Б. Принцип максимума давления в плоских течениях идеальной несжимаемой жидкости // Ученые записки ЦАГИ. 2016. Т. 47. № 6. С. 28-36.
6. Hopf E. Elementare Bemerkungen uber die Losungen partieller Differentialgleichungen zweiter Ordnung vom Elliptischen Typus // Sitzungsberichte der Preussischen Akademie der Wissenschaften. 1927. V. 19. P. 147-152.
7. Миранда К. Уравнения с частными производными эллиптического типа. М.: Издательство иностранной литературы, 1957.
8. Беспорточный А.И., Бурмистров А.Н., Сизых Г.Б. Вариант теоремы Хопфа // Труды МФТИ. 2016. Т. 8, № 1. С. 115-122.
9. Седов Л.И. Механика сплошной среды. Т. 2. М.: Наука, 1973.
10. Бэтчелор Дж. Введение в динамику жидкости. М.: Мир, 1973.
11. Голубкин В.Н., Сизых Г.Б. Принцип максимума функции Бернулли // Ученые записки ЦАГИ. 2015. Т. 46, № 5. С. 53-56.
References
1. Rowland H. On the motion of a perfect incompressible fluid when no bodies are present. American Journal of Mathematics. 1880. V. 3. P. 226-268.
2. Hamel G. Ein allgemeiner Satz uber den Druck bei der Bewegung volumbestandiger Flüssigkeiten. Monatshefte Math. Phys. 1936. V. 43. P. 345-363.
3. Truesdell C. Two Measures of Vorticity. J. Rational Mech. Anal. 1953. V. 2. P. 173-217.
4. Serrin J. Mathematical principles of classical fluid mechanics. Moscow: Publ. IL, 1963. (in Russian).
5. Golubkin V.N., Kovalev V.P., Sizykh G.B. Maximum principle for pressure in ideal incompressible fluid flows. TsAGI Science Journal. 2016. V. 47, I. 6. P. 28-36. (in Russian).
6. Hopf E. Elementary remarks on the solutions of partial differential equations of second order of Elliptic type. Meeting reports of the Prussian Academy of Sciences. 1927. V. 19. P. 147-152. (in German).
7. Miranda C. Partial Differential Equations of Elliptic Type. Moscow: Publ. IL, 1957. (in Russian).
8. Besportochnyy A.I., Burmistrov A.N., Sizykh G.B. Version of the Hopf theorem. Proceedings of MIPT. 2016. V. 8, N 1. P. 115-122. (in Russian).
9. Sedov L.I. A cource in continuum mechanics. V. 2. Moscow: Nauka, 1973. (in Russian).
10. Batchelor G.K. An Introduction to Fluid Dynamics. Moscow: Mir, 1973. (in Russian).
11. Golubkin V.N., Sizykh G.B. Maximum principle for Bernoulli function. TsAGI Science Journal. 2015. V. 46, N 5. P. 53-56. (in Russian).
Поступила в редакцию 05.06.2017
