计算电磁 on wuqq 的 Blog

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 方程”，也会学到一整套局部矩阵、全局组装和边界条件处理方法；但如果只停留在这一层，就很难真正回答下面这些问题：

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

现代观点的统一

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

0047 - 二维各向异性介质电磁散射有限元：参数对应、弱形式与二维三维误差分析

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

这篇笔记在做什么

这一篇整理的是二维各向异性介质散射问题。原始会话里先从更一般的 Maxwell 对应关系出发，随后又收束到“无磁、仅电各向异性”的情形；这里我把两部分都保留下来，方便你同时看到一般框架和具体约束后的推导结果。

文章后半部分还保留了二维与三维在离散误差、污染误差、自由度规模和求解策略上的系统对比。

二维各向异性介质电磁散射问题的有限元方法

从 Maxwell 方程到 Helmholtz 方程 · 完整参数对应 · 二维/三维误差分析

一、从 Maxwell 方程推导二维标量方程

1.1 问题设置

二维问题：所有场量与 $z$ 无关（$\partial/\partial z = 0$）。

各向异性介质参数（对角张量）：

$$ \bar{\bar{\varepsilon}}_r = \begin{pmatrix}\varepsilon_{xx}&0&0\\0&\varepsilon_{yy}&0\\0&0&\varepsilon_{zz}\end{pmatrix}, \quad \bar{\bar{\mu}}_r = \begin{pmatrix}\mu_{xx}&0&0\\0&\mu_{yy}&0\\0&0&\mu_{zz}\end{pmatrix} $$

时间约定：$e^{-j\omega t}$

Maxwell 方程：

$$ \nabla\times\mathbf{E} = j\omega\mu_0\bar{\bar{\mu}}_r\mathbf{H} $$$$ \nabla\times\mathbf{H} = -j\omega\varepsilon_0\bar{\bar{\varepsilon}}_r\mathbf{E} $$

1.2 TM 极化（$E_z$ 极化）

场分量：$\mathbf{E} = E_z\hat{z}$，$\mathbf{H} = H_x\hat{x} + H_y\hat{y}$

从 Faraday 定律：

$$ \nabla\times(E_z\hat{z}) = \frac{\partial E_z}{\partial y}\hat{x} - \frac{\partial E_z}{\partial x}\hat{y} = j\omega\mu_0(\mu_{xx}H_x\hat{x}+\mu_{yy}H_y\hat{y}) $$$$ H_x = \frac{1}{j\omega\mu_0\mu_{xx}}\frac{\partial E_z}{\partial y}, \quad H_y = \frac{-1}{j\omega\mu_0\mu_{yy}}\frac{\partial E_z}{\partial x} $$

从 Ampere 定律：

0044 - 三维各向异性介质散射有限元完整推导：四面体棱边元、ABC 边界与矩阵组装

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

这篇笔记在做什么

这一篇对应原始会话里最完整的一条三维电磁有限元推导链。我保留了四面体一阶 Nedelec 棱边元、$e^{-j\omega t}$ 约定、ABC 面矩阵、全局与局部映射、质量矩阵和代码框架的细致展开。

如果你关心的是“从 Maxwell 方程一路走到可实现的单元矩阵和组装算法”，这一篇就是主干正文。

三维各向异性介质散射问题有限元方法

——四面体棱边元 · ABC边界 · $e^{-j\omega t}$约定 · 完整推导

一、物理问题与几何设置

各向异性无磁电介质立方体，被平面波照射：

 E^inc (平面波)
↓ ↓ ↓
╔══════════════════════════╗
║ Ω₀ (自由空间, ε_r=I) ║ ← S_ABC (外边界, 球/长方体)
║ ║
║ ┌──────────┐ ║
║ │ Ω_d │ ║
║ │ 各向异性 │ ║
║ │ 介质体 │ ║
║ │ ε̄̄_r(张量)│ ║
║ └──────────┘ ║
║ S_d (介质面) ║
║ ║
╚══════════════════════════╝
μ_r = 1 (全空间无磁)

材料参数——各向异性介电张量：

0043 - 三维电磁散射有限元入门：从 PEC 立方体到弱形式、棱边元与全局编号

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

这篇笔记在做什么

这篇笔记从一份 Playground 导出会话整理而来。我保留了原始推导里的细节，重点把最简单的三维 PEC 立方体散射问题整理成一个可连续阅读的入门版本。

文章会按下面这条线展开：Maxwell 方程、散射场分解、ABC 截断、弱形式、棱边元、六面体单元、局部到全局编号、矩阵组装。

三维电磁场散射问题的有限元方法——以立方体为例

一、物理问题描述

考虑一个边长为 $a$ 的理想导体(PEC)立方体，置于自由空间中，被一个平面波照射。我们要求解散射场。

 入射平面波 E^inc
↓ ↓ ↓
___________
/ /|
/ / |
/___________/ | PEC 立方体
| | |
| | /
| | /
|___________|/
求：散射场 E^sca

二、理论推导

2.1 从 Maxwell 方程出发

时谐场（$e^{j\omega t}$ 约定）下的 Maxwell 方程：

$$ \nabla \times \mathbf{E} = -j\omega\mu\mathbf{H} $$$$ \nabla \times \mathbf{H} = j\omega\varepsilon\mathbf{E} + \mathbf{J} $$

消去 $\mathbf{H}$，得到矢量波动方程：

$$ \nabla \times \left(\frac{1}{\mu_r} \nabla \times \mathbf{E}\right) - k_0^2 \varepsilon_r \mathbf{E} = -jk_0 Z_0 \mathbf{J} $$

其中：

计算电磁 on wuqq 的 Blog

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

总结

目录

我的整理笔记

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

一套自洽的知识手册

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

传统观点的局限

现代观点的统一

0047 - 二维各向异性介质电磁散射有限元：参数对应、弱形式与二维三维误差分析

这篇笔记在做什么

二维各向异性介质电磁散射问题的有限元方法

从 Maxwell 方程到 Helmholtz 方程 · 完整参数对应 · 二维/三维误差分析

一、从 Maxwell 方程推导二维标量方程

1.1 问题设置

1.2 TM 极化（$E_z$ 极化）

0044 - 三维各向异性介质散射有限元完整推导：四面体棱边元、ABC 边界与矩阵组装

这篇笔记在做什么

三维各向异性介质散射问题有限元方法

——四面体棱边元 · ABC边界 · $e^{-j\omega t}$约定 · 完整推导

一、物理问题与几何设置

0043 - 三维电磁散射有限元入门：从 PEC 立方体到弱形式、棱边元与全局编号

这篇笔记在做什么

三维电磁场散射问题的有限元方法——以立方体为例

一、物理问题描述

二、理论推导

2.1 从 Maxwell 方程出发