Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

МОДЕЛИРОВАНИЕ СИЛОВОГО ГИДРОДИНАМИЧЕСКОГО ВОЗДЕЙСТВИЯ ВОЛН НА ОПОРЫ НАДВОДНЫХ КОНСТРУКЦИЙ

Сухинов А.И. 1 Никитина А.В. 2 Фоменко Н.А. 3 Тимофеева Е.Ф. 4 Проценко С.В. 5
1 Донской государственный технический университет
2 Научно-исследовательский институт многопроцессорных вычислительных систем им. А.В. Каляева Южный федеральный университет
3 Московский государственный технический университет им. Н.Э. Баумана
4 Северо-Кавказский федеральный университет
5 Таганрогский педагогический институт имени А.П. Чехова Ростовский государственный экономический университет (РИНХ)
Работа посвящена разработке математической модели волновых процессов, учитывающей гидродинамическое воздействие на различные надводные береговые сооружения при наличии поверхностных волн в мелководных водоемах. Проведен натурный эксперимент по измерению различных параметров распространения волны на мелководье. Разработана двумерная модель волновых гидродинамических процессов, описывающая поведение водной среды как в случае наличия надводной конструкции, установленной на дне водоема, так и в случае его отсутствия. На основе построенных алгоритмов разработан комплекс программ, предназначенный для моделирования распространения волновых гидродинамических процессов. Комплекс позволяет задавать форму и интенсивность источника колебаний, геометрию надводного объекта, а также учитывает затопления и осушения прибрежных участков.
волны
мелкая вода
затопления
надводные сооружения
силовое воздействие
1. Белоцерковский О.М., Гущин В.А., Щенников В.В. Метод расщепления в применении к решению задач динамики вязкой несжимаемой жидкости // Журнал вычислительной математики и математической физики. – 1975. –№ 15:1 – С. 197–207.
2. Васильев В.С., Сухинов А.И. Прецизионные двумерные модели мелких водоемов // Математическое моделирование. – 2003. – Т. 15, № 10. – С. 17–34.
3. Гущин В.А., Матюшин П.В. Математическое моделирование и визуализация трансформации вихревой структуры течения около сферы при увеличении степени стратификации жидкости // Журнал вычислительной математики и математической физики. – 2011. – Т. 51, № 2. – С. 268–281.
4. Коновалов А.Н. К теории попеременно-треугольного итерационного метода // Сибирский математический журнал. – 2002. – 43:3. – С. 552–572.
5. Самарский А.А. Теория разностных схем. – М.: Наука, 1989.
6. Самарский А.А. Методы решения сеточных уравнений / А.А. Самарский, Е.С. Николаев. – М.: Наука, 1978. – 592 с.
7. Сухинов А.И., Чистяков А.Е. Адаптивный модифицированный попеременно-треугольный итерационный метод для решения сеточных уравнений с несамосопряженным оператором // Математическое моделирование. – 2012. – Т. 24, № 1, – С. 3–20.
8. Сухинов А.И., Чистяков А.Е. Параллельная реализация трехмерной модели гидродинамики мелководных водоемов на супервычислительной системе // Вычислительные методы и программирование: Новые вычислительные технологии. – 2012. – Т.13. – С. 290–297.
9. Сухинов А.И., Чистяков А.Е., Алексеенко Е.В. Численная реализация трехмерной модели гидродинамики для мелководных водоемов на супервычислительной системе // Математическое моделирование. – 2011. – Т. 23, № 3. – С. 3–21.
10. Сухинов А.И., Чистяков А.Е., Проценко Е.А. Математическое моделирование транспорта наносов в прибрежных водных системах на многопроцессорной вычислительной системе // Вычислительные методы и программирование. – 2014. – Т. 15. – С. 610–620.
11. Сухинов А.И., Чистяков А.Е., Тимофеева Е.Ф., Шишеня А.В. Математическая модель расчета прибрежных волновых процессов // Математическое моделирование. – 2012. – Т. 24, № 8. – С. 32–44.

В результате сложного взаимодействия волн с рельефом дна, береговой линией и различными объектами наблюдается их рефракция, дифракция и изменения структуры береговой зоны. Возникает необходимость в построении математических моделей, учитывающих особенности рельефа дна прибрежной акватории и наличие надводных сооружений. Не приуменьшая целесообразности использования трехмерных моделей волновых гидродинамических процессов [8–11], следует подчеркнуть важность разработки двумерных моделей гидродинамики при исследовании воздействия волн на надводные объекты, на основе данных моделей при приемлемой точности можно производить оперативные численные расчеты, при условии ограниченности вычислительных ресурсов. Полученные численные результаты распространения гидродинамических колебаний позволят провести оценку воздействия волн на объекты, находящиеся в прибрежной зоне мелководных водоемов.

