均匀无源介质中的波动方程推导

在电磁场与电磁波理论的研究中,麦克斯韦方程组构成了经典电动力学的基石。然而,这些方程描述的是电场与磁场的一阶空间与时间导数之间的耦合关系,直接求解并不直观。为了揭示电磁场在空间中传播的本质,我们需要将这些一阶微分方程解耦,转化为描述电场或磁场单独行为的二阶偏微分方程——即波动方程。

本文将聚焦于总览级原理,系统推导在理想、均匀、无源的线性介质中电磁场的波动方程,并简要探讨其物理意义与全景应用。
在推导波动方程之前,必须明确所处物理环境的本构关系与基本假设。所谓“均匀无源介质”,通常具备以下数学和物理特征:

  • 无源(Source-free):介质内部不存在自由电荷与自由电流,即自由电荷密度 $\rho_f = 0$,传导电流密度 $\mathbf{J}_f = 0$。
  • 均匀(Homogeneous):介质的宏观电磁参数(介电常数 $\varepsilon$、磁导率 $\mu$、电导率 $\sigma$)在空间各点均保持常数,不随空间坐标的变化而改变。
  • 线性与各向同性:介质对电磁场的响应是线性的,且在各个方向上性质相同,因此 $\varepsilon$ 和 $\mu$ 均为标量常数。

在上述条件下,宏观介质中的麦克斯韦方程组(微分形式)可以简化为:

  1. 库伦定律(电场高斯定律):$\nabla \cdot \mathbf{D} = 0 \implies \nabla \cdot \mathbf{E} = 0$
  2. 磁场高斯定律:$\nabla \cdot \mathbf{B} = 0 \implies \nabla \cdot \mathbf{B} = 0$
  3. 法拉第电磁感应定律:$\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}$
  4. 安培-麦克斯韦定律:$\nabla \times \mathbf{H} = \mathbf{J}_f + \frac{\partial \mathbf{D}}{\partial t} \implies \nabla \times \mathbf{H} = \frac{\partial \mathbf{D}}{\partial t}$

同时,介质的本构关系为 $\mathbf{D} = \varepsilon \mathbf{E}$ 且 $\mathbf{B} = \mu \mathbf{H}$。

电场波动方程的推导

我们的目标是从法拉第定律和安培-麦克斯韦定律出发,消去磁场变量 $\mathbf{H}$(或 $\mathbf{B}$),仅留下电场 $\mathbf{E}$。

首先,对法拉第电磁感应定律两边取旋度(Curl):
$$\nabla \times (\nabla \times \mathbf{E}) = \nabla \times \left( -\frac{\partial \mathbf{B}}{\partial t} \right)$$

利用算子的线性性质,可以将时间偏导数与空间旋度算子交换位置:
$$\nabla \times (\nabla \times \mathbf{E}) = -\frac{\partial}{\partial t} (\nabla \times \mathbf{B})$$

根据本构关系,在均匀介质中 $\mathbf{B} = \mu \mathbf{H}$,其中 $\mu$ 为常数,故可提出:
$$\nabla \times (\nabla \times \mathbf{E}) = -\mu \frac{\partial}{\partial t} (\nabla \times \mathbf{H})$$

接下来,代入安培-麦克斯韦定律中的 $\nabla \times \mathbf{H} = \varepsilon \frac{\partial \mathbf{E}}{\partial t}$:
$$\nabla \times (\nabla \times \mathbf{E}) = -\mu \varepsilon \frac{\partial^2 \mathbf{E}}{\partial t^2}$$

在向量微积分中,有一个著名的矢量恒等式:
$$\nabla \times (\nabla \times \mathbf{E}) = \nabla(\nabla \cdot \mathbf{E}) - \nabla^2 \mathbf{E}$$

由于介质是无源的且电荷密度为零,根据电场高斯定律 $\nabla \cdot \mathbf{E} = 0$,矢量恒等式的第一项归零,简化为:
$$\nabla \times (\nabla \times \mathbf{E}) = -\nabla^2 \mathbf{E}$$

将此代入前式,即得到电场的齐次波动方程:
$$\nabla^2 \mathbf{E} - \mu \varepsilon \frac{\partial^2 \mathbf{E}}{\partial t^2} = 0$$

磁场波动方程的推导

采用完全对称的逻辑,我们可以推导出磁感应强度 $\mathbf{B}$(或磁场强度 $\mathbf{H}$)的波动方程。

从安培-麦克斯韦定律出发,两边取旋度:
$$\nabla \times (\nabla \times \mathbf{H}) = \nabla \times \left( \varepsilon \frac{\partial \mathbf{E}}{\partial t} \right)$$

交换空间与时间导数,并代入本构关系 $\mathbf{D} = \varepsilon \mathbf{E}$:
$$\nabla \times (\nabla \times \mathbf{H}) = \varepsilon \frac{\partial}{\partial t} (\nabla \times \mathbf{E})$$

再代入法拉第电磁感应定律 $\nabla \times \mathbf{E} = -\mu \frac{\partial \mathbf{H}}{\partial t}$:
$$\nabla \times (\nabla \times \mathbf{H}) = -\varepsilon \mu \frac{\partial^2 \mathbf{H}}{\partial t^2}$$

利用相同的矢量恒等式 $\nabla \times (\nabla \times \mathbf{H}) = \nabla(\nabla \cdot \mathbf{H}) - \nabla^2 \mathbf{H}$,并结合磁场高斯定律 $\nabla \cdot \mathbf{H} = 0$,左侧化简为 $-\nabla^2 \mathbf{H}$。最终得到磁场的齐次波动方程:
$$\nabla^2 \mathbf{H} - \mu \varepsilon \frac{\partial^2 \mathbf{H}}{\partial t^2} = 0$$

波动方程的物理意义与全景展望

通过上述推导,我们成功将麦克斯韦的一阶耦合方程组转化为解耦的二阶偏微分方程。这些方程揭示了电磁场演化的核心规律:

  • 波速的本质:方程中的系数项 $\mu \varepsilon$ 具有时间的平方除以空间平方的量纲。定义相速度 $v = \frac{1}{\sqrt{\mu \varepsilon}}$,在真空(自由空间)中则退化为光速 $c = \frac{1}{\sqrt{\mu_0 \varepsilon_0}}$。这在物理史上首次从理论上证明了光本身就是一种电磁波。
  • 统一框架的体现:电场波动方程与磁场波动方程具有完全相同的数学形式,表明电场与磁场在无源空间中对称地以波的形式向外传播,且两者相互激发。

在更广泛的电磁理论全景中,从均匀无源介质的理想波动方程出发,后续可以进一步推导出有源介质中的非齐次波动方程、边界条件下的波导传播特性,以及电磁场的能量传输(坡印廷定理)等高级课题。掌握这一基础推导,是理解复杂电磁环境与现代无线通信技术的必经之路。