CC BY
143
Вычислительные технологии
Том 13, Специальный выпуск 2, 2008
Моделирование динамики наносов в прибрежной зоне
И, С, Нуднер, В, В, Максимов, Ю.Б, Майоров, Д. В, Бабчик Научно-исследовательский центр 26 ЦНИИ Минобороны РФ, Санкт-Петербург, Россия е-таП: итах1тоу@таИ. ги, д.пис1пег@тад.1. ги
In this paper some theoretical and experimental results addressing the important problem of inshore sediment dynamics are presented. Main features of the three modem software packages solving this problem are considered. A detailed list of references (in Russian) relevant to this problem is given.
Введение
С проблемой заносимости акваторий морских портов и изменения береговой линии человечество знакомо с давних лет. Еще древние финикийцы при строительстве портов Тир и Силон ( X111 XII вв. до и. э.) стремились учитывать естественное движение наносов. В древнем Риме в ряде античных портов были построены сквозные молы для пропуска вдольбереговых наносов. Известно, что Леонардо да Винчи занимался лабораторным исследованием волнового движения. По-видимому, он первым из европейцев установил орбитальный характер движения жидких частиц в волне. Великий флорентиец имел ясное представление о существовании продольного перемещения береговых наносов и учитывал его при проектировании портов Чезаматико (Италия) и Гавр (Франция).
Как правило, решение проблемы заносимости акваторий, являющейся частью более общей проблемы динамики наносов, долгое время осуществлялось с помощью натурных наблюдений. Лишь в начале XX в. начали предприниматься попытки решения этой сложной задачи теоретическими методами. В Советском Союзе активное решение задач переноса донных отложений началось примерно с середины XX в. Под руководством В. П. Зенковича была создана первая научная школа по изучению геоморфологии берегов. К решению различных задач динамики наносов стали применяться как экс-нериментальные, так п теоретические методы [1-47]. В настоящее время удалось но-лучить теоретические решения ряда частных задач, подкрепив их экспериментальной проверкой.
Все основные изучаемые процессы взаимодействия волн, течений и движения наносов происходят в прибрежной зоне. Здесь происходит изменение параметров набегающих волн: начинает проявляться асимметрия движения. При подходе к берегу возникает рефракция волн, увеличение и укручение их гребней, что в большинстве случаев приводит к их обрушению (которое может происходить многократно). Волновое движение носит ярко выраженный турбулентный характер. Происходят процессы эрозии (размыва) и аккумуляции (накопления) берегового склона. Под действием волн, подходящих к
© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2008.
берегу под косым углом, и вдоль береговых течений появляется продольный транспорт наносов. Возникают прибрежные разрывные течения. Донные отложения, вовлекаемые в движение, имеют разнородный состав, от валунов (диаметром более 100 мм) до мелкого пелита (диаметром менее 0.001 мм). Промежуточный диапазон значений занимают галька, гравий, песок, алеврит.
При теоретическом решении задачи транспорта наносов делаются многочисленные упрощения (чтобы сделать задачу принципиально разрешимой) и вводятся дополнительные гипотезы (для замыкания получающихся систем уравнений). В немногочисленных случаях (например, [29-32]) удается получить удовлетворительное совпадение теоретических результатов с экспериментальными и натурными данными. Чаще же результаты теоретических моделей дают лишь качественное совпадение, различаясь по величине с экспериментальными данными в несколько раз. Указанное обстоятельство характеризует сложность и неизученность проблемы динамики наносов в целом, несмотря на предпринимаемые интенсивные попытки найти ее решение.
В данной работе кратко обсуждается современная классификация теоретических моделей, рассматривается методика проведения экспериментальных исследований заносимое™ акватории и проводится анализ ряда полученных результатов физического моделирования для акватории конкретного заданного объекта. Приводится подробное решение одной из важных задач геоморфологии берегов: задачи о переформировании береговой линии. Обсуждается три наиболее разработанных пакета программ, которые в принципе позволяют получить приближенное численное решение задачи о динамике наносов.
1. Классификация теоретических моделей динамики взвешенных наносов
В настоящее время существует три группы моделей, описывающих движение взвешенных наносов.
Первая группа моделей рассматривает движение твердой одиночной частицы. Размеры частицы полагаются малыми. Предполагается, что коэффициент сопротивления линейно зависит от относительной скорости частицы и окружающей ее жидкости. Учет нелинейности коэффициента сопротивления резко усложняет вычисления, которые приводят к громоздким выражениям и трудно алгоритмизуемы. Кроме того, остается неясным, как переносить результаты, полученные для одиночной частицы, на ансамбль частиц с учетом их взаимодействия между собой. Впрочем, имеются результаты, указывающие на то, что в большинстве случаев при достаточном удалении частиц друг от друга (один-два диаметра частиц) их взаимодействием можно пренебречь.
Вторая группа моделей вводит в рассмотрение систему из большого, но конечного числа частиц. Фактически они сводятся к начально-краевой задаче, записанной системой дифференциальных уравнений, которая описывает движение жидкости и частиц, находящихся в ней. Однако в настоящее время в рамках этих моделей рассматриваются лишь отдельные частные гипотезы о характере силового взаимодействия жидкости и твердых частиц. Кроме того, до сих пор остается неясным, как происходит взаимодействие движущихся частиц потока с твердыми границами. Поэтому существенное продвижение в решении задачи о движении наносов многочастичными моделями представляется проблематичным.
Третья группа моделей описывает движение потока жидкости со взвешенными наносами в виде двухфазных континуальных сред, взаимодействующих между собой. Уравнения для сред вводятся двумя способами: осреднением уравнений, выполняющихся для фаз, рассматриваемых в отдельности, или введением гипотез, получаемых из более общих соображений относительно возможного поведения двухфазной среды. Эта группа моделей в настоящее время применяется достаточно часто при построении вычислительных алгоритмов, несмотря на то, что можно указать большое количество ее слабых мест: корректное описание коэффициентов турбулентной диффузии; учет обратного воздействия взвеси на движение жидкости; большой произвол в выборе начальных и граничных условий; учет саморегуляции потока, содержащего взвесь, с дном и некоторые другие обстоятельства.
Фактически при выполнении расчетов динамики наносов сначала рассчитывают волновой режим в исследуемой акватории, а затем найденные гидродинамические характеристики вводят в уравнения для расчета концентрации. Это уравнения конвекции-диффузии. Следует отметить, что они справедливы для малых концентраций и частиц небольшого размера. (На практике производится расчет по каждой фракции взвеси отдельно и результаты осредняются в соответствии с теми или иными предположениями.) В случае больших изменений батиметрии и береговой линии волновой режим пересчитывается и процедура повторяется.
2. Экспериментальные исследования заносимости акватории объекта
Необходимым предварительным этапом исследований заносимости является определение волнового режима на акватории объекта. Оно выполнено на основе расчетно-теоретических и экспериментальных исследований. Рассмотрены все варианты компоновки оградительных сооружений. При их оценке основным требованием было уменьшение высоты волн на внутренней акватории до заданной величины. При определении защищенности акватории и выдаче рекомендаций по плановой компоновке принималась заданная обеспеченность расчетного шторма. Применялось моделирование с разными горизонтальным и вертикальным масштабами, что потребовало введения в экспериментальные результаты для высот волн коэффициентов искажения дифракции, получаемых расчетным путем. Экспериментальные исследования волнового режима проводились на неразмываемой пространственной модели объекта, выполненной на мелководной ступени большого опытового бассейна. Горизонтальный масштаб моделирования составил 1: 100, вертикальный — 1: 25. Это обстоятельство обусловлено возможностью размещения модели в бассейне и обеспечения необходимых параметров волнения. Волнение создавалось передвижным волнопродуктором балочного типа. Опыты проводились для наиболее волноопасного направления. Параметры волн регистрировались на подходе к акватории, в канале и на самой акватории. Полученные значения использовались далее для теоретической оценки заносимости и сравнения с данными экспериментов.
2.1. Задачи, состав и методика экспериментальных исследований
Основной задачей исследований являлась разработка инженерных решений по уменьшению заполнимости объекта наносами. Экспериментальные исследования были на-
правлены на изучение эффективности вариантов защиты акватории объекта от наносов, образующихся при комплексном воздействии на подводный склон волн и течений, и основывались на методе сравнительного анализа вариантов защиты.
Эксперименты проводились на пространственной установке, расположенной в мелководной ступени опытового бассейна, с использованием модели объекта, на которой предварительно были проведены исследования характеристик волнового режима, после покрытия рабочей части модели слоем песка толщиной 10 см со средним диаметром частиц Б0 = 0.22 мм (рис. 1). Песчаное основание модели подходного канала было покрыто цементным раствором. Тем самым создавалась возможность количественной оценки объема осевших в подходном канале наносов. Отметка воды, заполняющей установку, соответственно была увеличена на 10 см с целью сохранения глубин на входе в канал и составила 50 см.
За экспериментальный сезон проведено девять опытов (рис. 2). За основной принят вариант существующей в настоящее время компоновки акватории (опыт 1). Особое внимание в исследованиях обращено на варианты ограждения подходного канала парными параллельными молами. Было исследовано два случая: распределение скоростей течения воды по подводному склону и волнения, соответствующего северо-восточному румбу в натуре. Один случай распределения скоростей (опыт 6) аналогичен распределению скоростей по подводному склону (опыт 1). Другой (опыт 3) — с максимумом скоростей, смещенным на большие глубины, и числами Фруда Ег = 0.022 ... 0.036. В аналогичных условиях рассмотрены варианты с ограждением канала парными молами и усиленными подводными банкетами, развернутыми относительно оси канала на угол 25° (опыты 2 и 7).
Рис. 1. Компоновка акватории: 1 — моделируемая акватория; 2— волномеры; 3— волноиро-дуктор; 4 ^ выгородки из щитов; 5 — водобойный колодец; 6 — выгородки из железобетонных блоков
Рис. 2. Исследуемые варианты. Масштаб 1: 50
Ограждение канала парными молами увеличенной длины моделировалось в условиях течения с числами Фруда Гг = 0.023 на изобате 16 см и Гг = 0.035 — на изобате 24 см (опыт 4).
Волнение во всех вариантах имело северо-восточное направление.
Далее исследовался вариант с удлинением и разворотом южного мола в восточном направлении (опыт 5) и удлинением северного мола параллельно оси канала. Числа Фруда на тех же изобатах равны соответственно Гг = 0.025 и Гг = 0.031.
Опыты 8 и 9 — это варианты с разворотом оси канала на 30 и 10° относительно своего первоначального положения; продолжение южного мола имело ломаный профиль.
У голов молов были расположены подводные банкеты. Числа Фруда на упомянутых выше изобатах равны соответственно К- = 0.023 и К- = 0.02 ... 0.03.
Поскольку заносимость объекта обусловлена как вдольбереговым движением наносов, так и движением наносов по нормали к нему, моделирование их представляет собой задачу повышенной сложности. Задача решалась следующим образом. За основу было принято качественное сопоставление рассматриваемых вариантов защитных мероприятий. Каждый вариант сравнивался с базовым (существующая в настоящее время компоновка акватории) по количеству наносов, осевших на акватории и в подходном канале. Затем определялась деформация подводного склона на входе в канал.
Поскольку годовую заносимость натурного объекта моделировать практически невозможно, то выполнялось моделирование лишь части этой заносимости, по времени соответствующей длительности расчетного шторма. Используемый метод качественного сопоставления предполагает равенство значений времен проведения всех опытов, причем это время не должно превышать значения, при котором происходит предельное заполнение наносами подходного канала в опыте с существующей в настоящее время компоновкой объекта. Это время, оказавшееся равным 1.5 ч, и было принято в качестве времени проведения каждого опыта.
Параметры волн выбирались по результатам исследований волнового режима и были следующими: К = 7.5 см и Т = 1.33 с. Таким образом, число волн, воздей-
1.5
^т0с1 = ^/Ттоа = 1.5 ■ 3600/1.33 = 4100. В натуре среднее число волн в случае шторма с северо-восточного направления, повторяемостью один раз в 10 лет, составляет ^ = и/Тп.АЪ = (14 ... 16)3600/7.2 = (7 ... 8)103.
Течение на модели обеспечивалось с помощью установки, включающей насосную станцию, закачивающую слабосоленую воду с глубоководной ступени опытового бассейна в систему трубопроводов, и водовыпуски. При достижении расхода воды, подающейся па модель, максимального значения ^ = 0.36 м3/с включался волнопродуктор. Этот момент считался началом опыта.
В течение каждого опыта с помощью микровертушек осуществлялись промеры скоростей течения воды в намеченных створах, а также производилась запись волн четырьмя волномерами. После окончания опыта установка осушалась, производилась оценка заносимости подходного канала и акватории. Для этого после каждого шага опыта собирались наносы, осевшие в подходном канале и на акватории объекта, которые далее высушивались и взвешивались.
Для проведения линии деформации подводного склона в районе начала подходного канала визуально определялась область локального уменьшения глубин, обычно расположенная в начале подходного канала. Затем намечалась ось максимального уменьшения глубин (ось подводной косы), вдоль которой затем определялся профиль этого аккумулятивного образования.
2.2. Результаты экспериментальных исследований
В опыте 1 исследовалась заносимость подходного канала в варианте существующей в настоящее время компоновки объекта. Значения чисел Фруда, соответствующие средним скоростям на вертикалях, в этом опыте были несколько больше чисел Фруда для аналогичных вертикалей в натурных условиях. Во время проведения эксперимента на подводном склоне наблюдалось грядовое перемещение наносов, которое привело
к полной заносимостп подходного канала. Вес осевших в подходном канале наносов составил 314.7 кг.
В опытах 6 и 3 изучалось влияние удлинения южного мола на заносимость подходного канала. Результаты показали, что заносимость подходного канала существенно зависит от распределения скоростей по подводному склону. Если после проведения опыта 3 заносимость подходного канала практически отсутствовала, то после проведения опыта 6 вес осевших в канале наносов был равен 118.3 кг. Вместе с тем удлинение южного мола позволило снизить заносимость подходного мола на 62 % по сравнению с заносимостью подходного канала в варианте существующей в настоящее время компоновки объекта. Однако удлинение южного мола вызвало увеличение скоростей у головы этого мола, что способствовало образованию воронки размыва в этом районе и росту подводной косы. При этом локальное уменьшение глубин на входе в канал в опыте 6 оказалось больше, чем в опыте 3.
Таким образом, распределение скоростей по подводному склону сказывается и на размерах подводной косы, образующейся на входе в канал. Следует заметить, что в опыте 6 наблюдалось рифельное движение наносов вдоль подходного канала на акваторию, хотя поступление наносов на акваторию объекта за время опыта было незначительным.
С целью уменьшения заносимости подходного канала и уменьшения размеров подводной косы, образующейся на входе в канал, проведены опыты 2 и 7, в которых к удлиненным молам достраивались усиленные банкеты из каменной наброски, развернутые от канала относительно его оси на 25 %. Распределение скоростей по подводному склону в опыте 2 было аналогично распределению скоростей в опыте 3. Эпюра скоростей в опыте 7 была аналогична эпюре скоростей в опыте 6.
Результаты экспериментов показали, что вариант комплексной защиты (опыт 7) поз-
88
существующей в настоящее время компоновки объекта (опыт 1).
Распределение средних скоростей по подводному склону влияет на заносимость подходного канала и размеры аккумулятивных форм, образовавшихся за подводными банкетами. Сделанный вывод подтверждает то, что если в опыте 7 появилась подводная коса меньших размеров, чем в опыте 6, то в опыте 2 четко выраженного аккумулятивного образования не отмечено. В этих опытах наблюдалось рифельное движение наносов вдоль канала, однако заносимость акватории за время опыта была незначительной.
Дальнейшие эксперименты были посвящены вопросам уменьшения заносимости за счет увеличения длины оградительных сооружений и изменения их планового размещения (опыты 4, 5, 8 и 9).
В результате этих экспериментов установлено, что заносимость подходного канала и акватории после проведения этих опытов практически отсутствовала, однако у голов южных молов образовывались аккумулятивные формы в виде подводных кос, которые под разными углами пересекали ось подходного канала.
Наибольших размеров подводная коса образовывалась в опытах 4 и 5. Подводная коса наименьших размеров возникла после опыта 8.
Подводные косы, образовавшиеся в опытах 4 и 5, располагались в основном в пределах канала между южным и северным молами, т. е. перекрывали собой вход в канал.
В опытах 8 и 9 подводные банкеты предотвращали размыв головы южного мола; в опыте 4 размыв подводного склона наблюдался за банкетом. В опыте 5 сильному размыву подвергся участок разворота южного мола в восточном направлении (рис. 3).
а
Рис. 3. Опыт 5: а вид на модель объекта пород опытом; б вид на южный мол во время проведения опыта; в поело проведения опыта
Таким образом, результаты экспериментов позволяют сделать следующие выводы.
1. Усиленные каменные банкеты снижают размыв грунта у голов молов и уменьшают размеры подводных аккумулятивных образований в районе входа в канал.
2. Наибольшую защищенность от вызванной штормами заносимости наносами в восточной четверти подходного канала и акватории объекта обеспечивает развитие южного мола по вариантам 8 и 9.
3. Теоретическая модель изменения береговой линии
Лито динамические процессы, протекающие на участках побережья с прилегающей к ним акваторией, часто приводят к переформированию дна и линии берега. Поэтому при выборе мест строительства новых портов необходимо учитывать характер таких процессов. Несмотря на то что до возведения морского гидротехнического объекта проводятся определенные изыскания, а при проектировании предусматриваются необходимые меры защиты, проблемы, связанные с отложением наносов и размывом, решаются трудно и требуют значительных затрат при последующей эксплуатации комплексов морских гидротехнических сооружений. Поэтому для прикладной гидромеханики имеет большое значение исследование как самих процессов динамики наносов, так и процессов изменения береговой линии и форм профиля подводного берегового склона.
ЗЛ. Изменение береговой линии при отложении наносов
В естественных условиях дно моря представляет собой слой несвязного грунта из песка, гальки, ила или их смеси. При отсутствии волн можно говорить о поверхности раздела вода-грунт. При волнении возникает поток частиц из донного слоя через поверхность раздела в водный слой. В результате этого частицы грунта перемешиваются с водой. В водном слое частицы грунта находятся во взвешенном состоянии. Крупные, тяжелые частицы не переходят поверхность раздела.
В случае твердого недеформированного дна береговая линия задается уравнением
H(x, y) = 0, (1)
которое является неявным заданием функции y = f (x), определяющим линию берега. Ординаты поверхности раздела определяются функцией
z = —D(x, y, t), (2)
которая заранее не известна.
При волнении в одних местах береговой зоны происходит аккумуляция наносов, а в других — размыв грунта.
В случае нетвердого дна, верхняя граница которого является поверхностью раздела z = —D(x, y, t), линия берera y = n(x, t) будет определяться уравнением
D(x, y, t) = 0. (3)
Таким образом, имеется два слоя. Верхний, водный, со взвешенными в нем частицами (наносами) и нижний, донный, слой крупных и мелких частиц, насыщенный водой. Донный слой имеет сложную реологию. Задача исследования динамики как водного, так и донного слоев является не только сложной, но и не определенной.
В настоящее время математическое моделирование процесса воздействия волн на береговую зону сводят к составлению и решению уравнения, описывающего форму береговой линии y = t) с учетом соответствующей гипотезы о векторе потока наносов. Предполагают, что поперечный профиль берега принимает форму устойчивого равновесия. В этом случае поперечным перемещением наносов можно пренебречь. Это означает, что формирование (изменение) линии берега происходит только за счет вдольберегово-го потока наносов.
Рассмотрим сечение берега плоскостью x = £ и в нем область Б, ограниченную вертикалью в точке (£, n(£, t)). Полагаем, что ось x направлена вдоль берега, ось y — перпендикулярно в сторону моря, ось z — вертикально вверх. Система осей Oxyz — правая.
В момент времени площадь области Б определяется величиной
Б(£,t)= У H(£,y) dy, (4)
о
причем
% = и
Рассмотрим объем участка берега, равный Б(£, t)dx. За единицу времени этот объем изменится на величину
дБ , ,
—— ax at. dt
Такое изменение объема вызвано потоком наносов, равным
Qx = -Q cos а, (6)
направленным вдоль осп х.
Согласно закону баланса потока наносов имеем равенство
Н(х, = ЯхЦ+Лх. (7)
Из выражения (7) следует уравнение
дп д
Н{х,г))— = - — {дсоъа). (8)
Здесь Q — поток наносов вдоль линии берега; а — угол между касательной к линии берега и осью х; tg а = дп/дх. Расход наносов, формирующих линию берега, вызывается действием волн. Величина расхода определяется скоростью движения частиц наносов и зависит от очертаний линии берега. В формуле для расхода Q это обстоятельство учитывается зависимостью угла в между фронтом набегающей волны и касательной к линии берега [14, 15].
Пусть угол 7 определяет направление фронта волны (рис. 4). Тогда можно записать
в = 7 - а. (9)
/ у
Рис. 4. Схема к задаче об изменении береговой линии
Учитывая равенства
dQ _ ^a _ Vxx д/ _ 0Q
dx dx' dx 1 + ' дб
запишем уравнение для определения линии берега (8) в виде
дп
Н(х, Tj)— = \_Q'{7 — cos a + Q(7 — a) sin a]
1 + nL'
Представляя Q в гаде Q(0) = B sin 20, из уравнения (11) получим
дп
Н{х,г])-т^ = б|2 сое а [сое 27(1 - tg2 о) + 2sin27tga] +
+ вит [вт 27(1 — tg2 а) — 2 сое 27 tg а] —•
При малых значениях угла а уравнение (12) существенно упрощается:
(10)
cos а = cos(y — 0). (11)
(12)
(13)
Уравнение (13), как и более общее уравнение (12), нелинейное. Решать их нужно численными методами.
В случае крутого подводного склона можно считать, что H(x,y) = const. Тогда уравнение (13) будет линейным:
дп 2 д2п 2
, а =-— cos 27.
дt а дх2
2В ~Н
(14)
Для определения изменения береговой линии уравнение (14) решается при начальном условии
П(х, 0) = по(х), |х| < то, Ь> 0. (15)
Решение задачи (14) и (15) будет иметь вид
сю
n(x,t)= f G(xM)По(С) dt, G(xM)
1
л/47таН
exp
о* - е
АаЧ
(16)
Рис. 5. Береговая линия треугольной и прямоугольной бухт для трех моментов времени
Задавая функцию По (я), можно достаточно просто оп ределить п (х,£) — береговую линию.
Модельные примеры эволюции береговой линии бухт, имеющих первоначально треугольную
п (г)( |х| - 100, |х| < 100, По(х) = 1 0, |х| > 100,
и прямоугольную форму
по(х) =
{
-50, 0,
|х| < 10, \х\> 10,
(17)
(18)
под действием волн, приведены на рис. 5. В обоих случаях первая форма береговой линии (обозначена кружками) получена для некоторого момента времени, близкого к начальному. Вторая линия (крестики) была вычислена для момента времени, в 10 раз большего. И, наконец, третья линия (квадратики) рассчитывалась во время, в 10 раз большее времени вычисления второй линии. Из рисунков видно, что линейная модель в целом описывает процессы аккумуляции и абразии берегов, но не учитывает многих процессов, имеющих место при транспорте наносов.
3.2. Отложение наносов у мола
Рассмотрим задачу о защите акватории порта от наносов. Один из вариантов решения задачи — возведение мола поперек движения наносов. Мол представляет собой твердую вертикальную стенку и играет роль нанасозадерживающего сооружения (рис. 6). На плоскости хОу мол изображается в виде отрезка вдоль оси у длины Ь:
х = 0, 0 < у < Ь.
(19)
Пространство между линией берега и молом называют входящим углом. Здесь происходит отложение наносов, которые изменяют линию берега с течением времени.
Задачу о заполнении входящего угла рассматривают для двух стадий [15]. Граничным между стадиями является момент времени когда наносы достигают оголовка
Рис. 6. Схема к задаче о заносимости входящего угла
мола. Для первой стадии (0 < £ < Ь1) задача сводится к решению уравнения баланса наносов при условиях
П (х, 0) = по (х); (20)
дж = 0, х = 0, 0 < £ < £1, п(0, £1) = Ь П(х, £) ^ 0, х ^ то, £ > 0.
С момента времени £ = £1 условия задачи имеют вид
П(х,£1) = П1 (х) = п(х,£1 - 0); п(0, £) = Ь, £ < £1.
(21) (22)
(23)
(24)
Можно численно решить задачу о заполнении входящего угла при малых значениях а. В этом случае используются уравнение баланса наносов (13) и граничные условия
дп 1
—— = - tg 27, х = 0, г)(х, —> 0, х —> ос.
(25)
Рассмотрим частный случай задачи. Будем считать, что береговой склон крутой. Уравнение и условия задачи записываются следующим образом:
п = а2пхх, х > 0, 0 < £ < £1; п (х,£) = по (х), £ = 0; ?7Ж(0, = ^ = ж 0.
Решение удается получить в явном виде. Введем новую функцию
(26)
(27)
(28)
V = пх - М
Тогда функция V будет решением задачи
V = а2V,,,,, х > 0, 0 < £ < ,
(29)
Решение задачи имеет вид V 1
\J А-каЧ
exp
г (х-С)2] г (х + £)21
1 АаЧ \ — exp 1 АаЧ \
Vo(0 d£.
Учитывая, что
1
л/АттаЧ
exp
г (х-С)2] г (х + £)21
1 АаЧ \ — exp 1 АаЧ \
dC = Ф
x
\/АаН' '
(31)
(32)
Ф(ж) = —¡= / ехр(—a2) da, 'п J o
получим
Пх = ^
1 - ф
x
+
VAaH' у/АтгаН
exp
г (х-С)2] г (х + £)21
1 АаЧ \ — exp 1 АаЧ \
dx
dC
o
(33)
(34)
В результате интегрирования уравнение береговой линии примет следующую форму:
„2
r](x, t) = ¡л,{ 2а\ — exp (-1— х V п \Аа2tJ
1ф
x
2а\Д
+
г (Р-О2] г (.Р + О2]
1 АаЧ \ — exp 1 АаЧ \
oo
exp
x > 0, 0 < t < ti, n(0, ti) = L.
dx
dCdp
(35)
те
те
z
те
1
x те
4. Ведущие пакеты прикладных программ, доставляющие численное решение проблемы динамики наносов
Существует большое количество пакетов, позволяющих в той или иной мере численно решать проблему транспорта наносов. Мы остановимся на трех пакетах, на наш взгляд, наиболее полно подходящих к решению задачи. Следует отметить, что для успешного их использования необходимо наличие опыта в решении задач по прикладной гидромеханике и в подготовке исходных данных, выборе основных модулей, корректном задании начальных и граничных условий.
LITPACK. Второе название этого пакета — Интегрированная система моделирования и кинетики береговой линии. Разработчиком его является датская фирма DHI, Water and Environment. Он предназначен для моделирования перемещения несвязанных наносов при наличии волн и течения, литорального дрейфа, эволюции береговой линии и профиля пляжей. Пакет имеет графический интерфейс, полностью интегрированный в операционную систему Windows. LITPACK имеет два основных модуля: LITSTP — для вычисления перемещения наносов и LITDRIFT — для расчета вдоль
берегового перемещения наносов и литорального дрейфа. Три дополнительных модуля расширяют функции пакета: LITLINE позволяет рассчитывать изменение береговой линии, LITPROF — изменение поперечного профиля, LITTREN вычисляет отложение наносов в подводных желобах и траншеях. Все изменения можно выполнять в оперативном режиме. Пакет позволяет наносить на карту различные типы гидротехнических сооружений. Возможен динамический вывод промежуточных результатов.
SMS. Название этого пакета является сокращением от Suface-water Modeling System. Текущая версия 9.2. Разработчик — лаборатория окружающей среды (EMRL) из Bri-gham Young University, USA, в кооперации с Waterways Experiment Station американского Корпуса военных инженеров (US ACE-WES) и федеральным Управлением автомагистралей (FHWA). Пакет предназначен для моделирования и анализа движения воды, имеющей свободную поверхность. Имеет, как и предыдущий пакет, полностью графический, интегрированный в операционную систему Windows, интерфейс. Включает восемь основных модулей:
RMA2, HIVEL2D, FESWMS — расчет по двумерным моделям прибрежной циркуляции;
НМЛ 10. CH3D — расчет по трехмерным моделям прибрежной циркуляции;
ADCIRC, M2D, BOUSS2D — расчет по моделям океанической циркуляции;
CGWAVE — интегрирование уравнений рефракции-дифракции и BOOSS2D — интегрирование уравнений Буссинеска; STWAVE — спектральная модель расчета ветровых волн;
RMA4 и SED2D-WES — модели транспорта наносов;
HEC-RAS — одномерная прибрежная модель;
РТМ — модель, основанная на лагранжевом подходе.
Модуль RMA2 доставляет конечно-элементное решение уравнений Навье—Стокса для турбулентных потоков с осреднением по числу Рейнольдса. Трение определяется из уравнений Маннинга или Шези. Для расчета турбулентных характеристик используются коэффициенты турбулентной вязкости. Для решения уравнений используется процедура Галёркина. Модуль SED2D-WES позволяет находить решения как стационарных, так и нестационарных задач. Основное уравнение — уравнение конвекции— диффузии.
DELFT-3D. Полное название пакета — Интегрированная база трехмерного моделирования течений, транспорта наносов, волн, качества воды, морфологических и экологических изменений в прибрежных, речных, озерных областях и эстуариях. Разработчик — лаборатория Delft Hydraulics, Нидерланды. Основные модули, используемые при изучении динамики наносов: FLOW, SED, MOR. Так, FLOW — двумерная и трехмерная программы гидродинамического и транспортного моделирования нестационарного движения жидкости и переноса на криволинейных сетках. Гидродинамические характеристики, вычисленные в этом модуле, используются в качестве входных параметров в других модулях. Связь между модулями осуществляется через графический интерфейс пользователя. SED и MOR — модули для вычисления движения наносов, как связанных (глина), так и несвязанных (песок). Первый модуль используется при кратковременных расчетах, когда топографическими изменениями можно пренебречь. Второй — при долговременных расчетах, когда нужно учитывать эти изменения. Для расчета транспорта наносов применяется трехмерное уравнение конвекции-диффузии.
Заключение
Дана краткая характеристика теоретических моделей динамики наносов. Приведены некоторые результаты экспериментальных исследований, имеющих своей целью получить надежные оценки заносимости акватории заданного объекта и разработать рекомендации для снижения уровня заносимости до допустимых значений. Рассмотрено решение задач об изменении береговой линии при отложении наносов и об отложении наносов у мола. Приведены основные характеристики трех пакетов прикладных программ для расчета динамики наносов в прибрежной зоне. Приложен подробный список отечественных публикаций, посвященный динамике наносов.
Список литературы
[1] Айбулатов H.A. Исследование вдольберегового перемещения песчаных наносов в море. М.: Наука, 1966. 159 с.
[2] Айбулатов H.A. Динамика твердого вещества в шельфовой зоне. Л.: Гидрометеоиздат,
1990. 271 с.
[3] Айбулатов H.A., Артюхов Ю.В. Геоэкология шельфа и берегов Мирового океана. СПб.: Гидрометеоиздат, 1993. 304 с.
[4] Айбулатов H.A. и др. Геоэкология шельфа и берегов морей России. М.: Ноосфера, 2001. 428 с.
[5] Айбулатов H.A. Деятельность России в прибрежной зоне моря и проблемы экологии. М.: Наука, 2005. 364 с.
[6] Алешков Ю.З. Трансформация поверхностных гравитационных волн в жидкости при медленно изменяющейся глубине // Исследования по упругости и пластичности. 1980. № 13. С. 224-230.
[7] Алешков Ю.З. Теория волн на поверхности тяжелой жидкости. Л.: Изд-во ЛГУ, 1981. 196 с.
[8] Алешков Ю.З. Течения и волны в океане. СПб.: Изд-во Санкт-Петерб. ун-та, 1996. 228 с.
[9] Алешков Ю.З. Распространение примесей в океане // Вест. СПбГУ. Сер. 1. 2001. Вып. 1, № 1. С. 12-20.
[10] Алешков Ю.З. Изменение береговой линии // Петербургские фрагменты научной картины мира. СПб.: Изд-во РАН СПбНЦ, 2002. С. 19-27.
[11] Алешков Ю.З. Динамические процессы в прибрежной зоне моря// Вест. СПбГУ. Сер. 1. 2003. Вып. 1, № 1. С. 64-72.
[12] Анцыферов С.М., Косьян Р.Д. Взвешенные наносы в верхней части шельфа. М.: Наука, 1986. 223 с.
[13] Берега / П.А. Каплин, O.K. Леонтьев, С.А. Лукьянова, Л.Г. Никифоров. М.: Мысль,
1991. 479 с.
[14] Вечорек В.И. Об использовании мола в качестве наносозадерживающего сооружения // Гидромеханика. 1971. № 19. С. 35-39.
[15] Вечорек В.И. Переформирование волнами линии берега у поперечного наносоудержи-вающего сооружения // Гидротехническое строительство. 1981. № 1. С. 33-36.
[16] гришанин К.В. Динамика русловых потоков. Л.: Гидрометеоиздат, 1979. 311 с.
[17] Гришин H.H. Механика придонных наносов. М.: Наука, 1982. 160 с.
[18] Дебольский В.К., Зайдлер Р., Массель С. и др. Динамика русловых потоков и литодинамика прибрежной зоны моря. М.: Наука, 1994. 302 с.
[19] Динамические процессы береговой зоны моря / Под ред. Р. Д. Косьяна, И. С. Подымо-ва, Н.В. Пыхова. М.: Научный мир, 2003. 326 с.
[20] Жуковец A.M. О приближенном моделировании динамики песчаных наносов при волнении // Тр. координац. совещаний по гидротехнике. Л.: Энергия, 1971. Т. 66. С. 77-83.
[21] Зенкович В.П. Динамика и морфология морских берегов. Ч. 1. Волновые процессы. М.; Л.: Морской транспорт, 1946. 496 с.
[22] Зенкович В.П. Основы учения о развитии морских берегов. М.: Изд-во АН СССР, 1962. 710 с.
[23] Знаменская Н.С. Донные наносы и русловые процессы. Л.: Гидрометеоиздат, 1976. 191 с.
[24] Картвелишвили H.A. Потоки в недеформированных руслах. Л.: Гидрометеоиздат, 1973. 279 с.
[25] Косьян Р.Д., Пыхов Н.В. Гидрогенные перемещения осадков в береговой зоне моря. М.: Наука, 1991. 280 с.
[26] Кузнецов С.Ю., Сперанский Н.С. Расчет движения взвешенных частиц песка в монохроматической волне на мелководье // Океанология. 1985. Т. 25, № 5. С. 780-783.
[27] Леонтьев И.О. Динамика прибойной зоны. М.: Изд-во ПО АН СССР, 1989. 185 с.
[28] Леонтьев И.О. О возможности прогноза штормовых деформаций профиля подводного берегового склона // Океанология. 1991. Т. 31, № 5. С. 735-743.
[29] Леонтьев И.О. Расчет волн, течений и вдольбереговых потоков наносов в прибрежной зоне // Океанология. 1993. Т. 33, № 3. С. 429-434.
[30] Леонтьев И.О. Поперечный транспорт наносов и деформация профиля дна в прибрежной зоне // Океанология. 1994. Т. 34, № 2. С. 287-293.
[31] Леонтьев И.О. Пространственная модель прибрежной морфодинамики, обусловленной штормовым волнением // Океанология. 1995. Т. 35, № 2. С. 275-281.
[32] Leont'yev I.O. Numerical modelling of beach erosion during storm event // Coastal Engineering. 1996. Vol. 29. P. 187-200.
[33] Леонтьев И.О. О динамических изменениях профиля пляжа во время шторма // Океанология. 1997. Т. 37, № 1. С. 136-144.
[34] Leont'yev I.O. Modelling of morphological changes due to coastal structures // Coastal Engineering. 1999. Vol. 38. P. 143-166.
[35] Леонтьев O.K. Геоморфология морских берегов и дна. M.: Изд-во Моск. ун-та, 1955. 379 с.
[36] Леонтьев O.K. Основы геоморфологии морских берегов. М.: Изд-во Моск. ун-та, 1961. 418 с.
[37] Лонгинов В.В. Динамика береговой зоны бесприливных морей. М.: Изд-во АН СССР, 1963. 379 с.
[38] Лонгинов В.В. Очерки литодинамики океана. М.: Наука, 1973. 244 с.
[391 МАКАРОВ К.Н. Динамика берегов в зоне влияния продольных берегозащитных сооружений нетрадиционных конструкций // Гидротехническое строительство. 2002. № 6. С. 28-32.
[40] Пышкин Б.А., Максимчук B.JI., Цайтц Е.С. Исследование вдольберегового движения наносов на морях и водохранилищах. Киев: Наук, думка, 1967. 142 с.
[41] Пышкин Б.А., Цайтц Е.С., Сокольников Ю.М. Регулирование вдольберегового потока наносов. Киев: Наук, думка, 1972. 136 с.
[42] Пышкин Б.А. Динамика берегов водохранилищ. Киев: Наук, думка, 1973. 413 с.
[43] Сафьянов Г.А. Береговая зона океана в XX веке. М.: Мысль, 1978. 263 с.
[44] Сафьянов Г.А. Геоморфология морских берегов. М. 1996. 400 с.
[45] Селезов И.Т., Сидорчук В.Н., Яковлев В.В. Трансформация волн в прибрежной зоне шельфа. Киев: Наук, думка, 1983. 208 с.
[46] Сокольников Ю.Н., Вечорек В.И. Экспериментальное исследование процесса заполнения входящего угла мола // Тр. координац. совещаний по гидротехнике. Л.: Энергия, 1971. № 5. С. 114-117.
[47] Сокольников Ю.Н. Инженерная морфодинамика берегов и ее приложения. Киев: Наук, думка, 1976. 227 с.
Поступила в редакцию 14 марта 2008 г.
