总结

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

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

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

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

目录

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

我的整理笔记

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

一套自洽的知识手册

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

传统观点的局限

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

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

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

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

现代观点的统一

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

$$ \text{物理场} = \text{微分形式} \xrightarrow{\text{de Rham 复形}} \text{精确序列} \xrightarrow{\text{Whitney 映射}} \text{离散精确序列} \xrightarrow{\text{Galerkin}} \text{矩阵} $$

在这个框架下：

  • 节点元、棱边元、面元、体元不是彼此独立发明的元素族，而是同一代数结构在不同阶上的离散实现。
  • 棱边元没有伪解不是偶然现象，而是因为它保持了连续问题的拓扑结构。
  • 稳定性、收敛性和误差估计不再需要逐案猜测，而可以从统一的 Hilbert 复形理论中推出。

第一章：微分形式，电磁场的自然语言

1.1 为什么矢量分析不够好？

传统电磁学用矢量场描述场量，但“同样是矢量”的记法会掩盖物理本质上的差异。

物理量传统写法物理本质
电场 $\mathbf{E}$矢量场沿路径积分有意义：$\int_C \mathbf{E}\cdot d\mathbf{l}$
磁通密度 $\mathbf{B}$矢量场沿面积分有意义：$\int_S \mathbf{B}\cdot d\mathbf{S}$
电荷密度 $\rho$标量场沿体积分有意义：$\int_V \rho\,dV$

问题不在于矢量分析“错了”，而在于它把不同积分对象都写成了看起来类似的量。微分形式用“阶数”把这些差异直接编码进对象本身。

1.2 微分形式的定义

在三维空间 $\mathbb{R}^3$ 中：

0-形式：在点上求值

$$ f \in \Omega^0(\mathbb{R}^3) $$$$ \text{求值：} f(p) $$

1-形式：沿曲线积分

$$ \omega = E_x\,dx + E_y\,dy + E_z\,dz \in \Omega^1(\mathbb{R}^3) $$$$ \int_C \omega = \int_C (E_x\,dx + E_y\,dy + E_z\,dz) $$

电场 $\mathbf{E}$ 自然对应 1-形式，因为电压本质上是沿曲线的积分。

2-形式：沿曲面积分

$$ \sigma = B_x\,dy\wedge dz + B_y\,dz\wedge dx + B_z\,dx\wedge dy \in \Omega^2(\mathbb{R}^3) $$$$ \int_S \sigma = \int_S \mathbf{B}\cdot d\mathbf{S} $$

磁通密度 $\mathbf{B}$ 自然对应 2-形式，因为磁通是面积分。

3-形式：沿体积积分

$$ \tau = \rho\,dx\wedge dy\wedge dz \in \Omega^3(\mathbb{R}^3) $$$$ \int_V \tau = \int_V \rho\,dV $$

电荷密度 $\rho$ 自然对应 3-形式。

1.3 楔积与外导数

楔积满足：

$$ dx \wedge dy = -dy \wedge dx $$$$ dx \wedge dx = 0 $$

外导数 $d: \Omega^k \to \Omega^{k+1}$ 统一了梯度、旋度和散度：

$$ df = \frac{\partial f}{\partial x}dx + \frac{\partial f}{\partial y}dy + \frac{\partial f}{\partial z}dz $$$$ d(E_x dx + E_y dy + E_z dz) = \left(\frac{\partial E_z}{\partial y}-\frac{\partial E_y}{\partial z}\right)dy\wedge dz + \cdots $$$$ d(B_x\,dy\wedge dz + B_y\,dz\wedge dx + B_z\,dx\wedge dy) = \left(\frac{\partial B_x}{\partial x}+\frac{\partial B_y}{\partial y}+\frac{\partial B_z}{\partial z}\right)dx\wedge dy\wedge dz $$

统一表格如下：

矢量分析微分形式外导数映射
$\nabla f$$df$$d: \Omega^0 \to \Omega^1$
$\nabla\times\mathbf{E}$$dE$$d: \Omega^1 \to \Omega^2$
$\nabla\cdot\mathbf{B}$$dB$$d: \Omega^2 \to \Omega^3$

最关键的结构性质只有一句话：

$$ \boxed{d \circ d = 0} $$

它统一了：

