12
УДК 532.546
А. А. Юрова, А. В. Юров, И. В. Лукиных
ОБ ОДНОМ МЕТОДЕ ПОСТРОЕНИЯ ТОЧНЫХ РЕШЕНИЙ
УРАВНЕНИЙ ДВУМЕРНОЙ ГИДРОДИНАМИКИ НЕСЖИМАЕМОЙ ЖИДКОСТИ
Предложен простой алгебраический метод построения точных решений уравнений двумерной гидродинамики несжимаемой жидкости.
В случае невязкой жидкости задача сводится к последовательному решению трех линейных дифференциальных уравнений в частных производных, а в случае вязкой — к решению трех линейных дифференциальных уравнений в частных производных и одного обыкновенного дифференциального уравнения первого порядка.
We propose a simple algebraic method for constructing exact solutions of equations of two-dimensional hydrodynamics of an incompressible fluids.
The problem reduces to consecutively solution three linear partial differential equations for a nonviscous fluid and to solving three linear partial differential equations and one first-order ordinary differential equation for a viscous fluid.
Ключевые слова: двумерная несжимаемая жидкость, точные решения.
Key words: two-dimensional incompressible fluid, exact solutions.
Динамика несжимаемых вязких потоков жидкости описывается уравнениями Навье — Стокса (НС). В режиме больших чисел Рейнольдса возникает турбулентность — одна из важнейших нерешенных задач теоретической и математической физики. Результаты Като показывают, что двумерные (2D) уравнения НС глобально определены в C0 ([О, да ] Hs (R2)) при s > 2 и при 0 < T < да слабое решение уравнений 2DHC стремится к решению 20-уравнения Эйлера в C0 (0, T] Hs (R2 ))[l]. В свою очередь, уравнение 3DHC локально определено в C0 ([0,т]; Hs (r 2)), при s > 5/2 и аналогично предыдущему его слабые решения приближаются к решениями ЗО-уравнения Эйлера C0 (М; Hs (R3)), где Т определяется начальными данными (нормой) и внешними силами [2; 3]. Невязкий предел интенсивно исследовался в работах [4; 5]. Среди важнейших результатов можно отметить факт гамильтоновости 2D -уравнений Эйлера, доказанный Арнольдом [6], и результаты исследования симплектической структуры этих уравнений [7]. В работах [8] и [9] было найдено представление пары Лакса для 2D -уравнения Эйлера, записанное в эйлеровых переменных. Позднее было построено представление Лакса и в евклидовых переменных [10].
© Юрова А. А., Юров А. В., Лукиных И. В., 2013
Вестник Балтийского федерального университета им. И. Канта. 201З. Вып. 4. С. 12 — 17.
В настоящей работе мы описываем удивительно простой, но эффективный метод построения точных решений уравнений двумерной гидродинамики несжимаемой жидкости, применимый как к невязкой (уравнения Эйлера), так и вязкой жидкости (уравнения НС).
Рассмотрим двумерную невязкую жидкость. Две компоненты скорости Ух и Уу выражаются через функцию тока / = /(1,, х, у) по формулам
*х =/• =/ • (!)
ду дх
так что выполняется уравнение неразрывности й\уу = 0 . В этих переменных двумерное уравнение Эйлера принимает вид
дД/д/дД/д/дД/о ^
дt ду дх дх ду
где Д означает двумерный лапласиан [11]. Равенство (2) представляет собой нелинейное уравнение, которое не относится к числу так называемых интегрируемых (несмотря на наличие представления Лакса). Тем не менее для него можно развить процедуру построения точных решений.
Теорема. Пусть / — гармоническая в области & функция: Д / = 0, где Д — двумерный лапласиан. Пусть теперь F = , X, у) —
решение переопределенной системы линейных дифференциальных уравнений
дГ д / дГ д / дГ й 1п к
ДГ = кГ, ---= —-----------------------Г, (3)
дt дх ду ду дх &
где к ^) — некоторая функция времени. Тогда функция
/ = / + Г (4)
удовлетворяет в & уравнению (2), т. е.
дД/ д/1 дД/1 д/ дД/ о
дt ду ох дх ду
Справедливость теоремы проверяется прямым вычислением. В дальнейшем будем ссылаться на выражения (3) и (4) как процедуру «одевания» & = 2 уравнений Эйлера. Преобразование (4) похоже на преобразование Дарбу (ПД), применяемое в теории интегрируемых систем. Суть ПД заключается в нахождении решения пары Лакса / с заданным «затравочным» потенциалом (который, в свою очередь, является решением исследуемого нелинейного уравнения) и последующем использовании этого / для нахождения новых потенциалов.
13
14
Сродство ПД с описанным выше преобразованием очевидно. В самом деле, в качестве промежуточного шага необходимо решать два линейных уравнения с переменными коэффициентами (3), зависящими от гармонической функции /, которая может считаться «затравочным» потенциалом, поскольку удовлетворяет уравнению (2) и описывает плоское (стационарное, если / не зависит от ^) потенциальное течение. Новая функция тока /1 (4) будет уже описывать нестационарное течение с завихренностью.
Тем не менее это не преобразование Дарбу, поскольку линейная система (3) не является [Ь, А] — парой уравнения Эйлера (2). Условие совместимости уравнений (3) имеет виц
Ск
Ав- кв = — Г,
(5)
ду дГ ду дГ
где в =------------------. Используя (3) и факт гармоничности у,
дх ду ду дх
можно переписать условие (5) в другом виде:
С д2щ д2у| д2Г і С д2Г д2Г| д2у V
ч дх2 ду2 ) дхду
ду дх ) дхду 2
1 Ск
-Г.
(6)
Замечание. В систему (3) входит вспомогательная функция времени к = к (ї). Ее величина не фиксирована, а определяется формой гармонической функции у. В частности, не для всякой гармонической у система (3) совместна, т. е. не для всякой у существует к = к (ї) такая, что система (3) имеет решение. Сейчас мы покажем, этот элементарный подход неожиданно оказывается достаточно эффективным для построения точных решений двумерного уравнения Эйлера.
Рассмотрим гармоническую функцию в виде уУ = ХуЛуУу, где ЛУ — постоянная ( +1) х (У +1) матрица с элементами аік, а і и к пробегают от нуля до N . Вид матрицы ЛУ определяется путем подстановки уУ в уравнение А уУ = 0 . Ниже мы приводим несколько примеров матриц при различных N :
Л =
V а02
0а О о а011 , Л2 =
V а10 а11)
а01 а I 02
а11 0
0 0 ,
A =
a о о 1 a a02 a30^ aoo a01 a02 a30 a04
a10 a11 -3a30 a13 a1o a11 -3a30 a13 0
0 a - -3a03 0 0 , A4 = -a02 -3a03 -6"o4 0 0
0 J a - 1J 0 0 у a30 -a13 0 0 0
. a04 0 0 0 0 ,
Число свободных параметров M (N) (т. е. число
независимых эле-
ментов в матрицах AN) определяется следующим образом: если N нечетное число, то M(N) = 2(N +1), а если N четное число, то M (N) = 2 N +1. Такая функция тока служит лишь вспомогательным средством для вычисления поля скоростей. Поскольку после дифференцирования по пространственным переменным коэффициент a00 обращается в нуль, не теряя общности, можно положить a00 = 0, что мы и будем подразумевать в дальнейшем. Остальные коэффициенты следует рассматривать как функции времени a. = a. (t). Подставляя /N в систему дифференциальных уравнений (3), находим функцию F, после чего по формуле (4) вычисляем соответствующую / N.
В качестве простого примера положим a11 = const и a02 = const. Производя вышеуказанные вычисления, находим
a
a (t) = a0 cosat, b (t) =-— (a11 cosh at + a sinh at),
2a„,
15
где a0 — постоянная интегрирования, а a = яЦ + 4^ . Что же до c(t),
то удобно считать a10 (t)b(t) = a01 (t)a(t), что дает c = const, причем
можно без потери общности, выбрать c = 0 . Отметим, что решение (10) не сингулярно, а k (t) = -a2 (t) - b2 (t). Тем же самым методом можно
строить и решения, описывающие потенциальные течения.
Вышеописанная методика легко обобщается на случай двумерного течения несжимаемой вязкой жидкости, описываемой уравнением
дЛ/ д/ дЛ/ д/ дЛ/ 2
' = vЛ /,
дt
ду дх дх ду
(8)
где V — кинематическая вязкость. Пусть / — гармоническая в области П функция: А/ = 0. Пусть теперь ¥(, X, у) — решение переопределенной системы нелинейных дифференциальных уравнений
, ч дF д/ дF д/ дF , ч
ЛF = u (t)F, — +—-------------------------+ U(u)F,
дt дх ду ду дх
(9)
где u (t) удовлетворяет обыкновенному дифференциальному уравнению
1б
дu (t)
+ uU (u ) = vu2, (lO)
дt
а U (u) — произвольная функция от своего аргумента. Тогда функция /1 =/ + F удовлетворяет в D уравнению (8):
дЛ/ д/ дЛ/ д/ дЛ/ 2
----L +---------------------L-L-1-L = vЛ /
дt ду дx Dx ду
В качестве простейшего примера рассмотрим «одевание» на нулевом фоне / = 0 . Решая систему (9), находим
/ = Ke к(а + ^ (C1 sin (ax + Ьу) + C2 cos (ax + Ьу)),
где K,a,b,C12 — произвольные постоянные. Разумеется, это простое и
известное решение. Мы привели его только для демонстрации работоспособности метода.
Теперь рассмотрим случай N = 1 . По аналогии с уравнением (8) можно построить решение (9) вида
/1 = A(t) (у sina(t) + xcosa(t)) + # (t) sin g + #2 (t) cos g, (ll)
где g = R(xcosф + ysinф + f(t)), U(u) = vu; u = -R2 =const; A(t),a(t),f(t) — произвольные функции; ф — произвольная постоянная, а функции #1 (t) , #2 (t) ищутся из системы обыкновенных дифференциальных уравнений
d#d() =-vR2#1 + R - A (t)sin (-a(t))#
= -vR 2#2 - R f ^ - A (t )sin (--a(t ))V
dt
V
Разумеется, подобно уравнению (8) можно построить суперпозицию произвольного числа синусов и косинусов вместо выражения (11). Это наблюдение является следствием линейности уравнений (9) и линейности преобразуемой функции у/ по пространственным переменным. В частном случае df (t)/dt = A(t)sin( ф -a (t)) находим
/ = A(t) (y sina(t) + x cos a (t)) + e-VR2‘ (l (0) Sin g + £ (0) cos g) ,
где #12 — постоянные интегрирования.
Сравнивая уравнения (8) и (9), видим, что учет вязкости приводит, как следовало ожидать, к появлению дополнительного экспоненциального множителя описывающего диссипацию.
Таким образом, уравнения двумерной гидродинамики несжимаемой жидкости допускают простой алгебраический метод построения
точных решений. В данной работе мы ограничились демонстрацией простейших решений отдельно для случаев невязкой и вязкой жидкости. Не вызывает никаких сомнений то, что с помощью вышеописанной методики можно построить множество значительно более сложных и интересных с физической точки зрения решений.
Список литературы
1. Kato T. Global solutions of two-dimensional Navier-Stokes and Euler equations // Proc. Symp. Pure Math. 1986. Pt. 2. Vol. 45.
2. Kato T. Quasi-linear equations of evolution, with applications to partial differential equastions // Journal of Functional Analysis. 1972. Vol. 9. P. 296.
3. Kato T. Spectral theory and differential equastions // Lecture Notes in Mathematics / ed. W.N. Everitt. Berlin, 1975. Vol. 448. P. 25.
4. Constantin P., Wu J. The inviscid limit for non-smooth vortivity // Indiana University Mathematics Journal. 1996. № 1. P. 45.
5. Wu J. J. The Inviscid Limit of the Complex Ginzburg-Landau Equastion. //Differential Equestions. 1998. Vol. 142. P. 413.
6. Arnold V. I. Sur la Geometrie Differentielle des Groupes de Lie de Dimension Infinie et ses Applications a Lalhydrodinamique des Fluides Parfaits / / Ann. Inst. Fourier. Grenoble, 1966. Vol. 16. P. 319.
7. Marsden J. E. Lectures on mechanics. Navier - Stokes equastions: a mathematical analysis //Lond Math. Soc. Lect. Note. Ser. 174. Cambridge Univ. Press, 1992.
8. Fridlander S., Vishik M. Lax pair formulation for the Euler equation. // Phys. Lett. 1990. A. 148(6—7). P. 313—319.
9. Fridlander S., Vishik M. An inverse scattering treatment for the flow of an ideal fluid in two dimensions // Nonlinearity. 1993. Vol. 6. P. 231.
10. Li Y., Yurov A. Lax pairs and Darboux transformations for Euler equastions. // Studies In Applied Math. 2003. Vol. 111. P. 101—103.
11. Ландау Л. Д., Лифшиц Е. М. Теоретическая физика. М., 1986. Т. 6 : Гидродинамика.
12. Итс А. Р., Рыбин А. В., Салль М. А. К вопросу о точном интегрировании нелинейного уравнения Шредингера // ТМФ. 1988. № 1. C. 20, 74.
17
Об авторах
Алла Александровна Юрова — канд. физ.-мат. наук, доц., Калининградский государственный технический университет.
E-mail: yurov@freemail.ru
Артем Валерьянович Юров — д-р физ.-мат. наук, проф., Балтийский федеральный университета им. И. Канта, Калининград.
E-mail: artyom_yurov@mail.ru
Ирина Викторовна Лукиных — асп., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: luki777@mail.ru
About authors
Alla Yurova — PhD, ass. prof., Kaliningrad State Tehnical University.
E-mail: yurov@freemail.ru
Dr Artyom Yurov — prof., I. Kant Baltic Federal University, Kaliningrad.
E-mail: artyom_yurov@mail.ru
Irina Lukinyh — PHD stud., I. Kant Baltic Federal University, Kaliningrad. E-mail: luki777@mail.ru