磁気流体力学
 | この記事は検証可能な参考文献や出典が全く示されていないか、不十分です。 出典を追加して記事の信頼性向上にご協力ください。（2011年7月） |
磁気流体力学または磁性流体力学（英語：magnetohydrodynamics）とは、電導性の流体を扱うように拡張された流体力学であって、電磁流体力学とも呼ばれ、またしばしばmagneto-hydro-dynamicsの頭文字をとってMHDと称せられる。
目次
1 総説
2 磁気流体力学の仮定
3 有用な諸概念
3.1 磁場の圧力と張力
3.2 理想MHD
3.3 磁力線と流体との凍り付き
3.4 応用例 アルヴェーン波
4 関連項目
5 参考文献
6 外部リンク
総説
磁気流体力学の基本的アイデアは、電導性流体の中では流体の運動が磁場の変化をもたらして電流を誘起し、その電流と磁場との相互作用から流体への力を生じ、よって流体の運動自身が変化する、というものである。対象とする物質は主に液体金属（水銀など）とプラズマである。 そして基礎方程式として通常の流体力学の基礎方程式（ナビエ-ストークス方程式と連続の式）と電磁場のマクスウェルの方程式とを組み合わせて用いる。
磁気流体力学は 1942年に宇宙の諸現象研究の過程でハンス・アルヴェーンが発表した論文、すなわち今日アルヴェーン波として知られている磁場中電導性流体特有の波の存在を述べた論文から始まった。そしてアルヴェーン自身を含む多くの人々の研究により大きく発展し、今日では宇宙空間物理学研究や熱核融合研究の基礎として広く用いられている。アルヴェーンは「電磁流体力学の基礎研究、プラズマ物理学への応用」により1970年にノーベル物理学賞を受賞した。
磁気流体力学の仮定
磁気流体力学では、使用の実態に即して、近似として通常 次の仮定がなされる。
まず、ここで扱うのは電気伝導度の相当よい流体であるから、導体の電気力学に倣って変位電流を無視し、磁場 B と電流 j とはアンペールの法則で結び付けられているとする。 すなわち
rotB=μj{displaystyle {mbox{rot}}{boldsymbol {B}}=mu {boldsymbol {j}}} :（1）
ここで μ{displaystyle mu } は流体の透磁率で、定数と仮定されている。
ついで流体はほぼ中性とし、電荷を流体が運ぶことで生ずる対流電流は伝導電流と比較して小さいとして無視し、電流は伝導電流のみであるとする。そしてそれはオームの法則により定まるとする。すなわち
j=σ(E+v×B){displaystyle {boldsymbol {j}}=sigma left({boldsymbol {E}}+{boldsymbol {v}}times {boldsymbol {B}}right)} ：（2）
ここで σ{displaystyle sigma } は流体の電気伝導度である。ただし、この仮定は電荷密度ρe{displaystyle rho _{e}} を 0 とすることではない。実際、上記2つの式からE を求めてガウスの式 ρe=div(ϵE){displaystyle rho _{e}={mbox{div}}(epsilon {boldsymbol {E}})} に代入すれば 0 でない ρe{displaystyle rho _{e}} が求まる。ここでϵ{displaystyle epsilon } は流体の誘電率である。
また、これら2つの式を用いて、電磁誘導に関するマクスウェルの式 rotE+∂B/∂t=0{displaystyle {mbox{rot}}{mathit {E}}+partial {boldsymbol {B}}/partial t=0} から電場 E を消去すると、次の誘導方程式が得られる。
∂B∂t=rot(v×B)+1σμΔB{displaystyle {frac {partial {boldsymbol {B}}}{partial t}}={mbox{rot}}({boldsymbol {v}}times {boldsymbol {B}})+{frac {1}{sigma mu }}Delta {boldsymbol {B}}} ：（3）
ここでrot rotB=grad divB−ΔB=−ΔB{displaystyle {mbox{rot}} {mbox{rot}}{boldsymbol {B}}={mbox{grad}} {mbox{div}}{boldsymbol {B}}-Delta {boldsymbol {B}}=-Delta {boldsymbol {B}}} を用いた。また σ{displaystyle sigma } は定数と仮定した。
なおまた、運動方程式においても電荷密度が小さいので、電場による力は省略して、流体に及ぼす力は磁場による力 j×B{displaystyle {boldsymbol {j}}times {boldsymbol {B}}} のみであるとする。
こうして、運動方程式と方程式（1）および（3）により方程式系は変数 v{displaystyle {boldsymbol {v}}} と B{displaystyle {boldsymbol {B}}} とで閉じて、見かけ上、電場 E{displaystyle {boldsymbol {E}}} は完全に消去される。これが磁気流体力学の名の起こりである。しかし、電流を規定するという意味で電場の役割は本質的で、具体的問題では電場を考慮せずには解くことができないこともある。
有用な諸概念
磁場の圧力と張力
この節では都合により、ベクトル微分演算子 ∇{displaystyle nabla } を用いる記法で記述する。
流体に働く力 f=j×B{displaystyle {boldsymbol {f}}={boldsymbol {j}}times {boldsymbol {B}}} は j{displaystyle {boldsymbol {j}}} に（1）式を代入し、変形すると
f=−∇(B22μ)+1μ(B⋅∇)B{displaystyle {boldsymbol {f}}=-nabla left({frac {B^{2}}{2mu }}right)+{frac {1}{mu }}left({boldsymbol {B}}cdot nabla right){boldsymbol {B}}} ：（4）
となる。ここで右辺第1項は磁場の等方的圧力 （B2/2μ{displaystyle B^{2}/2mu }） による力と解釈出来る。また第2項はその点での磁場方向の単位ベクトルを b{displaystyle {boldsymbol {b}}}、その方向を z方向にとると
1μ(B⋅∇)B=∂∂z(B22μ)b+B2μ∂∂zb{displaystyle {frac {1}{mu }}left({boldsymbol {B}}cdot nabla right){boldsymbol {B}}={frac {partial }{partial z}}left({frac {B^{2}}{2mu }}right){boldsymbol {b}}+{frac {B^{2}}{mu }}{frac {partial }{partial z}}{boldsymbol {b}}} ：（5）
となり、第1項は磁場方向の張力で、（4）式第1項の z 成分と打ち消し合う。また第2項は磁力線が湾曲しているとき、
∂∂zb=1ReR{displaystyle {frac {partial }{partial z}}{boldsymbol {b}}={frac {1}{R}}{boldsymbol {e}}_{R}} ：（6）
（ここで R は磁力線のその点での曲率半径、eR{displaystyle {boldsymbol {e}}_{R}} はその点から曲率中心へ向かう単位ベクトル）となることから解る通り、大きさ （B2/μ{displaystyle B^{2}/mu }） の張力により流体を曲率中心の方向へ引っ張る力である。
かくして結局、流体に作用する力は
f=−∇⊥(B22μ)+B2μ1ReR{displaystyle {boldsymbol {f}}=-nabla _{perp }left({frac {B^{2}}{2mu }}right)+{frac {B^{2}}{mu }}{frac {1}{R}}{boldsymbol {e}}_{R}} ：（7）
（∇⊥{displaystyle nabla _{perp }} は磁力線に直角な2次元微分演算子）と書くことが出来て、磁力線に直角に働く大きさ (B2/2μ{displaystyle B^{2}/2mu }) の磁場の圧力と、磁力線の湾曲を引き延ばすように働く大きさ(B2/μ{displaystyle B^{2}/mu })の磁場の張力とからなる、と言うことが出来る。
理想MHD
MHDで、電気伝導度 σ{displaystyle sigma } を十分大きいとして σ→∞{displaystyle sigma rightarrow infty } とおいたものを 理想MHD(ideal MHD) といい、広い応用をもつ。 そこではオームの法則（2）は電場を E=−v×B{displaystyle {boldsymbol {E}}=-{boldsymbol {v}}times {boldsymbol {B}}} と定めるのに用いられ、電流 j{displaystyle {boldsymbol {j}}} はアンペールの法則（1）により定まるものとされる。そして誘導の方程式（3）は右辺の第2項が省略出来て
∂B∂t=rot(v×B){displaystyle {frac {partial {boldsymbol {B}}}{partial t}}={mbox{rot}}({boldsymbol {v}}times {boldsymbol {B}})} ：（8）
となる。理想MHDが使えるのは、式（3）右辺の第2項が第1項に比して十分小さく、（8）が正しくなる場合に限られる。
式（3）の右辺第1項と第2項との大きさの比は、流れを表す代表的な速度、大きさをそれぞれ U,L とすると、おおよそ
Rm=σμUL{displaystyle R_{m}=sigma mu UL} ：（9）
となり、磁気レイノルズ数( magnetic Reynolds number )と呼ばれる。そして Rm ≫ 1 が理想MHDの適用出来る条件である。通常の実験室内の電導性流体ではRm は高々1の程度の大きさの量であるが、宇宙プラズマでは L が非常に大きいため、また超高温プラズマでは σ{displaystyle sigma } が非常に大きいために、いずれも Rm →∞ の仮定が十分正しくなり、理想MHDが適用される。
磁力線と流体との凍り付き
理想MHDのもっとも顕著な特色は磁力線と流体との凍り付きである。ある時刻の磁力線はその時の磁場分布だけで定まり、異なる時刻の磁力線同士を関係付ける方法は一般には存在しない。しかし、理想MHDのもとでは方程式（8）のおかげで、ある点の磁力線はそこでの流体の速度 v{displaystyle {boldsymbol {v}}} で動く、すなわち磁力線は流体と一緒に動くとする扱いが許される。この現象を凍り付き(froze in)と言う。その結果、磁力線は流体により運ばれて時間とともに移動していく（対流）、もしくは磁力線は流体を凍り付かせて質量密度をもつ実体として運動する、という描像を画くことが出来る。
一方、もとの方程式（3）に戻れば、電気伝導度 σ{displaystyle sigma } が有限な場合は右辺第2項は磁場B{displaystyle {boldsymbol {B}}} の拡散を表す。すなわち、σ{displaystyle sigma } が有限な場合は磁力線と流体との凍り付きは完全ではなく、磁場は流体中を拡散する。そしてその緩和時間は tc=σμL2{displaystyle t_{c}=sigma mu L^{2}} で与えられる。
応用例 アルヴェーン波
一般に線密度 ρ{displaystyle rho }、張力 T をもつ糸には速さ T/ρ{displaystyle {sqrt {T/rho }}} の横波が伝播する。上により、磁力線は流体に凍り付いているから、単位断面積の磁力管を考えると、それは線密度 ρm{displaystyle rho _{m}}の流体の糸と見なせる（ρm{displaystyle rho _{m}} は流体の質量密度）。そしてそれに大きさ B2/μ{displaystyle B^{2}/mu } の張力がかかっているから、磁気流体中には磁力線に沿って速さ
vA=Bρmμ{displaystyle v_{A}={frac {B}{sqrt {rho _{m}mu }}}}
の横波が伝播する。これがアルヴェーン波である。
関連項目
- MHD発電
- プラズマ宇宙論
- ダイナモ理論
参考文献
今井功・桜井明『電磁流体力学』岩波講座 現代物理学V.H.、岩波書店、1959年、
外部リンク
Magnetohydrodynamics （英語） - スカラーペディア百科事典「磁気流体力学」の項目。
