Run-up of nonlinear monochromatic wave on a plane beach in presence of a tide

Cover Page

Abstract


The nonlinear problem of long wave run-up on a plane beach in a presence of a tide is solved within the shallow water theory using the Carrier-Greenspan approach. The exact solution of the nonlinear problem for wave run-up height is found as a function of the incident wave amplitude. Influence of tide on characteristics of wave run-up on a beach is studied.


ВВЕДЕНИЕ

Накат волн на побережье в большинстве районов Мирового океана происходит в присутствии прилива, высота которого может быть значительной. В случае цунами подход волны во время прилива (высокой воды) усугубляет катастрофу, и цунами проникает на большие расстояния в глубь побережья. Напротив, во время отлива (малая вода) даже сильное цунами может оказаться незаметным, и такие случаи известны из наблюдений. Именно поэтому в международной практике принято измерять характеристики цунами «без прилива», вычитая последний из мареограмм или просто уменьшая соответственно высоту волны в наблюдениях. Сам прилив в настоящее время хорошо предвычисляется, так что ошибки, связанные с приливными колебаниями, могут быть минимизированы. При расчете же общего уровня подъема воды колебания прилива складываются с колебаниями цунами по принципу линейной суперпозиции [3].

Между тем воздействие цунами на побережье с учетом прилива представляет собой сложную нелинейную задачу, и простой суперпозиции волн цунами и приливных волн недостаточно. Так, значительная роль прилива в процессе генерации волн цунами оползнями отмечается в работе [17]. Во время большой воды береговой склон разжижается, а во время отлива может потерять устойчивость и сползти в воду, генерируя при этом волну цунами. Взаимодействие волн с приливом влияет на экологию прибрежной зоны, в частности на смешение пресных грунтовых вод с соленой морской водой [11]. В численных расчетах трансформации волн цунами с приливом в одном из населенных пунктов Испании демонстрируется, что максимальная глубина потока воды на берегу очень сильно зависит от приливного уровня, в то время как максимальная высота наката менее подвержена такому влиянию [8]. Приливные колебания могут воздействовать на цунами и на подходе их к бухтам, усиливая или ослабляя в зависимости от фазы прилива [4]. Проникновение волн цунами в реки также зависит от уровня прилива, и численные эксперименты подчеркивают важную роль взаимодействия цунами с приливом и речным стоком [6, 9, 10]. Эти выводы, хорошо понятные с точки зрения геофизики, получены численно. Насколько нам известно, аналитические решения, описывающие влияние прилива на характеристики наката волн цунами на берег, отсутствуют, за исключением работы [7], в которой рассматривается накат волны на вертикальную стенку. Целью настоящей статьи является исследование влияния прилива на высоту наката длинной необрушенной волны на плоский откос в рамках одномерной нелинейной теории мелкой воды.

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ

Рассмотрим простейшую задачу наката длинных волн на плоский откос в присутствии приливных колебаний уровня воды. Анализ мы проведем в рамках нелинейной теории мелкой воды, уравнения которой запишем в классическом виде:

ηt+xh+ηu=0, (1)

ut+uux+gηx=0, (2)

где η — смещение уровня моря, u — усредненная по глубине скорость водного потока, g — ускорение свободного падения и h(x)=αx — переменная глубина бассейна, α — тангенс угла наклона берега к горизонту (рис. 1). Ось x направлена к берегу, и невозмущенный урез находится в точке х = 0. Граничным условием на берегу является равенство нулю полной глубины:

H(x,t)=h(x)+η(x,t)=0, (3)

которое определяет положение подвижного уреза x(t) со временем.

 

Рис. 1. Геометрия задачи.

 

Решение нелинейных уравнений (1)–(2) обычно находится с помощью преобразования Карриера-Гринспана [2, 5], сводящего систему (1)–(2) к линейному волновому уравнению для вспомогательной функции Ф(σ, λ) в новых переменных σ и λ:

2Φλ22Φσ21σΦσ=0, (4)

и все физические переменные находятся явно через данную функцию Ф(σ, λ):

η=12gΦλu2, u=1σΦσ, (5)

x=12αgΦλu2σ22, t=1αgλu. (6)

Такое аналитическое представление задачи связано с тем, что исходная система (1)–(2) в терминах полной глубины является гиперболической системой с постоянными коэффициентами, что позволяет использовать классическое преобразование годографа (Лежандра) (см. [2, 5]). Возникающая при этом волновая функция Ф не имеет физического смысла, но через нее полностью определяются физические переменные. Обратим внимание, что переменная λ пропорциональна времени, а σ пропорциональна полной глубине бассейна:

