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

Аннотация

Построен алгоритм 3D моделирования электромагнитного поля для произвольного типа источника возбуждения N-слойной среды с иерархическим проводящим включением, расположенным в J-ом слое. Построены алгоритмы 2-D моделирования для дифракции звука и линейно поляризованной поперечной упругой волны на включении иерархической структуры, расположенной в J-ом слое N-слойной упругой среды. Целью данной работы является разработка программного комплекса для проведения математического моделирования распространения упругих возмущений, зарождающихся в гипоцентре землетрясения, в гетерогенных средах. В качестве источника возмущения используется геофизическая модель очага землетрясения в виде подвижки по уже существующему разлому. Для описания динамического поведения среды используется система уравнений (гиперболического типа) механики твердого тела с явным выделением контактных границ неоднородностей. Она решается с помощью сеточно-характеристического метода на криволинейных структурных трехмерных сетках. Одной из особенностей используемого вычислительного алгоритма является его высокая масштабируемость по количеству вычислителей. Использование криволинейных сеток позволяет с высокой точностью описать большой класс расчетных геометрий. Авторами была проведена математическая постановка задачи, разработан программный комплекс и выполнена серия численных расчетов. В статье приведены результаты расчета распространения сейсмического возмущения в многослойной геологической среде и оценки сейсмостойкости наземного сооружения. Проведено исследование ускорения процесса расчета при увеличении количества вычислительных узлов.