Дозвуковой принцип максимума для неизоэнтропийных течений Текст научной статьи по специальности «Математика»

Ovil Aviation High Technologies

Vol. 20, No. 02, 2017

УДК 533.6.011

ДОЗВУКОВОЙ ПРИНЦИП МАКСИМУМА ДЛЯ НЕИЗОЭНТРОПИЙНЫХ ТЕЧЕНИЙ

Г.Б. СИЗЫХ1

1 Московский физико-технический институт (государственный университет), г. Москва, Россия

Дозвуковой принцип максимума справедлив для дозвуковых стационарных безвихревых течений газа. Согласно этому принципу, если модуль скорости не постоянен всюду, то его максимум достигается на границе и только на границе рассматриваемой области течения. Это свойство используется при разработке формы летательных аппаратов с максимальным критическим значением числа Маха: считается, что если в набегающем потоке и на поверхности обтекаемого тела местное число Маха меньше единицы, то в течении нет звуковых точек. Известное доказательство дозвукового принципа максимума существенным образом опирается на предположение о том, что во всей рассматриваемой области течения давление является функцией плотности. Для идеального (роль диффузии молекул пренебрежимо мала) совершенного (закон Менделеева - Клапейрона) газа давление является функцией плотности, если во всей рассматриваемой области течения энтропийная функция постоянна. Приведен пример дозвукового стационарного безвихревого течения газа, в котором энтропийная функция имеет различные значения на разных линиях тока, а давление не является функцией плотности. Применение дозвукового принципа максимума к такому течению было бы необоснованно. Приведенный пример показывает содержательность вопроса о месте точек максимума модуля скорости дозвуковых стационарных безвихревых неизоэнтропийных течений газа. Для выяснения закономерностей расположения таких точек был проведен анализ полных (без каких-либо упрощающих допущений) уравнений Эйлера в общем пространственном случае. Предложено новое доказательство дозвукового принципа максимума. Это доказательство не опирается на предположение об изоэнтропийности. Тем самым показано, что требование изоэнтропийности можно исключить из условий дозвукового принципа максимума. Дозвуковой принцип максимума оказывается верным и для неизоэнтропийных дозвуковых стационарных безвихревых течений идеального совершенного газа.

Ключевые слова: уравнения Эйлера, дозвуковой принцип максимума, безвихревые течения, изоэнтро-пийность.

ВВЕДЕНИЕ

В движущейся частице идеального (роль диффузии молекул пренебрежимо мала) совершенного (закон Менделеева - Клапейрона) газа давление p связано с плотностью р соотношением p = sp, где энтропийная функция s > 0 не меняется во время движения (здесь k > 1 -показатель адиабаты). Поэтому величина s одинакова во всех точках любой линии тока стационарного течения. В некоторых течениях энтропийная функция может быть различной на разных линиях тока. Если же величина s принимает одно и то же значение во всех точках течения, то такое течение называется изоэнтропийным. Здесь и далее предполагается, что во всей рассматриваемой области движется один и тот же газ, и показатель адиабаты k > 1 имеет одинаковое значение во всех точках. Это относится как к изоэнтропийным, так и неизоэнтропийным течениям. То есть в данной статье причиной неизоэнтропийности считается только различие s на разных линиях тока.

В статье [1] для стационарных дозвуковых безвихревых течений газа, в которых давление является функцией плотности, получен так называемый дозвуковой принцип максимума. Доказательство можно найти в [2, гл. II]. Согласно этому принципу, максимум модуля скорости, если он не постоянен всюду, достигается на границе и только на границе течения. В этой теореме в качестве границ течения могут выступать как поверхности обтекаемых тел, так и воображаемые поверхности, ограничивающие рассматриваемую область течения. Доказательство этого принципа существенным образом опирается на предположение о том, что во всей рассматриваемой области течения давление является функцией плотности. Для изоэнтропийных течений s = const во всех точках, и, следовательно, давление идеального совершенного газа

Vol. 20, No. 02, 2017

