毕奥-萨伐尔定律
理论与物理
毕奥-萨伐尔定律
老师，毕奥-萨伐尔定律是什么？
计算电流产生磁场的基本定律。相当于静电场的库仑定律的磁场版本。
$\mu_0 = 4\pi \times 10^{-7}$ H/m（真空磁导率）。微小电流元$Id\mathbf{l}$在距离$r$的点处产生的磁通密度$d\mathbf{B}$。
和库仑定律很像，但区别在于有向量积（叉积）。
是的。磁场方向由电流与位置向量的右手定则决定。库仑力是径向的，而磁场是切向的。
典型磁场的解析解
| 电流分布 | 磁通密度 $B$ |
|---|---|
| 无限长直线电流 $I$ | $B = \mu_0 I / (2\pi r)$ |
| 圆形线圈（中心） | $B = \mu_0 I / (2a)$，$a$: 线圈半径 |
| 螺线管（内部） | $B = \mu_0 n I$，$n$: 匝数/m |
| 亥姆霍兹线圈 | $B = (4/5)^{3/2} \mu_0 n I / a$ |
总结
- $d\mathbf{B} = \mu_0/(4\pi) \cdot I d\mathbf{l} \times \hat{r} / r^2$ — 电流→磁场的基本定律
- 右手定则 — 磁场方向是电流与位置的叉积
- 用解析解验证FEM — 螺线管、亥姆霍兹线圈
毕奥与萨伐尔——两人其实是共同研究的搭档
“毕奥-萨伐尔定律”冠以两人的名字，但据说他们共同实验的时间只有短短几周。让-巴蒂斯特·毕奥是实验物理学家，费利克斯·萨伐尔是医生兼物理学家，这对奇特的组合在1820年奥斯特发现后的几个月内就导出了电流与磁场的定量关系。另一方面，毕奥后来也因与安培围绕该成果的优先权激烈对立而闻名。科学史的背后，总是围绕着发现优先权的戏剧性故事。
各项的物理意义
- 电场项 $\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中的磁场分析
磁场也用FEM求解吗？
相当于静电场电位$\phi$的是矢量势$\mathbf{A}$呢。
但$\mathbf{A}$是矢量（3个分量），所以比标量$\phi$的自由度多。2D情况下只需$A_z$（1个分量）即可，效率高。3D则需要处理规范条件（$\nabla \cdot \mathbf{A} = 0$）。
边元（Nédélec元）
3D磁场FEM中，边元（edge element）是标准。不是将自由度分配给节点，而是分配给边，从而自然地满足$\mathbf{B}$法向分量的连续性。
总结
- $\nabla \times (\nu \nabla \times \mathbf{A}) = \mathbf{J}$ — 静磁场FEM的控制方程
- 2D中是$A_z$的标量问题 — 高效
- 3D中边元是标准 — 自然地处理规范条件
Coffee Break 闲谈
毕奥-萨伐尔数值积分的“奇点问题”——位置不同答案会爆炸
实现毕奥-萨伐尔定律数值计算时最大的陷阱是“当评估点离电流元太近时，分母趋近于零导致值发散”。实际工作中，不是将电流当作细导线处理，而是当作具有有限横截面积的实体模型处理，或者需要进行适当的偏移处理，使评估点远离电流元。特别是在细线圈的端部或折返部分进行评估时要格外注意，采用根据部位改变积分点密度的自适应积分是有效的。
边元（Nedelec元）
专用于电磁场分析的单元。自动保证切向分量的连续性，排除伪模式。3D高频分析的标准。
节点元
用于标量势公式化。在静磁场的标量势法或静电场分析中有效。
FEM vs BEM（边界元法）
FEM: 对应非线性材料·非均匀介质。BEM: 自然地处理无限域（开域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
用牛顿-拉弗森法处理B-H曲线的非线性。残差标准: $||R||/||R_0|| < 10^{-4}$ 是通用的。
频域分析
通过时间谐波假设归结为稳态问题。需要复数运算，但宽带特性需通过时域分析获取。
时域的时间步长
需要最高频率成分的1/20以下的时间步长。隐式时间积分中也可使用更大的步长，但需注意精度。
频域与时域的区分使用
频域分析类似于“调谐收音机到特定频率”——可以高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——可以再现包含所有频率成分的瞬态现象，但计算成本高。
3D磁场FEM中，边元（edge element）是标准。不是将自由度分配给节点，而是分配给边，从而自然地满足$\mathbf{B}$法向分量的连续性。
- $\nabla \times (\nu \nabla \times \mathbf{A}) = \mathbf{J}$ — 静磁场FEM的控制方程
- 2D中是$A_z$的标量问题 — 高效
- 3D中边元是标准 — 自然地处理规范条件
毕奥-萨伐尔数值积分的“奇点问题”——位置不同答案会爆炸
实现毕奥-萨伐尔定律数值计算时最大的陷阱是“当评估点离电流元太近时，分母趋近于零导致值发散”。实际工作中，不是将电流当作细导线处理，而是当作具有有限横截面积的实体模型处理，或者需要进行适当的偏移处理，使评估点远离电流元。特别是在细线圈的端部或折返部分进行评估时要格外注意，采用根据部位改变积分点密度的自适应积分是有效的。
边元（Nedelec元）
专用于电磁场分析的单元。自动保证切向分量的连续性，排除伪模式。3D高频分析的标准。
节点元
用于标量势公式化。在静磁场的标量势法或静电场分析中有效。
FEM vs BEM（边界元法）
FEM: 对应非线性材料·非均匀介质。BEM: 自然地处理无限域（开域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
用牛顿-拉弗森法处理B-H曲线的非线性。残差标准: $||R||/||R_0|| < 10^{-4}$ 是通用的。
频域分析
通过时间谐波假设归结为稳态问题。需要复数运算，但宽带特性需通过时域分析获取。
时域的时间步长
需要最高频率成分的1/20以下的时间步长。隐式时间积分中也可使用更大的步长，但需注意精度。
频域与时域的区分使用
频域分析类似于“调谐收音机到特定频率”——可以高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——可以再现包含所有频率成分的瞬态现象，但计算成本高。
实践指南
实务
线圈磁场设计、MRI磁体、电磁铁、继电器、传感器的磁场计算是主要应用。
检查清单
脑磁计（MEG）的传感器布局设计会用到毕奥-萨伐尔定律
用于测量大脑神经活动的脑磁计（MEG：脑磁图仪）是一种通过SQUID传感器检测神经细胞微弱电流（皮安量级）产生磁场的装置。其磁场计算就使用了毕奥-萨伐尔定律，数值积分对于优化传感器布局以最精确地定位信号源至关重要。静磁场分析技术支撑着尖端医疗设备的设计，所以最好摒弃“电磁学基础在现场用不上”这种先入为主的观念。
分析流程的比喻
电机的电磁场分析类似于“给吉他调音”。调整琴弦粗细（线圈匝数）和琴桥位置（磁铁布局），以引出最美妙的音色（高效的扭矩特性）。改变一个参数，整体平衡就会改变——所以参数化研究很重要。
初学者容易陷入的陷阱
“空气区域？为什么要用网格划分空气？”——这是几乎所有初次接触电磁场分析的人都会有的疑问。答案是“因为磁力线也会扩散到铁心之外”。如果将分析区域限制在紧贴铁心，无处可去的磁通会“撞上”边界壁面反射，产生实际中不可能出现的磁通集中。想象一下房间太小，球在墙壁上不断弹跳的状态。
边界条件的思考方式
远场边界条件看似不起眼但超级重要。需要在数值上表达“从这里开始是无限延伸的空间”。如果设置错误，磁通就会像碰到“看不见的墙”一样被弹回。