СЕМИНАР ПО МЕХАНИКЕ СПЛОШНЫХ СРЕД

Уважаемые коллеги!

В среду, 6 июня 2012 г., в кинозале Института механики МГУ в 12.00 состоится очередное заседание семинара по механике сплошных сред под руководством А.Г. Куликовского, В.П. Карликова и О.Э. Мельника.

Никитин Н.В.

НИИ механики МГУ им. М.В. Ломоносова, Москва

ЕЩЕ РАЗ О НИЖНЕМ КРИТИЧЕСКОМ ЧИСЛЕ РЕЙНОЛЬДСА В КРУГЛОЙ ТРУБЕ

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

Еще сам О. Рейнольдс заметил, что при малых Re (в районе Re=2000) течение в трубе является перемежающимся – участки “турбулентного” течения сменяются участками гладкого, ламинарного движения. Авторы работ [2,3] выделили два различающихся типа перемежающихся режимов: турбулентные пробки и турбулентные порывы. Первые увеличивают свою протяженность сносясь потоком и ведут к установившейся турбулентности, вторые либо затухают в процессе эволюции, либо делятся на несколько. С начала 2000-ых годов интерес к течению в круглой трубе вновь возродился в связи с обнаружением точных когерентных 3-мерных решений уравнений Навье-Стокса [4,5]. Эти решения имеют вид бегущих волн. Они неустойчивы, однако размерность отталкивающего многообразия в фазовом пространстве оказывается небольшой. Выдвинута гипотеза, что незатухающее турбулентное движение представляет собой случайное блуждание между различными когерентными состояниями.

Одновременно стали проводиться измерения и прямые расчеты характеристик турбулентных порывов – предвестников турбулентного течения в трубе. Обнаружено, что процесс затухания турбулентного порыва является Марковским (т.е. вероятность его затухания в каждый момент времени постоянна и не зависит от предыстории). Такой процесс аналогичен процессу распада радиоактивных веществ и характеризуется характерным временем жизни τ = τ (Re) (аналог периода полураспада радиоактивности). Характерное время жизни турбулентного порыва растет при увеличении числа Рейнольдса. Если при некотором Re оно обращается в бесконечность, то это число Рейнольдса и есть нижнее критическое. С точки зрения динамических систем этому соответствует переход от хаотического седла к хаотическому (странному) аттрактору. В первых исследованиях этого вопроса действительно был сделан вывод об обращении τ в бесконечность при некотором конечном Re [6-8].

Более тщательные измерения и вычисления [9-12] показали, однако, что τ растет при увеличении Re, оставаясь конечным. Даже был сделан вывод, что турбулентность в трубе принципиально затухающий процесс. Выходом из создавшегося противоречия может служить учет другого явления – стремления турбулентных порывов к делению. Оказалось [13], что процесс деления турбулентного порыва также является Марковским и характеризуется своим характерным временем, уменьшающимся при росте числа Рейнольдса. Число Рейнольдса Re=2040, при котором два конкурирующих процесса – процесс затухания и процесс деления уравновешивают друг друга является точкой статистического фазового перехода и может считаться нижним критическим числом Рейнольдса.

В лекции будет дан обзор описанных исследований и представлены некоторые визуализации соответствующих процессов, выполненные автором.

Работа выполнена при финансовой поддержке РФФИ, проект № 11-01-00088-а, вычисления проводились на суперкомпьютерном комплексе МГУ “Чебышев”.

.

1. O. Reynolds, Philos. Trans. R. Soc. London Ser. A 174, 935 (1883).
2. I. J. Wygnanski, F. H. Champagne, J. Fluid Mech. 59, 281 (1973).
3. I. Wygnanski, M. Sokolov, D. Fridman, J. Fluid Mech. 69, 283 (1975).
4. H. Faisst, B. Eckhardt, Phys, Rev. Lett. 91, 224502 (2003).
5. H. Wedin, R. R. Kerswell, J. Fluid Mech. 508, 333 (2004).
6. H. Faisst, B. Eckhardt, J. Fluid Mech. 504, 343 (2004).
7. J. Peixinho, T. Mullin, Phys. Rev. Lett. 96, 094501 (2006).
8. A. P. Willis, R. R. Kerswell, Phys. Rev. Lett. 98, 014501 (2007).
9. B. Hof, J. Westerweel, T. M. Schneider, B. Eckhardt, Nature 443, 59 (2006).
10. B. Hof, A. de Lozar, D. J. Kuik, J. Westerweel, Phys. Rev. Lett. 101, 214501 (2008).
11. D. Kuik, C. Poelma, J. Westerweel, J. Fluid Mech. 645, 529 (2010).
12. M. Avila, A. Willis, B. Hof, J. Fluid Mech. 646, 127 (2010).
13. K. Avila, D. Moxey, A. de Lozar, M. Avila, D. Barkley, B. Hof, Science 333, 192 (2011).