ИСТИНА |
Войти в систему Регистрация |
|
ИПМех РАН |
||
Данная работа служит основой для моделирования зон субдукции. Для этого решаются уравнения Навье-Стокса для вязко-пластической среды. В качестве базы, для численного решения уравнений Навье-Стокса были выбраны методики, предложенные в работах исследователя Т.В.Гери [1] и использование алгоритма SIMPLE на разнесённой (шахматной) сетке. Плюсом подхода является то, что при решении уравнений моментов, мы имеем известные значения давлений в соседних узлах, что избавляет нас от необходимости искать приближенное значение давлений в узлах сетки граничащих с узлами, в которых определены значения скоростей. Алгоритм заключается в следующем: сначала решаются уравнения моментов Х и Y, каждое относительно своих узлов сетки, при этом давления полагаются известными с предыдущего шага. Затем, решается уравнение сохранения массы, тем самым уточняется новое значение давления, которое снова подставляется в уравнения моментов и ищется решение уравнений моментов, которое будет удовлетворять закону сохранения массы. Для решения поставленной задачи написана первая версия кода программы на языке C++. Выбор языка определяется его низким уровнем, что должно обеспечивать более быстрый просчёт данных. Распределение плотностей выбрано в соответствии с моделью emc stw-105. Значения для температур взяты с сайта Аризонского института геофизики. Величина изменения ускорения свободного падения рассчитывается с использованием данных о плотностях. Плотность полагается зависящей от давления и температуры одновременно. Вязкость описывается согласно уравнению из статьи [2]. Обсуждаются первые результаты выполненных расчетов.