[論文レビュー] Comment on the paper Li-Shi Luo, Wei Liao, Xingwang Chen, Yan Peng and Wei Zhang, Numerics of the lattice Boltzmann method: Effects of collision models on the lattice Boltzmann simulations, Physical Review E 83, 056710 (2011)

I. V. Karlin, Sauro Succi|arXiv (Cornell University)|Jul 15, 2011

Lattice Boltzmann Simulation StudiesEngineering被引用数 16

ひとこと要約

このコメントは、2011年のPhysical Review E論文が、エントロピー則に基づく格子ボルツマン方程式（ELBE）が標準的なLBM（LBGK）モデルに劣っていると主張しているのを反論するものである。著者らは、ELBEとしてラベル付けされたシミュレーションが実際には、一定の緩和時間を持つ標準的なLBM（LBGK）スキームであったことを示している。これは、エントロピーに基づく自己適応的緩和パラメータを用いる真のELBEとは異なる。その結果、比較は誤りであり、循環的であるため、ELBEの安定性と性能に関する論文の結論は根拠を失う。

ABSTRACT

Critical comments on the entropic lattice Boltzmann equation (ELBE), by Li-Shi Luo, Wei Liao, Xingwang Chen, Yan Peng and Wei Zhang, Physical Review E 83, 056710 (2011), are based on simulations which make use of a model that, despite being called ELBE by the authors, is in fact fully equivalent to the standard lattice Bhatnagar-Gross-Krook equation. As a result, the conclusion of Luo et al on ELBE is circular, hence devoid of scientific bearing.

研究の動機と目的

• Luoら（2011）の科学的誤解を是正すること。同論文はELBEがLBMに劣ると主張しているが、これは誤りである。
• ELBEスキームが、Hの定理に従うエントロピーに基づく自己適応的緩和パラメータを用いる点で、LBGKとは本質的に異なることの明確化。
• Luoらが真のELBEをシミュレートしておらず、代わりに一定の緩和時間を使用していたことの証明。このため、比較は無効である。
• ELBEの真の価値は、直接数値シミュレーション（DNS）においてではなく、解像度が不足する状況での小スケール安定化（可変粘性を介して）にあり、それゆえにその有用性が発揮されることの強調。
• ELBEにおける可変粘性が非物理的であるという主張に反論し、このような有効粘性は計算流体力学で一般的であり、高レイノルズ数流れの解析に不可欠であると主張する。

提案手法

• Luoらがシミュレーションで一定の緩和時間τを使用していたことの特定。これは、エントロピー条件によって動的に決定される緩和パラメータαに依存するELBEのコアなメカニズムに反する。
• 正しいELBEの定式化の再確認：f_i(x+v_i,t+1) - f_i(x,t) = αβ(f^eq_i - f_i) として、人口が進化する。ここでαは、H(αf^eq + (1-α)f) = H(f) というエントロピー条件を解くことで決定され、Hの定理が保証される。
• ELBEとLBGKの唯一の差異は、平衡分布関数と緩和パラメータの性質にある。LBGKでは緩和パラメータが固定されているが、ELBEでは自己適応的である。
• 人口が平衡に近い（解像度が高い）状況では、ELBEはα = 2のときLBM（LBGK）に還元され、DNS領域ではLBM（LBGK）と同一となる。
• Luoらが用いたエントロピー平衡の多項式近似は、エントロピー条件（2）が満たされないため、真のELBE挙動を回復できないことの指摘。
• Keatingら、Spasovら、Geerdink & Hoekstra、Yasuda & Satofukaによる独立したELBE実装が、エントロピー条件と自己適応的緩和の正しく適用されていることを確認している。

実験結果

リサーチクエスチョン

• RQ1 エントロピー則に基づく格子ボルツマン方程式（ELBE）で一定の緩和時間τを用いた場合、標準的なLBM（LBGK）モデルと同等の数値的安定性を示すか？
• RQ2 緩和パラメータを固定した場合、ELBEスキームはLBM（LBGK）モデルと本質的に異なるか？
• RQ3 エントロピー条件（2）を満たさないシミュレーションがELBEとしてラベル付けされた場合、科学的に妥当であるか？
• RQ4 ELBEが「最も劣っている」という主張がなぜ科学的に正当化されないのか？
• RQ5 ELBEにおける可変粘性は非物理的であり、計算効率を著しく損なうものなのか？その結果、ELBEそのものが無効とされるのか？

主な発見

• Luoら（2011）がELBEとラベル付けしたシミュレーションは、ELBEではなく、一定の緩和時間τを持つ標準的なLBM（LBGK）スキームである。このため、比較は無効である。
• Luoら自身の論文内で、彼らが可変緩和時間でのELBEのテストを行わなかったと認めている。これは、彼らの実装が真のELBEではなかったことを裏付ける。
• エントロピー条件（2）はELBEのHの定理と小スケール安定化の根幹を成す。この条件が欠落している場合、そのシミュレーションはELBEとは見なせない。
• エントロピー平衡とその多項式近似との差はO(u^4)であり、低マッハ数流れにおける格子ボルツマン法の切り捨て誤差の範囲内にある。
• ELBEの真の強みは、解像度が不足する状況において、必要な場所と必要なときにのみ可変粘性を発動させる自己適応的安定化能力にある。一方、LBGKはこのような状況で失敗する。
• 非線形エントロピー条件を解く計算コストは、高レイノルズ数で必要なグリッド解像度が低下する恩恵により上回られ、実用的応用においてELBEは全体としてより効率的である。