Ovil Aviation High Technologies

р — 8рк является функцией плотности. Поэтому в формулировке дозвукового принципа максимума для идеального совершенного газа требование о том, чтобы давление являлось функцией плотности, можно заменить требованием изоэнтропийности.

Применение оператора ротора к (векторному) уравнению Эйлера [3] для безвихревого течения приводит к равенству нулю векторного произведения двух градиентов: [Ур, Ур] —0. Это равенство в случае Ур ^ 0 в (открытой) области означает существование функции р — р(р), а в случае Ур = 0 - постоянство плотности. В первом случае справедлив дозвуковой принцип максимума. Второй случай соответствует безвихревому течению несжимаемой жидкости, и принцип максимума модуля скорости вытекает из теории гармонических функций. В обоих случаях дозвуковой принцип максимума оказывается верным без дополнительного предположения об изоэнтропийности. Отсюда напрашивается вывод о том, что требование изо-энтропийности в формулировке дозвукового принципа максимума является необязательным. Но такой вывод был бы неверным. Это связано с тем, что возможны неизоэнтропийные безвихревые течения, в которых условие Ур ^ 0 нарушено не в (открытых) областях, а в некоторых внутренних точках, линиях или поверхностях. В таких течениях множество точек, в которых Ур — 0, не является открытым, и поэтому невозможно применение теории гармонических функций. Применение же дозвукового принципа максимума также может оказаться невозможным, поскольку давление не будет являться функцией плотности. Поясним это на примере течения газа между двумя соосными цилиндрами. Первый цилиндр имеет радиус г0, второй - радиус 3г0. Расположим декартову прямоугольную систему координат Оххх2хъ так, чтобы ось Ох3 совпадала с осью цилиндров. Обозначим У,У,У - компоненты скорости газа. Непосредственной проверкой можно убедиться, что уравнениям Эйлера [3], уравнению неразрывности

[3] и условию сохранения энтропийной функции р / рк вдоль линий тока удовлетворяет набор функций У1,У2,У3, р и р , вычисленных по формулам

У1 — У0г0Х2 / Г2 > У2 — ~У0Г0Х1 / Г2 > У3 — У0 ,

р — р0 (4г02 + (г - 2г0 )2 ) / г02, р — 5р0 V2 + р0Г02 [п(г / г ) + (4г0/ г) - (4г02 / г2 )], ¡2 2

где г — Лх + х , а положительные константы г0, р0, У0 имеют размерность длины, плотности и скорости соответственно. Предполагается, что координаты Х'х2'хз имеют размерность длины. Эти формулы описывают течение, в котором частицы жидкости движутся по цилиндрическим поверхностям с общей осью симметрии Ох3. Скорость с компонентами имеет нулевую

завихренность. Первое слагаемое в выражении для давления (5р0У0 ) подобрано таким образом,

что в зазоре между цилиндрами (г0 < г < 3г0) выполнено неравенство -^У2 + У2 + У3 кр/р . То есть течение является дозвуковым.

Несложно проверить, что во всем зазоре между цилиндрами (г0 < г < 3г0) давление строго монотонно возрастает с ростом радиуса г . Из выражения для плотности видно, что она строго монотонно убывает с ростом радиуса г на участке г0 < г < 2г0. И наоборот, плотность строго

монотонно возрастает с ростом радиуса г на участке 2г0 < г < 3г0 . Поэтому по отдельности на участках г0 < г < 2г0 и 2г0 < г < 3г0 давление является функцией плотности. Но во всем зазоре между цилиндрами (г0 < г < 3г0) давление не является функцией плотности. Следовательно, к

Сivil Aviation High Technologies

Vol. 20, No. 02, 2017

этому стационарному дозвуковому безвихревому течению нельзя применить дозвуковой принцип максимума.

