ПРЕОБРАЗОВАНИЯ ЭКВИВАЛЕНТНОСТИ УРАВНЕНИЙ ДВИЖЕНИЯ ДВУМЕРНОГО СЛОЯ ИДЕАЛЬНОЙ ЖИДКОСТИ

Полный текст

Описание системы уравнений. Многие исследователи нелинейных уравнений давно используют групповой анализ дифференциальных уравнений. Так как принципы инвариантности закладываются при выводе уравнений, применение группового анализа особенно эффективно и плодотворно в области фундаментальной механики и физики [1-4]. В последнее время особо актуальным стало исследование задач математического моделирования разнообразных физических процессов. Главной целью группового анализа является построение точных решений дифференциальных уравнений. На пути построения этих решений важной частью задачи является нахождение преобразований эквивалентности для заданной системы уравнений. В данной работе получены преобразования эквивалентности для системы уравнений, описывающей двумерное движение иде- альной жидкости, записанных в модифицированных переменных. Рассмотрим систему уравнений движения двумерного слоя идеальной жидкости по наклонной плоскости. Уравнения запишем в декартовой системе координат так, чтобы ось z была ортогональна к подложке, а ось x направлена в сторону действия скатывающей силы. Жидкость занимает область Ω = {( x, z): -ж < x < +ж, 0 < z < H ( x, t )} , где t - время; H - толщина слоя жидкости. В уравнения движения жидкости входят компоненты вектора скорости (u, w), p - давление (скатывающие силы заменой переменных можно включить в давление). Слой жидкости имеет твердую подложку при z = 0 и свободную границу при z = H ( x, t ) (см. рисунок). 46 Математика, механика, информатика -h^] = H 2[-px +^ξ ], (1) Решение определяющих уравнений. Поставим задачу группового анализа для уравнений (1)-(4). Необходимо найти преобразования эквивалентности для этой системы уравнений. Преобразования эквивалентности - это преобразования, которые сохраняют структуру исходной системы уравнений [1]. Введем индексные обозначения: x1 = x, x2 =ξ , x3 = t , u1 = u , u2 = w . u = p, u = H . Будем считать, что если f(x1,x2,x3) - некоторая функция, то fr = Cf fj = C2 f i, j = 1, 2,3. Движение двумерного слоя жидкости толщины H по наклонной плоскости Рассмотрим систему уравнений движения двумерного слоя идеальной жидкости по наклонной плоскости. Уравнения запишем в декартовой системе координат так, чтобы ось z была ортогональна к подложке, а ось x направлена в сторону действия скатывающей силы. Жидкость занимает область Ω = {(x, z) : -ж < x < +ж, 0 < z < H(x,t)} , где t - время; H - толщина слоя жидкости. В уравнения движения жидкости входят компоненты вектора скорости (u,w), p - давление (скатывающие силы заменой переменных можно включить в давление). Слой жидкости имеет твердую подложку при z = 0 и свободную границу при z = H ( x, t ) (см. рисунок). В дальнейшем вводятся новые координаты — — z t = t, x = x, ξ =-, H ( x, t) модифицированные компоненты вектора скорости и давление [5]: u = uH, w = w - uξHX , p = p. Это позволяет фиксировать границы области по ξ , Ω = {(x,ξ) : -ж < x < +ж, 0 < ξ < 1}. Система уравнений в новых переменных записывается в виде (черту над переменными опускаем) ρ[Hut - Ht (u + ξUξ ) + uux + wurq dx1 dx1 dxJ Система уравнений (1)-(4) перепишется в индексных обозначениях так: T7 1 и / 1 , 2 1ч. 11. 21 Hu3 - H3(u + xu2) + u u1 + uu2 - -—H1 (u1)2 + H2u3 - x2HH1u3 = 0, (5) H Hu32 + x2(u1 H1 + uH13 -—u1 H1H3) -H —x H^u2 +--H1 (u + x u2)) + H +u щ + u x (— Hu1 +--H11 u--—) + 1 H 11 H 11 (H1)2 +u (u~2 +--H-[(u + x ^)) + u~2 — 0, H (6) (7) (8) Инфинитезимальный оператор, допускаемый системой (5)-(8), ищем в виде [1] u\ + u2 — 0, u24 = 0. „ w д α C X = ξ1-+ ηα-, dxi Cua ' (9) где i, α = 1, 2, 3, 4, а координаты оператора ξ1, ηα зависят от всех зависимых и независимых переменных (по всем повторяющимся индексам производится суммирование). Поскольку в систему уравнений входят первые производные, то для формирования определяющих уравнений нужно продолжить оператор X на первые производные: p[H(w + S)t -Z1Ht(w + Ξ)ξ + u(w + S)x + w(w+Ξ)ξ] = -p, (2) ux + wZ= 0, (3) Ηξ= 0, (4) где ρ - плотность жидкости (ее можно положить равной единице); S = (Zu1Hx)/H ; нижние индексы обозначают дифференцирование по x, ξ , t функций u , w , p . Функция H(x, /) и ее производные существенно входят в уравнения (1)-(4), которые образуют замкнутую систему относительно u , w , p . X=X+cf due ζβ =d^ + u?^L· - uß dx1 дηß dxn - uj δξ_ dx1 • + ui δξ_ Cun где β, l = 1, 2, 3 (суммирование по повторяющимся индексам), координаты ζβ зависят от всех x', uα и производных u“. Далее, действуя оператором X на уравнения (8), (7), (6), (5), из критерия инвариантности [1] получаем определяющие уравнения, которые расщепляем по независимым переменным. При этом переходим на многообразие, задаваемое исходными 47 Вестник СибГАУ. № 4(50). 2013 уравнениями, т. е. u\, u^ , u^ выражаем через оставшиеся элементы и подставляем в определяющие уравнения. Используем также четвертое уравнение u4 = 0. После достаточно трудоемких выкладок при расщеплении определяющих уравнений относительно независимых переменных получаем следующее представление координат оператора (9): ξ1 = B1 x1 + C1, ξ2 = (B1 - A4)x2, ξ3 = B3 x3 + C3, η1 = (B + A4 - B3)u', η2 = (B - B3 )u 2, η3 = ( B1 - B3 )u 3 + D( x3), η4 = A4u4 , где A1, A4, B1, B3, C1, C3 - постоянные; D(x3) -произвольная функция. Выпишем операторы преобразования эквивалентности: X =A 1 Cx1 ’ X = — 2 Cx3 ’ 1 д 2 д 1 д 2 д ,-3 д X3 = x —■—+ x —— + u —■—+ u —— + 2u —~ , Cx Cx Cu Cu Cu X 4 = x^^ - u1 A - u 2^- - 2u 3^-, Cx Cu Cu Cu v 2 д 1 д 4 д X 5 —— x —T + u ——+ u -— , Cx Cu Cu X6(D(x3)) = D(x3)-^ . Cu На основании теоремы Ли для найденных операторов ставится во взаимно однозначное соответствие группа преобразований переменных x1 , x2 , x3 , u1 , u2, u3, h. Тем самым доказали, что преобразования эквивалентности для системы (5)-(8) имеют структуру бесконечномерной группы преобразований: G3ks :{x1 = b3b4x1 + b ; x2 = b3b4b5-1x2 ; 3 3 1 122 x = b4x + b2; u = b3b5u ; u = b3u ; u3 = b32u3 +T(x3); H = b5 H}, (10) где b1, b2, b3, b4, b5 - групповые параметры (10), Ψ( x3) - произвольная функция. Заметим, что для H допускается только растяжение с одновре- 21 менным растяжением для x и u . Автор благодарит своего научного руководителя В. К. Андреева за постановку задачи.

Об авторах

Дарья Александровна Краснова

Институт вычислительного моделирования СО РАН

Email: krasnova-d@mail.ru
аспирант Российская Федерация, 660036, Красноярск, Академгородок, 50/44

Список литературы

Дополнительные файлы