电磁场方程及其特征值问题高效高精度数值方法

结题报告

项目介绍

AI项目解读

基本信息

批准号：
11871092
项目类别：
面上项目
资助金额：
55.0万
负责人：
张智民
依托单位：
北京计算科学研究中心
学科分类：
A0504.微分方程数值解
结题年份：
2022
批准年份：
2018
项目状态：
已结题
起止时间：
2019-01-01 至2022-12-31

项目参与者：
安静；王丽修；张倩；
关键词：
误差估计电磁场方程谱元超收敛协调有限元

项目摘要

Electromagnetic equations have important applications in many fields such as target recognition and imaging, stealth and anti-stealth, telemetry and remote sensing, underground target detection, microwave device design, etc. Nédélec constructed a class of finite elements belonging to H(curl) , which starts a new era for solving electromagnetic equations by so-called "edge elements". In this project we plan to: 1) construct a series of H^2(curl) finite elements, based on which, all electromagnetic equations with order no more than 4 can be descriticized by conforming finite element methods. Then, we explore superconvergence properties of H(curl) and H^2(curl) elements, and extend Polynomial Preserving Recovery to electromagnetic equations. 2) propose a spectral Galerkin method based on a dimension reduction scheme for the eigenvalue problem of the electromagnetic field in the isotropic medium on the polar/spherical geometry, and a spectral element method for the eigenvalue problem of electromagnetic field in an anisotropic medium on general domains. Finally, we provide error estimation for every numerical method and prove the efficiency of the algorithm by theory and numerical experiments.

电磁场方程在目标识别与成像、隐身与反隐身、遥测遥感、地下目标探测、微波器件的设计等很多领域都有重要的应用。 Nédélec 构造了一系列H(curl)有限元，开启了棱元求解电磁场方程的新时代。本项目的研究目标有两个: 1) 构造一系列H^2(curl)有限元从而使得不高于4阶的电磁场方程均可用协调有限元方法离散; 进而，探究H(curl)和H^2(curl) 元的超收敛性质，同时将保多项式恢复型算子(Polynomial Preserving Recovery)推广到电磁场方程的数值解中。2) 提出极/球几何体上各向同性介质电磁场特征值问题基于降维格式的一种高效高精度的谱Galerkin方法，对于一般区域上各向异性介质电磁场特征值问题提出一种有效的谱元方法。最后, 我们将对提出的数值方法给出理论分析, 进行误差估计，并通过数值试验来验证算法的有效性。

结项摘要

电磁场方程在很多科学领域有着重要的应用。1980年代 Nédélec 构造了一系列H(curl)-协调有限元，开启了棱元求解电磁场方程的新时代。随着科学技术的发展，数值求解含有4阶旋度算子的电磁场方程提上了日程，因而需要构造H(curl2)-协调有限元。本项目成功地构造出一系列H(curl2)-协调有限元，包括二维三角形单元、矩形单元、一般四边形单元，和三维四面体单元、立方体单元上的H(curl2)-协调有限元。这些单元的构造基于de Rham 复型，因而是系统性的几族任意阶方法，包括最低阶的6个自由度的三角形单元、8个自由度的四边形单元、以及18个自由度的四面体单元。从而使得不高于4阶的电磁场方程均可用协调有限元方法有效离散求解。我们还探究了H(curl2)-协调元的收敛性和超收敛性质，同时将保多项式恢复型算子(Polynomial Preserving Recovery)推广到电磁场方程的数值解中。本项目还提出并构造了圆盘、椭圆域、以及三维球几何体上各向同性介质电磁场特征值问题基于降维格式的一种高效高精度的Galerkin多项式谱方法，并将其应用于数值求解三维球域中Maxwell传输特征值问题，得到了具有几何收敛阶（谱精度）的数值逼近。对于这个新的方法，我们建立了数学理论，进行了误差估计，数值实验显示了算法的有效性。这是目前解三维球域中Maxwell传输特征值问题唯一的具有谱精度的数值方法。