电磁场问题的弱形式表述
在电磁场数值分析中，直接求解偏微分方程（PDEs）的“强形式”（Strong Form）往往面临巨大的挑战，尤其是在处理复杂几何结构、不连续介质参数或非光滑边界条件时。为了能够利用有限元法（FEM）等数值技术进行计算，我们必须将物理方程转化为“弱形式”（Weak Form）。
弱形式不仅降低了对解的连续性（光滑度）的要求，还为处理边界条件和复杂边界提供了自然的数学框架。本文将深入探讨电磁场问题从强形式到弱形式的转化逻辑、数学原理及具体推导过程。
在电磁学中，麦克斯韦方程组通常以微分方程的形式给出，这就是强形式。强形式要求解函数在定义域内必须具有足够阶数的连续导数。例如，对于二阶偏微分方程，解函数必须是二阶可微的。然而，在实际工程问题中，介质的电磁参数（如相对介电常数 $\epsilon_r$ 或磁导率 $\mu_r$）往往在界面处发生突变，这导致电磁场在界面处导数不连续，强形式在这些点上不再成立。
弱形式（又称变分形式或积分形式）通过引入“测试函数”（Test Function）并利用积分运算，将微分阶数降低。其核心思想是：不再要求方程在每一个点上都严格成立，而是要求方程在积分意义下对所有的测试函数都成立。这种方法允许解函数处于更宽泛的函数空间中（如索伯列夫空间 Sobolev Space），从而能够处理物理上更真实的不连续情况。
弱形式的通用推导步骤
将一个电磁场强形式方程转化为弱形式，通常遵循以下四个标准步骤：
- 加权（Weighting）：将强形式方程乘以一个任意的、属于特定函数空间的测试函数 $v$（或 $\mathbf{E}_t$）。
- 积分（Integration）：在整个计算域 $\Omega$ 上对加权后的方程进行积分。
- 分部积分（Integration by Parts）：利用格林公式（Green's Identities）或散度定理，将高阶导数项转化为低阶导数项与边界积分项。这一步是降低解函数光滑度要求的关键。
- 边界条件处理（Boundary Conditions）：将边界积分项与物理边界条件（Dirichlet 或 Neumann 条件）相结合。
实例推导：时谐电磁场方程
为了更直观地理解这一过程，我们以时谐（Time-Harmonic）电磁场中的矢量波动方程为例。假设在无源区域，电场 $\mathbf{E}$ 满足以下强形式方程：
$$\nabla \times (\mu^{-1} \nabla \times \mathbf{E}) - \omega^2 \epsilon \mathbf{E} = 0$$
其中，$\mu$ 是磁导率，$\epsilon$ 是介电常数，$\omega$ 是角频率。
第一步：加权与积分
引入一个测试函数 $\mathbf{E}_t$，该函数与 $\mathbf{E}$ 具有相同的数学性质（例如属于相同的向量空间）。我们将方程两端乘以 $\mathbf{E}_t$ 并对区域 $\Omega$ 进行体积分：
$$\int_{\Omega} \left[ \nabla \times (\mu^{-1} \nabla \times \mathbf{E}) \right] \cdot \mathbf{E}t , d\Omega - \int{\Omega} \omega^2 \epsilon \mathbf{E} \cdot \mathbf{E}_t , d\Omega = 0$$
第二步：利用向量恒等式进行分部积分
为了降低 $\nabla \times$ 的阶数，我们利用向量分析中的分部积分公式：
$$\int_{\Omega} (\nabla \times \mathbf{A}) \cdot \mathbf{B} , d\Omega = \int_{\Omega} \mathbf{A} \cdot (\nabla \times \mathbf{B}) , d\Omega + \int_{\partial \Omega} (\mathbf{A} \times \mathbf{n}) \cdot \mathbf{B} , d\Gamma$$
其中 $\mathbf{n}$ 是边界 $\partial \Omega$ 的单位外法向量。
令 $\mathbf{A} = \mu^{-1} \nabla \times \mathbf{E}$，$\mathbf{B} = \mathbf{E}_t$，代入上式：
$$\int_{\Omega} (\mu^{-1} \nabla \times \mathbf{E}) \cdot (\nabla \times \mathbf{E}t) , d\Omega + \int{\partial \Omega} (\mu^{-1} \nabla \times \mathbf{E} \times \mathbf{n}) \cdot \mathbf{E}t , d\Gamma - \int{\Omega} \omega^2 \epsilon \mathbf{E} \cdot \mathbf{E}_t , d\Omega = 0$$
注意到 $\mu^{-1} \nabla \times \mathbf{E}$ 与磁场 $\mathbf{H}$ 有关（根据 $\nabla \times \mathbf{E} = -j\omega\mu\mathbf{H}$），边界项实际上涉及到了磁场的切向分量。
第三步：整理得到弱形式
整理上述各项，得到最终的弱形式表述：
寻找 $\mathbf{E} \in \mathcal{V}$，使得对于所有的 $\mathbf{E}_t \in \mathcal{V}0$，均满足：
$$\int{\Omega} (\mu^{-1} \nabla \times \mathbf{E}) \cdot (\nabla \times \mathbf{E}t) , d\Omega - \int{\Omega} \omega^2 \epsilon \mathbf{E} \cdot \mathbf{E}t , d\Omega = - \int{\partial \Omega} (\mathbf{n} \times \mathbf{H}) \cdot \mathbf{E}_t , d\Gamma$$
这里，左侧构成了双线性形式 $a(\mathbf{E}, \mathbf{E}_t)$，右侧构成了线性泛函 $L(\mathbf{E}_t)$。
边界条件的引入
在弱形式中，边界条件被分为两类：
- 本质边界条件（Essential / Dirichlet Boundary Conditions）：
直接作用于解函数本身的条件，例如在完美导体（PEC）表面 $\mathbf{E} \times \mathbf{n} = 0$。在弱形式推导中，我们通过限制测试函数空间 $\mathcal{V}_0$ 的取值，使得在边界上测试函数必须满足该条件。 - 自然边界条件（Natural / Neumann Boundary Conditions）：
通过分部积分产生的边界积分项自动引入的条件。例如，在磁场边界上给定的磁场切向分量 $\mathbf{n} \times \mathbf{H} = \mathbf{M}_s$。在弱形式中，这个 $\mathbf{M}_s$ 会直接进入方程右侧的线性泛函项中。
总结
弱形式表述是连接物理连续介质力学与数值计算（如有限元法）的桥梁。通过将微分方程转化为积分方程，我们获得了以下优势：
- 降低正则性要求：允许解在介质界面处存在导数不连续。
- 几何灵活性：积分形式对复杂、不规则的边界具有极强的适应性。
- 统一框架：将边界条件自然地融入方程体系，使得数值实现更加系统化。
掌握弱形式的推导，是进行电磁场高级数值仿真研究的基石。