Научная статья на тему 'Капиллярно-гравитационная устойчивость в слое жидкости, подогреваемом снизу, при наличии плоского течения Куэтта'

Капиллярно-гравитационная устойчивость в слое жидкости, подогреваемом снизу, при наличии плоского течения Куэтта Текст научной статьи по специальности «Физика»

CC BY
41
11
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Лапшин В. Б., Будников А. А., Сидоренко А. В.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Капиллярно-гравитационная устойчивость в слое жидкости, подогреваемом снизу, при наличии плоского течения Куэтта»

Капиллярно-гравитационная устойчивость в слое жидкости, подогреваемом снизу, при наличии плоского течения Куэтта

Лапшин В.Б. (1), Будников А.А. (2), Сидоренко А.В. (arturik@altavista.com) (3)

(1) Государственный Океанографический Институт, (2) Московский Государственный Университет, (3) Московский Физико-Технический Институт

Процессы, происходящие в приповерхостом слое океана, являются одим из определяющих параметров глобальной климатической системы. Физическо-химические свойства поверхностного микрослоя океана и процессы, протекающие в непосредственной близости от границы раздела оказывают большое влияние на потоки тепла, соли и газа между океаном и атмосферой.

Работа посвящена исследованию капиллярно-гравитационной устойчивости в слое жидкости, подогреваемом снизу при наличии плоского течения Куэтта. Данная задача является попыткой моделирования приповерхностного сантиметрового слоя океана при умеренном ветре.

Плоское течение Куэтта с линейным профилем скорости было одной из первых моделей теоретических исследований гидродинамической устойчивости [1]. Рэлей писал о необычайной сложности проблемы устойчивости плоского течения Куэтта. Ламб [2] отмечал, что подавляющее количество авторов склоняются к мнению об устойчивости плоского течения Куэтта, но этот вывод «настолько же возможен, насколько трудно его продемонстрировать».

К настоящему времени, основные теоретические результаты, связанные с устойчивостью течения Куэтта, изложены в целом ряде публикаций [см. напр. 3].

Дополнительные сложности возникают при исследовании устойчивости горизонтального сдвигового течения в слое жидкости подогреваемого (охлаждаемого) снизу (сверху). В монографии Д. Джозефа рассматривается гравитационная устойчивость плоского течения Куэтта в слое жидкости, подогреваемом снизу [3]. Критическое число Рэлея, полученное автором, составляет Re = 1708. При числах Релея, меньших критического числа, возмущения плоского течения Куэтта, не зависящие в начальный момент времени от горизонтальной координаты х, всегда затухают. Более того, этот критерий является необходимым и достаточным условием

устойчивости.

В работах [4, 5] рассматривается устойчивость слоя жидкости со свободной поверхностью, нагреваемого снизу, относительно периодического возмущения. Вместе с этим нам не известны публикации, в которых был бы изложен энергетический анализ капиллярно-гравитационной устойчивости в слое жидкости, подогреваемом снизу, при наличии плоского течения Куэтта. Анализу этого комбинированного течения и посвящена данная статья.

Постановка задачи

Рассмотрим уравнения движения горизонтального слоя вязкой несжимаемой ньютоновской жидкости (рис.1) толщиной d, с плотностью р и вязкостью д, ускорение свободного падения коэффициент термического расширения у, коэффициент термодиффузии к. Пространство над жидкостью заполнено нейтральным газом с пренебрежимо малыми вязкостью и плотностью. Нижняя граница слоя находится на твердой поверхности, температура которой задана.

Газ

11

2 = 1+ лМ

о

X, и

Рис.1 Разрез возмущенного слоя жидкости.

Используется декартова система координат, начало которой находиться на нижней границе слоя жидкости. Нижняя граница жидкости является твердой. Верхняя граница является свободной и задается уравнением Б^): z =1+п(хД). Начальный

АТ

профиль температуры имеет вид Тн = Т0 - ъ. Начальный профиль скорости ^ъ)^ = 0

ъ

U0 • ^ . Скорость на верхней границе ^ё) =

Безразмерный масштаб для скорости основного течения ^ = 1 м/с. В данной постановке задачи используется приближение Буссинеска.

Масштабы обезразмеривания возмущения скорости V, температуры 9, времени ^

да

давления р и поверхностного натяжения а = а0 - ^Т •Т - соответственно к/ё, АТ, ё2/к,

дк/ё2, а0, где АТ (>0) - разность температур между нижней и верхней границей, а0 -значение поверхностного натяжения на свободной поверхности.

Возмущения основного состояния такой системы в безразмерном виде удовлетворяют следующим уравнениям:

+ VjVi,j I = + Re0ki - Рг -1и - Рг -1w ^ ,