σ=2gH(x,t), (7)

так что точка σ = 0 соответствует подвижному урезу. В результате, волновое уравнение (4) решается на полуоси σ ≥ 0, и граничное условие на подвижном урезе (3) трансформируется в

Φ(σ=0,λ)<, Φσ(λ,σ=0)=0, (8)

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

Преобразование Карриера-Гринспана неоднократно использовалось в литературе для анализа наката длинных волн на плоский откос и получения разнообразных характеристик наката. Одним из свойств полученного Карриером и Гринспаном решения для плоского откоса является возможность использования двухшагового подхода для изучения свойств движущегося уреза [2]: на первом этапе решается линейная задача и определяется вертикальное смещение уровня воды на неподвижном урезе (х = 0): R(t) и скорость частиц в этой точке U(t), а затем находятся «истинные» нелинейные вертикальные смещения подвижного уреза (σ = 0): r(t) и скорость движения подвижного уреза u(t) по формулам:

r(t)=Rt+ugαu22g, u(t)=Ut+ugα. (9)

Более того, удается записать критерий обрушения волны, под которым понимается градиентная катастрофа — обращения в бесконечность du/dt:

Br=max(d2R/dt2)gα2=1. (10)

РЕЗУЛЬТАТЫ ЧИСЛЕННЫХ РАСЧЕТОВ

Мы выполнили серию численных расчетов динамики волн на фоне прилива по формулам (4)–(6), где в качестве функции Ф(σ, λ) взята сумма приливной и «цунамной» компонент функции:

Φ=A0J0ΩσgαsinΩλgα+AJ0ωσgαsinωλgα+φ, (11)

и φ представляет собой сдвиг фаз между ними. Здесь A0 и А — амплитудные характеристики прилива и цунами (высоты прилива и цунами рассчитываются по формуле (5)), Ω — частота прилива, ω — частота цунами, а J0 — функция Бесселя первого рода. В качестве примера взята волна цунами с периодом 10 минут и амплитудой наката 8.9 м на фоне прилива с периодом 12 часов и амплитудой 3 м. Тангенс угла берегового склона α равен 0.01. Параметр обрушения для этого случая, рассчитанный по формуле (10), равен Br = 0.99, что характеризует сильно нелинейные волны перед самым обрушением. Сдвиг фаз φ = 0. Из-за большой разницы в периодах прилив фактически представляет собой «статическое» смещение уровня воды в соответствующей фазе с очень малой скоростью течения, однако здесь мы не будем использовать это естественное физическое предположение.

Колебания уровня воды и скорость движения уреза показаны на рис. 2–4. Можно сказать, что колебания уровня воды на урезе (рис. 2) представляют собой суперпозицию прилива и волн цунами. Это подтверждает и рис. 3, на котором проведено сравнение динамики цунами в условиях максимального и минимального приливов. Скорость же движения уреза полностью определяется волнами цунами. Нелинейные эффекты, связанные с фазой прилива, проявляются в незначительном изменении скорости волн на берегу и достаточно малы (рис. 4). Так, из рис. 4 видно, что при отливе слегка повышаются скорости наката волн цунами, в то время как при максимальном приливе слегка увеличиваются скорости отката.

 

Рис. 2. Колебания волн с периодом 10 минут и амплитудой наката 8.9 м на берегу на фоне прилива с периодом 12 часов и амплитудой 3 м. Тангенс угла берегового склона α = 0.01. Сдвиг фаз φ = 0.

 

Рис. 3. Сравнение колебаний наката волн при максимальном приливе (сплошная линия) и отливе (штриховая линия) для условий рис. 2.

 

Рис. 4. Сравнение скоростей наката волн при максимальном приливе (сплошная линия) и отливе (штриховая линия) для условий рис. 2.

 

ЗАКЛЮЧЕНИЕ

Таким образом, нами рассмотрен накат длинной волны на плоский откос в присутствии прилива. Исследование проводилось в рамках нелинейной теории мелкой воды с помощью преобразования Карриера-Гринспана. Расчеты подтверждают, что за время добегания волны до берега прилив можно считать однородным в пространстве. Даже в случае сильно нелинейного режима, когда длинные волны близки к обрушению, колебания уровня воды на берегу представляют собой почти линейную суперпозицию прилива и волн цунами. Скорость же движения уреза полностью определяется волнами цунами. Нелинейные эффекты, связанные с фазой прилива, проявляются в малом изменении скорости волн на берегу; при отливе скорости наката длинных волн цунами повышаются, а при приливе растет скорость отката волны от берега. Эти выводы обусловлены сильной идеализацией задачи (необрушенные волны, плоский откос, простирающийся на берег, одномерное распространение). При этом очевидно, что береговой откос не является постоянным и может перейти в любой другой, например в почти пологий пляж. В этом случае волна цунами во время высокой воды может очень далеко проникнуть на побережье в отличие от низкой воды. Однако для такого откоса, как и для пространственной задачи, не удается получить аналитического решения, которое могло бы использоваться в качестве тестового в численных расчетах.

