С. В. Борисёнок
КАЧЕСТВЕННАЯ МОДЕЛЬ СПЕКТРАЛЬНОЙ ФУНКЦИИ ВЯЗКОСТИ ДЛЯ ОДНОМЕРНОГО СПЕКТРА ЭНЕРГИИ РАЗВИТОЙ ИЗОТРОПНОЙ ЗАТУХАЮЩЕЙ ТУРБУЛЕНТНОСТИ
Обсуждается общая структура одномерного спектра энергии развитой затухающей изотропной турбулентности в рамках так называемых качественных моделей (англ. toy models). Уравнение спектрального баланса энергии представлено в безразмерной форме и замкнуто с помощью модифицированной схемы Гейзенберга. В результате вычислена спектральная функция полной вязкости и продемонстрирован ее универсальный характер.
104
Адекватное теоретическое описание систем с сильно развитой изотропной турбулентностью все еще остается одной из актуальных проблем теоретической физики. При этом спектры энергии — парный коррелятор скорости и функция переноса — играют при изучении турбулентности особую роль. Имеется достаточный экспериментальный материал, касающийся как их общего вида в инерционном интервале и за его пределами, так и деталей процесса переноса энергии по спектру. По этой причине любая модель турбулентности должна сравнивать полученные теоретически энергетические спектры с экспериментальными, что является важным аргументом в пользу ее применимости. Однако не все характеристики спектров на нынешнем этапе допускают экспериментальное исследование, так что желательно сравнивать новые теоретические модели с простыми качественными моделями развитой турбулентности, заведомо более «грубыми», то есть не учитывающими всех деталей явления, зато введенными на основе простых и легко формализуемых математически предположений. Такие полуфеноменологические модели, получившие английское название toy models, имеют в теории турбулентности (на фоне ее сложных нелинейных математических уравнений) весьма широкое применение.
В данной статье обсуждаются качественные модели (toy models) энергетических спектров трехмерной развитой изотропной затухающей турбулентности, которая (теоретически) должна описываться во всем многообразии ее свойств математической моделью, имеющей в основе уравнение Навье - Стокса. В принципе, из того же уравнения должны выводиться и спектры энергии поля скоростей, описывающие как общие процессы накачки и диссипации энергии в системе, так и детали переноса энергии по спектру. Практические трудности построения общей теории турбулентных спектров связаны со сложным существенно нелинейным характером уравнения Навье - Стокса. Поэтому все реально используемые модели, из которых получаются сопоставляемые с экспериментом результаты, основаны на той или иной схеме замыкания, в основу которой положена определенная физическая гипотеза. В случае конструирования одномерного спектра энергии E это означает, что соответствующая функция переноса T энергии по спектру должна быть выражена через сам спектр E из каких-либо внешних по отношению к уравнению Навье - Стокса соображений.
Подчеркнем еще раз, что все используемые практически схемы замыкания основаны не на точном решении уравнения Навье - Стокса (то есть нахождении общего вида функции переноса), а на какой-либо гипотезе замыкания, так что в этом смысле нет принципиальной разницы между построенной на простых соображениях качественной и реальной рабочей моделью турбулентности. Среди наиболее успешных схем замыкания выделим модель максимальной хаотичности поля скорости [1], введенную на основе модифицированной стохастической модели Уайлда.
Материал данной статьи относится только к трехмерной затухающей однородной изотропной турбулентности. Особо подчеркнем, что случай двумерной турбулентности здесь не обсуждается из-за более сложного характера законов сохранения (и числа сохраняющихся величин!), что приводит и к принципиально иному устройству энергетических спектров.
105
В статье описана общая структура одномерного энергетического спектра развитой турбулентности, изложены основы феноменологической теории Колмогорова, обосновывающей универсальность спектров энергии развитой турбулентности в инерционном интервале волновых чисел и в диссипационной области. Сама колмогоровская теория не замыкает уравнение спектрального баланса энергии, поэтому в статье обсуждается выбор схемы замыкания этого уравнения. Прежде чем перейти собственно к качественным моделям, мы приводим краткий обзор стохастических схем замыкания, знание основ которых существенно для понимания проблем, связанных с расчетом энергетических спектров. Затем обсуждаются качественные схемы замыкания уравнения баланса (toy models). За основу выбрана модифицированная качественная модель Гейзенберга как наиболее приемлемая для описания экспериментальных спектров турбулентности. Сконструированная схема замыкания пригодна для всего диапазона волновых чисел. В работе уравнение спектрального баланса представлено в безразмерной форме, после чего вычислена спектральная функция вязкости, описывающая детали переноса энергии по спектру и включающая кинематический коэффициент турбулентной вязкости.
Структура одномерного энергетического спектра развитой изотропной турбулентности
Одномерный энергетический спектр развитой изотропной турбулентности
Векторное поле скорости <Pi(x,t) несжимаемой вязкой жидкости описывается уравнением Навье - Стокса [2]
a << + (< ■д < = -—д P + vA<, (l)
Р
дополненным условием несжимаемости
д < = о.
Здесь р — плотность жидкости, P — давление, v — коэффициент кинематической вязкости, А — оператор Лапласа. Уравнение (l) представляет собой второй закон Ньютона, записанный в системе координат, движущейся вместе с жидкостью, а нахождение его решения относится к числу наиболее сложных задач классической гидродинамики. Однако существует возможность переформулировать задачу, перейдя от уравнения (l) к одномерному уравнению спектрального баланса в представлении Фурье волновых чисел к [3]. Ниже подобное описание принято как основное.
Для одномерного спектра энергии Е(к) развитой затухающей изотропной турбулентности уравнение спектрального баланса принимает вид
д 2
—Е (к, t) = -2vk2 Е (к, t) + T (к, t). (2)
Функция Т(к,г) называется одномерной функцией переноса энергии по спектру, причем на нее наложено условие
| (¡к ■ Т (к, г) = 0.
0
Введем полную энергию пульсаций
да
е(г) = | (к ■ Е (к, г) (3)
0
и мощность диссипации на единицу массы
да
ё (г) = -1 (к ■ 2ук2 Е (к, г). (4)
0
Тогда от уравнения (2) можно перейти к его проинтегрированной форме:
^ е(г) = -ё(г). (5)
дг
Соотношение (5) называется интегральной формой уравнения баланса энергии. В дальнейшем для краткости явная зависимость спектральных функций от времени г указываться не будет.
Феноменологическая теория Колмогорова
Преимущество описания системы на языке соотношения баланса (2) заключается в физической наглядности последнего. Согласно феноменологической теории Колмогорова, режим затухающей турбулентности характеризуется двумя основными масштабами: энергия исходно содержится в системе на малых волновых числах к (то есть в крупных вихрях) порядка т = Ь~ , где Ь — внешний (интегральный) масштаб турбулентности, и рассеивается в области больших к порядка к( (область диссипации). Между этими двумя областями лежит широкий инерционный интервал, в котором реализуется лишь перенос энергии по спектру. Подобное устройство спектра отвечает дроблению крупных вихрей и их последующему затуханию под действием вязкого трения. Колмогоров сформулировал две гипотезы:
1. Вид спектра энергии Е(к) в инерционном интервале и в области больших волновых чисел к определяется только двумя параметрами — уи е .
2. В самом инерционном интервале он определяется только накачкой, которая, в соответствии с интегральной формой уравнения баланса (5), равна мощности диссипации е .
Феноменологическая теория Колмогорова демонстрирует тем самым универсальность спектров энергии Е(к) развитой турбулентности в инерционном интервале волновых чисел и в диссипационной области. Это означает, что при-
107
веденные спектры Е(к) £ _2/3к//3 являются универсальной для всех турбулентных систем функцией безразмерного аргумента к/ка. Для широкого класса турбулентных систем экспериментальные данные подтверждают это заключение теории Колмогорова [3]. Что касается энергосодержащей области волновых чисел, то здесь приведенные спектры не ложатся на единую кривую, образуя широкий «веер». В этом отражается чувствительность длинноволновой части спектра к деталям накачки энергии в пульсационную составляющую, что обусловлено граничными условиями и взаимодействием с усредненным сечением. Однако для определенного класса течений можно рассчитывать на универсальное поведение и в данной области спектра. Впервые возможность подобной универсальности обсуждалась в работе [4]. Для этого должно выполняться условие локальной однородности, означающее малость интегрального масштаба турбулентности по сравнению с характерным расстоянием, на котором заметно меняются статистические характеристики пульсаций. Такое условие выполняется, в частности, для затухающей за решеткой турбулентности.
По теории Колмогорова, отношение двух основных масштабов турбулентности (т и ка) характеризуется числом Рейнольдса Яе — безразмерной комбинацией внешнего масштаба Ь, вязкости V и характерной скорости набегающего потока V, связанной с диссипацией соотношением £ = V3 / Ь . Тогда
к
^ = (Яе)3/4, (6)
т
и, таким образом, на основании выражений (3) и (4) приходим к равенству
К =
V V у
(7)
Иными словами, число Рейнольдса определяет ширину инерционного интервала. Развитой турбулентности соответствуют большие числа Рейнольдса Яе порядка 10 .
Спектр энергии Е(к) можно представить в виде
Е (к) = СК£2Пк ~5/3Е
к к
V т ка У
(8)
где Ск — безразмерная постоянная (константа Колмогорова), а Е — безразмерная скейлинговая функция от безразмерных же аргументов, на которую в инерционном интервале налагается требование Е = 1 (поскольку в нем спектр определяется только диссипацией).
Для нахождения спектра Е(к) из уравнения баланса в виде (2) или (5) необходимо указать схему замыкания, то есть выразить функцию Т(к), описывающую детали переноса энергии по спектру, через Е(к). (Феноменологическая теория Колмогорова сама по себе схемы замыкания не содержит!) Выбор подходящей модели для функции переноса позволяет рассчитать полный спектр энергии и тем самым выявить механизмы накачивания и диссипации энергии в
108
турбулентной системе. Проблема эта до конца не решена, большинство существующих на текущий момент схем замыкания имеют существенный недостаток — отсутствие универсальности во всем диапазоне волновых чисел, в том числе в области диссипации, полученные же универсальные решения, хорошо согласующиеся с экспериментом, основываются на сложных математических моделях (см., например, работу [5]).
Для любой качественной модели процесса достаточно указать простой вид подобного замыкания, по возможности доставляющий универсальный вид спектра энергии во всем диапазоне волновых чисел, то есть для вихрей любого масштаба.
Выбор схемы замыкания уравнения баланса энергии Стохастические модели турбулентности
Успех феноменологической теории Колмогорова для развитой турбулентности стимулировал построение статистических моделей. Применительно к модели изотропной турбулентности построение статистического ансамбля впервые было осуществлено Уайлдом в теперь уже классической работе [6]. Его
подход основывался на введении поперечной (д .f. = 0 ) случайной стохастической внешней силы f, действующей на единичную массу жидкости. Тогда уравнение (1) представляется равенством
1
дf + f д)р =--дP + vApi + f. (1а)
Р
Случайная внешняя сила в выражении (1а) моделирует стохастический характер поля скоростей ф1. Она считается гауссово распределенной, с нулевым средним и коррелятором в виде белого шума. Подобное устройство случайной силы (максимальная хаотичность ее поля) соответствует максимуму ее информационной энтропии. Случайная сила тем самым имитирует компенсирующую диссипацию накачку энергии в систему от взаимодействия с крупномасштабными пульсациями.
Предложенная Уайлдом стохастическая модель до настоящего времени является основной при изучении развитой турбулентности. Тем не менее, она встречается с определенными трудностями: искусственность введения случайной силы, сложность обобщения модели на неоднородный случай, техническая громоздкость осуществления принципа максимальной хаотичности. Кроме того, в ней парный коррелятор диссипации энергии
([дP, + д](x) [д,Pm + дP ](x h) * const • Ix - x 1
имеет неверное степенное поведение. Модель Уайлда доставляет значение / = 0, тогда как эксперимент демонстрирует / * 0,2.
Альтернативой модели уайлдовской выступает модель максимальной хаотичности [1], в основу которой положен спектральный баланс пульсацион-ной энергии. Основное предположение этой модели о характере эволюции системы состоит в том, что на начальной стадии затухания турбулентности скорость убывания энергии относительно мала, так что в каждый момент времени функция распределения пульсаций по-прежнему будет иметь квазистационарный вид.
Задача ставится аналогично и при рассмотрении турбулентного потока за решеткой при постоянной средней скорости потока V вдоль оси г. Такая задача стационарна, а роль времени I в ней играет расстояние 2 от решетки. Предполагая при этом локальную однородность системы, связанную с относительной малостью изменения вдоль г характеристик пульсаций на расстояниях порядка внешнего масштаба турбулентности Ь (иногда называемого также корреляционным радиусом), баланс энергии можно записать с заменой д( на Vдг.
Модель максимальной хаотичности позволяет написать замкнутую систему уравнений и решать ее, пользуясь условием максимума энтропии. При этом получить ответ в общем виде затруднительно. Обычно с помощью квантово-полевого по происхождению аппарата ренорм-группы и техники в-разложения ищется автомодельное решение в области волновых чисел из инерционного интервала и меньших.
Именно в модели максимальной хаотичности впервые теоретически рассчитан и обоснован спектр энергии в области всех волновых чисел [4].
Качественные схемы замыкания
Наиболее простой качественной схемой, замыкающей уравнение спектрального баланса (2), является модель, предложенная в 1948 году Гейзенбер-гом [3]. Проинтегрируем уравнение (2) по малым волновым числам к и представим поток энергии по спектру
к
](к) = | Т (д^д
0
(убывающий до нуля при к ^ да) в виде, аналогичном вкладу (4) от вязкости:
к к ]{к) = |Т(д)dq = -2^ (к)| д2Е(д)dq . (9)
0 0
Введенная в выражении (9) функция Ук(к) называется кинематическим коэффициентом турбулентной вязкости. В схеме Гейзенберга ее можно найти, не решая точно исходное уравнение (2), которое эквивалентно уравнению Навье - Стокса (1) для несжимаемой вязкой жидкости, а комбинируя тк(к) из соображений размерности из всех величин, от которых она может зависеть, то есть из Е(к) и к. Гейзенберг предложил простую модель
V (к) =
3
Я
ёд
(10)
с безразмерным подгоночным параметром у Легко убедиться, что расстановка пределов интегрирования в выражении (10) приводит к правильному асимптотическому поведению 1к(к), а именно к вязкому затуханию крупных вихрей. На самом деле функция у(к) зависит статистически от всех возмущений с разными волновыми числами к, так что она не может быть выражена через одни только Е(к) и к. В этом упрощенном предположении и проявляется главная особенность качественной модели турбулентного спектра.
Модель (9)—(10) доставляет верную асимптотику Е(к) в области диссипации. При этом левая часть проинтегрированного по к уравнения (9), связанная с накачкой энергии в систему, набирается, напротив, в энергосодержащей области (напомним, что для развитой турбулентности между этими областями лежит широкий инерционный интервал), поэтому с хорошей точностью
— [Е(д) ёя
= -е
(11)
Обозначение е, в отличие от (4), без черты сверху подчеркивает, что мы полагаем ее в нашем приближении константой. Соотношение (11) представляет собой закон сохранения энергии для установившегося турбулентного режима: заданная мощность энергии е, рассеиваемой на больших к, равна мощности энергии, накачиваемой в систему на малых к.
Собирая выражения (2) и (9)—(11) вместе, приходим к замкнутому интегральному уравнению на функцию Е(к):
к к е = 2 [у + ук (к ) ] ] я 2 Е (д) ёя = 21 д2 Е ( д)ёд •
0 0
Для его решения введем функцию
к
Н (к) = | д2Е(д)ёд ,
Е ( р )
ёр
(12)
(13)
поделим на Н(к) левую и правую части (5) (это допустимо, поскольку тривиальное решение Н(к) = 0 физически некорректно) и возьмем производную по к. Тогда получаем, что
Н (к) =
С 1 2
• к-4 + С
2
V е
(14)
У
где С — постоянная интегрирования. Соответственно,
Е (к) = ^ • к -5/3
С 1 2 ^ + Ск4
(15)
к
СО
0
Мы получили решение Баса—Чандрасекара (1949) [3] в модели Гейзенберга. Асимптотики этого решения суть к ~5/3 в инерционном интервале и к в области диссипации, что полностью соответствует эксперименту. Задача, однако, не решена до конца: энергосодержащая область, порождающая еще один размерный параметр — внешний масштаб Ь, в модели (10) не учтена. Следовательно, для описания всего спектра, включая малые волновые числа, необходимо расширение схемы.
Легко увидеть, что модель Гейзенберга для кинематического коэффициента турбулентной вязкости не пригодна для энергосодержащей области, хотя бы в силу неаналитичности выражения (10) при к ^ 0. Кроме того, в области малых волновых чисел в схему замыкания должен войти еще один параметр
внешний масштаб Ь (или, что удобнее, обратная к нему величина т = Ь"1). С учетом этого в работе [7] была предложена простая модификация модели Гейзенберга, а именно: в знаменателе интеграла в выражении (10) под знаком кор-
3 2 2 3/2
ня была сделана замена д ^ (д +т ) для продолжения свойств аналитичности по к вплоть до малых масштабов. В силу малости т, что следует из выражения (6), модифицированная модель Гейзенберга переходит асимптотически в выражение (10) на больших волновых числах. Кроме того, скейлинговая функция Г в равенство (8) была представлена в виде произведения
Г
'к кЛ
V т к*
= Г
к
V т у
к
V кл у
(16)
причем в инерционном интервале Г1 = Г2 =1, а Г2 совпадает со скейлинговой функцией Г модели Гейзенберга. В результате в работе [7] была найдена
У= 0,31.
Интересно отметить, что даже такая грубая модель приводит к одномерному спектру энергии, хорошо совпадающему с экспериментальным (несколько хуже, чем более точная схема замыкания [5]). Это доказывает практическую эффективность предложенной Гейзенбергом схемы замыкания, поскольку один из критериев приемлемой качественной модели - возможность ее расширения на весь диапазон волновых чисел.
Альтернативная схема замыкания (независимо от Гейзенберга) была предложена для функции переноса Коважным (1948) [3]:
№ = Еъ,2(к)кт
с соответствующей подгоночной константой ^соу.. Таким образом, Коважный предложил моделировать из размерных соображений не только собственно турбулентную добавку к вязкости Ук(к), а функцию у(к) целиком. Однако схема замыкания Коважного, как и более сложные модели Стюарта—Таунсенда (1951), Огуры—Миякоды (1953) и прочие (подробнее см. в работе [3]), не обладают должной степенью универсальности. Они не могут быть непрерывно продолжены на весь интервал волновых чисел, не доставляют правильные асимптотики, так что схему Гейзенберга и ее расширения следует считать лучшей качественной моделью из существующих.
Безразмерная форма уравнения спектрального баланса и вычисление спектральной функции полной вязкости
Представим уравнение баланса
к
е= 2 [у + ук (к) ]] q2 Е ^) dq (17)
0
в безразмерном виде. Для этого введем безразмерную функцию / — такую, что
е (к) = Гпк -5/7 (к / кл ),
и перейдем к безразмерной переменной х = к/к^ Учитывая соотношение (7), перепишем выражение (17) в форме
1 = 2В(х)\ у1П/(у^у, (18)
0
где введено обозначение для безразмерной спектральной функции полной вязкости:
В-1 + ^ . (19)
V
В частности, при В = 1 турбулентный вклад в вязкость отсутствует. Для модели Гейзенберга (10) имеем:
да
В(х) = 1 + г\ у~7/71/2 (у^у . (20)
х
Любая качественная модель в таком подходе принимает вид
В(х) = 1 + ГЪ(хЛх)1
где у— подгоночная постоянная, а Ъ^0 при х^-да. На функцию7х) наложены условия:
7(х)^Ск при х ^ 0 и7х) ^С1х"16/3 при х ^ да. (21)
Из сопоставления (15) с (8) видно, что
/ о Л2/3
4
С =
к -4/3 2/3
3 у
V 9У )
(22)
В модифицированной модели Гейзенберга спектральная функция вязкости принимает соответственно форму
(с 2 Л1/3 да ,-.1/2
В( х) = 1 +
Г Г (х) х
V 9 ) \ х5'6 (х2 + Ке—3/2 )3/4 ^ (23)
где в соответствии с выражением (8) сделана подстановка / (х) = Ск¥ (х) и учтено соотношение (22). Скейлинговая функция ^(х) найдена в работе [7], где она имеет вид
F ( х) = х ) х ); (х • Яе34 )П/3
х) =
(2,0 + (х • Яе34 )2 )• 1 +1,06-(1 + (х • Яе34 )2)
(24)
х) =
^3 V
V 2 С,
Отметим, что в работе [7] использовался альтернативный безразмерный аргумент X = к / т = х • Яе34 вместо х = к / кл, но для исследования функции спектральной вязкости он менее удобен. Таким образом, спектр В(х) содержит лишь один параметр - число Рейнольдса Яе, тогда как у, и соответственно СК, находятся путем численного решения уравнения спектрального баланса. Как показано в работе [7], у= 0,31.
Исследуем поведение нормированных спектральных функций при различных значениях чисел Рейнольдса. На рис. 1 приведен график спектральной функции для типичного в случае развитой турбулентности значения Яе = 100 при у= 0,31.
На рис. 2 представлено семейство образующих поверхность в координатах (х, Яе) спектральных кривых, соответствующих значению у= 0,31 и различным числам Рейнольдса.
Рис. 1. Спектральная функция В(х) при Яе = 100
Рис. 2. Зависимость спектральной функции В(х) от х и Яе
4/3
Из графиков спектральной функции В(х) видно, что она сводится, собственно, к вязкому вкладу при больших значениях х (в области диссипации), тогда как в энергосодержащей области В(х) >> 1, то есть в ней существенную роль играет именно кинематическая вязкость Ук(к). По мере возрастания числа Рейнольдса кинематическая вязкость становится все более существенной в области инерционного интервала (что хорошо видно на рис. 2). Таким образом, спектральная функция В(х) демонстрирует правильное поведение в зависимости от возрастания роли эффектов вязкости в турбулентном спектре.
Обратимся теперь к безразмерному потоку энергии по спектру. Для его получения выразим функциюу(к) через скейлинговую функцию ¥. Введем
у ( х) = Ж = -2 .Нк)
V
(25)
Из выражения (19) следует vK/v = В(х) - 1; кроме того, с учетом равенства (7),
V Н (к)
ск I Упр (у) Лу.
так что безразмерный поток имеет вид
х
3(х) = -2Ск (В(х) -1)\у1^(у)ф.
(26)
На рис. 3 этот поток построен для Re = 100 и у= 0,31. Напомним, что 3 = 0 (то есть В = 1) соответствует случай отсутствия турбулентной вязкости.
Поверхность 3(хДе) приведена на рис. 4 при у = 0,31. Для удобства визуального восприятия она построена для абсолютной величины 3 (так как сам поток 3 < 0).
Рис. 3. Безразмерный поток энергии по спектру 3(х) при Re = 100
Рис. 4. Зависимость потока 3(х) от х и Re
0
0
На последнем графике видно, что процесс переноса энергии наиболее интенсивно осуществляется на малых х, в энергосодержащей области. При этом набирающийся на больших числах Рейнольдса максимум абсолютной величины функции переноса довольно острый и устойчивый по волновому числу х (порядка 0,4) при различных числах Рейнольдса. Последнее означает, что, несмотря на возрастание роли вязкостных эффектов, кардинального изменения схемы переноса энергии по спектру не происходит.
В действительности, что не учтено в модифицированной качественной модели Гейзенберга, турбулентные возмущения набираются статистически на всех волновых числах к.
* * *
Из сказанного следует, что:
1. Существенная нелинейность модели развитой турбулентности делает невозможным аналитическое решение уравнения баланса энергии (2) во всем диапазоне волновых чисел. Однако условие развитости турбулентности (большие значения чисел Рейнольдса и, следовательно, большая разнесенность энер-госодержащей и диссипационной областей) позволяет не только получить правильные асимптотики в простой качественной модели, но и числено решить уравнение (2), используя закон сохранения энергии и метод анализа размерностей.
2. Качественная схема замыкания позволяет на основе найденной скей-линговой функции ^(х) вывести ряд других спектральных функций, таких как функция переноса Т(к) (для которой трудно предложить схему ее непосредственного экспериментального измерения) или более удобную, по нашему мнению, спектральную функцию полной вязкости В(х).
3. Спектральная функция полной вязкости В(х) в модифицированной модели Гейзенберга имеет универсальный характер (то есть охватывает весь диапазон волновых чисел к, параметрически зависит от числа Рейнольдса Яе) и обнаруживает правильное качественное поведение при повышении роли эффектов вязкости в спектре.
4. Любая реальная модель турбулентного энергетического спектра должна обладать теми же качественными свойствами спектров, что и модифицированная модель Гейзенберга, так что в условиях отсутствия удовлетворительных экспериментальных данных о полном спектре Е(х) (на практике измеряется его продольная составляющая) в качестве грубой оценки пригодности можно использовать спектры, полученные в рамках качественной модели.
БИБЛИОГРАФИЧЕСКИЕ ССЫЛКИ
1. Аджемян Л. Ц., НалимовМ. Ю. Теор. и мат. физ. Т. 91. 1992. № 2. С. 294-308.
2. Ландау Л. Д., Лифшиц Е. М. Гидродинамика. М., 1986.
3. Монин А. С., Яглом А. М. Статистическая гидромеханика. Т. 2. СПб., 1996.
4. Аджемян Л. Ц., Богданов С. Р., Сыщиков Ю. В. Вестник ЛГУ. № 10. 1982. С. 76-79.
5. Борисенок С. В. Вестник СПбГУ. Сер. 4. Вып. 2. 1994 (№ 11). С. 70-76.
6. WyldH. D. Ann. Phys. 14. 1961. P. 143-165.
7. Борисенок С. В., Новиков С. В. Вестник СПбГУ. Сер. 4. Вып. 3. 1999. (№ 18). С. 12-21.
S. Borisenok
TOY MODEL OF VISCOSITY SPECTRUM FUNCTION FOR ONE-DIMENSIONAL ENERGY SPECTRUM OF DEVELOPED ISOTROPIC DECAYING TURBULENCE
The general structure of the one-dimensional energy spectrum for developed decaying isotropic turbulence in the frame of the so-called «toy models» is discussed. The energy balance equation is represented in a dimensionless form with the modified Heisenberg closure scheme. As a result the full viscosity spectrum function is calculated and its universal character is demonstrated.