ABC on wuqq 的 Blog

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} $$

其中：