プラズマのモデリングのソースを表示

'''プラズマのモデリング''' または '''プラズマのモデル化''' とは、[[プラズマ]]の状態を記述する[[運動方程式]]を解く事を指す。
[[電磁場]]を記述する[[マクスウェルの方程式|マクスウェル方程式]]、もしくは静電場を記述する[[ポアソン方程式]]も連成して解かれる場合が多い。
プラズマは多数の粒子([[電子]]、[[イオン (化学)|イオン]]、電気的中性な[[分子]]、[[ラジカル (化学)|ラジカル]]等)から構成される。
それらをモデル化する方法にはいくつかの種類があり、それぞれ利点と制限を持っている<ref>
{{cite journal 
 | author1 = L. L. Alves et. al.
 | title = Foundations of modelling of nonequilibrium low-temperature plasmas
 | journal = Plasma Sources Sci. Technol.
 | volume = 27
 | year = 2018
 | issue = 2
 | page = 023002
 | doi = 10.1088/1361-6595/aaa86d
 }}
</ref>。
[[Image:A Comparison Chart For Modeling Plasma2.png|400px|right|Chart for modeling plasma]]
{{Clear|right}}
==試験粒子モデル==
試験粒子モデルは、印加された[[電場]](電界)と[[磁場]](磁界)の中を電子とイオンが独立的に動くものとしてプラズマを記述する。ここで電場と磁場は電子やイオンに対して[[運動の第3法則|外力]]として作用し、荷電粒子間の相互作用は考慮されない。そのため、この手法は[[セルフコンシステント]]ではない。これらの理想的な状態にある電子やイオンは '''試験粒子''' ('''テストパーティクル''') と呼ばれる<ref>
{{Cite book|和書
|editor = 日本機械学会
|title = 原子・分子モデルを用いる数値シミュレーション
|year = 1996
|publisher = コロナ社
|series = コンピュータアナリシスシリーズ
|issue = 7
|isbn = 4-339-04141-6
}}
</ref>。

試験粒子は[[ローレンツ力]]に従って旋回運動をするが、多くの場合、比較的速い[[円運動]]と旋回中心の比較的遅い[[ドリフト (プラズマ物理学)|ドリフト]]の重ね合わせとして扱うことが可能である。

==運動論モデル==
運動論モデルは、プラズマを記述する最も基礎的な方法であり、粒子の位置 <math>\vec{x}</math> と速度 <math>\vec{v}</math> の6次元[[位相空間 (物理学)|位相空間]] <math>(\vec{x},\vec{v})</math> における '''分布関数''' <math>f(\vec{x},\vec{v},t)</math> でプラズマを記述する。ここで、位置 <math>\vec{x}</math> と速度 <math>\vec{v}</math> は独立変数である。分布関数の[[基礎方程式|支配方程式]]として、次のいずれかが解かれる。

* [[ボルツマン方程式]]
粒子間相互作用を2体衝突で近似した衝突項 (ボルツマンの衝突項) が置かれる。衝突項には複雑な多重積分が含まれており、近似された衝突モデルか[[モンテカルロ法]]が使用される。長距離にまで影響が及ぶ[[クーロンの法則|クーロン相互作用]]を衝突項で考慮すると、厳密には不正確となる。そこで、荷電粒子間の相互作用や多体衝突効果を正確に考慮する場合は、集団運動が作る電磁場の作用として扱う<ref name="Lieberman">
{{cite book |和書
|author1= M. A. Lieberman
|author2= A. J. Lichtenberg
|others= 堀勝 監修、佐藤久明 訳
|publisher= 丸善出版
|title= プラズマ/プロセスの原理
|edition= 第2版
|year= 2010
|month= 1
|origtitle= Principles of Plasma Discharges and Materials Processing
|origyear= 2005
|isbn= 978-4-621-08223-2
}}
</ref><ref name="IEEJ">
{{cite book |和書
|author= 関口忠
|series= 電気学会大学講座
|publisher= オーム社
|title= プラズマ工学
|year= 1997
|isbn= 4-88686-220-9
}}
</ref>。

