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

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

Чистяков А.Е. 1 Хачунц Д.С. 1 Тимофеева Е.Ф. 2 Фоменко Н.А. 3 Moscow State Technical University named after N.E. Bauman, Moscow 4
1 Научно-исследовательский институт многопроцессорных вычислительных систем имени академика А.В. Каляева, ГОУ ВО «Южный федеральный университет»
2 ФГАОУ ВПО «Северо-Кавказский федеральный университет»
3 Московский государственный технический университет имени Н.Э. Баумана, Москва
4 Инженерно-технологическая академия, ГОУ ВО «Южный федеральный университет»
Рассматривается непрерывная математическая модель расчета выхода волны на берег и ее разрушения, учитывающая такие физические параметры, как турбулентный обмен, сложная геометрия дна и береговой линии, трение о дно. Для этой модели построен эффективный алгоритм ее реализации, учитывающий динамическое изменение уровня возвышения жидкости. Разработан метод построения сеток с динамически изменяющейся геометрией расчетной области. Построен алгоритм решения задачи диффузии-конвекции на основе явных регуляризованных схем и выполнена его параллельная реализация. Разработан программный комплекс, предназначенный для визуализации двумерных полей скоростей движения водной среды в случае математического моделирования наката и обрушения волны на берег, а также для прогнозирования волновых процессов в прибрежных акваториях, на основе которого получены поля течений водной среды, согласующиеся с реальным физическим процессом.
гидродинамика
математическая модель
волновые гидродинамические процессы
система уравнений навье – стокса
рельеф дна
численный эксперимент
1. Васильев В.С., Сухинов А.И. Прецизионные двумерные модели мелких водоемов// Математическое моделирование. – 2003. – Т. 15, № 10. – С. 17.
2. Самарский А.А. Теория разностных схем. – М.: Наука, 1989. – 616 с.
3. Самарский А.А., Николаев Е.С. Методы решения сеточных уравнений. – М.: Наука, 1978. – 532 с.
4. Сухинов А.И., Никитина А.В., Чистяков А.Е. Моделирование сценария биологической реабилитации Азовского моря// Математическое моделирование. – 2012. – Т. 24, № 9. – С. 3–21.
5. Сухинов А.И., Никитина А.В., Чистяков А.Е., Семенов И.С. Математическое моделирование условий формирования заморов в мелководных водоемах на многопроцессорной вычислительной системе // Вычислительные методы и программирование: новые вычислительные технологии. – 2013. – Т. 14, № 1. – С. 103–112.
6. Сухинов А.И., Хачунц Д.С., Чистяков А.Е. Математическая модель распространения примеси в приземном слое атмосферы прибрежной зоны и ее программная реализация // Журнал вычислительной математики и математической физики. – 2015. – Т. 55, № 7. – С. 1238.
7. Сухинов А.И., Хачунц Д.С., Чистяков А.Е. Параллельные алгоритмы для прогноза состояния воздушной среды рекреационных зон прибрежных систем Вестник компьютерных и информационных технологий. – 2015. – № 5 (131). – С. 55–60.
8. Сухинов А.И., Чистяков А.Е., Проценко Е.А. Математическое моделирование транспорта наносов в прибрежной зоне мелководных водоемов // Математическое моделирование. – 2013. – Т. 25, № 12. – С. 65–82.
9. Сухинов А.И., Чистяков А.Е., Тимофеева Е.Ф., Шишеня А.В. Математическая модель расчета прибрежных волновых процессов // Математическое моделирование. – 2012. – Т. 24, № 8. – С. 32–44.
10. Сухинов А.И., Чистяков А.Е., Фоменко Н.А. Методика построения разностных схем для задачи диффузии-конвекции-реакции, учитывающих степень заполненности контрольных ячеек // Известия Южного федерального университета. Технические науки. – 2013. – № 4 (141). – С. 87–98.
11. Четверушкин Б.Н. Пределы детализации и формулировка моделей уравнений сплошных сред // Математическое моделирование. – 2012. – Т. 24, № 11. – С. 33–52.

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

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

