[論文レビュー] Beyond pressureless gas dynamics: Quadrature-based velocity moment models

Christophe Chalons, Damien Kah+1 more|arXiv (Cornell University)|Nov 12, 2010

Gas Dynamics and Kinetic Theory参考文献 18被引用数 42

ひとこと要約

本稿は、無限Knudsen数極限における運動方程式の四次モーメント系を、二点則による積分近似を用いて拡張し、滑らかでない解および特異解を捉えることができる四モーメント系を提案する。エントロピー測度解の厳密な枠組みを構築し、キネティックスキームを用いた検証を通じて収束性と粒子軌道の重なりや一般化された$δ$-ショックの適切な取り扱いを示している。

ABSTRACT

Following the seminal work of F. Bouchut on zero pressure gas dynamics which has been extensively used for gas particle-flows, the present contribution investigates quadrature-based velocity moments models for kinetic equations in the framework of the infinite Knudsen number limit, that is, for dilute clouds of small particles where the collision or coalescence probability asymptotically approaches zero. Such models define a hierarchy based on the number of moments and associated quadrature nodes, the first level of which leads to pressureless gas dynamics. We focus in particular on the four moment model where the flux closure is provided by a two-node quadrature in the velocity phase space and provide the right framework for studying both smooth and singular solutions. The link with both the kinetic underlying equation as well as with zero pressure gas dynamics is provided and we define the notion of measure solutions as well as the mathematical structure of the resulting system of four PDEs. We exhibit a family of entropies and entropy fluxes and define the notion of entropic solution. We study the Riemann problem and provide a series of entropic solutions in particular cases. This leads to a rigorous link with the possibility of the system of macroscopic PDEs to allow particle trajectory crossing (PTC) in the framework of smooth solutions. Generalized $\\delta$-choc solutions resulting from Riemann problem are also investigated. Finally, using a kinetic scheme proposed in the literature without mathematical background in several areas, we validate such a numerical approach in the framework of both smooth and singular solutions.

研究の動機と目的

• 速度位相空間における二点則による積分近似を用いた四モーメント系を導入することで、圧力なしガス力学を超えて、希薄な粒子流れにおける複雑な位相空間ダイナミクスを捉える四次モーメント系の拡張を図ること。
• キネティック方程式に由来する弱超臨界保存則系の文脈において、測度解とエントロピー解の数学的枠組みを確立すること。
• 速度が一致する場合を含め、モーメント空間の境界における系の挙動を調査し、数値的実現可能性を保証すること。
• キネティックスキームを用いて数値的手法を検証し、滑らかでない解および特異解の両方を的確に捉える収束性と精度を示すこと。
• 粒子の軌道が重なることや複雑な位相空間ダイナミクスを有する粒子を含む流れのモデル化のための理論的・数値的基盤を提供すること。

提案手法

• 速度分布関数の二点則近似に基づく四モーメント系を定式化し、分布関数をデルタ関数の和として表現する。
• 系を弱超臨界保存則系として定義し、実現可能性条件およびモーメント空間の境界を含む数学的構造を特徴付ける。
• 系のダイナミクスと整合する一連のエントロピーとエントロピー流束を構成することで、エントロピー測度解の概念を導入する。
• 積分ノードと重みに基づくキネティックスキームを適用し、実現可能性と特異性の適切な取り扱いを保証して系を数値的に解く。
• 構造化グリッド上の$L_1$ノルムを用いた誤差ノルムを用いて、リーマン問題の解と収束性の検証を実施する。
• $q/e^{3/2}$などの量のモニタリングを通じて、モーメント空間の内部から境界への遷移を分析し、錐体保存を確認する。

実験結果

リサーチクエスチョン

• RQ1 積分近似モーメントモデルは、圧力なしガス力学を超えて、希薄な粒子流れにおける複雑な位相空間ダイナミクスをどのように拡張できるか？
• RQ2 四モーメント系の数学的構造は何か？また、滑らかでない解および一般化された$δ$-ショックを含む特異解をどのように扱えるか？
• RQ3 この系におけるエントロピー測度解はどのように定義され、特徴づけられるか？また、物理的整合性を保証する役割は何か？
• RQ4 キネティックスキームは実現可能性を保持し、モーメント空間の境界における粒子軌道の重なりや遷移を的確にシミュレートできるか？
• RQ5 数値スキームの収束性はいかほどか？また、リーマン問題の解析解を捉える性能はどの程度か？

主な発見

• 四モーメント積分近似モデルは、滑らかな解と一般化された$δ$-ショックの両方を的確に捉えており、粒子軌道の重なりを含む複雑なダイナミクスのモデル化が可能であることを示している。
• キネティックスキームは、400セルから3200セルへのグリッドの細分化に伴う誤差低減を確認したところ、$L_1$ノルムにおける収束率が約0.5であることが判明した。
• スキームは実現可能性錐体を保持し、モーメント空間の境界でも数値的安定性を維持しており、$q/e^{3/2}$の挙動と不自然な振動の不在によって裏付けられている。
• 系はエントロピーとエントロピー流束の族を支持しており、これにより、元のキネティック方程式と整合するエントロピー測度解を定義できる。
• 系は単一速度極限において圧力なしガス力学に帰着し、既存の枠組みと整合性があることが確認された。
• 数値的検証により、スキームは特異解や軌道の重なりを含むリーマン問題の解析解を的確に再現していることが示された。