基于谱元法的低频平面波电磁场数值模拟
来源期刊:桂林理工大学学报2020年第4期
论文作者:万文武 李长伟 熊彬 高磊 宋卫文
关键词:谱元法;GLL多项式;正演;精度;
摘 要:谱元法是近年来发展成熟的一种数值方法,具有高精度和指数收敛特性。基于谱元法实现了大地电磁的正演模拟,从麦克斯韦方程出发得到电场矢量Helmholtz方程,使用Gauss-Lobatto-Legendre多项式作为基函数离散Helmholtz方程,采用稳定的双共轭梯度算法求解大型稀疏线性方程组得到电磁场分布。通过与国际标准模型COMMEMI2D-0、 COMMEMI2D-1的计算对比,验证了算法的有效性和准确性,对几个典型的地电模型进行模拟计算,表明了算法具有有限元的灵活性和谱方法高精度的特性,是一种高效实用的电磁正演计算方法。
万文武,李长伟,熊彬,高磊,宋卫文
摘 要:谱元法是近年来发展成熟的一种数值方法,具有高精度和指数收敛特性。基于谱元法实现了大地电磁的正演模拟,从麦克斯韦方程出发得到电场矢量Helmholtz方程,使用Gauss-Lobatto-Legendre多项式作为基函数离散Helmholtz方程,采用稳定的双共轭梯度算法求解大型稀疏线性方程组得到电磁场分布。通过与国际标准模型COMMEMI2D-0、 COMMEMI2D-1的计算对比,验证了算法的有效性和准确性,对几个典型的地电模型进行模拟计算,表明了算法具有有限元的灵活性和谱方法高精度的特性,是一种高效实用的电磁正演计算方法。
关键词:谱元法;GLL多项式;正演;精度;