Источник финансирования. Работа выполнена при финансовой поддержке гранта PUT1378, грантов РФФИ (№ 17-05-00067 и № 18-05-80019), а также гранта Президента Российской Федерации для государственной поддержки ведущих научных школ Российской Федерации НШ-2685.2018.5.

I. I. Didenkulova

Tallinn University of Technology; Nizhny Novgorod State Technical University n.a. R.E. Alekseev

Author for correspondence.
Email: dii@appl.sci-nnov.ru

Estonia, Tallinn; Nizhny Novgorod, Russia

Department of Marine Systems

E. N. Pelinovsky

Nizhny Novgorod State Technical University n.a. R.E. Alekseev; Federal Research Center "Institute of Applied Physics of the Russian Academy of Sciences"; National Research University — Higher School of Economics; Special Research Bureau for Automation of Marine Researches

Email: pelinovsky@appl.sci-nnov.ru

Russian Federation, Nizhny Novgorod; Nizhny Novgorod; Moscow; Yuzhno-Sakhalinsk

  1. Куликов Е. А., Рабинович А. Б., Файн И. В. и др. Генерация цунами оползнями на тихоокеанском побережье Северной Америки и роль приливов в этом процессе // Океанология. 1998. Т. 38. C. 361-367.
  2. Пелиновский Е. Н. Гидродинамика волн цунами. Нижний Новгород: ИПФ РАН, 1996. 276 с.
  3. Шевченко Г. В. Влияние прилива на формирование опасных подъемов уровня на побережье Охотского моря при цунами и штормовых нагонах // Океанология. 2017. Т. 57. № 5. C. 690-701.
  4. Behera M. R., Murali K., Sundar V. Effect of the tidal currents at the amphidromes on the characteristics of an N-wave-type tsunami // Proceedings of the Institution of Mechanical Engineers Part M: Journal of Engineering for the Maritime Environment. 2011. V. 225. P. 43-59.
  5. Carrier G. F., Greenspan H. P. Water waves of finite amplitude on a sloping beach // J. Fluid Mech. 1958. V. 4. P. 97-109.
  6. Chanson H. Tidal Bores, Aegir, Eagre, Mascaret, Pororoca: Theory and Observations. Singapore: World Scientific. 2011. 220 p.
  7. Didenkulova I., Pelinovsky E. Tsunami wave run-up on a vertical wall in tidal environment // Pure Appl. Geophys. 2018. V. 175. P. 1387-1391.
  8. Lima V. V., Miranda J. M., Baptista M. A. et al. Impact of a 1755-like tsunami in Huelva, Spain // Natural Hazards and Earth System Sciences. 2010. V. 10. P. 139-148.
  9. Tolkova E., Tanaka H., Roh M. Tsunami observations in rivers from a perspective of tsunami interaction with tide and riverine flow // Pure and Applied Geophysics. 2015. V. 172. P. 953-968.
  10. Tolkova E. Tsunami Propagation in Tidal Rivers. Springer, 2018. 131 p.
  11. Xin P., Robinson C., Li L., Barry D., Bakhtyar R. Effects of wave forcing on a subterranean estuary // Water Resources Research. 2010. V. 46. P. W12505.

Supplementary files

Supplementary Files Action
1. Fig. 1. The geometry of the problem. View (826KB) Indexing metadata
2. Fig. 2. Oscillations of waves with a period of 10 minutes and a coastal amplitude of 8.9 m on the shore against a tide with a period of 12 hours and an amplitude of 3 m. The tangent of the angle of the coastal slope is α = 0.01. The phase shift is φ = 0. View (875KB) Indexing metadata
3. Fig. 3. Comparison of fluctuations in the wave run at maximum tide (solid line) and low tide (dashed line) for the conditions of Fig. 2. View (939KB) Indexing metadata
4. Fig. 4. Comparison of the speeds of the wave run-up at the maximum tide (solid line) and low tide (dashed line) for the conditions of Fig. 2. View (940KB) Indexing metadata

Views

Abstract - 16

PDF (Russian) - 28

Cited-By


PlumX

Refbacks

  • There are currently no refbacks.

Copyright (c) 2019 Russian academy of sciences

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies