计算物理 ›› 2002, Vol. 19 ›› Issue (6): 537-543.

• 论文 • 上一篇    下一篇

用边有限元方法计算磁偶极子的三维电磁响应

沈金松   

  1. 石油大学(北京)盆地与油藏研究中心, 北京 102200
  • 收稿日期:2001-08-03 修回日期:2002-03-01 出版日期:2002-11-25 发布日期:2002-11-25
  • 作者简介:沈金松(1964-),男,江苏吴江,副教授,博士,主要从事地球物理测井理论、方法及应用的研究.
  • 基金资助:
    留学回国人员科研启动基金(J0201)

MODELING 3-D ELECTROMAGNETIC RESPONSES OF THE MAGNETIC DIPOLE IN FREQUENCY DOMAIN BY USING THE EDGE FINITE ELEMENT METHOD

SHEN Jin-song   

  1. Basin and Reservoir Research Center, Petroleum University, Beijing 102200, China
  • Received:2001-08-03 Revised:2002-03-01 Online:2002-11-25 Published:2002-11-25

摘要: 用边有限元基函数导出了麦克斯韦(Maxwell)方程的有限元关系式,计算了地下三维介质中磁偶极子的电磁场响应.将场分量定义在有限单元的边上,解决了结点有限元方法中场切向分量不连续的矛盾,保证了源除外的所有单元内有旋无散的特性.将总场分离成背景场和二次场,使该方法适用于任何方向的磁偶极子源.通过模拟算例分析了7种Krylov子空间迭代算法以及不完全乔累斯基分解预处理手段在解大型线性代数方程组中的计算效率和收敛特性.对比结果表明,施加不完全乔累斯基分解作预处理的广义乘积型双共轭梯度算法GPBiCG (Pbicg)收敛最快,是三维复杂介质电磁响应数值模拟的首选算法.

关键词: 边有限元方法, 磁偶极子, 双共轭梯度迭代, 不完全乔累斯基分解

Abstract: The finite element formulation of the Maxwell equation is derived by using the edge finite element method and is used to simulate the 3-D electromagnetic responses to the subspace medium.Through the definition of the electric field in the edge of each cell,the problem of the discontinuity of the tangential components of the electric field between different cells in node finite element method is solved and this method also keeps the property of the non-divergence and finite rotation.The separation of the total electric field into the background and the secondary fields makes it possible that the technique proposed in the paper is suitable for the magnetic dipole in any direction.The properties of the convergence of seven Krylov subspace iterative solvers and the effectiveness of the incomplete Cholesky decomposition pre-conditioner are analyzed through two simulation examples.Comparisons of different solvers show that the generalized productive bi-conjugate gradient iterative solver processed with pre-conditioner has the best convergence and it is the best selection in 3-D electromagnetic modeling.

Key words: edge finite element method, magnetic dipole, bi-conjugate gradient iterative solver, incomplete Cholesky decomposition

中图分类号: