CC BY
10
НАУКИ О ЗЕМЛЕ / SCIENCE ABOUT THE EARTH
DOI: https://doi.org/10.23670/IRJ.2017.58.016 Булатов В.В.1, Владимиров Ю.В.2
1ORCID: 0000-0002-4390-4013, Доктор физико-математических наук, 2ORCID: 0000-0002-2139-5078, Кандидат физико-математических наук, Институт проблем механики им.А.Ю.Ишлинского РАН (Москва) Работа выполнена при частичной поддержке РФФИ (проект 17-01-00010) ВОЛНОВАЯ ДИНАМИКА ЭКСПОНЕНЦИАЛЬНО СТРАТИФИЦИРОВАННЫХ СРЕД
Аннотация
В статье рассмотрены вопросы, связанные с динамикой внутренних гравитационных волн в экспоненциально стратифицированных средах бесконечной глубины. Обсуждаются адекватность и пределы применимости используемой модели стратифицированной среды, а также основные качественные особенности динамики плоских волновых пакетов внутренних гравитационных волн. Показано что как минимум в двух отношениях, внутренние гравитационные волны в безграничной экспоненциально стратифицированной среде качественно отличаются акустических и электромагнитных колебаний: частоты колебаний распространяющихся внутренних волн не могут превосходить частоту плавучести, и внутренние гравитационные волны распространяются не по нормали к волновым фронтам.
Ключевые слова: внутренние гравитационные волны, динамика стратифицированной среды, частота плавучести, волновой фронт.
Bulatov V.V.1, Vladimirov Yu.V.2
1ORCID: 0000-0002-4390-4013, PhD in Physics and Mathematics, 2ORCID: 0000-0002-2139-5078, PhD in Physics and Mathematics, Institute for Problems in Mechanics RAS (Moscow) The work was partially supported by the RFBR (project 17-01-00010) WAVE DYNAMICS OF EXPONENTIALLY STRATIFIED MEDIUM
Abstract
The problems of internal gravity waves dynamics in unbounded exponentially stratified medium are considered. We discuss the relevance and limits of applicability of the used stratified medium model. The basic qualitative features of the internal gravity waves plane packets dynamics are obtained. It is shown that internal gravity waves in an unbounded exponentially stratified medium qualitatively differ from acoustic and electromagnetic oscillations. The frequency of propagating internal waves oscillations can not exceed the buoyancy frequency. Internal gravity waves propagate not along the normal to the wave fronts.
Keywords: internal gravity waves, stratified medium dynamics, buoyancy frequency, wave front.
В настоящей работе рассмотрены особенности волновой динамки внутренних гравитационных волн в безграничной экспоненциально стратифицированной среде. Вначале сделаем несколько замечаний об адекватности моделей стратифицированных сред, в которой распространяются колебания. Простейшей моделью природной стратифицированной среды является безграничная по вертикали экспоненциальная среда, т. е. среда, в которой невозмущенная плотность р (z) экспоненциально зависит от z : р (z) = exp(—qz). Для такой среды
линеаризованная система уравнений внутренних волн сводится к системе с постоянными коэффициентами, и её решение можно получить методом Фурье, хотя вычисление асимптотики волнового поля вдали от источника возмущений, равномерной по координатам точки наблюдения представляет собой достаточно сложную, но вполне решаемую аналитическую задачу [1], [2], [3], [4], [5], [6], [7]. Однако соответствие качественных свойств дальних полей в экспоненциальной безграничной по вертикали среды, полученных на основе асимптотических аналитических оценок, реальной волновой картине в океане требует, по всей видимости, дополнительного физического обоснования. Как правило, в океанических условиях распространение внутренних волн происходит в слое конечной толщины (толщины слоя скачка плотности), т.е. порядка 50-100 метров, и лишь в океане, где имеются два термоклина -сезонный и главный, толщина слоя имеет порядок одного километра и более. В то же время, как правило, представляет интерес распространение волн на дистанциях порядка десятков километров, т.е. на дистанциях, превышающих на один или два порядка толщину стратифицированного слоя. Ясно, что в этой ситуации заменять слой конечной толщины на безграничную по вертикали среду неправомерно. Дальнее поле в безграничной по вертикали среде и в слое конечной толщины имеют два качественно общих свойства: в обеих задачах поле при удалении от источника убывает, впрочем, характер убывания различен, и в обеих задачах имеются области, где поле имеет осциллирующий, волнообразный характер, однако эти пространственные области существенно различаются [1 ], [2], [3], [4]. При решении полной системы нелинейных уравнений гидродинамики, и в дальнейшем отщеплении линейной части решения в процессе численного расчета зон перемешивания, турбулентных следов, различных вихревых образований и возбуждаемых ими внутренних гравитационных волн модель экспоненциальной безграничной среды может оказаться физически вполне адекватной, т.к. основные гидрофизические процессы локализованы в относительно ограниченной пространственной области. Поэтому если найти из этих рассмотрений некоторое распределение источников возмущений, аппроксимирующее, например, турбулентный след или иные локализованные начальные неволновые возмущения, то далее можно рассчитывать по этим численным данным дальние волновые поля, уже учитывая конечную толщину стратифицированного слоя и реальное распределение невозмущенной плотности. Кроме того, модель бесконечной экспоненциально стратифицированной среды, по-видимому, является достаточно адекватной для описания полей внутренних гравитационных волн, возбуждаемых в искусственных
стратифицированных средах, в том числе в лабораторных условиях [5], [6], [7].
Рассмотрим далее плоские волновые пакеты внутренних гравитационные волн в экспоненциально стратифицированной среде. Самый простой вид уравнения внутренних гравитационных волн приобретают в случае экспоненциально стратифицированной жидкости, в которой невозмущенная плотность - экпоненциальная функция глубины: р (z) = exp(-qz). В такой среде N2 (z) = gq = const и основное уравнение, описывающее внутренние гравитационные волны (в приближении Буссинеска), становится уравнением с постоянными коэффициентами
LW = 0 (1)
Т 5 2 Л ХТ2Л Л 5 2 5 2 А Л 5 2
L = —rA,+ N2 А, А = —- + —-, А =А +—-
512 3 5x2 5y2 3 5z2
и система линеаризованная уравнений гидродинамик имеет вид
5U, 5p Л 5U2 5^ 5W 5p
—1 + —= 0, —1 + — = 0, -+ —+ pg = 0 (2)
5t 5x 5t 5y 5t 5z
5P-N2W = 0, 5Ui+5Ui+5W = 0
5t g 5x 5y 5z
Простейшие решения этих уравнений - это плоские волны, с зависимостью от t, Г = (x, y, z) в форме exp 1ф = exp i(-rat + kr), где k = (а, P, у) - волновой вектор. Повернем оси x, y так, чтобы Р обратилось в
5U 5p _ч ^
нуль, то есть чтобы рассматриваемое поле не зависело от y (тогда, в частности, -- =--— = 0). Полагая
5t 5y
U2 = 0, U = U10 exp(iO), W = W0 exp(iO), p = p0 exp(iO), p = p0 exp(iO), получим из (1) дисперсионное уравнение
ю2(а2 + у2) = N2 а2 (3)
а из (2) - систему уравнений относительно U10, , po, Р0
iN2
®U10 = а p0, «W0 = у p0 + igp0, ®P„ =--W0, aüj0 +у W0 = 0 (4)
g
Условие разрешимости системы (4) - выполнение дисперсионного соотношения (3). Положим а = kcOS 0 у = ksin0; |0 < —. Ограничиваясь для определенности волнами, распространяющимися в положительном
направлении оси x , получим ( U - постоянный амплитудный множитель)
^ хт Л
Ю = ,-= = NCOS 0 (5А)
л/а^г^2
U = Usin 0 expi[-©t + k(xcos 0 + zsin 0] W = -Ucos 0 expi[-©t + k(xcos 0 + zsin 0]
p = NU sin 0 expi[—®t + k(xcos 0 + zsin 0] (5В)
k
p = iNUexpi[-rat + k(xcos 0 + zsin 0] g
Прокомментируем выписанные соотношения. Частота колебаний ю не зависит от модуля k волнового вектора, определяется лишь его направлением и не превосходит частоты плавучести N. Максимальное значение Ю приобретает при 0 = 0 , в этом случае: p = U = 0, W = —U exp( —iNt + ikx ),
iNU , ч Л
р =--exp(—iNt + ikx), то есть волновой вектор направлен по горизонтали, а среда испытывает лишь
g
—
вертикальные колебания. При 0 —^ i — в пределе Ю = 0, и мы получаем стационарное движение, в котором:
W = 0 , U = Uexp(+ikz) , р = iNU exp(+ikz), p = ± NUexp(+ikz) , то есть имеется сдвиговое
g k
горизонтальное течение с гармонической зависимостью от z и возмущение плотности (с соответствующим ему возмущением давления) с той же зависимостью от z . Заметим, что для волнового уравнения и системы уравнений
Максвелла дисперсионное уравнение имеет вид Ю = c|k|, то есть ю, наоборот, не зависит от направления волнового
вектора и определяется лишь его длиной. Из последнего уравнения в (4) следует, что вектор скорости
(U10 exp(i9),U20 exp(i9)) перпендикулярен волновому вектору (a,Р), то есть что вектор скорости направлен
вдоль волновых фронтов - линий постоянной фазы Ф(^ x, z) = ®t — ax — yz = const. Поэтому, казалось бы,
внутренние гравитационные волны, подобно уравнениям Максвелла (где для плоской волны векторы E и H также направлены вдоль волновых фронтов), являются поперечными волнами. Однако, такое заключение преждевременно.
Согласно общей теории колебаний в средах с дисперсией, вектор групповой скорости, то есть скорости движения волновых пакетов с высокочастотным заполнением, определяется соотношением
дю(а, у) _ дю(а, у) да ' 2у ду
где ю(а, у) вычисляются из дисперсионного соотношения (5А)
дЮ(^ у) тт
U 2а = -Г- ' U 2у = ~ (6)
дю _ д да да
Na
Nsin2 9 дю _ Nsin0cos0
k ду k
Откуда видно, что направление вектора групповой скорости совпадает с направлением вектора (U1(), W0). То же
направление имеет среднее за период T =- значение вектора плотности потока энергии (pUj, pW). При
Ю
вычислении (pU^p^W) следует учесть, что физический смысл имеют не комплексные функции (5В), а их вещественные части
и = IU sin 9 cos(—at + T) (7)
W = — U| cos 9 cos(—rat + T), p = NU sin 9 cos(—rat + T)
k
где T = arg U + k(x cos 9 + z sin 9) и при вычислении потока энергии надо в качестве U, W, p брать
выражения (7). К аналогичному направлению распространения можно прийти, если искать решения этих уравнения лучевым методом (методом геометрической оптики), то есть строить для внутренних волн аналог ВКБ решения для волнового уравнения или уравнений Максвелла. Проведенный анализ показывает, что как минимум в двух отношениях, внутренние гравитационные волны в безграничной экспоненциально стратифицированной среде качественно отличаются акустических и электромагнитных колебаний. Во-первых, частоты колебаний распространяющихся внутренних волн не могут превосходить частоту плавучести N, и во-вторых, внутренние гравитационные волны распространяются не по нормали к волновым фронтам.
Список литературы / References
1. Bulatov V.V. Internal gravity waves: theory and applications / V.V.Bulatov, Yu.V.Vladimirov. - Moscow: Nauka Publishers, 2007. - 304 p.
2. Булатов В.В. Динамика негармонических волновых пакетов в стратифицированных средах / В.В.Булатов, Ю.В.Владимиров. - М.: Наука, 2010. - 470 с.
3. Bulatov V.V. Wave dynamics of stratified mediums / V.V.Bulatov, Yu.V.Vladimirov. - Moscow: Nauka Publishers, 2012. - 584 p.
4. Булатов В.В. Волны в стратифицированных средах / В.В.Булатов, Ю.В.Владимиров. - М.: Наука, 2015. - 735 с.
5. Sutherland B.R. Internal gravity waves / B.R.Sutherland. - Cambridge: Cambridge University Press, 2010. - 394 p.
6. Miropol'skii Yu.Z. Dynamics of internal gravity waves in the ocean /Yu.Z.Miropol'skii, O.V.Shishkina. - Boston: Kluwer Academic Publishers, 2001. - 421 p
7.Pedlosky J. Waves in the ocean and atmosphere: introduction to wave dynamics / J.Pedlosky. - Berlin-Heidelberg: Springer, 2010. - 260 p.
Список литературы на английском языке / References in English
1. Bulatov V.V. Internal gravity waves: theory and applications / V.V.Bulatov, Yu.V.Vladimirov. - Moscow: Nauka Publishers, 2007. - 304 p.
2. Bulatov V.V. Dinamika negarmonicheskih volnovih paketov v stratificirovannih sredah [Non-harmonic wave packets dynamics in stratified medium] / V.V.Bulatov, Yu.V.Vladimirov. - M.: Nauka, 2010. - 470 p. [in Russian]
3. Bulatov V.V. Wave dynamics of stratified mediums / V.V.Bulatov, Yu.V.Vladimirov. - Moscow: Nauka Publishers, 2012. - 584 p.
4. Bulatov V.V. Volni v stratificirovannih sredah [Waves in stratified medium] / V.V.Bulatov, Yu.V.Vladimirov. - M.: Nauka, 2015. - 735 p. [in Russian]
5. Sutherland B.R. Internal gravity waves / B.R.Sutherland. - Cambridge: Cambridge University Press, 2010. - 394 p.
6. Miropol'skii Yu.Z. Dynamics of internal gravity waves in the ocean /Yu.Z.Miropol'skii, O.V.Shishkina. - Boston: Kluwer Academic Publishers, 2001. - 421 p
7.Pedlosky J. Waves in the ocean and atmosphere: introduction to wave dynamics / J.Pedlosky. - Berlin-Heidelberg: Springer, 2010. - 260 p.
DOI: https://doi.org/10.23670/IRJ.2017.58.117 Гайдукова Е.В.1, Джалалванд А.2, Синкпеун Л.3
1ORCID: 0000-0002-3547-5538, кандидат технических наук, доцент, 2ORCID: 0000-0003-2311-5065, аспирант, 3ORCID: 0000-0001-9297-2904, аспирант, Российский государственный гидрометеорологический университет (РГГМУ), г. Санкт-Петербург АПРОБАЦИЯ МЕТОДИКИ ДОЛГОСРОЧНОГО ПРОГНОЗИРОВАНИЯ СТАТИСТИЧЕСКИХ ХАРАКТЕРИСТИК НА РЕЧНЫХ БАССЕЙНАХ ИРАНА
Аннотация
В статье рассматривается возможность применения методики оценки последствий изменения климата к речным бассейнам Ирана. В качестве исследуемых взяты бассейны рек Сефидруд и Карун. Апробация методики произведена на ретроспективном материале с выделением маловодных и многоводных фаз водности с целью имитации условий изменения климата. В основе прогнозной методики лежит уравнение Фоккера-Планка-Колмогорова. Получено, что для надежных оценок статистических характеристик речного стока рек Ирана необходимо использовать региональные зависимости для параметров прогнозного уравнения.
Ключевые слова: гидрология, прогнозирование, разностно-интегральные кривые, ретроспективные данные, изменение климата.
Gaidukova E.V.1, Jalalvand A.2, Sinkpehoun L.3
1ORCID: 0000-0002-3547-5538, PhD in Engineering, Associate professor,
2ORCID: 0000-0003-2311-5065, postgraduate student,
3ORCID: 0000-0001-9297-2904, postgraduate student, Russian State Hydrometeorological University (RSHU), Saint-Petersburg
APPROBATION OF METHODS OF LONG-TERM FORECASTING OF STATISTICAL CHARACTERISTICS
ON RIVER BASINS OF IRAN
Abstract
The article considers the possibility of applying the methodology for assessing consequences of the climate change to the river basins of Iran. The paper investigates Sefid-Rud and Karun river basins. Approbation of the technique was performed with the help of retrospective material with the identification of low-water and high-water content phases with the aim of simulating climate change conditions. The forecast method is based on the Fokker-Planck-Kolmogorov equation. According to the results, reliable estimates of statistical characteristics of the river flow in Iran, it is necessary to use regional dependencies for the parameters of the forecast equation.
Keywords: hydrology, forecasting, difference integral curves, retrospective data, climate change.
Введение
В настоящее время факт изменения климата подтверждается инструментальными данными [1]. Возникает задача оценки последствий изменения климата на различные отрасли экономики. Существует несколько подходов в этом направлении (см. например, [2], [3], [4], [5]), в том числе методика, разработанная в Российском государственном гидрометеорологическом университете (РГГМУ) [6].
Целью данного исследования является апробация методики РГГМУ на речных бассейнах Ирана на ретроспективном материале, т. е. составление поверочных прогнозов статистических характеристик речного стока.
Исходные данные
В качестве исследуемых речных бассейнов были взяты водосборы рек Сефидруд и Карун, которые располагаются на севере и юго-западе Ирана. Расположение гидрологических постов в бассейнах Сефидруд и Карун представлено на рис. 1. Были отобраны посты с незарегулированным стоком и замыкающие зональные площади водосборов. Всего получилось 28 постов из 56 имеющихся. Продолжительность рядов составляет в среднем 39 лет, большинство рядов начинаются с 1968 года и заканчиваются 2006 годом.
