[論文レビュー] Purified Two-Relaxation-Time Lattice Boltzmann Method: Removing Ghost Modes from TRT for Enhanced Stability

Yuan Yu, Yaolong Yu+4 more|arXiv (Cornell University)|Feb 6, 2026

Lattice Boltzmann Simulation Studies被引用数 0

ひとこと要約

この論文は Hermite ベースの正則化が TRT-RLB において数学的に TRT 演算子からゴーストモードを除去することと等価であり、Purified TRT (P-TRT) を実現して TRT-RLB の安定性と同等を、計算コストをはるかに低く抑えつつ再現できることを示し、D2Q9 および D3Q19 格子で実証する。

ABSTRACT

The two-relaxation-time (TRT) lattice Boltzmann model is widely adopted for its simplicity and tunable boundary accuracy. However, its collision operator relaxes the full symmetric non-equilibrium component, implicitly retaining non-hydrodynamic ghost modes that degrade stability at high Reynolds numbers. In this work, we establish a rigorous connection between ghost-mode filtering and regularization within the TRT framework. By decomposing the discrete velocity space into hydrodynamic and non-hydrodynamic subspaces, we prove that the TRT-regularized lattice Boltzmann (TRT-RLB) model is mathematically equivalent to the standard TRT model with ghost modes explicitly removed. This equivalence holds exactly for D2Q9 and D3Q19 lattices, where the symmetric and antisymmetric subspaces are completely spanned by the physically relevant Hermite modes and identifiable ghost modes. Based on this finding, we propose the Purified TRT (P-TRT) model, which achieves regularization-level stability through simple algebraic ghost-mode subtraction rather than expensive tensor projections. For D2Q9, the non-equilibrium collision cost is reduced from 180 to 52 floating-point operations per node, a 71% reduction. Linear stability analysis in moment space further reveals that the P-TRT operator annihilates the ghost eigenvalue, proving its spectral radius is bounded above by that of standard TRT and that stability is governed exclusively by hydrodynamic modes. Benchmarks including the double shear layer at Re up to 10^7, Taylor--Green vortex decay, force-driven Poiseuille flow, and creeping flow past a square cylinder confirm that P-TRT preserves the stability, second-order accuracy, and zero-slip boundary properties of TRT-RLB while retaining the simplicity of the TRT family.

研究の動機と目的

• Reynolds 数が高い場合の TRT ライトボルツマン法の安定性制限を動機づけ、対処する。
• ゴーストモードのフィルタリングを TRT フレームワーク内の正則化と結びつける。
• Purified TRT (P-TRT) を導入し、代数的なゴーストモードの減算だけで TRT-RLB の性能に匹敵させる。
• ベンチマーク流れにおいて計算効率の向上と精度・境界条件の維持を示す。

提案手法

• TRT-RLB とその Hermite ベースの非平衡モーメント正則化をレビューする。
• 速度空間を対称性パリティによって対称部分空間（流体力学的+ ゴースト）と反対称部分空間に直交射影で分解する。
• D2Q9 および D3Q19 において二次 Hermite 射影とゴースト減算の完全な同等性を証明する。
• 計算コストの高いテンソル射影を避け、代数的にゴーストモードを除去して P-TRT を提案する。
• ゴーストモードが P-TRT のスペクトル半径を標準の TRT を超えて不安定にしないことを示す線形安定性分析を提供する。
• FLOPs を削減しつつ安定性を保ち、精度を検証する benchmark テストを提示する。

実験結果

リサーチクエスチョン

• RQ1 TRT でのゴーストモードを完全な Hermite 射影なしで明示的に除去して物理を保てるか。
• RQ2 TRT-RLB における Hermite ベース正則化は TRT 内のゴーストモードフィルタリングと異なるのか。
• RQ3 Purified TRT (P-TRT) は TRT-RLB の安定性と精度を低コストで保持できるか。
• RQ4 D2Q9, D3Q19 などの一般的な格子において P-TRT は TRT-RLB と比較してどの程度計算節約になるか。

主な発見

• P-TRT は D2Q9 および D3Q19 に対して TRT-RLB と数学的に等価である。対称部分空間におけるゴーストモードの完全な除去と既存モードを含む反対称部分空間の除去が達成されるため。
• ゴーストモードの消去により P-TRT のスペクトル半径は標準の TRT を上回らないことが保証され、ゴースト固有値はゼロに固定される。
• D2Q9 では非平衡衝突の総コストが元の TRT-RLB の 180 FLOPs から P-TRT の 52 FLOPs に低下し、71%の削減となる。
• P-TRT は安定性、二次精度、ゼロ滑り境界特性を TRT-RLB と同等に維持しつつ、計算オーバーヘッドを大幅に削減する。
• 双重せん断層（Re が最大 1e7 まで）、減衰する Taylor–Green渦、力駆動ポアソイル流、スクエア状円筒を越える流体のクリープ流などのベンチマークが P-TRT の安定性と精度を検証する。