Приведенный пример показывает существенность требования изоэнтропийности для дозвукового принципа максимума. Однако это связано лишь с особенностями «старого» доказательства [1, 2]. В данной статье показано, что дозвуковой принцип максимума остается верным и для неизоэнтропийных течений, то есть для течений, в которых энтропийная функция может принимать разные значения на разных линиях тока. Тем самым показано, что требование изо-энтропийности является несущественным и может быть исключено из условий дозвукового принципа максимума.

В данной статье рассматриваются области течения, в которых нет разрывов аэродинамических параметров. Достаточная для дальнейшего исследования гладкость всех аэродинамических параметров в таких областях предполагается естественным свойством течений газа, что является одной из общепризнанных гипотез теории непрерывной сплошной среды.

СЛУЧАЙ ОТСУТСТВИЯ ТОЧЕК ТОРМОЖЕНИЯ

Рассмотрим стационарное безвихревое ( rot V = 0, где V - скорость) течение идеального совершенного газа в общем пространственном случае. Такое течение может быть неизоэнтро-пийным. В силу возможной неизоэнтропийности нельзя считать, что полная энтальпия

\ = spk-1k/(k-1) + V2/2 > 0

одинакова на разных линиях тока. Но, независимо от этого, энтропийная функция s не меняется вдоль линий тока. Поэтому из уравнений Эйлера [3] следует, что полная энтальпия h0 сохраняется вдоль каждой линии тока.

В данном разделе будем считать, что в рассматриваемой области течения нет точек торможения, то есть что модуль скорости V = V и, следовательно, параметр Чаплыгина

т = V / sj2h всюду отличны от нуля. Поэтому в каждой точке существует единичный вектор, касательный к линии тока е = V / V. Тогда уравнение неразрывности div(pV) = div(^j2h()рте) = 0, с учетом (е ■ V/z()) = 0, можно представить в виде

div е + (е ■ V) 1п(рт)= 0.

Подставим в это равенство плотность, полученную из выражения для полной энтальпии:

_L J_

р = (1 -т2)k-1 (h (k -1)/ sk) k-1.

Учитывая, что (e ■ V(/z0 (к — 1) / sk)) = 0, получим

dive -2т2(k-\y\l-T2y\e -Ч)1ш + (e-V)lm = 0.

Далее, поскольку (e ■ V) In т = (e ■ V) In (r^2h0 ) =(e ■ V)ln V, имеем

dive = (az2-\)(\-т2)~\ё-УЯ), a = (k + l)/(k-l)>l, A = lnV. (1)

Vol. 20, No. 02, 2017

Civil Aviation High Technologies

Из условия отсутствия завихренности rot(Ve) = 0 находим

rote=[e,VÁ], (2)

Докажем, что величина л — lnV не может достигать максимального значения во внутренней точке течения, если только не является константой. Вычисляя ротор от обеих частей (2) и используя формулу для ротора векторного произведения, получим

rot rot е = (VA- V)e-(e-V) УЛ+ е АЛ- (div e)VÁ ,

где А - оператор Лапласа. Поскольку производная единичного вектора е по любому направлению перпендикулярна е (и, следовательно, (ё-(УЛ-У)ё)-0), последнее уравнение после скалярного умножения на вектор е дает

(е • rot rot е) = -(е- ((е ■ V) УЛ))+ АЛ-(е- У Л) div е .

Преобразуем левую часть этого равенства, используя известное векторное тождество У div е = rot rote + Ае . После перестановки слагаемых получим:

АЛ - (е- У div е)-(е- ((е ■ У) У Л)) - (е ■ У Л) div ё = -(е-Ае).

Второе слагаемое в левой части выразим с помощью (1), учитывая при этом, что в силу (ё-УЬъ) = 0 верно равенство (е ■ У Л) = (е-У1пт). Имеем

Al- (ост2 -1)(1 - г2 )-1 (е ■ У(е ■ У Л)) - 2т2 (а -1)(1 - г2 )"2 (е ■ У Л)2 -

