Coulomb规范下地电磁场的自适应有限元模拟的理论分析

被引：26

作者：

汤井田

任政勇

化希瑞

机构：

[1] 中南大学信息物理工程学院

来源：

地球物理学报 | 2007年 / 05期

关键词：

地电磁场; Coulomb规范; 矢量势-标量势; 后验误差; 自适应有限元;

D O I：

暂无

中图分类号：

P631.3 [电法勘探];

学科分类号：

0818 ; 081801 ; 081802 ;

摘要：

地电磁场的直接求解法存在伪解现象,且电磁场分量在界面上的不连续性与节点型有限元的基本要求矛盾.本文将Coulomb规范下磁矢量势-电标量势与自适应有限元相结合,提出了地球物理电磁场计算的快速、高精度方法.首先从地电磁场一般边值问题出发,给出了Coulomb规范下磁矢量势-电标量势的公式系统,分析了求解域内势的连续性.采用Galerkin加权余值法推导出积分弱解形式和Delaunay非结构化四面体单元时Hierarchal基函数的有限元方程.基于超收敛恢复技术,提出了适用于电磁场的后验误差估计方法,阐述了地电磁场自适应计算的策略及迭代算法,分析了计算时间消耗和误差收敛性质,表明本文方法可以用最优的计算资源得到呈拟指数收敛到准确解的数值结果,从而为后续的数值计算奠定了理论基础.

引用

页码：1584 / 1594

页数：11

共 10 条

[1] 任意地球物理模型的三角形和四面体有限单元剖分 [J].

汤井田 ;

任政勇 ;

化希瑞 .

地球物理学进展, 2006, (04) :1272-1280

[2] 电导率分块均匀的瞬变电磁2.5维有限元数值模拟 [J].

熊彬 ;

罗延钟 .

地球物理学报, 2006, (02) :590-597

[3] 二维有限元网格自动自适应生成 [J].

刘荣 ;

郑永红 ;

游亚戈 .

计算机辅助设计与图形学学报, 2005, (11) :79-84

[4] 复杂介质有限元法2.5维可控源音频大地电磁法数值模拟 [J].

底青云 ;

Martyn Unsworth ;

王妙月 .

地球物理学报, 2004, (04) :723-730

[5] 大地电磁法三维交错采样有限差分数值模拟 [J].

谭捍东 ;

余钦范 ;

John Booker ;

魏文博 .

地球物理学报, 2003, (05) :705-711

[6] 有限元网格生成方法研究的新进展 [J].

关振群 ;

宋超 ;

顾元宪 ;

隋晓峰 .

计算机辅助设计与图形学学报, 2003, (01) :1-14

[7] 大地电磁二次函数插值的有限元法正演模拟 [J].

史明娟 ;

徐世浙 ;

刘斌 .

地球物理学报, 1997, (03) :421-430

[8]

Grid generation and optimization based on centroidal Voronoi tessellations[J] . Qiang Du,Max Gunzburger.Applied Mathematics and Computation . 2002 (2)

[9]

Superconvergence and the superconvergent patch recovery[J] . O.C. Zienkiewicz,J.Z. Zhu.Finite Elements in Analysis & Design . 1995 (1)

[10] THE P-VERSION OF THE FINITE-ELEMENT METHOD [J].

BABUSKA, I ;

SZABO, BA ;

KATZ, IN .

SIAM JOURNAL ON NUMERICAL ANALYSIS, 1981, 18 (03) :515-545

← 1 →