– уравнения Навье – Стокса:

chistyak01.wmf (1)

chistyak02.wmf (2)

– уравнение неразрывности:

chistyak03.wmf (3)

Уравнение (3) в случае несжимаемой жидкости примет вид

chistyak04.wmf (4)

Уравнения (1)–(4) рассматриваются при следующих граничных условиях:

– на дне области:

chistyak05.wmf chistyak06.wmf

chistyak07.wmf (5)

– на поверхности жидкости:

chistyak08.wmf chistyak09.wmf

chistyak10.wmf chistyak11.wmf (6)

– на боковой границе:

chistyak12.wmf chistyak13.wmf

chistyak14.wmf (7)

где V = {u, v} – вектор скорости движения водной среды; P – давление; ?, ? – коэффициенты турбулентного обмена по горизонтальному и вертикальному направлениям соответственно; g – ускорение свободного падения; ? – плотность жидкости; ?x, ?y – составляющие тангенциального напряжения на дне жидкости; ? – поток вектора скорости через боковую поверхность; L – расстояние от поверхности жидкости до дна (глубина жидкости с учетом возвышения уровня) на боковой границе.

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

Расчетная область представляет собой прямоугольник. Для численной реализации дискретной математической модели поставленной задачи волновой гидродинамики вводится равномерная сетка [2, 3]:

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

Применяется непосредственная аппроксимация. В результате аппроксимации уравнений (1), (2) по временной переменной, вводя промежуточный слой n + ? согласно MAC-методу и расщепляя уравнения по физическим процессам, уравнения приводятся к следующему виду [1, 5, 9]:

chistyak16.wmf

chistyak17.wmf (8)

chistyak18.wmf chistyak19.wmf (9)

Продифференцировав уравнения (9) и преобразовав с учетом уравнения неразрывности (4), получим

chistyak20.wmf (10)

Расчет задач гидродинамики по данному методу «поправок к давлению» осуществляется в три этапа. На первом этапе на основе уравнений (8) считается поле скоростей. На втором этапе рассчитывается давление по уравнению (10). На третьем этапе из соотношений (9) уточняется поле скоростей по давлению.

Программная реализация

Для демонстрации визуализации процесса выхода и обрушения волны на берег с учетом особенностей рельефа дна прибрежной акватории и суши было разработано экспериментальное программное обеспечение «Waves» на базе ЭВМ. Программа «Waves» предназначена для построения двумерных полей скоростей движения водной среды в случае математического моделирования наката и обрушения волны на берег, а также для прогнозирования возможных ситуаций, связанных с волновыми процессами в прибрежных акваториях, с целью предопределения строительства сооружений и использования конкретного участка береговой линии. Разработанное экспериментальное программное обеспечение на базе ЭВМ предназначено для математического моделирования и демонстрации визуализации процесса выхода и обрушения волны на берег с учетом особенностей рельефа дна прибрежной акватории и суши.

Программный компонент «Waves» включает в себя следующие блоки: управляющий блок (в данном блоке содержится цикл по временной координате и вызываются функции: расчет поля скорости без учета давления, расчет давления, расчет двумерного поля скорости, расчет заполненности ячеек водной средой и функции ввода-вывода данных); блок ввода начальных данных для расчета течений и давления (задаются начальные распределения поля скорости и давления, а также маски граничных условий и начальные значения заполненности ячеек); блок построения сеточных уравнений для поля скорости без учета давления в соответствии с конечно-объемной схемой (считаются и записываются в массив коэффициенты и правая часть соответствующего сеточного уравнения, представленного в канонической форме); блок построения сеточных уравнений для поля давления; блок расчета заполненности ячеек [4, 10]; блок расчета поля скорости с учетом давления (результатом работы данного блока является расчет значений поля двумерного вектора скоростей на следующем временном слое); блок расчета пятидиагональных сеточных уравнений адаптивным попеременно-треугольным методом скорейшего спуска; блок вывода значений поля скоростей, давления и заполненности ячеек.

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

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