$$ \nabla\times(\nabla f) = 0 $$$$ \nabla\cdot(\nabla\times\mathbf{E}) = 0 $$

1.4 Maxwell 方程的微分形式版本

传统矢量形式：

$$ \nabla\times\mathbf{E} = -\frac{\partial\mathbf{B}}{\partial t}, \quad \nabla\cdot\mathbf{B} = 0 $$$$ \nabla\times\mathbf{H} = \mathbf{J} + \frac{\partial\mathbf{D}}{\partial t}, \quad \nabla\cdot\mathbf{D} = \rho $$

微分形式写成：

$$ E \in \Omega^1, \quad H \in \Omega^1 $$$$ B \in \Omega^2, \quad D \in \Omega^2 $$$$ J \in \Omega^2, \quad \rho \in \Omega^3 $$

于是 Maxwell 方程变成：

$$ \boxed{dE = -\frac{\partial B}{\partial t}, \quad dB = 0} $$$$ \boxed{dH = J + \frac{\partial D}{\partial t}, \quad dD = \rho} $$

其中 $dB = 0$ 意味着 $B$ 是闭形式。在单连通区域中，由 Poincare 引理，存在 $A \in \Omega^1$ 使得：

$$ B = dA $$

这就是矢量势的自然来源。

1.5 Hodge 星算子与材料关系

在微分形式语言中，$E$ 是 1-形式，$D$ 是 2-形式，所以物质关系不再是“同类型对象乘一个系数”，而要借助 Hodge 星算子：

$$ \star: \Omega^k \to \Omega^{3-k} $$

材料关系写成：

$$ D = \star_\varepsilon E, \quad B = \star_\mu H $$

对于各向异性介质，这个结构可以写成：

$$ \star_\varepsilon (E_x\,dx + E_y\,dy + E_z\,dz) = \varepsilon_{xx}E_x\,dy\wedge dz + \varepsilon_{yy}E_y\,dz\wedge dx + \varepsilon_{zz}E_z\,dx\wedge dy $$

这一步非常关键：各向异性不是“方程里某个系数变了”，而是 Hodge 星算子的结构变了。于是：

  • 外导数 $d$ 对应刚度部分，体现拓扑结构；
  • Hodge 星 $\star_\varepsilon$ 对应质量部分，体现材料和度量结构。

第二章：de Rham 复形，有限元的拓扑骨架

2.1 连续 de Rham 复形

外导数把不同阶微分形式串成一个复形：

$$ \boxed{ 0 \xrightarrow{} \Omega^0 \xrightarrow{d_0} \Omega^1 \xrightarrow{d_1} \Omega^2 \xrightarrow{d_2} \Omega^3 \xrightarrow{} 0 } $$

精确性意味着：

$$ \text{Im}(d_k) = \text{Ker}(d_{k+1}) $$

用 Sobolev 空间写成：

$$ H^1(\Omega) \xrightarrow{\nabla} H(\text{curl};\Omega) \xrightarrow{\nabla\times} H(\text{div};\Omega) \xrightarrow{\nabla\cdot} L^2(\Omega) $$
空间微分形式矢量分析物理量连续性
$H^1$0-形式标量场电位 $\phi$函数值连续
$H(\text{curl})$1-形式矢量场电场 $\mathbf{E}$切向连续
$H(\text{div})$2-形式矢量场磁通密度 $\mathbf{B}$法向连续
$L^2$3-形式标量场电荷密度 $\rho$无连续性要求

2.2 为什么精确序列重要？

连续物理定律本身就是精确序列的表现。

例如无磁荷：

$$ \nabla\cdot\mathbf{B} = 0 $$

意味着：

$$ \mathbf{B} \in \text{Ker}(\nabla\cdot) = \text{Im}(\nabla\times) $$

所以存在 $\mathbf{A}$ 使得：

$$ \mathbf{B} = \nabla\times\mathbf{A} $$

同样，若 $\nabla\times\mathbf{E} = 0$，则在适当区域内：

$$ \mathbf{E} = \nabla\phi $$

2.3 伪解的根源

如果离散空间不保持这个精确关系，就会出现离散核空间比真实核空间更大的情况。多出来的那部分元素在连续物理中没有对应对象，但在离散方程中却可能被当成“解”，这就是伪解。

更准确地说：

