泊松方程静电场
泊松方程静电场的理论基础
泊松方程
老师，静电场的泊松方程是什么？
高斯定律＋电场与电位关系的结合结果:
对于一般介质:
静电场FEM的支配方程就是这个。给定电荷分布$\rho_v$时，求解电位$\phi$。
这与结构的平衡方程$\nabla \cdot \sigma + f = 0$很相似呀。
数学上是同一类椭圆型偏微分方程。结构的杨氏模量→介电常数$\varepsilon$、外力→电荷密度$\rho_v$、位移→电位$\phi$，简单替换即可。
边界条件
| 种类 | 数学形式 | 物理含义 |
|---|---|---|
| Dirichlet | $\phi = \phi_0$ | 电极（电位固定） |
| Neumann | $\partial\phi/\partial n = 0$ | 对称面、绝缘面 |
| 混合 | $\varepsilon \partial\phi/\partial n = \sigma_s$ | 表面电荷密度 |
总结
锂离子电池设计中泊松方程大显身手的原因
在锂离子电池的充放电仿真中，需要同时求解电解质中的锂离子浓度分布和电位分布。控制电位分布的正是泊松方程 $\nabla^2 \phi = -\rho/\varepsilon$。电池内部的正极、负极和电解质的介电常数差异很大，边界条件的设定会直接影响精度。实际上，松下等电池制造商在使用CAE进行"电极多孔结构最优化"时，就是通过求解泊松方程获得的电位分布来推估局部电流密度，进而识别容易劣化的位置。这个不起眼的二阶微分方程正在背后支撑着最新的电池技术。
泊松方程静电场的数值计算手法
FEM离散化
用Galerkin方法进行弱形式化，按单元进行离散化:
$[B]$是形状函数的梯度矩阵。这是结构解析中$[B]^T [D] [B]$的$[D]$替换为$\varepsilon [I]$后的形式。
二维三角形单元的话$B$是常数（一阶单元），可以手算呀。
正是。CST三角形单元的静电场版是教科书的练习题最好的素材。单元刚度矩阵是$\varepsilon \cdot A_e / (4A_e) \cdot [b_i b_j + c_i c_j]$。
总结
用FDM求解泊松方程——5点星差分的直观理解
用有限差分法（FDM）对泊松方程进行离散化，二维情况下得到"注目格点的电位＝周围4点电位的平均＋电荷密度项"的形式。这就是"5点星"差分式，直观上意思是"无电荷处，该点电位是邻近点的平均值"。这个性质与SOR（逐次超松弛法）等迭代解法的速度密切相关，选择不同的超松弛系数 $\omega$ 可以让收敛速度快一倍，也可能导致发散。虽然现在FEM占主流，但很多工程师仍用FDM实现泊松方程来进行教学或概念验证，它是理解矩阵"稀疏性"和数值技巧的最好入门素材。
泊松方程静电场的实务应用
实务
高压套管、电缆接头、半导体掺杂分布是典型的泊松方程问题。
检查清单
电动汽车充电系统绝缘设计与泊松方程
快速充电器（CHAdeMO、CCS）所用的高压连接器绝缘设计必须进行泊松方程电位分布解析。充电时数百安培的电流流过，连接器内部某些地方的电位梯度会变得很陡，这些地方就是绝缘材料的最薄弱点。设计人员通过泊松方程的解获得"电场集中的角部和阶梯处"，再通过倒角加工或增加绝缘厚度来应对。实际的设计流程中，先用泊松方程得到电位分布，再将其作为输入用于热解析，经常形成联合分析。一个方程在多个设计阶段中发挥作用呢。
泊松方程静电场的软件比较
工具
所有静电场求解器都求解泊松方程。不需要特别选择工具。COMSOL、Maxwell、FEMM任选其一。
泊松方程求解器的"多重网格法"为什么这么快
商用工具的泊松方程求解器中很多采用了"多重网格法"。速度快的原因在于，通常的迭代解法对网格细度相应的"短波长误差"能迅速消除，但"长波长误差"的消除很慢。多重网格法通过交替使用粗网格和细网格，无论误差波长如何都能高效消除。计算量相对问题规模N的复杂度达到 $O(N)$ 这种惊人的性能。这就是在COMSOL或ANSYS中"大规模静电解析意外地快速完成"的原因——都是求解器设计的功劳。
泊松方程静电场的先端研究
先端
用神经网络求解泊松方程——Physics-Informed NN的冲击
从2019年左右，"Physics-Informed Neural Network（PINN）"这种方法备受关注。这种方法是在神经网络的损失函数中加入支配方程（包括泊松方程）的残差，通过学习来求解边值问题。不用生成网格、可在任意点预测，是主要优势。但收敛慢、精度保证困难等问题使得"完全替代FEM还言之过早"是现状。不过，将来在电池和半导体器件仿真中进行"大量设计参数的高速扫描"时，PINN可能会成为主流。那一天的到来也许不远。
泊松方程静电场的故障排除
故障
"求解器不收敛"——泊松方程中常见的陷阱
泊松方程的FEM解析中"迭代解法不收敛"是常见的故障。最典型的原因是"边界条件设置错误"。如果没有设置迪利克雷边界（电位固定），联立方程的系数矩阵就会变成奇异矩阵，解无法确定。当电位的绝对值无所谓、只有电场有意义时，做法是"设一个参考电位点来接地"。还有一个常见的陷阱是介电常数的单位错误——本该输入相对介电常数 $\varepsilon_r$ 的地方直接填入真空介电常数 $\varepsilon_0 = 8.85\times10^{-12}$ F/m，结果电位数值偏离12个数量级。这是老牌工程师都踩过的坑。
价值
详细
错误