国内外计算电磁学研究的回顾与展望.docxVIP

国内外计算电磁学研究的回顾与展望 1 广泛而广泛的组织学术活动 计算电磁学是一门综合磁体理论、数值计算方法和计算机软件技术的新兴学科。他的研究历史可以追溯到半个多世纪前（例如，在20世纪40年代提出了有限的差分法），但他的增长是近30年的。众所周知,计算机技术的持续进步和日益普及、计算数学和软件技术的快速发展、电气装备的工业需求和控制技术的广泛应用是促进计算电磁学研究的关键因素。与此同时,广泛而有组织的学术活动不仅加速了研究进程,也是计算电磁学积极发展的标志。其中,在国际学术交流方面,COMPUMAG会议(Conference on Computation of Electromagnetic Fields)是这一领域最重要的世界性会议之一。该会议由1974年在英国卢瑟福阿普尔顿国家实验室召开的国际指导委员会发起,自1976年起至今,大约两年一次,分别在英、法、美、德、意大利、奥地利、日本、巴西等国家举行了共14届会议。经过提前5年的申办与准备,第15届COMPUMAG会议于2005年6月在中国沈阳举行。这也是这一系列会议首次在中国举行。本文拟从第15届COMPUMAG会议的概述开始,讨论近年来国内外计算电磁学研究的若干重要成果。 2 会议论文形式 2005年6月26日至30日,第15届COMPUMAG会议在中国沈阳成功召开。实际到会代表354人,来自30个国家和地区,其中国内代表73人,国外代表281人。会议共收到以2页短文形式提交的论文680篇,经过世界各地的178名评审专家的网上评审,共接受了508篇论文刊登在会议论文(短文)集上。参加交流的论文数超过了这一系列会议历史上最好的法国会议(COMPUMAG 2001,Evian,提交论文651篇,接受的会议交流论文427篇)。目前参加交流的论文正以扩展的全文形式(邀请报告论文为6页,其它为4页)参加第二阶段评审,预计将有约250篇刊登在IEEE的磁学会刊上。各专题会的主题与上届会议(COMPUMAG 2003,Saratoga Springs,USA)类似,表1列出了专题内容与所交流的文章数。 3 计算公式的发展 近30年来,计算电磁学领域已经取得了重大的科学技术进步。从二维到三维,从线性到非线性,从单一电磁场问题到电磁场与电路系统或与其它物理场的耦合问题,从正问题到包括优化技术在内的逆问题,计算能力有了飞跃的提高。在第15届COMPUMAG会议的开幕大会上,有两个邀请报告分别从国际范围和我国国内的研究情况出发综述了计算电磁学的近期发展,本节拟以这两个报告为线索,说明这一领域的若干重要研究成果。 有限元法是应用最为广泛的数值方法,下面列出的前4个重要成果均与有限元法密切相关。 3.1 基于大型工程问题的有限元分析 ICCG法即不完全乔累斯基分解-共轭梯度法(Incomplete Cholesky Conjugate Gradient Method),是一种预处理共轭梯度法。该方法最初由J.A.Meijerink和V.der Vorst提出,用于求解有限元分析中的大型稀疏对称代数方程组,目前已广泛用于大型工程问题的有限元计算,并出现了一些改进形式。这一方法的突出优点是,其计算次数随nlogn增长,而不是像高斯消去法及其派生方法那样随n3增长,n为未知数个数;其存储量需求独立于系数矩阵的带宽,也大大小于高斯消去法。ICCG法的应用使得在1980年以前很难进行的许多三维大型工程问题的有限元计算得以实现,但目前研究者仍试图改进这一方法,因为包括涡流问题、耦合问题在内的复杂工程问题其离散化方程的系数矩阵非正定,也并不满足原始ICCG法要求的“M矩阵”条件,从而在应用ICCG法时收敛性变坏。 3.2 基于自适应生成的三角形 对于一组给定的节点,怎样的网格才是最好的三角形网格?许多研究者认为,迪朗尼三角化方法最适合完成这一任务。该方法由俄国数学家Delaunay提出(1934),Cendes(1983)首先将这一方法应用于二维有限元网格自适应生成。这一方法能使所生成三角形的最小角度之和极大化,也就是所生成的三角形在一定条件下最接近于等边三角形,于是能避免细长畸形单元的产生,这正是有限元分析所要求的。迪朗尼法已从二维推广到三维,此时需要将三角形转变成四面体,将二维情况下的外接圆转变成外接球。 3.3 确定磁储能边界 开域问题的边界在无穷远处,求解区域无限大,在有限元计算中无法实现。最通用的经验式权宜之计是截断法,即人为地设定一个远离场源的边界,认为在该边界上电磁场已衰减到零,将无限域截断成有限域。为了确定这一边界的设定是否合理,通常需要以重要场域内的总体参数(例如电感,电容或电磁储能)在不同边界设定下的变化为判据,经过多次试算,最后设定一个满足给定误差要求的相对较小的区域作为求解域,显然