收藏本站
《湘潭大学》 2012年
收藏 | 手机打开
二维码
手机客户端打开本文

几种求解H(curl)与H(grad)型偏微分方程有限元离散系统的多水平快速算法

刘春梅  
【摘要】:H(curl)型和H(grad)型偏微分方程组(PDEs)在电磁场计算和固体力学等领域具有广泛应用.本文针对频域电磁场PML逼近方程组、频域电磁场隐身超常材料模型和不定时谐Maxwell方程组等几类重要的H(curl)型偏微分方程组的棱有限元离散系统,以及Possion方程和线弹性力学方程组等几类重要的H(grad)型偏微分方程的间断Galerkin (DG)有限元离散系统和二次协调有限元离散系统,研究了相应的多水平快速算法.主要内容如下: 针对一类二维电磁场散射问题PML逼近方程的对称不定线性棱有限元离散系统,首次利用迭代两网格法的思想为其构造了一种迭代两网格法,该方法本质地将对原复杂离散系统的求解转化为三个更为简单的辅助系统的求解.同时为其中的两个辅助系统分别构造了高效预条件子,以及基于相应预条件子的PMINRES方法和PCG方法.数值实验表明,与经典的迭代法(如PMINRES法)相比,我们设计的迭代两网格法在迭代次数和CPU时间上均具有明显的优势. 通过将完美匹配层(PML)方法、电磁场理论以及隐身超常材料理论等相结合,给出了一种二维频域电磁场隐身超常材料模型,数值实验检验了模型的正确性.进一步,为一类电磁场散射问题PML逼近方程组和上述隐身超常材料模型方程组的线性棱有限元离散系统,结合HX预条件方法,构造了一种Schur补型的多水平预条件子,数值实验验证了算法的有效性. 针对不定时谐Maxwell方程组的棱有限元离散系统,构造了一种基于两水平方法和HX预条件子的多水平加性预条件子,本质性地将细网格空间上原不定问题的求解转化为粗网格空间和细网格上curl算子核空间上原不定问题的求解,以及在细网格空间上的对称正定问题的求解.当粗网格尺寸充分小时,证明了基于该预条件子的PGMRES法的一致收敛性.数值实验验证了理论结果的正确性. 针对一类二阶椭圆问题的一般对称内罚间断有限元离散系统,通过将DG有限元空间分解为“高频”部分和线性协调元空间,为其构造一种新的稳定的空间分解.利用该空间分解,首先为DG有限元离散系统设计了多水平预条件子.基于辅助空间预条件理论,我们证明了该预条件系统的条件数是一致有界的,同时为其中的线性元系统设计了BPX预条件子,从而得到一种多水平加性预条件子.接着还为DG有限元离散系统设计了一种两水平迭代法,并证明该迭代法是一致收敛的.数值实验验证了上述理论的正确性. 针对平面弹性力学问题二次有限元离散系统,给出了一种基于二分网格的局部多层网格(LMG)法.通过将二次有限元空间分解成‘高频’部分和线性元空间,结合二分网格和插值算子的特性,证明了该空间分解的稳定性和强Cauchy-Schwarz不等式成立,进而得到了LMG算法的一致收敛结果.数值实验验证了理论的正确性.
【学位授予单位】:湘潭大学
【学位级别】:博士
【学位授予年份】:2012
【分类号】:O241.82

知网文化
【参考文献】
中国期刊全文数据库 前7条
1 梁力,林韵梅;有限元网格修正的自适应分析及其应用[J];工程力学;1995年02期
2 王建华;线弹性有限元的自适应加密与多重网格法求解[J];河海大学学报;1994年03期
3 ;TWO-GRID DISCRETIZATION SCHEMES OF THE NONCONFORMING FEM FOR EIGENVALUE PROBLEMS[J];Journal of Computational Mathematics;2009年06期
4 杨强;翟明杰;周维垣;;三维弹性自适应有限元及其在拱坝分析中的应用[J];水力发电学报;2006年03期
5 王建华,李佳川;自适应网格加密和多重网格求解严重病态方程[J];上海力学;1996年01期
6 陈竹昌,王建华,王卫中;自适应多重网格有限元求解应力集中问题[J];同济大学学报(自然科学版);1994年02期
7 黄云清,陈艳萍;解非线性奇异两点边值问题有限元的一种分层迭代法[J];湘潭大学自然科学学报;1994年01期
中国博士学位论文全文数据库 前1条
1 钟柳强;求解两类Maxwell方程组棱元离散系统的快速算法和自适应方法[D];湘潭大学;2009年
中国硕士学位论文全文数据库 前1条
1 刘春梅;求解平面弹性问题的自适应有限元法及其多层网格法[D];湘潭大学;2009年
【共引文献】
中国期刊全文数据库 前10条
1 任仪;聂在平;赵延文;;基于高阶叠层基函数的加速迭代求解方法[J];电波科学学报;2008年01期
2 尹文禄;邓聪;赵菲;柴舜连;毛钧杰;;高阶矢量有限元方法实现及关键问题[J];电波科学学报;2009年02期
3 尹文禄;杨虎;肖科;柴舜连;毛钧杰;;高阶曲线矢量有限元方法实现及关键问题[J];电波科学学报;2011年04期
4 樊明武;颜威利;唐任远;;电磁场数值计算的发展、应用与前景(第五届国际电磁场计算会议评述)[J];电工技术学报;1987年02期
5 张秀敏,苑津莎,程志光;电磁场数值分析中棱单元矢量插值函数的研究[J];电工技术学报;2003年02期
6 张秀敏,苑津莎,徐永生,程志光;基于工程损耗模型的棱边有限元与节点有限元的算法比较[J];电工技术学报;2003年03期
7 孙佳佳;吴刚;;高阶矢量基函数在腔体本征值问题中的应用[J];电子科技;2012年02期
8 盛新庆;徐善驾;;一种系统构造三角形杂交边缘元空间的新方法[J];电子学报;1998年06期
9 任仪;聂在平;赵延文;马文敏;;矩量法中阻抗矩阵的稀疏化研究[J];电子学报;2007年12期
10 ;Meshing effects of the 3-D FEM numerical modeling in seismo-electromagnetics:An application in selectivity of seismic electric signal (SES)[J];Earthquake Science;2011年06期
中国重要会议论文全文数据库 前8条
1 刘春梅;钟柳强;舒适;肖映雄;;弹性力学问题自适应有限元方法及其局部多层网格法[A];中国计算力学大会'2010(CCCM2010)暨第八届南方计算力学学术会议(SCCM8)论文集[C];2010年
2 梁力;林韵梅;;瞬动态问题有限元网格单元修正的自适应分析[A];第五届全国结构工程学术会议论文集(第一卷)[C];1996年
3 王超;殷红成;黄培康;;四边参数曲面层叠高阶基函数的构造及分析[A];2007年全国微波毫米波会议论文集(上册)[C];2007年
4 Douglas N.Arnold;;Differential complexes and numerical stability[A];Structure-preserving Algorithms 2003--Proceedings of CCAST(World Laboratory) Workshop[C];2003年
5 R.Hiptmair;;Finite elements in computational electromagnetism[A];Structure-preserving Algorithms 2003--Proceedings of CCAST(World Laboratory) Workshop[C];2003年
6 R.Hiptrnair;;DISCRETE HODGE-OPERATORS: AN ALGEBRAIC PERSPECTIVE[A];Structure-preserving Algorithms 2003--Proceedings of CCAST(World Laboratory) Workshop[C];2003年
7 R.HIPTMAIR;;CANONICAL CONSTRUCTION OF FINITE ELEMENTS[A];Structure-preserving Algorithms 2003--Proceedings of CCAST(World Laboratory) Workshop[C];2003年
8 汪文帅;张怀;李小凡;;间断的Galerkin方法在地震波场数值模拟中的应用概述[A];中国科学院地质与地球物理研究所2013年度(第13届)学术论文汇编——地球深部结构与过程研究室[C];2014年
中国博士学位论文全文数据库 前10条
1 任志刚;预条件算法及在电磁场数值模拟中的应用[D];电子科技大学;2010年
2 黄卓红;PDE离散方程组和鞍点问题的预处理方法[D];电子科技大学;2010年
3 尹文禄;高阶矢量有限元方法在电磁领域中的研究及应用[D];国防科学技术大学;2010年
4 刘洋;非标准混合元方法分析及数值模拟[D];内蒙古大学;2011年
5 汪波;求解时域麦克斯韦方程组的间断伽略金方法[D];湖南师范大学;2011年
6 朱剑;复杂电磁问题的有限元、边界积分及混合算法的快速分析技术[D];南京理工大学;2011年
7 李胜坤;鞍点问题和Sylvester型矩阵方程(组)的数值解法研究[D];电子科技大学;2011年
8 王振;制备MgO晶体电弧炉的建模研究[D];大连理工大学;2011年
9 戴培良;工程力学中的有限元方法及其误差估计[D];南京航空航天大学;2003年
10 何小祥;基于有限元法的混合技术及其在复杂电大腔体电磁分析中的应用[D];南京航空航天大学;2003年
中国硕士学位论文全文数据库 前10条
1 王瑞英;H(curl)空间的各向异性分析[D];郑州大学;2010年
2 王娅娜;平面弹性问题协调的矩形混合元[D];郑州大学;2010年
3 孙刚;外圈旋转轴承多点加载装置研究[D];哈尔滨工程大学;2010年
4 魏宁;定常非线性薛定谔方程的两网格有限元方法[D];湘潭大学;2010年
5 叶珍;微波管输入输出窗及任意结构腔体有限元理论与CAD技术[D];电子科技大学;2011年
6 郑瑞瑞;Sobolev方程的两类变网格混合元方法[D];山东大学;2011年
7 罗平;地下水污染模型与麦克斯韦方程的数值方法[D];山东大学;2011年
8 周洁;求解一类Maxwell特征值问题的两网格法和自适应方法[D];湘潭大学;2011年
9 彭妮娜;一种基于GIF图像轮廓曲线的磨光算法[D];湘潭大学;2011年
10 刘鸣;特种工况轴承试验台的研制[D];哈尔滨工程大学;2011年
【二级参考文献】
中国期刊全文数据库 前10条
1 梁力,林韵梅;有限单元法的非结构化网格及其自动生成[J];东北大学学报;1994年03期
2 梁力,林韵梅;有限元网格修正的自适应分析及其应用[J];工程力学;1995年02期
3 杨一都;;特征值问题迭代伽略金法与Rayleigh商加速[J];工程数学学报;2008年03期
4 陈艳萍,黄云清;奇异非对称两点边值问题的有限元解的整体高精度[J];高等学校计算数学学报;1994年03期
5 ;An Adaptive Uniaxial Perfectly Matched Layer Method for Time-Harmonic Scattering Problems[J];Numerical Mathematics:Theory,Methods and Applications;2008年02期
6 王建华;线弹性有限元的自适应加密与多重网格法求解[J];河海大学学报;1994年03期
7 ;Uniform convergence of multigrid V-cycle on adaptively refined finite element meshes for second order elliptic problems[J];Science in China(Series A:Mathematics);2006年10期
8 ;The order-preserving convergence for spectral approximation of self-adjoint completely continuous operators[J];Science in China(Series A:Mathematics);2008年07期
9 周张义;李芾;;P型自适应有限元法及其在机车车辆强度分析中的应用[J];机车电传动;2007年02期
10 肖尚宏,王石刚,余俊,周济;二维自适应有限元网格生成[J];计算机工程;1993年05期
【相似文献】
中国期刊全文数据库 前10条
1 刘兴平;广义预条件迭代方法[J];计算数学;1990年01期
2 卫加宁,王伟沧,皮新明,章社生;一类虚拟边界预条件多重网格并行算法[J];武汉理工大学学报(交通科学与工程版);2001年01期
3 黄文彬,许传炬;Poisson方程谱元法的一个有限元预条件分析[J];厦门大学学报(自然科学版);2003年04期
4 孙丽英;线性方程组的AOR预条件迭代法的两个性质[J];广东教育学院学报;2003年02期
5 庄伟芬,卢琳璋;(I+S_(max))预条件Gauss-Seidel迭代法进一步探索[J];厦门大学学报(自然科学版);2004年S1期
6 朱振宇,裴江云,吕小林,刘洪,李幼铭;基于预条件共轭梯度法的盲源反褶积方法(英文)[J];Petroleum Science;2004年03期
7 李铮,邵新慧,李长军;矩阵B~TA~(-1)B的特征值估计及预条件处理[J];东北大学学报(自然科学版);2005年06期
8 雷刚;畅大为;王慧勤;;预条件[I+C(α)]加速2PPJ型方法的收敛性[J];江汉大学学报(自然科学版);2006年01期
9 卫加宁,章社生,Yakup Paker;区域分裂并行计算中预条件迭代若干性质[J];武汉理工大学学报(交通科学与工程版);2001年02期
10 雷刚;;两类预条件后SOR迭代法收敛性的比较[J];江南大学学报(自然科学版);2011年02期
中国重要会议论文全文数据库 前10条
1 梅金顺;刘洪;;预条件方程组及其应用[A];中国科学院地质与地球物理研究所二○○四学术论文汇编·第三卷(油气·矿产·水资源)[C];2004年
2 刘喜武;刘洪;;实现稀疏反褶积的预条件双共轭梯度法[A];中国科学院地质与地球物理研究所二○○三学术论文汇编·第四卷(油气资源)[C];2003年
3 龙桂华;李小凡;张美根;;粘弹性声波波形反演中的预条件最速下降法研究[A];中国地球物理学会第二十四届年会论文集[C];2008年
4 李正光;吴柏生;;结构单元和结点删除后的一种重分析方法[A];第十二届全国结构工程学术会议论文集第Ⅰ册[C];2003年
5 李振华;王彦飞;杨长春;;一种能提高叠前偏移剖面分辨率的预条件正则化反演方法[A];中国地球物理·2009[C];2009年
6 谷同祥;刘兴平;;预条件多搜索方向共轭梯度方法[A];中国工程物理研究院科技年报(2003)[C];2003年
7 孙树立;陈璞;;求解多右端向量方程组的块共轭梯度法及其相关研究进展[A];计算爆炸力学进展[C];2006年
8 胡玉生;;基于GMRES的边界元一有限元混合法在EMC中的应用[A];第17届全国电磁兼容学术会议论文集[C];2007年
9 吴建平;李晓梅;;块三对角线性方程组不完全分解预条件的一种并行化方法[A];全国计算物理学会第六届年会和学术交流会论文摘要集[C];2007年
10 牛臻弋;徐金平;;PO-MoM方程迭代求解中的稀疏化递归Cholesky分解预条件技术[A];2005'全国微波毫米波会议论文集(第二册)[C];2006年
中国重要报纸全文数据库 前1条
1 张未名 作者单位:中国科学院上海微系统与信息技术研究院;环保要与发展阶段相适应[N];中国经济导报;2005年
中国博士学位论文全文数据库 前10条
1 索奇峰;预条件可压流算法在桥梁风工程中的应用研究[D];西南交通大学;2012年
2 李正光;结构布局修改静力重分析的预条件共轭梯度法[D];吉林大学;2004年
3 吴建平;稀疏线性代数方程组迭代法中的预处理技术研究[D];中国人民解放军国防科学技术大学;2002年
4 张振宇;数值线性代数中的若干问题[D];复旦大学;2003年
5 单润红;Toeplitz系统的高效预处理技术和分布式并行算法研究[D];国防科学技术大学;2004年
6 王振营;人口迁移的规律[D];中国人民大学;1993年
7 秦梅;一类奇异方程组的求解和扰动分析[D];复旦大学;2006年
8 陈(日午);关于真实感服装生成算法的研究[D];西北工业大学;2000年
9 刘福体;线性方程组迭代求解及相关问题的研究[D];电子科技大学;2009年
10 张胜;不定椭圆问题有限体积法的算法研究[D];复旦大学;2003年
中国硕士学位论文全文数据库 前10条
1 白文军;线性方程组的预条件迭代法和误差估计[D];陕西师范大学;2010年
2 王炜;线性方程组新的预条件迭代法的研究[D];太原理工大学;2012年
3 孟国艳;SAOR-AL预条件共轭梯度法[D];太原理工大学;2003年
4 刘翠;离散正弦变换与TOEPLITZ系统的快速算法与并行算法[D];国防科学技术大学;2004年
5 薛秋芳;一类矩阵的AOR迭代收剑性分析及两种预条件迭代的收敛性定理[D];陕西师范大学;2005年
6 柳卫东;一类预条件迭代法研究[D];陕西师范大学;2008年
7 杨勇;大型线性方程组不完全分解预条件的研究[D];电子科技大学;2013年
8 魏洪川;超大规模集成电路中寄生电容的提取和快速解法研究[D];清华大学;2003年
9 赵丹;两类矩阵的预条件迭代法的收敛性[D];扬州大学;2009年
10 卢晓平;精化块共轭梯度法及其应用[D];复旦大学;2007年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026