ビオ・サバールの法則のソースを表示

{{出典の明記|date=2021年11月}}
'''ビオ・サバールの法則'''（ビオ・サバールのほうそく、{{lang-en-short|Biot–Savart law}}）とは[[電流]]の存在によってその周りに生じる[[磁場]]を計算する為の[[電磁気学]]における法則である。この法則は静電場に対する[[クーロンの法則]]に対応する。

この法則によって磁場は距離、方向、およびその電流の大きさなどに依存することが論じられる。この法則は静的な近似の元では[[アンペールの法則]]および[[磁場に対するガウスの法則]]と同等である。

[[1820年]]に[[フランス]]の物理学者[[ジャン＝バティスト・ビオ]]と[[フェリックス・サヴァール]]によって発見された。
<!-- 先人の意志を尊重して磁場Hのまま記述しているが、他言語版では磁束密度Bで記述しているものが多い
    ↑Hに対して成り立つ式なのだからHで書くのが当然である。Bで書くと磁化Mがあらわれる (H = B/μ_0 - M)
 -->

==概要==
微小な長さの'''電流要素''' ''I'' d'''''l''''' によって '''''r''''' 離れた位置に作られる微小な磁場 d'''''H''''' は
:<math>\mathrm{d}\boldsymbol{H} = \frac{I \mathrm{d}\boldsymbol{l} \times \boldsymbol{r}}{4 \pi r^3}</math>

で表される。ここで ''r'' := |'''''r'''''| である。

