泊松方程（静电场）
理论与物理
泊松方程
老师，静电场的泊松方程是什么？
高斯定律 + 电场与电位关系的组合结果：
一般电介质中：
静电场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 \cdot (\varepsilon \nabla \phi) = -\rho_v$ — 静电场FEM的控制方程
- 椭圆型PDE — 与结构FEM相同的数学结构
- Dirichlet = 电极，Neumann = 对称/绝缘
Coffee Break 闲谈
泊松方程在锂离子电池设计中活跃的原因
锂离子电池充放电模拟中，需要同时求解电解质中的锂离子浓度分布和电位分布。支配该电位分布的正是泊松方程 $\nabla^2 \phi = -\rho/\varepsilon$。电池内部正极、负极、电解质的介电常数差异很大，边界条件的设定左右着精度。实际上，松下等电池制造商在进行“电极多孔结构优化”的CAE应用时，会基于泊松方程求解得到的电位分布来估算局部电流密度，从而识别易发生劣化的位置。这个不起眼的二阶微分方程，在背后支撑着最新的电池技术。
各项的物理意义
- 电场项 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：法拉第电磁感应定律。随时间变化的磁通密度产生电动势。【日常例子】自行车发电机（发电机）通过旋转磁铁使附近线圈产生电压——磁场随时间变化会感应出电场，这是该定律的直接应用。IH电磁炉也基于相同原理，高频磁场的变化在锅底感应出涡流，通过焦耳热加热。
- 磁场项 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：安培-麦克斯韦定律。电流和位移电流产生磁场。【日常例子】电线通电时周围产生磁场——这就是安培定律。电磁铁根据此原理工作，通过线圈通电产生强磁场。智能手机的扬声器也应用了电流→磁场→振膜的力这一原理。高频（GHz频段天线等）情况下，位移电流 $\partial D/\partial t$ 不可忽略，用于描述电磁波辐射。
- 高斯定律 $\nabla \cdot \mathbf{D} = \rho_v$：表明电荷是电通量的发散源。【日常例子】用垫板摩擦头发会产生静电使头发竖起——带电的垫板（电荷）放射状地发出电力线，对轻的头发施加力。电容器设计中，电极间的电场分布根据此定律计算。ESD（静电放电）对策也以基于高斯定律的电场分析为基础。
- 磁通守恒 $\nabla \cdot \mathbf{B} = 0$：表示不存在磁单极子。【日常例子】将条形磁铁切成两半也无法得到只有N极或只有S极的磁铁——N极和S极总是成对存在。这意味着磁力线描绘的是“没有起点和终点的闭合回路”。在数值分析中，为了满足此条件，采用矢量势 $\mathbf{B} = \nabla \times \mathbf{A}$ 这种公式化方法，自动保证磁通守恒。
假设条件与适用范围
- 线性材料假设：磁导率、介电常数不依赖于磁场、电场强度（饱和区域需要非线性B-H曲线）
- 准静态近似（低频）：位移电流项可忽略（$\omega \varepsilon \ll \sigma$）。涡流分析中常用
- 2D假设（截面分析）：电流方向均匀且可忽略端部效应时有效
- 各向同性假设：各向异性材料（如硅钢板的轧制方向等）需要定义方向特性
- 不适用的案例：等离子体（电离气体）、超导体、非线性光学材料需要额外的本构关系
| 类型 | 数学形式 | 物理意义 |
|---|---|---|
| Dirichlet | $\phi = \phi_0$ | 电极（电位固定） |
| Neumann | $\partial\phi/\partial n = 0$ | 对称面、绝缘面 |
| 混合 | $\varepsilon \partial\phi/\partial n = \sigma_s$ | 表面电荷密度 |
- $\nabla \cdot (\varepsilon \nabla \phi) = -\rho_v$ — 静电场FEM的控制方程
- 椭圆型PDE — 与结构FEM相同的数学结构
- Dirichlet = 电极，Neumann = 对称/绝缘
泊松方程在锂离子电池设计中活跃的原因
锂离子电池充放电模拟中，需要同时求解电解质中的锂离子浓度分布和电位分布。支配该电位分布的正是泊松方程 $\nabla^2 \phi = -\rho/\varepsilon$。电池内部正极、负极、电解质的介电常数差异很大，边界条件的设定左右着精度。实际上，松下等电池制造商在进行“电极多孔结构优化”的CAE应用时，会基于泊松方程求解得到的电位分布来估算局部电流密度，从而识别易发生劣化的位置。这个不起眼的二阶微分方程，在背后支撑着最新的电池技术。
各项的物理意义
- 电场项 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：法拉第电磁感应定律。随时间变化的磁通密度产生电动势。【日常例子】自行车发电机（发电机）通过旋转磁铁使附近线圈产生电压——磁场随时间变化会感应出电场，这是该定律的直接应用。IH电磁炉也基于相同原理，高频磁场的变化在锅底感应出涡流，通过焦耳热加热。
- 磁场项 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：安培-麦克斯韦定律。电流和位移电流产生磁场。【日常例子】电线通电时周围产生磁场——这就是安培定律。电磁铁根据此原理工作，通过线圈通电产生强磁场。智能手机的扬声器也应用了电流→磁场→振膜的力这一原理。高频（GHz频段天线等）情况下，位移电流 $\partial D/\partial t$ 不可忽略，用于描述电磁波辐射。
- 高斯定律 $\nabla \cdot \mathbf{D} = \rho_v$：表明电荷是电通量的发散源。【日常例子】用垫板摩擦头发会产生静电使头发竖起——带电的垫板（电荷）放射状地发出电力线，对轻的头发施加力。电容器设计中，电极间的电场分布根据此定律计算。ESD（静电放电）对策也以基于高斯定律的电场分析为基础。
- 磁通守恒 $\nabla \cdot \mathbf{B} = 0$：表示不存在磁单极子。【日常例子】将条形磁铁切成两半也无法得到只有N极或只有S极的磁铁——N极和S极总是成对存在。这意味着磁力线描绘的是“没有起点和终点的闭合回路”。在数值分析中，为了满足此条件，采用矢量势 $\mathbf{B} = \nabla \times \mathbf{A}$ 这种公式化方法，自动保证磁通守恒。
假设条件与适用范围
- 线性材料假设：磁导率、介电常数不依赖于磁场、电场强度（饱和区域需要非线性B-H曲线）
- 准静态近似（低频）：位移电流项可忽略（$\omega \varepsilon \ll \sigma$）。涡流分析中常用
- 2D假设（截面分析）：电流方向均匀且可忽略端部效应时有效
- 各向同性假设：各向异性材料（如硅钢板的轧制方向等）需要定义方向特性
- 不适用的案例：等离子体（电离气体）、超导体、非线性光学材料需要额外的本构关系
数值解法与实现
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]$。
总结
- 刚度矩阵$[K] = \int \varepsilon [B]^T[B] d\Omega$ — 与结构分析同形
- 一次三角形单元可手算 — 最适合FEM教学
用FDM求解泊松方程——五点星差分的直觉
用有限差分法（FDM）离散化泊松方程，在2D中会得到“关注网格点的电位 = 周围4点电位的平均值 + 电荷密度项”的形式。这就是“五点星”差分式，直观上意味着“如果没有电荷，该点的电位就是相邻点的平均值”。这个性质直接关系到SOR（逐次超松弛法）这种迭代解法的速度，超松弛因子 $\omega$ 的选择不同，收敛速度可能快一倍，也可能发散。即使在FEM主流的今天，仍有许多工程师为了教学或概念验证而用FDM实现泊松方程，它是利用矩阵“稀疏性”的数值技术的绝佳入门题材。
边单元（Nedelec元）
专用于电磁场分析的元素。自动保证切向分量的连续性，排除伪模式。是3D高频分析的标准。
节点单元
用于标量势公式化。在静磁场的标量势法或静电场分析中有效。
FEM vs BEM（边界元法）
FEM: 对应非线性材料、非均匀介质。BEM: 可自然处理无限区域（开区域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
用牛顿-拉弗森法处理B-H曲线的非线性。残差标准: $||R||/||R_0|| < 10^{-4}$ 是通用的。
频域分析
通过时间谐波假设归结为稳态问题。需要复数运算，但宽带特性需通过时域分析获取。
时域的时间步长
需要最高频率成分的1/20以下的时间步长。隐式时间积分中也可使用更大的步长，但需注意精度。
频域与时域的使用区分
频域分析类似于“将收音机调到特定频率”——可以高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——可以再现包含所有频率成分的瞬态现象，但计算成本高。
实践指南
实务
高压套管、电缆接头、半导体的掺杂分布是典型的泊松方程问题。
检查清单
- [ ] 空间电荷$\rho_v$的分布是否正确（半导体中是施主/受主浓度）
- [ ] 是否考虑了介电常数的各向异性、非线性
- [ ] 电极与接地的边界条件是否正确
- [ ] 是否确认了网格收敛性（电场是电位的微分，精度会降低一阶）
EV充电系统的绝缘设计与泊松方程
快速充电器（CHAdeMO、CCS）中使用的高压连接器的绝缘设计，离不开基于泊松方程求解的电位分布分析。充电时流过数百安培的电流，连接器内部的电位梯度在某些位置会变得陡峭，那里就成为绝缘材料的薄弱点。设计者根据泊松方程的解识别“电场集中的拐角和台阶”，并通过倒角加工或增加绝缘材料厚度来应对。实际的设计流程中，常常是先通过泊松方程得出电位分布，再将结果用作热分析的输入，进行耦合分析。一个方程在多个设计阶段都发挥着作用呢。
分析流程的比喻
电机的电磁场分析感觉上接近“给吉他调音”。调整琴弦粗细（线圈匝数）和琴桥位置（磁铁配置），引出最美的音色（高效的扭矩特性）。改变一个参数，整体平衡就会改变——所以参数化研究很重要。
初学者容易陷入的陷阱
“空气区域？为什么要用网格划分空气？”——这是几乎所有初次接触电磁场分析的人都会产生的疑问。答案是“因为磁力线也会扩散到铁芯之外”。如果分析区域紧贴铁芯，无处可去的磁通会“撞上”壁面反射，产生实际中不可能出现的磁通集中。想象一下房间太小，球在墙上不断弹跳的状态。
边界条件的思考方式
远方的边界条件虽然不起眼但超级重要。需要在数值上表达“从这里开始是无限延伸的空间”。如果设置错误，磁通就会像撞上“看不见的墙”一样被弹回。
软件比较
工具
所有静电场求解器都求解泊松方程。无需特别的工具选型。COMSOL、Maxwell、FEMM 哪一个都可以。
泊松方程求解器“多重网格法”速度快的真正原因
商用工具的泊松方程求解器很多都采用了“多重网格法”。为什么快呢？因为通常的迭代解法对于与网格细度相应的“短波误差”能快速消除，但“长波误差”消除得非常慢。多重网格法通过交替使用粗网格和细网格，可以高效地消除任何波长的误差。实现了计算量相对于问题规模N达到 $O(N)$ 的惊人性能。COMSOL或ANSYS中“大规模静电场分析意外地快速完成”，多亏了求解器的这种巧妙设计。
选型时最重要的三个问题
- “要解什么”：泊松方程（静电场）所需的物理模型、单元类型是否支持。例如，流体中是否有LES支持，结构中接触、大变形的对应能力会成为差异点。
- “谁来使用”：新手团队适合GUI充实的工具，有经验者适合脚本驱动的灵活工具。类似于汽车的AT车（GUI）和MT车（脚本）的区别。
- “要扩展到什么程度”：基于对未来分析规模扩大（HPC支持）、向其他部门扩展、与其他工具联动的前瞻性选择，有助于长期降低成本。
尖端技术
尖端
- 泊松-玻尔兹曼方程 — 电解质溶液中的离子分布。生物分子模拟
- 半导体器件模拟 — 泊松 + 漂移扩散方程的耦合（Sentaurus, Silvaco）
- 等离子体模拟 — 泊松 + 粒子法（PIC法）
EV充电系统的绝缘设计与泊松方程
快速充电器（CHAdeMO、CCS）中使用的高压连接器的绝缘设计，离不开基于泊松方程求解的电位分布分析。充电时流过数百安培的电流，连接器内部的电位梯度在某些位置会变得陡峭，那里就成为绝缘材料的薄弱点。设计者根据泊松方程的解识别“电场集中的拐角和台阶”，并通过倒角加工或增加绝缘材料厚度来应对。实际的设计流程中，常常是先通过泊松方程得出电位分布，再将结果用作热分析的输入，进行耦合分析。一个方程在多个设计阶段都发挥着作用呢。
分析流程的比喻
电机的电磁场分析感觉上接近“给吉他调音”。调整琴弦粗细（线圈匝数）和琴桥位置（磁铁配置），引出最美的音色（高效的扭矩特性）。改变一个参数，整体平衡就会改变——所以参数化研究很重要。
初学者容易陷入的陷阱
“空气区域？为什么要用网格划分空气？”——这是几乎所有初次接触电磁场分析的人都会产生的疑问。答案是“因为磁力线也会扩散到铁芯之外”。如果分析区域紧贴铁芯，无处可去的磁通会“撞上”壁面反射，产生实际中不可能出现的磁通集中。想象一下房间太小，球在墙上不断弹跳的状态。
边界条件的思考方式
远方的边界条件虽然不起眼但超级重要。需要在数值上表达“从这里开始是无限延伸的空间”。如果设置错误，磁通就会像撞上“看不见的墙”一样被弹回。
なった
詳しく
報告