[論文レビュー] Transport formulas for multi-component plasmas within the effective potential theory framework

Grigory Kagan, Scott Baalrud|arXiv (Cornell University)|Nov 29, 2016

Dust and Plasma Wave Phenomena参考文献 3被引用数 8

ひとこと要約

本稿では、有効ポテンシャル理論（EPT）を用いて、多数体効果を一般化されたクーロン対数を介して修正された二体衝突に還元することで、結合プラズマの高精度な流体モデル化を可能にする、コン act で解析的に取り扱いやすい多成分プラズマの輸送式を提示する。主な貢献は、弱い結合極限では完全な解析的表現が得られ、強い結合領域では数値ルーチンを備えた実用的な輸送係数計算フレームワークを提供することにある。

ABSTRACT

The recently proposed effective potential theory [Phys. Rev. Lett. 110, 235001 (2013)] allows evaluating transport in coupled plasmas with the well-developed formalisms for systems with binary collisions. To facilitate practical implementation of this concept in fluid models of multi-component plasmas, compact expressions for the transport coefficients in terms the generalized Coulomb logarithms are summarized from existing prescriptions. For weakly coupled plasmas, characterized by Debye-shielded Coulomb interaction potential, expressions become fully analytical. In coupled plasmas the generalized Coulomb logarithms need to be evaluated numerically. Routines implementing the described formalisms are included as supplemental material.

研究の動機と目的

• 複数のイオン種と電子を含む多成分プラズマへの有効ポテンシャル理論（EPT）の拡張。
• 結合プラズマの流体モデルにおける輸送係数を計算するための実用的で計算効率の良いフレームワークの提供。
• 分子動力学シミュレーションと流体モデルの間のギャップを埋めるために、修正されたポテンシャルを用いた二体衝突形式を用いて正確な輸送予測を可能にする。
• 弱い結合および強い結合プラズマの両方に対応する実装可能な数倉ルーチンの提供により、既存のプラズマシミュレーションワークフローへの統合を促進する。

提案手法

• 結合プラズマにおける多数体効果を、修正された相互作用ポテンシャルを用いて有効な二体衝突に還元する有効ポテンシャル理論（EPT）の応用。
• 輸送係数（拡散、粘性、熱フラックス）を、有効ポテンシャルの影響を内蔵した一般化されたクーロン対数 Ξαβ(l,k) を用いて表現。
• チャップマン＝エンスコーグ法およびグラド法を用いて線形化されたボルツマン方程式の解を得る。速度空間における直交多項式展開を用い、さまざまな近似レベル（ξ = 2, 3）での輸送係数を導出。
• デバイ遮蔽クーロンポテンシャルに基づく一般化されたクーロン対数を用いて、弱い結合極限における輸送係数の完全な解析的表現を導出。
• 事前に計算された Ξαβ(l,k) 値を用いて、強い結合プラズマにおける輸送係数を評価するためのMATLABにおける数値ルーチンを提供。
• 一成分プラズマおよび二成分イオンプラズマにおける動的摩擦、熱拡散、バロ拡散の既報結果を再現するための例題スクリプトを含む。

実験結果

リサーチクエスチョン

• RQ1 多数体効果が支配的となる、特に強い結合領域における多成分プラズマの輸送係数をどのように高精度にモデル化できるか。
• RQ2 有効ポテンシャル理論（EPT）は、複数のイオン種を含む結合プラズマにおいて、二体衝突形式の有効性をどの程度まで拡張できるか。
• RQ3 多成分プラズマの弱い結合極限における輸送係数（拡散、粘性、熱フラックス）の解析的表現は何か。
• RQ4 一般化されたクーロン対数 Ξαβ(l,k) はどのように計算され、熱拡散や動的摩擦といった既知の輸送現象をどのように回復できるか。
• RQ5 結合プラズマの流体モデルにおける効率的かつ高精度な輸送計算を実現するための数値実装戦略は何か。

主な発見

• 有効ポテンシャル理論により、多数体物理学を修正されたポテンシャルを用いた有効な二体衝突に還元することで、多成分プラズマにおける輸送係数の高精度な予測が可能になる。
• 弱い結合極限では、デバイ遮蔽クーロンポテンシャルに基づく一般化されたクーロン対数を用いて、すべての輸送係数が完全に解析的に表現可能である。
• 本形式を用いて計算された動的摩擦係数 $ A_{lh} $ および熱拡散比 $ k_T^{(i)} $ は、Kaganら（2014年）の報告結果を、イオン質量比および濃度比の異なる条件下でも正確に再現する。
• バロ拡散比 $ k_p $、電気拡散比 $ k_E $、電子温度依存の熱拡散比 $ k_T^{(e)} $ は、導出された行列表現から一貫して回復可能である。
• 数値ルーチン 'heat diffusion.m'、'viscosity.m'、'heat diffusion weak.m'、'viscosity weak.m' は、Kagan ら（2016年）が示したOCPおよびDT混合物の輸送結果を正確に再現した。
• 水素の1成分プラズマ（OCP）の一般化されたクーロン対数 Ξαβ(l,k) は、'xi lk ocp.mat' に格納されており、同じ電荷を有する種の相対濃度に依存しないため、水素同素体の混合物に対しても再利用可能である。