ビオ・サバールの法則
ビオ・サバールの法則の理論基礎
ビオ・サバールの法則
先生、ビオ・サバールの法則って何ですか？
電流が作る磁界を計算する基本法則。静電場のクーロンの法則に対応する磁場版。
$\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$ |
まとめ
ビオとサバール——2人は実は共同研究の相棒だった
「ビオ・サバールの法則」は2人の名前を冠しているが、彼らが共同で実験したのはわずか数週間という説もある。ジャン＝バティスト・ビオは実験物理学者、フェリックス・サバールは医師兼物理学者という異色のコンビで、1820年にエルステッドの発見から数ヶ月のうちに電流と磁場の定量的な関係を導いた。一方、ビオは後にその功績をめぐってアンペールと激しく対立したことでも知られている。科学史の裏側には、発見の優先権をめぐるドラマが常にある。
ビオ・サバールの法則の数値計算手法
FEMでの磁場解析
磁場もFEMで解くんですか？
直接ビオ・サバールを積分するのではなく、ベクトルポテンシャル$\mathbf{A}$を導入する。
$\mathbf{A}$が満たす方程式:
$\mathbf{J}$: 電流密度。これが静磁場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）が標準。節点要素ではなく辺にDOFを割り当てることで、$\mathbf{B}$の法線成分の連続性を自然に満たす。
まとめ
ビオ・サバール数値積分の「特異点問題」——場所によって答えが爆発する
ビオ・サバールの法則を数値実装するとき最大の罠は「評価点が電流素片に近すぎると分母がゼロに近づいて値が発散する」こと。実務では電流を細いワイヤとして扱うのではなく、有限の断面積をもつソリッドモデルとして扱うか、評価点を電流要素から適切に離すオフセット処理が必要だ。特に細いコイルの端部や折り返し部分での評価は要注意で、数値積分点の密度を部位ごとに変えるアダプティブ積分が効果的なんだ。
ビオ・サバールの法則の実務適用
実務
コイルの磁界設計、MRI磁石、電磁石、リレー、センサの磁界計算が主な適用。
チェックリスト
脳磁計（MEG）のセンサー配置設計にビオ・サバールが使われる
脳の神経活動を計測する脳磁計（MEG：磁気脳波計）は、神経細胞の微弱な電流（ピコアンペアオーダー）が作り出す磁場をSQUIDセンサーで検出する装置だ。その磁場計算にはビオ・サバールの法則が使われていて、センサーをどう配置すれば最も精度よく信号源を特定できるかを最適化するのに数値積分が欠かせない。静磁場解析の技術が最先端医療機器の設計を支えているわけで、「電磁気の基礎なんて現場で使わない」という思い込みは捨てた方がいい。
ビオ・サバールの法則のソフトウェア比較
ツール
| ツール | 特徴 |
|---|---|
| JMAG | モーター・アクチュエータの業界標準。日本発 |
| Ansys Maxwell | 自動適応メッシュ。3D静磁場 |
| COMSOL AC/DC | マルチフィジックス連成 |
| FEMM | 2D静磁場。無償。教育・初期検討 |
| Cedrat Flux | 電力機器向け |
ビオ・サバール計算に強い無償ツール——PyCOMSOLとBField3D
商用ツール不要でビオ・サバール積分を実行したいなら、PythonベースのBField3DやmagpylibといったOSSが実用になる。これらはコイル形状をNumpyの配列で定義すれば任意点の磁場ベクトルを返してくれる。簡単なソレノイドやヘルムホルツコイルの設計検討なら、ANSYSのライセンス料を使わずとも十分。ただし鉄心や磁性材料が絡む問題では透磁率の非線形性を扱えないので、最終的にはFEMツールに渡す前の「概算見積もり」ツールとして使うのが賢い使い分けだ。
ビオ・サバールの法則の先端研究
先端
宇宙の磁場計算もビオ・サバールで——ITER核融合炉の超伝導コイル設計
フランスで建設中の核融合実験炉ITER（国際熱核融合実験炉）のプラズマ閉じ込め磁場設計は、ビオ・サバールの法則の巨大な応用例だ。超伝導コイル18基が発生する磁束密度は最大5.3テスラ、コイルに流れる電流は68,000アンペア。任意点での磁場を求めるためにビオ・サバール積分を数値的に解くのだが、コイル形状が複雑な3次元曲線なので、線積分を数万点に離散化して計算する。先端研究ではGPUを使った並列化で計算時間を数百倍短縮している。
ビオ・サバールの法則のトラブル対応
トラブル
「磁場の方向が逆になった」——右手則の見落としで起きる古典トラブル
ビオ・サバール積分を実装したとき「磁場の符号が逆」になるミスは驚くほど多い。原因の大半は電流ベクトルdlと変位ベクトルrの外積の順序ミス——コードに書いた r × dl を dl × r と逆にしてしまうだけで磁場が反転する。チェック方法は単純で、無限長直線電流の解析解（B = μ₀I/2πr）と比較して符号と大きさの両方を検証すること。この検証テストを実装の最初にパスさせる習慣を付ければ、複雑な3次元コイルに入ってからのデバッグが格段に楽になる。
なった
詳しく
報告