磁気ベクトルポテンシャル
理論と物理
ベクトルポテンシャル$\mathbf{A}$
先生、磁気ベクトルポテンシャルって何ですか？
磁束密度$\mathbf{B}$の「ポテンシャル」。$\mathbf{B}$は発散がゼロ（$\nabla \cdot \mathbf{B} = 0$）なので、必ずあるベクトル場$\mathbf{A}$の回転として表せる:
$$ \mathbf{B} = \nabla \times \mathbf{A} $$
静電場の電位$\phi$に対応するんですね。
対応関係:
先生、磁気ベクトルポテンシャルって何ですか？
磁束密度$\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の最大の利点。
まとめ
$\mathbf{A}$は一意に決まらない（任意のスカラー場$\psi$の勾配を加えても$\mathbf{B}$は変わらない）。一意性のためにゲージ条件を課す:
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ソルバの基礎となっている。
各項の物理的意味
- 電場項 $\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定式化
弱形式:
$\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ゲージは行列の条件数を改善し、大規模問題での収束性向上に貢献する。ゲージ固定の選択が計算の安定性と効率に直接影響する重要な数値技術だ。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
実践ガイド
実務
ほぼ全ての磁場FEM解析がA法（ベクトルポテンシャル法）に基づいている。JMAG、Maxwell、COMSOL、FEMMの内部ソルバーは全てA法。
チェックリスト
「A定式化と磁束密度の精度差」——エッジ要素vs節点要素の選択
磁場FEMでAを離散化するとき、節点要素（スカラー自由度）かエッジ要素（Nedelec要素）かの選択が精度を左右する。節点要素は接線成分の連続性を保証しないため、異なる透磁率の界面でBの法線成分の誤差が生じる。エッジ要素は接線成分の連続性を自然に満たし、界面での数値的ジャンプを抑えられる。現代の商用ソルバはほぼエッジ要素を採用しているが、ライセンスコスト・メモリの観点から節点要素を選ぶカスタム実装も存在する。精度が重要な永久磁石・高透磁率材料の界面ではエッジ要素が標準的な選択だ。
解析フローのたとえ
モータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ（コイル巻数）とブリッジの位置（磁石配置）を調整して、最も美しい音色（効率の良いトルク特性）を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。
初心者が陥りやすい落とし穴
「空気領域？ なんで空気をメッシュで切るの？」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。
境界条件の考え方
遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。
ソフトウェア比較
ツール
全磁場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定式化を採用しており、流体-電磁気連成問題の研究に活用されている。
選定で最も重要な3つの問い
- 「何を解くか」：磁気ベクトルポテンシャルに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端
AV定式化の拡張——時間高調波解析と複素ベクトルポテンシャル
交流磁場解析ではAを複素量として扱う「時間高調波（周波数領域）定式化」が広く使われる。A=Re(Â exp(jωt))として実部・虚部を同時に解くことで、定常周期解析が1解で完了する。ただし非線形BH曲線を持つ鉄心は複素Aでは厳密に扱えないため、「等価損失抵抗」近似か「時間ステッピング（FDTD的）」との使い分けが必要だ。多周波成分（インバータ高調波）を扱う場合は複数周波数の重ね合わせ計算や非線形時間ステッピングが必要で、計算コストは周波数数に比例して増大する。
トラブルシューティング
トラブル
「ゲージ未固定でFEMが収束しない」——典型的な初心者の罠
磁気ベクトルポテンシャルA定式化でゲージを固定しないと、解が一意に定まらず反復ソルバが収束しないかゆっくりしか収束しない。この問題は「特異行列（正則でない系）」として現れ、ソルバのエラーメッセージが「インデフィニット行列」や「ゼロ固有値検出」となる。対策：①ソルバのゲージ設定を確認する（多くは「Automatic Gauge」設定）、②境界条件でAの法線成分または接線成分を一か所以上固定する。商用ソルバはゲージ固定を自動処理するが、カスタムFEMコードでは意識的なゲージ実装が必要で、見落としがちな重要ポイントだ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——磁気ベクトルポテンシャルの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告