- (ё ■ ((ё ■ V) va» - (ё ■ у л) сИуё =-{ё-Аё). (з)

Рассмотрим множитель (е - У(е - УЛ)) во втором слагаемом левой части этого уравнения. Используя известную формулу для градиента скалярного произведения, имеем

У(ё ■ У Л) = (ё ■ У)УЛ + (У Л- У)ё + [ё, rotyX\ + [У Л, rote] =

= (ё ■ у)ул+(у л- у) ё+[у л, rot ё].

Как уже было замечено выше, производная единичного вектора е по любому направлению перпендикулярна е , и поэтому (е -{УЛ- V)e ) = 0. Далее, используя (2), получаем, что

[у л, rot e¡ =[ул, [ё, у л]\ = (у л ■ул)ё-(ё ■ ул) у л,

и, следовательно,

(ё ■ У(ё ■ у л)) = (ё ■ ((ё ■ V) vi))+(У л ■ у л) - (ё ■ у л)2.

Поэтому уравнение (3) можно записать в виде

АЛ- (1 + (ост2 -1)(1 - т2 )-1 )(е ■ У(е ■ У Л)) + (Ь-УЛ)=-(е- Ае) , (4)

Ovil Aviation High Technologies

Vol. 20, No. 02, 2017

где

b = -{ост2 -1)(1 -T2ylV A +{((.ат2 -1)(1 -T2yl —It2 (a -1 )(\ - т2 )~2)(e ■VA) - div e } e . (5)

Учитывая первое из равенств (1), упростим выражение (5):

Ъ =-(ат2 -l)(l-T2yWA —2т2(а — 1)(1 — т2у2(ё-УЛ)ё .

Выберем прямоугольную декартову систему координат Ох1х2х3. Пусть е1, е2, е3 и ЪХ,Ъ2,Ъ3 - координаты единичного вектора е и вектора b в этой системе. Для любых i = 1, 2, 3; к = 1, 2, 3 верно дифференциальное равенство

2 2 д х 2

/

= 2

de.

2

v dx. ! V i y

+ 2e

d2e

к

k 2~ д x.

Сложим девять таких равенств, соответствующих различным комбинациям i = 1, 2, 3; к = 1, 2, 3:

A (e2 + + e2) = 2(Ve )2 + 2( Ve2 )2 + 2(Ve3 )2 + 2e a e + 2e2 a e2 + 2e3 A e3

или

A(ef + e2+e3) = 2(Vex)2 + 2(Ve2)2 + 2(Ve3 f + 2(e ■ Ae).

Учитывая, что ex, e2, e3 - координаты единичного вектора, имеем A (e2 + e^ + e^) = A (1) = 0. Следовательно:

0 = 2(Welf + 2(Ve2)2 + 2(Ve3f + 2(e ■ Ae). Поэтому уравнение (4) может быть записано в виде

АЯ- (1 + (ост2 -1)(1- г2 )-1 )(е ■ V(e • VA))+ (Я • VA) = (V^ )2 + (Ve2 )2 + (Ve3 )2 .

Обозначим р = 1 + (от2 -1)(1 -т2)-1 и запишем левую часть через частные производные по координатам xt,x2,x3:

2, d2X

d2 2

(1 - ßel ) — + (1 -ße22) d-i + (1 - ße32 ) df - 2ßexe2

2, d2X

dx2

dx2

dx-2

dx1dx2

- 2 ße e

d 2 dx2dx3

- 2ße2e2

d2Ä dx dx

+

7 dX dX dX 2 2 2

+ bi — + b2 — + b3 — = (Vei)2 + (Ve2)2 + (Ve3)2.

dx2

dx^,

dx.

Vol. 20, No. 02, 2017

Civil Aviation High Technologies

Для исследования свойств решений уравнения (6) воспользуемся одним из следствий теоремы Хопфа [4, 5], приведенным в [6]. Сформулируем это следствие.

Пусть во всех точках ограниченной области G коэффициенты аи, а12, а22, а13, а23, а33 уравнения

2 2 2 2 2 2 д u д u д u д u д u д u ап--h а22--h а33--h 2а12--h 2а13--h 2а23--h

0x2 dx 2 dx 2 0x10x2 0x10x3 dx2 dx3

du du du + \-+ Ъ2-+ Ъ3-+ cu = f (7)

0xx дх2 0x3

являются коэффициентами положительно определенной квадратичной формы A. Пусть, далее, для любой точки M(x^, x2, x3) е G существуют а = а(xx, x2, x3 ) > 0 и

S = S(xiyx2,x3) > 0 такие, что замкнутый шар U(M,а) с центром в точке M и радиусом

а целиком лежит в области G и в нем все коэффициенты уравнения (7) ограничены и выпол-

2

няется неравенство detA > S. Тогда в случае c = 0, f > 0, если решение u е C (G) уравнения (7) непрерывно в замкнутой области G, то u < maxu во всей области G. При этом, если

0G

u Ф const в G, то равенство u = maxu возможно только на 0G.

0G

Если уравнение (6) используется для поиска решения уравнений Эйлера, следствием которых оно является, нельзя считать, что коэффициенты при производных являются известными функциями координат. Однако если уравнение (6) используется для выяснения свойств его решения в предположении, что такое решение существует, то эти коэффициенты можно считать некоторыми известными функциями координат. Такой подход (считать коэффициенты известными функциями) неоднократно использовался в доказательствах, например, в [2, 7-10]. Рассмотрим уравнение (6) как уравнение вида (7). Т. е. будем считать коэффициенты при производных Л в левой части (6) заданными функциями координат x>x2,x3.

Поскольку местная скорость газа равна скорости звука при т = 1/ л/а , то во всей области дозвукового течения без точек торможения выполняется неравенство 0 < т < 1/ л/а . Поэтому 0 < Р < 1, и использование критерия Сильвестра показывает, что коэффициенты при вторых производных Л в уравнении (6) являются коэффициентами положительно определенной квадратичной формы. В силу непрерывности параметров течения, величина Р отграничена от единицы. Поэтому существует такое число S > 0, что во всех точках течения выполнено условие

detA = 1 — Р > S. Кроме того, величина 1 — т отграничена от нуля, и все коэффициенты как при вторых, так и при первых производных Л в уравнении (6) ограничены. Правая часть (6) неотрицательна, что в приведенной выше формулировке следствия теоремы Хопфа соответствует требованию f > 0. Следовательно, уравнение (6) и его решение удовлетворяют приведенному выше следствию теоремы Хопфа. Поэтому если величина Л не постоянна, то ее максимум (и, следовательно, максимум модуля скорости) достигается на границе и только на границе области течения. Таким образом, дозвуковой принцип максимума доказан для неизоэнтропийного течения без точек торможения.

Научный Вестник МГТУ ГА_Том 20, № 02, 2017

Civil Aviation High Technologies Vol. 20, No. 02, 2017

СЛУЧАЙ СУЩЕСТВОВАНИЯ ТОЧЕК ТОРМОЖЕНИЯ

Теперь рассмотрим дозвуковое стационарное безвихревое течение идеального совершенного газа, в котором могут быть точки торможения. (Течение может быть неизоэнтропий-ным.) Дозвуковой принцип максимума будет доказан, если из предположения о достижении максимума V во внутренней точке ограниченной области течения последует, что величина V одинакова во всех точках этой области.

Допустим, что в некоторой внутренней точке P ограниченной области течения G модуль скорости достигает своего максимального значения Vmax . Достаточно рассмотреть

случай Vmax > 0 .

Точка P по предположению является внутренней, то есть существует такая окрестность U точки P, в которой V > Vmax / 2. Эта окрестность является областью (открытым связным множеством), в которой нет точек торможения. Поэтому непустым будет множество G' œ G, определяемое как объединение всех областей, в которых V > Vmax / 2 и содержащих точку P.

Это объединение является областью, и в этой области, согласно доказанному выше варианту дозвукового принципа максимума (при отсутствии точек торможения), скорость всюду

равна Vmax .

Допустим, что область G" не совпадает с областью G . Тогда существует точка P еdG являющаяся внутренней точкой G. В силу непрерывности скорости V(P ) = Vmax . Следовательно, существует окрестность U' œ G точки P", в которой V > Vmax / 2 . Но тогда объединение G" = G' u U' является областью, содержащей точку P, и в которой V > Vmax / 2. Область G"

содержит в себе область G и не совпадает с ней. Это противоречит определению области G . Полученное противоречие показывает, что G" = G. Итак, доказано следующее утверждение, в котором отсутствует требование изоэнтропийности.

Дозвуковой принцип максимума. В дозвуковом стационарном безвихревом течении идеального совершенного газа модуль скорости, если не является константой во всей ограниченной области течения G, достигает своего максимального значения на границе и только на границе области G.

ЗАКЛЮЧЕНИЕ

Известное ранее доказательство дозвукового принципа максимума [1, 2] не позволяло исключить изоэнтропийность из условий этого принципа. Приведено новое доказательство. В этом доказательстве свойство изоэнтропийности не используется. То есть требование изоэнт-ропийности является несущественным для справедливости дозвукового принципа максимума. Таким образом, если модуль скорости (дозвукового стационарного безвихревого течения идеального совершенного газа) не постоянен всюду, то его максимум достигается на границе и только на границе рассматриваемой области течения. Полученный результат расширяет область применения дозвукового принципа максимума для качественного анализа течений.

СПИСОК ЛИТЕРАТУРЫ

1. Shiffman M. On the Existence of Subsonic Flows of a Compressible Fluid. J. Ration. And Analysis. 1952, vol. 1, pp. 605-652.

2. Берс Л. Математические вопросы дозвуковой и околозвуковой аэродинамики. М.: ИЛ, 1961. 208 с.

3. Ламб Г. Гидродинамика. М.-Л.: ОГИЗ ГИТТЛ, 1947. 256 с.

Vol. 20, No. 02, 2017

Civil Aviation High Technologies

4. Hopf E. Elementare Bemerkungen uber die Losungen partieller Differentialgleichungen zweiter Ordnung vom Elliptischen Typus. Sitzungsberichte der Preussischen Akademie der Wissenschaften. 1927, vol. 19, pp. 147-152.

5. Миранда К. Уравнения с частными производными эллиптического типа. М.: Издательство иностранной литературы, 1957. 256 с.

6. Беспорточный А.И., Бурмистров А.Н., Сизых Г.Б. Вариант теоремы Хопфа // ТРУДЫ МФТИ. 2016. Т. 8, N 1. С. 115-122.

7. Gilbarg D., Shiffman M. On Bodies Achieving Extreme Value of the Critical Mach Number. I. J. Ration. And Analysis. 1954, vol. 3, no. 2, pp. 209-230.

8. Бурмистров А.Н., Ковалёв В.П., Сизых Г.Б. Принцип максимума для решения уравнения эллиптического типа с неограниченными коэффициентами // ТРУДЫ МФТИ. 2014. Т. 6, № 4. С. 97-102.

9. Сизых Г.Б. Признак наличия точки торможения в плоском безвихревом течении идеального газа // ТРУДЫ МФТИ. 2015. Т. 7, № 2 (26). С. 108-112.

10. Голубкин В.Н., Сизых Г.Б. Принцип максимума функции Бернулли // Ученые записки ЦАГИ. 2015. Т. 46, N 5. С. 53-56.

СВЕДЕНИЯ ОБ АВТОРЕ

Сизых Григорий Борисович, кандидат физико-математических наук, доцент МФТИ, [email protected].

MAXIMUM PRINCIPLE FOR SUBSONIC FLOW WITH VARIABLE ENTROPY

Grigory B. Sizykh1

1Moscow Institute of Physics and Technology (State University), Moscow, Russia

ABSTRACT

Maximum principle for subsonic flow is fair for stationary irrotational subsonic gas flows. According to this principle, if the value of the velocity is not constant everywhere, then its maximum is achieved on the boundary and only on the boundary of the considered domain. This property is used when designing form of an aircraft with a maximum critical value of the Mach number: it is believed that if the local Mach number is less than unit in the incoming flow and on the body surface, then the Mach number is less then unit in all points of flow. The known proof of maximum principle for subsonic flow is based on the assumption that in the whole considered area of the flow the pressure is a function of density. For the ideal and perfect gas (the role of diffusion is negligible, and the Mendeleev-Clapeyron law is fulfilled), the pressure is a function of density if entropy is constant in the entire considered area of the flow. Shows an example of a stationary subsonic irrotational flow, in which the entropy has different values on different stream lines, and the pressure is not a function of density. The application of the maximum principle for subsonic flow with respect to such a flow would be unreasonable. This example shows the relevance of the question about the place of the points of maximum value of the velocity, if the entropy is not a constant. To clarify the regularities of the location of these points, was performed the analysis of the complete Euler equations (without any simplifying assumptions) in 3-D case. The new proof of the maximum principle for subsonic flow was proposed. This proof does not rely on the assumption that the pressure is a function of density. Thus, it is shown that the maximum principle for subsonic flow is true for stationary subsonic irrotational flows of ideal perfect gas with variable entropy.

Key words: Euler equations, subsonic maximum principle, irrotational flow, flow with variable entropy.

REFERENCES

1. Shiffman M. On the Existence of Subsonic Flows of a Compressible Fluid. J. Ration. And Analysis. 1952, vol. 1, pp. 605-652.

Civil Aviation High Technologies

Vol. 20, No. 02, 2017

2. Bers L. Mathematical Aspects of Subsonic and Transonic Gas Dynamics. Surveys in Applied Mathematics. Vol. 3. New York, John Wiley & Sons Inc. 1958. 164 p.

3. Lamb H. Hydrodynamics. Cambridge University Press. 1895. 604 p.

4. Hopf E. Elementare Bemerkungen uber die Losungen partieller Differentialgleichungen zweiter Ordnung vom Elliptischen Typus. Sitzungsberichte der Preussischen Akademie der Wissenschaften. 1927, vol. 19, pp. 147-152.

5. Miranda C. Partial Differential Equations of Elliptic Type. Berlin-Heidelberg-New York, Springer-Verlag. 1970. 370 p.

6. Besportochnyy A.I., Burmistrov A.N., Sizykh G.B. Variant teoremy Hopfa [Version of the Hopf Theorem]. Trudy MFTI [Proceedings of MIPT], 2016, vol. 8, no. 1, pp. 115-122. (in Russian)

7. Gilbarg D., Shiffman M. On Bodies Achieving Extreme Value of the Critical Mach Number. I.J. Ration. And Analysis. 1954, vol. 3, no. 2, pp. 209-230.

8. Burmistrov A.N., Kovalev V.P., Sizykh G.B. Printsip maksimuma dlia resheniia uravne-niia ellipticheskogo tipa neogranichennymi koeffitsientami [Maximum Principle for the Solution of an Elliptic Equation With Unbounded Coefficients]. Trudy MFTI [Proceedings of MIPT], 2014, vol. 6, no. 4, pp. 97-102. (in Russian)

9. Sizykh G.B. Priznak nalichiia tochki tormozheniia v ploskom bezvihrevom techenii ideal-nogo gaza [Sign of the Presence of Braking Points in a Planar Irrotational Flow of an Ideal Gas]. Trudy MFTI [Proceedings of MIPT], 2015, vol. 7, no. 2, pp. 108-112. (in Russian)

10. Golubkin V.N., Sizykh G.B. Maximum Principle for Bernoulli Function. TsAGI Science Journal. 2015, vol. 46, issue 5, pp. 485-490.

INFORMATION ABOUT THE AUTHOR

Grigory B. Sizykh, PhD, Associate Professor of Moscow Institute of Physics and Technology, o1o2o3@yandex .ru.