* [[運動論的方程式#ブラソフ方程式|ブラゾフ方程式]]
ボルツマン方程式において衝突項を0と置いた式である。近距離の粒子間相互作用の影響が小さい系に対して適しており、[[天体物理学]]等の分野ではよく使用される。

* [[フォッカー・プランク方程式]]
衝突による小角散乱を[[マルコフ過程]]で扱った衝突項 (フォッカー・プランクの衝突項) が置かれる。この衝突項は[[移動現象論|輸送現象]]や[[緩和現象]]の扱いに適している。

分布関数より得られる[[電荷]]と[[電流]]は、マクスウェル方程式を介して電磁場をセルフコンシスタントに定める。

数値シミュレーションを実施する上では、分布関数を6次元位相空間で離散化する必要がある。この離散化には[[粒子法|粒子的な手法]]または[[連続体力学|連続体]]的な手法が使用される<ref>
{{Cite journal|和書
|author1 = 内藤裕志
|author2 = 佐竹真介
|title = 粒子シミュレーションのコーディング技法
|date = 2013-04-25
|journal = プラズマ・核融合学会誌
|volume = 89
|issue = 4
|pages = 245-260
|issn = 09187928
|naid = 110009604189
}}
</ref><ref>
{{Cite journal|和書
|author = 真壁利明
|title = 反応性非平衡プラズマのモデリング
|year = 1991
|journal = 応用物理学会
|volume = 60
|issue = 7
|pages = 663-673
|url=https://doi.org/10.11470/oubutsu1932.60.663
|CRID=1390001204595199872
|doi=10.11470/oubutsu1932.60.663
}}
</ref>。プラズマ解析に使用される代表的な粒子法には[[PIC法|PIC]]-MCC法がある<ref>
{{Cite journal 
 | author = Birdsall, C.K.
 | title = Particle-in-cell charged-particle simulations, plus Monte Carlo collisions with neutral atoms, PIC-MCC 
 | journal = IEEE Transactions on Plasma Science
 | volume = 19
 | issue = 2
 | year = 1991
 | pages = 65&ndash;85
 | doi = 10.1109/27.106800
 | issn = 0093-3813
 | bibcode = 1991ITPS...19...65B 
}}
</ref>。

==ジャイロ運動論モデル==
ジャイロ運動論モデルは、粒子の円運動をジャイロ周期に渡って平均化する。このモデルはバックグラウンドに強い磁場を伴う系に対して適切であり、トカマクプラズマが持つ不安定性のシミュレーションに対して広く使用されている<ref>
{{Cite journal|和書
|author = 洲鎌英雄
|title = ジャイロ運動論
|year = 2003
|journal = プラズマ・核融合学会誌
|volume = 79
|issue = 2
|pages = 107-120
|url=https://doi.org/10.1585/jspf.79.107
|doi=10.1585/jspf.79.107
|CRID=1390001206513415552
}}
</ref>。また最近では天体物理学における[[オーロラ (代表的なトピック)|オーロラ]]解析等でも利用されている。

==流体モデル==
運動論モデルはプラズマの物理現象を正確に記述するが、多くの[[計算複雑性理論|複雑性]] (および数値シミュレーションでは多くの計算負荷) を伴う。流体モデルでは '''巨視量''' ([[密度]]、[[流速|平均速度]]、平均エネルギー等の、分布関数の速度空間における[[積率]])を基礎変数にしてプラズマを記述することで、計算をより単純にする。これらの巨視量に対する方程式は、いわゆる流体方程式と呼ばれ、[[運動論的方程式]]の速度空間における積率を取る事で得られる。

流体方程式は、[[輸送係数]] ([[移動度]]や[[拡散係数]]等) の定め無しでは閉じられていない。輸送係数を定めるためには、分布関数をいくつかある候補の中から選択して仮定する必要がある。[[熱力学的平衡|熱平衡]]が実現している系に対しては、[[マクスウェル分布]]は1つの良い選択である。しかし非平衡性が強くなるプラズマの系では、粒子の速度分布がマクスウェル分布から大きく変化する事が知られており、マクスウェル分布を仮定するといくつかの矛盾を抱える事がある<ref name="Lieberman" /><ref name="IEEJ" /><ref>
{{cite journal 
 | author1 = J. van Dijk et. al.
 | title = Plasma modelling and numerical simulation
 | journal = J. Phys. D: Appl. Phys.
 | volume = 42
 | year = 2009
 | issue = 19
 | page = 190301
 | doi = 10.1088/0022-3727/42/19/190301
 }}
</ref><ref>
{{Cite book|和書
|editor = 電気学会・マイクロ波プラズマ調査専門委員会
|title = マイクロ波プラズマの技術
|year = 2003
|publisher = オーム社
|isbn = 4-274-19710-7
}}
</ref>。

流体方程式と電磁場方程式との結合系を考える理論体系は、[[磁気流体力学|電磁流体力学]]として知られている。

==運動論と流体のハイブリッドモデル==
ハイブリッドモデルでは、系の粒子種のいくつかを流体として、その他を運動論的に扱う。流体で扱う粒子種と運動論的に扱う粒子種、および電磁場がセルフコンシステントに解かれるように構成される必要がある。

==プラズマ解析専用のシミュレータ製品==
===運動論ベース===
*[https://www.northropgrumman.com/space/pic-code-software/ LSP] - ''www.northropgrumman.com''
*[https://www.northropgrumman.com/space/pic-code-software/ Magic] - ''www.northropgrumman.com''
*[https://www.wavefront.co.jp/CAE/particle-plus/ Particle-PLUS] - ''www.wavefront.co.jp''

===流体ベース===
*[http://esgeetech.com/products/vizglow/ VizGlow] - ''esgeetech.com''
*[http://esgeetech.com/products/vizspark-thermal-plasma-modeling/ VizSpark] - ''esgeetech.com''

また、上記に挙げたプラズマ解析専用の[[シミュレーション|シミュレータ]]の他に、[[数値流体力学#著名な数値流体力学ソフトウェア|汎用的な流体解析シミュレータ]]にプラズマ解析機能が含まれるものもある。

== 教科書 ==
*Francis F. Chen, 'Introduction to Plasma Physics and Controlled Fusion', 3rd edition, Springer (2015), {{ISBN2|978-3-319-22308-7}}
*Nicholas Krall and Alvin Trivelpiece, 'Principles of Plasma Physics', San Francisco Press (1986), {{ISBN2|978-0-911302-58-5}}
*M. A. Lieberman and A. J. Lichtenberg「プラズマ/プロセスの原理」第2版、堀勝 監修、佐藤久明 訳、丸善出版 (2010年)、原著 'Principles of Plasma Discharges and Materials Processing'、{{ISBN2|978-4-621-08223-2}}
*宮本健郎「核融合のためのプラズマ物理」『NIFS-PROC-80』シリーズ、核融合科学研究所 (2010年)、{{ISSN|0915-6348}}
*関口忠「プラズマ工学」『電気学会大学講座』シリーズ、オーム社 (1997年)、{{ISBN2|4-88686-220-9}}

== 脚注 ==
{{Reflist}}

== 関連項目 ==
*[[プラズマ物理]]
*[[磁気流体力学]]
*[[マルチフィジックス]]
*[[計算科学]]
*[[数値解析]]
*[[数値流体力学]]

{{DEFAULTSORT:ふらすまのもてりんく}}
[[Category:プラズマ]]
[[Category:科学モデリング]]
[[Category:CAE]]
[[Category:計算物理学]]
[[Category:数値流体力学]]