Постановка задачи волновой гидродинамики

Для построения двумерной математической модели движения водной среды использовалась трехмерная гидростатическая модель. Исходными уравнениями гидродинамики являются:

– система уравнений Навье – Стокса

Sukhinov01.wmf

Sukhinov02.wmf (1)

– уравнение неразрывности для несжимаемой жидкости

Sukhinov03.wmf (2)

– уравнение гидростатики

Sukhinov04.wmf (3)

Система уравнений (1)–(3) рассматривается при следующих граничных условиях: на дне условие непроницаемости и трения [2–5]

Sukhinov05.wmf Sukhinov06.wmf Vn = 0,

на поверхности задается подъем уровня и ветровые напряжения

Sukhinov07.wmf Sukhinov08.wmf Sukhinov09.wmf

на боковых границах условие скольжения без трения

Sukhinov10.wmf Sukhinov11.wmf Sukhinov12.wmf

где ξ – функция подъема уровня (функция возвышения); V = {u, v, w} – вектор скорости движения водной среды; P – давление; μ, η – коэффициенты турбулентного обмена по горизонтальному и вертикальному направлениям соответственно; g – ускорение свободного падения; ρ – плотность жидкости; τx, τy – тангенциальное напряжение на дне жидкости.

Уравнение расчета давления в гидростатическом приближении при наличии надводной конструкции, установленной на дне водоема:

Sukhinov13.wmf

Доопределим функцию возвышения в случае наличия на поверхности надводного тела:

Sukhinov14.wmf

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

В уравнение неразрывности (2) в гидростатическом случае запишем в виде

Sukhinov15.wmf (4)

где θ = min(χ, ξ); H = h + θ; h – глубина водоема.

Система уравнений (1) с учетом гидростатического приближения в двумерном случае имеет вид

Sukhinov16.wmf

Sukhinov17.wmf (5)

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

Дискретная модель гидродинамики

Область моделирования вписана в прямоугольник и покрыта равномерной прямоугольной расчетной сеткой ω = ωt×ωx×ωy:

Sukhinov18.wmf

Sukhinov19.wmf

Sukhinov20.wmf

где n, i, j – индексы по временной координате и пространственным координатным направлениям Ox, Oy соответственно; ht, hx, hy – шаги по временной координате и пространственным координатным направлениям Ox, Oy соответственно; Nt, Nx, Ny – количество узлов по временной координате и пространственным координатным направлениям Ox, Oy соответственно; lt, lx, ly – длина расчетной области по временной координате и пространственным координатным направлениям Ox, Oy соответственно.

Расчетные ячейки представляют собой прямоугольники, они могут быть заполненными, частично заполненными или пустыми [11]. Поле скоростей и давление рассчитываются в вершинах ячейки. Центры ячеек и узлы разнесены на hx/2 и hy/2 по координатам x и y соответственно. Обозначим через oi,j «заполненность» ячейки (i, j). Вершинами ячейки (i, j) являются узлы (i, j), (i – 1, j), (i, j – 1), (i – 1, j – 1). Вводятся коэффициенты q0, q1, q2, q3, q4, описывающие «заполненность» областей, находящихся в окрестности ячейки. Коэффициенты qm можно вычислить по формулам

Sukhinov21.wmf

Sukhinov22.wmf

Sukhinov23.wmf Sukhinov24.wmf

Sukhinov25.wmf

Sukhinov26.wmf (6)

Рассмотрим двумерную модель движения водной среды, представленную уравнениями (4), (5). Воспользуемся схемами расщепления по физическим процессам [1–4] для системы (5). При этом решения задачи находим вначале на некоем промежуточном временном слое:

Sukhinov27.wmf

Sukhinov28.wmf (7)

затем на следующем временном слое

Sukhinov29.wmf Sukhinov30.wmf (8)

Для решения задачи (8) вычислим функцию возвышения уровня, дифференцируя первое уравнение системы по переменной x, второе по переменной y, сложив их, имеем

Sukhinov31.wmf

Полученное выражение с учетом уравнения неразрывности примет вид

Sukhinov32.wmf (9)