где ki = (0,1), mi = (1,0) (1.1)

д0

+ Vi0,i = V2© + w (1.2)

V* = 0, (1.3)

Ту = - рбу + ц(ui,j +Щ,0.

где И - безразмерная скорость основного течения, V - безразмерное скорость возмущения, 0 - безразмерное возмущение температуры, w - безразмерная вертикальная составляющая возмущения скорости, р - возмущенная составляющая давления, 5у - символ Кронекера.

V уАТдё3

Рг = к - число Прандтля, Re = ^ - число Релея.

Баланс напряжений на поверхности:

ОуП = С [К(п)(а(0) + МСп) - 0(Л + ^ 5п2)]ш - М(0,к - пжИкЪ ,

при ъ = 1 + п (14)

Кинематическое условие на поверхности: П = Nvini , при ъ = 1 + п, (1.5)

Уравнение для потока тепла на верхней границе: 1 - N

0дщ = N , при ъ = 1 + п, (1.6)

На нижней границе:

V = 0 = 0, при z = 0, (1.7)

где К(п) = N , N = л/1 + Пх2 , П = (-Пх , 1) NN, * = (1, Пх) NN , а(0) - безразмерный коэффициент поверхностного натяжения,

- Ш ^ ^

М = - число Марангони, С = ^ d - капиллярное число,

О = - число Бонда, а0 - значение коэффициента поверхностного натяжения на

свободной поверхности.

Мы будем предполагать, что каждый из параметров параметры у^ 0, р и п периодичен по оси х с периодом Х0. Кроме требования периодичности и достаточной гладкости у^ 0, р и п произвольны, в частности, они имеют произвольную амплитуду. Блок-схема решения данной задачи представлена на рис.2. Энергетический функционал типа Ляпунова для задачи (1) имеет вид

1 1 л б2 Г1Г 2 уАТп3^ 1

е = <2 рг -1и№> + я,1<2> + ^ |2л2 + JN ,

где использованы обозначения

Х01 + п(хД) S(t) Х0

<...> = J ... dzdx, J ... = ^ ... dS = J ... •Ndx [7,8].

0 0 S(t) 0 0

Первый член в уравнении для энергетического функционала представляет собой механическую энергию, второй член может быть назван балансом энтропии. Оба этих члена тесно связаны с перемещением поверхности п. Третий член учитывает кинематическое условие на верхней границе и является уравнением эволюции поверхности.

> 0, > 0 - параметры связи. Параметры связи являются свободными и выбираются так, чтобы оптимизировать предел устойчивости [6]. В нашем случае для удобства выбирается = О/С.[4,5] Достаточное условие устойчивости (1) дает вариационная задача

Рис.2 Блок-схема энергетического анализа начально-краевой задачи.

dE

+ <2pvi,i> + Р^ п) = 0, где 2р(x,z,t) и в - множители Лагранжа, которые вводятся х

для учета уравнения непрерывности уц = 0 и условия периодичности возмущения ^ п

х

dx = 0.

Поскольку для природных условий С ~ 10 , можно использовать этот параметр для разложения полученного уравнения. Производя разложение членов при вариациях скорости возмущения дщ, температуры 50 и множителей Лагранжа др, дв, в нулевом приближении имеем:

др д2у д2у 52w 1 Зи^)

- дх + 2дХ2 + д"2 + дХяТ - РГ - 1 w -д^- = 0

дх дх га охга га

др 3^ 3\ ^ „ 1 Зи^)

- + + ХТ + + Иф - РГ - 1 У = 0 га Зz Зх охо^ Зz

У2ф + ^ = 0

ду + 5w = 0

Зх Зz

у = w = ф = 0, при z = 0 w = 0, при z = 1

дф 11 Зу

Зz = 2 дх , при z = 1

ду дw 1 1 дф & + & = - 2 М Ж дх, при z = 1

(2.1)

(2.2) (2.3)

(2.4)

(2.5) (2.6)

(2.7)

(2.8)

И = 1 ^ (7 + 7),

(2.9)

2 7

Где введено обозначение ф = л/яё 70.

Решение уравнений (2) находиться методом разделения переменных в виде fi(z)exp(iax), где индекс i = 1, 2, 3, 4 относиться соответственно к p, у, w, и 0.

Чай + 2^а)^3 +£3" + iaf4' - РГ f4 = 0

-fl' + 2Й4'' + (ia)2f4 + iafз' + И£г - РГ^^ЦЗ^^^ *3 = 0

(ia)2f2 + Г + Й4И = 0

iafз + Й4' = 0

f3 = - icf4' = ¿«-Л' + f2"')

1 2 f4 = - ^((ia)2f2 + f2'')

