УДК 532.529.2
Структура индуцированных диффузией течений на клине с искривленными гранями
© 2016 Н.Ф. Димитриева*, Ю.Д. Чашечкин**
*Институт гидромеханики НАН Украины, Киев, Украина E-mail: [email protected]
**Институт проблем механики им. А.Ю. Ишлинского РАН, Москва, Россия E-mail: [email protected]
Поступила в редакцию 02.02.2016 г.
Разработаны математическая модель и метод численной реализации, позволяющие одновременно изучать все элементы внутренне многомасштабных стратифицированных течений. С использованием суперкомпьютерных систем впервые рассчитаны и визуализированы индуцированные диффузией течения на симметричном неподвижном клине с искривленными гранями. Показано влияние краевых эффектов и кривизны боковой поверхности клина на структуру течения. Дефицит давления в примыкающих струйных течениях формирует интегральную силу, вызывающую самодвижение свободного клина вдоль горизонта нейтральной плавучести в устойчиво стратифицированной среде в направлении вершины.
Ключевые слова: стратифицированная жидкость, индуцированные диффузией течения, самодвижение.
Введение. Растворенные в жидкости вещества и взвешенные частицы под действием массовых сил (гравитационных или центробежных) распределяются неравномерно и формируют устойчивую стратификацию, поле плотности в которой определяется профилями температуры, концентрации примеси и давления в соответствии с видом уравнения состояния. Неравновесная среда с молекулярными потоками стратифицирующих компонент находится в состоянии покоя, только если градиенты плотности параллельны линии действия силы тяжести. Прерывание молекулярного потока на непроницаемых границах произвольной формы формирует течения, индуцированные диффузией, которые начали теоретически изучать в начале 40-х годов прошлого века [1] и несколько позднее стали моделировать в лаборатории [2].
Интерес к изучению градиентных течений сохраняется в силу их распространенности в окружающей среде. В обычных условиях возмущения сосредоточены в тонких слоях у непроницаемых поверхностей и достигают штормовых значений при формировании больших градиентов температуры в атмосфере на крутых склонах горных систем [3] и вблизи ледников [4]. В последние годы особое внимание уделяется изучению влияния течений, индуцированных диффузией, на динамику и структуру процессов в водохранилищах, озерах, морях [5] и океанах с учетом влияния глобального вращения [6]. В ряде работ учитывались сложность уравнения состояния морской среды и наличие дополнительной стратифицирующей компоненты с собственным коэффициентом диффузии («differential diffusion») [7].
В конце прошлого века наряду с аналитическими исследованиями стационарных течений [1, 2, 5] стали разрабатываться модели процесса установления течений, индуцированных диффузией на различных препятствиях -МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 3 2016 77
наклонной и горизонтальной пластине [8, 9], цилиндре, сфере, а также стали проводиться лабораторные исследования с помощью высокочувствительных теневых приборов. В экспериментах у полюсов препятствий кроме больших вихрей наблюдались «розетки» нестационарных диссипативно-гравитацион-ных волн, которые с удалением от препятствия трансформировались в протяженные прослойки [10].
Интегральное силовое действие конвергентных течений, которое отсутствует на симметричных препятствиях (сфере, цилиндре или горизонтальной пластине), принимает конечное значение на наклонной пластине и других несимметричных относительно линии действия силы тяжести телах. Возникающие градиенты давления достаточно велики и способны вызывать самодвижение свободных тел нейтральной плавучести («diffusion fish» по терминологии из статьи [11]), играющее важную роль в динамике гидросферы. Расчет установившихся течений на плоском неподвижном клине выполнен в работе [12]. Практический интерес представляет изучение механизма формирования пропульсивных сил на телах различной геометрии.
В данной работе впервые выполнен расчет динамики установления и структуры двумерного течения, индуцированного диффузией на выпуклом или вогнутом симметричном клине.
Постановка задачи и метод решения. Поскольку скорости исследуемых течений малы по сравнению со скоростью звука, математическая модель изучаемых физических процессов базируется на системе дифференциальных балансных уравнений механики неоднородных многокомпонентных жидкостей в приближении Буссинеска и пренебрежении эффектами сжимаемости [13], которая включает уравнения неразрывности, баланса импульса (Навье - Сто-кса), вещества и замыкающее алгебраическое уравнение состояния p(S (z)):
Р = Ро (exp(- z/ Л)+s ), (1)
divv = 0, (2)
rs 1
— + (vV)v =--VP + vAv - sg, (3)
dt ро
ds v
— + vVs = ks As , (4)
dt s Л
где S = S0 (z) + s - полная соленость, включающая коэффициент солевого сжатия, в которой выделены исходная стационарная стратификация и возмущение S(x, z, t), р0 - плотность на нулевом уровне (горизонте нейтральной плавучести); v = (vx,vz) - скорость жидкости; Р - давление за вычетом гидростатического; v и ks - постоянные коэффициенты кинематической вязкости и диффузии соли; t - время; g - ускорение свободного падения; V и A -операторы Гамильтона и Лапласа; Л = р / | Vр |, N = -yfg/Л , Tb = 2п/N -масштаб, частота и период плавучести; ось z направлена вертикально вверх. Член с вертикальной компонентой скорости течения vz в формуле (4) харак-
78
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 3 2016
теризует конвективным перенос солености, нарушающим исходную стратификацию.
В начальный момент времени t = 0 в покоящуюся непрерывно стратифицированную жидкость горизонтально помещается непроницаемый клин выбранной формы. На поверхности тела Е задается условие прилипания для скорости и непротекания для вещества, на бесконечном удалении - затухание всех возмущений:
V, " = 0, у| = 0, V, 4 = 0, (5)
' и<0 ' II ' ' I х,' 4 '
д8 1 дг д"
= ---+
дп I Л дп дп
= 0, (6)
I
где п - внешняя нормаль к поверхности препятствия.
В задаче представлен набор размерных параметров (у = 10-6м2/с, К = 1,41-10-9 м2/с, gz = 9,8 м/с2, N = 1 с-1), которые формируют характерные масштабы времени t = Ть, а также линейные макромасштабы (плавучести Л = 9,8 м, размера препятствия Ь = 0,1 м) и микромасштабы (5УЫ =-,]у/N = 10-3 м и ЗК" =^¡кs|N = 4 -10-5 м). Характерный масштаб скорости задается отношением базовых масштабов иуы = л!ш = 10-3 м/с, иК" = 4 -10-5 м/с.
Макромасштабы Л и Ь характеризуют исходную стратификацию и геометрию течения. Микромасштабы диссипативной природы З^ и определяют поперечные размеры тонкоструктурных компонентов. Широкий диапазон значений масштабов длины указывает на сложность внутренней структуры стратифицированного течения.
Исследования решений линеаризованных фундаментальных уравнений и результаты лабораторного моделирования показали, что крупномасштабные элементы течений (волны и вихри) характеризуются регулярно возмущенными компонентами полного решения. Тонкоструктурные элементы течений, которые визуализируются во всем диапазоне параметров изучаемых процессов, описываются членами обширного семейства сингулярно возмущенных компонент. Все компоненты сосуществуют и активно взаимодействуют между собой, формируя эволюционирующую тонкую структуру [14].
Отношения масштабов задают традиционные безразмерные комплексы -шкалу плавучести С = Л/ Ь , числа Рейнольдса Яе = |у| Ь /V , Пекле
Ре = |у|Ь/К", Шмидта 8с = у/К" и Фруда Бг = |у|/Ж , где |у| - величина характерной скорости течения, которая при рассмотрении задачи течения, индуцированного диффузией на неподвижном препятствии, принимается равной иК".
В нелинейной постановке система (1) - (6) позволяет одновременно изучать все элементы течений в рамках единого описания в естественных физических переменных без привлечения дополнительных констант и связей. В
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 3 2016 79
силу многомасштабности явлений и нелинейности уравнений одним из основных инструментов анализа таких эволюционирующих процессов служит численное моделирование.
Дискретизация расчетной области осуществлялась в открытой интегрируемой платформе SALOME (www.salome-platform.org), а также с использованием встроенных утилит blockMesh, topoSet и refineMesh открытого пакета OpenFOAM (www.openfoam.com) с учетом собственных масштабов. Расчетная область представляет собой прямоугольник, в котором горизонтально расположен симметричный клин длиной L = 10 см и высотой основания h = 2 см. Боковые стороны тела - дуги окружностей, радиусы и положения центров которых в декартовой системе координат с центром на вершине клина приведены в табл. 1.
Т а б л и ц а 1
Геометрия клина
Координаты центра, см
Отклонение от грани хорды AR, см
(-7,6; 126,7) (17,6; 125,7)
126,9 125,7
-0,1 0,1
(9,2; 41,4) (5,8; 126,7)
43,2 10,0
0,3 1,4
С учетом формы тела построены блочно-структурированные гексаэд-ральные расчетные сетки с совмещением линий на границах блоков. Для разрешения тонкоструктурных компонентов течения на минимальном микромасштабе SKs размещались нескольких ячеек расчетной сетки в областях с большими значениями градиентов. Проведенные тестовые расчеты на более грубой сетке выявили вблизи препятствия локальные нефизические осцилляции решения, которые при продолжительных вычислениях приводят к накоплению погрешности и остановке счета.
Численное решение строилось методом конечных объемов в открытом пакете OpenFOAM. Для учета эффектов стратификации и диффузии стандартная численная модель icoFoam, реализующая нестационарные течения однородной жидкости, была дополнена новыми переменными ( р и s) и соответствующими уравнениями (1), (4), а также новыми вспомогательными параметрами (N, Л , Ks, g и др. [15]).
Граничное условие для возмущения солености (6) реализовано с помощью расширенной утилиты funkySetBoundaryField, которая позволяет задавать аналитические выражения для физических переменных. Для проведения последующего анализа с помощью утилит vorticity, funkySetFields, probes, forces и forceCoeffs вычислялись отсутствующие в стандартных программах
80
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 3 2016
ОрвпЕОАЫ информативные физические величины - завихренность и темп ее бароклинной генерации, скорость диссипации механической энергии.
Результаты и обсуждение. Сложную многомасштабную структуру индуцированных диффузией течений иллюстрирует поле горизонтальной компоненты градиента возмущения солености (рис. 1). Положительные значения визуализируемой величины обозначены зеленым цветом, отрицательные -синим. Между изолиниями одинаковый диапазон изменения величины. Структура изображения согласуется с данными теневой визуализации течений в лабораторном бассейне [11].
Рис. 1. Поле горизонтальной компоненты градиента возмущения солености дя/дх на клине
с искривленными гранями (Ь = 10 см, к = 2 см, Ть = 6,28 с, т = Г/Т4 = 16): а - Ая = -0,1 см; б - Ая = 0,3 см; в - Ая = 1,4 см
Около угловых точек клина формируются дополнительные тонкоструктурные элементы - диссипативно-гравитационные волны, в которых максимумы продольной компоненты градиента возмущения солености достигали
|дя/дх| = 4 -10-2 (возмущения солености я имели порядок 10-5). С удалением от препятствия градиент возмущения солености резко падает и составляет около 10-6 на расстоянии 5 см по горизонтали и 0,5 см по вертикали от угловой точки.
МОРЖОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 3 2016
81
Структура поля горизонтальной компоненты градиента возмущения солености существенно зависит от знака кривизны боковой поверхности: чем острее вершины у основания клина, тем ярче выражены визуализируемые пучки знакопеременных полос (рис. 1, а). В то же время максимальные значения дх слабо зависят от кривизны граней клина.
I ' I тах А А
Для выпуклого клина (рис. 1, б), у которого угол между основанием и боковой гранью приближается к 90°, пучок тонкоструктурных элементов расплывается.
Для клина, максимальная высота которого смещена по отношению к основанию, картина течений перестраивается (рис. 1, в) и становится подобной наблюдаемой на цилиндре [8].
г, см 2-
0
-2
а
г, см 2
0
-2
6
2, СМ 2
0
-2
в
г, см 20-
-2-10 -5 0 5 10 15 дт, см
г
Рис. 2. Поля горизонтальной Ух (а, б) и вертикальной Уг (в, г) компонент скорости течений на клине (Ь = 10 см, к = 2 см, Ть = 6,28 с, т = Г/Г4 = 16): а, в - АК = -0,1 см; б, г - АК = 0,3 см
Непроницаемое препятствие, погруженное в неподвижную устойчиво стратифицированную жидкость, формирует сложную систему течений,
82 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 3 2016
включающую тонкие главные струи вдоль наклонных сторон с примыкающими компенсационными противотечениями (рис. 2). С удалением от препятствия ширина слоев увеличивается, а максимальная скорость падает.
В распределении горизонтальной компоненты скорости ух толщина главной струи примыкающего течения растет с уменьшением локального угла наклона касательной к горизонту нейтральной плавучести. Соответственно для вогнутого клина толщина уменьшается от 0,8 см у острой вершины до 0,2 см вблизи основания (рис. 2, а), а для выпуклого клина, наоборот, увеличивается (рис. 2, б).
В поле вертикальной компоненты скорости у2 (рис. 2, в, г) главные струи (слои с максимальным значением скорости) отрываются от поверхности в угловых точках у основания клина и формируют систему диссипативно-гравитационных волн. В середине основания клина располагается протяженная невозмущенная область, так как вертикальная поверхность не нарушает диффузионный перенос примеси.
V , мкм/с
-1
-2
vz, мкм/с 1,00,5-
0,0
-0,5
-1,0
-1
-------2
----3
.....4
\ I ; I
-1
1
■ I
I
о 6
1 2, СМ
Рис. 3. Профили горизонтальной Ух (а) и вертикальной (б) компонент скорости у вершины (кривые 1, 2) и основания (кривые 3, 4) клина (Ь = 10 см, к = 2 см, Ть = 6,28 с, т = 1/Ть = 16): кривые 1, 3 - Ая = -0,1 см; кривые 2, 4 - Ая = 0,3 см
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 3 2016
83
Распределения компонент скорости течения вблизи вершин клина представлены на рис. 3. Около острой вершины, точки бифуркации течения, наблюдается скачкообразное изменение вертикальной компоненты скорости (кривые 1, 2 на рис. 3, б). Значения компонент скорости зависят от угла наклона грани клина к оси х.
В окрестности угловых точек основания клина формируется сложная знакопеременная структура, в которой вертикальная компонента скорости Уг становится сравнимой по абсолютному значению с горизонтальной ух (кривые 3, 4 на рис. 3, б).
У острой вершины клина, а также в тонком течении вдоль его боковых сторон возмущение давления убывает (рис. 4), его значение зависит от кривизны боковых сторон. При этом толщина протяженной области дефицита давления на вогнутом клине при Ак = -0,1 см составляет около 3 см (рис. 4, а), что значительно больше, чем на выпуклом клине: при Ак = 0,3 см толщина области 2,0 см (рис. 4, б), при Ак = 1,4 см толщина области 1,5 см (рис. 4, в). Для выпуклого клина на рис. 4, в дефицит давления наблюдается у всех вершин, что обусловлено изменением угла наклона граней.
Рис. 4. Поля возмущений давления течений на клине (Ь =10 см, к = 2 см, ТЪ = 6,28 с, т = ИТЪ = 16): а - АК = -0,1 см; б - АК = 0,3 см; в - АК = 1,4 см
84
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 3 2016
Разность давлений создает горизонтальную компоненту интегральной силы Рх, толкающую горизонтальный клин в направлении вершины, что согласуется с наблюдениями [11, 12]. Рассчитанные значения Рх приведены в табл. 2 при т = t|Tъ = 16 .
Т а б л и ц а 2
Зависимость горизонтальной компоненты результирующей силы, действующей на клин, от искривления его поверхностей (т = t|Tъ = 16)
Ак, см -0,1 0,0 0,1 0,3 1,4
Р , Н 8,4-Ю-6 8,3-Ю-6 8,0-Ю-6 -6,1-10-6 -2,2-10-6
Вертикальные компоненты результирующей силы Рг, вычисленные отдельно для верхней и нижней граней, превосходят ¥х на два порядка. Вследствие антиколлинеарности они компенсируют действие друг друга.
Заключение. На основе фундаментальной системы уравнений разработаны математическая модель и программа построения численного решения, позволяющие одновременно изучать все элементы внутренне многомасштабных стратифицированных течений без привлечения дополнительных гипотез и связей. Критерии построения и оптимизации сеток сформулированы с учетом масштабного анализа системы уравнений. Модель впервые использована для расчета и визуализации течений, индуцированных диффузией на симметричном неподвижном клине с искривленными гранями.
Показано влияние краевых эффектов и кривизны поверхностей на структуру течения. Дефицит давления на клине, возникающий вследствие перестройки исходной стратификации, формирует интегральную силу, способную вызывать самодвижение свободного клина, как и других несимметричных тел, вдоль горизонта нейтральной плавучести. Тело в устойчиво стратифицированной среде перемещается в направлении своей вершины, что подтверждается результатами выполненных экспериментов. Величина возмущения давления зависит от кривизны боковых сторон клина.
Тонкие струйные течения жидкости, формирующиеся вдоль каждой из сторон клина, порождают в областях схождения с угловых точек внутренние волны, образующие тонкие прослойки.
Результаты расчетов течений, индуцированных диффузией, могут использоваться в качестве начальных условий при решении задач обтекания препятствий стратифицированной морской средой.
Вычисления проводились в ^еЪ-лаборатории итНиВ (www.unihub.ru) и на суперкомпьютере ВЦ МГУ «Ломоносов» (www.parallel.ru). Работа выполнена в ИПМех РАН при финансовой поддержке РФФИ (проект 15-01-09235) и ОЭММПУ РАН (Программа ГУ-4-12 «Динамика формирования и взаимодействия волн и вихрей в сплошных средах»).
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 3 2016
85
СПИСОК ЛИТЕРАТУРЫ
1. Прандтль Л. Гидроаэромеханика. - М.: Изд-во иностр. лит., 1949. - 488 с.
2. Phillips O.M. On flows induced by diffusion in a stably stratified fluid // Deep-Sea Res. -1970. - 17. - P. 435 - 443.
3. Shapiro A., Fedorovich E. A boundary-layer scaling for turbulent katabatic flow // Bound.-Lay. Meteorol. - 2014. - 153, Issue 1. - P. 1 - 17.
4. Oerlemans J., van Pelt W.J.J. A model study of Abrahamsenbreen, a surging glacier in northern Spitsbergen // The Cryosph. - 2015. - 9, № 2. - P. 767 - 779.
5. Зырянов В.Н., Лапина Л.Э. Склоновые течения в морях, озерах и водохранилищах, обусловленные диффузионными эффектами // Водные ресурсы. - 2012. - 39, № 3. -С. 292 - 303.
6. Garrett C., MacCready P., Rhines P.B. Boundary mixing and arrested Ekman layers: rotating, stratified flow near a sloping boundary // Ann. Rev. Fluid Mech. - 1993. - 25. - P. 291 - 323.
7. GargettA.E. Differential diffusion: an oceanographic primer // Progr. Oceanogr. - 2003. - 56. - P. 559 - 570.
8. Linden P.F., Weber J.E. The formation of layers in a double-diffusive system with a sloping boundary // J. Fluid Mech. - 1977. - 81. - P. 757 - 773.
9. Чашечкин Ю.Д., Загуменный Я.В. Течения непрерывно стратифицированной жидкости, индуцированные прерыванием диффузионного переноса неподвижной пластиной // Морской гидрофизический журнал. - 2012. - № 5. - C. 3 - 22.
10. Zagumennyi Ia.V., Chashechkin Yu.D. The structure of convective flows driven by density variations in a continuously stratified fluid // Physica Scripta. - 2013. - 155. - article id. 014034.
11. Mercier M.J., Ardekani F.M., Allshouse M.R. et al. Self-propulsion of immersed object via natural convection // Phys. Rev. Lett. - 2014. - 112, № 20. - P. 204501(5).
12. PageM.A. Propelled by diffusion // Nature Physics. - 2010. - 6. - P. 486 - 487.
13. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. VI. Гидродинамика. - М.: Наука, 1986. - 736 с.
14. Чашечкин Ю.Д. Дифференциальная механика жидкостей: согласованные аналитические, численные и лабораторные модели стратифицированных течений // Вестник МГТУ им. Н.Э. Баумана. Сер. «Естественные науки». - 2014. - № 6. - C. 67 - 95.
15. Димитриева Н.Ф., Загуменный Я.В. Численное моделирование стратифицированных течений с использованием OpenFOAM // Тр. ИСП РАН. - 2014. - 26, № 5. - С. 187 -200.
Structure of the diffusion-induced flows on the wedge with the curved surfaces
N.F. Dimitrieva, Yu.D. Chashechkin
Institute of Hydromechanics, National Academy of Sciences of Ukraine, Kiev, Ukraine
e-mail: [email protected] A.Yu. Ishlinsky Institute for Problems in Mechanics, Russian Academy of Sciences, Moscow, Russia
e-mail: [email protected]
The mathematical model and the numerical implementation method permitting to study simultaneously all the elements of the internally multi-scale stratified flows are developed. For the first time, the diffusion-induced flows on a fixed symmetrical wedge with the curved surfaces are calculated and visualized using the supercomputer systems. Influence both of the edge effects and the wedge side surface curvature upon the flow structure is shown. In a stably stratified environment, the integral force inducing self-motion of a free wedge along the neutral buoyancy horizon towards the top is formed by pressure deficit in the adjacent jet streams.
Keywords: stratified fluid, diffusion-induced flows, self-motion.
86 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 3 2016