间断Galerkin有限元地震数值模拟算法的自适应实现方法

摘要

一种间断Galerkin有限元地震数值模拟算法的自适应实现方法，包括以下步骤：一、定义网格剖分密度，按网格密度对计算区域进行非结构网格剖分；二、确定每个网格单元上的多项式展开阶数p；三、对相邻的两个网格单元解耦，选择合适的数值流通量格式作为相邻单元间波场交换方式；四、对解耦方程数值化，获得每个网格单元上的空间离散方程；五、求相邻两个单元公共边界上波场函数之间的内积，获得边值映射矩阵；六、选择合适的时间积分格式；七、利用边值映射矩阵对不同阶数的相邻单元波场进行映射，计算空间离散部分，施加震源条件，更新波场分量，获得模拟结果。