f1 = ^ ^.....+ 2f2''' - a2f2' + iV1 ^ (f2'' - a2f2) }

Далее, приводя полученные уравнения к одной из функций £(ъ), а именно к ^(ъ), получим:

Г2(6) - 3а^(4) + 3а4^2(2) - а2^2 + - 2iaА(a2f2(1) - Г2(3)) = 0,

л т, 1дИ(ъ) где А = Рг - - сдвиговое число. дъ

Общее решение этого уравнения имеет вид £г(ъ) = Е С exp(qiz), i = 1, 2, ... 6, где qi - корни характеристического уравнения

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

- а2)3 + ^2а2 + 2iaАq(q2 - а2) = 0. Соотношение между критическими числами находиться из граничных условий: ЕС = 0

1 Еф2С = 0

1

ia^

E(qi3 - a2qi )Ci = 0

Ш Z(qi2 - a2)exp(qi)Ci = 0

1 2 2 2 1 M 2

Ш (^(qi - a )qi Ciexp(qi)) - 2 a ^Ciexp(qi) = 0

1 M 1 2 2

2 ^Re ш (zqi(qi - a )exp(qi)Ci) - 2qiCiexp(qD = 0.

Полученная связь имеет вид:

ME (R, A) = min max M(X,a, VRe, A), a X

Где a волновое число в направлении х. На рис.3 представлен пример седловой точки на плоскости М( X,a ).

На рис.4 изображены полученные нейтральные кривые при фиксированных числах А в координатах Re, M и a. При числе А = 0, т.е. при отсутствии плоского течения Куэтта полученная нейтральная кривая в плоскости (Re, M) совпадает с нейтральной кривой, приведенной в работах [4, 5].

Рис.3. Зависимость критического числа Марангони от волнового числа а и параметра связи X, при Яе = 700 и А = 10.

Рис.4. Зависимость критических чисел Марангони от критического числа Релея при различных волновых числах а. В плоскости (Re, M) приведены проекции 1', 2', 3', 4', 5' нейтральных кривых 1, 2, 3, 4, 5 при сдвиговых числах 60, 50, 30, 10 и 0 соответственно. Область устойчивости с ростом А уменьшается.

При отсутствии капиллярной конвекции, т. е. при числе Марангони равном нулю, и линейном вертикальном профиле скорости плоского течения Куэтта, критическое число Рэлея = 1708, приведенное в монографии Джозефа [3], получается при волновом числе a = 0.75.

Основные выводы:

1. Впервые проведен анализ капиллярно-гравитационной устойчивости слоя жидкости, подогреваемой снизу, при наличии плоского течения Куэтта. Рассчитывались зависимости критического числа Марангони от числа Релея при заданных значениях сдвигового числа А. Диапазоны значений безразмерных параметров начально-краевой задачи выбирались характерными для приповерхностного слоя океана.

2. В результате решения минимаксной задачи

проведен детальный анализ зависимости критического чисела Марангони от параметра связи 7 и волнового числа a. Поверхность М^Д) имеет вид седла (рис.3). Однако, при волновом числе a ~ 3 на этой поверхности возникает дополнительный максимум, не исчезающий при изменении масштаба сетки. Причем положение этого максимума

3. Построено семейство нейтральных кривых М = M(Re) при различных числах А

4. Показано, что с ростом сдвигового числа область устойчивости уменьшается (рис.4). Так число Марангони при фиксированном числе Рэлея уменьшается в три раза при изменении числа А от 0 до 60.

Литература:

1. Гольдштик М.А., Штерн В.Н. Гидродинамическая устойчивость и турбулентность. -Новосибирск: Наука. 1977. -366с.

2. Ламб Г. Гидродинамика. -М.-Л., ОГИЗ - Гостехиздат. 1947. -928с.

3. Джозеф Д.Д. Устойчивость движений жидкости: пер. С англ. -М.: Мир, 1981. -638с.

ME (R, A) = min max M(X ,a, VRe, A), a X

E

зависит от числа A=Pr 1

_ 1dU(z)

Pr ' в диапазоне от 0 до 100 (рис.4).

4. Сastillo J.L., Velarde M.G. Buoyancy - thermocapillary instability: the role of interfracial deformation in one- and two-component fluid layers heated from below or above. J. Fluid. Mech., 1982. V.125, P.463 - 474.

5. Davis S.H., Homsy G.M. Energy stability theory for free-surface problems: buoyancy -thermocapillary layers. J. Fluid. Mech., 1980, v.98, P.553 - 572.

6. Joseph D.D. On the stability of the Boussinesq equations, Arch. Rat. Mech. Anal., 1965, V.20, No.59, P.564 - 576.

i Надоели баннеры? Вы всегда можете отключить рекламу.