CC BY
43
2013 Механика № 3
УДК 532.5
С.Н. Аристов, Д.В. Князев
Институт механики сплошных сред УрО РАН, Пермь, Россия
АВТОМОДЕЛЬНЫЕ ЛОКАЛИЗОВАННЫЕ КОНВЕКТИВНЫЕ СТРУКТУРЫ
Рассмотрена задача о конвективном течении в слое вязкой жидкости, вызываемом ее локальным нагревом. Поиск решения задачи осуществлялся в рамках класса точных решений уравнений термогравитационной конвекции, обобщающего известный класс решений уравнений Навье-Стокса, к которому относятся вихри Бюргерса и Салливана. Для единичного числа Пран-дтля найдены два семейства автомодельных решений задачи, позволившие описать эволюцию двух различных типов радиально-локализованных вихрей. В обоих случаях радиальная компонента скорости на большом расстоянии от оси симметрии вихрей убывает обратно пропорционально радиусу, в то время как вертикальная составляющая скорости и температура в первом случае затухают как квадрат расстояния от оси, а во втором - экспоненциально. Для азимутальной скорости получено отдельное линейное уравнение с коэффициентами, зависящими от функции тока меридионального течения. В силу автомодельности это уравнение допускает частные решения с разделяющимися переменными, суперпозиция которых дает возможность описать перенос момента импульса (циркуляции, если она отлична от нуля) от бесконечности к центру вихря, а также проследить эволюцию произвольного локализованного начального возмущения азимутальной скорости. Под действием вихревой и тепловой диффузии рассмотренные вихревые образования затухают со временем. Полученные точные решения являются простыми, обозримыми моделями локализованных конвективных вихрей и ранее известны не были.
Ключевые слова: конвекция, уравнения Обербека-Буссинеска, точные решения, локализованные вихри.
S.N. Aristov, D.V. Knyazev
Institute of Continuous Media Mechanics UB RAS, Perm, Russian Federation
SELF-SIMILAR LOCALIZED CONVECTIVE STRUCTURES
The problem of convective flow in a layer of viscous fluid, caused by local heating, is considered. The problem solution was sought within the class of exact solutions of the thermo-gravitational convection generalizing well-known class of solutions of the Navier-Stokes equations involved the Burger's and Sullivan's vortexes. The two families of self-similar solutions of the problem for a unit Prandtl's number are found, it allow to describe the evolution of two different types of radial - localized vortices. In both cases, the radial component of the velocity at a distance from the axis of symmetry of the vortices is inversely proportional to the radius, while the vertical component of the velocity and the temperature in the first case decayed as the square of the distance from the axis, and the second - exponentially. The separate linear equation for the azimuthal velocity with coefficients depending on the current function of meridian flow is obtained. By virtue of the similarity equation admits particular solu-
tions with separated variables, the composition of which allows us to describe the transfer of angular momentum (circulation, if it is different from zero) from infinity to the center of the vortex, as well as to trace the evolution of an arbitrary localized initial perturbation of azimuthal velocity. These vortex formations are decay in the time due to the vortex and thermal diffusion. The Found exact solutions are the simple, usefully models of localized convective vortices and previously were not known.
Keywords: convection, Oberbeck - Boussinesq equations, exact solutions, localized vortexes.
Введение
Движение тяжелой жидкости, возникающее в результате неоднородности распределения плотности, обусловленного неравномерным нагревом среды, называется термогравитационной конвекцией. Простейшим объектом при изучении конвекции является слой жидкости с твердыми или свободными горизонтальными границами. Необходимым условием равновесия такого слоя является вертикальность однородного градиента давления. Конвекция в слое не может возникнуть до тех пор, пока число Рэлея, характеризующее баланс между силами вязкого трения и плавучести (Архимеда), не достигнет критического значения. Например, для слоя с твердыми границами критическое значение числа Рэлея Ra = 1707,762 [1]. Положительность безразмерного комплекса Рэлея означает, что температура жидкости уменьшается с высотой (неустойчивая стратификация), в то время как при подогреве сверху (устойчивая стратификация) плоского горизонтального слоя с твердыми границами конвекция никогда не возникает.
Возникнув, конвективное течение стремится перемешать жидкость так, чтобы выровнять температуру или перейти в новое равновесное состояние, при этом в качестве источника движения используется в том числе запас энергии, содержащийся в первоначальной равновесной стратификации. В связи с этим представляет интерес проблема использования энергетического потенциала термически стратифицированной среды для генерации и поддержания пространственно-локализованных вихревых структур, наблюдаемых в лабораторных [2] и естественных условиях (смерчи, «пыльные дьяволы» и пр.). Далее рассматриваются точные решения уравнений гидродинамики в приближении Обербека-Буссинеска, моделирующие такого рода структуры.
1. Постановка и общий анализ задачи
Пусть в тяжелой жидкости, находящейся над непроницаемой плоскостью z = 0, имеется равновесный градиент температуры VT = -A,e z. В начальный момент времени t = 0 в слой вносится ради-
ально локализованное осесимметричное возмущение, в центре которого вертикальный градиент температуры V2Т = -А0ег. Требуется описать возникающее конвективное течение в зависимости от соотношения величин равновесного и возмущающего градиентов температуры и типа локализованных возмущений гидродинамических и тепловых полей.
В цилиндрической системе координат (ег ,еф ,е г) осесимметричные уравнения Обербека-Буссинеска записываются как
'2 ( д дгуг д2уг ^
+ V
д г
дуг дуг ду
■ + + у
ді
д г
д г
ф др
-----?- = -^- + ч
Уф
г
дг дг д г2
ду
^ ф ді
■ + у.
ду
^ ф
д г
■+у.
ду
^ ф
д г
(
= \
д дгуф д\
дг дг д г:
дуг дуг ду
■ + Уг^г^ + у
= др
ді ' 'г дг '2 дг д г + У
(1 д дуг д2уг^
----г—- +-------
г д г д г д 2 2
V
+ Т,
дТ дТ дТ
-----+ у-------+ у ----= X
д і г д г 2 д г
(1 д дТ д 2Т}
-----г-----1-----
г д г д г д г 2
1 дгуг ду7 л —= 0.
г дг дг
Граничные условия, выражающие требования аналитичности скорости на оси возмущения [3] и механического равновесия на бесконечном удалении от нее, имеют вид
= А0;
_ дуг п дТ
г = 0: уг = уФ= о.-дТ = о.-дг
дР
г уг =уф = уг =0— = —№г.
і=0
(1)
Здесь уг ,уф,у7 - скорость жидкости; Т, Р - отклонения температуры
и давления, отнесенного к средней плотности р0, от их средних значений, соответствующих р0; в - коэффициент теплового расширения; g -модуль ускорения свободного падения; V - коэффициент кинематической вязкости жидкости.
Уравнения Обербека-Буссинеска допускают класс точных решений, обобщающий известное изотермическое решение [4, 5]
V = -
^ ( )
V- =^r-u(x,^, У<р =
у/Х
ди д х
Ж (т, х) + 2
(т, х) + 1у(т, х )^,
, вТ = -[©(т,х) + 20(т,х)] ,
(2)
Р = [р0 (х) - 42Р1 (т,х) - 22 2Р2 (х)] ,
где т = .
2^8^^I, х = г2, 2 = 3^V22 - безразмерные время
и пространственные переменные.
Представление (2) можно рассматривать как первые члены разложения гидродинамических величин и температуры в ряд Тейлора по поперечной координате 2, справедливое для тонкого слоя. Обрыв ряда оправдан точным удовлетворением (2) уравнениям Обербека-Буси-неска. В то же время это разложение не позволяет в полной мере учесть краевые условия на горизонтальных границах слоя.
Далее принимается, что V = Ж = & = рх = 0, тогда уравнения На-вье-Стокса в приближении Буссинеска редуцируются к системе
д 2 и
д и ди д и
д
д хдт+ и д х 2 д х д х 2 д х
V д х2 у
0-0,
д0 д0 0ди = 2 д I д0) дт + и дх дх Рг дх ^х дху,
д¥ д¥ „ д2У + и= 2 х—-
дт дх дх2
с граничными условиями, следующими из (1),
х = 0: и = V = 0, 0(О,О) = 0 0 =в А031
(3)
(4)
■го: •
ди 0.^ 0, 0, 0^0, = вА,зГ
V х V 8
д х ' ’Тх
Здесь Рг = у/х - число Прандтля; 0,, 0 0 - безразмерные равновесный и возмущенный градиенты температуры. Величины Р0 и Р2 определяются следующим образом:
ди и2 1 Г
Р° = их ~ 2х ~ 2 ]
V ] 2 д и
х ) х дт
Первые два уравнения (3) образуют замкнутую подсистему, которая в случае совпадения характерных вязкого и теплового диссипативных времен (Рг = 1) преобразованием
ди и а а
0 = а(т)дх + 0», -Й7 = Q,--Г, а(0) =
д х
4
-а0
(5)
сводится к одному уравнению с граничными условиями, следующими из (4) (последнее уравнение (3) и краевые условия к нему остаются неизменными):
д 2и + и д 2и ди ди = ^ д д 2и ^
д хдт+ и д х 2 д х д х д х
д х
х = 0: и = 0
ди , д х
ди
а д и 4 д х.
00 -0, ди ^ и
0 * • х^0, -=^0.
дх л/х
(6)
(7)
т=0 ^0
Безразмерный параметр а0/(00 -0,) характеризует соотношение
между начальным импульсом возмущения и создаваемой им (возмущением) дополнительной силой Архимеда на оси вихря.
В отсутствие фоновой стратификации 0, = 0 второе уравнение
(5) имеет решение а(т) = 4а0 (а0т + 4) 1, что позволяет ввести в (6), (7)
автомодельную переменную <^ = х(т + 4/ а0) 1 и переписать последнюю
систему в виде краевой задачи для обыкновенного дифференциального уравнения
2(,и") = (и -£)и"- и'и'- 2м', 0
Е,= 0: и = 0, и' = 4—2; ^^го: и' ^ 0, а п
->0.
(8)
(9)
Впредь штрихом обозначается производная по переменной £. Естественно полагать, что неизвестная V, подчиняющаяся последнему уравнению (3) с соответствующими граничными условиями (4), является функцией безразмерного времени т и переменной £. Тогда для на-
и
хождения азимутальной составляющей поля скорости имеем линейную начально-краевую задачу
т+0;)-з7=2^-Щ2+'£-ит-щ;,
V
V(o£) = Vo (); £= 0: V = 0; £^го: ^ 0,
решение которой можно искать в виде суперпозиции частных решений с разделяющимися переменными, что приводит к выражению
V (ТЛ) = Го | ехр
П+Р
{(к)
2 •> 2к
о
ЛП + +ЕГкКк (5)[т+аг
к Ф0
(10)
Постоянные Гк являются коэффициентами разложения начального возмущения V0 (5) по собственным функциям спектральной задачи (первое слагаемое в (10) соответствует собственному значению к = 0)
25 К1+\5- и (5)] Кк + ккк = °> К (0) = 0. 5т
К
(5) =
5-~° >15
2. Точные автомодельные решения
= 0.
(11)
Уравнение (8) обладает следующим решением (его стационарный изотермический аналог рассматривался в [6]):
и = а + Ь£ +
(Ь + 2)Ь = 0, (с - 12Л)с = 0, \[а + 4 - ё (Ь - 1)]с = 0, Ьс = 0. С учетом краевых условий (9) имеем
и
(5) = -15, а = -4, Ь = 0, с = 96, ё = 8, 00 = -3а2 <0.
5 + 8’
В результате, принимая во внимание (10), получаем следующие выражения для гидродинамических величин:
(12)
в(т +4-600 )2 ( + 8)2’
где кк удовлетворяют системе
Например, при к = 2, имеем: ^2 =£(£ + 8) 5ехр(-£/2).
Автомодельность вида х/ (т + т0) свойственна диффундирующим
течениям, что проявляется в дрейфе максимума радиальной скорости от оси вихря. В то же время имеет место радиальная адвекция по направлению к оси симметрии потока и подъем жидкости вверх от подстилающей непроницаемой плоскости г = 0 на фоне устойчивой стратификации возмущения при нулевой равновесной стратификации 0, = 0. Вдали от оси симметрии горизонтальное адвективное течение, аналогичное потенциальному вихрестоку, преобладает над вертикальным, так как (,уф ,Уг)~ (-г _1,Г 0 г _1, гг -2 ). В конечном счете тепловая
и вихревая диффузия берет верх над конвергентной адвекцией, и вихрь (12) со временем затухает. Следует отметить, что тепловое возмущение (12) локализовано «сильнее» гидродинамического и быстрее исчезает со временем.
Еще одним решением краевой задачи (8), (9), существующим при а 0 = 2 -00 , является (к этому типу относятся стационарные изотермические вихри Бюргерса [7] и Салливана [8] и их нестационарные аналоги [9, 10]) следующее:
С его помощью вновь удается построить пример точного решения уравнений конвекции, описывающего радиально-локализованные вихри:
Здесь базисные функции для разложения начального возмущения азимутальной скорости удовлетворяют задаче о собственных функциях и собственных значениях
Вихрь (13), подобно предыдущему случаю, затухает со временем вследствие тепловой и вихревой диффузии, преодолевающей конвергентную адвекцию. Вдали от оси горизонтальное течение по-прежнему представляет собой потенциальный вихресток, а температура и вертикальная составляющая скорости обращаются в ноль экспоненциально.
В силу сложности организации экспериментов количество работ, направленных на лабораторное изучение нестационарных вихрей, невелико. В [11] выполнено сравнение вихрей Бюргерса [7], Салливана [8] и Беллами-Найтса [10], относящихся к рассматриваемому здесь классу точных решений (2), с результатами лабораторного моделирования эволюции вихрей, образующихся вследствие локального отбора (аналог конвективной подъемной силы) со свободной поверхности жидкости, заполняющей вращающуюся кювету (аналог циркуляции на бесконечности). В результате обнаружено хорошее качественное, а в ряде случаев и количественное совпадение экспериментальных и теоретических профилей азимутальной скорости в последовательные моменты времени.
вТ = 40 0 2
вТ = 40 0 2
Таким образом, получены два примера точных автомодельных решений уравнений гидродинамики в приближении Обербека-Буси-неска, описывающих радиально-локализованные конвективные вихри, эволюционирующие в слое без температурной стратификации.
Работа выполнена в рамках Программы совместных проектов СО, УрО и ДВО РАН (проект 12-С-1-1006) и Российского фонда фундаментальных исследований (№ 12-01-00023 а)
Библиографический список
1. Гершуни Г.З., Жуховицкий Е.М. Конвективная устойчивость несжимаемой жидкости. - М.: Наука, 1972. - 392 с.
2. Богатырев Г.П. Возбуждение циклонического вихря или лабораторная модель тропического циклона // Письма в ЖЭТФ. - 1990. -Т. 51. - Вып. 11. - С. 557-559.
3. Гольдштик М. А., Штерн В.Н., Яворский Н.И. Вязкие течения с парадоксальными свойствами. - Новосибирск: Наука, 1989. - 336 с.
4. Prager S. Spiral flow in a stationary porous pipe // Phys. Fluids -1964. - Vol. 7. - P. 907-908.
5. Аристов С. Н. Стационарный цилиндрический вихрь в вязкой жидкости // Докл. АН. - 2001. - Т. 377, № 4. - С. 477-480.
6. Бурдэ Г.И. О движении жидкости вблизи растягивающегося кругового цилиндра // Прикл. математика и мех. - 1989. - Т. 53. -Вып. 4. - С. 343-345.
7. Burgers J.M. A mathematical model illustrating the theory of turbulence // Adv. Appl. Mech. - 1948. - No. 1. - P. 171-199.
8. Sullivan R.D. A two-cell solution of the Navier-Stokes equations // J. Aero/Space Sci. - 1959. - Vol. 26. - P. 767-768.
9. Rott N. On the viscous core of a line vortex // ZAMP. - 1958. -Vol. 9b. - P. 543-553.
10. Bellamy-Knights P.G. An unsteady two-cell vortex solution of the Navier-Stokes equations // J. Fluid Mech. - 1970. - Vol. 41. - Part 3. -P. 673-687.
11. Huang S.-L., Chen H.-S., Chu C.-C., Chang C.-C. On the transition process of a swirling vortex generated in a rotating thank // Exp. Fluids. -2008. - Vol. 45. - P. 267-282.
References
1. Gershuny G.Z., Zkuhovicky E.M. Konvektivnaya ustoychivost neszhimaemoy zhidkosti [The convective stability of the incompressible fluid]. Moscow: Nauka, 1972, 392 p.
2. Bogatyrev G.P. Excitation of cyclonic vortex or laboratory model of tropical cyclone. JEPT Lett, 1990, vol. 51, no. 11, pp. 630-633.
3. Goldshtik M.A., Shtern V.N., Yavorsky N.I. Vyazkie techeniya s paradoksalnimi svoyctvami [Viscose flows with paradoxical properties]. Novosibirsk: Nauka, 1989, 336 p.
4. Prager S. Spiral flow in a stationary porous pipe. Phys. Fluids, 1964, vol. 7, pp. 907-908.
5. Aristov S.N. A stationary cylindrical vortex in a viscous fluid. Dok-lady Physics, 2001, vol. 46, iss. 4, pp. 251-253.
6. Burde H.I. On the motion of fluid near a stretching circular cylinder. J. of Appl. Math. and Mech, 1989, vol. 53, iss. 2, pp. 271-273.
7. Burgers J.M. A mathematical model illustrating the theory of turbulence. Adv. Appl. Mech, 1948, no. 1, pp. 171-199.
8. Sullivan R.D. A two-cell solution of the Navier-Stokes equations. J. Aero/Space Sci., 1959, vol. 26, pp. 767-768.
9. Rott N. On the viscous core of a line vortex. ZAMP, 1958, vol. 9b, pp. 543 - 553.
10. Bellamy-Knights P.G. An unsteady two-cell vortex solution of the Navier-Stokes equations. J. Fluid Mech., 1970, vol. 41, part 3, pp. 673-687.
11. Huang S.-L., Chen H.-S., Chu C.-C., Chang C.-C. On the transition process of a swirling vortex generated in a rotating thank. Exp. Fluids, 2008, vol. 45, pp. 267-282.
Об авторах
Аристов Сергей Николаевич (Пермь, Россия) - доктор физикоматематических наук, главный научный сотрудник Института механики сплошных сред УрО РАН (614013 Россия, г. Пермь, ул. Академика Королева, 1, e-mail: [email protected]).
Князев Денис Вячеславович (Пермь, Россия) - кандидат физикоматематических наук, научный сотрудник ИМСС УрО РАН (614013, Россия, г. Пермь, ул. Академика Королева, 1, e-mail: [email protected]).
About the authors
Aristov Sergey Nikolaevich (Perm, Russian Federation) - Doctor of Physical and Mathematical Sciences, general scientist, ICMM UB RAS (1, Academic Korolev st., 614013, Perm, Russian Federation, e-mail: [email protected]).
Knyazev Denis Vyacheslavovich (Perm, Russian Federation) - Ph.D. in Physical and Mathematical Sciences, scientist, ICMM UB RAS (1, Academic Korolev st., 614013, Perm, Russian Federation, e-mail: [email protected]).
Получено 11.04.2013
