Рассеяние звука на турбулентных флуктуациях давления и энтропии Текст научной статьи по специальности «Физика»

Рассеяние звука на турбулентных флуктуациях давления и энтропии

Е. В. Юшков, В. П. Юшков

Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра физики атмосферы. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.

E-mail: yushkov@phys.msu.ru

Статья поступила 10.02.2011, подписана в печать 29.06.2011

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

Ключевые слова: турбулентность, энтропия, давление, звуковые волны, сжимаемость, бароклинность.

УДК: 504.35. PACS: 47.27.nb, 43.28.+h, 92.60.hk.

Введение

Турбулентность в пограничном слое атмосферы изучена достаточно подробно. Классические монографии [1-5], сотни научных статей, лабораторное и численное моделирование, множество натурных экспериментов, кажется, охватили все вопросы, которые могут быть исследованы. Однако новые измерения порождают и новые вопросы, которые можно обобщить следующим образом: какие приближения или предположения, заложенные в основание теории турбулентности, выполняются в наименьшей степени в новых задачах и к чему приведет замена этих приближений на, возможно, более адекватные. Настоящая работа в значительной степени повторяет вывод нескольких классических формул, лишь добавляя к рассмотрению новые члены, которые могут быть значимы при определенных условиях или которые дают уточняющие поправки к классическим формулам. Одновременно эти новые поправки могут изменять интерпретацию экспериментальных результатов и потому порождать необходимость дополнительных исследований.

Хорошо известно, что A.M. Обухов еще в 1949 г. в работе [7] рассмотрел вопрос о турбулентных флуктуациях давления в несжимаемой жидкости при условии р = const. Однако в атмосфере приближение постоянной плотности не является адекватным, в отличие от приближения несжимаемости (div« = 0) для достаточно медленных турбулентных флуктуаций, поэтому для атмосферы в формулу Обухова следует внести коррективы. Несмотря на эту очень важную работу, до сих пор традиционно полагается, что турбулентные флуктуации давления в атмосфере малы, а флуктуации плотности связаны лишь с флуктуациями температуры. Так, например, в 1953 г. В. И. Татарским была получена формула, широко используемая в настоящее время и описывающая рассеяние звука в турбулентной атмосфере [5], и в ней отсутствуют турбулентные

флуктуации давления. Приближения, использованные Татарским, являются общепринятыми в настоящее время, и большинство последующих работ использует эти приближения (см., например, книгу [6]).

Повторение вывода формулы Татарского показывает, что использованные приближения проще для анализа и, возможно, в большинстве практических задач флуктуации давления действительно дают малый вклад в формулу рассеяния звуковых волн «слышимого» диапазона. Однако этот дополнительный анализ оказывается полезен в методическом плане и открывает новые механизмы, влияющие на рассеяние акустических волн с большей длиной волны — сотни метров и более (акустический инфразвук). Экспериментальные измерения индикатрисы рассеяния звуковых волн в турбулентной атмосфере, проведенные М. А. Каллистар-товой [9], подтвердили формулу Татарского, однако в иных метеорологических условиях, например ночью, особенно при сильной устойчивой стратификации, эти дополнительные члены могут давать и более значимый вклад в рассеяние инфразвуковых волн даже и не очень низкой частоты. К тому же точность измерений при проведении такого эксперимента в настоящее время может быть значительно увеличена.

Принцип рассеяния звуковых волн на турбулентных неоднородностях лежит в основе, например, акустического зондирования атмосферного пограничного слоя (АПС), которое активно развивается уже более 30 лет [11]. Известно, что интенсивность рассеяния, а также высота слоя перемешивания, измеряемые акустическим способом, особенно в условиях устойчивой стратификации, не всегда совпадают с данными других измерений как интенсивности [12], так и высоты пограничного слоя [13]. Для объяснения этих различий привлекаются дополнительные предположения, справедливость которых порой трудно проверить. В этом смысле настоящая работа указывает на еще одну возможную причину.

1. Флуктуации давления и энтропии в индикатрисе рассеяния звуковых волн

