电荷守恒是电磁学最基本的原理之一，它在数学上可以用 连续性方程（continuity equation）来表述。连续性方程把电荷的时空分布与电流的流动联系起来，体现了电荷既不会凭空产生，也不会凭空消失的物理事实。下面从数学推导、物理意义以及实际应用三个层面系统阐述电荷守恒方程的连续性形式。

1.1 基本假设

• 体电荷密度 (\rho(\mathbf{r},t))：单位体积内的电荷量，单位 (\mathrm{C/m^3})。
• 电流密度 (\mathbf{J}(\mathbf{r},t))：单位面积、单位时间通过的电荷量，单位 (\mathrm{A/m^2})。

1.2 控制体与守恒原理

取任意闭合体积 (V)（控制体），其边界为闭合曲面 (S)。在时间 ([t, t+\Delta t]) 内，进入或离开该体积的电荷量应等于体积内电荷的增减：

[
\underbrace{\int_{V}\rho(\mathbf{r},t+\Delta t),dV}_{\text{末时刻电荷}}

\underbrace{\int_{V}\rho(\mathbf{r},t),dV}_{\text{初时刻电荷}}

-\underbrace{\oint_{S}\mathbf{J}\cdot d\mathbf{S},\Delta t}_{\text{净流出电荷}}
]

两边同除以 (\Delta t)，并令 (\Delta t\to0)，得到

[
\frac{d}{dt}\int_{V}\rho,dV = -\oint_{S}\mathbf{J}\cdot d\mathbf{S}.
]

1.3 使用散度定理

根据散度定理

[
\oint_{S}\mathbf{J}\cdot d\mathbf{S}= \int_{V}\nabla!\cdot!\mathbf{J},dV,
]

于是

[
\int_{V}\frac{\partial\rho}{\partial t},dV

• \int_{V}\nabla!\cdot!\mathbf{J},dV = 0.
  ]

因为体积 (V) 任意，积分式只能在积分子内部恒等，于是得到 连续性方程的局部形式：

[
\boxed{\displaystyle \frac{\partial\rho}{\partial t} + \nabla!\cdot!\mathbf{J}=0 } \tag{1}
]

这就是电荷守恒的微分表达式。

2. 物理图像与直观解释

• (\partial\rho/\partial t) 表示某点电荷密度随时间的变化率。如果该点电荷密度在增大，则该项为正。
• (\nabla!\cdot!\mathbf{J}) 是电流密度的散度，描述电流在该点的“汇聚”或“发散”。若电流从该点向外流出，则散度为正，表示电荷在该点被抽走。
• 方程 (1) 说明：电荷密度的增加只能由电流的汇聚导致，电荷密度的减少只能由电流的发散导致，两者恰好相抵，保证总电荷守恒。

2.1 一维示例

在一根均匀导线中，仅考虑 (x) 方向的电流 (J_x(x,t)) 与线电荷密度 (\lambda(x,t))（(\rho) 换成 (\lambda)），连续性方程化为

[
\frac{\partial \lambda}{\partial t} + \frac{\partial J_x}{\partial x}=0.
]

若在某段导线两端分别接入电源，使得左端电流为 (I_0)，右端电流为 (I_0+\Delta I)，则 (\partial J_x/\partial x = \Delta I/L)（(L) 为段长），对应的 (\partial\lambda/\partial t) 为正，说明该段电荷密度随时间上升，形成瞬时的电荷堆积。

3. 与麦克斯韦方程组的联系

连续性方程并非独立于麦克斯韦方程组，而是 由麦克斯韦方程组的两条基本方程（高斯定律和安培-麦克斯韦定律）相容性条件推出的。

• 高斯定律（电场）
  [
  \nabla!\cdot!\mathbf{E}= \frac{\rho}{\varepsilon_0}
  ]
• 安培-麦克斯韦定律（磁场）
  [
  \nabla!\times!\mathbf{B}= \mu_0\mathbf{J}+ \mu_0\varepsilon_0\frac{\partial\mathbf{E}}{\partial t}
  ]

对安培-麦克斯韦定律取散度，利用向量恒等式 (\nabla!\cdot!(\nabla!\times!\mathbf{B})=0)，得到

[
0 = \mu_0\nabla!\cdot!\mathbf{J}+ \mu_0\varepsilon_0\frac{\partial}{\partial t}\bigl(\nabla!\cdot!\mathbf{E}\bigr)
= \mu_0\nabla!\cdot!\mathbf{J}+ \mu_0\varepsilon_0\frac{\partial}{\partial t}\left(\frac{\rho}{\varepsilon_0}\right).
]

整理即得

[
\frac{\partial\rho}{\partial t}+ \nabla!\cdot!\mathbf{J}=0,
]

这正是连续性方程 (1)。因此 电荷守恒是麦克斯韦方程组内部的一致性要求，不需要额外的假设。

4. 实际应用案例

4.1 电容充放电过程

在平行板电容器中，板间的电场随时间变化，导致位移电流 (\mathbf{J}_D=\varepsilon_0\partial\mathbf{E}/\partial t)。若仅考虑导体内部的导电电流 (\mathbf{J})，连续性方程写成

[
\frac{\partial\rho}{\partial t}+ \nabla!\cdot!(\mathbf{J}+\mathbf{J}_D)=0.
]

这说明 位移电流同样参与电荷守恒，是解释“电容器内部没有真实电荷流动却仍然满足电流连续性”的关键。

4.2 天线辐射中的电荷收敛

在振荡天线的驱动点，电流随时间呈正弦变化。连续性方程告诉我们，在电流节点（电流为零）处必然出现 电荷积聚（(\partial\rho/\partial t\neq0)），这正是天线辐射的源头——随时间变化的电荷分布产生电磁波。

5. 常见误区与注意事项

• 误区 1：连续性方程只适用于稳态
  实际上，方程 (1) 对任意时变情形均成立，稳态（(\partial\rho/\partial t=0)）只是其特例，此时 (\nabla!\cdot!\mathbf{J}=0)，即电流是无源的闭合回路。

• 误区 2：电荷守恒可以被“破坏”
  在经典电磁学框架下，电荷守恒是不可违背的。如果在数值仿真或实验中出现不满足 (1) 的现象，往往是 网格离散、边界条件或测量误差 所致，需要检查模型的完整性。

• 注意事项

  • 在使用 有限体积法 或 有限元法 进行电磁数值求解时，必须在离散方程中显式保留连续性约束，以避免出现“伪电荷”。
  • 当涉及 等离子体 或 半导体 中的产生/复合过程时，连续性方程会加入源项 (G - R)（产生率减复合率），形式为
    [
    \frac{\partial\rho}{\partial t}+ \nabla!\cdot!\mathbf{J}= G - R,
    ]
    这仍然是守恒原理的扩展形式。

6. 小结

• 连续性方程 (\displaystyle \frac{\partial\rho}{\partial t}+ \nabla!\cdot!\mathbf{J}=0) 是电荷守恒的局部数学表达，直接来源于控制体的守恒原理并通过散度定理得到。
• 它在物理上描述了电荷密度的时空变化与电流的散度之间的精确平衡，是 麦克斯韦方程组内部的一致性条件。
• 在电容充放电、天线辐射、等离子体动力学等实际问题中，连续性方程提供了分析电荷分布演化、验证数值模型以及设计电磁装置的基本工具。

掌握连续性方程的推导与应用，不仅有助于深入理解电磁场的本质，也为后续学习更高级的电磁波、量子电动力学奠定坚实的理论基础。