基于有限体积法的二维大地电磁各向异性数值模拟          下载免费PDF全文

王宁  汤井田  任政勇  肖晓  皇祥宇 《地球物理学报》2019,62(10):3912-3922

为了计算带任意地形的各向异性介质中二维大地电磁响应,本文在非结构化网格的基础上,采用有限体积法,开发了二维大地电磁各向异性正演模拟的新算法.首先,从Maxwell方程出发,推导二维各向异性介质中大地电磁场的边值问题;然后,采用三角网格自动生成技术对求解区域进行非结构化网格剖分,进而构建节点中心控制体积单元,利用有限体积方法,得到求解边值问题的大型稀疏线性方程组;最后,利用Pardiso精确地计算了大地电磁响应值.三个各向异性模型的计算结果表明,本文开发的有限体积算法,不仅能够高精度求解带任意地形的大地电磁电导率各向异性问题,而且对于同一模型,该方法的计算消耗和精度都与有限单元法相当.因此,有限体积法是处理电磁法各向异性问题的一种有效方法.  相似文献   

三维正交非结构网格数值模型中物理流场的定义方法与改进          下载免费PDF全文

刘晓辉  董礼先  Ralph T. Cheng 《海洋与湖沼》2010,41(4):621-627

使用Casulli等(2000)的数值方法建立了基于正交非结构网格的U型渠道的三维近岸正压水动力数值模型,检验了模型使用ELM(Eulerian-Lagrangian Method)方法处理平流项时,物理流场的定义对模拟结果的作用。数值试验和结果分析说明,使用ELM方法计算平流项所需的物理流场对计算结果影响明显,而目前常用的两种物理流场定义方法都存在一定不足——水位计算值对流速变化的反映不灵敏或者会产生小扰动等。针对常用物理流场定义方法的不足之处,作者在Casulli方法的基础上利用对多边形切向流速分量进行平均的方法改进了物理流场定义,抑制了小扰动。  相似文献