磁矢量势
磁矢量势
理论与物理
矢量势$\mathbf{A}$
老师，磁矢量势是什么？
磁通密度$\mathbf{B}$的“势”。由于$\mathbf{B}$的散度为零（$\nabla \cdot \mathbf{B} = 0$），它总可以表示为某个矢量场$\mathbf{A}$的旋度：
$$ \mathbf{B} = \nabla \times \mathbf{A} $$
这对应于静电场中的电势$\phi$呢。
对应关系：
静电场 静磁场 $\mathbf{E} = -\nabla\phi$ $\mathbf{B} = \nabla \times \mathbf{A}$ $\phi$: 标量 $\mathbf{A}$: 矢量（3分量） 泊松方程: $\nabla^2\phi = -\rho/\varepsilon$ $\nabla^2\mathbf{A} = -\mu\mathbf{J}$（库仑规范）
规范条件
$\mathbf{A}$不是唯一确定的（加上任意标量场$\psi$的梯度，$\mathbf{B}$也不会改变）。为了保证唯一性，施加规范条件：
- 库仑规范: $\nabla \cdot \mathbf{A} = 0$ — 静磁场的标准
- 洛伦兹规范: $\nabla \cdot \mathbf{A} + \mu\varepsilon \partial\phi/\partial t = 0$ — 动态问题
2D情况下的$A_z$
对于2D问题（$z$方向均匀），$\mathbf{A} = A_z(x,y) \hat{z}$，变为标量问题。磁力线是$A_z$的等高线。这是2D磁场FEM的最大优势。
总结
- $\mathbf{B} = \nabla \times \mathbf{A}$ — 磁场FEM的基础
- 通过规范条件保证唯一性 — 库仑规范是标准
- 2D中是$A_z$的标量问题 — 磁力线 = $A_z$的等高线
Coffee Break 闲谈
磁矢量势——“没有物理实体的量”支撑FEM的原因
磁矢量势A（B=rot A）也被称为电场·磁场的“影子量”，虽然无法直接测量，但在计算中不可或缺。在FEM磁场分析中，将A而非B作为未知数的原因是，自动满足 div B=0（磁通连续性）比处理 rot B=μ₀J 的数值计算更为重要。在A公式化中，div B=0 作为 rot(rot A)=0 被自动保证。1970年代，Oliver C. Zienkiewicz 等人将A公式化应用于FEM，如今已成为所有主流磁场FEM求解器的基础。
各项的物理意义
- 电场项 $\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假设（截面分析）：电流方向均匀且可忽略端部效应时有效
- 各向同性假设：对于各向异性材料（如硅钢板的轧制方向等），需要定义方向相关的特性
- 不适用的情形：等离子体（电离气体）、超导体、非线性光学材料需要额外的本构关系
老师，磁矢量势是什么？
磁通密度$\mathbf{B}$的“势”。由于$\mathbf{B}$的散度为零（$\nabla \cdot \mathbf{B} = 0$），它总可以表示为某个矢量场$\mathbf{A}$的旋度：
$$ \mathbf{B} = \nabla \times \mathbf{A} $$
这对应于静电场中的电势$\phi$呢。
对应关系：
| 静电场 | 静磁场 |
|---|---|
| $\mathbf{E} = -\nabla\phi$ | $\mathbf{B} = \nabla \times \mathbf{A}$ |
| $\phi$: 标量 | $\mathbf{A}$: 矢量（3分量） |
| 泊松方程: $\nabla^2\phi = -\rho/\varepsilon$ | $\nabla^2\mathbf{A} = -\mu\mathbf{J}$（库仑规范） |
$\mathbf{A}$不是唯一确定的（加上任意标量场$\psi$的梯度，$\mathbf{B}$也不会改变）。为了保证唯一性，施加规范条件：
- 库仑规范: $\nabla \cdot \mathbf{A} = 0$ — 静磁场的标准
- 洛伦兹规范: $\nabla \cdot \mathbf{A} + \mu\varepsilon \partial\phi/\partial t = 0$ — 动态问题
2D情况下的$A_z$
对于2D问题（$z$方向均匀），$\mathbf{A} = A_z(x,y) \hat{z}$，变为标量问题。磁力线是$A_z$的等高线。这是2D磁场FEM的最大优势。
总结
- $\mathbf{B} = \nabla \times \mathbf{A}$ — 磁场FEM的基础
- 通过规范条件保证唯一性 — 库仑规范是标准
- 2D中是$A_z$的标量问题 — 磁力线 = $A_z$的等高线
Coffee Break 闲谈
磁矢量势——“没有物理实体的量”支撑FEM的原因
磁矢量势A（B=rot A）也被称为电场·磁场的“影子量”，虽然无法直接测量，但在计算中不可或缺。在FEM磁场分析中，将A而非B作为未知数的原因是，自动满足 div B=0（磁通连续性）比处理 rot B=μ₀J 的数值计算更为重要。在A公式化中，div B=0 作为 rot(rot A)=0 被自动保证。1970年代，Oliver C. Zienkiewicz 等人将A公式化应用于FEM，如今已成为所有主流磁场FEM求解器的基础。
各项的物理意义
- 电场项 $\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假设（截面分析）：电流方向均匀且可忽略端部效应时有效
- 各向同性假设：对于各向异性材料（如硅钢板的轧制方向等），需要定义方向相关的特性
- 不适用的情形：等离子体（电离气体）、超导体、非线性光学材料需要额外的本构关系
数值解法与实现
A法的FEM公式化
弱形式：
$$ \int_\Omega \nu (\nabla \times \mathbf{A}) \cdot (\nabla \times \mathbf{N}_i) \, d\Omega = \int_\Omega \mathbf{J} \cdot \mathbf{N}_i \, d\Omega $$
$\mathbf{N}_i$: 边单元的形函数（Nédélec单元）。
2D公式化
2D中，用节点单元离散化$A_z$：
$$ \int_\Omega \nu \nabla A_z \cdot \nabla N_i \, d\Omega = \int_\Omega J_z N_i \, d\Omega $$
与静电场的泊松方程形式完全相同。$\varepsilon \to \nu$、$\phi \to A_z$、$\rho \to J_z$。
总结
- 3D: 边单元（Nédélec） — 自然地处理规范条件
- 2D: 节点单元处理$A_z$ — 与泊松方程同形
Coffee Break 闲谈
规范固定——AV公式化中解决“非唯一性”的数值技术
磁矢量势公
この記事の評価
ご回答ありがとうございます！
参考に
なった
なった
もっと
詳しく
詳しく
誤りを
報告
報告