$$ \text{伪解} = \text{离散核空间中不该出现的多余元素} $$

如果离散层面仍然满足

$$ \text{Ker}(d_h^k) = \text{Im}(d_h^{k-1}) $$

则伪解被结构性排除。

2.4 Sobolev 空间定义

$$ H^1(\Omega) = \{u \in L^2(\Omega) : \nabla u \in [L^2(\Omega)]^3\} $$$$ \|u\|_{H^1} = \left(\int_\Omega |u|^2 + |\nabla u|^2\,dV\right)^{1/2} $$$$ H(\text{curl};\Omega) = \{\mathbf{u} \in [L^2(\Omega)]^3 : \nabla\times\mathbf{u} \in [L^2(\Omega)]^3\} $$$$ \|\mathbf{u}\|_{H(\text{curl})} = \left(\int_\Omega |\mathbf{u}|^2 + |\nabla\times\mathbf{u}|^2\,dV\right)^{1/2} $$$$ H(\text{div};\Omega) = \{\mathbf{u} \in [L^2(\Omega)]^3 : \nabla\cdot\mathbf{u} \in L^2(\Omega)\} $$

嵌入关系为：

$$ H^1 \hookrightarrow H(\text{curl}) \hookrightarrow L^2 $$

但 $H^1$ 和 $H(\text{curl})$ 并不相同。$H(\text{curl})$ 只要求切向连续，而不要求整体函数值连续。

第三章：离散 de Rham 复形与 Whitney 形式

3.1 离散化的核心原则

一个好的离散化，首先要保持精确序列。也就是说，要构造离散空间 $V_h^k$ 使得下面的图交换：

$$ \begin{array}{ccccccc} H^1 & \xrightarrow{\nabla} & H(\text{curl}) & \xrightarrow{\nabla\times} & H(\text{div}) & \xrightarrow{\nabla\cdot} & L^2 \\ \downarrow\pi_h^0 & & \downarrow\pi_h^1 & & \downarrow\pi_h^2 & & \downarrow\pi_h^3 \\ V_h^0 & \xrightarrow{\nabla} & V_h^1 & \xrightarrow{\nabla\times} & V_h^2 & \xrightarrow{\nabla\cdot} & V_h^3 \end{array} $$

交换性意味着：

$$ \pi_h^{k+1} \circ d = d \circ \pi_h^k $$

“先微分再投影”和“先投影再微分”必须一致。

3.2 Whitney 形式

0-形式离散：节点元

自由度位置是节点，基函数是重心坐标：

$$ W_i^0 = \lambda_i $$

自由度是节点值：

$$ \sigma_i^0(f) = f(v_i) $$

1-形式离散：棱边元

自由度位置是棱边，基函数是 Whitney 1-形式：

$$ \boxed{W_{ij}^1 = \lambda_i\,d\lambda_j - \lambda_j\,d\lambda_i} $$

用矢量形式写就是：

$$ \mathbf{W}_{ij}^1 = \lambda_i\nabla\lambda_j - \lambda_j\nabla\lambda_i $$

这正是一阶 Nedelec 棱边元。

自由度是沿棱边的线积分：

$$ \sigma_{ij}^1(\omega) = \int_{e_{ij}} \omega $$

离散函数只要求切向连续。

2-形式离散：面元

自由度位置是面，基函数是 Whitney 2-形式：

$$ W_{ijk}^2 = 2(\lambda_i\,d\lambda_j\wedge d\lambda_k + \lambda_j\,d\lambda_k\wedge d\lambda_i + \lambda_k\,d\lambda_i\wedge d\lambda_j) $$

对应的矢量形式可以理解为 Raviart-Thomas / 面元结构，法向连续。

3-形式离散：体元

自由度位置是体，通常只要求体积分意义下的离散，不要求跨单元连续。

3.3 Whitney 形式保持精确序列

Whitney 复形在适当条件下保持离散精确性：

$$ V_h^0 \xrightarrow{\nabla} V_h^1 \xrightarrow{\nabla\times} V_h^2 \xrightarrow{\nabla\cdot} V_h^3 $$

例如：

$$ \nabla\times(\lambda_i\nabla\lambda_j - \lambda_j\nabla\lambda_i) = 2\nabla\lambda_i\times\nabla\lambda_j $$

从而离散梯度落在离散 1-形式空间中，离散旋度又落在离散 2-形式空间中。

3.4 统一表格

$$ \boxed{ \begin{array}{c|c|c|c|c} k\text{-形式} & V_h^k & \text{DOF位置} & \text{基函数} & \text{连续性} \\\hline 0 & \text{节点元} & \text{节点} & \lambda_i & C^0 \\ 1 & \text{棱边元} & \text{棱边} & \lambda_i\nabla\lambda_j-\lambda_j\nabla\lambda_i & \text{切向连续} \\ 2 & \text{面元} & \text{面} & \lambda_i\nabla\lambda_j\times\nabla\lambda_k+\cdots & \text{法向连续} \\ 3 & \text{体元} & \text{体} & \text{常数或不连续} & \text{无连续性} \end{array} } $$

这说明节点元、棱边元、面元、体元不是彼此独立的发明，而是同一个离散复形在不同阶上的四个层次。

第四章：变分框架与 Hodge-Laplace 问题

4.1 内积结构

在 $k$-形式空间上，自然的 $L^2$ 型内积写成：

$$ (\omega, \eta)_{L^2} = \int_\Omega \omega \wedge \star\eta $$

对于电磁场对应的 1-形式，这就退化成熟悉的形式：

$$ (\mathbf{u}, \mathbf{v})_{L^2} = \int_\Omega \mathbf{u}\cdot\mathbf{v}\,dV $$$$ (\nabla\times\mathbf{u}, \nabla\times\mathbf{v})_{L^2} = \int_\Omega (\nabla\times\mathbf{u})\cdot(\nabla\times\mathbf{v})\,dV $$

4.2 Hodge-Laplace 算子

$k$-形式上的 Hodge-Laplace 算子为：

$$ \Delta_k = d_{k-1}\delta_k + \delta_{k+1}d_k $$

其中 $\delta = \star d \star$ 是余外导数。

对于 1-形式，可把它理解为：

$$ \Delta_1\mathbf{E} = \nabla\times(\nabla\times\mathbf{E}) - \nabla(\nabla\cdot\mathbf{E}) $$

在适当约束下，Maxwell 波动方程正可以视为 Hodge-Laplace 问题的一个具体实现。

4.3 电磁问题的现代表述

设 $E, T \in H(\text{curl};\Omega)$，则双线性形式可以写成：

$$ a(T, E) = \underbrace{(dT, dE)_{L^2}}_{\text{拓扑项}} - k_0^2\underbrace{(T, \star_\varepsilon E)_{L^2}}_{\text{材料项}} + \underbrace{\text{边界项}}_{\text{ABC/PML}} $$

对应关系可以总结为：

项数学来源物理来源依赖拓扑？依赖材料？
$(dT,dE)$外导数 $d$传播结构是否
$(T,\star_\varepsilon E)$Hodge 星介质响应否是
边界项边界上的 $\star$ 与迹算子ABC / PML / DtN否一般与边界模型有关

这正是“刚度矩阵主要对应拓扑结构，质量矩阵体现材料结构”的现代解释。

4.4 Galerkin 离散

取 $E_h, T_h \in V_h^1$，并写成：

$$ E_h = \sum_{I=1}^{N_E} E_I\,W_I^1 $$

则矩阵方程为：

$$ [A]\{E\} = \{b\} $$

其中

$$ A_{JI} = \underbrace{(dW_J^1, dW_I^1)_{L^2}}_{S_{JI}} - k_0^2\underbrace{(W_J^1, \star_\varepsilon W_I^1)_{L^2}}_{T_{JI}} + \underbrace{\text{边界项}}_{B_{JI}} $$

第五章：有限元外微积分 FEEC

5.1 Arnold-Falk-Winther 框架

有限元外微积分（Finite Element Exterior Calculus, FEEC）由 Arnold、Falk 和 Winther 在 2000 年代系统建立。它的核心思想不是为每个 PDE 单独发明有限元，而是从离散复形出发，自动得到稳定的离散方案。

传统方法更像“逐案试验”：

Maxwell 方程 -> 试节点元 -> 发现伪解 -> 发明棱边元
Stokes 方程  -> 试不同速度-压力配对 -> 找 inf-sup 稳定元
弹性问题    -> 试位移元 -> 发现 locking -> 发明修正方法