После разработки программы «Waves» были проведены тестирования ее на соответствие уже имеющимся результатам, которые были получены другими научными работами и опытным путем. Анализ показывает, что приближенная математическая модель выхода волны на берег и ее разрушение в прибрежной зоне водоема с определенной точностью согласуется с существующими данными.

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

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

pic_58.tif

Динамика изменения профиля уровня возвышения жидкости в случае наличия искусственного препятствия (волнореза)

chistyak15.wmf

Данные о времени счета, ускорении и эффективности

Количество ядер

Время, с

Ускорение

Эффективность

1

1215

1

100

2

623,1

1,95

97,5

4

319,7

3,8

95

8

164,2

7,4

92,5

16

85,6

14,2

88,8

32

46,2

26,3

82,2

64

25

48,6

75,9

Параллельная реализация задачи. Для параллельной реализации поставленной задачи будем использовать регуляризованные по Б.Н. Четверушкину явные схемы [11], при этом уравнение диффузии-конвекции [6–8]

chistyak21.wmf

с учетом слагаемых регуляризаторов запишется в виде [11]:

chistyak22.wmf

Следует отметить, что условием устойчивости для задачи диффузии-конвекции является chistyak23.wmf а в случае использования слагаемых регуляризаторов данное ограничение запишется в виде

chistyak24.wmf

Обратим внимание на то, что полученные с учетом ограничения минимальных масштабов дополнительные члены, по сути, не меняют решений, использующих классические постановки. Они выступают лишь в качестве физически обоснованных регуляризаторов, сглаживающих физически необоснованные эффекты, получающиеся при численном решении. В этой связи конкретное значение величин, входящих в виде коэффициентов в регуляризирующие члены, уже не играет особой роли. Важно, чтобы они находились в правильном диапазоне лишь по порядку величины. Этот факт существенно облегчает практическое использование указанного подхода. При параллельной реализации использованы методы декомпозиции сеточных областей для вычислительно трудоемких задач диффузии-конвекции, учитывающие архитектуру и параметры многопроцессорной вычислительной системы РОЦ НИТ ЮФУ. Пиковая производительность МВС составляет 18,8 TFlops. В качестве вычислительных узлов используются 128 однотипных 16-ядерных Blade-серверов HP ProLiant BL685c, каждый из которых оснащен четырьмя 4-ядерными процессорами AMD Opteron 8356 2,3 GHz и оперативной памятью в объеме 32 ГБ.

Результаты использования многопроцессорных технологий для решения уравнения диффузии-конвекции на основе явных регуляризованных схем приведены в таблице.

Заключение

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

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

Работа выполнена при частичной поддержке Задания №2014/174 в рамках базовой части государственного задания Минобрнауки России.


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

Чистяков А.Е., Хачунц Д.С., Тимофеева Е.Ф., Фоменко Н.А., Moscow State Technical University named after N.E. Bauman, Moscow ПРОГРАММНАЯ РЕАЛИЗАЦИЯ ДИСКРЕТНОЙ МАТЕМАТИЧЕСКОЙ МОДЕЛИ РАСЧЕТА ПРИБРЕЖНЫХ ВОЛНОВЫХ ПРОЦЕССОВ НА ОСНОВЕ РЕГУЛЯРИЗОВАННЫХ ПО Б.Н. ЧЕТВЕРУШКИНУ ЯВНЫХ СХЕМ НА ВЫЧИСЛИТЕЛЬНЫХ СИСТЕМАХ С МАССОВЫМ ПАРАЛЛЕЛИЗМОМ // Фундаментальные исследования. – 2015. – № 12-3. – С. 540-544;
URL: http://www.fundamental-research.ru/ru/article/view?id=39577 (дата обращения: 16.07.2019).

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

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