Проделаем еще раз вывод формулы Татарского, добавляя те члены, которые отсутствовали в исходном анализе в силу принятых упрощений, но которые могут иметь тот же порядок, по крайней мере теоретически, что и оставленные члены уравнений. Распространение в турбулентной среде звуковых волн с малой длиной волны (в инерционном интервале турбулентности) описывается упрощенными уравнениями гидродинамики [14]. Пренебрегая пока диссипативными процессами, влиянием силы тяжести и сил плавучести:

dV 1

dp at

div(pV)=0.

(1) (2)

Будем рассматривать скорость, плотность, давление и температуру как сумму средней, турбулентной и акустической составляющих: Vi = Vqi + щ + , где £г- — скорость адиабатических движений (звуковая), щ — скорость турбулентных флуктуаций, V0i — среднее течение, например в пограничном слое (i — компоненты скоростей); р = ро +ра +р', ро — гидростатическое давление, ра — акустическое давление, р' — турбулентная составляющая. Для упрощения выкладок, поскольку нас интересует в первую очередь влияние турбулентных флуктуаций давления на рассеяние звука, положим V0i =О (отсутствие доплеровского сдвига в рассеянии) и ро = const, Т0 = const (рассеяние в небольшом объеме, нет зависимости интенсивности рассеяния от температурной стратификации атмосферы — рефракции акустических волн).

Считая воздух идеальным газом, запишем выражение для энтропии: S = Св In р - Ср In р, где Св и Ср — удельные теплоемкости при постоянном объеме и давлении (Ср/Св =7). Уравнение адиабатичности для акустической компоненты

(3)

замыкает систему уравнений (1), (2), если считать поля турбулентных флуктуаций заданными. Из (3) следует

din р 1 d\np

dt 7 dt

1 dP А- Т/

7p dt

(4)

(5)

(6)

Выразим теперь термодинамическую температуру Т через потенциальную температуру в и давление, поскольку потенциальная температура не меняется в звуковой волне и ва = 0, а в = во + в'. Тогда