FEEC 则是：

识别 PDE 对应的复形
-> 选择保持交换图的离散子复形
-> 稳定性与误差分析统一完成

5.2 抽象 Hilbert 复形

一般 Hilbert 复形写作：

$$ \boxed{ (W, d): \quad W^0 \xrightarrow{d^0} W^1 \xrightarrow{d^1} W^2 \xrightarrow{d^2} \cdots \xrightarrow{d^{n-1}} W^n } $$

满足：

$$ d^{k+1}\circ d^k = 0 $$

对三维电磁学而言，可以取：

$$ W^0 = H^1(\Omega), \quad W^1 = H(\text{curl};\Omega), \quad W^2 = H(\text{div};\Omega), \quad W^3 = L^2(\Omega) $$$$ d^0 = \nabla, \quad d^1 = \nabla\times, \quad d^2 = \nabla\cdot $$

5.3 抽象 Hodge-Laplace 问题

在 Hilbert 复形上，$k$-阶 Hodge-Laplace 问题可以写成混合形式：

$$ (\sigma, \tau) - (d\tau, u) = 0, \quad \forall\,\tau \in W^{k-1} $$$$ (d\sigma, v) + (du, dv) + (v, p) = (f, v), \quad \forall\,v \in W^k $$$$ (u, q) = 0, \quad \forall\,q \in \mathfrak{H}^k $$

这里 $\mathfrak{H}^k$ 是调和形式空间。

这个混合框架的意义在于：它把核空间和调和空间明确地纳入方程结构中，所以不会在低频极限、零特征值附近或拓扑非平凡区域里失去控制。

5.4 离散子复形与共链投影

选取有限维子空间 $V_h^k \subset W^k$，要求：

$$ V_h^0 \xrightarrow{d^0} V_h^1 \xrightarrow{d^1} V_h^2 \xrightarrow{d^2} V_h^3 $$

仍构成复形。

再引入投影：

$$ \pi_h^k: W^k \to V_h^k $$

要求满足：

$$ \boxed{\pi_h^{k+1}\circ d^k = d^k \circ \pi_h^k} $$

这就是交换图性质，也是 FEEC 中最重要的结构条件。

5.5 FEEC 的核心定理

在满足交换性、有界性和逼近性的前提下，离散 Hodge-Laplace 问题满足类似最佳逼近的误差估计：

$$ \|\sigma - \sigma_h\| + \|u-u_h\| + \|p-p_h\| \leq C \inf_{(\tau_h, v_h, q_h)} \left(\|\sigma-\tau_h\| + \|u-v_h\| + \|p-q_h\|\right) $$

对电磁 1-形式问题，这直接给出：

$$ \|\mathbf{E} - \mathbf{E}_h\|_{H(\text{curl})} \leq C\,\inf_{\mathbf{v}_h\in V_h^1}\|\mathbf{E}-\mathbf{v}_h\|_{H(\text{curl})} $$

再结合 Whitney / Nedelec 插值误差估计，就可以推出标准收敛结论，而不必为每个离散方案重新搭一遍理论。

第六章：Helmholtz 分解与数值困难

6.1 连续 Helmholtz 分解

对任意 $\mathbf{u} \in H(\text{curl};\Omega)$，存在分解：

$$ \boxed{\mathbf{u} = \nabla\phi + \mathbf{z}} $$

其中 $\nabla\phi$ 是无旋部分，$\mathbf{z}$ 是其补空间中的旋度部分。

对电场来说，这意味着：

$$ \mathbf{E} = -\nabla\phi + \mathbf{z} $$

静电分量和辐射分量在函数空间层面是分开的。

6.2 离散 Helmholtz 分解

在棱边元空间中，也有类似结构：

$$ V_h^1 = \nabla V_h^0 \oplus Z_h $$

其中 $\nabla V_h^0$ 对应离散梯度场，$Z_h$ 对应其离散正交补。

由于交换图成立，这个分解和连续问题是对齐的，而不是“拍脑袋分块”。

6.3 低频崩溃

时谐问题写成：

$$ (\nabla\times\mathbf{T}, \nabla\times\mathbf{E}) - k_0^2(\mathbf{T}, \varepsilon\mathbf{E}) = (\mathbf{T}, \mathbf{f}) $$

