Научный журнал
Международный журнал прикладных и фундаментальных исследований

ISSN 1996-3955
ИФ РИНЦ = 0,570

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

Шишеня А.В. 1 Кузнецова И.Ю. 2 Сухинов А.И. 3
1 ООО «CVisionLab»
2 Институт компьютерных технологий и информационной безопасности Южного федерального университета
3 Донской государственный технический университет
В статье рассматривается построение, исследование и программная реализация математической модели гидродинамики со свободной поверхностью и изменяемой геометрией береговой линии на основе гиперболизированных уравнений Навье – Стокса. Выполнен асимптотический анализ времени работы алгоритма моделирования гидродинамики на основе гиперболизированных уравнений и сравнение со временем работы алгоритма моделирования на основе классических уравнений (уравнений параболического типа) гидродинамики на основе неявных схем. Основываясь на этом анализе, получена оценка числа узлов, при котором гиперболизированная модель вычислительно более эффективна по сравнению с классической моделью. Выполнена параллельная реализация алгоритма расчета гиперболизированной модели гидродинамики, проведены численные эксперименты и получены значения ускорения и эффективности параллельного алгоритма. Показано, что за счет регулярной структуры хранения данных в памяти и малого числа условных переходов параллельный алгоритм достигает суперлинейного ускорения. В качестве тестового эксперимента для проверки корректности модели выполнено численное моделирование штормового нагона в Таганрогском заливе, произошедшего 24–25 сентября 2014 года. Полученная в результате моделирования зависимость уровня воды вблизи Таганрогского полуострова соответствует реально измеренным значениям.
гидродинамика
мелководный водоем
параллельное программирование
1. Ландау Л.Д., Лифшиц В.М. Гидродинамика. – М.: Наука, 1988. – 733 с.
2. Морозов Д.Н., Трапезникова М.А., Четверушкин Б.Н., Чурбанова Н.Г. Использование явных схем для моделирования процесса двухфазной фильтрации // Математическое моделирование. – 2011. – № 23 (7). – С. 52–60.
3. Д’Асчензо Н., Савельев В.И., Четверушкин Б.Н. Об одном алгоритме решения параболических и эллиптических уравнений // Журнал вычислительной математики и математической физики. – 2015. – № 55 (8). – С. 1320–1328.
4. Иванова К.Н., Четверушкин Б.Н., Чурбанова Н.Г. Квазигазодинамическая система уравнений и уравнения Навье – Стокса // Математическое моделирование. – 2004. – № 16 (4). – С. 98–104.
5. Шишеня А.В., Сухинов А.И., Кузнецова И.Ю. Моделирование штормового нагона в мелководном водоеме на основе прецизионных 3D моделей гидродинамики // Фундаментальные исследования. – 2017. – № 12–2. – С. 335–340.
6. Четверушкин Б.Н. Минимальные размеры в задачах механики сплошной среды // Математическое моделирование. – 2005. – № 17 (4). – С. 27–39.
7. Четверушкин Б.Н., Морозов Д.Н., Трапезникова М.А., Чурбанова Н.Г., Шильников Е.В. Об одной явной схеме для решения задач фильтрации // Математическое моделирование. – 2010. – № 22 (4). – С. 99–109.
8. Сухинов А.И., Шишеня А.В. Улучшение оценки параметра 1 попеременно-треугольного итерационного метода с априорной информацией // Известия ЮФУ. Технические науки. – 2010. – № 6 (107). – С. 7–15.
9. Шишеня А.В., Сухинов А.И. Построение регуляризированной математической модели гидродинамики и сгонно-нагонных явлений в мелководных водоемах и её параллельная реализация на параллельной вычислительной системе с общей памятью // Известия ЮФУ. Технические науки. – 2014. – № 12 (161). – С. 219–230.

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

При моделировании процессов, происходящих в мелководных водоемах, классическими являются модели, основанные на уравнениях Навье – Стокса (уравнениях параболического типа) [1]. Аппроксимация такой модели с помощью неявных схем обладает достаточным запасом устойчивости и налагает обременительные ограничения на шаг по времени как в случае явной схемы. В ряде работ [2–4] удалось ослабить это ограничение за счет гиперболизации уравнений модели. Кроме того, в случае большого числа узлов расчетной сетки такая гиперболизированная явная схема оказывается вычислительно более выгодной при параллельной программной реализации благодаря сокращению количества передач данных при решении системы линейных алгебраических уравнений (СЛАУ).

