CC BY
34
УДК 532.516
В. П. Ковалёв1, Е. Ю. Просвиряков2, Г. Б. Сизых1
1 Московский физико-технический институт (государственный университет) 2Институт машиноведения УрО РАН
Получение примеров точных решений уравнений Навье—Стокса для винтовых течений методом суммирования скоростей
В рамках уравнений Навье-Стокса рассмотрены винтовые течения вязкой несжимаемой жидкости в потенциальном поле внешних сил. Винтовыми в статье называются течения, в которых завихренность параллельна скорости, а отношение величины завихренности к величине скорости постоянно по пространству и по времени. Показано, что можно суммировать векторные поля скоростей двух винтовых решений, если у этих решений отношение величины завихренности к величине скорости одинаково. В результате получается поле скорости некоторого «нового» винтового решения уравнений Навье-Стокса. Такой метод получения новых решений назван «методом векторного суммирования». В качестве «исходных» точных решений уравнений Навье-Стокса рассмотрены, в частности, решения Громеки-Бельтрами-Тркала (то есть решения, которые получены методом Тркала из решений Громеки-Бельтрами). Метод векторного суммирования позволяет складывать скорости таких решений Громеки-Бельтрами-Тркала, у которых оси симметрии не совпадают. В результате получается неосесимметричное точное решение. Приведен пример нового точного решения, полученного методом векторного суммирования трех известных точных решений. Два из этих трех известных решений - это два разных решения Громеки-Бельтрами-Тркала, у которых оси симметрии не совпадают. Третье решение - это решение уравнений Навье-Стокса, полученное методом Тркала из ABC-решения (Arnold-Beltrami-Childress). Полученное новое (суммарное) решение имеет более сложную структуру, чем три «исходных» решения. Это новое решение не периодично. У него нет оси симметрии, нет плоскости симметрии и нет центра симметрии. Точные решения, полученные методом векторного суммирования, могут использоваться для тестирования численных алгоритмов и компьютерных программ.
Ключевые слова: точные решения уравнений Навье-Стокса, винтовые решения уравнений Навье-Стокса, точные решения уравнений Эйлера, вязкая несжимаемая жидкость, осесимметричные течения.
V. P. Kovalev1, E. Yu. Prosviryakov2, G. B. Sizykh1
1 Moscow Institute of Physics and Technology (State University) 2Institute of Engineering Science, UB RAS, Ekaterinburg, Russia
Obtaining examples of exact solutions of the Navier-Stokes equations for helical flows by the method of summation of velocities
Corkscrewlike flows of a viscous incompressible fluid in the potential field of external forces are considered in terms of the Navier-Stokes equations. (In this paper, flows are considered to be corkscrewlike when their vorticity is parallel to the velocity and the vorticity-to-velocity ratio is constant in space and time.) It is shown that if this ratio is equal for two corkscrew-type solutions, the vector sum of velocities for these solutions is the velocity for some corkscrew-type solution of the Navier-Stokes equations. This property is proposed to be used for obtaining new exact solutions. It is demonstrated that the vector summation yields solutions with new properties, which are not inherent to original solutions). In particular, the nonaxisymmetric exact solution can be obtained by the summation of axisymmetric Gromeka-Beltrami-Trkal solutions, i. e. solutions obtained by Trkal's method
from Gromeka-Beltrami solutions). Exact solutions obtained by vector summation can be used to test numerical algorithms and computer programs.
Key words: exact solutions of the Navier-Stokes equations, helical solutions of the Navier-Stokes equations, exact solutions of the Euler equations, incompressible viscous fluid, axially symmetric flow.
1. Введение
Недавно была обобщена теорема Томсона [1] о сохранении циркуляции по контуру, состоящему из одних и тех же частиц идеальной баротропной жидкости. Для небаротропной и вязкой жидкости найдена скорость движения воображаемых частиц. При этом движении сохраняется циркуляция скорости жидкости по контуру, состоящему из воображаемых движущихся частиц [2].
Этот новый закон сохранения еще не заложен в основу численных методов. Поэтому в работе [2] предложено использовать его для проверки численных расчетов. Численные решения, как правило, получаются с некоторой ошибкой - с некоторым отклонением от точного решения. Точное выполнение закона сохранения на численном решении невозможно. Для выбраковки неприемлемых решений необходимо решить вопрос о величине ошибки, превышение которой свидетельствует о нарушении закона сохранения.
Вопросы подобного рода еще не изучались. В качестве первого шага предполагается провести ряд численных экспериментов для выяснения точности, с которой выполняется закон сохранения [2] в случае, когда численное решение имеет некоторую «приемлемую» точность. В рамках такого эксперимента требуется вычислять как ошибку выполнения закона сохранения, так и отклонение численного решения от точного решения. Для этого требуется запас различных точных решений. Следует заметить, что некоторые точные решения трудно реализовать в эксперименте или наблюдать в реальной обстановке. Это связано со сложностью создания граничных условий (которые должны быть взяты из самих решений). Но если речь идет о верификации, то есть о проверке точности численного алгоритма, не обязательно иметь техническую возможность реализовать краевые условия в эксперименте. Большинство численных алгоритмов работают с любыми граничными условиями. Поэтому для упомянутых выше численных экспериментов можно использовать и такие точные решения, для которых трудно создать реальные граничные условия, и которые до сих пор считались практически непригодными. К ним относятся точные нестационарные решения уравнений Навье-Стокса, полученные методом Тркала из стационарных винтовых решений уравнений Эйлера. Винтовыми в настоящей статье называются решения, у которых векторы скорости и завихренности связаны постоянным по пространству и времени скалярным коэффициентом. Тркал в статье [3] в общем пространственном случае показал, что для несжимаемой жидкости скорость стационарного винтового решения уравнений Эйлера, умноженная на exp(-ai), где t - время, при надлежащем выборе параметра a будет удовлетворять уравнениям Навье-Стокса. Считается [4], что этот результат Тркала является логическим продолжением идей, высказанных Тейлором и Кальдонаццо. В качестве стационарных винтовых решений уравнений Эйлера можно использовать, например, ABC-решение или осесимметричные решения Громеки-Бельтрами [5, 6], которые были получены независимо Громекой и Бельтрами в 1881 и в 1889 годах соответственно. (Текст [5] можно найти в [7].) В этих решениях проявляется сложное взаимодействие движений жидкости в разных направлениях. Однако решения Громеки-Бельтрами имеют ось симметрии, а ABC-решение обладает свойством периодичности: оно состоит из одинаковых кубических ячеек.
Отметим, что известны и другие винтовые решения - точные решения Богоявленского, обладающие кристаллографической симметрией [8-10]. Этот класс течений описывает со-литонные течения на сфере, которые в общем случае являются квазипериодическими по
координатам [8-10]. Кроме того, в недавней работе [11] приведены новые геометрические алгоритмы построения винтовых векторных полей. Это позволило построить новые классы точных винтовых решений уравнений Эйлера [12].
Представляется, что результаты численных экспериментов по выяснению «приемлемой» точности будут более информативными, если в запас точных решений включить точные решения с более сложной структурой, в которой нет симметрии и периодичности. Для этого в данной работе предложено использовать метод получения таких сложных решений, названный методом векторного суммирования. Этот метод был известен ранее для винтовых решений уравнений Эйлера. Основная научная новизна данной работы состоит в обосновании возможности использовать этот метод для винтовых решений уравнений Навье-Стокса.
С теоретической точки зрения метод векторного суммирования очень прост. Но его практическая реализация связана с громоздкими формулами, которые невозможно найти в одной книге. Учитывая это, в четвертом разделе и в Приложениях приведены справочные сведения, достаточные для компьютерного представления новых точных решений, и приведен пример представления одного из таких решений. Это сделано для того, чтобы данную статью можно было использовать как пособие по построению новых точных решений с целью верификации компьютерных программ.
Для знакомства с теоретической частью данной статьи можно опустить чтение четвертого раздела и Приложений.
2. Уравнения движения
Рассмотрим течение ньютоновской вязкой несжимаемой жидкости в потенциальном поле массовых сил. Будем использовать следующие безразмерные переменные: V - скорость; О = rot V - завихренность; р - давление, отнесенное к плотности; П - потенциал объемных сил; Re - число Рейнольдса, t - время. С учетом равенства AV = — rot О, верного для несжимаемой жидкости, движение описывается уравнениями Навье-Стокса в форме [1]:
V + ^ rot О + О х V = —V
ot Re
2
V
Р +
div V = 0. (1)
В случае ^ V = 0 при нулевой вязкости эти уравнения упрощаются и принимают вид
уравнений Эйлера для несжимаемой жидкости [13]:
•2
О х V = —V
V
р + —+ п
, div V = 0. (2)
Здесь и далее знак волны над символом означает, что речь идет о решении уравнений Эйлера.
В данной работе уравнения (1) и (2) будут рассматриваться в ограниченной замкнутой области. Все выкладки и рассуждения будут проводиться в предположении, что компоненты скорости, давление и потенциал объемных сил имеют непрерывные вторые производные по пространственным координатам и непрерывные первые производные по времени. Чтобы далее не обсуждать вопросы гладкости функций, заметим, что это требование для компонент скорости и для давления будет выполнено в полученных ниже точных решениях при условии, если оно будет выполнено для потенциала объемных сил.
3. Получение точных решений уравнений Навье-Стокса методом Тркала
Говорят, что скорость V удовлетворяет уравнениям Навье-Стокса (уравнениям Эйлера), если существует давление р, составляющее вместе с этой скоростью решение уравнений Навье-Стокса (уравнений Эйлера). Соленоидальность (ё1у V = 0) является необходимым условием того, чтобы скорость удовлетворяла уравнениям Навье-Стокса. Но этого
недостаточно. Если взять произвольную соленоидальную скорость, то может оказаться невозможным «подобрать» соответствующее давление р. Поэтому не всякая соленоидаль-ная скорость удовлетворяет уравнениям Навье-Стокса. Аналогичная ситуация имеет место и для уравнений Эйлера.
В данной работе исследуются винтовые решения. Это решения, в которых
О = rot V = k V, к = const. (3)
(Здесь равенство к = const означает постоянство и в пространстве и во времени.) Скорость, удовлетворяющую условию (3), будем называть винтовой.
Заметим, что винтовые решения являются подклассом течений Бельтрами, в которых Ох V = 0, а величина к, вообще говоря, зависит от пространственных координат и времени.
Для винтовых течений уравнения Навье-Стокса (1) и уравнения Эйлера для стационарной скорости (2) имеют соответственно вид
д 1
—V + — rot О = -V
ot Re
2
V
Р
, div V = 0, (4)
0=
' V2 "
Р + ^ + П
, div V = 0. (5)
Из уравнения (5) следует важный для поиска точных винтовых решений вывод, позволяющий при поиске скорости точного решения исключить из рассмотрения давление и потенциал объемных сил. Для винтовых решений поиск стационарной скорости, удовлетворяющей уравнениям Эйлера для несжимаемой жидкости, сводится к поиску стационарной скорости V, удовлетворяющей уравнению (3) и уравнению неразрывности div V = 0. Если такая скорость найдена, то давление может быть вычислено по формуле
V2
р = Ро - ~Y - п, (6)
где ро - произвольная константа, обеспечивающая неотрицательность давления р в рассматриваемой области течения. Выбор такой константы ро всегда возможен в ограниченной замкнутой области в силу требований к гладкости, сформулированных в конце второго раздела. Существование давления (6), удовлетворяющего уравнениям (5), означает, что скорость V удовлетворяет уравнениям Эйлера. Такие решения ("V,р) будем называть стационарными винтовыми решениями уравнений Эйлера.
Метод Тркала. Пусть р, V и О = rotV = k V - давление, скорость и завихренность винтового решения уравнений Эйлера со стационарной скоростью. Рассмотрим зависящую от времени скорость V = Vexp(-tk2/Re). Эта скорость соленоидальна, поскольку div V = exp(-tk2/Re) div V = 0. Завихренность этой скорости
О = rotV = rot V exp(-ifc2/Re) = к V exp(-ifc2/Re) = к V.
Следовательно, для скорости V выполнено условие (3) и rot О = к2 V. Поэтому левая часть (4) равна
д 1 ^2 Ъ2
—V +— rot О = - — Vexp(-ifc2/Re) + — V = 0. Ot Re Re ' Re
Правая часть (4) будет равна нулю, если давление вычислить по формуле
V2 V2
p = Ро - — - П = ро - — exp(-2i к2/Re) - П. (7)
Таким образом, (V,p) будет винтовым решением уравнений Навье-Стокса. Такой способ получения решений уравнений Навье-Стокса предложен Тркалом в [3].
и
Обратим внимание на отличие давлений р и р. При постоянном во времени поле П давление р непостоянно во времени и вычисляется по формуле (7), а давление р не зависит от времени и вычисляется по формуле (6).
Итак, каждому стационарному винтовому решению уравнений Эйлера соответствует нестационарное винтовое решение уравнений Навье-Стокса. Учитывая это, будем искать скорость V стационарных винтовых решений уравнений Эйлера. И не будем каждый раз напоминать, что скорость V = V exp(-tk2/Re) и давление, вычисленное по формуле (7), будут решениями уравнений Навье-Стокса.
4. ABC-решение и решения Громеки—Бельтрами
Напомним известные способы суммирования стационарных винтовых решений на примере ABC-решений и решений Громеки-Бельтрами.
ABC-решение. Пусть Oxyz - прямоугольная декартова система координат; i, j, k - базисные векторы этой системы. И пусть А, В, С, к - произвольные константы. Непосредственной проверкой можно убедиться, что если
то скорость V = u i + v j + w k будет удовлетворять уравнениям div V = 0 и rot V = к V. То есть V будет скоростью винтового решения уравнений Эйлера. В виде (8) это решение предложено в [14]. Однако частные случаи ABC-решения были известны значительно раньше. Например, в 1919 году (для случая А = 1, В = С = 0) оно было приведено в статье [3] в качестве примера винтового течения.
Суммирование ABC-решений. Различным наборам чисел А, В, С, к соответствуют различные винтовые решения. Сумма решений, соответствующих наборам А\, В\, С\,к\ и А2, В2, С2,к2, в случае к\ = к2 дает новое винтовое решение, соответствующее набору А\ + ^2, В\ + В2, С\ + С2, к\. Если же к\ = к2, то суммарная скорость не будет винтовой скоростью и в общем случае не будет удовлетворять уравнениям Эйлера.
Решения Громеки—Бельтрами. Введем цилиндрическую систему координат г, (р, х с началом в точке В и будем рассматривать симметричные относительно оси В% течения. «Необычное» использование греческих букв В и % объясняется тем, что в следующем разделе эта система координат будет рассматриваться одновременно с прямоугольной декартовой системой координат Oxyz, оси которой могут не совпадать с осью В%.
Параметры осесимметричного течения не зависят от окружной координаты (р. Пусть er, e<^, ex - правая тройка единичных векторов в радиальном, окружном и осевом направлениях соответственно. Вектор скорости закрученного (Vv = 0) течения имеет вид: V = Vrer + V^e^ + Vxex, где функции Vr, Vv, Vx зависят от г, х и не зависят от ip. Пусть к - произвольная константа, а функция ф = ф(г, х) удовлетворяет уравнению
В работах [5, 6] показано, что для произвольного потенциала П = П(г, %) набор функций
u = A sin(fcz) + С cos(ky), v = В sin(fc^) + A cos(fcz), w = С sin(fcy) + В cos(kx),
(8)
(9)
(10)
(11)
р = ро - 0.5 (Vr2 + V^2 + Vx2) - П
является винтовым решением уравнений Эйлера, а функция ф = ф(г, х) постоянна вдоль линий тока и поэтому называется функцией тока.
Таким образом, любое решение уравнения типа (9) «порождает» точное винтовое решение уравнений Эйлера для несжимаемой жидкости.
Все решения уравнения (9) до сих пор не найдены. Громека [5] и Бельтрами [6] искали решения в виде произведения ф = ф(г,х) = R(r)X(x) (метод разделения переменных). Но в их статьях представлены не все решения, которые можно получить таким методом. Во-первых, они не рассматривали решения, которые имеют особенность на оси симметрии. Во-вторых, статьи Громеки и Бельтрами опубликованы в то время, когда еще не были исследованы цилиндрические функции мнимого аргумента (Macdonald functions). Поэтому в статьях Громеки и Бельтрами отсутствуют решения, которые выражаются с использованием таких функций.
Приведем все известные в настоящее время решения (9), которые можно получить методом разделения переменных. Все эти решения следует называть решениями Громеки-Бельтрами.
Для каждого фиксированного числа к = 0 решением ф = ф(г, х) уравнения (9) является произвольная линейная комбинация функций вида
rl\(pir) sin(miх), rl\(pir) cos(miх), ш2 = к2 + р\;
rK1(p2r) sin(m2х), rK1(p2r) cos(m2х), т2 = к2 + р2,; rJi(pir) sin(m3х), rJi(p3r) cos(m3х), тЗ, = к2 - рЗ; rY1(p4г) sin(m4х), rY1(p4r) cos(m4х), = к2 - pi;
rJ1(p5r) exp(m5х), ^Y1(p5r) exp(m5х), m2i = р5 - к2; sin(fc%); cos(fc%); г2 sin(fc%); г2 cos(fc%); r%J1(|fc| r); r %Y1(|fc| г), (13)
где Р1 р2 - произвольные положительные константы, а константы рз ,р4, р§ положительны и выбраны так, что р2 < к2, р\ < к2, р2 ^ к2. J1, Y1 - функции Бесселя и Неймана первого порядка; I1 и К1 - цилиндрические функции мнимого аргумента первого порядка.
Некоторые из приведенных выше решений Громеки-Бельтрами имеют особенность на оси симметрии. Они могут использоваться для тестирования программ расчета течений в зазоре между двумя соосными поверхностями вращения. Такие течения не содержат ось симметрии.
Суммирование решений Громеки—Бельтрами с общей осью симметрии.
Приведем пример применения известного метода суммирования винтовых решений с общей осью симметрии. Этот метод указан в работах [5, 6] и основан на возможности сложения различных решений уравнения (9), полученных для одинаковых значений коэффициента к. Рассмотрим линейную комбинацию двух функций тока для случая к = 2: ф' = гЛ(^3г) sin(x) и ф" = г J1(\/5r) exp(-%). Соответствующие скорости V = v;er + V^ev + V^ex и V = VГ'er + ev + V"ex вычисляются по формулам (10) -(12):
1 д ib' vl(r,x) = 1 о^Ф' = M^3r) cos (x),v£(r,x) = -к f = -2^(^3г) sin(X),
1 д
v£(r,x) = --Q-Ti>' = -\/3Зо(^3т) sin(x); (14)
V?(r,X) = -Л(Лг) exp(-x),V£(r,x) = -2^(Лг) exp(-x),
V"(r,x) = -V5MV5r) exp(-x), (15)
где ,1о - функция Бесселя нулевого порядка. На рис. 1 показаны в осевой плоскости линии тока радиально-осевых скоростей Vrx = Vler + V^ex и Vr% = V"er + ex каждого из этих течений.
' 3 X ' ' X
Рис. 1. Линии тока радиально-осевых скоростей в осевой полуплоскости
Оба поля скорости V и V имеют общую ось симметрии и одинаковый коэффициент к = 2. В силу линейности уравнения (9) сумма скоростей V = V + V соответствует функции ф = ф1 + ф" = rJl(V3r) 8ш(%) + г Jl(V5r) ехр(—\) и поэтому является решением уравнений Эйлера. Линии тока суммарного течения изображены в двух разных масштабах на рис. 2. (Этот пример и рис. 1 и 2 заимствованы из статьи [15]).
10 1 13+ г^ со os ю 15 го 15 as зе ^^
Рис. 2. Линии тока суммарного течения в осевой полуплоскости в разных масштаба
Этот пример показывает, что различные линейные комбинации скоростей стационарных винтовых решений уравнений Эйлера с общей осью симметрии и с одинаковым коэффициентом к могут давать различные и интересные сложные течения. Еще раз заметим, что возможность получать новые решения с помощью такого суммирования давно известна [5, 6]. При этом все суммарные решения остаются осесимметричными.
5. Векторное суммирование решений
Рассмотрим две разные прямоугольные декартовы системы координат Oxyz и Оxyz. Пусть i,j,k и i, j, k - базисные векторы этих систем координат. Пусть, далее, А, В, C,k -произвольные скалярные константы. Запишем ABC-решение в системе Oxyz:
u = Asin(kz) + С cos(ky), v = В sin(kx) + A cos(kz), w = С sin(ky) + В cos(kx).
Вектор скорости V = у, i + v j + w k может быть представлен через базисные векторы i,j,k, то есть может быть представлен в виде: V = u i+vj + w k. В Приложении 1 приведены
формулы для расчета компонент u, v, w через координаты х, у, z. В общем случае точки О и О могут не совпадать, а тройки базисных векторов i,j,k и i, j, к могут иметь произвольное взаимное расположение. Поэтому в системе Oxyz скорость V = и i + v j + w k является «сдвинутым и повернутым» ABC-решением. В Приложении 1 показано, что в общем случае взаимного расположения систем координат Oxyz и Оxyz невозможно подобрать константы А, В, С, к так, чтобы компоненты и, v, w оказались представлены в виде (8). То есть в системе координат Oxyz скорость V = и i + vj + w k в общем случае не будет являться ABC-решением (поскольку в общем случае ее невозможно записать в виде (8)). Поэтому формально можно считать, что «сдвинутое и повернутое» ABC-решение является новым точным решением для тестирования численных алгоритмов. Однако по своей физической сущности это решение остается ABC-решением. В частности, оно, как и «классическое» ABC-решение, состоит из периодически расположенных одинаковых ячеек с одинаковыми картинами течения.
Аналогичная ситуация имеет место для решений Громеки-Бельтрами. Для тестирования численного алгоритма можно «замаскировать» осесимметричность решения Громеки-Бельтрами. Достаточно расположить ось симметрии решения Громеки-Бельтрами так, чтобы она не совпадала с осями прямоугольной декартовой системы координат Oxyz. А затем выразить компоненты скорости и давление через х, у, z. Техника такого представления описана в Приложении 1. Компьютерная программа не сможет воспользоваться осесимметричностью решения.
У такого подхода есть один существенный недостаток. Дело в том, что ось симметрии хотя и «замаскирована», но все-таки существует. Поэтому в течении могут отсутствовать процессы, которые имеют место в общем пространственном случае. В этой связи заметим, что симметрия иногда существенно упрощает свойства течения. Так, например, в статье [16] для осесимметричных течений вязкой жидкости удалость получить аналоги теорем Гельмгольца о вихрях, которые в общем пространственном случае верны только для идеальной жидкости.
Предложим способ получения несимметричных и непериодичных (по пространственным координатам) решений. То есть непериодичных течений, у которых нет оси симметрии, нет плоскости симметрии и нет центра симметрии. Как было показано выше, разные ABC-решения с одинаковым коэффициентом к в уравнении (3) можно складывать и получать в результате новое ABC-решение. Это эквивалентно сложению векторов скорости ABC -решений, записанных в одной системе координат Oxyz. Аналогичным свойством обладают решения Громеки-Бельтрами (с одинаковым коэффициентом к). Получение функции тока нового решения Громеки-Бельтрами сложением двух функций тока ф' и ф" эквивалентно сложению векторов скорости решений Громеки-Бельтрами, имеющих общую ось симметрии. Оказывается, что такая возможность сложения векторов является следствием более общего свойства, присущего всем винтовым течениям.
Пусть V i и V2 - две соленоидальные винтовые скорости с одинаковым коэффициентом к:
div V1 = div V2 = 0, rot V1 = кV1, rot V2 = кV2, к = const.
Например, Vi может быть скоростью ABC-решения, а V2 - скоростью одного из решений Громеки-Бельтрами. (У этих решений должны быть равны коэффициенты к.)
Рассмотрим векторную линейную комбинацию V = aVi + bV2, где a и b - произвольные скалярные константы. Она соленоидальна:
div( aV 1 + bV 2) = a div V1 + b div V 2 = 0,
или
divV = 0. (16)
Завихренность этой скорости равна
rot (aV 1 + bV2) = rot (aV 1) + rot (bV2) = a к V1 + bk V2 = к V,
или
rot V = k V.
(17)
Как показано в третьем разделе (абзац после формулы (5)), условий (16) и (17) достаточно, чтобы скорость V являлась скоростью винтового решения уравнений Эйлера. Поэтому любая векторная линейная комбинация скоростей винтовых решений уравнений Эйлера с одинаковым коэффициентом к в уравнении (3) является скоростью винтового решения уравнений Эйлера с тем же коэффициентом к.
В частности, это означает возможность суммирования скоростей решений Громеки-Бельтрами со скоростями ABC-решений, а также возможность суммирования скоростей нескольких решений Громеки-Бельтрами с несовпадающими осями симметрии. Представляется естественным назвать такой метод получения новых винтовых решений методом векторного суммирования. Заметим, что метод векторного суммирования «автоматически» распространяется на решения уравнений Навье-Стокса, полученные методом Тркала.
Несмотря на то, что метод векторного суммирования давно известен для винтовых течений идеальной несжимаемой жидкости [17], возможность его использования для вязких течений не была замечена ранее и представляется впервые в данной работе. Так, например, в статьях [18, 19] обсуждалась возможность суммирования различных решений уравнений Навье-Стокса. Авторы этих статей указывают лишь на возможность умножения на константу. Другими словами, они заметили, что V = aVi будет скоростью винтового решения уравнений Навье-Стокса, если Vi является скоростью винтового решения уравнений Навье-Стокса. Однако «полноценная» линейная комбинация V = aV i + bV2 винтовых решений в работах [18, 19] не исследована.
Для практического использования метода векторного суммирования можно записывать разные винтовые решения уравнений Эйлера в одной прямоугольной декартовой системе координат Oxyz. И затем складывать компоненты скоростей этих решений.
Продемонстрируем возможности метода векторного суммирования на примере. Сложим компоненты скоростей трех решений с одинаковым коэффициентом к = 2 в одной прямоугольной декартовой системе координат Oxyz.
Первое решение. Это решение Громеки-Бельтрами, соответствующее функции ф', рассмотренной выше в четвертом разделе. Ось симметрии В'этого решения совместим с осью Oz так, как это показано на рис. 3. (Единичный вектор ех/, направленный вдоль оси В'%', в базисе i,j,k основной системы координат представляется в виде: ех' =0 ■ i + 0 ■ j + 1 ■ к.) В цилиндрической системе координат В'г'<р'х' функция тока ф' имеет вид ф1 = ф'(г', х') = г' Ji(V3r') sin(x').
3 У
Рис. 3. Расположение осей трех решений в системе координат Охуг
Второе решение. Это решение Громеки-Бельтрами, соответствующее функции ф", рассмотренной в четвертом разделе. Зададим еще одну цилиндрическую систему координат В"г"<р"~х". Точку В" расположим так, чтобы в системе Охуг она имела координаты (2, 0, 0).
Единичный вектор ex», направленный вдоль оси В''х", в базисе i,j,k основной системы координат зададим следующим образом: ex» = 1/л/3 ■ i + 1/л/3 ■ j + 1/л/3 ■ k. Положение оси В''х" показано на рис. 3. В цилиндрической системе координат В''г''(р"х!' функция тока ф" имеет вид: ф" = ф''(г'',х") = г'' Ji(V5r'') exp(-У).
Третье решение. Это «сдвинутое и повернутое» ABC-решение, которое является «классическим» ABC-решением с параметрами А = 1, В = 2, (7 = 3, к = 2 в системе координат Oxyz. Систему координат Oxyz расположим следующим образом. Точку О расположим так, чтобы в системе Oxyz она имела координаты (3, 3, 0). Оси системы Оxyz получим вращением осей системы Oxyz вокруг оси Оу на угол к/4, как показано на рис. 3.
Компоненты скоростей этих трех винтовых решений в системе Oxyz обозначим соответственно u', v', w', и", v'', w'' и и, v, w . Эти три скорости являются соленоидальными винтовыми скоростями с одинаковым коэффициентом к = 2. Поэтому их сумма
Vs = (и' + и" + и) i + (v' + v'' + v) j + (w' + w" + w)k (18)
является скоростью винтового решения уравнений Эйлера с тем же коэффициентом к. В Приложении 2 получены формулы для вычисления и', v', w', и", v", w'' и и, v, w через координаты x,y,z. Представлять отдельно громоздкие формулы для величин (и' + и" + и) , (v' + v'' + v) и (w' + w'' + w) представляется излишним. Поэтому будем считать, что в Приложении 2 получены формулы, выражающие компоненты скорости (18) в системе Oxyz, через координаты x,y,z.
Для произвольного потенциала внешних массовых сил П = П(ж, у, z, t) давление
р = Р0 - П - Vs/2 (19)
и скорость (18) вместе будут составлять винтовое решение уравнений Эйлера. При этом скорость
V s exp(-4í/Re) (20)
и давление
V S
P = Po exp(-8í /Re) - П (21)
вместе будут составлять винтовое решение уравнений Навье-Стокса.
Таким образом, формулы Приложения 2 вместе с формулами (18) и (19) дают новое винтовое решение уравнений Эйлера. А вместе с формулами (18), (20) и (21) - новое винтовое решение уравнений Навье-Стокса. Это новое решение не периодично. У него нет оси симметрии, нет плоскости симметрии и нет центра симметрии.
6. Заключение
Показано, что линейная комбинация скоростей нестационарных решений уравнений Навье-Стокса, полученных методом Тркала, с одинаковым коэффициентом к в уравнении (3) является скоростью некоторого нестационарного решения уравнений Навье-Стокса. Такой метод получения новых винтовых решений предложено называть методом векторного суммирования. В качестве примера приведено новое точное решение, представляющее собой сумму двух решений Громеки-Бельтрами с непараллельными и непересекающимися осями симметрии и одного «сдвинутого и повернутого» ABC-решения. Методом Тркала (с помощью формул (20) и (21)) это новое стационарное винтовое решение Эйлера легко пересчитывается в новое нестационарное решение уравнений Навье-Стокса. Это решение не периодично. У него нет оси симметрии, нет плоскости симметрии и нет центра симметрии. То есть представлено решение, обладающее качественно новой и более сложной картиной течения.
Предложенный метод значительно обогащает запас точных решений уравнений Эйлера и точных решений уравнений Навье-Стокса со сложной топологией течения. Кроме того,
некоторые новые решения могут использоваться для описания неосесимметричных волн произвольной амплитуды в рамках недиссипативных уравнений магнитной гидродинамики [20].
Такие решения пригодны не только для теоретических исследований, но и для проверки методов численного интегрирования уравнений сплошной среды.
7. Приложение 1. Представление решений в прямоугольной декартовой системе координат
Основная система координат и матрица перехода. Пусть 1^к - базисные векторы системы Охуг. И пусть эти векторы имеют такую взаимную ориентацию, что выполнено векторное равенство к = [1 х ^. Эту систему координат будем называть основной системой координат или, для краткости, основной системой. В точке О, которая в основной системе Охуг имеет координаты (х0,у0, г0), расположим другую прямоугольную декартову систему координат Охуг. Оси этой системы координат могут быть не параллельны осям системы Охуг. Их направления определяются матрицей перехода:
5 =
Sil S12 S13 «21 S22 S23 «31 «32 «33
Столбцы матрицы перехода состоят из компонент единичных базисных векторов 1, ^ к системы Охуг в базисе 1^,к. Например, j = • 1 + «22 • j + «32 • к. В трехмерном случае
семейство таких матриц является трехпараметрическим. С учетом требования к = все возможные варианты матрицы 5 могут быть представлены в виде
i х j
cos ß cos 7 sin a sin ß cos 7 — cos a sin 7 cos a sin ß cos 7 + sin a sin 7 cos a sin 7 sin a sin ß sin 7 + cos a cos 7 cos a sin ß sin 7 — sin a cos 7
(П.1)
— sin ß
sin a cos ß
cos a cos ß
где a, ß, 7 - произвольные константы. (Например, s23 = cos a sin ß sin 7 — sin a cos 7.)
Таким образом, при фиксированной основной системе координат Охуг, задавая числа х0,у0,г0 и a, ß, 7, можно задавать различные системы координат Охуг. Числа х0,у0,г0 определяют положение точки О, а числа а, ß, 7 - направление осей системы Охуг.
Любая точка в пространстве имеет координаты как в основной системе Охуг, так и в системе Охуг. Эти координаты связаны между собой через транспонированную матрицу перехода:
х\ í Su S21 «31 \ / х — хо \
У I = I «12 «22 «32 II У — Уо I . (П.2)
z ) V «13 «23 «33 / V 2 — Z0 J
Перечисленные выше сведения из аналитической геометрии позволяют представлять «сдвинутые и повернутые» известные решения уравнений Эйлера в основной системе координат.
ABC-решение. Представим «сдвинутое и повернутое» ABC-решение в системе координат Охуг. Пусть А, В, С, к - произвольные скалярные константы. Согласно (8), скорость V = u i + V j + w k будет скоростью ABC-решения в системе Охуг, если представить u,v,w в виде
и = Á sin(M) + С cos(ky), v = В sin(fcá) + Á cos(fci), w = С sin(fcy) + В cos(kx).
и \ i 1 t Su S12 S13
V I =I S21 S22 S23
W ) \ S31 S32 S33
Вектор скорости этого ABC-решения V = и i + v j + w k может быть представлен через базисные векторы i,j,k, то есть может быть представлен в виде: V = и i + vj + w k, где
A sin(fcá) + С cos(fcy) \ В s'rn(kx) + Ácos^á) I . (П.3)
С sin(fcy) + В cos(fc;7) )
Последовательное применение формул (П.2) и (П.3) позволяет представить компоненты u,v,w соленоидальной винтовой скорости в базисе i,j,k через координаты x,y,z. Это скорость «сдвинутого» и «повернутого» ABC-решения. Задавая различное расположение системы координат Oxyz, можно получать различные «сдвинутые» и «повернутые» ABC -решения. Хотя по своей физической сущности эти решения являются ABC-решениями, их в общем случае невозможно представить в системе координат Oxyz в «классическом» виде (8). Например, для случая а = 0, [3 = 0, 7 = п/3, х0 = yo = z0 = 0 матрица S будет иметь вид
f 1/2 -^И2 0 5 = I У3И2 1/2 0
V 0 0 1
а компонента скорости и, при значении к = 1, вычисленная по формуле (П.3), будет равна: и = 0.5 Asin(z) + Ccos(£V^2 - У/2) - ^Л3Ё sin(^/2 + у^ЭИ2)-^34cos(z) . Эта компонента и зависит от координаты х, что невозможно для «классического» ABC-решения.
Решение Громеки—Бельтрами. Получим выражения компонент скорости решения Громеки-Бельтрами в основной системе координат Oxyz. Чтобы отличать эти выражения от полученных выше выражений компонент скорости для ABC-решения, будем использовать обозначения со штрихом. (Это значит, в частности, что используемые ниже обозначения, подобные V¿, не являются обозначением производной.) Чтобы задать решение Громеки-Бельтрами, во-первых, нужно выбрать положение цилиндрической системы координат В'г'^'х', которое определяется положением точки В' и единичным вектором ех/, направленным вдоль оси В'. Это можно сделать, задав координаты х'0,у'0,70 точки В' в системе координат Oxyz и задав компоненты d\, d2, d'3 единичного вектора е%/ = d1 ■ i + d'2 ■ j + d'3 ■ k в базисе i,j,k основной системы координат. Во-вторых, нужно выбрать функцию тока ф' = ф'(г,х) одного из решений Громеки-Бельтрами. То есть выбрать некоторую линейную комбинацию функций вида (13), соответствующую некоторому выбранному числу к = 0.
Пусть числа х'о,у'о,z'o, dx, d<2, d3 и функция тока ф' = ф' (г,х), соответствующая числу к = 0, выбраны. Рассмотрим вспомогательную прямоугольную декартову систему координат О'х'77, у которой начало О1 совпадает с началом В' системы В'r'ip'x', а ось О'7 совпадает с осью В'х'. Очевидно, что эти условия однозначно определяют положение оси О'7, но положения осей О'х' и О'у' определяются этими условиями с точностью до поворота вокруг оси О'7. Эта неоднозначность выражается в том, что условию совпадения осей О'7 и В'х' соответствуют различные матрицы перехода от базиса i, j, k системы Oxyz к базису i', j', к' системы О'х'у'7:
s' 11 s' 12 s' 13
/ / /
s 21 s 22 s 23
s' 31 S' 32 s' 33
S' =
Действительно, представим эту матрицу перехода Б' в виде (П.1), но вместо символов а, Р, 7 для обозначения констант будем использовать символы а', ¡3', 7'. Условие совпадения осей О'7 и в'х' означает равенство к' = ех/ = ^ ■ 1 + ■ j + ■ к, которое (согласно определению матрицы перехода) однозначно определяет только коэффициенты третьего
столбца матрицы Б':
сов Ы вш $ осе 7' + вш Ы вш 7' = $ 1, сое а' вш вш — вш о! сое 7' = &2, сов а' сов $ = й' 3.
Поскольку ех' = • 1 + ¿2 • j + ¿3 • к - единичный вектор (то есть (^)2 + (б!2)2 + (^'3)2 = 1), то эта система из трех уравнений имеет бесконечное количество решений. Различным решениям будут соответствовать, вообще говоря, различные положения осей О'' х' и О'у' при заданном положении оси О' х'. Выберем одно из этих решений:
а =
( — агсвш(^'2 ), &3 ^ 0; \ ж + агсвш(^'2 ), $3 < 0;
агсвш
Р' =
(й' 1 />/1 — д!22), <!'3 ^ 0, й'22 = 1; агсвш(^'^\/1 — 22), й'3 < 0, й'22 = 1;
0,
^ 22 = 1;
(П.4)
у = 0.
Соответствующая этому набору констант а', /3', ^ матрица перехода имеет вид
Б' =
5 11 « 12 « 13 ^ 21 ^ 22 ^ 23
31 32 ^ 33
сов вш а' вш сое а' вш 0 сое а' — вш а'
— вш вш а' сое сов а' сов
(П.5)
Транспонированная матрица перехода З'т позволяет вычислить координаты любой точки в системе О'х'у'х' через координаты этой же точки в основной системе Охуг:
х' \ / 11 в' 21 в' 31
У' ] | = | 12 22 32
г' I 1 \ в' 13 8' 23 8' 33
X — X О
У — У' о 2 — г' о
(П.6)
Теперь запишем выражения для компонент скорости одного из решений Громеки-Бельтрами в системе координат О'х'уЧ'. Пусть функции VI = VI (г,х), V|p = У^(г,х), = ^х(г, X) получены по формулам (10) - (12) для выбран-
ной функции тока ф' = ф'(г, х) и выбранного числа к = 0. Тогда скорость V' = VI(г',х') ег/ + У^(г',х') е^' + У^(г',х') ех/ будет скоростью одного из решений Громеки-Бельтрами.
Если выбрано решение Громеки-Бельтрами, не имеющее особенностей на оси симметрии, то скорость V' = VIег/ + е^ + ^ех' в базисе 1', j', к' записывается в виде: V' = и' 1'+у'j' + и:' к'. При этом функции и', у', ь)' вычисляются через координаты х',у',г' следующим образом. Сначала вычисляются координаты г' и :
'х'2 + у'2, х'
г .
(П.7)
Затем функции и', у', ь)' вычисляются по формуле
(х^ ' г (г',х') — у'У' ^(г^х' ))/г' (у'V 'г (г',х') + Х^' <р(г',х'))/г' V' Х^^х') 0 0
V' Х(0,У)
, г' = 0;
г' = 0.
(П.8)
г
Наконец, вектор скорости "V' = VIег/ + У^е^ + ^ех/ = и' 1' + 7j' + го' к' может быть представлен в базисе 1, j, к в виде: "V = и' 1 + V' j + и)' к, где
7 \ / .в' 11 .в' 12 в' / 7 \
7 I = I в'21 в'22 в'23 II 7 I . (П.9)
Ы' ) \ в'31 в'32 8'33 ) \ 7 /
Последовательное применение формул (П.4) - (П.9) позволяет представить компоненты и',ь',и)' соленоидальной винтовой скорости решения Громеки-Бельтрами в базисе 1, j, к через координаты х,у,г.
Интересно отметить одно свойство формул (П.4) - (П.9). Как сказано выше, выбор констант а', ¡3', 7', обеспечивающий совпадение осей О'7 и в'х', неоднозначен. Формулы (П.4) дают лишь один из возможных вариантов. При другом выборе констант а', /3', 7', обеспечивающим совпадение осей О'7 и в'х', матрица перехода (П.5), вообще говоря, изменится. Но окончательный результат последовательного применения формул (П.5) - (П.9) не изменится. Без доказательства заметим, что это свойство следует из наличия осевой симметрии у решения Громеки-Бельтрами.
Геометрический смысл констант а, 7. Иногда ориентацию осей прямоугольной декартовой системы координат Оху7 относительно основной системы Охуг удобно задавать как результат последовательного выполнения некоторых поворотов. В этом случае для вычисления матрицы перехода полезно использовать геометрический смысл констант а, Р, 7, которые являются углами некоторых поворотов. Матрица перехода не зависит от положения точки О. Поэтому, для простоты, рассмотрим случай О = О и покажем, какие именно повороты приводят к системе Охуг.
Первый поворот - это поворот системы Охуг вокруг оси Ог на угол 7. При этом положительным считается угол такого поворота, при котором винт с правой резьбой продвигается в положительном направлении оси Ог. В результате получается промежуточная система Ох1У1г1.
Второй поворот - это поворот промежуточной системы Ох1У1г1 вокруг (новой) оси Оу1 на угол р. При этом положительным считается угол такого поворота, при котором винт с правой резьбой продвигается в положительном направлении оси Оу1. В результате получается вторая промежуточная система Ох2У2%2.
И, наконец, третий поворот - это поворот второй промежуточной системы Ох2У2%2 вокруг оси Ох2 на угол а. При этом положительным считается угол такого поворота, при котором винт с правой резьбой продвигается в положительном направлении оси Ох2.
В результате получается система Охуг. При использовании этого геометрического свойства констант а, 7 следует иметь в виду, что в общем случае приведенный выше порядок выполнения поворотов имеет существенное значение.
8. Приложение 2. Примеры представления точных решений в прямоугольной декартовой системе координат
В пятом разделе данной работы предложено просуммировать скорости трех винтовых решений. Получим выражения компонент скорости для каждого из этих решений в системе Охуг по отдельности.
Первое решение. Для получения выражений компонент скорости первого решения воспользуемся обозначениями Приложения 1. Согласно описанию первого решения, начало цилиндрической системы координат в'г'ф'х' совпадает с началом системы Охуг. То есть все три числа х'0,у'о, г'0, определяющие положение точки в', равны нулю. Если использовать обозначения Приложения 1, то условие ех/ =0-1+0 ^ + 1 ■к, задающее направление оси в'х', означает, что = 0, = 0, в!3 = 1. Отсюда по формулам (П.4) получаем а' = = 7' = 0,
и формула (П.5) дает следующую матрицу перехода:
Б' =
8 11 8 12 8 13 ^ 21 ^ 22 ^ 23 «'31 ^ 32 в' 33
1 0 0 0 1 0 001
(П.10)
Для первого решения выбрана функция тока ф' = ф'(г',х') = г'Л(л/3г') 8т(%'). Она задает скорость V' = VI(г',х') ег' + У^(г',х') е^' + У^(г',х') ех' винтового решения, где функции VI = VI(г',х'), Уф = Ур(г',х'), У^ = Ух(г',^) вычисляются по формулам (14). Например: У^(г',х') = — 231(у/Ъг') 8т(%'). Последовательное применение формул (П.10), (П.6), (П.7), (14), (П.8), (П.9) позволяет представить столбец компонент скорости решения Громеки - Бельтрами, в системе координат Охуг:
и у'
[х31(у/3г1) сов (г) + 2у31(у/3г') 8т(г))/г' {уЫл/З г') сов (г) — 2Х31(Л/ЗГ') 8т(г))/У
г') 8ш(^) 0 0
— л/3 8ш(^)
, г' = у/х2 + у2 = 0;
у/х2 + у2 = 0.
(П.11)
Формула (П.11) позволяет вычислить компоненты и',и'' ,4)' скорости первого решения через переменные х,у,г.
Второе решение. Этому решению соответствует другая цилиндрическая система координат В''г''(р"х''. Далее все обозначения с двумя штрихами будут иметь тот же смысл, какой имеют эти же обозначения с одним штрихом в Приложении 1. По условию, точка В'' в системе Охуг имеет координаты х'0 = 2, у'0 = 0, г'0 = 0, а положительное направление оси В''х" задается единичным вектором ех« = 1/ л/3 ■ 1 + 1/ у/3 ■ j + 1/ л/3 ■ к. Поэтому
4' = 1/у/3, 4 = 1/у/3, 4 = 1/у/3.
Отсюда по формулам (П.4) получаем а" = — агс8т(1/у/3), (П.5) дает следующую матрицу перехода:
[3'' = к/4, 7'' = 0, и формула
Б'' =
„// „// „//
5 11 5 12 5 13
„// „// „//
5 21 5 22 5 23
„// „// „//
5 31 5 32 5 33
1//2 —1//6 1 //3 0 /2//3 1 //3 —1 //2 —1 //6 1 //3
Для второго решения выбрана функция тока ф" = ф"(г'',х") = т"З1(у/Ъг") ехр(—%"). Соответствующие ей функции V" = V"(г'',х"), У'/ = У'/(г'',х"), У" = Ух (г'',х") вычисляются по формулам (15). Матрица перехода Б'' не так проста, как матрица (П.10). Поэтому столбец компонент скорости выбранного решения Громеки-Бельтрами в системе координат Охуг удобнее представить не в виде одной громоздкой формулы, а как результат последовательного применения нескольких формул:
х
11/2 0 —1 //2 = | —1 //6 /2//3 —1 //6
1//3 1 //3 1 //3
ж — 2
У
г
(П.12)
п
а
г'' = у/х1'2 + у''2, = ¿'',
(П.13)
и
w
= <
-x"Ji(y/5r") exp(-x") + 2y"Ji(vV/) exp(-x"))/r" -y"Ji(/5r") exp(-') - 2x" Ji(/5r") exp(-x"))/г" -/5 Jq (/5 r'') exp(-X")
r" = 0;
-л/б exp(-x")
1//2 -1//6 1 //3 o /2//3 i//3 -i//2 -i/V6 i//3
и
го'
(П.14)
(П.15)
Последовательное применение формул (П.12) - (П.15) позволяет вычислить компоненты 7', -и", и)'' скорости второго решения через переменные х,у,г.
Третье решение. Согласно описанию третьего решения, оси системы Оху7 получаются вращением осей системы Охуг вокруг оси Оу на угол к/4 в направлении, показанном на рис. 3. Воспользуемся геометрическим смыслом констант а, 7, описанным в конце Приложения 1. Учитывая правило для определения положительного направления поворота, получаем, что системе Оху7 соответствует последовательность поворотов с углами 7 = 0, @ = —к/4, а = 0. Поэтому матрица перехода, вычисленная по формуле (П.1), имеет вид
5 =
1//2 0 -1//2
0 1 1 / /2 o
0
1 //2
Точка О задается координатами х0 =3, у0 = 3, г0 третьего решения имеет вид
/ 1 //2 0 1 //2
У I = I 0 1 0
7 / \ —1//2 0 1/л/2
= 0. Поэтому формула (П.2) для
х - 3 У - 3 z
(П.16)
В системе координат Охуг третье решение является «классическим» ABC-решением с параметрами Á = 1, В = 2, С = 3, к = 2. Поэтому формула (П.3) имеет вид
и
V I =
W
1//2 0 -1//2
0 1 1 /л/2 0
0
1/V2
sin(2o) + 3cos(2y) 2 sin(2£) + cos(2i) 3 sin(2y) + 2 cos(2á)
(П.17)
Таким образом, компоненты скорости третьего решения в системе координат Охуг вычисляются через переменные х,у,г в результате последовательного применения формул (П.16) и (П.17).
и
0
г
V
Литература
1. Бэтчелор Дж. Введение в динамику жидкости. М.: Мир, 1973.
2. Голубкин В.Н., Марков В.В., Сизых Г.Б. Интегральный инвариант уравнений движения вязкого газа // ПММ. 2015. Т. 79, вып. 6. С. 808-816.
3. Trkal V. Poznamka k hydrodynamice vazkych tekutin // Casopis pro pestovani matematiky a fysiky (Praha). 1919. V. 48, I. 3. P. 302-311.
4. Неменьи П.Ф. Современное развитие обратных и полуобратных методов в механике сплошной среды // Проблемы механики / под ред. Р. Мизеса и Т. М. Кармана. М.: ИЛ, 1955. С. 234-257.
5. Громека И.С. Некоторые случаи движения несжимаемой жидкости // Ученые записки Казанского университета. 1881.
6. Beltrami E. Considerazioni idrodinamiche // Rend. Inst. Lombardo Acad. Sci. Lett. 1889. V. 22. P. 122-131.
7. Громека И.С. Некоторые случаи движения несжимаемой жидкости. Собрание сочинений. М., 1952. С. 76-148.
8. Bogoyavlenskij O.I. Method of symmetry transforms for ideal MHD equilibrium equations. The legacy of the inverse scattering transform in applied mathematics // Contemp. Math. Providence, RJ, Amer. Math. Soc. 2002. V. 301. P. 195-218.
9. Bogoyavlenskij O.I. Exact solutions to the Navier-Stokes equations // Comptes Rendus Math. Acad. Sci. Soc. R. Canada. 2002. V. 24, N 4. P. 138-143.
10. Bogoyavlenskij O.I. Exact solutions to the Navier-Stokes equations and viscous MHD equations // Phys. Lett., Ser. A. 2003. V. 307, N 5-6. P. 281-286.
11. Vereshchagin V.P., Subbotin Y.N., Chernykh N.I. On the mechanics of helical flows in an ideal incompressible nonviscous continuous medium // Proceedings of the Steklov Institute of Mathematics. 2014. V. 284, N 1. P. 159-174.
12. Vereshchagin V.P., Subbotin Y.N., Chernykh N.I. Some solutions of continuum equations for an incompressible viscous medium // Proceedings of the Steklov Institute of Mathematics. 2014. V. 287, N 1. P. 208-223.
13. Ламб Г. Гидродинамика. М.: ОГИЗ. ГИТТЛ, 1947.
14. Arnold V.I. Sur la topologie des ecoulements stationnaires des fluides parfaits // C. R. Acfd. Sci. Paris. 1965. V. 261, N 1. P. 17-20.
15. Ковалёв В.П., Сизых Г.Б. Осесимметричные винтовые течения идеальной жидкости // ТРУДЫ МФТИ. 2016. Т. 8, № 3. С. 171-179.
16. Сизых Г.Б. Эволюция завихренности в закрученных осесимметричных течениях вязкой несжимаемой жидкости // Ученые записки ЦАГИ. 2015. Т. 46, № 3. С. 14-20.
17. Васильев О.Ф. Основы механики винтовых и циркуляционных потоков. М.: Госэнергоиздат, 1958.
18. Jonas P. Trkalian flow fields // Institute of Thermomechanics AS CR, v.v.i., Praha. 2008.
19. Lakhtakia A., Trkal V. Beltrami fields and Trkalian flows // Czech. J. Phys. 1994. V. 44, I. 2. P. 89-96.
20. Dritschel D.G. Generalized helical Beltrami flows in hydrodynamics and magnetohydrodynamics // Journal of Fluid Mechanics. 2006. V. 222. P. 525-541.
References
1. Batchelor G.K. An Introduction to Fluid Dynamics. Moscow: Mir, 1973. (in Russian).
2. Golubkin V.N., Markov V.V., Sizykh G.B. The integral invariant of the equations of motion of a viscous gas. J. Appl. Math. Mech. 2015. V. 79, I. 6. P. 808-816. (in Russian).
3. Trkal V. A note on the hydrodynamics of viscous fluids (translated by I. Gregora). Czech. J. Phys. 1994. V. 44. P. 97-106.
4. Nemenyi P.F. Recent developments in inverse and semi-Inverse methods in the mechanics of Continua. Moscow: IL, 1955. P. 234-257. (in Russian).
5. Gromeka I.S. Some cases of the motion of an incompressible fluid. Uchenie zapiski Kazanskogo universiteta. 1881. (in Russian).
6. Beltrami E. Considerazioni idrodinamiche. Rend. Inst. Lombardo Acad. Sci. Lett. 1889. V. 22. P. 122-131.
7. Gromeka I.S. Some cases of motion of incompressible fluid. Moscow: Akad. nauk, 1952. P. 76-148. (in Russian).
8. Bogoyavlenskij O.I. Method of symmetry transforms for ideal MHD equilibrium equations. The legacy of the inverse scattering transform in applied mathematics. Contemp. Math. Providence, RJ, Amer. Math. Soc. 2002. V. 301. P. 195-218.
9. Bogoyavlenskij O.I. Exact solutions to the Navier-Stokes equations. Comptes Rendus Math. Acad. Sci. Soc. R. Canada. 2002. V. 24, N 4. P. 138-143.
10. Bogoyavlenskij O.I. Exact solutions to the Navier-Stokes equations and viscous MHD equations. Phys. Lett. Ser. A. 2003. V. 307, N 5-6. P. 281-286.
11. Vereshchagin V.P., Subbotin Y.N., Chernykh N.I. On the mechanics of helical flows in an ideal incompressible nonviscous continuous medium. Proceedings of the Steklov Institute of Mathematics. 2014. V. 284, N 1. P. 159-174.
12. Vereshchagin V.P., Subbotin Y.N., Chernykh N.I. Some solutions of continuum equations for an incompressible viscous medium. Proceedings of the Steklov Institute of Mathematics. 2014. V. 287, N 1. P. 208-223.
13. Lamb H. Hydrodynamics. Moscow: OGIZ GITTL, 1947.
14. Arnold V.I. Sur la topologie des ecoulements stationnaires des fluides parfaits. C. R. Acfd. Sci. Paris. 1965. V. 261, N 1. P. 17-20.
15. Kovalev V.P., Sizykh G.B. Corkscrew axisymmetric flow of an ideal fluid. Proceedings of MIPT. 2016. V. 8, N 3. P. 171-179. (in Russian).
16. 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).
17. Vasil'ev O.F. Fundamentals of mechanics of screw and flow patterns. Moscow, 1958. (in Russian).
18. Jonas P. Trkalian flow fields. Colloquium Fluid Dynamics 2008. Institute of Thermomechanics AS CR. v.v.i. Prague. October 22-24, 2008.
19. Lakhtakia A., Trkal V. Beltrami fields and Trkalian flows. Czech. J. Phys. 1994. V. 44, I. 2. P. 89-96.
20. Dritschel D. G. Generalized helical Beltrami flows in hydrodynamics and magnetohydrodynamics. Journal of Fluid Mechanics. 2006. V. 222. P. 525-541.
Поступила в редакцию 06.03.2017