当 $k_0 \to 0$ 时，质量项衰减，梯度空间上的控制只剩下很小的 $k_0^2$ 级项，于是系统的条件数迅速恶化。

大致上：

$$ \kappa([A]) \sim \frac{O(1)}{k_0^2} $$

这就是低频崩溃的结构原因。

6.4 Tree-Cotree 分解

在网格图论结构中，可以把棱边分成树与余树：

$$ \text{棱边集} = \text{tree} \cup \text{cotree} $$

它们分别对应离散梯度部分和旋度部分。这个分解的本质仍然是离散复形的维数分解，而不是纯粹的程序技巧。

第七章：误差分析的现代观点

7.1 抽象误差理论

对不定问题，经典 Cea 引理需要用 Schatz 类型论证推广。若双线性形式满足连续性和 Gårding 不等式，并且离散空间满足适当逼近性质，则可得到：

$$ \|\mathbf{E}-\mathbf{E}_h\|_V \leq C(1+\Gamma_h)\inf_{\mathbf{v}_h\in V_h}\|\mathbf{E}-\mathbf{v}_h\|_V $$

这说明离散解依然是一个稳定的近最佳逼近。

7.2 Whitney 1-形式的插值误差

对一阶棱边元，设 $\mathbf{E}$ 足够光滑，则有：

$$ \|\mathbf{E}-\pi_h^1\mathbf{E}\|_{L^2(T)} \leq C\,h_T^{s+1}|\mathbf{E}|_{H^{s+1}(T)} $$$$ \|\nabla\times(\mathbf{E}-\pi_h^1\mathbf{E})\|_{L^2(T)} \leq C\,h_T^{s+1}|\nabla\times\mathbf{E}|_{H^{s+1}(T)} $$

于是：

$$ \|\mathbf{E}-\pi_h^1\mathbf{E}\|_{H(\text{curl},T)} \leq C\,h_T\left(|\mathbf{E}|_{H^1(T)} + |\nabla\times\mathbf{E}|_{H^1(T)}\right) $$

7.3 与节点元插值的对比

节点元在 $H^1$ 空间里的误差估计是经典的，而棱边元在 $H(\text{curl})$ 空间中还额外涉及 $\nabla\times\mathbf{E}$ 的正则性。

量节点元棱边元
自然空间$H^1$$H(\text{curl})$
连续性要求值连续切向连续
额外结构无需要控制旋度
对材料界面更自然？否是

材料界面上通常只满足：

$$ \hat{n}\times(\mathbf{E}_1-\mathbf{E}_2)=0 $$

但法向分量不一定连续，所以 $\mathbf{E}$ 往往不属于 $[H^1]^3$，而只属于更自然的 $H(\text{curl})$。这正是棱边元优于节点元的函数空间原因。

7.4 对偶论证与 $L^2$ 收敛

在光滑凸域中，Aubin-Nitsche 型对偶论证可能给出更高的 $L^2$ 收敛阶；但在一般电磁问题里，边角、非凸域和材料界面会破坏所需的正则性，于是常见情形下棱边元的 $L^2$ 误差并不会像标量椭圆问题那样自然获得额外一阶超收敛。

这也是电磁问题和经典 Poisson 问题在误差分析上最容易被混淆、但实际差别很大的地方。

第八章：散射问题的函数空间框架

8.1 无界域与辐射条件

散射问题定义在无界域上，函数空间不仅要包含局部可积与旋度可积性，还要嵌入辐射条件。

Silver-Muller 条件写成：

$$ \lim_{r\to\infty}(\nabla\times\mathbf{E}^{sca})\times\hat{r} - jk_0\mathbf{E}^{sca} = o(r^{-1}) $$

它的作用不是“再加一个边界方程”，而是从所有形式解里选出唯一的外向传播解。

8.2 截断、ABC 与 DtN

把无界域截断成有界计算域后，需要在外边界上引入等价条件。精确的边界条件是 DtN 映射：

$$ \hat{n}\times(\nabla\times\mathbf{E}^{sca}) = \mathcal{T}(\hat{n}\times\mathbf{E}^{sca}\times\hat{n}) $$

但它是非局部算子，实现代价高。

一阶 ABC：

