[論文レビュー] Implementation of on-site velocity boundary conditions for D3Q19 lattice Boltzmann

Martin Hecht, Jens Harting|TU/e Research Portal (Eindhoven University of Technology)|Nov 27, 2008

Lattice Boltzmann Simulation Studies被引用数 10

ひとこと要約

本稿では、三次元D3Q19格子ボルツマン法に局所的で現地適用可能な速度境界条件を提示し、格子の向きとは無関係に任意のフラックス方向を可能にする。この手法は2次精度を保証し、数値的スリップを排除し、速度をゼロに設定することで自然に滑らかでない条件を含む。複数のテストケースによる検証により、複雑な幾何形状においても高い精度と頑健性が確認された。

ABSTRACT

On-site boundary conditions are often desired for lattice Boltzmann simulations of fluid flow in complex geometries such as porous media or microfluidic devices. The possibility to specify the exact position of the boundary, independent of other simulation parameters, simplifies the analysis of the system. For practical applications it should allow to freely specify the direction of the flux, and it should be straight forward to implement in three dimensions. Furthermore, especially for parallelized solvers it is of great advantage if the boundary condition can be applied locally, involving only information available on the current lattice site. We meet this need by describing in detail how to transfer the approach suggested by Zou and He to a D3Q19 lattice. The boundary condition acts locally, is independent of the details of the relaxation process during collision and contains no artificial slip. In particular, the case of an on-site no-slip boundary condition is naturally included. We test the boundary condition in several setups and confirm that it is capable to accurately model the velocity field up to second order and does not contain any numerical slip.

研究の動機と目的

• 3次元D3Q19格子ボルツマン法に一般化され、局所的に適用可能な速度境界条件を開発すること。これにより、任意の流入方向を許容する。
• 数値的スリップを排除し、速度場再構築における2次精度を確保すること。
• BGK緩和時間に依存しない境界条件を提供することで、多孔質媒体やマイクロ流体デバイスのシミュレーションにおける正確な壁位置決めを可能にすること。
• ZouとHeの手法をD2Q9からD3Q19に一般化し、非垂直フラックスに対応する横方向運動量補正を含む。
• 壁が直交座標平面と一致しない複雑な幾何形状（例：マイクロミキサー、多孔質材料など）においても正確なシミュレーションを可能にすること。

提案手法

• D2Q9用に開発されたZouとHe（1997）の手法をD3Q19格子に適応し、任意の流入方向を想定したフラックス境界条件を導出する。
• 境界ノードにおける分布関数を、既知の密度と速度成分を用いて局所的に再構築する局所的衝突ステップを採用する。
• 非平衡成分に対してバンプバック則を適用し、運動量バランスを維持するために、項 $ N^{x}_{y}, N^{x}_{z}, N^{y}_{x}, N^{y}_{z} $ を用いた横方向運動量補正を追加する。
• 速度成分と隣接する分布関数に基づいて、入口・出口境界における密度および分布関数の式を導出する。
• 速度成分と局所的な分布値に依存する補正を組み込むことで、平衡分布関数と整合性を保つ。
• 入口・出口境界の両方に対応する対称的定式化を採用し、各面（例：x=0, x=max, y=0, y=max）に対して異なる方程式を定義することで、局所的適用性を確保する。

実験結果

リサーチクエスチョン

• RQ1 ZouとHeの境界条件アプローチを、任意のフラックス方向を想定した3次元D3Q19格子ボルツマンシミュレーションに一般化できるか？
• RQ2 提案手法は、数値的スリップが存在しない状態で、速度場再構築において2次精度を達成するか？
• RQ3 現在の格子節点からの情報のみを用いて局所的に適用可能であり、大規模シミュレーションにおける効率的な並列処理を可能にするか？
• RQ4 壁が格子軸と一致しない複雑な幾何形状において、この手法はどのように性能を発揮するか？
• RQ5 速度をゼロに設定することで、滑らかでない境界条件が自然に含まれるか？また、BGK緩和時間に依存しないか？

主な発見

• 数値収束テストにより、提案された境界条件が速度場再構築において2次精度を達成していることが確認された。
• 分布関数再構築における適切な横方向運動量補正のおかげで、非垂直フラックス方向に対しても数値的スリップが存在しない。
• 現在の格子節点からの情報のみを必要とするため、局所的適用が可能であり、大規模シミュレーションにおける効率的な並列処理を支援する。
• 速度ベクトルをゼロに設定することで、滑らかでない境界条件が自然に含まれ、壁の位置がBGK緩和時間に依存しない。
• チャネルの任意の方向への流れや多孔質媒体を含む複数の設定での検証により、多様な幾何形状において頑健性と精度が確認された。
• 境界における分布関数の導出式（例：式51–54, 65–68）は、格子ボルツマン平衡関数および運動量保存則と整合性を持つ。