什么是低频崩溃（Low‑Frequency Breakdown）问题及其解决
电磁学的发展始终与实际工程需求紧密绑定。早期，远距离通信的迫切需求推动了海底电缆、电报系统的研发，但很快发现，铜质电缆的有限电导率会导致信号在低频传输中严重畸变，而当时仅靠经验化的电报方程无法精准描述这一现象。随着实验数据的积累与理论研究的深入，麦克斯韦方程组成为解决这类问题的核心工具。然而，在低频、准静态乃至直流极限场景下，这套经典方程组却遭遇了低频崩溃这一关键难题：数值求解过程中方程奇异、解不唯一，无法得到稳定可靠的结果。低频崩溃问题不仅制约了低频电磁仿真、电路 - 电磁场联合分析、静磁学建模等工程应用，也暴露了麦克斯韦方程组在低频极限下的理论缺陷。
麦克斯韦方程组与低频极限特征
1.1 时频域麦克斯韦方程组
其中，前两式为旋度方程，描述电磁场的涡旋耦合；后两式为散度方程，描述场源的分布特征。
1.2 低频极限下的方程退化
低频场景对应角频率ω->0，此时ω²量级远小于ω量级，频域旋度方程中含ω²的项可近似忽略。将本构关系代入旋度方程，消去B与H，可推导出电场满足的波动方程：
低频极限下，方程退化为：
该退化方程是低频崩溃问题的直接载体，其数学特性决定了低频场景下的求解困境。
低频崩溃问题的本质：旋度算子的零空间
2.1 旋度算子的核心恒等式
低频崩溃的根源在于旋度算子存在零空间，这一特性可由矢量微积分基本恒等式直接推导，任意标量场梯度的旋度恒为零，即：
其中ψ为任意光滑标量场。该恒等式表明：无旋场（标量场梯度）无法通过旋度算子区分，旋度方程无法约束电场的无旋分量。
2.2 解的不唯一性：低频崩溃的数学表现
对于低频退化方程
假设E0是方程的一组解，根据旋度恒等式，对任意标量场ψ，
也必然满足方程：
这意味着低频退化方程存在无穷多组解，解不唯一。
从数值求解角度看，通过有限差分、有限元等离散方法将连续方程转化为矩阵方程时，方程的零空间会导致系数矩阵奇异、行列式为零、矩阵不可逆，无法通过常规数值方法求解，直接造成仿真发散、结果失真，即低频崩溃。
2.3 物理层面的直观理解
从物理本质看，电场可分解为两部分：涡旋电场（由磁场变化感应产生，无散有旋）与静电场（由电荷分布产生，无旋有散）。
高频场景：电磁波传播主导，涡旋电场占主导，旋度方程可有效约束电场分布，解唯一稳定；
低频 / 准静态场景：静电场分量占主导，而旋度方程仅能约束涡旋电场，无法约束无旋的静电场分量，导致电场分布不确定，对应数学上的解不唯一。
低频崩溃的解决方案：A-Φ 势法
矢量势A：描述磁场与电场的涡旋分量（对应 
）；标势Φ：描述电场的无旋静电分量（对应-▽Φ）。
该分解完全覆盖电场、磁场的全部物理信息，无冗余、无遗漏，从根本上解决了旋度方程无法约束无旋分量的问题。
3.2 洛伦兹规范与势函数的控制方程
为保证势函数方程的唯一性，需引入洛伦兹规范：
将势函数分解式代入时域麦克斯韦方程组，结合洛伦兹规范，可推导出均匀介质中势函数满足的波动方程：
这两个方程分别为标势波动方程与矢量势波动方程，均以 ** 拉普拉斯算子▽²** 为核心。
3.3 低频极限：退化为泊松方程，消除零空间
低频极限下ω->0，时域偏导数趋向于0，势函数波动方程退化为泊松方程
泊松方程的核心优势在于：拉普拉斯算子▽²无零空间、负定且可逆。
从数学特性看，拉普拉斯算子不存在非零函数满足▽²f=0（齐次方程仅零解），对应的离散矩阵非奇异、行列式非零、可逆，可通过常规数值方法稳定求解，从根本上消除了低频崩溃。
3.4 非均匀介质中的推广与数值实现
对于非均匀介质，A-Φ 势法可进一步推广，广义洛伦兹规范下的势函数方程为：
其中，
（α为任意非零常数），适配介质参数的空间变化。
在数值实现层面，离散外微积分（DEC）是适配 A-Φ 势法的高效工具，通过构建原网格与对偶网格，分别承载矢量势、标势的离散自由度，天然规避低频崩溃问题。其核心挑战在于霍奇星算子的对角化—— 当前离散形式下该算子非对角，求逆成本较高，而 FDTD 方法中天然对角化的特性，为后续优化提供了方向。
[1] W. C. Chew, B. Zhang, and J. Zhu, “Some Selected Unsolved Problems in Classical and Quantum Electromagnetics,” PIER, vol. 180, pp. 79–87, 2024, doi: 10.2528/PIER24072910.
