Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

MATHEMATICAL MODELING OF WAVE PROCESSES AND TRANSPORT OF BOTTOM MATERIALS TAKING INTO ACCOUNT COASTAL STRUCTURES IN COASTAL WATER AREAS

Chistyakov A.E. 1 Protsenko E.A. 2 Timofeeva E.F. 3 Nikolich N.A. 4
1 Don State Technical University
2 Taganrog Institute named after A.P. Chekhov Rostov State Economic University
3 North Caucasus Federal University
4 Russian State University of Oil and Gas (National Research University) named after I.M. Gubkin
The necessity of numerical modeling of the mechanisms of formation of the field of suspension and sediment flow determines the urgency of the problem under study. Under modern conditions, the process of developing shallow coastal zones, with constructive transformation and changes in coastal reliefs, is becoming increasingly important. Assessment of the hydrodynamic impact on coastal infrastructure facilities located in the coastal waters is an important task at present. To determine the dynamics of the phenomena occurring in the coastal water areas and to predict possible interference in the ecosystem, the need to build mathematical models of the transport processes of matter in the coastal waters taking into account the presence of coastal structures is actualized. In accordance with the goal, it is necessary to analyze existing models of sediment transport in aquatic ecosystems and develop a set of interrelated models of three-dimensional wave processes and sediment transport designed to simulate hydrodynamic wave processes in the presence of various shore structures.
sediment transport
waves
shallow water
flooding
surface structures
and power action

Последние десятилетия характеризуются интенсивным использованием прибрежных районов морей и океанов в экономической и хозяйственной сферах деятельности страны. Изучение гидродинамических процессов прибрежных акваторий связано с исследованием влияния волновых процессов, генерируемых в открытом море или в самой прибрежной зоне водоема. Движение волн может привести к негативным результатам, оказывающим влияние на эксплуатацию как самой прибрежной зоны, так и непосредственно объектов прибрежной инфраструктуры: к трансформации донной поверхности, происходящей вследствие подъема донных отложений, к абразии, то есть процессу разрушения волнами и прибоем берегов различных водных систем.

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

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

Характерной особенностью прибрежных акваторий является значительное влияние донной поверхности на волновые процессы, что затрудняет исследование приливных явлений в прибрежных районах морей и устьях рек. Влияние волновых процессов на прибрежную зону может иметь амбивалентный характер: волновые процессы могут оказывать существенное влияние на аккумуляцию и абразию береговой зоны водоема и непосредственно на прибрежные конструкции [2].

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

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

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

Математическая модель транспорта донного материала

Исходными уравнениями математической модели транспорта донных веществ являются уравнения, описывающие переформирование прибрежного рельефа [3]:

chist01.wmf

chist02.wmf

chist03.wmf (1)

где Н – глубина дна, отсчитываемая от невозмущенной поверхности водоема; ε – пористость грунта; chist04.wmf – расход наносов, chist05.wmf; x, y – горизонтальные декартовы координаты; τb – касательное напряжение на дне; τbc – критическое значение касательного напряжения, при котором начинается перемещение наносов, g – ускорение свободного падения, ρ – плотность жидкости, А и β – безразмерные постоянные (А равно 19,5, β равно 3), ω – частота волны, d – характеристика осадков.

Принимая во внимание ограничения на касательные напряжения на дне расчетной области, система уравнений (1) может быть представлена в виде [4]:

chist06.wmf, (2)

chist07.wmf, (3)

где chist08.wmf – функция Хэвисайда.

Уравнение (2) дополняется начальным и граничным условиями:

chist09.wmf (4)

chist10.wmf (5)

chist1.tif

Рис. 1. Начальный рельеф дна

chist2.tif

Рис. 2. Геометрия дна через 80 минут

Для расчета тангенциальных напряжений в модели транспорта наносов необходима информация о поле вектора скорости на дне водоема [5].

Рис. 1–2 иллюстрируют результаты численных экспериментов моделирования динамики донной поверхности. При этом размер расчетной сетки равен 200×200, шаг по пространственным переменным 0,1 м, шаг по времени 0,01 с, скорость движения ветра 5 м/с и направлена слева направо.

Математическая модель волновой гидродинамики

Для математического описания задачи волновой динамики жидкости использована модель [6], представленная на рис. 3.

chist11.wmf, (6)

chist12.wmf, (7)

chist13.wmf. (8)

chist14.wmf, chist15.wmf,

chist16.wmf, chist17.wmf, (9)

chist18.wmf, chist19.wmf,

chist20.wmf, (10)

chist21.wmf, chist22.wmf,

chist23.wmf, (11)

где chist24.wmf – вектор скорости движения водной среды; Р – давление; μ, η – коэффициенты турбулентного обмена по горизонтальному и вертикальному направлениям соответственно; g – ускорение свободного падения; ρ – плотность жидкости; τx, τy – составляющие тангенциального напряжения на дне жидкости; П – поток вектора скорости через боковую поверхность; L – расстояние от поверхности жидкости до дна (глубина жидкости с учетом возвышения уровня) на боковой границе. Система координат выбрана таким образом, что ось Ox совмещена с поверхностью невозмущенной жидкости и направлена в сторону берега, ось Oy – вертикально вниз.

Для численной реализации дискретной модели задачи волновой гидродинамики введена сетка:

chist25.wmf,

где τ – шаг по времени, hx, hy – шаги по пространству, Nt – верхняя граница по времени, Nx, Ny – границы по пространству.

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

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

Для разработки математической модели движения водной среды использована трехмерная гидростатическая модель, представленная на рис. 5.

chist3.tif

Рис. 3. Исходные уравнения математической модели волновой гидродинамики

chist4a.tif chist4b.tif

Рис. 4. Динамика профиля уровня возвышения жидкости при наличии прибрежной конструкции

chist5.tif

Рис. 5. Исходные уравнения математической модели волновых процессов, учитывающей наличие надводных береговых сооружений

chist26.wmf,

chist27.wmf, (12)

chist28.wmf, (13)

chist29.wmf. (14)

chist31.wmf, chist32.wmf (15)

chist32.wmf chist33.wmf, (16)

chist34.wmf, (17)

где ξ – функция подъема уровня (функция возвышения), chist35.wmf – вектор скорости движения водной среды, Р – давление,

μ, η – коэффициенты турбулентного обмена по горизонтальному и вертикальному направлениям соответственно, g – ускорение свободного падения, ρ – плотность жидкости, τx, τy – тангенциальное напряжение на дне жидкости, оси Ox, Oy имеют горизонтальное направление.

chist6.tif

Рис. 6. Функция возвышения уровня при обтекании водной средой прибрежного объекта

Для расчета давления в гидростатическом приближении при наличии надводного объекта, установленного на дне водоема, используем уравнение

chist36.wmf

Функция возвышения в этом случае примет вид

chist37.wmf

где χ – функция, описывающая геометрию дна.

Взаимосвязанный комплекс программ был использован для расчета гидродинамического воздействия волновых процессов на прибрежный объект. Размеры подводной части конструкции: 5 м×10 м×20 см. Объект установлен на шести опорах. Зона моделирования имеет размеры 50 м×50 м×1 м (рис. 6).

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

Заключение

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

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

Работа выполнена при поддержке РНФ (проект № 17–11–01286).