含点电荷与固定电势导体的静电势泊松方程数值求解方法咨询

架构老林

2026-4-22

嗨，看起来你已经把问题的物理背景理得很清楚了——这其实是经典的点电荷+等势导体的静电场问题，数值求解完全可行，我给你梳理几个常用的思路和实操步骤：

首先先明确一下咱们要解决的核心方程与边界条件，方便后续对应数值方法：

导体外部区域的泊松方程：$$\nabla^2\phi = k\delta(\mathbf{r})$$
边界条件：
1. 导体边界$\partial\alpha$上：$\phi = \phi_0$
2. 无穷远边界（数值计算中需做截断处理）：$\phi \to 0$

一、有限差分法（FDM）——最容易上手的入门方案

适合规则形状的计算区域（比如长方体、球体截断域），步骤非常直观：

截断计算区域：因为无穷远没法直接计算，取一个足够大的立方体/球体区域$\Omega$，把点电荷和导体都完整包含在内，$\Omega$的外边界直接设为$\phi=0$（只要$\Omega$足够大，这个近似的误差可以忽略不计）。
网格离散化：把$\Omega$划分成均匀的网格（三维就是一个个小立方体单元），尽量让点电荷所在的原点正好落在某个网格点上，能大幅简化后续源项处理。
处理点电荷的δ函数源：如果原点是网格点，那么在该点的离散方程中直接加入源项$k/\Delta V$（$\Delta V$是单个网格单元的体积）；如果原点在网格内部，可以把源项平均分配给周围的网格点，但还是建议调整网格让原点落在节点上。
离散泊松方程：对每个非导体内部、非点源的网格点，用二阶中心差分近似拉普拉斯算子，比如三维中的形式：
$$\frac{\phi_{i+1,j,k} - 2\phi_{i,j,k} + \phi_{i-1,j,k}}{\Delta x^2} + \frac{\phi_{i,j+1,k} - 2\phi_{i,j,k} + \phi_{i,j-1,k}}{\Delta y^2} + \frac{\phi_{i,j,k+1} - 2\phi_{i,j,k} + \phi_{i,j,k-1}}{\Delta z^2} = 0$$
导体内部的网格点直接固定为$\phi_0$，不需要参与后续求解。
求解线性方程组：把所有网格点的方程整理成$\mathbf{A}\phi = \mathbf{b}$的线性方程组形式，用迭代法（比如Gauss-Seidel、SOR迭代）或者直接法求解，迭代法更适合大网格的场景。

因为导体是等势体，边界元法不需要处理整个空间，只需要离散导体的边界，计算量会小很多：

利用格林函数分解解：外部区域的电势可以拆成两部分：点电荷的原始电势，加上导体表面感应电荷产生的电势。设导体表面的感应电荷面密度为$\sigma(\mathbf{r}')$，那么外部电势的表达式为：
$$\Phi(\mathbf{r}) = -\frac{k}{4\pi|\mathbf{r}|} + \int_{\partial\alpha} \frac{\sigma(\mathbf{r}')}{4\pi|\mathbf{r}-\mathbf{r}'|} dS'$$
代入边界条件建立方程：在导体边界$\partial\alpha$上，$\Phi(\mathbf{r}) = \phi_0$，所以对边界上的每个点$\mathbf{r} \in \partial\alpha$，都有：
$$\phi_0 = -\frac{k}{4\pi|\mathbf{r}|} + \int_{\partial\alpha} \frac{\sigma(\mathbf{r}')}{4\pi|\mathbf{r}-\mathbf{r}'|} dS'$$
离散边界求解感应电荷：把导体边界$\partial\alpha$划分成若干小面元（比如三角形面元），假设每个面元上的$\sigma$是常数，把积分方程转化为线性方程组，求解出每个面元的感应电荷密度。
计算任意外部点的电势：得到$\sigma$后，代入上面的电势表达式，就能计算出外部任意位置的$\Phi(\mathbf{r})$了。

如果导体的形状很不规则，有限元法的网格适配性会更强：

备注：内容来源于stack exchange，提问作者Adrien Amour