磁気ベクトルポテンシャル
磁気ベクトルポテンシャルの理論基礎
ベクトルポテンシャル$\mathbf{A}$
先生、磁気ベクトルポテンシャルって何ですか？
磁束密度$\mathbf{B}$の「ポテンシャル」。$\mathbf{B}$は発散がゼロ（$\nabla \cdot \mathbf{B} = 0$）なので、必ずあるベクトル場$\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の最大の利点。
まとめ
磁気ベクトルポテンシャル——「物理的実体のない量」がFEMを支える理由
磁気ベクトルポテンシャルA（B=rot A）は電場・磁場の「影の量」とも呼ばれ、直接測定はできないが計算では欠かせない。FEM磁場解析でBではなくAを未知数とする理由は、rot B=μ₀Jの数値的処理より div B=0（磁束の連続性）の自動満足が重要だからだ。A定式化では div B=0がrot(rot A)=0として自動的に保証される。1970年代にOliver C. ZienkiewiczらがFEMにA定式化を適用し、今や全ての主要磁場FEMソルバの基礎となっている。
磁気ベクトルポテンシャルの数値計算手法
A法のFEM定式化
弱形式:
$\mathbf{N}_i$: 辺要素の形状関数（Nédélec要素）。
2D定式化
2Dでは節点要素で$A_z$を離散化:
静電場のポアソン方程式と全く同じ形。$\varepsilon \to \nu$、$\phi \to A_z$、$\rho \to J_z$。
まとめ
ゲージ固定——AV定式化での「非一意性」を解決する数値技術
磁気ベクトルポテンシャルAはゲージ変換（A→A+grad φ）に対して不変でないため、一意に決まらない「ゲージの自由度」がある。FEMでAを解くためにはゲージを固定する必要があり、主な手法が「クーロンゲージ（div A=0）」と「木-余木（Tree-Cotree）ゲージ」だ。クーロンゲージはペナルティ法やラグランジュ乗数法で実装され、多くの商用ソルバで採用されている。Tree-Cotreeゲージは行列の条件数を改善し、大規模問題での収束性向上に貢献する。ゲージ固定の選択が計算の安定性と効率に直接影響する重要な数値技術だ。
磁気ベクトルポテンシャルの実務適用
実務
ほぼ全ての磁場FEM解析がA法（ベクトルポテンシャル法）に基づいている。JMAG、Maxwell、COMSOL、FEMMの内部ソルバーは全てA法。
チェックリスト
「A定式化と磁束密度の精度差」——エッジ要素vs節点要素の選択
磁場FEMでAを離散化するとき、節点要素（スカラー自由度）かエッジ要素（Nedelec要素）かの選択が精度を左右する。節点要素は接線成分の連続性を保証しないため、異なる透磁率の界面でBの法線成分の誤差が生じる。エッジ要素は接線成分の連続性を自然に満たし、界面での数値的ジャンプを抑えられる。現代の商用ソルバはほぼエッジ要素を採用しているが、ライセンスコスト・メモリの観点から節点要素を選ぶカスタム実装も存在する。精度が重要な永久磁石・高透磁率材料の界面ではエッジ要素が標準的な選択だ。
磁気ベクトルポテンシャルのソフトウェア比較
ツール
全磁場FEMソルバーがA法を実装。ユーザーが意識する必要はほとんどない。FEMMは$A_z$を直接可視化できるので教育的に有用。
ベクトルポテンシャル定式化ツール——ANSYS Maxwell vs COMSOL AC/DC
磁気ベクトルポテンシャルA（またはAV）定式化を採用する主要ツールはANSYS Maxwell（A-φ定式化）、COMSOL AC/DC Module（A定式化、エッジ要素標準装備）、Elmer FEM（オープンソース、Tree-Cotreeゲージ実装済）だ。Maxwell は大規模3D問題のAMGソルバ（代数的マルチグリッド）最適化が進んでいる。COMSOLはLiveLinkでMATLABと連携し、カスタムゲージや独自材料モデルの実装が容易だ。OpenFOAMのelectromagneticsモジュールもA定式化を採用しており、流体-電磁気連成問題の研究に活用されている。
磁気ベクトルポテンシャルの先端研究
先端
AV定式化の拡張——時間高調波解析と複素ベクトルポテンシャル
交流磁場解析ではAを複素量として扱う「時間高調波（周波数領域）定式化」が広く使われる。A=Re(Â exp(jωt))として実部・虚部を同時に解くことで、定常周期解析が1解で完了する。ただし非線形BH曲線を持つ鉄心は複素Aでは厳密に扱えないため、「等価損失抵抗」近似か「時間ステッピング（FDTD的）」との使い分けが必要だ。多周波成分（インバータ高調波）を扱う場合は複数周波数の重ね合わせ計算や非線形時間ステッピングが必要で、計算コストは周波数数に比例して増大する。
磁気ベクトルポテンシャルのトラブル対応
トラブル
「ゲージ未固定でFEMが収束しない」——典型的な初心者の罠
磁気ベクトルポテンシャルA定式化でゲージを固定しないと、解が一意に定まらず反復ソルバが収束しないかゆっくりしか収束しない。この問題は「特異行列（正則でない系）」として現れ、ソルバのエラーメッセージが「インデフィニット行列」や「ゼロ固有値検出」となる。対策：①ソルバのゲージ設定を確認する（多くは「Automatic Gauge」設定）、②境界条件でAの法線成分または接線成分を一か所以上固定する。商用ソルバはゲージ固定を自動処理するが、カスタムFEMコードでは意識的なゲージ実装が必要で、見落としがちな重要ポイントだ。
関連トピック
なった
詳しく
報告