Слоистые градиентные стационарные течения вертикально завихренной вязкой несжимаемой жидкости Е.Ю. Просвиряков evgen_pros@mail.ru Уральский государственный университет путей сообщения (Екатеринбург) Институт машиноведения УрО РАН (Екатеринбург) Аннотация Найдено точное решение, описывающее стационарное течение за- вихренной вязкой несжимаемой жидкости. Приведенное решение системы Навье-Стокса имеет характерную особенность при опи- сании движущейся жидкости — справедливо при учете не только вязких сил, но и сил инерции. Учет сил инерции приводит к по- явлению застойных точек в слое жидкости, противотечений и су- ществованию толщин слоя, при которых касательные напряжения на нижней границе обращаются в нуль. Показана роль давления в формировании числа застойных точек в слое жидкости. Ключевые слова: точное решение; противотечение; течение Пу- азёйля; слоистое течение; застойная точка. 1 Введение При изучении крупномасштабных течений в изотермических уравнениях Навье-Стокса и их конвективных аналогах используется гидростатическое приближение. Применение «наивных» уравнений океанических течений является спорным ввиду их заведомой приближенности. Использование уравнений термодиффу- зии при моделировании течений жидкости в океане основано на удержании вертикальной скорости во всех уравнениях движения кроме уравнения сохранения импульсов на ось, относительно которой предлагается использовать гидростатическое приближение. Применение слоистых течений для описания крупномас- штабных течений затруднительно, поскольку для неизотермического океана система Обербека-Буссинеска является переопределенной. Очевидно, что выходом из этой ситуации может служить интегрирование пол- ной системы уравнений Навье-Стокса, но это не облегчит, если не усложнит описание движения жидкости на первых порах. Наиболее перспективным представляется путь, при котором будут найдены точные ре- шения для нелинейных уравнений движения жидкости, описывающих слоистые течения, а затем будет проведено сравнение с классическими («наивными») уравнениями теории океана и с полной системой уравнений Навье-Стокса. Для этой цели можно использовать известный класс точных решений, в котором скорости линейны по части координат [1, 2, 3, 4, 5, 6]. К первым точным решениям указанного класса движений относятся классические течения Куэтта [7] и Экмана [8]. Течение Куэтта [7] описывает однонаправленное слоистое течение вязкой несжимаемой жидко- c by the paper’s authors. Copying permitted for private and academic purposes. Copyright ° In: G.A. Timofeeva, A.V. Martynenko (eds.): Proceedings of 3rd Russian Conference "Mathematical Modeling and Information Technologies" (MMIT 2016), Yekaterinburg, Russia, 16-Nov-2016, published at http://ceur-ws.org 164 сти при отсутствии вращения и постоянном давлении. В то время, как решение Экмана [8] справедливо для описания противотечений вращающихся масс жидкости. Взяв за основу эти два классических точных реше- ния, Стоммел [9] и Чарни [10] предложили математические модели, которые описывали бы экваториальные противотечения. Как было показано Стоммелом [9], использование линейных определяющих соотношений (линеаризованных уравнений Навье-Стокса) не пригодно для описания движений жидкости на экваторе. Точное решение Чарни [10], качественно описывающее нелинейные эффекты генерации противотечений, имеет очень узкую область применимости. Этот недостаток решения Чарни можно объяснить сильным сглаживанием нелинейных слагаемых уравнений Навье-Стокса и отсутствием вертикальной компоненты завихренности в модельном течении. Косвенным подтверждением этих результатов можно указать непри- годность использования точного решения Куэтта для описания противотечений. Противотечения в слое жидкости для течения Куэтта возможны при учете скольжения жидкости на нижней границе. В этом случае противотечения могут наблюдаться только для несжимаемых реологических сред. В статьях [11, 12, 13] были предложены математические модели, которые посредством точных реше- ний уравнений Навье-Стокса объясняют генерацию противотечений в жидкости и усиление скоростей по сравнению с заданными граничными возмущениями. Характерной особенностью полученных результатов в [11, 12, 13] является демонстрация существования противотечений в слоистых течениях при постоян- ном давлении. Было показано, что существование застойных точек и противотечений возможно благодаря наличию вертикального вихря в жидкости. Далее интерес представляют слоистые течения, у которых дав- ление будет влиять на количество застойных точек. В данной статье дается ответ на этот вопрос в рамках одного точного решения уравнений Навье-Стокса. 2 Постановка задачи Уравнения Навье-Стокса, описывающие слоистые течения вязкой несжимаемой жидкости, в прямоуголь- ной декартовой системе координат Oxyz записываются следующим образом [2]: µ 2 ¶ ∂Vx ∂Vx ∂Vx ∂P ∂ Vx ∂ 2 Vx ∂ 2 Vx + Vx + Vy =− +ν + + , ∂t ∂x ∂y ∂x ∂x2 ∂y 2 ∂z 2 µ ¶ ∂Vy ∂Vy ∂Vy ∂P ∂ 2 Vy ∂ 2 Vy ∂ 2 Vy + Vx + Vy =− +ν + + , ∂t ∂x ∂y ∂y ∂x2 ∂y 2 ∂z 2 ∂P = 0, ∂z ∂Vx ∂Vy + = 0. (1) ∂x ∂y p Здесь (Vx ; Vy ; 0) — вектор скорости слоистого течения жидкости, P = — давление p, деленное на посто- ρ янную плотность жидкости ρ; ν — коэффициент кинематической (молекулярной) вязкости. Система нелинейных уравнений в частных производных (1) является переопределенной, поскольку для определения трех неизвестных функций Vx , Vy , P в ней содержится четыре уравнения. Первые три урав- нения в (1) определяют баланс сил инерции, градиента приведенного давления P, имеющего размерность энергии, и внутреннего трения в жидкости. Отметим, что для уравнения импульсов, спроецированного на ось Oz, выполняется условие гидростатического равновесия. Для системы (1) при постоянном давлении (P = const) в статье [11] Vx = U (z, t) + u1 (z, t) x + u2 (z, t) y, Vy = V (z, t) + v1 (z, t) x + v2 (z, t) y. (2) Поле скоростей (2) линейно по координатам x и y, которые далее будем называть горизонтальными (про- дольными). Соответственно координату z будем называть вертикальной (поперечной). Ограничимся далее исследованием стационарного точного решения специального вида (2) системы (1). Рассмотрим следующую структуру поля скоростей: Vx = U (z) + u2 (z) y, 165 Vy = V (z) . (3) Очевидно, что при описании стационарных движений вязкой несжимаемой жидкости неизвестные функции U, V, u2 определяются только вертикальной координатой. В дальнейшем аргумент z, от которого зависят скорости U, V и градиент u2 скорости Vx указываться не будет. Решение (3) будем использовать при описании движения вязкой несжимаемой в бесконечно протяжен- ном слое в горизонтальном направлении, определяемом декартовым произведением множеств (x; y; z) = (−∞; +∞) × (−∞; +∞) × [0; h] . Уравнения z = 0 и z = h определяют границы слоя в поперечном направлении, следовательно, h — толщина слоя жидкости. Подставим выражения (3) в уравнения (1), получим систему уравнений вида µ 2 ¶ ∂P ∂ U ∂ 2 u2 V u2 = − +ν + y , ∂x ∂z 2 ∂z 2 ∂P ∂2V =ν 2, ∂y ∂z ∂P = 0. (4) ∂z Исследуя структуру последней системы уравнений, получим, что давление описывается линейной формой P = P0 (z) + xP1 + yP2 = S + g (h − z) + xP1 + yP2 . (5) В выражении (5) S, P1 и P2 — постоянные величины; g — ускорение свободного падения; h — толщина слоя. Подставляя в систему (4) выражение, описывающее давление (5), получим, приравненивая коэффици- енты при одинаковой степени координаты y, уравнения для определения гидродинамических полей: d2 u2 d2 V 2 = 0, ν 2 = P2 , dz dz d2 U ν = V u2 + P1 . (6) dz 2 3 Вычисление точного решения Последовательное интегрирование уравнений (6) позволяет получить общее решение этой системы z2 u2 = c1 z + c2 , V = P2 + c3 z + c4 , 2ν · ¸ 1 z5 z4 z3 z2 P1 z 2 U= c1 P2 + (2c1 c3 + c2 P2 ) + (c1 c4 + c2 c3 ) + c2 c4 + + c5 z + c6 . (7) ν 40ν 24 6 2 ν 2 Для вычисления постоянных интегрирования ci , где i = 1, 6, определяющих структуру решений (7), сформулируем граничные условия. На нижней границе, описываемой уравнением плоскости z = 0, выпол- няются условия прилипания: Vx = Vy = 0, которые из-за структуры класса (3) записываются следующим образом: U = V = 0, u = 0. (8) На верхней границе z = h заданы скорости Vx = W cos ϕ + yΩ, Vy = W sin ϕ, 166 что эквивалентно условиям U = W cos ϕ, V = W sin ϕ, u = Ω. (9) В (9) W — значение скорости (модуль вектора) жидкости на поверхности z = h слоя жидкости, а угол ϕ — это направление этой скорости, относительно двумерной прямоугольной декартовой системы координат Oxy; Ω — величина градиента скорости, совпадающего с точностью до знака со значением вертикальной ∂Vy ∂Vx компоненты завихренности на верхней границе слоя жидкости: Ωz = − = −Ω. ∂x ∂y Воспользовавшись значениями скоростей на границах (8) и (9), получим частное решение, описывающее течение вязкой несжимаемой жидкости, индуцируемое неоднородным распределением поля скоростей: µ ¶ z P2 h2 ³ z ´2 P 2 h2 z u=Ω , V = + W sin ϕ − , h 2ν h 2ν h µ ¶ P2 Ωh4 ³ z ´5 ΩW sin ϕh2 P2 Ωh4 ³ z ´4 P1 h2 ³ z ´2 U= + − + + 40ν 2 h 12ν 24ν 2 h 2ν h µ ¶ ΩW sin ϕh2 P2 Ωh4 P1 h2 z + W cos ϕ − + − . (10) 12ν 60ν 2 2ν h 4 Исследование и обсуждение свойств гидродинамических полей При анализе движения жидкости из соображения подобия удобно записывать уравнения Навье-Стокса в безразмерном виде. При точном интегрировании уравнений движения можно приводить к безразмерному виду не только уравнения, но и решения. Далее будет реализован второй путь при анализе поля скоростей. Выражения (10) в безразмерном виде записываются следующим образом: µ ¶ δ 2 Ga 2 δ 2 Ga VY = Z + sin ϕ − Z, 2Re 2Re µ ¶ δ 4 GaRo 5 Taδ 2 sin ϕ δ 4 GaRo δ 2 Ga 2 UX = Z + − Z4 + a Z + 80 24 48 2Re µ ¶ Taδ 2 sin ϕ δ 4 GaRo δ 2 Ga 1 + cos ϕ − + −a Z+ ZY. (11) 24 120 2Re 2Ro Wl P2 l 3 W Re 2Ωl2 Здесь Re = , Ga = , Ro = = , Ta = — число Рейнольдса, Галилея, Россби, моди- ν ν2 2Ωl Ta ν h фицированное число Тейлора второго рода соответственно, δ = — отношение характерных размеров; l P1 z x y a= — отношение градиентов давления; Z = , X = , Y = — безразмерные координаты. Отметим, P2 h l l что приведение к безразмерной форме (11) решений (10) осуществлялось посредством нормировки скоро- стей на число W. Отметим, что при приведении решений (11) к безразмерному виду для горизонтальных координат используется конечный характерный масштаб. Это позволяет использовать формулы (11) для описания течения в замкнутых объемах, задавая на боковых границах y = −l (Y = −1) и y = l (Y = 1) условия проницаемости. Рассмотрим для начала случай, когда W = 0, Ω = 0 на плоскости z = h. В этом случае жидкость покоится на обеих границах и формулы (11) определяют комбинацию однонаправленных классических течений Пуазёйля: δ 2 Ga 2 δ 2 Ga δ 2 Ga aδ 2 Ga 2 aδ 2 Ga aδ 2 Ga Vx = Z − Z= Z (Z − 1) , Vy = Z − Z= Z (Z − 1) . 2 2 2 2 2 2 Отметим, что для представления скоростей Vx и Vy в безразмерном виде, нужно в размерных выражениях ν (10) провести нормировку скоростей на число . Очевидно, что при выполнении условий прилипания на l обеих границах слоя, поле скоростей в жидкости не расслаивается, но достигает экстремального значения в середине слоя. Характер экстремума определяется знаком коэффициента, стоящего перед слагаемым 167 Z 2 . Иными словами, поскольку решения являются непрерывно дифференцируемыми функциями, в жид- кости не существуют застойные (критические) точки, которые характеризуют наличие противотечений, а усиление скоростей может иметь место. Аналогично можно получить комбинацию классических течений Куэтта, полагая P1 = P2 = 0, Ω = 0 в формулах (10). Если в (11) потребовать равенства нулю Ga, a, или в решениях (10) положить P1 = P2 = 0, то имеем неоднородное течение Куэтта, изученное подробно в статьях [11, 12, 13]. В [11, 12, 13] было пока- зано, что при P1 = P2 = 0 и отсутствии вращения в жидкости может существовать одна застойная точка. Данные решения могут описывать течения жидкости в экваториальной зоне океана. Точные решения (11) при Ga = a позволяют уточнить механизм образования противотечений в жидкости. Они образуются из-за нелинейных эффектов в жидкости, которые проявляются при наличии неоднородного распределения ско- ростей (ветра) на одной из границ даже при постоянном давлении. Таким образом, градиентное течение не является основополагающим при объяснении противотечений в экваториальной зоне океана, но оно может влиять на количество застойных точек в потоке вязкой несжимаемой жидкости. Исследуем далее количество нулей скоростей VX и VY . Скорость VY может быть записана в виде: · 2 µ ¶¸ δ Ga δ 2 Ga VY = Z Z + sin ϕ − . 2Re 2Re Очевидно, что скорость VY может принимать только одно нулевое значение на области определения Z ∈ 1 2Re sin ϕ [0; 1] в точке Z = − . 2 δ 2 Ga Перейдем к исследованию скорости VY . Однородное слагаемое U скорости VY допускает мультиплика- тивное представление: · 4 µ ¶ δ GaRo 4 Taδ 2 sin ϕ δ 4 GaRo δ 2 Ga U =Z Z + − Z3 + a Z+ 80 24 48 2Re µ ¶¸ Taδ 2 sin ϕ δ 4 GaRo δ 2 Ga ¡ ¢ + cos ϕ − + −a = Zf (Z) = Z AZ 4 − BZ 3 − CZ − D . 24 120 2Re Знак минус перед коэффициентами A, B, C и D поставлен из технических соображений. Исследование количества нулей функции f на отрезке [0; 1] сведем к исследованию функционального уравнения f1 = f2 : AZ 4 − BZ 3 = CZ + D. Здесь f1 = AZ 4 − BZ 3 (кривая 4 на рис. 1), f2 = CZ + D (прямые 1, 2, 3, 5, 6, 7 на рис. 1). Очевидно, что при A = B = 0 справедливо уравнение f = f2 = −CZ − D = 0, которое имеет не более чем одно решение (прямые 5, 6, на рис. 1). Только один корень может существовать у уравнения f = 0, если A = 0 или B = 0, когда соответственно имеем одно из уравнений −BZ 3 = CZ + D, или AZ 4 = CZ + D. В данных ситуациях уравнение f = 0 имеет только единственное решение, поскольку f1 и f2 являются монотонными функциями, которые могут пересекаться только в одной точке или не пересекаться вовсе (прямые 1, 2, 5, 6, 7 на рис. 1). Отметим, что при C = 0, уравнение f = 0 может иметь решения внутри отрезка [0; 1] (прямая 3 на рис. 1). Прямые 3 и 6 (рис. 1) при фиксированном значении углового коэффици- ента C являются предельными прямыми. Если при параллельном переносе прямая принадлежит полосе, ограниченной прямыми 3 и 6 на рис. 1, то многочлен f имеет ровно один корень на отрезке [0; 1] . Отметим, что прямая 6 (рис. 1) характеризует существование кратных корней Z = 0 и Z = 1. Установим далее ограничения на коэффициенты A, B, C и D, при которых возможно два нуля у функции B B f (рис 2). Функция f1 имеет единственный корень Z = при ∈ (0; 1) , что эквивалентно разбиению A A функции на два участка монотонности (кривая 4 на рис. 2). На одном участке функция f1 монотонно 168 Рис. 1: Графики функций f1 и f2 , иллюстрирующие существование одного нуля у многочлена f B убывает, а на другом — возрастает. Это означает, что для любой комбинации ∈ (0; 1) существуют такие A числа C и D, что уравнение f = 0 имеет 2 корня (прямые 3 и 5 на рис. 2). При A > 0 параметр C должен быть положительным, а при A < 0 — отрицательным. Следовательно, если числа A и C имеют разные знаки, то уравнение f = 0 имеет не более одного корня (прямая 1 на рис. 2). При C = 0 в этом случае у уравнения f = 0 всегда будет не больше одного корня (кривая 4 на рис. 2). Единственный корень будет существовать при касании прямой 6 кривой 4 (рис. 2). Таким образом, многочлен f не может иметь на отрезке [0; 1] три корня. Иными словами, задание неоднородного распределения скоростей и постоянных градиентов давления по толщине слоя жидкости, образует трехслойную структуру течения (рис. 3). В средней области II на рис. 3 скорость будет направлена в противоположную сторону относительно скоростей в областях I и III, примыкающих к границам слоя (рис. 3). Следовательно, учет градиентного течения позволяет использовать трехслойные модели океана (кусочно-непрерывно стратифицированы) при проведении теоретических исследований. Исследование числа нулей у полиномиального поля скоростей можно проводить другими методами. Можно использовать теоремы о локализации нулей многочленов. Например, теорему Рауса-Гурвица и их аналоги, в том числе и частотные критерии. Для многочленов невысокой степени исследование пробле- мы Рауса-Гурвица можно заменить изучением спектральных свойств. Для существования одного корня функция должна принимать экстремальное значение хотя бы в одной точке, но это значение должно быть противоположного знака с граничной скоростью. Аналогичные рассуждения справедливы при существо- вании двух застойных точек. Для этого функция должна иметь хотя бы два экстремума со значениями противоположного знака. Решения (11) определяют касательные напряжения: µ ¶ δ 2 Ga δ 2 Ga τY Z = τZY = Z + sin ϕ − , Re 2Re µ ¶ δ 4 GaRo 4 Taδ 2 sin ϕ δ 4 GaRo 2δ 2 Ga τXZ = τZX = Z + − Z3 + a Z+ 16 6 12 2Re Taδ 2 sin ϕ δ 4 GaRo δ 2 Ga 1 + cos ϕ − + −a + Y. 24 120 2Re 2Ro Напряжение τY Z , очевидно, обращается в нуль не более, чем в одной точке на области определения. Од- нородное слагаемое τXZ может принимать нулевые значения на [0; 1] не более чем в трех точках. Это можно 169 Рис. 2: Графики функций f1 и f2 , иллюстрирующие существование двух нулей у многочлена f показать, используя алгоритм приведенный выше. Отсюда вытекает очень важное свойство обобщенного неоднородного течения Пуазёйля, заключающееся в том, что существует такое значение Z ∈ [0; 1] , при котором обе компоненты тензора касательных напряжений равны нулю: τY Z = τXZ = 0. Иными словами, в жидкости существует такой слой, в котором сила трения тождественно равна нулю. В частности, трение может отсутствовать при Z = 0. При исследовании структуры течения вязкой несжимаемой жидкости (11), получим, что движение во всем пространстве является вихревых. Вычисляя компоненты завихренности, получим ∂Vy ∂V Ωx = − =− , ∂z ∂z ∂Vx ∂U ∂u Ωy = = +y , ∂z ∂z ∂z ∂Vy ∂Vx Ωz = − = −u. ∂x ∂y Z Ранее было показано, что Ωz = − . Следовательно, завихренность не может принять нулевое значение. 2Ro Таким образом, рассматриваемое движение жидкости не может локально и глобально вырождаться в по- тенциальное. Кроме того, решения (11) формируют спиральную структуру течения, для количественного описания которой используется спиральность поля скорости. Вклад в генерацию спиральности осуществ- ляется вертикальной и горизонтальной компонентами движения. Ввиду этого, полученные решения будут полезны в качестве тестовых при исследовании свойств спиральной турбулентности и проведении числен- ного моделирования течений диссипативной жидкости. Метки ссылок должны начинаться с фамилии первого автора [?]. Список литературы [1] S. N. Aristov , A. D. Polyanin. New classes of exact solutions of Euler equations. Doklady Physics, 53(3): 166-171, 2008. [2] S. N. Aristov , E. Y. Prosviryakov. A New Class of Exact Solutions for Three-Dimensional Thermal Diffusion Equations. Theoretical Foundations of Chemical Engineering, 50(3): 286-293, 2016. 170 Рис. 3: Профиль скорости Vx при наличии двух застойных точек при фиксированном значении y [3] P. G. Drazin and N. Riley. The Navier–Stokes Equations: A Classification of Flows and Exact Solutions. Cambridge: Cambridge Univ. Press, 2006. [4] V. V. Kuznetsov, V. V. Pukhnachev A New Family of Exact Solutions of Navier-Stokes Equations. Doklady Physics, 54(3): 126-130, 2009. [5] C. C. Lin. Note on a class of exact solutions in magneto-hydrodynamics. Arch. Rational Mech. Anal., 1: 391-395, 1958. [6] S. V. Meleshko. A particular class of partially invariant solutions of the Navier–Stokes equations. Nonlinear Dynamic, 36(1): 47-68, 2004. [7] M. Couette. Études sur le frottement des liquides. Ann. Chim. Phys., 21: 433-510, 1890. [8] V. W. Ekman. On the influence of the earth’s rotation on ocean currents. Ark. Mat. Astron. Fys., 2: 1-5, 1905. [9] H. Stommel. Wind-drift near the equator. Deep-Sea Res., 6(4): 298-302, 1960. [10] T. G. Charney. Nonlinear theory of a wind-driven homogenous layer near the equator. Deep Sea Res., 6(4): 303-310, 1960. [11] S. N. Aristov , E. Y. Prosviryakov. Inhomogeneous Couette flow. Rus. J. Nonlin. Dyn., 10(2): 177-182, 2014. (In Russian) [12] S. N. Aristov , E. Y. Prosviryakov. Large-Scale Flows of Viscous Incompressible Vortical Fluid. Russian Aeronautics, 58(4): 413-418, 2015. [13] S. N. Aristov , E. Y. Prosviryakov. Unsteady Layered Vortical Fluid Flows. Fluid Dynamics, 51(2): 148-154, 2016. 171 Layered gradient stationary flow vertically swirling viscous incompressible fluid Eugenii Yu. Prosviryakov Institute of Engineering Science, The Ural Branch of Russian Academy of Sciences (Yekaterinburg, Russia) Ural State University of Railway Transport (Yekaterinburg, Russia) Abstract. An exact solution describing the flow of a vortical viscous incompressible fluid is derived. The solution of the Navier-Stokes equation possesses a characteristic feature in describing a fluid in motion, namely, it holds true when not only viscous but also inertia forces are taken into account. Taking the inertia forces into account leads to the appearance of stagnation points in a fluid layer and counterflows, as well as the existence of layer thicknesses at which the tangent stresses vanish on the lower boundary. Pressure role in the formation of stagnant fluid layer at the points is discussed. Keywords: exact solution, countercurrent, Hagen-Poiseuille flow, laminar flow, stagnation point. 172