Дискретные аналоги операторов конвективного Sukhinov33.wmf и диффузионного Sukhinov34.wmf переноса, учитывающие частичную заполненность ячеек, в случае граничных условий третьего рода Sukhinov35.wmf запишем в следующем виде:

Sukhinov36.wmf

Sukhinov37.wmf

Полученные сеточные уравнения решены на основе адаптивного модифицированного попеременно-треугольного итерационного метода вариационного типа [4, 6, 7].

Измерение параметров волновых процессов на основе натурных наблюдений

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

Sukhinov38.wmf

где [xn] – заданная последовательность чисел с постоянным шагом дискретизации размерности N.

С учетом обозначений: a = Rex, b = Imx, тригонометрическая функция, проходящая через точки {n, xn}, примет вид

Sukhinov39.wmf

На рис. 1 точками отмечены значения спектра волны, линиями представлены функции, распределенные по нормальному и логнормальному законам и имеющие математические ожидания и дисперсии, соответствующие реальным натурным данным.

На рис. 2 красным представлена функция, распределенная по логнормальному закону, синим – функция, распределенная по нормальному закону, черным – реальные натурные значения функции возвышения уровня.

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

Волновые процессы по большей части можно описать тремя величинами: математическим ожиданием (период волны), дисперсией и амплитудой нормального или логнормального распределения компонент спектра. Эти величины получены при обработке данных натурного эксперимента и используются в качестве граничных условий для математических моделей волновых гидродинамических процессов.

pic_51.tif

Рис. 1. Спектр функции возвышения уровня: 1 – нормальное распределение; 2 – логнормальное распределение

pic_52.wmf

Рис. 2. Сопоставление результатов измерений с нормальным и логнормальным распределением

№ п/п

Глубина, см

Период волны, с

Средняя высота волны, см

Максимальное значение высоты волны, см

Дисперсия функции возвышения уровня

Корреляция с нормальным распределением

Корреляция с логнормальным распределением

1

12,734

3,181

1,434

3,266

3,384

0,67622403

0,72818161

2

21,657

3,187

2,216

5,127

2,875

0,71970734

0,75497854

3

34,296

3,257

2,673

6,673

2,587

0,76756352

0,80809736

4

47,696

3,208

2,903

7,278

2,373

0,80434285

0,81516631

5

50,221

3,238

3,408

8,779

2,465

0,80072646

0,82234947

6

56,95

3,323

3,423

10,05

2,539

0,82520735

0,83499856

7

58,256

3,094

3,538

13,742

2,468

0,70451786

0,75010325

8

75,284

3,482

3,595

12,716

2,317

0,80464887

0,82816629

9

83,353

3,056

4,472

14,647

2,498

0,7677805

0,80442466

10

123,251

3,23

4,671

15,749

2,327

0,78716382

0,82809779

 

Результаты численных экспериментов

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

В качестве примера практического использования проблемно-ориентированного комплекса программ решается задача расчета гидродинамического воздействия волн опоры на сооружения. Размеры надводной конструкции: ширина 5 м, длина 10 м, глубина погружения 20 см. Конструкция установлена на дне водоема при помощи шести опор. Выделенный участок моделирования имеет размеры 50 на 50 м и глубину 1 м. Источник возмущений задается на некотором расстоянии от надводного объекта. В начальный момент времени жидкость находится в состоянии покоя. Требуется определить последующее движение водной среды при наличии на поверхности надводного объекта и силовые гидродинамические нагрузки на опоры сооружения. Для решения поставленной задачи использована сетка размерами 100×100, шаг по времени равен 0,01 с.

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

pic_53.tif

Рис. 3. Движение водной среды

pic_54.tif

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

pic_55.tif

Рис. 5. Результаты расчета по распространению волновых колебаний

Заключение

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

Работа выполнена при финансовой поддержке РФФИ по проектам № 16-37-00129, № 15-07-08626, № 15-01-08619.


Библиографическая ссылка

Сухинов А.И., Никитина А.В., Фоменко Н.А., Тимофеева Е.Ф., Проценко С.В. МОДЕЛИРОВАНИЕ СИЛОВОГО ГИДРОДИНАМИЧЕСКОГО ВОЗДЕЙСТВИЯ ВОЛН НА ОПОРЫ НАДВОДНЫХ КОНСТРУКЦИЙ // Фундаментальные исследования. – 2016. – № 12-4. – С. 777-783;
URL: https://fundamental-research.ru/ru/article/view?id=41168 (дата обращения: 29.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674