CC BY
35
УДК 532.517,4.:536.35
МОДЕЛИРОВАНИЕ КАСКАДНЫХ ПРОЦЕССОВ В КОНВЕКТИВНОЙ ТУРБУЛЕНТНОСТИ ПРИ ЭКСТРЕМАЛЬНЫХ ЗНАЧЕНИЯХ ЧИСЛА ПРАНДТЛЯ
Ложкин С.А., Фрик П.Г. ( Пермь )
Abstract
A shell model of fully developed turbulent convection is considered. In contrast to traditional shell models this model includes the nonlocal interactions that become of vitaI importance in the case of extreme value of Prandtl number. We consider cascade processes in convective turbulence under very low and very large Prandtl number. New kinds of spectral intervals are described.
Каскадные модели появились как попытка описать процессы спектрального переноса в развитой турбулентности с помощью небольшого числа переменных U„. каждая из которых является коллективной характеристикой амплитуд
пульсаций поля скорости в целой октаве волновых чисел 2" <!&!< 2"н1 .Уравнения для ¿/„должны воспроизводить "базовые" свойства уравнений Навье-Стокса, а именно, иметь по меньшей мере те же интегралы движения и тот же вид нелинейности. Каскадные модели использовались рядом авторов начиная с 70-х годов [1,2.5.7,9,13] . но особенно активно их стали исследовать в 90-х [14,18,19,21] после того, как было показано, что они воспроизводят поведение структурных функций высших порядков в реальной турбулентности [23] . Наиболее активно исследуется модель, называемая GOY (Gledzer-Ohkitani-Yamada), которая записывается для комплексных переменных Uп в виде
d,u, - ¡Кiauixui2 + bu;jj;+1 + cu:_tu;_2; - Re ' k]un +/„, (i)
где Re - число Рейнольдса, кк=к02" - волновое число, а коэффициенты
о,ft,(.выбираются так, чтобы обеспечить сохранение энергии Е = ^ |i7„|2 в пределе Re —> œ, fn —> 0 (например, а = 1, b = -1/4, с = -1/8 ).
Следует указать два основных качественных отличия каскадных моделей от уравнений Навье-Стокса. Во-первых, это полная потеря информации о пространственном распределении исследуемых случайных полей. Во-вторых, уравнения типа (1) описывают только локальные (в /с-пространстве) взаимодействия. Первый недостаток является следствием низкой размерности динамической системы, второго можно избежать и оставаясь в рамках каскадных моделей. Впервые модель с нелокальными взаимодействиями была построена в работе [7] для двумерной турбулентности. Модель для конвективной турбулентности описана в [8]. Модификация этой модели используется в данной работе для рассмотрения конвективной турбулентности при экстремальных значениях числа Прандтля. когда учет нелокальных взаимодействий становится принципиальным.
Турбулентная конвекция. Рассматриваются течения, описываемые в рамках приближения Буссинеска для термогравитационной конвекции несжимаемой жидкости. Уравнения движения запишем в безразмерной форме
д,и+(и-Ч)и = -УР + СГ1/2\и+еТ, (2)
д,Т + й -ЧТ = а~'Сг~У2АТ, (3)
У-Й = 0, (4)
где й - скорость, Р - давление, Т - температура, е - единичный вектор вдоль вертикальной оси, Сг = gfiL}T'v'2 - число Грассгофа.о = у/х - число Прандтля. V -кинематическая вязкость, / - температуропроводность. В качестве единицы длины выбран макромасштаб Ь, единицы температуры - характерная для этого масштаба разность температуры Г*, единицы скорости V = (^^ЬТ")У1 и единицы времени г=/УК. При выбранной единице скорости число Грассгофа просто связано с числом Рейнольдса Сг = V21}\'2 - Яе‘.
В случае, когда единственной причиной движения являются неоднородности температуры, а малые значения вязкости и температуропроводности позволяют развиться инерционному интервалу, простые соображения размерности согласно теории Колмогорова позволяют предсказать спектральные законы для пульсаций скорости и температуры. Спектральный поток кинетической энергии, который в случае колмогоровской турбулентности пронизывает все масштабы и равен скорости диссипации энергии е, в конвективной турбулентности перестает быть константой, гак как силы плавучести привносят энергию в каждый масштаб. Зато сохраняющейся величиной становится поток энергии пульсаций температуры е; , который на
масштабе / есть бГ,2/?, (характерное время переноса /, - //81/). Постоянство ст дает для характерных пульсаций скорости 8У, и пульсаций температуры 87", на расстоянии 1(Ь» I» 1У,1Х)
ЪТ2Щ ~ 1. (5)
Динамика конвективного интервала определяется помясштабным балансом нелинейного члена и члена, описывающего архимедовы силы в уравнении для скорости, что позволяет написать соотношение
5^2//~57;. (6)
Комбинируя (6) и (5), легко получить для инерционного интервала оценки
ЬУ,~1Ъ!\ 67(7)
соответствующие спектральным законам
Еу{к)~к^\ (8)
Ет{к)~к-7/\ (9)
Эти законы, впервые полученные Обуховым [6] и Болджиано [5] для стратифицированной среды, многократно воспроизводились и обсуждались на различном уровне строгости [3,4,16,22]. Отметим, что исключение из определяющих параметров спектрального потока энергии £ делает результаты (8)-(9) не зависящими от размерности пространства. В обычной же гидродинамической турбулентности переход от трехмерного пространства к двухмерному связан с появлением дополнительного интеграла движения (энстрофии) и приводит к качественному изменению характера каскадных процессов.
Появления инерционного интервала с законами (8-9) можно ожидать, когда масштабы вязкой диссипации и тепловой диффузии имеют один порядок малости,
т.е. число Прандтля о~1. Однако и в пределе бесконечно малых V и % интервал (8-9) не растет неограниченно. Поступление кинетической энергии за счет архимедовых сил падает по мере уменьшения масштаба как /4' и на масштабе I = (масштабе Болджиано) возникает ситуация, когда обычный
каскадный перенос становится эффективнее термогравитационной конвекции и в трехмерном течении устанавливается колмогоровский спектр. Цель данной работы состоит в рассмо фении специфики каскадных процессов вблизи масштаба Болджиано в двумерной турбулентности, а также каскадных процессов при очень низких и очень высоких значениях числа Прандтля.
Каскадная модель. Для двумерной турбулентной конвекции каскадная модель была построена в работе [8] и имела вид
4и„ = £ 7; т, и-и> ~Яе~11и» + /Г»0»■ (10)
т.!
4®, = - к'-(аЯе)-1@„, (11)
т.!
где /%, =^,2", Тпт1 =2* Т0Нпт1 = 2Л'Я0а значения элементов для центральных частей матриц Т0л1, и Ная) приведены в таблицах. Структура матриц следует из разбиения пространства волновых векторов на октавы и из требования сохранения кинетической энергии Еу = Х„1^»Г ’ энстрофии
I |2 . I |2
О = ^ \2пи,г\ и энергии пульсаций температуры ЕГ - ¿,и|2ч©и| •
Т
1 \ т -4 -3 -2 -1 0 1 2 3
4 0.155
3 0.242
2 0.431
1 -0.0088 -0.0257 -0.0796 -0.269
0
-1 0.0032 0.0096 0.0269
.2 і 1
/ \ т -4 -3 -2 -1 0 1 2 3
4 -0.0537
3 -0.0941 -1.493
і -0.125 -0.720 -0.153
1 1 -0.0058 -0.0145 -0.0374 -0.0996 -0.221 -0.365 -0.145
о!
-і 0.0018 0.00468 0.0125 0.0277 0.045'Г 0.0181
0.00196 о.оиз 0.00239
-3 0.000184 0.00291 0.00030
-4 0.000013 0.00073 0.00004
В данной работе используется модифицированная модель вида (10-11). Во-первых, рассматриваются комплексные переменные, использование которых существенно снижает время интегрирования, необходимое для получения устойчивых статистических характеристик [17]. Во-вторых, в матрице НпШ
рассматриваются только члены, описывающие генерацию неоднородностей температуры крупномасштабным полем скорости (строки / = ±1, т < 0), и диагонали т-п и 1-т, которые очевидно доминируют над соответствующими боковыми столбцами. Тогда, с учетом связей между элементами матрицы, следующих из законов сохранения, можно записать
- ;2 ^Ta-j nijU -UK_jUn+i
- k2n Re~' U„+Fn@„,
-8íC/41©;+1)+
(13)
+- 2yí/;+y©;+I)) - Rey'®,,,
Параметр У фиксирует наиболее далекие взаимодействия (при J =1 система возвращается к стандартному виду каскадных уравнений, описывающих только локальные взаимодействия).
Умеренные числа Прандтля (ст~ 1). Спектральные законы вида (8-9) получались в моделях как двух-, так и трехмерной конвективной турбулентности при решении каскадных уравнений для относительно небольшого диапазона масштабов (порядка десяти) [8,9,3]. Рассмотрим эволюцию спектров двумерной турбулентной конвекции при очень больших числах Грассгофа, когда большой интервал значений волновых чисел позволяет проследить за формированием спектров по обе стороны от масштаба Болджиано. Система уравнений (12-13) для случая, когда число Прандтля равно единице, а число Грассгофа Сг = 1024 (что соответствует Яе = 1012), интегрировалась методом Рунге-Кутта четвертого порядка с постоянным шагом по времени для 0 < п < 30. Равномерный нагрев на макромасштабе моделировался путем поддержания стационарного значения модуля переменной |©0| = 1.
В отличие от трехмерного случая, в двумерной гидродинамической турбулентности существование инерционного интервала с прямым каскадом энергии невозможно. Это обстоятельство препятствует установлению стационарного распределения энергии по спектру. Процесс передачи энергии к мелкомасштабному движению блокируется на масштабе /„, вправо от которого формируется инерционный интервал переноса энстрофии. Влево от /„ развивается интервал обратного переноса энергии к крупным масштабам со спектральным законом "-5/3", причем граница этого интервала продвигается ьлево по мере накопления системой энергии.
Стационарной ситуации удается добиться путем введения дополнительной диссипации кинетической энергии на больших масштабах (в уравнение для Vп дописывается член вида -уип , так называемое линейное трение, обычно используемое и при прямых численных экспериментах с двумерной турбулентностью). На рис.1 показаны осредненные по времени значения энергии пульсаций скорости и температуры в отдельных октавах Е,-„ и ЕТп. Проведены
линии, соответствующие степенным законам для спектров Еу(к) и Ет(к). При рассматривании рисунков следует помнить, что показатель степени для величины
£„(£„) на единицу меньше, чем для самого спектра Е(к), что связано с принятым делением оси волновых векторов на октавы. Границы различных интервалов более четко выражены в спектре пульсаций скорости. В спектре пульсаций температуры переходы размытые и степенные участки не столь ярко выражены.
1с^ Л"„/ 1с^2
О 5 И) 1 Г> 20 25 п
Рис. 1
Любопытно отметить, что описанный численный эксперимент есть первая удачная попытка получения стационарного интервала обратного каскада энергии в каскадных моделях двумерной турбулентности. Инверсный каскад с законом, близким к "-5/3", регистрировался в каскадных моделях двумерной гидродинамической турбулентности в узком интервале волновых чисел [7], но все попытки получения стационарного интервала с постоянным потоком энергии не дали положительного результата - после накопления некоторой энергии в крупных масштабах (малые п) устанавливался режим равнораспределения энергии между переменными Un с соответствующим спектром Ег(к)~-к~' [10,17] Рассмотренная нами ситуация отличается способом подвода энергии, а именно, отсутствием какой-либс. искусственной силы. возбуждающей турбулентность на промежуточных масштабах.
Малые числа Прандтля (о<< 1) приводят к возникновению инерционнодиффузионного интервала в спектре пульсаций температуры. Он возникает в масштабах, на которых сохраняется обычный инерционный интервал в поле скорости, но существенна тепловая диффузия. Из сопоставления соответствующих членов уравнения (4) bVfiTeL~' ~ ЪТ,ГХ и колмогоровской оценки для пульсаций скорости в инерционном интервале получается [12]
E.ik; А1’1. (14)
В этом случае результат зависит, естественно, от динамики инерционного интервала и в случае рассмотрения двухмерной турбулентности, в которой реализуется инерционный интервал переноса энстрофии со спектральным законом Е(к) - к~' . мы придем к закону [8]
Е, iк ) '■ к ’. (15)
При cto.il быстром затухании энергии пульсаций трудно рассчитывать на формирование протяженного интервала. Это подтверждают и результаты
численного счета для случая о = 10 8 . приведенные на рис.2. где не удается выделить интервала с постоянным степенным законом.
log Ап / log 2
0 5 10 15 20 25
Рис. 2
На рисунке даны результаты счета с различными значениями параметра J. В конвективном интервале отличие невелико, так как здесь доминируют локальные взаимодействия полей скорости и температуры. Различия хорошо видны в мелкомасштабной части спектра, где формируется инерционный интервал переноса энстрофии. Известно, что каскадные модели испытывают проблемы с описанием каскада энстрофии, выражающиеся в том, что поток энстрофии в них слаб в сравнении с пульсациями энстрофии в отдельном масштабе и падает с ростом п, а спектры энергии не следуют единому степенному закону [10,17]. Увеличение J усиливает поток энстрофии и приводит к появлению протяженного интервала, в котором спектр кинетической энергии следует степенному закону с наклоном Е(к)~ к“10/3. Интересно отметить, что именно такой наклон спектра был получен при исследовании двумерной турбулентности с помощью иерархической модели, в которой число переменных растет как 22" по мере роста волнового числа кп- 2”, что позволяет в отличие от каскадных моделей учитывать и пространственную неоднородность турбулентного течения [11].
Большие числа Прандтля (а» I) способствуют формированию вязкоконвективного интервала, в котором соответствующие масштабы поля скорости подавлены вязкостью, но остается спектральный поток пульсаций температуры, поддерживаемый лишь крупномасштабным полем скорости. Поскольку диффузия тепла происходит на существенно меньи их масштабах, поток энергии пульсаций температуры остается постоянным, но в оценке (5) характерное время переноса определяется крупномасштабными пульсациями скорости и может считаться для этого интервала постоянным.Это дает оценку 57J2 ~ corist, что равносильно спектральному закону, впервые полученному Бэтчелором [12],
Ет(к) ~ к~'. (16)
Между конвективным (обуховским) и вязко-конвективным интервалами можно ожидать появления промежуточного интервала, в котором вязкий член становится весомее нелинейного, но остается существенной сила плавучести. Тогда баланс архимедовых и вязких сил приводит к соотношению
&8Т,~8У,Г2.
которое вместе с (5) однозначно определяет спектральные законы
Е„(к)~к-и\ Ет(к)~к~1,ъ. (18)
Очевидно, что такой интервал может установиться лишь в узком диапазоне волновых чисел, так как быстрое ослабление пульсаций скорости стимулирует развитие вязко-конвективного интервала (16). Однако баланс (17) может соблюдаться и в вязко-конвективном интервале (когда доминирует перемешивание крупномасштабным голем скорости). Тогда спектр пульсаций температуры будет следовать закону ' -1" (16), а соотношение (17) приводит к спектру
Ег(к)~к-\ (20)
Результаты вычислений для случая о = 106, Ке = 10 представлены на рис.З и показывают, что интервал, в котором устанавливаются законы (16-19), может быть достаточно протяженным. Можно видеть, что увеличение J приводит к растяжению интервала (16), но практически не влияет на распределение энергии пульсаций скорости.
1од Л’,, / <0^2
1П
L\-(J = 1) EtU = 1) Ev(J 5) br(J = Ь)
in j-
i
I
-20 j-
'•--- i--5
К
4
Pnc. ‘A
Работа выполнена при финансовой поддержке Российского фонда
фундаментальных исследований (грант 94-01-0095!-а).
Библиографический список
1. Гдедзер Е.Б Система гидродинамического типа, допускающая два квадратичных интеграла движения //ДАН СССР. 1973. Т. 209. № 5.
2. Зимин В.Д. Иерархическая модель турбулентности // Изв. АН СССР. Физика атмосферь! и океана. 1981. Т. 17, № 12. С. 1265-1273.
3. Зимин В.Д., Фрик П.Г. Турбулентная конвекция. М.: Наука, 1988. 175с.
4. Мокин A.C., Яглом А.М. Статистическая гидромеханика. М.: Наука, 1967. Ч. 2. 720с
5. Деснянский В.H.. Новиков Е.А. Моделирование каскадных процессов в
турбулентных течениях // Прикл. математика и механика. 1974. Т. 38, № 3. С.
507-513.
6. Обухов А.М.//ДАН СССР. 1959. Т. 125, №4, С. 1246.
7. Фрик Г1.Г. Иерархическая модель двумерной турбулентности // Магнитная
гидродинамика. 1983 Т. 19. № 1. С. 60-66.
8. Фрик П.Г. Моделирование каскадных процессов в двумерной турбулентной конвекции // Прикл. мех. и тех. физ. 1986. Т. 27, № 2. С. 71-79.
9. Шайдурова Т.Е. Иерархическая модель турбулентной конвекции. Препринт. Свердловск: ИМСС, 1986. 32с.
10.Aurell Е., Bofetta G., Crisanti A., Frick P., Paladin G., Vuipiani A. Statistical mechanics of shell models for 2D-turbulence // Phys. Rev. E. 1994. Vol. 50. P. 47US-4715.
11.Aurell E., Frick P., Shaidurov V. Hierarchical tree-model of 2D-t urbuler.ce // Physica D. 1994. Vol. 72. P. 95-109.
12.Batchelor G.K., Howells I.D., Townsend A,A. // J.Fluid Mech. 1959. Vol. 5, N
1.P.134.
13.Bell T.L., Nelkin M. Time-dependent scaling relations and cascad model of turbulence//J. Fluid Mech. 1978. Vol. 88, N 2. P. 369-393.
14.Biferale L.. Lambert R,, Lima R. & Paladin G. Transition tc chaos in a shell model of turbulence// Physica D. 1995. Vol. 80, P. 105.
15.Bolgiano R. //J.Geophys.Res. 1959. Vol. 46, P. 2226.
16.Brandenburg A. Energy spectra in a model for convective turbulence //
Phys.Rev.Let. 1992. Vol. 69, N 4. P. 605-608.
17.Frick P. & Aurell E. On the spectral laws in shell-models of 2D-turbulence II Europhys. Lett. 1993. Vol. 24, N 9. P. 725-730.
18.Frick P.. Dubrulle B., Babiano A. Scale invariance in a class of shell-models // Phys. Rev. E. 1995. Vol. 51, N 6, P. 5582-5593.
19.Jensen M.H., Paladin G. & Vuipiani A. Ashell model for turbulent advection of
passive scalar fields//Phys. Rev. A. 1991. Vol. 43, P. 798.
20.Kraichnan R.H. Inertial rangesin two-dimensional turbulence // Phys. of Fluids. 1967. Vol. 10. N 7. P. 1417.
21.Pisarenko D., Biferale L., Courvoisier D.. Frisch U. & Vergassola M. Further resuits on multifractality in shell models // Phys. Fluids A. 1993. Vol. 5, N 10. P. 2533-2538
22.Yakhot V. // Phys. Rev. Lett. 1992. Vol. 69. P. 769.
23 Yamada M. & Ohkitani K. The constant of motion and inertial subrange spectrum in fully-developed model turbulence// Phys. Lett. 1988. Vol. A134. P. ¡65-169.
