[論文レビュー] Analysis and boundary condition of the lattice Boltzmann BGK model with two velocity components

Xiaoyi He, Qisu Zou|arXiv (Cornell University)|Jul 22, 1995

Lattice Boltzmann Simulation Studies参考文献 21被引用数 15

ひとこと要約

本稿では、2次元格子ボルツマンBGKモデルに2つの速度成分を適用し、パイプ流れとクーレット流れの解析的解を用いて、ナビエ＝ストークス方程式の2次精度の有限差分近似を導出することで、理論的裏付けのある境界条件を提案する。この手法は運動量交換の一貫性を保証し、境界条件が適切に適用されれば、速度プロファイルにおいて機械精度に達する結果を得る。

ABSTRACT

In this paper, we study the two dimensional lattice Boltzmann BGK model (LBGK) by analytically solving a simple flow in a 2~-D channel. The flow is driven by the movement of upper boundary with vertical injection fluid at the porous boundaries. The velocity profile is shown to satisfy a second-order finite-difference form of the simplified incompressible Navier-Stokes equation. With the analysis, different boundary conditions can be studied theoretically. A momentum exchange principle is also revealed at the boundaries. A general boundary condition for any given velocity boundary is proposed based on the analysis.

研究の動機と目的

• 格子ボルツマンBGKモデルにおける流体力学的境界条件に理論的根拠が欠如している問題に対処すること。
• 注入とせん断力によって駆動される、多孔質境界と移動壁を有する2次元チャネル流れにおける速度プロファイルを分析すること。
• 任意の指定された速度境界に対して、2次精度と正しい運動量交換を保証する一般境界条件を導出すること。
• LBGK方程式の解析的解に基づく数学的に厳密な枠組みを提供することで、既存の境界スキームを拡張すること。
• 提案された境界条件を、解析的解との比較と誤差解析を通じて検証すること。

提案手法

• 2次元チャネル内における定常かつ非圧縮性流れに対して、LBGK方程式を解析的に解く。速度成分はxおよびy方向に存在する。
• 定常状態および密度一定の仮定の下で、LBGKモデルからナビエ＝ストークス方程式の2次精度の有限差分形を導出する。
• 境界速度の一貫性を補正するバウンスバック則を適用し、境界における運動量交換原理を導入する。
• 境界ノードにおける既知の密度および速度に基づいて分布関数を調整することで、正しい巨視的速度を強制する一般境界条件を用いる。
• 境界条件によって定まる係数を有する指数関数的解の線形結合を用いて、速度プロファイルの解を求める。
• 指定された境界速度と計算された速度との差に起因する境界誤差を定量化するための補正項（ε₀, εₙ）を導入する。

実験結果

リサーチクエスチョン

• RQ1 2次元チャネル流れにおける2速度成分付き格子ボルツマンBGKモデルは、ナビエ＝ストークス方程式をどのように近似するか？
• RQ2 標準的なバウンスバック境界条件が、LBGKモデルに導入する誤差の性質は何か？
• RQ3 任意の指定された速度境界に対して2次精度を保証する一般境界条件を導出可能か？
• RQ4 2速度成分付きのLBGKモデルにおいて、境界における運動量交換原理はどのように現れるか？
• RQ5 LBGK法における境界条件を体系的に評価・改善するための解析的枠組みとして、どのようなものがあるか？

主な発見

• 2次元LBGKモデルにおける速度プロファイルは、非圧縮性ナビエ＝ストークス方程式の2次精度有限差分近似を満たす。
• 境界速度が正しく強制されれば、平面パイプ流れにおいて提案された境界条件は機械精度に達する。
• 速度プロファイルの解は、解析的解の2次近似であり、誤差はレイノルズ数および格子間隔によって上限が定まる。
• 補正項ε₀およびεₙは、指定された境界速度と計算された速度との差に起因する境界誤差を定量化するものであり、境界条件が満たされた場合には消える。
• 一般境界条件は3次元15速度モデルへ拡張可能であり、法線速度がゼロの場合には既存のスキーム（例：Maierら）に還元される。
• 解析により、境界における運動量交換原理が明らかとなり、巨視的速度と分布関数との一貫性が保証される。