Формулировка непрерывной модели

В работе [5] была рассмотрена модель гидродинамики мелководных водоемов с уточненными граничными условиями. Недостатком описанной модели, в основе которой лежат уравнения параболического типа, является то, что в отличие от реально протекающего процесса, модель не способна описывать распространение значений физических величин (таких как скорость движения среды, давление и др.) с конечной скоростью. Данный недостаток отсутствует в уравнениях гиперболического типа, при этом формируется фронт волны, скорость движения которого зависит от коэффициента при второй производной по времени, а сама вторая производная принимается в качестве регуляризирующего слагаемого, сглаживающего решение и позволяющего ослабить ограничения на шаг по времени [6].

Распространение возмущений в водоеме возможно благодаря сжимаемости жидкости (зависимость плотности от давления), поэтому полагаем

hih01.wmf, (1)

где β – коэффициент сжимаемости жидкости (для воды hih02.wmf.

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

hih03.wmf, (2)

hih04.wmf, (3)

hih05.wmf (4)

где τ – регуляризирующий параметр, определяющий время между взаимодействиями, по порядку равное времени, за которое акустическая волна проходит ячейку расчетной сетки [7]:

hih06.wmf, (5)

где c – скорость звука в воде, hih07.wmf – безразмерный параметр.

Граничное условие для давления на свободной границе для уравнений (2)–(4) имеет вид

hih08.wmf

hih09.wmf. (6)

Дискретизация

При расщеплении по физическим процессам система (2)–(4) примет вид

hih10.wmf,

hih11.wmf,

hih12.wmf.

Граничное условие (6) примет вид

hih13.wmf.

Аппроксимация по пространственным направлениям выполнена с помощью интегро-интерполяционного метода с частичной заполненностью [6]. Условия устойчивости полученной разностной схемы:

hih14.wmf, hih15.wmf, hih16.wmf,

hih17.wmf, hih18.wmf,

hih19.wmf, hih20.wmf.

Оценим величину шага по времени, допустимого при решении практических задач. Примем hih21.wmf м, hih22.wmf, тогда hih23.wmf.

Аналитическое сравнение вычислительной сложности алгоритмов

Выполним оценку количества полных проходов по узлам расчетной сетки, необходимого для получения решения в конкретный момент времени для алгоритмов на основе классических уравнений Навье – Стокса и системы гиперболизированных уравнений гидродинамики. Пусть размер сетки N1× N2×N3, размер расчетной области l1×l2×l3, u, v, w – три компоненты поля скорости.

Предположим, что решение необходимо получить в момент времени hih24.wmf. В случае неявной схемы это означает, что нужно выполнить только один шаг по времени и решить уравнения типа диффузии-конвекции-реакции. При использовании попеременно-треугольного метода (ПТМ) количество внутренних шагов по времени асимптотически пропорционально максимальному числу узлов сетки по одному из направлений [8]. Учтем также, что в нашем случае эта константа пропорциональна отношению шага расчетной сетки по горизонтальному направлению к шагу по вертикальному направлению. Таким образом, количество шагов ПТМ совместно с методами сопряженных и бисопряженных направлений асимптотически равно [9] hih25.wmf.

При численном решении задач гидродинамики мелководных водоемов на основе трехмерной системы уравнений Навье – Стокса ограничение на шаг по времени определяется вертикальными составляющими скорости и шага по пространству. В случае применения модели на основе гиперболизированной системы уравнений гидродинамики, количество шагов по времени, требуемых для получения решения в момент времени T асимптотически равно hih26.wmf.

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

hih27.wmf. (7)

где hih28.wmf – скорость звука в воде.

В случае Азовского моря lx = 380000 м, lz = 17 м, характерная вертикальная скорость w = 0,01 м, тогда модель на основе гиперболизированных уравнений вычислительно более эффективна по сравнению с моделью на основе уравнений Навье – Стокса параболического типа при hih29.wmf. Такое соотношение количества узлов расчетной сетки вполне достижимо, то есть оба метода обладают примерно одинаковой применимостью. Далее будет показано, что на практике параллельная реализация метода на основе гиперболизированных уравнений позволяет получить значительно более эффективный алгоритм расчета модели.

Численные эксперименты

Проверку адекватности построенной модели на основе гиперболизированных уравнений реально протекающему физическому процессу, а также сравнение прогнозируемых гидродинамических параметров и быстродействие разработанных алгоритмов прогнозирования на основе классических (параболических) и гиперболизированных уравнений Навье – Стокса выполним на основе данных о затоплении прибрежных регионов Азовского моря из-за штормового нагона, которое произошло 24–25 сентября 2014 г. в районе порта г. Таганрога. Скорость и направление ветра приведены в таблице [9].

Данные о скорости, направлении ветра и уровне воды в районе порта г. Таганрога, 24–25 сентября 2014 г.

Момент времени

24 сентября 2014

25 сентября 2014

00:00

04:00

08:00

12:00

16:00

20:00

00:00

04:00

08:00

Направление ветра

Ю-В

Ю-В

Ю

Ю-Ю-З

Ю-Ю-З

З

Ю-З

Ю-З

Ю-З

Скорость ветра, м/с

17

18

27

33

37

33

28

22

18

Уровень воды, см

12

24

66

160

380

420

390

320

280

Результаты численного моделирования представлены на рис. 1. В качестве начального приближения использовались результаты расчета течений в Азовском море в условиях отсутствия ветра. При численном моделировании нагонных явлений предполагалось, что возрастание скорости ветра 24 сентября 2014 г. от 0 м/с до 17 м/с (первый столбец таблицы) происходило линейным образом в течение 5 часов. Из графиков на рис. 1 видно, что результаты расчетов на грубой сетке достаточно точно повторяют результаты, полученные на подробной сетке, а значения функции возвышения уровня располагаются достаточно близко по отношению к наблюдаемым значениям.

hih1.tif

Рис. 1. Функция возвышения уровня в районе порта г. Таганрога (пунктирная линия – значения, рассчитанные на сетке с шагами по горизонтальным направлениям – 1003 м и 1013 м, тонкая линия – значения, рассчитанные на сетке с шагами по горизонтальным направлениям – 503 м и 510 м; шаг сетки по вертикальному направлению – 0,1 м; ломаная линия – натурные данные)

Также выполнена программная реализация математической модели гидродинамики и сгонно-нагонных явлений в мелководных водоемах на основе регуляризированных уравнений Навье – Стокса на параллельной вычислительной системе с общей памятью с помощью библиотеки OpenMP [9]. Зависимость ускорения от количества процессоров представлена на рис. 2, а эффективности – на рис. 3.

hih2.tif

Рис. 2. Ускорение параллельной программы

hih3.tif

Рис. 3. Эффективность параллельной программы

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

1. Явная схема расчета позволяет значительно сократить число синхронизаций параллельных потоков;

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

3. Применение интегро-интерполяционного метода с частичной заполненностью позволило использовать структурированные регулярные сетки и регулярные структуры данных для их хранения, что повышает эффективность работы КЭШ-памяти процессора.

Выводы

В работе рассмотрена математическая модель гидродинамики со свободной поверхностью и изменяемой геометрией береговой линии на основе гиперболизированных уравнений Навье – Стокса. Введение регуляризирующего параметра τ позволяет увеличить шаг по времени при численном решении уравнений с помощью явных схем и повысить эффективность параллельной программной реализации алгоритма по сравнению с моделью на основе классических уравнений (уравнений параболического типа) Навье – Стокса.

Проведено аналитическое сравнение вычислительной сложности алгоритмов на основе классических уравнений Навье – Стокса и системы гиперболизированных уравнений. Показано, что при выполнении условия (7) гиперболизированная модель гидродинамики вычислительно предпочтительнее классической модели на основе уравнений Навье – Стокса параболического типа. Другими словами, классическая модель предпочтительнее при малой скорости течений и грубой расчетной сетке, а при увеличении количества узлов сетки (при измельчении сетки) или при возрастании скорости течений лучше использовать регуляризированную модель.

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

Работа выполнена при финансовой поддержке РФФИ по проекту № 15-07-08408.


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

Шишеня А.В., Кузнецова И.Ю., Сухинов А.И. РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ ЗАТОПЛЕНИЯ ПРИБРЕЖНЫХ РАЙОНОВ НА ОСНОВЕ ГИПЕРБОЛИЗИРОВАННЫХ УРАВНЕНИЙ ГИДРОДИНАМИКИ // Международный журнал прикладных и фундаментальных исследований. – 2018. – № 5-2. – С. 306-311;
URL: https://applied-research.ru/ru/article/view?id=12260 (дата обращения: 15.09.2019).

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

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