耦合岩浆/地幔动力学模型的分块预条件子分析
数值分析
2016-07-08 v3 计算工程、金融与科学
地球物理
摘要
本文考虑岩浆动力学方程有限元离散的迭代求解。简化形式的岩浆动力学方程与斯托克斯方程具有某些共同特征。因此，我们为岩浆动力学构造、分析并数值测试了一种 Elman、Silvester 和 Wathen 型分块预条件子。我们从分析上证明并从数值上展示了该预条件子的最优性。所呈现的分析揭示了预条件子对岩浆动力学方程中参数的依赖性，这些参数可能影响迭代线性求解器的收敛性。通过一系列二维和三维非结构化网格上的数值算例——从简单的说明性问题到俯冲带类几何体上的大型问题——验证了该分析。用于重现所有数值算例的计算机代码可作为支撑材料免费获取。
引用
@article{arxiv.1311.6372,
title = {Analysis of block-preconditioners for models of coupled magma/mantle dynamics},
author = {Sander Rhebergen and Garth N. Wells and Richard F. Katz and Andrew J. Wathen},
journal= {arXiv preprint arXiv:1311.6372},
year = {2016}
}