$$ \hat{n}\times(\nabla\times\mathbf{E}^{sca}) = jk_0\hat{n}\times(\hat{n}\times\mathbf{E}^{sca}) $$

是对精确 DtN 的低阶近似。

8.3 PML 的数学本质

PML 本质上是复坐标拉伸。写成复坐标变换后，外行波在 PML 区域里指数衰减。对 FEEC 而言，这件事最自然的解释是：

  • 外导数 $d$ 不变，因为拓扑不变；
  • Hodge 星被修改，因为度量结构被复拉伸改变。

也就是说，PML 在现代框架里更像是“改 Hodge 星”，而不是“发明一个额外 PDE”。

第九章：高阶方法与 hp 自适应

9.1 高阶 Whitney 形式

一阶棱边元属于最低阶空间。更高阶时，可以使用 Nedelec 第一类或第二类高阶空间，它们依然保持离散复形与交换图结构。

常见的理解方式是：

  • 第一类更紧凑；
  • 第二类更富多项式自由度；
  • 两者都保持精确序列，因此都不会引入拓扑错误。

9.2 hp 自适应

对光滑区域，提高多项式阶数 $p$ 往往更高效；对奇异点、尖角和材料界面附近，局部减小网格尺度 $h$ 更重要。

因此理想策略不是只做 $h$ 加密，也不是只做 $p$ 提升，而是：

$$ \text{奇异处做 } h\text{-加密，光滑处做 } p\text{-提升} $$

在解析解或近解析解区域，$hp$ 方法可达到接近指数型的误差衰减；而在奇异性主导区域，合理的局部加密仍然是关键。

第十章：从理论到计算的桥梁

10.1 理论与代码结构的对应

现代数学理论                          代码实现
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
de Rham 复形                         节点 -> 棱边 -> 面 -> 体 的拓扑结构
Whitney 1-形式                       Nedelec 棱边基函数
外导数 d                             curl 项、拓扑项、刚度结构
Hodge 星 ★_ε                         材料项、质量矩阵、各向异性
交换图投影 π_h^k                     全局-局部映射与 DOF 一致性
Galerkin 投影                        单元矩阵 -> 全局组装
Helmholtz 分解                        tree-cotree / 低频稳定化

10.2 几个“为什么”的终极回答

为什么用棱边元而不是节点元？

因为电场是 1-形式，天然属于 $H(\text{curl})$，需要的是切向连续。棱边元正好实现这个要求；节点元给了过强的值连续约束，破坏了正确的离散核空间结构。

为什么刚度矩阵不依赖材料，而质量矩阵依赖材料？

因为刚度矩阵来自外导数 $d$，它体现的是拓扑和微分结构；质量矩阵来自 Hodge 星 $\star_\varepsilon$，它体现的是度量和材料结构。

为什么组装时需要方向符号？

因为 1-形式的自由度本质上是沿棱边的积分，而积分对象有定向。全局方向和局部方向不一致时，符号修正是微分形式定向性在离散层面的直接体现。

为什么 ABC 有误差，而 PML 在连续层面可做到无反射？

因为 ABC 是对精确 DtN 映射的截断近似，PML 则是复坐标拉伸后的等价模型。前者天然带有截断误差，后者在连续层面可以做到精确匹配。

10.3 为什么高频问题更难？

因为 Helmholtz / Maxwell 时谐问题在高频下同时叠加了三类困难：

  • 双线性形式的不定性更强，稳定常数随波数变坏。
  • 相位误差会在传播过程中累积，形成污染误差，而不只是局部插值误差。
  • 每波长需要足够多的自由度，频率越高，网格越细、系统越大、求解越难。

所以“高频更难”不只是因为未知量更多，而是因为连续问题的稳定性、离散问题的相位精度、线性系统的条件数和求解成本都在一起恶化。

结语

如果把整篇笔记再压缩成一句话，那就是：

$$ \boxed{ \text{三维电磁散射有限元的核心，不是某个单独公式，而是拓扑结构、度量结构与离散结构三者的统一。} } $$

在这个统一视角下，微分形式提供语言，de Rham 复形提供骨架，Whitney 形式提供离散对象，FEEC 提供稳定性与误差理论，最终才落到我们熟悉的棱边元、局部矩阵、全局组装和数值求解。