论文部分内容阅读
从上个世纪六十年代以来,计算电磁学得到了迅猛的发展,并得到了广泛的应用。经过四十多年众多计算电磁学和计算数学领域科学工作者的共同努力,开发了多种有效的数值计算方法求解各种电磁场问题。这些方法主要有时域有限差分方法、有限元方法、有限体积方法和矩量法等。这些方法在各种电磁场问题中得到了广泛的应用。然而这些方法又各自有自己的缺点,比如:时域有限差分方法并不适合计算几何结构复杂的问题;有限元方法和矩量法由于单元之间互相依赖,因此这两种方法并不适合并行计算;有限体积方法不适合使用高阶基函数。为了利用这些方法的优点,同时又尽量回避这些传统方法的缺点,本论文研究了一类较新的数值计算方法,这类方法是间断伽辽金有限元方法,并实现了其中一种基于节点的间断伽辽金有限元方法。间断伽辽金有限元方法是一种基于有限体积方法和有限元方法而发展出来的一种方法。我们研究了如何使用节点间断伽辽金有限元方法求解电磁场谐振腔问题、电磁波传播问题、电磁目标散射问题、不确定形状电磁目标随机散射等问题。电磁场谐振腔问题是电磁场领域的基本问题之一。在研究运用一维到三维的节点间断伽辽金有限元方法求解Maxwell方程组时,我们首先用该方法求解了电磁场谐振腔问题。电磁场传播问题也是电磁场领域的基本问题之一。我们首先研究了如何使用节点间断伽辽金有限元方法求解电磁波在自由空间的传播,然后研究了使用该方法求解电磁波在有耗介质中的传播问题。接着,我们重点研究了使用节点间断伽辽金有限元方法求解电磁场散射问题。电磁场散射问题在电磁场领域的研究具有非常重要的意义。自从第二次世界大战雷达发明以来,不确定形状目标散射的计算一直是电磁场研究领域的热点问题,其对于军事目标的隐身和反隐身具有不言而喻的重要意义。首先,我们研究了二维金属圆柱的散射问题,然后研究了二维金属方柱的散射问题。最后,我们在已有开源网格生成程序基础上开发了计算不确定形状目标散射问题所需要的复杂动态网格生成程序,并结合稀疏网格积分方法和随机配置方法研究了不确定形状目标散射问题。稀疏网格积分方法(Sparse Grid Method)是在求解高维积分问题时大大减少计算量的一种有效方法。随机配置方法(Stochastic Collocation Method)则为求解随机问题提供了一种和蒙特-卡洛方法(Monte-Carlo Method)一样容易实现,但是其比蒙特-卡洛方法收敛更快。为了提高计算速度,我们提出了基于Richardson外插方法的自适应稀疏网络积分方法。其核心思想是在稀疏网格积分方法中从第三层起使用Richardson外插方法,由当前层和上一层的计算结果来预测下一层的结果。如果预测的结果与当前层的结果之间的误差在设计的误差范围之内,则停止计算,反之,则继续计算下一层,从而加速求解不确定目标随机散射问题。当所求解问题的维数越高,节省的计算时间就越多。随着现代科学研究的不断发展,科学研究和工程要解决的问题规模越来越大。科学计算和仿真是理论研究和实验研究之外的第三种重要的方法。与此同时,人们对求解大规模的电磁场问题需求也日夜增长。单台个人计算机或工作站的计算能力通常满足不了计算需求。因此,并行计算是求解大规模电磁场问题的唯一途径。间断伽辽金有限元方法由于单元之间通过数值通量进行耦合,只有单元边界上的值与相邻单元有联系,所以特别容易实现并行计算。我们研究了在超级计算机环境下三维的间断伽辽金有限元方法实现。在实现中,我们采用了开源的基于消息传递接口的并行编程函数库,以及使用了ParMETIS快速并行网格划分函数库进行负载均衡。间断伽辽金有限元方法在空间离散时通常采用高阶的多项式基函数。因此,为了使整个计算格式保持高阶精度,间断伽辽金有限元方法中采取的时间离散方法通常也采用高阶的常微分方程求解方法,比如:两阶的蛙跳格式,两到四阶的龙格-库塔方法,多步的预测器与校正器方法等。但是这些方法都是基于解析或半解析方法推导的。因此,这些方法的精度是固定的。要想推导同样的更高精度格式,非常困难乃至不可能。本文研究并实现了一种完全基于数值方法构造的时间离散格式。这种新格式,在形式上和传统的预测器与校正器方法一样。但是,不同的是,新格式不是采用多项式函数近似常微分方程的解,而是采用指数函数来近似常微分方程的解。所有的这些研究为进一步推广间断伽辽金有限元方法在计算电磁学领域的应用提供了较好的参考价值。