[論文レビュー] On the numerical solution of a variable-coefficient Burgers equation arising in granular segregation

Ivan C. Christov|arXiv (Cornell University)|Jun 30, 2017

Granular flow and fluidized beds参考文献 2被引用数 3

ひとこと要約

本稿では、乾燥二成分混合物における粒子分離をモデル化する可変係数バーゲルス方程式を解くために、2次精度で強く陰的なクランク・ニコルソン数値スキームを提示する。この手法は粒子フラックスの非線形境界条件を扱うことができ、正確なキント解との比較により2次収束性と良好な保存性を示し、粒子動力学シミュレーションから得られる運動的応力プロファイルに基づく定数係数および可変係数ケースへの応用も行っている。

ABSTRACT

We study a variable-coefficient Burgers equation arising in the modelling of segregation of dry bidisperse granular mixtures. The equation is subject to nonlinear boundary conditions for the particle flux. We construct a strongly implicit Crank--Nicolson type of numerical scheme for the latter equation. The scheme is benchmarked against a standard exact solution of kink type, showing second-order of accuracy and good discrete conservation properties. Two segregation problems considered in the literature are then solved and discussed. The first is the case of a linear kinetic stress profile, which renders the governing equation of constant-coefficient type, while the second is the case of a variable kinetic stress profile based on an expression fit to particle dynamics simulation data.

研究の動機と目的

• 二成分混合物における粒状分離に由来する可変係数バーゲルス方程式を正確に解くための堅牢な数値スキームの開発。
• 支配方程式における粒子フラックスの非線形境界条件の取り扱い。
• 数値解法における2次精度および離散的保存性の確保。
• 正確なキント型解との比較によるスキームの妥当性検証。
• 2つの物理的状況への適用：粒子動力学シミュレーションから得られる運動的応力プロファイルに基づく線形および可変係数プロファイル。

提案手法

• 可変係数バーゲルス方程式に対して、強い陰性のクランク・ニコルソン型有限差分スキームを構築する。
• 粒子フラックスにおける非線形境界条件を処理できるようにスキームを定式化する。
• 正確なキント型解析解を用いた妥当性検証により、2次収束性が確認された。
• 物理的整合性に不可欠な離散的保存性を保持する。
• 2つの物理的状況をシミュレートした：線形運動的応力プロファイル（定数係数ケース）と、粒子動力学シミュレーションから導出された可変プロファイル。
• 数値結果を既存の文献と比較し、精度および物理的妥当性を評価した。

実験結果

リサーチクエスチョン

• RQ1 粒状分離における非線形フラックス境界条件を有する可変係数バーゲルス方程式を正確に解くための数値スキームは、どのように設計できるか？
• RQ2 提案された数値スキームの収束次数および保存性の挙動は何か？
• RQ3 線形運動的応力プロファイルを仮定した定数係数ケースに適用した場合、モデルはどのように動作するか？
• RQ4 粒子動力学シミュレーションから導出された可変運動的応力プロファイルを仮定した場合、モデルはどのように振る舞うか？
• RQ5 可変係数方程式を数値的に解くことにより、分離ダイナミクスに関するどのような知見が得られるか？

主な発見

• 提案されたクランク・ニコルソンスキームは、正確なキント解との収束テストにより、空間および時間方向で2次精度を達成していることが確認された。
• スキームは良好な離散的保存性を示し、粒子フラックスのモデリングに不可欠な物理的不変量を保持している。
• 線形運動的応力ケースでは、既知の定数係数バーゲルス方程式の挙動と一致する。
• 粒子動力学シミュレーションデータへのフィットに基づく可変係数ケースでは、定数係数モデルでは捉えきれない特徴的な分離ダイナミクスが明らかになった。
• 数値結果は先行研究と整合的であり、モデルの物理的妥当性およびスキームの信頼性が裏付けられた。