本文翻译自SeisSol官网(https://seissol.readthedocs.io/en/latest/index.html)
SeisSol
SeisSol是一种基于任意高阶精确导数不连续Galerkin方法(ADER-DG)的波传播和动态破裂模拟软件包。
SeisSol模拟软件的特点如下:
- 在时间和空间上使用任意高的近似阶;
- 使用四面体网格来逼近复杂的三维模型几何形状(断层和地形)和快速模型生成;
- 使用各向同性弹性、粘弹性和粘塑性材料来逼近真实的地下地质特性;
- 并行地理信息输入(ASAGI);
- 生成可靠和足够精确的合成地震图或其他地震数据集;
History
软件包 SeisSol ( http://www.SeisSol.org/)允许真实地模拟在复杂地球结构中传播的三维地震波场,这种地震波场是由一个有限的动态地震震源产生的,该震源受描述几何复杂断层上的断层应力和滑动之间关系的本构规律支配。SeisSol 是一种基于 Kaser 和 Dumbser (2006)提出的 ADER-DG 方法的高精度不连续 Galerkin 有限元求解器,能够精确地模拟断层摩擦破坏与地震波在传播过程中产生的最小色散误差,而波长本质上是消散的,去除了网格中无法解析的频率,而不影响更长的物理意义波长。De la Puente et al. (2009)和 Pelties et al. (2012)引入黎曼解来处理不连续的断层滑动条件和实现地震波的传播动力学。
该软件最近被证明在当前和未来的高性能计算基础设施上具有高度可伸缩性。在1992年兰德斯7.2级地震的开创性模拟中,它在世界上一些最大的超级计算机上达到了几千兆次/秒的性能。利用非平面地震断层结构模拟了工程频带(0ー10hz)的高精度破裂演化和地面综合振动。2017年初,SeisSol 进行了迄今为止时间最长、规模最大的动态破裂模拟,得益于当时的逐步推进(Uphoff 等人,2017年) ,解决了2004年苏门答腊-安达曼地震,包括复杂的扇形断层几何形状。该报纸获得了著名的美国国际超级计算大会协会的“最佳论文奖”。地震结果表明,承认几何复杂性、实际的断层性质和速度模型不仅影响震源动态,而且对合成的地面震动也有重要影响。该软件包以开源发行版的形式提供给社区使用( https://github.com/seissol/seissol )。