電流がある程度の幅をもって流れているとき（すなわち、太さ無限小の線でなく領域 ''V'' を占めているとき）、[[電流密度]] '''''j''''' を使った積分形で書く必要がある：
:<math>\boldsymbol{H} =\int_V \frac{\boldsymbol{j}(\boldsymbol{r}') \times (\boldsymbol{r} - \boldsymbol{r}')}{4 \pi |\boldsymbol{r} - \boldsymbol{r}'|^3} \mathrm{d}^3 \boldsymbol{r}'</math>
なお上式では左辺の磁場 '''''H''''' は微小量ではない。

この法則は積分を実行して初めて有効な値が出る、すなわち実験的検証が間接的にならざるを得ない欠点がある。

==歴史==
{{Expand-section|date=March 2011}}
[[1820年]][[4月]]、[[デンマーク]]の物理学者[[ハンス・クリスティアン・エルステッド]]は[[コペンハーゲン大学]]での講義中、電気回路をいじっていた時近くにあった方位磁石が北ではない方角を指し示していることに気が付き、電流と磁場の関係について数か月の研究の末、電流の磁気作用を発表した。
これを受け、ジャン・バティスタ・ビオとフェリックス・サバールは共同で実験を行い、この法則を発表するに至った。さらにこの数ヵ月後には[[フランソワ・アラゴー]]が電磁石の原理を、[[アンドレ・マリー・アンペール]]がアンペールの法則を発見している。これらの功績がエルステッドの発見から僅か一年以内のことであったのは驚くべきことである。

さらに3年後の1823年にスタージャンが実際に電磁石を作成し、24年にアラゴーは回転磁気を発見している。この[[1820年代|1820年からの数年間]]は科学史上重要な期間である。

==その他の形式==
===均一な電流===
電流''I'' が如何なる点においても一定の場合磁場'''''H''''' は、
:<math>\boldsymbol{H} = \frac{I}{4 \pi} \int \frac{\mathrm{d}\boldsymbol{l} \times \boldsymbol{r}}{r^3}</math>
となる。

===等速度運動する点電荷===<!--E, B は力場である-->
点電荷''q'' が一定の速度'''''v''''' で運動しているとき、[[特殊相対性理論]]と[[マクスウェルの方程式]]より以下の電束密度と磁場が与えられる<ref>Griffiths, David J. (1998). Introduction to Electrodynamics (3rd ed.). Prentice Hall. pp. 222–224, 435–440. ISBN 0-13-805326-X</ref>。
:<math>\begin{align}
\boldsymbol{D} &= \frac{q}{4\pi} \frac{1-\beta^2}{(1-\beta^2\sin^2\theta)^{3/2}}\frac{\boldsymbol{r}}{r^3}\\
\boldsymbol{H} &= \boldsymbol{v} \times \boldsymbol{D}
\end{align}</math>

ただし、''&beta;'' = ''v'' / ''c''、''&theta;'' は '''''v''''' と'''''r''''' のなす角であり、''c'' は光速度である。

''v'' が ''c'' に対して十分に小さい（<math>v^2 \ll c^2</math>）ときは、近似的に
:<math>\begin{align}
\boldsymbol{D} &=\frac{q}{4\pi}\ \frac{\boldsymbol{r}}{r^3} \\
\boldsymbol{H} &=\frac{q \boldsymbol{v}}{4\pi} \times \frac{\boldsymbol{r}}{r^3}
\end{align}</math>
と表すことができる。

これらの電束密度と磁場に関する式は、'''点電荷に対するビオ・サバールの法則'''と呼ばれ、[[1888年]]に[[オリヴァー・ヘヴィサイド]]によって導かれた。

==計算例==
{{Expand-section|date=March 2011}}

===無限に長い直線電流の周りの磁場===
[[画像:ビオ・サバール.jpg|thumb|right]]
ビオ・サバールの法則は[[積分]]することによりアンペールの法則の磁場と一致する。
例えば無限に長い直線電流であれば、図より、
:<math>r = \frac{a}{\sin \theta}, \qquad s = -\frac{a}{\tan \theta}</math>
したがって
:<math>\frac{\mathrm{d}s}{\mathrm{d}\theta} = \frac{a}{\sin^2 \theta}</math>
となるから、ビオ・サバールの法則を積分して、
:<math>\begin{align}
H &= \frac{1}{4 \pi}\int_{-\infty}^{\infty} \frac{I\sin \theta}{r^2}\, \mathrm{d}s\\
&=\frac{1}{4 \pi}\int_{0}^{\pi} \frac{I\sin \theta}{r^2}\frac{a}{\sin^2 \theta}\, \mathrm{d}\theta = \frac{I}{4 \pi}\int_{0}^{\pi}\frac{\sin \theta}{a}\, \mathrm{d}\theta = \frac{I}{4 \pi a}\left[ -\cos \theta \right]_0^\pi\\
&=\frac{I}{2 \pi a}
\end{align}</math>
を得る。これは[[アンペールの法則]]の磁場の大きさと一致する。

===立体角を用いた解析===
[[File:BiSa-Amp.PNG|thumb|right|300px|閉経路C1から回路C2を俯瞰する立体角を&Omega;とする。]]
以下のようにしてもビオ・サバールの法則からアンペールの法則が成り立つことを示すことができる<ref>{{cite|和書|title=朝倉物理学選書2 電磁気学 |author=伊藤敏夫 |ISBN=978-4-254-13757-6}}</ref>。

閉回路C<sub>1</sub>上の点Pから回路C<sub>2</sub>を俯瞰する立体角を ''&Omega;'' とする。ここで回路C<sub>1</sub>上を点Pから微小距離 d'''''s''''' だけ移動した点をP&prime;とすると、点P&prime;から回路C<sub>2</sub>を俯瞰する立体角''&Omega;''+d''&Omega;'' は、−d'''''s''''' だけ平行移動された回路を俯瞰する立体角と等しい。

このとき、回路上の微小長さ d'''''s'''''&prime; と平行移動した微小距離 −d'''''s''''' によって作られる面の面素ベクトル d'''''S''''' は

:<math>\mathrm{d}\boldsymbol{S} = \mathrm{d}\boldsymbol{s'} \times (-\mathrm{d}\boldsymbol{s}) = \mathrm{d}\boldsymbol{s} \times \mathrm{d}\boldsymbol{s'}</math>

であるが、点Pから回路上の微小長さ d'''''s'''''&prime; へ向かうベクトルを '''''r''''' とすると、点Pから面素 d'''''S''''' を見る立体角は
:<math>\frac{\mathrm{d}\boldsymbol{S} \cdot \boldsymbol{r}}{r^3} = \frac{(\mathrm{d}\boldsymbol{s} \times \mathrm{d}\boldsymbol{s'}) \cdot \boldsymbol{r}}{r^3}
= \frac{(\mathrm{d}\boldsymbol{s'} \times \boldsymbol{r})}{r^3} \cdot \mathrm{d}\boldsymbol{s}
</math>

と表すことができる。
これを d''s''&prime; に関して回路一周分[[線積分]]すれば立体角の変化d''&Omega;''を得ることができる。
:<math>\mathrm{d}\Omega\ = \left(\oint_{C_2} \frac{\mathrm{d}\boldsymbol{s'} \times \boldsymbol{r}}{r^3}\right) \cdot \mathrm{d}\boldsymbol{s}
</math>

回路上の微小電流要素が点Pに作る磁場はビオ・サバールの法則を積分して、
:<math>\boldsymbol{H} = \oint_{C_2} \frac{1}{4\pi} \frac{I\mathrm{d}\boldsymbol{s'} \times (-\boldsymbol{r})}{r^3}</math>

と得られるが、この両辺に d'''''s''''' を内積で乗じ、先の式を代入すると
:<math>\boldsymbol{H} \cdot \mathrm{d}\boldsymbol{s} = -\frac{I}{4\pi} \mathrm{d}\Omega</math>
の関係が得られる。
点Pが閉曲線C<sub>1</sub>上を一周するような''&Omega;''の変化は
:<math> \oint_{C_1} \mathrm{d}\Omega = -4\pi</math>
であるので、
:<math>\oint_{C_1} \boldsymbol{H} \cdot \mathrm{d}\boldsymbol{s} = I</math>
とする、アンペールの法則そのものが導かれる。

この場合、C<sub>1</sub>で囲む領域 ''D''（図赤色）の面積を ''S''<sub>1</sub>とすると、面S<sub>1</sub>に対する電流面密度の大きさ''j'' は
:<math>j = \frac{I}{S_1}</math>
となるが、例えば閉曲線C<sub>1</sub>が1周する間に回路C<sub>2</sub>が3周するような場合には電流面密度の大きさは 3''j'' 、閉曲線C<sub>1</sub>が2周する間に回路C<sub>2</sub>が1周するような場合には電流面密度の大きさは ''j''/2 である。このことを考慮すれば、
:<math>\oint_{C_1} \boldsymbol{H} \cdot \mathrm{d}\boldsymbol{s} = \int_{D} \boldsymbol{j} \cdot \mathrm{d}\boldsymbol{S_1}</math>
と書くことができる。

===円形電流の中心付近に於ける磁場===
[[ファイル:Bios2.PNG|thumb|right|500px|中心に生じる磁場]]
アンペールの法則を使った場合では求めることが難しい場合も、ビオ・サバールの法則を用いることで簡易に計算できる場合がある。例えば円形電流の中心付近に発生する磁場を求める場合がそうである。まず、右図のような半径 ''a'' の円周上P点に存在する電流 ''I'' によって、中心Oに生じる磁場について考える。

d'''''s''''' と '''''r''''' の為す角度を ''&phi;'' とおくと、図より
:<math>\mathrm{d}\boldsymbol{s} \perp \boldsymbol{r}</math>
となり、また
:<math>| \boldsymbol{r} | = a </math>
であるので、
:<math>\mathrm{d}H = \frac{1}{4 \pi} \frac{I\sin \phi\,\!}{a^2}\, \mathrm{d}s = \frac{1}{4 \pi} \frac{I\, \mathrm{d}s}{a^2}\, </math>
である。これを円周上で積分して、
:<math>H = \frac{1}{4 \pi} \frac{I}{a^2}\,\int_{0}^{2\pi a} \mathrm{d}s = \frac{1}{4 \pi} \frac{2\pi aI}{a^2}
= \frac{I}{2a} </math>
となる。

[[ファイル:Biot-savart.PNG|thumb|right|500px|中心よりzだけずれた位置に生じる磁場]]
次に、右図のようなOより面に垂直に ''z'' だけずれた位置Qに生じる磁場について考える。図より、
:<math>| \boldsymbol{r} | = \sqrt{a^2 + z^2}, \quad \angle \mathrm{OPQ} = \alpha </math>
である。

d'''''H''''' はビオ・サバールの法則より d'''''s''''' と '''''r''''' に垂直で、面に平行な成分 <math>\mathrm{d}H_\parallel = \mathrm{d}H\sin \alpha</math> は対称性により円周上を積分すると 0 になってしまうので、面に垂直な成分 <math>\mathrm{d}H_\perp = \mathrm{d}H\cos \alpha</math> のみを考えればよい。
:<math>\mathrm{d}H_\perp (z) = \mathrm{d}H(z)\cos \alpha = \frac{1}{4 \pi} \frac{I\mathrm{d}s}{a^2 + z^2} \frac{a}{\sqrt{a^2 + z^2}}</math>

ここで、d''s'' = ''a'' d''&theta;'' であることを用いて、
:<math>\begin{align}
H_\perp (z) &= \frac{1}{4 \pi} \frac{Ia^2}{(a^2 + z^2)^{3/2}} \int_{0}^{2\pi} \mathrm{d}\theta\\
&= \frac{1}{4 \pi} \frac{2\pi\,Ia^2}{(a^2 + z^2)^{3/2}} = \frac{Ia^2}{2(a^2 + z^2)^{3/2}} \end{align}</math>

ここで、''z'' = 0 とすれば、円の中心部に生じている磁場''H''<sub>O</sub> が得られる。即ち、
:<math>H_\mathrm{O} = H_\perp (0) = \frac{I}{2a} </math>
であり、これは先ほど求めたものに一致する。

==発散と回転==
{{Expand-section|date=March 2011}}
<!--- スペイン語版記事を翻訳する予定 --->
===発散===
ビオ・サバールの法則の両辺の[[発散 (ベクトル解析)|発散]]を取る。
:<math>\operatorname{div} \boldsymbol{H} = \frac{1}{4 \pi} \int_V \operatorname{div} \left(\boldsymbol{j} \times \frac{\boldsymbol{r}}{r^3}\right) \mathrm{d}^3 \boldsymbol{r}'</math>

ここで、ベクトル解析の恒等式より
:<math>\operatorname{div} \left(\boldsymbol{j} \times \frac{\boldsymbol{r}}{r^3}\right) = -\boldsymbol{j} \cdot \operatorname{rot} \frac{\boldsymbol{r}}{r^3}</math>

また、
:<math>\frac{\boldsymbol{r}}{r^3} = -\operatorname{grad} \frac{1}{r}</math>
なので、これを代入すると、
:<math>\operatorname{div} \boldsymbol{H} = \frac{1}{4 \pi} \int_V \boldsymbol{j} \cdot \operatorname{rot} \left(\operatorname{grad} \frac{1}{r}\right) \mathrm{d}^3 \boldsymbol{r}'
= 0</math>
を得る。これは、磁場に対するガウスの法則より導かれる結果に等しい。

===回転===
ビオ・サバールの法則の両辺の[[回転 (ベクトル解析)|回転]]を取る。
:<math>\operatorname{rot} \boldsymbol{H} = \frac{1}{4 \pi} \int_V \operatorname{rot} \left(\boldsymbol{j} \times \frac{\boldsymbol{r}}{r^3}\right) \mathrm{d}^3 \boldsymbol{r}'</math>

ここで、ベクトル解析の恒等式より
:<math>\operatorname{rot} \left(\boldsymbol{j} \times \frac{\boldsymbol{r}}{r^3}\right) = \left(\operatorname{div} \frac{\boldsymbol{r}}{r^3}\right) \boldsymbol{j} - (\operatorname{div} \boldsymbol{j}) \frac{\boldsymbol{r}}{r^3}</math>

また、
:<math>\operatorname{div} \frac{\boldsymbol{r}}{r^3}= 4\pi \delta (r) </math>
なので、
:<math>\begin{align}\operatorname{rot} \boldsymbol{H} &= \frac{1}{4 \pi} \int_V 4\pi \delta (r) \cdot \boldsymbol{j} \mathrm{d}^3 \boldsymbol{r}'\\
&= \boldsymbol{j}\end{align}</math>
が得られる。これはアンペールの法則そのものである。

ただし、この書き換えは静磁場でのみ有効であることに留意しなければならない。

==ベクトルポテンシャル==
{{Expand-section|date=March 2011}}
静磁場で、[[ベクトルポテンシャル]]が
:<math>\boldsymbol{B} = \operatorname{rot}\boldsymbol{A}</math>
と定義出来るとき、
:<math>\frac{\boldsymbol{r}-\boldsymbol{r}'}{|\boldsymbol{r}-\boldsymbol{r}'|^3} = -\operatorname{grad} \frac{1}{|\boldsymbol{r}-\boldsymbol{r}'|}</math>
であるので、ベクトルの恒等式、
:<math>\operatorname{rot} \frac {\boldsymbol{j_0(\boldsymbol{r}')}}{|\boldsymbol{r}-\boldsymbol{r}'|} = -\boldsymbol{j}_0(\boldsymbol{r}') \times \operatorname{grad} \frac {1}{|\boldsymbol{r}-\boldsymbol{r}'|}</math>

からビオ・サバールの法則はベクトルポテンシャル '''''A''''' によって
:<math>\boldsymbol{A} =\frac{\mu_0}{4 \pi} \int_V \frac{\boldsymbol{j}_0(\boldsymbol{r}')}{|\boldsymbol{r} - \boldsymbol{r}'|} \mathrm{d}^3 \boldsymbol{r}'</math>
と書き換えられることになる。ここで、<math>\mu_0</math>は[[真空の透磁率]]、<math>\boldsymbol{j}_0=\boldsymbol{j} + \boldsymbol{j}_\boldsymbol{M}</math>は[[電流密度]]+磁化電流密度である。

また、ビオ・サバールの法則は静電場におけるクーロンの法則に対応するものであるが、同様に電場の[[スカラーポテンシャル]] ''&phi;''
:<math>\boldsymbol{E} = -\operatorname{grad} \phi\,\!</math>
によって静電場におけるクーロンの法則、あるいは[[ガウスの法則]]を書き換えると
:<math>\phi\,\! =\frac{1}{4 \pi \varepsilon_0} \int_V \frac{\rho_0(\boldsymbol{r}')}{|\boldsymbol{r} - \boldsymbol{r}'|} \mathrm{d}^3 \boldsymbol{r}'</math>
となり対称性を見ることができる。ここで<math>\varepsilon_0</math>は[[真空の誘電率]]、<math>\rho_0=\rho + \rho_\boldsymbol{P}</math>は[[電荷密度]]+分極電荷密度である。

== マクスウェル方程式からの導出 ==
もちろん、電磁気学の法則なので[[マクスウェルの方程式|マクスウェル方程式]]から導出することができる。前節において電場と磁場のそれぞれがスカラーポテンシャルとベクトルポテンシャルによって記述され、しかもかたや電荷密度、かたや電流密度を距離の逆比で重み付けして積分するという形式になっていることをみたが、その対称性の説明もできる。

まずクーロンの法則をマクスウェル方程式から導出するプロセスを考える。これは静電場仮定でのポテンシャルと電場の関係

:<math> 
 \boldsymbol{E} = -\operatorname{grad} \phi\,\!
</math>

を用いて、ガウスの法則

:<math> 
 \nabla\cdot\boldsymbol{D}=\rho
</math>

を変形すると

:<math> 
 \nabla^{2}\phi=-\frac{\rho}{\varepsilon_0}
</math>

が得られる。これを<math>\phi</math>に対する微分方程式だと考えると、[[グリーン関数法]]によって解くことができる。すなわち、方程式の両辺をフーリエ変換して、各モードに対するウェイトの方程式として読み替え、ウェイトがわかったところでフーリエ逆変換によって解を構成する。解は前節の通りで

:<math>
 \phi\,\! =\frac{1}{4 \pi \varepsilon_0} \int_V \frac{\rho(\boldsymbol{r}')}{|\boldsymbol{r} - \boldsymbol{r}'|} \mathrm{d}^3 \boldsymbol{r}'
</math>

となる。本題のビオ・サバールの法則の場合はアンペール・マクスウェルの法則

:<math> 
 \operatorname{rot}\boldsymbol{H}=\boldsymbol{j}+\frac{\partial}{\partial t}\boldsymbol{D}
</math>

をやはり静電場仮定で右辺第二項を無視した上で、ポテンシャル表示

:<math>
 \boldsymbol{H} = \frac{1}{\mu_0} \operatorname{rot}\boldsymbol{A}
</math>

を適用する。すると

:<math> 
 \operatorname{rot}\operatorname{rot}\boldsymbol{A}=\mu_0\boldsymbol{j}
</math>

となるが、<math> 
 \operatorname{rot}\operatorname{rot}\boldsymbol{A}=\operatorname{grad}\operatorname{div}\boldsymbol{A}-\nabla^{2}\boldsymbol{A}
</math>および、クーロンゲージ<math> 
 \operatorname{div}\boldsymbol{A}=0
</math>を適用することで

:<math> 
 \nabla^{2}\boldsymbol{A}=-\mu_0 \boldsymbol{j}
</math>

となる。今度はベクトルになっているが、演算子がスカラーなので各ベクトル成分に対して独立に同じ方程式を立てているに過ぎず、また、式の形上もスカラーポテンシャルの時と係数を除いて同じなので同じ方法を使って構成すれば

:<math>
 \boldsymbol{A} =\frac{\mu_0}{4 \pi} \int_V \frac{\boldsymbol{j}(\boldsymbol{r}')}{|\boldsymbol{r} - \boldsymbol{r}'|} \mathrm{d}^3 \boldsymbol{r}'
</math>

を得ることになる。これに回転を取ればビオ・サバールの法則が得られる。

以上のことをまとめると、
* 考えているベクトル場<math> 
 \boldsymbol{H}
</math>が何らかのベクトル場<math> 
 \boldsymbol{A}
</math>の回転で記述できる(<math> 
 \operatorname{div}\boldsymbol{H}=0
</math>)
* ベクトル場<math> 
 \boldsymbol{H}
</math>の回転がベクトル場<math> 
 \boldsymbol{j}
</math>に比例する
* 空間次元が3であるとき(グリーン関数法の結果がこのような重み付けになる上で必要)
ビオ・サバールの法則が(比例係数を除いて)得られることになる。

== 流体力学 ==
前節で説明したように、数式上の性質さえ共通すればビオ・サバールの法則を得ることになる。それは具体的には[[流体力学]]が好例である。

流体力学における[[渦度]]は流速の回転として定義されている。

:<math>\begin{align} \boldsymbol{\Omega} &= \operatorname{rot} \boldsymbol{v} = \nabla \times \boldsymbol{v} \\
&= \begin{pmatrix}\dfrac{\partial v_z}{\partial y} - \dfrac{\partial v_y}{\partial z},& \dfrac{\partial v_x}{\partial z} - \dfrac{\partial v_z}{\partial x},& \dfrac{\partial v_y}{\partial x} - \dfrac{\partial v_x}{\partial y} \end{pmatrix}\end{align}</math>

したがって、渦度の具体的な場<math> 
 \boldsymbol{\omega}
</math>を知ることができ(それ自身なので比例係数1で比例しているとみなす)、<math> 
 \operatorname{div}\boldsymbol{v}=0
</math>を言える場合には、(もちろん空間次元3の領域で問題を扱うので)<math> 
 \boldsymbol{\omega}
</math>からベクトル場<math> 
 \boldsymbol{v}
</math>をビオ・サバールの法則によって記述できる。

具体的には

:<math> 
 \boldsymbol{v}=\frac{1}{4\pi}\int_{V}\mathrm{d}V \boldsymbol{\omega}\times\frac{\boldsymbol{r}}{r^{3}}
</math>

という結果を得ることができる。

上記結果を渦糸に適用してその微小部分の寄与を取り出すことで次のように記述されるビオ・サバールの法則をみる<ref>{{cite|和書 |author=小池勝 |title=流体機械工学 |publisher=コロナ社 |year=2009 |isbn=978-4-339-04474-4 |page=29}}</ref>。
:<math>dv = \frac{\Gamma\sin\theta \mathrm{d}s}{4\pi r^2}</math>
ただし変数の意味を以下のように読み替える。
* ''dv'' ：観測点で誘導される速度
* &Gamma;：[[渦糸]]まわりの[[循環 (流体力学)|循環]]
* ''ds'' ：渦糸の微小部分
* &theta;：''ds'' の方向とそこから観測点を結ぶ直線とのなす角
* ''r'' ：''ds'' と観測点の距離

==脚注==
{{脚注ヘルプ}}
=== 出典 ===
<references />

==関連項目==
* [[物理学]]
* [[電磁気学の年表]]
* [[気象庁地磁気観測所]] - 当法則が[[直流電化]]による悪影響の根拠となった。
* [[電磁ポテンシャル]]

{{Physics-stub}}

{{電磁気学}}
{{Normdaten}}

{{DEFAULTSORT:ひおさはあるのほうそく}} 
[[category:電磁気学]]
[[category:自然科学の法則]]
[[Category:物理学のエポニム]]