4. 泊松方程与拉普拉斯方程

前面几节我们从库仑定律和叠加原理出发，得到了电位与电荷分布之间的积分关系。但在实际工程中，电荷分布往往不是已知量，真正已知的是导体上的电位、导体表面的电荷总量或不同介质分界面上的衔接条件。这时，把静电场规律写成微分方程形式，再配合边界条件求解，才是更通用的做法。本节要建立的正是这套“方程 + 边界条件”的求解框架。

一、从高斯定理到泊松方程

静电场的基本微分方程来自高斯定理的微分形式。在真空中，电位移矢量与电场强度的关系为 D = ε₀E，而 E = −∇φ（负号表示电场指向电位降低的方向，这是全文统一的符号约定）。代入高斯定理 ∇·D = ρ，得到：

∇·(ε₀(−∇φ)) = ρ
⇒ ∇²φ = −ρ/ε₀        （真空中泊松方程）

当讨论区域内没有自由电荷，即 ρ = 0 时，方程退化为 ∇²φ = 0，这就是拉普拉斯方程。因此可以这样理解：泊松方程描述“有源区”的电位分布，拉普拉斯方程描述“无源区”的电位分布，两者是同一物理规律在不同条件下的两种形式。

在直角坐标系中，拉普拉斯算子展开为 ∂²φ/∂x² + ∂²φ/∂y² + ∂²φ/∂z²。若场沿某一方向不变（如无限长带电圆柱），可退化为二维形式，只保留两个坐标的二阶偏导，这也是工程计算中常用的简化手段。

二、不同媒质中的形式

当空间充满均匀、线性、各向同性的电介质时，D = εE，ε 为常数，可以提到算子外面，于是泊松方程写成 ∇²φ = −ρ/ε。与真空形式相比，只是把 ε₀ 换成了介质的介电常数 ε = ε_rε₀。这说明介电常数越大，同样体电荷密度产生的电位弯曲越平缓，电场被削弱。

如果媒质不均匀，ε 随位置变化，就不能再把它提出算子，方程变为 ∇·(ε∇φ) = −ρ。此时在两种介质的分界面上，电位连续（φ₁ = φ₂），而电位移的法向分量满足 ε₁E₁ₙ − ε₂E₂ₙ = σ_s（σ_s 为分界面自由面电荷密度）；若分界面上无自由电荷，则 ε₁∂φ₁/∂n = ε₂∂φ₂/∂n。这两条衔接条件是把不同区域的解“缝合”起来的关键，缺一不可。

三、边值问题的分类

仅有微分方程不足以确定唯一解，必须补充边界条件。按边界上给定的物理量不同，静电场边值问题分为三类：

• 第一类（狄利克雷问题）：已知边界上的电位值 φ|_S = f。例如已知两平行导体板的电位分别为 0 V 和 100 V。
• 第二类（诺伊曼问题）：已知边界上的电位法向导数，等价于已知表面电荷密度 ∂φ/∂n|_S = g。例如已知导体表面电荷分布。
• 第三类（混合问题）：一部分边界给电位，另一部分给法向导数。实际工程中导体与介质混合边界最常见这种情况。

唯一性定理保证：只要满足泊松（或拉普拉斯）方程并满足给定边界条件，解就是唯一的。这条定理非常重要，它允许我们用电像法、分离变量法等“猜解 + 验证”的手段求解，而不必担心漏解。

四、一个可动手的求解示例

设两块无限大平行导体板相距 d = 10 mm，板间为空气（ε ≈ ε₀），左板电位 0 V，右板电位 100 V，板间无体电荷。取 x 轴垂直板面、左板为 x = 0，则问题为一维拉普拉斯方程：

d²φ/dx² = 0，  0 < x < d
φ(0) = 0，φ(d) = 100 V

积分两次得 φ = C₁x + C₂，代入边界条件得 C₂ = 0，C₁ = 100/d = 10⁴ V/m。于是 φ(x) = 10⁴x（V），电场 E = −dφ/dx = −10⁴ V/m，方向由高电位板指向低电位板，与负号约定一致。若在板间插入 ε_r = 4 的介质且保持两板电位不变，由于无自由体电荷，电位分布仍为线性，但介质中的电场强度降为原来的 1/4，导体表面电荷密度相应改变——这正是介质削弱电场的定量体现。

若把右板改为“已知表面电荷密度 σ_s”，则边界条件变为第二类，此时电位分布需由 σ_s 反推，解的形式仍为线性但斜率由 σ_s/ε 决定。可见边界条件类型直接决定了解的确定方式。

五、求解流程与常见易错点

工程上处理静电场边值问题，通常按以下顺序推进：先判断求解区域内是否有自由体电荷，决定用泊松还是拉普拉斯方程；再划分媒质区域，写出各区的方程与界面衔接条件；然后根据导体电位或电荷情况确定边界条件类型；最后选择分离变量、电像法或数值方法求解并验证。

掌握本节建立的“方程—衔接—边界”框架后，后续的分离变量法、镜像法和有限差分法都只是在这一框架下选择不同的求解工具而已。