0049 - 三维电磁散射有限元方法的现代数学观点：微分形式、de Rham 复形与 FEEC

Wed, 15 Apr 2026 14:20:00 +0800

总结

这篇笔记把三维电磁散射有限元放回更现代的数学框架里理解。核心结论不是“棱边元更适合 Maxwell 方程”这么一句经验话，而是：电场天然是 1-形式，离散空间必须保持 de Rham 复形的拓扑结构；Whitney 1-形式也就是 Nedelec 棱边元，正是这个结构在离散层面的自然实现。

从这个角度看，很多熟悉但分散的知识都被统一起来了：为什么节点元会出现伪解，为什么刚度矩阵主要来自外导数而与材料无关，为什么质量矩阵体现各向异性，为什么需要交换图投影，为什么 FEEC 能提供统一的稳定性与误差理论，为什么低频崩溃、高频污染、PML、hp 自适应都可以放进同一个框架里讨论。

如果只记三句话，可以记成：

Maxwell 方程的离散化首先是一个拓扑问题，其次才是一个数值线性代数问题。
棱边元不是经验修补，而是 1-形式离散化的自然结果。
FEEC 的价值在于：把“一个个问题分别发明元素”的做法，提升为“从复形结构自动生成稳定离散化”的统一方法。

第零章：为什么需要现代观点？
第一章：微分形式，电磁场的自然语言
第二章：de Rham 复形，有限元的拓扑骨架
第三章：离散 de Rham 复形与 Whitney 形式
第四章：变分框架与 Hodge-Laplace 问题
第五章：有限元外微积分 FEEC
第六章：Helmholtz 分解与数值困难
第七章：误差分析的现代观点
第八章：散射问题的函数空间框架
第九章：高阶方法与 hp 自适应
第十章：从理论到计算的桥梁

我的整理笔记

三维电磁散射有限元方法的现代数学观点

一套自洽的知识手册

第零章：为什么需要现代观点？

传统观点的局限

传统教科书常把有限元过程讲成一条计算链：

$$ \text{Maxwell方程} \xrightarrow{\text{弱形式}} \text{变分问题} \xrightarrow{\text{选基函数}} \text{矩阵方程} \xrightarrow{\text{求解}} \text{数值解} $$

你会学到“棱边元比节点元更适合 Maxwell 方程”，也会学到一整套局部矩阵、全局组装和边界条件处理方法；但如果只停留在这一层，就很难真正回答下面这些问题：

为什么电场偏偏要用棱边元，而不是节点元？
为什么节点元会出现伪解？
为什么刚度矩阵和质量矩阵在物理含义上差别这么大？
为什么有些稳定性结论看上去像是“经验总结”，但其实背后有统一结构？

现代观点的统一

现代数学给出的是一个统一框架：

0045 - 棱边元有限元里的全局编号与局部编号：从弱形式离散到稀疏矩阵组装

Wed, 15 Apr 2026 12:18:00 +0800

这篇笔记在做什么

这篇笔记专门抽出原始会话里最难也最容易混乱的部分：全局基函数如何代入连续弱形式，为什么绝大多数单元积分自动为零，局部基函数怎样带着方向符号回到全局矩阵。

我保留了推导里的每一步筛选与变换，让“全局到局部、局部回到全局”不再只是口号，而是一条完整可验证的计算链。

从全局到局部，再从局部到全局

——弱形式离散化的完整推导

一、出发点：连续弱形式

求 $\mathbf{E}^{sca} \in H(\text{curl};\Omega)$，使得对所有 $\mathbf{T} \in H(\text{curl};\Omega)$：

$$ \begin{aligned} &\underbrace{\int_\Omega (\nabla\times\mathbf{T})\cdot(\nabla\times\mathbf{E}^{sca})\,dV}_{a_1(\mathbf{T},\mathbf{E}^{sca})} \\ &\quad - k_0^2 \underbrace{\int_\Omega \mathbf{T}\cdot\bar{\bar{\varepsilon}}_r\mathbf{E}^{sca}\,dV}_{a_2(\mathbf{T},\mathbf{E}^{sca})} \\ &\quad - jk_0 \underbrace{\oint_{S_{ABC}}(\hat{n}\times\mathbf{T})\cdot(\hat{n}\times\mathbf{E}^{sca})\,dS}_{a_3(\mathbf{T},\mathbf{E}^{sca})} \\ &= \underbrace{k_0^2\int_{\Omega_d}\mathbf{T}\cdot(\bar{\bar{\varepsilon}}_r-\bar{\bar{I}})\mathbf{E}^{inc}\,dV}_{F(\mathbf{T})} \end{aligned} $$

简记为：

$$ \boxed{a(\mathbf{T}, \mathbf{E}^{sca}) = F(\mathbf{T}), \quad \forall\,\mathbf{T} \in H(\text{curl};\Omega)} $$

其中：

$$ a(\mathbf{T}, \mathbf{E}) = a_1(\mathbf{T},\mathbf{E}) - k_0^2 a_2(\mathbf{T},\mathbf{E}) + a_3(\mathbf{T},\mathbf{E}) $$

注意 $a_3$ 已经包含了 ABC 边界条件（自然边界条件），它来自弱形式推导中的面积分项。

二、全局编号下的场逼近

2.1 全局棱边编号与全局基函数

设网格共有 $N_E$ 条全局棱边，编号为 $I = 1, 2, \ldots, N_E$。

Hcurl on wuqq 的 Blog