无限元-谱元混合法在2.5维引力位计算中的应用

引力位是地球自由振荡数值模拟中不可或缺的一部分,也是重力异常研究中需要计算的对象。由于引力位满足无界的泊松-拉普拉斯方程,其在无穷远边界处为零,这对数值模拟造成了困扰。针对此,采用无限元-谱元混合法,直接对无穷远边界进行拟合,不再对边界条件做近似;同时考虑到三维地球模型计算效率问题,采用2.5维控制方程;最后,通过数值试验,验证了方法的准确性。