К»Т = 7 I {йУ> где х = а 4 = Или

Ямг = |(1 + Г) где Г = в'/во- Аналогично

обозначим Р' = р'/ро и Р = ра/ро, в соответствии с обозначениями Татарского [5]. Отметим сразу, что Т в этой части работы описывает относительные

флуктуации потенциальной, а не термодинамической температуры.

Тогда из первого уравнения системы (5) получаем

d(pa+p') , , , ,w дщ

Jt = -lipo+pa +Р )(q- '

aeu

дх/'

(7)

или

д(Р + Р')

dt

(Щ + Ъ)

д(Р + Р')

= ^7

Ч I SI/ о

axi

дщ

дх,

■1

dxi

7 (Р + Р')

дщ dxi

дЪ

dxi

(8)

Линеаризуем уравнения относительно акустических переменных и отбросим члены, содержащие лишь турбулентные величины, так как они будут скомпенсированы в силу тех же уравнений движения (без акустических составляющих). В следующем разделе эти члены будут рассмотрены подробнее. Отметим лишь, что уравнения для турбулентных компонент замыкаются дополнительным условием сНуи = 0. Тогда

д(Р + Р') дР

■6

дР'

dCi

dt ^Ut dxi 1 dxi ^'dxi '' dxi'

.7P

№

(9)

Будем считать, что акустические величины зависят от времени посредством множителя ехр(—го>0 (волны в комплексном представлении), и период этих волн значительно меньше характерного времени изменения турбулентных флуктуаций такого же пространственного масштаба. Это означает, что в сумме акустической и турбулентной составляющих турбулентная часть описывает флуктуации, усредненные за время много больше периода акустических колебаний (для акустических колебаний инерционного масштаба длин волн

v ~ 1 кГц), поэтому дР

iojP = Щ-£— ОХ1

дР'

д t

•С

. Тогда

дР' _ ¡ю

о7Г РЩ-

(10)

Поступая аналогичным образом со вторым уравнением системы (6), получаем

dt

с,

41

7

Г)( 1

дхк

Р + Р')*-1

д(1+Р + Р')

дх.

(11)

Заносим скобку с давлением под производную и используем волновое представление для акустической компоненты. Получаем

дхк дхк) 7х дх1

(12)

а линеаризуя это уравнение относительно турбулентных флуктуаций, получим

-ги>£ + ик

дхь

■ik

дщ

дхь

— (1+7")

7>i

д{\ +нР + нР')

дх.

или

• i- "St , J- , Co I C0TI

ldxk

дщ

'dxk

7 dxi 7

,dP_

дх i'

(13)

(14)

Выше мы учли, что члены, содержащие только турбулентные компоненты, сократятся в силу уравнений движения (если акустические возмущения по амплитуде значительно меньше турбулентных).

Возьмем дивергенцию от (14):

■

Ш1Г~ = ОХ!

д_

дXI

Щ

дхь

дщ дхь

с2 с2 В

7 дх1

Т

Умножим (10) на ш: дР

дР'

дх,

■ Ш'уР'

дх.

дР

дх1/ (15)

(16)

—ь?Р = --Ь /^ц/—--Ь 7

оXI ОХ{

и подставим (14), отбрасывая величины второго порядка малости по и, Т', Р' (произведения, содержащие три переменные). Получим

с2 дР дР'

2л- дР иг, -иГР = шщ -т + г- +1^7

дх1 7 дх1 дх1

а подставляя (15) в (17):

с2 дР дР' д

'' Ох; + 7 ()Х; ()Х; +7 дх1

■С2 АР + с2 4-

дЪ

дх,

■с$Р'АР, (17)

9 Л • дР ^п „. , ■игР = —Ь —

ик

дхь

С

2°АР + с2°Ш-Л'

Т'^-)+С2Р'АР, (18)

или (еще раз используя (14) и отбрасывая члены следующего порядка малости)

2 Л г, 2п • дР с2 дРдР'

г'Сц д

ш дх,

дх1 7 дх1 дхг д2Р дР дщ \ 2 д дх{дхк дхк дхк ) 0 дх{

дх,

Это уравнение, правая часть которого выписана с точностью до членов первого порядка по турбулентным флуктуациям и, Т', Р' и описывает распространение звуковых волн в турбулентной атмосфере.

Найдем выражение для рассеянного поля, ориентируясь на выкладки Татарского [5]. Положим Р = Р0+Р5 где Р0 — падающая волна, а Р3 — рассеянное поле. Для слабого рассеяния (в линейном приближении) в правой части (22) можно положить Р = Ро. Так как ДРо + к2Ро = 0, получаем

АР8 + к2Р8 = -

д2

ш дх,дхь

т

' дхь

д_

дх1

,дРо дх.

Р'дяь.

7 ох1 0X1

(23)

Решение уравнения (23) имеет вид

йъг'

^Щг—г11

д2

ш дх,дхь

т

1 д.п

_д_

дх{

Т

дх.

1 дР0 дР'

7 дх{ дх{

■ Р'АР0

(24)

Пусть Ро(г') = Аое1кг . Сделаем приближение малости рассеивающего объема:

еЩг-г' I е

1кг

-1кпг'

(25)

где п — единичный вектор вдоль г (в направлении рассеянной волны). Преобразуем интеграл (24) по теореме -СдР'ДР. Гаусса, пренебрегая поверхностными интегралами, т.е.

считая поверхностные эффекты малыми по сравнению

(19)

с объемными. Тогда

Используем цепочку равенств

дР Щ-^— 4 ох1

ш дх.

ш дх.

1С0 д

д2р

дх,дхк

_д_

дхк

ик

1Сп

д2

Ш

дх,дхк

ик

дР

дх.

дР_

дXI

дР

■«¡т-

дхк

дР дщ

дхк дхк

дР дщ

дхк дхк

Р(г) =

г'Л0£'

{кг

4-лт

¿3 Г,е-Шш>

2 д2 ш дх,дх -,

щк/е1^^

ЩР)

дх.

, ч 1 яр/

(ТЧцеР" ) + -к^ + 1Р'кУкг V /7 дх.

А0е

1кг

ш дх.

иг

щ АР Н—=-Р

47ГГ

ё3г' еЖг>

кп,к,Т' ■

* 2 кп/к/п^

и,

(20)

7

-к(кщ

со

к + фдР'

(26)

отметим, что ДР + %-Р линейно по малым величинам

со

и, Т',Р', поэтому, отбрасывая выражения второго порядка малости, в итоге получаем

2 л г, 2п С2 дРдР'

Со АР + ш Р = —^ д—"о— 7 ох1 0X1

2с2 д2 ш дх,дхк

Щ

,21)

дР

дхк

или (используя дисперсионное соотношение ш2 = $с\)

АР + к2Р =

д2

ш дх,дхк

дР\ _ _д_

'дхк) дх1 _ \_дР_дР

7 дх{ дх1

Т

дР_

дх11

Р'АР.

(22)

где К = к-кп. Введем также дополнительный вектор К ц = кш+ (-)-№.

Как известно, средний вектор плотности потока энергии рассеянного излучения [5] :

5=^1гп[Р*УР} = \рос1п{Р*(г)Р(г)).

(27)

Предполагая, что мелкомасштабная турбулентность статистически изотропна, обозначим:

<7"(г,)Г(г2)> =Вт(р), (Г(г1)Р'(г2))=ВТР(р), <Г(г,Н(г2))= 0, (Р>(г1)Р>(г2))=ВР(р), (28)

(щ(п)щ(г2)) = с2Вц{р), {Р'{г1)щ{г2)) = 0,

где р = |г( — г21. Получаем

s_nPQC I у

еШр [й2(кп)2 + k2

V

Тогда S = п

(31)

Ро4 Ао

2 16тг2г2

V8tt3\k2®T{K){kn)2

+4£2Фг7(Ж){кп)2щп;+ФР{К){К'к)2^2ФТр(К){К'к)(Ал)].

(32)

Обозначим теперь через в угол рассеяния (чтобы обозначения совпадали с классической формулой Татарского). Тогда кп = !гсо$в и

щщФц{Щ =

cos

2 0

Е(К)

4 тгК2

(33)

a K'k = k2(cos6 + 1). Окончательно получаем

с 3^0

ж /г^ 2 л £(А") COS2 0 COS2 |

ФгСйОсоэ2^ —^-^

7Г

+ Ф/>(ЛГ)(соя0 + 7 - I)2 - 2®TP{K){cose + 7 - 1) cos в

(34)

Члены во второй строке дополняют формулу Татарского, но теперь описывает относительные флуктуации потенциальной, а не термодинамической температуры.

2. Турбулентные флуктуации давления и энтропии. Формулы Обухова

Вполне естественно, что уточненная формула Татарского совпадает с оригинальной в пренебрежении флук-туациями давления (Р' = 0), которые меньше температурных, если говорить об относительных величинах. Однако учет флуктуаций давления связывает формулу Татарского с формулой Обухова, описывающей связь флуктуаций давления с флуктуациями скорости [7]. Добавление же в уравнение связи новых членов качественно изменяет соотношение между флуктуациями скорости и флуктуациями термодинамических величин: они становятся относительно независимыми.

Чтобы показать необходимость такой независимости, повторим еще раз вывод A.M. Обуховым уравнения связи структурных функций термодинамических характеристик и скорости для модели локально изотропной

турбулентности, но с учетом сжимаемости воздуха. Запишем снова упрощенное уравнение Эйлера:

+ (К'к)2ВР - 2(К'к)(кп)ВТР\d3p. (29)

В традиционных обозначениях Фурье-образов корреляций

eiKpBT d5p = (2тг)3Фг(/0, ежрВц d5p = (2тг)3Фц{К).

(30)

Обозначим по аналогии (все флуктуации относительные)

eiKf>BP d^p = (2тг)3Ф/>(/0, eiKpBTp d^p = (2тг)3ФТР(К).

dV 1

■ (W)V=—Vp^V Ф;

(35)

Ы р

здесь Ф — потенциал внешних сил (геопотенциал). Влиянием силы Кориолиса на мелкомасштабные флуктуации пока, для простоты, пренебрегаем. Выразим силу градиента давления через две температуры:

1 ы

= \пв), (36)

или

дУ

~dt

(VV) V = —CpVT + TVS - \7Ф.

(37)

Даже если рассматривать только мелкомасштабные флуктуации и считать пока, что \7Ф = const = -CPVT0, то уже это уравнение показывает, что флуктуации температуры нельзя просто заменить флуктуациями потенциальной температуры (или энтропии), так как они связаны через флуктуации скоростей. Измерения турбулентности в пограничном слое [17] показывают, что интенсивности флуктуаций скорости и флуктуаций температуры относительно независимы, а это означает, что флуктуации термодинамической и потенциальной температуры имеют один масштаб, но не связаны функционально, иначе при отбрасывании одного из членов возникала бы одновременная корреляционная связь между интенсивностями флуктуаций двух других, подобная формуле Обухова [7].

В самом деле, если просто пренебречь флуктуациями энтропии (считать S = const), и полагая пока Vqi =0, для турбулентных флуктуаций скорости имеем:

^ + (и, \7)и = —CpVT. (38)

Здесь уже Т — абсолютные значения турбулентных флуктуаций термодинамической температуры. Если взять дивергенцию от правой и левой частей, то с учетом несжимаемости турбулентного движения divu = 0, получим

dvi dvk Ч °

Здесь по повторяющимся индексам подразумевается суммирование. Полагая, что поля турбулентных флуктуаций скоростей и температур статистически однородны и изотропны (в инерционном интервале масштабов), умножим оператор Лапласа от температуры в точке Mi на такое же выражение в точке М2 и усредним результат. Обозначим структурную функцию температуры Dj{p) = ((Т(М2) - Т(М\))2), так что A2DT(p) = -2А2(Т(М1)Т(М2)), где р — расстояние между точками наблюдений, которое, конечно, не следует путать с плотностью. Такое обозначение использовано здесь для сопоставления с формулой Обухова [7]. Пользуясь условием статистической однородности турбулентных полей, получаем

-СРАТ= .

дхь дх.

(39)

А2(Т(М{)Т(М2)) =

1 / dvi . „, . dvi . „, ч dvk

Дифференцирование в левой части проводится по координатам вектора, соединяющего точки М\ и М2. Тогда

д2£!г = ( тгт + --к-* ) От =

др4 р др3

2 дюк , 00; Д ,,,.

(41)

Полагая (аппроксимируя) случайное поле скоростей гауссовым, можно выразить четвертые моменты скорости через вторые [7], и

AzDT(p) = ^74

1 a2D'iS д2оа-у

С2 д£адС '

(42)

где

= ((va(M2) - {va{Ml)){vl3{M2) - (v^iMi))

= Dnn(p)SaP-

Dnn — Pit

(43)

Здесь — единичные орты, а Опп и Оц — продольная и поперечная структурные функции поля скоростей. Пользуясь формулами для поперечной и продольной структурной функции статистически однородного и изотропного поля [3], получаем связь структурных функций температуры и скорости:

От = ^Шр)}2,

(44)

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

Приближение S = const является довольно грубым. Независимость (функциональная) флуктуаций потенциальной и термодинамической температуры приводит к бароклинности турбулентных флуктуаций, что в свою очередь способствует генерации турбулентных вихрей. В самом деле, вектор бароклинности

В = rot

К

Р

Vp ) = Ср rot(^VF + TVS) = CpVT х VS,

(45)

поэтому флуктуации потенциальной температуры (или энтропии) нельзя сравнивать по амплитуде с флук-туациями термодинамической температуры (и отбрасывать малые компоненты), так как функциональная независимость этих двух температур является «источником» мелкомасштабной завихренности турбулентных течений в АПС (более крупномасштабная связана, очевидно, со средним течением и его вертикальным градиентом).

3. Энтропия и функция Ляпунова

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

пии атмосфере адиабатические флуктуации — звуковые волны — затухают не сразу. При этом (\75)2 принимает минимальное значение в статистически стационарном пределе, не равное нулю при генерации турбулентности из внешних источников. На это указывают и хорошо известные результаты анализа флуктуаций термодинамических характеристик в классической статистической физике [15]. Если считать, что выравнивание неоднородностей энтропии на малых масштабах происходит из-за турбулентного движения частиц жидкости, то можно предположить, что скорость изменения энтропии «частицы жидкости» единичной массы (в «физически бесконечно малом объеме») при отсутствии источников пропорциональна его лапласиану, поскольку это простейший оператор, который инвариантен к поворотам и сдвигам при условии локальной однородности и изотропности турбулентности. Коэффициент «диффузии» энтропии или сглаживания ее пространственных неоднородностей ниже обозначим для отличия от коэффициента температуропроводности как 7 (это обозначение использовано в работе [17], не следует путать с соотношением Ср/С0, введенным выше). Он может быть связан также и с кинематической вязкостью среды. Чтобы не путать пока термодинамическую энтропию с характеристикой турбулентности, выражающей экстремальные и инвариантные свойства турбулентного перемешивания, обозначим эту скалярную величину как ф. Смысл такого обозначения станет очевидным чуть ниже.

Уравнение для турбулентных микромасштабных флуктуаций скалярной характеристики ф тогда будет аналогично уравнению для флуктуаций температуры. В самом деле, если

(46)

то, учитывая несжимаемость жидкости и считая коэффициент «перемешивания» (диффузии) 7 = const, можно переписать уравнение (46) в виде

дф ~dt

дхк

фЩ~1

дф дхь

= 0.

(47)

Умножая теперь (47) на значение ф в той же точке и статистически усредняя ((•)), получим

,дф\ ФШ/

фд-^)^(фАф) = 0.

(48)

Используя свойство несжимаемости турбулентных флуктуаций скоростей, легко увидеть, что

дхк / Y~Kdxk Также очевидно, что д

Ф2щ

7 {фАф) = 7

,, дф. dxk

Поэтому, подставляя (49) и (50) в (48), получаем

(Ф\)

дф2 дхк

(49)

(50)

= 0.

(51)

Но, согласно теореме о некоррелированности однородного и изотропного случайного поля скоростей со скалярным полем, первый член в квадратных скобках равен нулю. Второй член также равен нулю в силу однородности флуктуации ф {(ф2) = const). Таким образом,

7<(W02). (52)

Поскольку уравнение (46) не включало источников поддержания флуктуаций ф, то и (52) описывает процесс затухания флуктуаций ф со временем. Эта формула может интерпретироваться и как определение ф: среди всех скалярных термодинамических характеристик турбулентной среды функция ф обладает свойством наиболее «быстрого» пространственного выравнивания, если ф имеет наибольший коэффициент «диффузии» 7. При наличии внешних источников турбулентного перемешивания, не обязательно тепловых, в стационарном пределе среднее значение (Чф)2 определяется интенсивностью этих источников, а величина ф2 принимает статистически стационарное значение в каждой пространственной точке.

Для случайного поля среднее значение ф2 может интерпретироваться и как плотность вероятности распределения частиц жидкости (одинаковой массы) в координатном пространстве, в том смысле, что для однородного случайного поля эта величина имеет постоянное значение, а для локально однородного (Чф)2 min по вероятности. Коэффициент «диффузии» или перемешивания 7 также характеризует время корреляции т турбулентных флуктуаций с пространственным масштабом А: т ~ . Общее распределение флуктуаций скалярной характеристики ф по спектру определяет характерное время корреляции турбулентных флуктуаций то, так что измерения через интервалы Д^»то могут рассматриваться как независимые. Функция ф является «почти детерминированной» функцией координат (т. е. функцией имеющей в основном низкочастотную и крупномасштабную динамическую изменчивость), а стационарность понимается как возможность рассматривать флуктуации ф в разные моменты времени реализациями стационарного случайного процесса и среднее по времени интерпретировать как среднее по ансамблю (эргодичность флуктуаций ф).

Заключение

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

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

Значимость добавленных членов в формуле Татарского быстро возрастает с ростом масштаба адиабатических колебаний, так как спектр турбулентных колебаний давления быстро падает с увеличением волнового числа [10]. Так, анизотропный градиент давления, связанный с геопотенциалом (см. формулу (37)), или термическая стратификация пограничного слоя и горизонтальные волны температуры/энтропии приводят к появлению новых добавочных членов в уравнении Эйлера для турбулентных пульсаций скорости и к завихренности бароклинного турбулентного движения [16].

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

Список литературы

1. Обухов A.M. Турбулентность и динамика атмосферы. М., 1981.

2. Монин A.C., Яглом A.M. Статистическая гидромеханика. Ч. 1. М„ 1965.

3. Монин A.C., Яглом A.M. Статистическая гидромеханика. Ч. 2. М„ 1967.

4. Ламли Дж.Л., Пановский Х.А. Структура атмосферной турбулентности. М., 1966.

5. Татарский В.И. Распространение волн в турбулентной атмосфере. М., 1967.

6. Осташев В.Е. Распространение звука в движущихся средах. М„ 1992.

7. Обухов A.M. И ДАН. 1949. 66, № 1. С. 17.

8. Татарский В.И. // ЖЭТФ. 1953. 25, № 1. С. 74.

9. Каллистратова М.А. // ДАН. 1959. 125, № 1. С. 69.

10. Голицын Г.С. II Изв. АН СССР. Сер. геофиз. 1964. № 8. С. 1253.

11. Kallistratoua М.А., Coulter R.L. 11 Meteorol. Atrnos. Physics. 2004. 85. P. 21.

12. Rusakou Y.S. 11 Meteorol. Z.. 2007. 16, N 4. P. 349.

13. Beyrich F., Weill A. 11 Bound. Layer Met. 1993. 63. P. 97.

14. Ландау Л.Д., Лифшиц Е.М. Гидродинамика М., 1988.

15. Климонтович Ю.Л. Статистическая физика. М., 1982.

16. Пономарев В.М., Чхетиани О.Г. 11 Изв. РАН. Физ. атм. океана. 2005. 41, № 5. С. 418.

17. Юшков В.П. II Вестн. Моск. ун-та. Физ. Астрон. 2011. № 3. С. 100.

The sound scattering on turbulent fluctuations of pressure and entropy E. V. Yushkov, V. P. Yushkov

Department of Physics of Atmosphere, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia. E-mail: yushkov@phys.msu.ru.

A well-known Tatarskij formula of the sound scattering in a turbulent atmosphere was analyzed. Appearance of additional terms in the classical formula is shown in the assumption that the acoustic fluctuations are adiabatic and the turbulent ones are incompressible. Comparison of the thermodynamic and potential temperature turbulent fluctuations shows that Obukhov formula, which describes the relationship of the pressure and velocity turbulent fluctuations, is changed in the atmosphere. Finally, a formula for entropy fluctuations is deduced by analogy with the formula for small-scale isotropic temperature fluctuations and the assumption that the thermodynamic entropy can be proportional to the informational one is discussed.

Keywords: turbulence, entropy, pressure, acoustic waves, compressibility, baroclinicity. PACS: 47.27.nb, 43.28.+h, 92.60.hk. Received 10 February 2011.

English version: Moscow University Physics Bulletin 6(2011).

Сведения об авторах

1. Юшков Егор Владиславович — аспирант; тел.: (495) 939-28-77, e-mail: yushkov@phys.msu.ru.

2. Юшков Владислав Пролетарьевич — канд. физ.-мат. наук, ст. науч. сотрудник; тел.: (495) 939-15-41, e-mail: yushkov@phys.msu.ru.