[論文レビュー] Entropy Decay Rates for Conservative Spectral Schemes Modeling Fokker-Planck-Landau Type Flows in the Mean Field Limit

Clark A. Pennie, Irene M. Gamba|arXiv (Cornell University)|Oct 7, 2019

Gas Dynamics and Kinetic Theory参考文献 21被引用数 4

ひとこと要約

本稿では、クーロン、マクスウェル、および硬い球体粒子相互作用をモデル化するフォッカー＝プランク＝ランダウ方程式におけるエントロピーの減衰をシミュレートするため、保存則を満たすハイブリッドスペクトル-不連続ガラーキンスキームを提示する。クーロン相互作用の場合は理論的な2/3乗則に近い減衰率を示す退化スペクトルを示す一方、マクスウェルおよび硬い球体の場合はスペクトルギャップを示し、解析的予測を確認するとともに、クーロン相互作用の推定値の鋭さを検証する。

ABSTRACT

The focus of this work is to create benchmark simulations of decay rates to statistical equilibrium in transport plasma models for Coulomb particle interactions given by a coupled Vlasov-Poisson Fokker-Planck-Landau equation, as well as with Maxwell type and hard sphere interactions. The qualitative decay to the equilibrium Maxwell-Boltzmann distribution through relative entropy is studied in detail for all three types of particle interactions by means of a conservative hybrid spectral and discontinuous Galerkin scheme adapted from previous work. More precisely, the Coulomb case shows that there is a degenerate spectrum, with a decay rate close to the law of two thirds predicted by upper estimates in a work of Strain and Guo in 2006, while the Maxwell type and hard sphere examples both exhibit a spectral gap as predicted by Desvillettes and Villani in 2000. Such decay rate behavior indicates that the analytical estimates for the Coulomb case is sharp while, still to this date, there is no analytical proof of sharp degenerate spectral behaviour for the Fokker-Planck-Landau operator. Simulations are presented, both for the space-homogeneous case of just particle potential interactions and the space-inhomogeneous case for the mean field coupling through the Poisson equation for total charges in periodic domains. New explicit derivations of spectral collisional weights are presented in the case of Maxwell type and hard sphere interactions and the stability of all three scenarios, including Coulomb interactions, is investigated.

研究の動機と目的

• 平均場極限におけるヴラソフ＝ポアソン・フォッカー＝プランク＝ランダウ方程式の長時間ダイナミクスをシミュレートする保存則を満たす数値スキームの開発。
• クーロン、マクスウェル、および硬い球体の3種類の粒子相互作用について、最大エントロピー分布への相対エントロピーの減衰率の調査。
• 特にクーロン相互作用における退化スペクトルの解析的エントロピー減衰推定値の鋭さの検証。
• 最小化による衝突インヴァリアントの保存により、長時間にわたる数値的精度と安定性を確保。
• 空間均一な結果を、周期的領域におけるポアソン結合付きの空間非均一系へと拡張。

提案手法

• スペクトル法による衝突演算子の評価とDG法による移流を組み合わせた保存則を満たすハイブリッドスペクトルおよび不連続ガラーキンスキームを採用。
• スペクトル法により、O(N³log N)の計算量で衝突演算子のフーリエ変換を計算し、O(N³)の補正ステップにより保存則を強制。
• ボルツマン方程式に関するGambaら（2017）の質量集中結果に基づき、速度空間のカットオフを適用。
• マクスウェル型および硬い球体相互作用のためのスペクトル衝突重みの明示的導出を提供。
• 全質量、運動量、エネルギーの保存を維持することで、長時間安定性と正しい平衡ダイナミクスを保証。
• 空間均一および空間非均一な設定の両方で数値的シミュレーションを実施し、自己自己場を生成するためのポアソン方程式への結合を含む。

実験結果

リサーチクエスチョン

• RQ1 保存則を満たすスペクトルスキームは、クーロン相互作用を有するフォッカー＝プランク＝ランダウ方程式において、エントロピーの平衡への減衰を正確に捉えられるか？
• RQ2 平均場極限において、クーロン、マクスウェル型、および硬い球体相互作用の間でエントロピー減衰率はどのように比較されるか？
• RQ3 クーロン相互作用の予測された2/3乗則の減衰率は数値的に確認され、解析的推定値の鋭さを示唆するか？
• RQ4 マクスウェル型および硬い球体相互作用は、DesvillettesとVillani（2000）の予測通りにエントロピー減衰にスペクトルギャップを示すか？
• RQ5 このスキームは長時間シミュレーションにおいてインヴァリアントの保存を維持でき、正しい平衡への到達を保証するか？

主な発見

• クーロン相互作用の場合は、理論的な2/3乗則に近い減衰率を示す退化スペクトルを示し、StrainとGuo（2006）の上界推定値の鋭さを支持する。
• マクスウェル型および硬い球体相互作用の場合は、シミュレーションによりエントロピー減衰にスペクトルギャップが存在することが確認され、DesvillettesとVillani（2000）の予測と整合的である。
• 保存則を満たすスペクトルスキームは、長時間にわたるシミュレーションにおいて質量、運動量、エネルギーの保存を維持し、マクスウェル＝ボルツマン平衡の正確な近似を可能にする。
• 空間均一および空間非均一な両ケースにおいて、ポアソン結合を含め、最大400時間単位にわたって安定かつ正確なダイナミクスを達成。
• マクスウェルおよび硬い球体相互作用のための明示的スペクトル重みが導出され、数値的にも検証され、スキームの適用可能性が向上。
• 数値結果から、フォッカー＝プランク＝ランダウ作用素の退化したスペクトル的挙動の解析的証明は未解決の問題のままであるが、クーロン相互作用におけるその鋭さは数値的に確認されている。