CC BY
73
Вычислительные технологии
Том 15, № 1, 2010
Численное моделирование поверхностных волн, возникающих при движении подводного оползня
по неровному дну*
Г. С. Хакимзянов Институт вычислительных технологий СО РАН, Новосибирск, Россия
e-mail: khak@ict.nsc.ru
Н. Ю. Шокина
Институт вычислительных технологий СО РАН, Новосибирск, Россия Section of Applied, Mathematics, University of Freiburg, Freiburg, Germany e-mail: shokina@mathematik.uni-freiburg.de
Выполнены расчеты волновых режимов при движении подводного оползня по криволинейному склону. Получены уравнения движения подводного оползня под действием сил тяжести, плавучести, трения и сопротивления воды, и исследована зависимость возникающей волновой картины от начального заглубления оползня, его размеров и кривизны склона. Моделирование поверхностных волн, генерируемых при движении оползня по неровному дну, проводилось в рамках нелинейной модели мелкой воды с использованием схемы предиктор-корректор на адаптивной сетке.
Ключевые слова: подводный оползень, неровное дно, уравнения мелкой воды, поверхностные волны, силы тяжести, плавучести, трения и сопротивления воды, численное моделирование, конечно-разностная схема, адаптивная сетка.
Введение
В естественных и искусственных водоемах подводный оползень может привести к образованию поверхностных волн больших амплитуд, опасных для береговых сооружений. В последние годы появилось много публикаций (см. библиографический список работы [1]), посвященных результатам вычислительных и лабораторных экспериментов, воспроизводящих механизм генерации поверхностных волн движением по плоскому подводному откосу недеформируемого твердого тела. В статье [2] показано, что замена реального оползня модельным недеформируемым твердым телом и выбор для него подходящего закона движения приводят к адекватному описанию волновых процессов в широком диапазоне изменения таких параметров как угол наклона плоского откоса, толщина и длина оползня, его начальное заглубление.
В работах [2-4] показано, что численное моделирование возникающих волновых процессов на основе "полных" математических моделей, учитывающих вертикальные перемещения жидкости, позволяет воспроизводить значения амплитуд генерируемых
*Работа выполнена при финансовой поддержке РФФИ (грант № 09-05-00294).
© ИВТ СО РАН, 2010.
волн, близкие к экспериментальным, В указанных работах исследовалась возможность использования различных приближенных математических моделей гидродинамики для описания оползневого механизма генерации поверхностных волн. Эти исследования показали, что общая картина возникающих волновых режимов приемлемо описывается даже с помощью модели мелкой воды, в которой вертикальные движения жидкости не принимаются во внимание.
Заметим, что в упомянутых работах рассматривалось движение только недеформи-руемого оползня и только по плоскому откосу, В настоящих исследованиях выполнено моделирование генерации поверхностных волн оползнем, движущимся по неровному склону. Предложено новое уравнение движения оползня, при выводе которого учитываются неровность склона, силы тяжести, плавучести, трения и сопротивления воды. Для плоского откоса полученное уравнение при пренебрежении силой трения совпадает с приведенным в [2].
Моделирование поверхностных волн, генерируемых при движении оползня по неровному дну, авторами осуществлялось в рамках нелинейной модели мелкой воды с учетом подвижности дна. Для расчетов использовалась схема предиктор-корректор [5] на адаптивной сетке, сохраняющая в линейном случае монотонность профилей численного решения. Схема опробована на тестовой задаче с известным аналитическим решением, описывающим генерацию поверхностных волн недеформируемым телом, движущимся с постоянной скоростью по горизонтальному дну. Для неровного дна выполнен анализ зависимостей характеристик волновых режимов от крутизны склона, начальной глубины залегания оползня, его длины и толщины,
1. Постановка задачи
Рассмотрим плоский слой жидкости, ограниченный сверху свободной поверхностью, снизу — непроницаемым подвижным дном. Предполагается, что жидкость находится в поле силы тяжести, является несжимаемой и невязкой. Пусть декартова система координат Охг выбрана так, что уравнение свободной поверхности покоящейся жидкости имеет вид г = 0 и все характеристики жидкости зависят только от переменных х, Для моделирования гидродинамических параметров изучаемых явлений будем использовать нелинейную модель теории мелкой воды, уравнения которой имеют вид
£ + §[ = С, *е(0,1), (1)
где £ — время, Ь — длина области, и — вектор решения, ¥ — вектор потоков,
и = ( Ни)> ¥ (и) = ( Ни2 4^/2 ) ' ° = ( дНЬ* ) '
и(х, £) — скороеть, Н = п + Ь — полная глубина, п(х, £) — отклонение свободной поверхности от невозмущенного уровня г = 0, уравнение
г = -к(х, £) = Ны(х) + (х, £)
задает форму подвижного дна, Н\л (х) и Н&\(х,£) — функции, определяющие его неподвижную и подвижную части, д — ускорение свободного падения.
На левой границе области предполагается наличие вертикальной непроницаемой стенки, правая граница считается открытой. На этих границах ставились условия непротекания и свободного прохода соответственно.
Предполагается, что в начальный момент времени жидкость покоится и известна функция г = к°1(ж), описывающая начальную форму оползня:
На (ж, 0) = (ж).
Форма и положение оползня при £ > 0 определяются законом его движения (выведен в следующем разделе).
Результаты [3] моделирования поверхностных волн, генерируемых оползнем, движущимся по плоскому откосу, дают основание утверждать, что приведенная математическая модель позволяет получить общее представление об исследуемом волновом процессе, В настоящей работе эта модель применяется для исследования поверхностных волн, возникающих при движении оползня по неровному склону, заданному однозначной функцией
г = кы (ж). (2)
2. Закон движения подводного оползня
Одним из главных факторов, определяющих особенности процесса волнообразования, обусловленного движением подводного оползня, является закон его движения. При моделировании оползня твердым телом такой закон задается для одной из характерных точек тела, например, для его центра масс, В работе [2] предложен закон движения оползней по плоскому откосу, В настоящем разделе дается вывод закона движения оползня по неровному подводному склону.
Предполагается, что в начальный момент времени оползень имеет конечные протяженность Ь и толщину Т, а его поверхность описывается функцией г = (ж) + к0(ж), где Н0 (ж) — заданная неотрицательная функция с конечным носителем (ж°, ж°) длины
Ь (вдоль оси Ож) и максимальным значением Т > 0 в некоторой точке промежутка (ж° ж0)
При £ > 0 оползень может начать движение по наклонному дну. Опишем применяемый в работе упрощенный способ моделирования оползневого процесса, основанный на следующих предположениях:
1 — при выводе закона движения оползня последний отождествляется в каждый момент времени с некоторой материальной точкой хс(£), имеющей абсциссу жс(£) и скользящей вдоль неровного основания (2) согласно закону движения материальной точки по плоской кривой, при этом жс(0) = ж0, где ж° € (ж0, ж°);
2 — положение оползня определяется исключительно значением абсциссы жс(£) выделенной точки хс(£), а его поверхность при £ > 0 описывается функцией
г = кы (ж) + Н81 (ж,£), (3)
где Н81 (ж,£) = (ж + ж° — жс(£)).
Таким образом, в момент времени £ оползень располагается на склоне между точками с абсциссами жД£) = ж° — ж° + жс(£) и жт (£) = ж° — ж° + жс(£) = ж\ (£) + Ь и при движении его длина (вдоль горизонтальной оси Ож) остается неизменной. Отметим также, что при движении оползня его поверхность деформируется в соответствии с встречающимися
неровностями дна, однако объем оползня V и его масса M = psiV остаются неизменными (здесь psi — плотность материала оползня).
Опишем далее закон движения оползня, точнее, укажем способ вычисления функции x = xc(t), поскольку именно эта функция участвует в определении поверхности (3) движущегося оползня.
Пусть s — длина дуги кривой (2), отсчитываемая от ее начала. По предположению
x
xs
x
s(x) = I у/1 + [/4(£) (4)
0
Предположим, что выделенная в начальный момент времени точка xc(0) соответствует параметру s = S0, а движущаяся точка xc(t) с абсциссой xc(t) — параметру s = S(t). Тогда закон движения s = S(t) первоначадьно покоящейся материальной точки xc(t) вдоль кривой (2) имеет следующий вид:
шЗД = FT (t), S (0) = So, ¿(0) = 0, (5)
где S = d2S/dt2, S = dS/dt, m — имеющая размерность массы, FT(t) —
t
точку xc(t),
mM массу воды Cwpw V, т. e. m = M + Cwpw V = (psi + Cwpw) V, где Cw — коэффициент присоединенной массы, pw — плотность воды, принимаемая в расчетах равной единице,
FT
материальную точку xc(t), В вертикальном направлении на каждый элемент оползня dx
В произвольной точке x G (xi(t), xr(t)) касательная составляющая этих двух противоположно направленных сил к кривой (2) вычисляется по формуле
fg(x, t) = (psi - pw) hi(x,t)Wg sin 0(x) dx,
где W — ширина оползня в направдении, перпендикулярном плоскости xOz, 0(x) — локальный угол наклона дна, при этом
h (x)
9(х) = — arctan 1г'ы(х), sin 9(х) =-- bt
y/i + [K tW]2
Таким образом, суммарное действие сил тяжести и плавучести на оползень конечного размера определяется силой
Xr(t)
Fg (t) = (psi - pw ) Wg y hi (x, t) sin 0(x) dx, (6)
Xl(t)
действующей в касательном направлении па материальную точку xc(t), движущуюся вдоль кривой (2), Сила (6) заставляет оползень ускоряться, в то время как сопротивление воды и трение оползня о дно способствуют замедлению движения. Касательная
составляющая силы сопротивления воды направлена против движения оползня, пропорциональна наибольшей площади TW поперечного сечения оползня и равна
Fr(t) = ~^CdPwTW(S)2, (7)
где S? — скорость движения точки x0(t) вдоль кривой (2), Cd — коэффициент сопротивления воды. Отметим, что сила сопротивления (7) является пассивной, т, е, при отсутствии движения она исчезает.
Движение оползня тормозит также сила трения, направленная по касательной к кривой (2) в сторону, обратную направлению движения оползня. Эта сила тоже пассивна, поскольку при отсутствии движения она не может изменить местоположение оползня и ее влияние на положение оползня будет проявляться только при его движении (S? = 0), Вначале рассмотрим элемент оползня длины dx, для которого силу трения будем определять по нормальной реакции N(x,t), действующей на выделенный элемент со стороны дна:
f (x,t) = -Cf N (x,t), (8)
где Cf — динамический коэффициент трения (коэффициент трения скольжения при
движении), В рассматриваемой упрощенной модели оползневого процесса будем счи-
Cf
ента трения скольжения при покое, который будем определять путем указания угла трения 6*, т.е. такого порогового значения угла плоского откоса, при превышении которого оползень начинает соскальзывать по плоскому откосу. Этот же угол трения будет использоваться для криволинейных откосов. Из приведенных предположений вытекает следующая формула для определения динамического коэффициента трения: Cf = tan 6*.
Нормальную силу N(x, t), действующую на элемент dx оползня, определим как сумму нормальной составляющей силы тяжести и центробежной силы:
N(x,t) = (psi - pw) К\(x,t)Wg cos 6(x) dx + (psi - pw) hsi(x,t)WK(x) (s?(t)) dx, (9) где K(x) — знакоопределенная кривизна основания дна:
K(x) = hl(x)f^l + [h>bt(x)f
Для плоского откоса K (x) = 0, в силу чего выражение (9) не содержит второго слагаемого, Из формул (8), (9) видно также, что в сравнении с плоским откосом в рассматриваемом случае неровность дна способствует уменьшению силы трения (по абсолютной величине) на выпуклых вверх участках дна (h'¿t (x) < 0) и, наоборот, возрастанию — на выпуклых вниз участках дна (h^ (x) > 0),
С учетом (9) получаем, что на движущийся оползень будет действовать следующая сила трения (приложенная в точке xc(t) и касающаяся в этой точке кривой (2)):
Ffr(t) = -Cf (psi - Pw) W J hi(x,t)
Xl(t)
Xr (t)
g cos 6(x) + K (x) (S(t)~
dx. (10)
Равнодействующая сил (6), (7), (10) и принимается за величину фигурирующую в формуле (5), Таким образом, с учетом принятых допущений задача (5) запишется в следующем виде:
(рsi + Cwpw) VS = (psi — pw) Wgli - W
(Psi ~ Pw) h + 7^CdpwT
(S?)2, (11)
S (0) = So, ¿(0) = 0, (12)
где
xr (t) Xr(t)
m=J ь (x,t)H - C/ cos = c/j Ь (x,t)K(x)dx,
x¡(t) Xl(t)
а искомая величина xc(t) связана с решением S(t) задачи (11), (12) соотношением (4):
Xc{t)
S(t)= I yJi + iKM)]2 С
o
Заметим, что уравнение (11) справедливо только для движущегося оползня. Поскольку оползень в начальный момент времени покоился, то этим уравнением можно воспользоваться при t > 0 только в том случае, если оползень начнет движение. Из (11) следует, что достаточным условием возможности сдвига оползня в начальный момент времени является положительность интеграла Д при t = 0:
I1 = J h°(x) [sin Í(x) — Cf cos Í(x)] dx > 0.
Xo xi
В силу нелинейности уравнения (11) для поиска его решения необходимо использовать численные методы. Для этого целесообразно переписать задачу (11), (12) в виде задачи Коши для системы обыкновенных дифференциальных уравнений первого порядка:
1 vi г>с(0) = 0, (13)
(7 + C»)^c=(7-l)^i-
(7"1)/2 + уТ
¿ = Б (0) = Бс, (14)
где 7 = р^/рад — отношение плотностей материала оползня и воды, 7 > 1, Решение данной задачи ищется до момента остановки оползня, т, е, до такого момента времени, при котором впервые значение ус скорости оползня обратится в пуль (с заданной точностью),
3. Модельная акватория и модельный оползень
Для исследования влияния неровности рельефа дна на поверхностные волны, порождаемые движением подводного оползня, был выбран модельный рельеф (рис, 1) с ненулевой кривизной, заданный с помощью гладкой монотонно убывающей функции (2) вида
7 / \ + Л— — Л— г . ., . .
г = пы\х) =----1----1ап11 [с(х - , (15)
z, м Q} град
Рис. 1. Модельный рельеф z = hbt(x) (i) и локальные углы наклона рельефа 9(x) (2)
где h+ h- — глубина дна соответственно в правой и левой бесконечно удаленной точке;
2 tan во 1 ho - h+
с = ---—: во — максимальным vro.n наклона рельефа: £ = —ш --— — точка
h- — h+ 2c h- — h0
перегиба; h0 — глубина в точке x = 0. В расчетах использовались следующие значения параметров: h+ = —1400 м, h- = — 5 м, в0 = 6o, h0 = —10 м. Длина области составляла 70 000 м.
Начаньпая форма модельного оползня задавалась по формуле
Т
h^(x) = { 2
, 2n(x — x0)
1 + eos ' v
|x Х0| — b/2, ^g-j
ь
0, |х - х0| > Ь/2
где х0 Т, Ь — заданные абсцисса вершины оползня при £ = 0, толщина оползня и его длина (вдоль оси Ох) соответственно, при этом х0 = х0 — Ь/2 х0 = х0 + Ь/2 = х0 + Ь. В расчетах толщина оползня Т была выбрана равной 25 м, длина Ь — 5000 м. В начальный момент времени вершина оползня находилась над точкой дна, которой соответствует глубина 4 = —200 м, поэтому
х° = £ - - ln (~ 12644 м. 2c \ h- — dcJ
Указанные значения наряду со значениями констант, входящих в закон движения (7 = 1.5 = 1, С = 1, 0* = 1°), будут использоваться в вычислительных экспериментах "по умолчанию", если иное не оговорено специально.
4. Результаты вычислительных экспериментов
Дня расчетов поверхностных волн, генерируемых движением оползня, использовалась схема предиктор-корректор па адаптивной сетке |5|, аппроксимирующая уравнения
мелкой воды (1), записанные в новых координатах (q,t):
(Ju)t + (f - XtU), = G.
Здесь J(q, t) — якобиан невырожденного преобразования координат
x = x(q,t), x(0,t) = 0, x(1,t) = L, (17)
t
зок [0, 1] на отрезок [0, L], Для построения адаптивной сетки использовался метод экви-раепределения [6], суть которого заключается в том, что из множества всевозможных отображений вида (17) выбирается то, которое является решением начально-краевой задачи
(w(x,t)xq)q = ^ q £ (0, 1), t> 0, , ,
х(0, t) = 0, x(1, t) = L, x(q, 0) = x0(q), 1 '
где в _ положительный параметр, подбираемый экспериментальным путем в целях уменьшения осцплляцпй траекторий узлов подвижной сетки, w(x,t) — заданная управляющая функция,
4.1. Тестовая задача о генерации поверхностных волн при движении тела по горизонтальному дну
В качестве тестовой рассмотрим задачу о генерации поверхностных волн твердым телом заданной формы (16), движущимся с постоянной скоростью vc > 0 по горизонтальному дну hbt(x) = — h°, Таким образом, в этой задаче функция, описывающая подвижное дно, изменяется со временем по следующему закону:
h(x, t) = h0 — h° (x + x° — xc(t)),
где xc(t) = x° +vct, В рассматриваемой тестовой задаче полная глубина H°(x) и скорость жидкости u°(x) в начальный момент времени подбирались таким образом, чтобы задача Коши для уравнений (1) с этими начальными данными имела решение в виде волны,
vc
u(x,t) = u°(x — vct), H (x,t) = H°(x — vct). (19)
Подставляя данное решение в уравнения (1) и требуя, чтобы на бесконечности жидкость покоилась, получаем, что начальные функции должны удовлетворять соотношениям
Н0(х) [щ(х) - vc] = -vcho, ^ + 9 [Но(х) - h(x, 0)] = (20)
Откуда следует кубическое уравнение для определения функции H°(x):
#o-#o(W) + |)+| = 0, (21)
где Н0 = Н0/h0, h = h/ho, vc = vc/л/gho- Уравнение (21) может не иметь положительных корней либо иметь один или два положительных корня, соответствующих докри-тичеекому и сверхкритическому течению. Это зависит от заданных значений скорости
движения тела ус и его толщины Т. Здесь рассматривается случай докритического течения, который реализуется при выполнении условия докритичности течения над вершиной движущегося тела:
ЗгГ - Щ Т < 1-----,
2
где Т = Т/Л,0, На рис, 2, а показана верхняя граница области докритического течения. Задавая значения Т > 0 и ус > 0, соответствующие этой области, находим больший положительный корень уравнения (21), по первой из форму:: (20) определяем начальную скорость и0(х), а затем по формулам (19) находим точное решение уравнений (1),
На рис, 3, а показаны графики свободной границы дня точного и численного решений в момент времени £ = 2000 с, к которому поверхностная волна в форме впадины прошла расстояние, примерно равное восьми своим длинам. В расчетах значения параметров х0, Ь брались такими же, как указано в разд. 3, скорость ус движения тела полагалась равной 20 м/с. Линии 1 на рис. 3, а получены при к0 = 200 м, Т = 25 м (на рис, 2, а этим данным соответствует точка 1). Видно, что несмотря па использование в расчетах сетки с небольшим количеством узлов N = 100, имеет место достаточно хорошая точность численного решения. Причиной этого является применение адаптивной сетки, сгущающейся в окрестности впадины (см, рис, 3, б) и отслеживающей движение волны. Отметим, что при использовании равномерной неподвижной сетки сравнимая точность достигается при на порядок большом количество узлов.
Линии 2 па рис, 3, а и точка 2 па рис, 2, а соответствуют движению тона но горизонтальному дну глубины к0 = 400 м. Видно, что двукратное увеличение глубины ведет к уменьшению амплитуды волны почти в три раза. Нелинейному уменьшению амплитуды волны способствует и уменьшение толщины тела (см, линии 3 па рис, 3, а, и точку 3 па рис, 2, а, полученные при к0 = 200 м, Т = 12.5 м). Нелинейный характер зависимостей минимального значения Птш во впадине волны от входных данных демонстрируется на рис, 2, б, где показаны графики изменения амплитуды волны Птш ПРИ варьировании
а б
Рис. 2. Граница области докритичеекшх) течения в плоскости параметров Т и ус (а) и графики зависимостей (Т) (!), ("Сс) (2), ^т;п(М (3) (б)
а б
Рис. 3. Графики точншх) (сплошные .линии) и чиеленших) (штрих) решений в момент времени £ = 2000 с (а) и траектории узлов адаптивной сетки (б), N = 100
одного из трех параметров Т, ус ют и Л,0 и неизменных значениях двух других, которые полагались равными Т = 25 м, ус = 20 м/с и Л,0 = 200 м. Заметим, что амплитуда генерируемой волны не зависит от длины тела Ь, влияющей лишь на длину волны.
Хорошее согласование полученных численных результатов с точным решением тестовой задачи позволяет надеяться па достоверность приведенных ниже расчетных данных дня задачи о генерации воли оползнем, движущимся по неровному дну.
4.2. Общие характеристики волновых режимов,
порождаемых движением оползня по неровному склону
Перейдем к рассмотрению особенностей волнообразования при движении оползня по неровному склону (15). Неровное дно оказывает влияние пе только па характер трапе-формации воли при их распространении, по и па закон движения оползня. На рис. 4, а показаны траектория движения вершины оползня (линия 1), перемещающегося под действием сил тяжести, выталкивания, трепия и сопротивления окружающей жидкости, и график изменения скорости движения оползня (линия 2). На профиле рельефа (линия 3) два крайних кружка изображают начальное (t = 0 с, xc = 12 644 м) и конечное (t = 1370 с, xc = 37 233 м) положения оползня, а средний — положение (xc = 22 624 м) оползня в момент (t = 512 с) достижения максимальной скорости (vc ~ 37.5 м/с) его движения, при этом в данный момент времени оползень находится ниже точки перегиба (£ = 18 673 м) профиля дна. Как видно из графика скорости, промежуток времени, соответствующий разгону оползня (Ait = 512 с), существенно меньше промежутка торможения (A2t = 858 с). Примерно в таком же отношении находятся пути, пройденные оползнем па участках разгона и торможения. По графику скорости также видно, что остановка оползня происходит резко (при ненулевом отрицательном значении ускорения) .
На рис, 4, б" изображены профили свободной границы в моменты времени, близкие к моменту достижения оползнем максимальной скорости. Точками здесь показаны положения вершины оползня в эти моменты времени на траектории ее движения (на линии х = хс(£)), Из приведенных графиков следует, что максимальное отклонение свободной границы от певозмущеппого состояния достигается точно над вершиной оползня и образовавшаяся здесь "впадина" движется вместе с оползнем, при этом в отличие от тестовой задачи амплитуда "впадины" со временем меняется. Перед оползнем формируется волна повышения, которая толкается оползнем в сторону увеличения глубин и иостеиенно отделяется от него с образованием на заднем фронте волны понижения (линия 3 па рис, 4,6), Слова от "виадииы", сопровождающей оползень, также формируется волна понижения, которая движется в сторону берега.
Указанные особенности оползневого механизма волнообразования хорошо иллюстрируются рис, 5, а, где изображена поверхность г = п(х,£). Видно, что волна понижения, двигающаяся к берегу, отражается от расположенной здесь вертикальной стенки в виде волны повышения, уходящей вправо. Рассмотрение динамики свободной границы на более продолжительные промежутки времени показывает, что колебания уровня воды на .новой стенке продолжаются и поело отражения набегающей волны понижения, однако максимальный вертикальный заилеек происходит во время взаимодействия со стенкой именно этой волны. Видно также, что при торможении оползня амплитуда идущей над ним "впадины" монотонно уменьшается и при остановке оползня формируются две противоположно направленные волны с небольшими амплитудами.
На рис. 5, б представлены траектории узлов адаптивной сетки, которая, как и в рассмотренной выше тестовой задаче, сгущалась в зонах максимальных значений вони-чины |п(х, ¿)|. Отчетливо видно, что равномерная в начальный момент сетка начинает перестраиваться с ростом времени в неравномерную, автоматически отслеживающую движущиеся волны повышения и понижения. На рисунке штриховой .пинией нанесена траектория движения вершины оползня, согласно которой при больших скоростях его
Рис. 4. Графики функций х = хс(Ь) (1), V = vc(t) (2) и г = (х) (3) (а) и графики свободной границы в моменты времени Ь = 300 с (1), 500 с (2) и 700 с (3) (б)
Рис. 5. Динамика свободной поверхности (а) и траектории узлов адаптивной сетки (б)
движения сетка имеет сильное сгущение в окрестности этой траектории и адаптируется к движущейся "впадине". При уменьшении амплитуды "впадины" сгущение рассасывается и сетка адаптируется к другим волнам, имеющим преобладающую амплитуду.
4.3. Зависимость волновых процессов от размеров оползня
Дня изучения влияния на волновые процессы таких характеристик как крутизна склона, длина, толщина и начальное заглубление оползня были выполнены вычислительные эксперименты, в которых один из параметров варьировался, а значения всех остальных выбирались такими же, как указано в раздело 3.
Вначале рассмотрим зависимость генерируемого оползнем волнового режима от толщины оползня Т. Поскольку объем модельного оползня (16) определяется по формуле V = РРТ6/2, то величина Т в качестве множителя будет входить во все слагаемые уравнения (13). Следовательно, закон движения оползня не будет меняться при изменении Т
Т
На рис. 6, б показаны мареограммы, записанные в двух мареографных точках. Мареограф А расположен на берегу и фиксирует высоту свободной поверхности на вертикальной стенке (в точке х = 0). Подвижный мареограф В перемещается вместе с оползнем, имеет абсциссу х = хе(£) и измеряет глубину "впадины" над оползнем. Видно, что эта
Т
ствии с решением тестовой задачи (см. поведение линии 1 иа рис. 2,6). Отметим, что
Т
ной величины на участке разгона, причем это происходит гораздо раньше момента достижения оползнем максимальной скорости. Рост глубины "впадины" при возрастании толщины оползня способствует увеличению амплитуды волны понижения, приходящей к берегу и, как следствие, к росту вертикального заилеска (см. мареограммы А), при
Т
а б
Рис. 6. Влияние толщины оползня Т на характеристики движения оползня и генерируемых им волн; Ь = 5000 м, (1с = —200 м, 0о = 6°; а — графики функций ж = хс(£), V = -ис(£) и г = (ж); ................. мареограммы А и В; Т = 5 м (1), 15 м (2), 25 м (3), 35 м (4)
Ь
волн; Т = 25 м, (1с = —200 м, 0о = 6°; а — графики функций х = жс(£), V = vc(t) и г = ^(х); ^................. мареограммы А и В; Ь = 3000 м (1), 5000 м (2), 7000 м (3), 9000 м (4)
Ь
Т
Ь
с большой скоростью и, как следствие, проходит большее расстояние по склону, чем менее протяженный (рис. 7, а, здесь старт оползня показан точкой).
4.4. Зависимость волновых процессов от начального заглубления оползня и крутизны склона
На рис, 8, а маркерами указаны точки старта оползней с различными величинами начального заглубления йс. Видно, что при любых ¿с точки остановки оползней практически не различаются, по вместе с тем заметно различаются времена движений и максимальные значения скоростей: дня "мелководных" оползней они больше, чем дня "глу-
а б
Рис. 8. Влияние начального заглубления оползня <1с на характеристики движения оползня и генерируемых им волн; Т = 25 м, Ь = 5000 м, во = 6°; а — графики функций ж = хс(£), V = vc(t) иг = к\л (х); б — мареограммы А и В; ¿с = —200 м (!), —300 м (2), —400 м (3), —500 м (4)
Рис. 9. Влияние крутизны склона на характеристики движения оползня и генерируемых им волн; Т = 25 м, Ь = 5000 м, = — 200 м; а — графики функций х = хс(£), V = vc(t) и г = Нл(х);
во = 6° 9 12 15°
бо ко водных". Соответственно h амплитуды генерируемых волн больше для оползней, стартующих с меньшей глубины (рис, 8, б), что качественно согласуется с результатами решения тестовой задачи (см, поведение линии 3 на рис, 2, б).
На рис, 9, a показаны графики траекторий оползней и их скоростей при разных значениях угла наклона 0о профиля дна в точке перегиба. Поскольку начальное заглубление оползней одинаково, то на более крутых склонах точки старта оползней (обозначены маркерами) располагаются ближе к берегу, чем на менее крутых. Большая крутизна склона приводит к тому, что развиваемая оползнем скорость принимает большие значения, а участок разгона и время движения — меньшие, чем эти величины при более пологом склоне, С ростом крутизны склона возрастают и заплески на стенку (см, ма-реограммы А на рис, 9, б).
5. Заключение
В представленной работе получены уравнения, описывающие движение подводного оползня по криволинейному склону под действием сил тяжести, плавучести, трения и сопротивления воды, С использованием предложенного закона движения выполнены расчеты волновых режимов при различных значениях начального заглубления модельного оползня, его размеров и крутизны модельного склона. Результаты численного моделирования подтвердили предположение о том, что наибольшую опасность для береговых сооружений представляют волны, генерируемые движением оползней больших размеров вблизи берега на крутых склонах, В дальнейшем полученный закон движения оползня будет использован для расчета генерируемых поверхностных волн на основе нелинейно-дисперсионных и полных гидродинамических моделей.
Список литературы
[1] Шокин Ю.И., Чуваров Л.Б. О подходах к численному моделированию оползневого механизма генерации волн цунами // Вычисл. технологии. 2006. Т. 11, Спец. выпуск, часть 2. С. 100-111.
[2] Watts Р., 1мamura F., Grilli S. T. Comparing model simulations of three benchmark tsunami generation cases // Science of Tsunami Hazards. 2000. Vol. 18, No. 2. P. 107-123.
[3] Елецкий C.B., Майоров Ю.Б., Максимов В.В. и др. Моделирование генерации поверхностных волн перемещением фрагмента дна по береговому склону // Вычисл. технологии. 2004. Т. 9, Спец. выпуск, часть 2. С. 194-206.
[4] Watts P., Grilli S.T., Kirby J.T. et al. Landslide tsunami case studies using a Boussinesq model and a fully nonlinear tsunami generation model // Natural Hazards and Earth System Sciences. 2003. Vol. 3, No. 5. P. 391-402.
[5] Шокин Ю. И., Хакимзянов Г. С. Схема предиктор-корректор, сохраняющая гидравлический скачок // Вычисл. технологии. 2006. Т. 11, Спец. выпуск, часть 2. С. 92-99.
[6] Численное моделирование течений жидкости с поверхностными волнами / Г.С. Хакимзянов, Ю.И. Шокин, В.Б. Барахнин, Н.Ю. Шокина. Новосибирск: Изд-во СО РАН, 2001. 394 с.
Поступила в редакцию 20 января 2010 г.
