CC BY
31
ISSN: 2310-7081 (online), 1991-8615 (print) d http://doi.org/10.14498/vsgtu1527
Математическое моделирование, численные методы и комплексы программ
УДК 517.958:532.51
Крупномасштабная слоистая стационарная конвекция вязкой несжимаемой жидкости под действием касательных напряжений на верхней границе. Исследование поля скоростей
Н. В. Бурмашева1'2, Е. Ю. Просвиряков2
1 Уральский федеральный университет им. первого Президента России Б. Н. Ельцина, Россия, 620002, Екатеринбург, ул. Мира, 19.
2 Институт машиноведения УрО РАН, Россия, 620049, Екатеринбург, ул. Комсомольская, 34.
Аннотация
Рассматривается точное решение задачи об определении конвективных движений в слоистых крупномасштабных течениях вязкой несжимаемой жидкости в стационарном случае. Показано, что получаемая задача является, во-первых, переопределенной и, во-вторых, нелинейной (за счет наличия членов конвективной производной в уравнении теплопроводности). Также показано, что выбором класса решений систему можно избавить от переопределенности, а уточнением ряда краевых условий свести задачу к задаче исследования термокапиллярной конвекции (конвекции Бенара-Марангони). Далее определяются условия появления противотечений. Исследуется их возможное количество. Кроме того, проводится анализ исследуемого течения на предмет наличия безвихревых областей. Показано, что при некоторых сочетаниях параметров системы возможны смены направления вихря.
Ключевые слова: слоистое течение, противотечения, застойные точки, точное решение.
Получение: 20 января 2017 г. / Исправление: 6 марта 2017 г. / Принятие: 13 марта 2017 г. / Публикация онлайн: 19 мая 2017 г.
Научная статья
3 ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Бурмашева Н. В., Просвиряков Е. Ю. Крупномасштабная слоистая стационарная конвекция вязкой несжимаемой жидкости под действием касательных напряжений на верхней границе. Исследование поля скоростей // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2017. Т. 21, № 1. С. 180-196. doi: 10.14498/vsgtu1527. Сведения об авторах
Наталья Владимировна Бурмашева А http://orcid.org/0000-0003-4711-1894 кандидат технических наук; доцент; институт математики и компьютерных наук, каф. механики и математического моделирования1; научный сотрудник; лаб. механики деформаций2; e-mail:[email protected]
Введение. Известно, что в потоке жидкости могут существовать точки, где нулевое значение принимает одна или несколько компонент вектора скорости. Такие точки называются застойными, или критическими [1]. Рядом с такими точками, как правило, образуются застойные зоны, то есть области с обратным течением. Данное утверждение справедливо для решений уравнений Навье—Стокса и их модификаций, являющихся непрерывными функциями координат. Другие точные решения, обладающие застойными точками, приведены в [2-4].
Противотечения возникают, например, в следе за цилиндром при небольших числах Рейнольдса и Фруда в жидкости с небольшим периодом плавучести [5]. Экспериментальные исследования показали, что противотечения наблюдаются в поверхностных, подповерхностных и глубинных слоях практически во всех районах Мирового океана [6-9].
Проблема учета и исследования застойных зон оказывается актуальной не только в задачах, связанных с обтеканием движущихся в жидкой или газообразной среде тел применительно, в частности, к судам и летательным аппаратам. Она возникает, например, при использовании скважин в нефтедобыче [10], когда необходимо определить характер течения в окрестности точки изменения направления потока. Данная проблема также имеет место при обтекании ветром зданий и является определяющей для расчета выбросов вредных веществ техногенных систем.
Теоретическое изучение этих и многих других процессов, связанных с исследованием противотечений, в настоящее время далеко от завершения. В связи с этим отыскание новых точных решений уравнений гидродинамики, описывающих появление застойных зон в жидкости, является актуальной зада-
Первым изотермическим точным решением, посвященным исследованию застойных точек, является течение Хименца [11]. Исследование противотечений в океане было начато Экманом [12]. В [12] было представлено точное решение, описывающее течение вязкой несжимаемой жидкости при балансе вязких и кориолисовых сил. Давление в этом случае полагалось однородным по горизонтальным координатам.
На основе течений Экмана [12], Чарни [13] и Стоммела [14] были предложены различные математические модели описания экваториальных противотечений. Дальнейшие исследования в этой области опирались на вычислительный эксперимент. Отметим, что в [15-20] построены широкие классы точных решений для описания жидкостей при вращении.
Несмотря на большое число работ, в которых учитывается планетарная завихренность, особый интерес представляют течения жидкости в отсутствие поля сил Кориолиса. В этом случае для описания течений можно использовать «бесконечно замкнутый» горизонтальный слой. В [21-23] приведены точные решения крупномасштабных изотермических течений при постоянном давлении и конвективных движений различной природы. В статьях [24-35] приведены точные решения для конвективных термокапиллярных течений,
Евгений Юрьевич Просвиряков Э http://orcid.org/0000-0002-2349-7801
кандидат физико-математических наук; ведущий научный сотрудник; лаб. прикладной
механики2; e-mail: [email protected]
которые являются обобщениями течения Остроумова—Бириха [36,37].
В данной работе мы рассматриваем движение вязкой несжимаемой жидкости под действием заданных на верхней границе тангенциальных сил, не являющихся результатом капиллярного эффекта.
1. Постановка задачи. Рассмотрим задачу описания конвективных движений в слоистых крупномасштабных течениях вязкой несжимаемой жидкости в стационарном случае. Наиболее распространенным подходом в таких задачах является использование уравнений тепловой конвекции в приближении Буссинеска, предполагающего аппроксимацию плотности как линейную функцию температуры. С математической точки зрения данная модель [38] состоит из
- уравнения Навье—Стокса
(V • Ч)У = -ЧР + VДУ + д/ЗТк, (1)
- уравнения теплопроводности
V • ЧТ = хДТ, (2)
- уравнения несжимаемости
V ■ V = 0. (3)
Здесь V(x, y, z) = {Vx, Vy, Vz} — вектор скорости, причем Vz = 0, так как рассматривается слоистое течение; P — отклонение давления от гидростатического, деленное на постоянную среднюю плотность р жидкости; T — отклонение от средней температуры; v, % — коэффициенты кинематической вязкости и температуропроводности жидкости соответственно; k — орт оси z, направленной вертикально вверх; V — оператор Гамильтона; Д — двумерный оператор Лапласа.
Заметим, что данная система суть система из пяти скалярных уравнений относительно четырех неизвестных (компонент скорости VX, Vy и физических полей T и P). То есть система переопределена. Следовательно, чтобы решать такую систему, надо либо убедиться в непротиворечивости системы условий, либо искать такое решение, на котором одно из имеющихся условий удовлетворится тождественно.
Кроме того, отметим, что члены конвективной производной из уравнения Навье-Стокса и уравнения несжимаемости тождественно равны нулю, но они сохраняются в уравнении теплопроводности. Поэтому, строго говоря, решается нелинейная задача.
2. Условие разрешимости переопределенной системы. Для того чтобы разрешить переопределенную исходную систему (1)-(3) предлагается пойти вторым из предложенных выше путей и искать решение в классе, предложенном в [20,24]:
Vx = u(z), Vy = v(z). (4)
В этом случае уравнение неразрывности (3) обращается в тождество:
dVx + OVy + dVz = du(z) + dvM =0 dx dy dz dx dy
Далее будем искать давление и температуру в виде специальных линейных форм по горизонтальным координатам:
Т = То(г) + ВДж + :вду, Р — Ро (г) + Р^х + Р^у. (5)
Подставим выбранный класс (4), (5) в уравнение Навье—Стокса (1):
т ВУХ т , дУу т ЯУг дР
Ух + Уу ^ + Уг = - — + V ДУх, дх ду дг дх
«(г) ^ «(г) ^+0 ^ — - дХ + V ДУх, дх ду дг дх
д / д2 д2 д2 \
- (ро(г)+р1(г)х+р2ш = + ^ + Ы «(г),
откуда согласно методу неопределенных коэффициентов получаем
Р1 = V«".
Здесь и всюду далее штрихом обозначено дифференцирование по переменной г. Аналогично из второго уравнения Навье—Стокса (относительно компоненты Уу) имеем
Р2 = V«".
Из третьего уравнения сохранения импульсов в жидкости имеем
дУХ дУу дУг дР Ух ^ + Уу ^ + Уг = - т- + V ДУг + двТ дх ду дг дг
«(г) ^ + «(г) ^^+0 £ = - £ + двТ,
£=дет
или
д
— (Ро(г) + Р1(г)х + Р2(г)^ = дв(То(г) + ВДх + Т2(г)у),
откуда согласно методу неопределенных коэффициентов получаем три скалярных соотношения:
дР-
—^ = двТ-, г = 0,1,2. дг
И, наконец, подставим выбранный класс решений (4), (5) в уравнение теплопроводности (2):
дТ удТ у
'дж у ду г дг \ дх2 ду2 дг2 /'
дТо(г) + Т1(г)х + Т2(г)у + дТо(г) + ВДх + Т2(г)у =
« м +« м —
дх ду
= х( дх ++) сад+ад*+ад^
■Т + ^Т2 = хТО'(^) + хТ'(г)ж + хТВДу
Отсюда согласно методу неопределенных коэффициентов имеем
■Т + ^Т2 = хТО'(^), ВД = ТВД = 0.
Таким образом, окончательно получаем нелинейную систему уравнений:
Т1' = о, Т2' = о.
дР1
= двТ1,
дР2
vu// = Р1, = Р2
хТо' = иТ1 + ^Т2,
= двТ2
дРо
= двТо,
(6)
(7)
(8) (9)
(10)
3. Краевая задача. Система обыкновенных дифференциальных уравнений (6)-(10) является системой тринадцатого порядка. Следовательно, чтобы корректно определить краевую задачу, необходимо иметь тринадцать условий, заданных на границе слоя, которые будут использованы для определения констант, возникающих в ходе интегрирования рассматриваемой системы.
Будем далее считать, что нижняя граница слоя жидкости является абсолютно твердой и неподвижной. Верхнюю границу будем полагать свободной. В качестве краевых условий рассмотрим следующие:
- пусть на нижней границе (г = 0) слоя жидкости выполняется условие прилипания, а температура задается функцией
Т = А* + Ву;
- на верхней границе (г = Л) действует постоянное давление 5, а температура задается функцией
Т = 0 + С* + Бу;
- кроме того, на свободной поверхности г = Л заданы напряжения:
= ^ = ^
Таким образом, приходим к следующей системе граничных условий:
То (0) = 0, То (Л) = 0, Ро (Л) = 5,
и(0) = -(0) = 0, Т1(0) = А, Т1 (Л) = С, Р1(Л) = 0,
(Л) =
Т2(0) = В, Т2(Л) = Б, Р2(Л) = 0,
(П) — «2.
Без ограничения общности можно считать 5 — 0, тем самым ведя отсчет приведенного давления от уровня, задаваемого на верхней границе.
4. Решение системы уравнений. Рассмотрим частный случай поставленной выше краевой задачи — задание температурного возмущения только на верхней границе (А — В — 0). Начнем последовательно интегрировать основную систему уравнений с уравнения (6). С учетом соответствующих граничных условий получаем
г г
Т1 — С-, Т2 — Я-. п п
Далее полученные выражения подставляем в уравнения (7) и находим градиенты Р1, Р2 приведенного давления:
г
2 п2 г2 п2
Р1 — две——, Р2 — дв^
2п 2 2п
Затем из уравнений (8) находим компоненты вектора скорости:
«М — 8(^(2П2г - ^). «(г) — 8(^(2П2г - ^)■
Далее найденные выражения компонент скорости и градиентов температуры подставляем в уравнение (9), откуда
т (г) г(-П3(С«1 + Р6) + 12ХП^ + № + Р6)г3) + То(-) — ^ +
СРдвг(-16П6 + 21П2г4 - 5г6) + 840хП2^ ,
и, наконец, из уравнения (10) определим составляющую Ро давления:
Ро (г) — шХПП +
2вд(С«1 + Р«2)г5 в2СРд2(41П8 - 64П6г2 + 28П2г6 - 5г8) + 120хпП + 6720хЛ2^ .
5. Переход к термокапиллярной конвекции. Рассматриваемая краевая задача может быть сведена к задаче исследования термокапиллярной конвекции (конвекции Бенара—Марангони), описанной в [24]. Для этого достаточно уточнить условия на нижней (твердой) и верхней (свободной) поверхностях, положив
То (П)—0, ^ (П)— «1 — -аЗД, п^ (П)— «2 — -аТ2(П).
6. Исследование компонент скорости. Функции компонент скорости u, v имеют ноль в начале координат (множитель z выделяется явно). Изучим вопрос о том, сколько еще нулей может быть у этих функций. Этот вопрос является довольно принципиальным, поскольку наличие нулей у этих функций говорит о существовании противотечений (расслоении) в конвективном потоке вязкой несжимаемой жидкости.
Пронормируем найденные в п. 4 выражения компонент u(z), v(z) вектора скорости V на величину C, приведя их к безразмерному виду. Для этого введем новые параметры системы:
А* = D/C, S = h/l,
где h — характерный вертикальный размер слоя, а l — характерный горизонтальный размер слоя. Кроме этого, перейдем к безразмерной координате z ^ z/h. Для этого разделим u(z), v(z) на gß13/v. В результате получим
zS3 zS3
u(z) = ^ (A*(2z - z3) + 8W1), v(z) = ^ (2z - z3 + 8W2) . 8 8
Здесь
W = ^
двСпЛ2
— числа Вебера, посчитанные для значений 6, г = 1, 2. Очевидно, что
^1/^2 = 6/6.
Отметим также, что после обезразмеривания координата г теперь меняется в диапазоне [0,1].
Рассмотрим сначала функцию и(г). Случай Д* = 0 опустим из-за тривиальности. Поэтому далее представим функцию и(г) в мультипликативном виде:
"М = ^ (2. - + ^)
и введем в рассмотрение функцию
■1(г) = 2г - г3 + .
Эта функция является полиномом третьей степени относительно г.
Чтобы определить, сколько нулей у функции "1(г), рассмотрим вспомогательную функцию /(г) = 2г — г3, график которой представлен на рис. 1.
Отметим, у функции /(г) на [0,1] есть всего один ноль — точка .о = 0. Кроме этого, глобальный максимум функции / на интересующем нас отрезке [0,1] есть точка г* = л/2/3, значение функции в которой равно
^ - ^ - 1 /тах = /(г*) = з^з > 1.
А на правой границе г = 1 она принимает значение /1 = /(1) = 1.
О 0.2 0.4 0.6 0.8 1.0
z
Рис. 1. График функции f (z) [Figure 1. The graph of f (z)]
Производя параллельный перенос графика этой функции вдоль вертикальной оси, то есть рассматривая семейство функций вида fc(z) = f (г) + с, можно добиться того, что функция ^(г) пересечет ось г на отрезке [0,1] один или два раза, или же вообще не будет пересекать ось г (не имея при этом общих точек или только единственную общую точку — точку касания). Различные положения функции ^ представлены на рис. 2.
Таким образом, возвращаясь к функции щ, получаем следующие оценки:
1) -
8Wi
Д*
< 1 —единственный ноль, одна смена знака; 8Wi
2) 1 = ——--два различных нуля, одна смена знака;
3) 1 < -
8Wi _
3^3
*
<
два различных нуля, две смены знака;
4^/2 8Wi 3^3
*
один ноль, функция ui(z) не меняет знак;
Рис. 2. Семейство функций fc(z) = f (z) + c для c = —f (1) + 0.05 (линия 1), для c = —(f (z*) + f (1))/2 (линия 2), для c = —f (z*) (линия 3), для c = —f (z*) — 0.05 (линия 4)
[Figure 2. The family of functions fc(z) = f (z) + c for c = —f (1) + 0.05 (line 1), for c = —(f (z* ) + f (1))/2 (line 2), for c = —f (z*) (line 3), for c = —f (z*) — 0.05 (line 4)]
4^2 8^! „ 5) < ——--нулей нет, функция иця) не меняет знак.
В соответствии с каждым из перечисленных случаев для функции ч(г) можем получить следующие картины (без ограничения общности взято 5 = 1, Д* = 1) (см. рис. 3).
s -0.02
Рис. 3. Семейство функций u(z) = ^Г" (2z - z3 + iW*) Лля Wi = -(f С1) - a05) (линия 1)
для Wi = -(f (z*) + f (1))/2 (линия 2), для Wi = -f (z*) (линия 3), для Wi = -(f (z*) + 0.05)
(линия 4), когда S = 1, A* = 1
(2z - z3 + 8W1) for Wi = -(f (1) - 0.05)
[Figure 3. The family of functions u(z)
(line 1), for Wi = -(f (z*) + f (1))/2 (line 2), for Wi = -f (z*) (line 3), for Wi = -(f (z*) + 0.05) (line 4), when S = 1, A* = 1]
Проводя аналогичный анализ для функции vi (z), можно получить следующую систему условий, определяющую расположение этой кривой относительно оси z:
1) — 8W2 < 1 —единственный ноль, одна смена знака;
2) 1 = —8W2 — два различных нуля, одна смена знака; 4—2
1 < —8W2 < -—^ —два различных нуля, две смены знака;
4V2
3—3
4—2
3—3
3^3
= —8W2 — один ноль, функция v1(z) не меняет знак;
< —8W2 — нулей нет, функция vi(z) не меняет знак.
Отметим, что если вернуться к изначальным параметрам задачи, то получим
Wi
Д* = gfiDnh2'
W2 =
g^Cr/h2'
Таким образом, наличие противотечений в слое жидкости, их число и положение застойных точек зависят от значения комбинаций параметров д^щь?
а дРСф2 '
7. Годограф скорости. Особый интерес в этой ситуации представляет изучение траектории, описываемой концом вектора-годографа скорости. Чисто теоретически может быть не более 25 различных форм годографа (5 вариантов для кривой п(х) и 5 вариантов для кривой у(г)), но есть ситуации, которые порождают один и тот же годограф (например, при Д* = 1
в случаях Wi = W2 = —1 и Wi = W2 = — годографом будет биссектриса первой четверти). Кроме того, не следует забывать про связь между параметрами задачи (Wi^2 = W2^i), которая приводит к дополнительным ограничениям и, следовательно, к снижению числа принципиально разных траекторий, описываемых концом годографа.
Итак, если функция ui(z) такова, что имеет единственный ноль и меняет знак один раз, то годограф имеет вид, представленный на рис. 4, a.
Если функция ui(z) имеет два нуля и меняет знак один раз, то годограф имеет вид, представленный на рис. 4, b.
Если функция ui(z) имеет два нуля, но меняет знак два раза, то годограф имеет вид, представленный на рис. 4, с.
Если функция ui(z) имеет один ноль, но не меняет знак, то годограф имеет вид, представленный на рис. 4, d.
Если у функции ui(z) нулей нет (функция не меняет знак), то годограф имеет вид, представленный на рис. 4, е.
8. Исследование завихренности и спиральности. Выпишем основную количественную меру завихренности — псевдовектор вихря или завихренность:
О = rot V =
i j k
д d d
dx dy dz
u v 0
или покоординатно:
дг А3 ди Л3
а* = - ^ = - + г - г3), Оу = ди = у (2т + Д*(г - г3)), О = 0.
Особый интерес здесь вызывают комбинации параметров, при которых завихренность обращается в ноль хотя бы в нескольких точках. Принципиально таких случая два.
Случай первый. При Д* = ^ = = 0 завихренность на интересующем нас интервале обращается в ноль в точках г = 0 и г = 1. В этом случае
¿3
и(г) = 0, г(г) = — г2(2 - г2).
8
То есть противотечений в слое не возникает (поскольку нули функции г(г) лежат вне интервала (0,1)). Жидкость в слое движется в одном строго зафиксированном направлении и имеет место «волновой» эффект изменения модуля скорости движения жидкости по высоте г рассматриваемого слоя. Случай второй. При Д* = 0, Д*^2 = Так как = последнее
условие равносильно условию Д^2 = В этом случае псевдовектор вихря становится нулевым в точках, являющихся корнями уравнения
+ г - г3 = 0.
Число действительных корней этого уравнения зависит от ^2. Экстремум функции 2^2 + г - г3 равен и достигается он в точке г** = .
СО
о
Рис. 4. Вид годографа
Ъ
в плоскости и(г) —
Таким образом, если W2 > 0 или W2 < , то у приведенного выше уравнения корней на отрезке [0,1] нет. Если W2 = — , то существует единственный действительный корень на отрезке [0,1] —точка z**. Если — 373 ^ W2 ^ 0, то таких корня два —один левее z**, другой правее. Кроме того,
u(z) = (2z — z3 + 8W) = Д* z83 (2z — z3 + 8W2) = ДЧ^,
то есть имеет место линейная зависимость между компонентами u и v вектора скорости V. В этом случае годограф скорости описывает своим концом прямую, проходящую через начало координат, тангенс угла наклона которой равен 1/Д*. А значит, если координата z такова, что удовлетворяет уравнению 2W2 + z — z3 = 0, то направление движения жидкости на этой высоте слоя оказывается заранее известным.
Таким образом, в зависимости от сочетаний конкретных значений параметров задачи можно попасть в ситуации, когда вихрь меняет свое направление.
Исследуем теперь спиральность [39]. По определению, спиральность есть число
H = V ■ О = V ■ rot V = — u^ + vdu = — (Wi — Д*W2)z2(—2 + 3z2).
dz dz 8
Очевидно, что спиральность оказывается равной нулю при Wi — Д*W2 = 0, а это есть одно из условий обращения завихренности в ноль. Это выглядит логичным, учитывая, что завихренность Q входит множителем в определение спиральности H.
Однако есть два случая, когда спиральность обращается в ноль, а завихренность может и не быть нулевой — это случай z = 0 и z = л^/2/3 (нас интересуют только точки слоя 0 ^ z ^ 1). В точке z = 0 вектор скорости V оказывается равен нулю (условие прилипания), а компоненты завихренности определяются выражениями
= —W2^3, Пу = Wi£3.
Если Wi = W2 = 0, то завихренность наряду со спиральностью становится нулевой при z = 0. Аналогичной подстановкой точки z = л^/2/3 легко можно показать, что если Wi = — 376, то обращаются в ноль u(^/2/3) и Qy, а если
W2 = — 3"76, то также в ноль обращаются v(^/2/3) и Qx.
Таким образом, даже если спиральность стала нулевой, то в зависимости от значений констант в граничных условиях завихренность может в ноль и не обращаться.
Конкурирующие интересы. Мы заявляем, что у нас нет конфликта интересов в отношении авторства и публикации этой статьи.
Авторская ответственность. Мы несем полную ответственность за предоставление окончательной рукописи в печать. Каждый из нас одобрил окончательную версию рукописи.
Финансирование. Работа выполнена при поддержке фонда содействия развитию
малых форм предприятий в научно-технической сфере (программа УМНИК).
Библиографический список
1. Аристов С. Н., Князев Д. В., Полянин А. Д. Точные решения уравнений Навье—Сток-са с линейной зависимостью компонент скорости от двух пространственный переменных// Теоретические основы химической технологии, 2009. Т. 43, №5. С. 547-566.
2. Dorrepaal J. M. An exact solution of the Navier-Stokes equation which describes nonorthogonal stagnation-point flow in two dimensions // Journal of Fluid Mechanics, 1986. vol.163, no. 1. pp. 141-147. doi: 10.1017/s0022112086002240.
3. Stuart J. T. The viscous flow near a stagnation point when the external flow has uniform vorticity// Journal of the Aerospace Sciences, 1959. vol.26, no. 2. pp. 124-125. doi: 10. 2514/8.7963.
4. Riesco-Chueca P., de la Mora J. F. Brownian motion far from equilibrium: a hypersonic approach// Journal of Fluid Mechanics, 1990. vol.214. pp. 639-663. doi: 10.1017/ S0022112090000301.
5. Гущин В. А., Рождественская Т. И. Численное исследование явлений, возникающих вблизи кругового цилиндра в течениях стратифицированных жидкостей с небольшими периодами плавучести // Прикладная механика и техническая физика, 2011. Т. 52, №6. С. 69-76.
6. Штокман В. Б. Экваториальные противотечения в океанах. Ленинград: Гидрометео-издат, 1948. 156 с.
7. Коротаев Г. К., Михайлова Э. Н., Шапиро Н. Б. Теория экваториальных противотечений в Мировом океане. Киев: Наук. думка, 1986. 208 с.
8. Коротаев Г. К. Теоретическое моделирование синоптической изменчивости океана. Киев: Наук. думка, 1988. 160 с.
9. Бондаренко А. Л. Крупномасштабные течения и долгопериодные волны Мирового океана. М.: Институт водных проблем РАН, 2011. 163 с.
10. Яруллин А. Р. Результаты экспериментальных исследований двухфазного расслоенного потока в горизонтальном стволе скважины со знакопеременной траекторией // Каротажник, 2014. №9(243). С. 63-71.
11. Hiemenz K. Die Grenzschicht an einem in den gleichförmigen Flüssigkeitsstrom eingetauchten geraden Kreiszylinder // Dingler's Politech. J, 1911. vol. 326. pp. 321-324, http: //dingler.culture.hu-berlin.de/article/pj326/ar326114.
12. Ekman V. W. On the Influence of the Earth's Rotation on Ocean-Currents // Ark. Mat. Astron. Fys., 1905. vol.2, no. 11. pp. 1-52, http://jhir.library.jhu.edu/handle/1774. 2/33989.
13. Charney J. G. Non-linear theory of a wind-driven homogeneous layer near the equator // Deep Sea Research, 1960. vol.6. pp. 303-310. doi: 10.1016/0146-6313(59)90089-9.
14. Stommel H. Wind-drift near the equator // Deep Sea Research, 1960. vol. 6. pp. 298-302. doi:10.1016/0146-6313(59)90088-7.
15. Аристов С. Н., Фрик П. Г. Нелинейные эффекты влияния экмановского слоя на динамику крупномасштабных вихрей в «мелкой воде» // Журнал прикладной механики и технической физики, 1991. Т. 32, №2. С. 49-54.
16. Ингель Л. Х., Аристов С. Н. Класс точных решений нелинейных задач о термических циркуляциях, связанных с объемным тепловыделением в атмосфере // Тр. Ин-та эксперим. метеорол., 1996. №27(162). С. 142-157.
17. Аристов С. Н., Шварц К. Г. Вихревые течения адвективной природы во вращающемся слое жидкости. Пермь: ПГУ, 2006. 155 с.
18. Аристов С. Н., Шварц К. Г. Вихревые течения в тонких слоях жидкости. Киров: ВятГУ, 2011. 207 с.
19. Аристов С. Н., Шварц К. Г. Адвективное течение во вращающейся жидкой пленке // Прикладная механика и техническая физика, 2016. Т. 57, №1. С. 216-223. ((од.: 10. 15372/РМТР20160121.
20. Аристов С. Н., Просвиряков Е. Ю. Новый класс точных решений уравнений термодиффузии// Теоретические основы химической технологии, 2016. Т. 50, №3. С. 294-301. ао1:10.7868/80040357116030027.
21. Аристов С. Н., Просвиряков Е. Ю. Неоднородные течения Куэтта// Нелинейная динам., 2014. Т. 10, №2. С. 177-182. аод.: 10.20537/па1402004.
22. Аристов С. Н., Просвиряков Е. Ю. Крупномасштабные течения завихренной вязкой несжимаемой жидкости // Известия высших учебных заведений. Авиационная техника, 2015. №4. С. 50-54.
23. Аристов С. Н., Просвиряков Е. Ю. Нестационарные слоистые течения завихренной жидкости// Изв. РАН. Механика жидкости и газа, 2016. №2. С. 25-31. ((од.: 10.7868/ 80568528116020055.
24. Аристов С. Н., Просвиряков Е. Ю. О слоистых течениях плоской свободной конвекции// Нелинейная динам., 2013. Т. 9, №4. С. 651-657. аод.: 10.20537/па1304004.
25. Аристов С. Н., Просвиряков Е. Ю., Спевак Л. Ф. Нестационарная конвекция Бе-нара-Марангони слоистых течений вязкой несжимаемой жидкости // Теоретические основы химической технологии, 2016. Т. 50, №2. С. 137-146. ((од.: 10.7868/ 80040357116020019.
26. Аристов С. Н., Просвиряков Е. Ю., Спевак Л. Ф. Нестационарная слоистая тепловая и концентрационная конвекция Марангони вязкой несжимаемой жидкости // Вычислительная механика сплошных сред, 2015. Т. 8, №4. С. 445-456. аоД: 10.7242/1999-6691/ 2015.8.4.38.
27. Андреев В. К. Решения Бириха уравнений конвекции и некоторые его обобщения: Препринт ИВМ СО РАН № 1-10. Красноярск, 2010. 68 с.
28. Андреев В. К., Бекежанова В. Б. Устойчивость неизотермических жидкостей (обзор) // Прикладная механика и техническая физика, 2013. Т. 54, №2. С. 3-20.
29. Андреев В. К., Степанова И. В. Однонаправленные течения бинарных смесей в модели Обербека-Буссинеска// Изв. РАН. Механика жидкости и газа, 2016. №2. С. 13-24. аоД: 10.7868/80568528116020043.
30. Гончарова О. Н., Кабов О. А. Гравитационно-термокапиллярная конвекция в горизонтальном слое при спутном потоке газа// Докл. РАН, 2011. Т. 426, №2. С. 183-188.
31. Гончарова О. Н., Резанова Е. В. Пример точного решения стационарной задачи о двухслойных течениях с испарением на границе раздела // Прикладная механика и техническая физика, 2014. Т. 55, №2. С. 68-79.
32. Бирих Р. В., Пухначев В. В. Осевое конвективное течение во вращающейся трубе с продольным градиентом температуры// Докл. РАН, 2011. Т. 436, №3. С. 323-327.
33. Бирих Р. В., Пухначев В. В., Фроловская О. А. Конвективное течение в горизонтальном канале с неньютоновской реологией поверхности при нестационарном продольном градиенте температуры // Изв. РАН. Механика жидкости и газа, 2015. №1. С. 192-198.
34. Пухначев В. В. Нестационарные аналоги решения Бириха // Известия АлтГУ, 2011. № 1-2. С. 62-69.
35. Рыжков И. И. Термодиффузия в смесях: уравнения, симметрии, решения и их устойчивость. Новосибирск: СО РАН, 2013. 200 с.
36. Остроумов Г. А. Свободная конвекция в условиях внутренней задачи. М.: Гостехиздат, 1952. 256 с.
37. Бирих Р. В. О термокапиллярной конвекции в горизонтальном слое жидкости // Прикладная механика и техническая физика, 1966. Т. 7, №3. С. 69-72.
38. Гершуни Г. З., Жуховицкий Е. М. Конвективная устойчивость несжимаемой жидкости. М.: Наука, 1972. 392 с.
39. Арнольд В. И. Математические методы классической механики. М.: Едиториал УРСС, 2003. 416 с.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2017, vol. 21, no. 1, pp. 180-196
d http://doi.org/10.14498/vsgtu1527
ISSN: 2310-7081 (online), 1991-8615 (print)
MSC: 76F02, 76M45; 76F45, 76R05, 76U05
A large-scale layered stationary convection
of an incompressible viscous fluid under the action of shear
stresses at the upper boundary. Velocity field investigation
N. V. Burmasheva1'2, E. Yu. Prosviryakov2
1 Ural Federal University named after the First President of Russia B. N. Yeltsin,
19, Mira st., Ekaterinburg, 620002, Russian Federation.
2 Institute of Engineering Science, Urals Branch, Russian Academy of Sciences,
34, Komsomolskaya st., Ekaterinburg, 620049, Russian Federation.
Abstract
The exact solution of the definition of convective motions in a layered large-scale flows of a viscous incompressible fluid in a steady case is considered. It was shown that the received problem is, firstly, overdetermined and, secondly, a nonlinear (due to the presence of members of a convective derivative in a heat conduction equation). Also it was shown that the solution class choice can eliminate the override, and the specification of a boundary conditions can reduce the problem to the study of a thermal capillary convection (convection Benard-Marangoni). Then conditions of the counterflow appearance are defined, and their possible amount is investigated. In addition, the analysis of the nonvortex region in the test flow is made. And it was shown that under certain combinations of system parameters the vortex can change the direction.
Keywords: layered flow, counterflow, stagnant point, exact solution.
Received: 20th January, 2017 / Revised: 6th March, 2017 / Accepted: 13th March, 2017 / First online: 19th May, 2017
Competing interests. We declare that we have no conflicts of interest in the authorship and publication of this article.
Authors' responsibilities. We are fully responsible for submitting the final manuscript in print. Each of us has approved the final version of manuscript.
Funding. This work was supported by the Foundation for Assistance to Small Innovative Enterprises in Science and Technology (the UMNIK program).
Research Article
3 ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Burmasheva N. V., Prosviryakov E. Yu. A large-scale layered stationary convection of an incompressible viscous fluid under the action of shear stresses at the upper boundary. Velocity field investigation, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2017, vol. 21, no. 1, pp. 180-196. doi: 10.14498/vsgtu1527 (In Russian). Authors' Details:
Natalya V. Burmasheva http://orcid.org/0000-0003-4711-1894
Cand. Techn. Sci.; Associate Professor; Institute of Mathematics and Computer Science, Dept. of Mechanics and Mathematical Modeling1; Researcher; Laboratory of Mechanics of Deformations 2; e-mail: [email protected]
Evgeny Yu. Prosviryakov http://orcid.org/0000-0002-2349-7801
Cand. Phys. & Math. Sci.; Leading Researcher; Laboratory of Applied Mechanics2;
e-mail: [email protected]
References
1. Aristov S. N., Knyazev D. V., Polyanin A. D. Exact solutions of the Navier-Stokes equations with the linear dependence of velocity components on two space variables, Theoretical Foundations of Chemical Engineering, 2009, vol.43, no. 5, pp. 642-662. doi: 10.1134/ S0040579509050066.
2. Dorrepaal J. M. An exact solution of the Navier-Stokes equation which describes nonorthogonal stagnation-point flow in two dimensions, Journal of Fluid Mechanics, 1986, vol.163, no. 1, pp. 141-147. doi: 10.1017/s0022112086002240.
3. Stuart J. T. The viscous flow near a stagnation point when the external flow has uniform vorticity, Journal of the Aerospace Sciences, 1959, vol. 26, no. 2, pp. 124-125. doi: 10.2514/ 8.7963.
4. Riesco-Chueca P., de la Mora J. F. Brownian motion far from equilibrium: a hypersonic approach, Journal of Fluid Mechanics, 1990, vol.214, pp. 639-663. doi: 10.1017/ S0022112090000301.
5. Gushchin V. A., Rozhdestvenskaya T. I. Numerical study of the effects occurring near a circular cylinder in stratified fluid flows with short buoyancy periods, J. Appl. Mech. Tech. Phys., 2011, vol.52, no. 6, pp. 905-911. doi: 10.1134/S0021894411060083.
6. Shtokman V. B. Ekvatorial'nye protivotecheniia v okeanakh [Equatorial Countercurrents in the Oceans]. Leningrad, Gidrometeoizdat, 1948, 156 pp. (In Russian)
7. Korotaev G. K., Mikhailova E. N., Shapiro N. B. Teoriia ekvatorial'nykh protivotechenii v Mirovom okeane [Theory of Equatorial Countercurrents in the World Ocean]. Kiev, Nauk. dumka, 1986, 208 pp. (In Russian)
8. Korotaev G. K. Teoreticheskoe modelirovanie sinopticheskoi izmenchivosti okeana [Theoretical Simulation of the Synoptic Variability of the Ocean]. Kiev, Nauk. dumka, 1988, 160 pp. (In Russian)
9. Bondarenko A. L. Krupnomasshtabnye techeniia i dolgoperiodnye volny Mirovogo okeana [The large-scale flows and long-period waves of the World Ocean]. Moscow, 2011, 163 pp. (In Russian)
10. Yarullin A. R. The results of experimental studies on a two-phase flow bundle in horizontal wells with a variable-sign trajectory, Karotazhnik, 2014, no. 9(243), pp. 63-71 (In Russian).
11. Hiemenz K. Die Grenzschicht an einem in den gleichförmigen Flüssigkeitsstrom eingetauchten geraden Kreiszylinder, Dingler's Politech. J, 1911, vol.326, pp. 321-324, http: //dingler.culture.hu-berlin.de/article/pj326/ar326114.
12. Ekman V. W. On the Influence of the Earth's Rotation on Ocean-Currents, Ark. Mat. Astron. Fys., 1905, vol.2, no. 11, pp. 1-52, http://jhir.library.jhu.edu/handle/1774. 2/33989.
13. Charney J. G. Non-linear theory of a wind-driven homogeneous layer near the equator, Deep Sea Research, 1960, vol.6, pp. 303-310. doi: 10.1016/0146-6313(59)90089-9.
14. Stommel H. Wind-drift near the equator, Deep Sea Research, 1960, vol. 6, pp. 298-302. doi:10.1016/0146-6313(59)90088-7.
15. Aristov S. N., Frik P. G. Nonlinear effects of the Ekman layer on the dynamics of large-scale eddies in shallow water, J. Appl. Mech. Tech. Phys., 1991, vol.32, no. 2, pp. 189-194. doi: 10.1007/BF00858033.
16. Ingel' L. Kh., Aristov S. N. The class of exact solutions of nonlinear problems on thermal circulation associated with volumetric heat release in the atmosphere, Tr. In-ta eksperim. meteorol., 1996, no. 27(162), pp. 142-157 (In Russian).
17. Aristov S. N., Shvarts K. G. Vikhrevye techeniia advektivnoi prirody vo vrashchaiushchemsia sloe zhidkosti [Vortical Flows of Advective Nature in a Rotating Fluid Layer]. Perm, Perm State Univ., 2006, 155 pp. (In Russian)
18. Aristov S. N., Shvarts K. G. Vikhrevye techeniia v tonkikh sloiakh zhidkosti [Vortical Flows in Thin Fluid Layers]. Kirov, Vyatka State Univ., 2011, 207 pp. (In Russian)
19. Aristov S. N., Shvarts K. G. Advective flow in a rotating liquid film, J. Appl. Mech. Tech. Phys., 2016, vol.57, no. 1, pp. 188-194. doi: 10.1134/S0021894416010211.
20. Aristov S. N., Prosviryakov E. Yu. A new class of exact solutions for three-dimensional thermal diffusion equations, Theor. Found. Chem. Eng., 2016, vol.50, no. 3, pp. 286-293. doi:10.1134/S0040579516030027.
21. Aristov S. N., Prosviryakov E. Yu. Inhomogeneous Couette flow, Nelin. Dinam., 2014, vol.10, no. 2, pp. 177-182 (In Russian). doi: 10.20537/nd1402004.
22. Aristov S. N., Prosviryakov E. Yu. Large-scale flows of viscous incompressible vortical fluid, Russ. Aeronaut., 2015, vol.58, no. 4, pp. 413-418. doi: 10.3103/S1068799815040091.
23. Aristov S. N., Prosviryakov E. Yu. Unsteady layered vortical fluid flows, Fluid Dyn., 2016, vol.51, no. 2, pp. 148-154. doi: 10.1134/S0015462816020034.
24. Aristov S. N., Prosviryakov E. Yu. On laminar flows of planar free convection, Nelin. Dinam., 2013, vol.9, no. 4, pp. 651-657 (In Russian). doi: 10.20537/nd1304004.
25. Aristov S. N., Prosviryakov E. Yu., Spevak L. F. Unsteady-state Bénard-Marangoni convection in layered viscous incompressible flows, Theor. Found. Chem. Eng., 2016, vol.50, no. 2, pp. 132-141. doi: 10.1134/S0040579516020019.
26. Aristov S. N., Prosviryakov E. Yu., Spevak L. F. Nonstationary laminar thermal and solutal Marangoni convection of a viscous fluid, Vychislitel'naia mekhanika sploshnykh sred [Computational Continuum Mechanics], 2015, vol.8, no. 4, pp. 445-456 (In Russian). doi: 10.7242/1999-6691/2015.8.4.38.
27. Andreev V. K. Resheniia Birikha uravnenii konvektsii i nekotorye ego obobshcheniia [Birikh Solutions to Convection Equations and Some of its Extensions], The ICM SB RAS Preprint no. 1-10. Krasnoyarsk, 2010, 68 pp. (In Russian)
28. Andreev V. K., Bekezhanova V. B. Stability of non-isothermal fluids (Review), J. Appl. Mech. Tech. Phys., 2013, vol.54, no. 2, pp. 171-184. doi: 10.1134/S0021894413020016.
29. Andreev V. K., Stepanova I. V. Unidirectional flows of binary mixtures within the framework of the Oberbeck-Boussinesq model, Fluid Dyn., 2016, vol.51, no. 2, pp. 136-147. doi: 10. 1134/S0015462816020022.
30. Goncharova O. N., Kabov O. A. Gravitational-thermocapillary convection of fluid in the horizontal layer in co-current gas flow, Dokl. Phys., 2009, vol. 54, no. 5, pp. 242-247. doi: 10. 1134/S1028335809050061.
31. Goncharova O. N., Rezanova E. V. Example of an exact solution of the stationary problem of two-layer flows with evaporation at the interface, J. Appl. Mech. Tech. Phys., 2014, vol. 55, no. 2, pp. 247-257. doi: 10.1134/S0021894414020072.
32. Birikh R. V., Pukhnachev V. V. An axial convective flow in a rotating tube with a longitudinal temperature gradient, Dokl. Phys., 2011, vol.56, no. 1, pp. 47-52. doi: 10.1134/ S1028335811010095.
33. Birikh R. V., Pukhnachev V. V., Frolovskaya O. A. Convective flow in a horizontal channel with non-Newtonian surface rheology under time-dependent longitudinal temperature gradient, Fluid Dyn., 2015, vol.50, no. 1, pp. 173-179. doi: 10.1134/S0015462815010172.
34. Pukhnachev V. V. Non-stationary Analogues of the Birikh Solution, Izvestiia AltGU, 2011, no. 1-2, pp. 62-69 (In Russian).
35. Ryzhkov I. I. Termodiffuziia v smesiakh: uravneniia, simmetrii, resheniia i ikh ustoichivost' [Thermodiffusion in Mixtures: Equations, Symmetries, Solutions and its Stability]. Novosibirsk, 2013, 200 pp. (In Russian)
36. Ostroumov G. A. Free convection under the condition of the internal problem, NACA Technical Memorandum 1407. Washington, National Advisory Committee for Aeronautics, 1958.
37. Birikh R. V. Thermocapillary convection in a horizontal layer of liquid, J. Appl. Mech. Tech. Phys., 1966, vol.7, no. 3, pp. 43-44. doi: 10.1007/BF00914697.
38. Gershuni G. Z., Zhukhovitskii E. M. Convective Stability of Incompressible Fluids, Israel Program for Scientific Translations. Jerusalem, Keter Publishing House, 1976, 330 pp.
39. Arnold V. I. Mathematical methods of classical mechanics, Graduate Texts in Mathematics, vol.60. New York, Springer-Verlag, 1978, x+462 pp. doi: 10.1007/978-1-4757-1693-1 ; 1989, xvi+519 pp. doi: 10.1007/978-1-4757-2063-1.
