[論文レビュー] Numerical Simulations of Particulate Suspensions via a Discretized Boltzmann Equation Part II. Numerical Results

Anthony J. C. Ladd|arXiv (Cornell University)|Jul 1, 1993

Lattice Boltzmann Simulation Studies被引用数 9

ひとこと要約

本稿では、粒子数に線形にスケーリングする離散化ボルツマン方程式法を用いて、コロイド系における流体力学的相互作用の効率的モデリングを可能にする、粒子系の数値シミュレーションを提示する。この手法は、ブラウン運動を伴う・なしの両方のクリーピングフロー、および有限レイノルズ数の流れを正確にシミュレートでき、1024粒子まで安定したダイナミクスを達成し、輸送係数と短時間ダイナミクスが検証済みである。

ABSTRACT

A new and very general technique for simulating solid-fluid suspensions has been described in a previous paper (Part I); the most important feature of the new method is that the computational cost scales with the number of particles. In this paper (Part II), extensive numerical tests of the method are described; for creeping flows, both with and without Brownian motion, and at finite Reynolds numbers. Hydrodynamic interactions, transport coefficients, and the short-time dynamics of random dispersions of up to 1024 colloidal particles have been simulated.

研究の動機と目的

• 粒子数に比例する線形な計算スケーリングを実現する、新規な数値手法を、粒子系のシミュレーションに検証すること。
• 熱揺らぎを伴う・なしの両方のクリーピングフロー条件下におけるコロイド系の流体力学的相互作用を調査すること。
• この手法を有限レイノルズ数の流れに拡張し、複雑な系のダイナミクスを的確に捉えられるかを評価すること。
• 最大1024粒子の不規則な分散系における輸送係数を計算し、短時間ダイナミクスを分析すること。
• ソフトマターおよび流体力学研究分野における固体-流体系のシミュレーションに、計算的に効率的で汎用性の高いフレームワークを提供すること。

提案手法

• 流体-粒子相互作用をモデル化するため、離散化された速度空間を用いた格子ボルツマン法を採用する。
• 規則的な格子上でバウンスバック則を用いて、粒子表面でのスリップなし境界条件を実装する。
• 格子ボルツマン方程式に確率的揺らぎを組み込み、ブラウン運動をランダム力項としてモデル化する。
• 衝突とストリーミングステップを伴うD3Q19格子上で離散ボルツマン方程式を解き、流体の流れと粒子の運動をシミュレートする。
• 個々の粒子を追跡し、流体場を介してその流体力学的相互作用を計算する、粒子解像度のアプローチを採用する。
• 質量および運動量の保存を保証するために、時間積分アルゴリズムを用いて系を離散時間ステップごとに進める。

実験結果

リサーチクエスチョン

• RQ1 提案手法は、低レイノルズ数におけるコロイド系の流体力学的相互作用を正確にシミュレートできるか？
• RQ2 この手法は、系内での粒子ダイナミクスに及ぼすブラウン運動の効果をどの程度的確に捉えられるか？
• RQ3 有限レイノルズ数において、この手法はどの程度安定的かつ正確に保たれるか？
• RQ4 シミュレーションで予測される輸送係数（例：自己拡散率、粘度）は何か？理論的予測と比較するとどうなるか？
• RQ5 計算コストは粒子数にどの程度依存するか？また、大規模系（例：1024粒子）を効率的に扱えるか？

主な発見

• 計算コストが粒子数に線形にスケーリングすることを達成し、最大1024個のコロイド粒子のシミュレーションが可能となった。
• クリーピングフロー下では流体力学的相互作用が正確に捉えられ、予測された自己拡散率が理論的予測と一致した。
• 確率的力項の導入により、ブラウン運動の効果が的確に再現され、熱揺らぎモデルの妥当性が検証された。
• 有限レイノルズ数のシミュレーションでは、安定的かつ物理的に整合性のある流れ場と粒子軌道が得られた。
• 不規則な分散系の短時間ダイナミクス、特に粒子の拡散および相関関数が、理論的モデルと良好に一致した。
• 低〜中程度のレイノルズ数の範囲にわたり、本手法は強健性と正確性を示した。