多网格谱元法及其在高性能计算中的应用*

针对传统谱元法在每个单元内只能存在单一均匀介质,应用在复杂非均匀介质的波传播模拟中可能造成极大计算规模的问题,发展了多网格谱元法。该方法在谱元法单元内引入独立的辅助网格,用于精细描述单元内的介质和外力分布变化,在较稀疏的主网格上进行波场的求解。基于声波和弹性波方程推导了多网格谱元法公式,并对几种典型模型进行了波场的数值模拟。与传统谱元法的对比结果表明,此算法在复杂非均匀介质的弹性波传播模拟中可以利用较少的网格点数达到不低于传统算法的精度。此外,实现了并行化的多网格谱元法,获得了较好的并行效率。