[論文レビュー] Lattice-Boltzmann Simulations of Fluid Flows in MEMS

Xiaobo Nie, Gary D. Doolen+1 more|arXiv (Cornell University)|Jun 11, 1998

Lattice Boltzmann Simulation Studies被引用数 23

ひとこと要約

本稿では、格子ボルツマン法（LBM）がマイクロエレクトロメカニカルシステム（MEMS）におけるマイクロスケールの流れを正確にシミュレートできることを示しており、マイクロチャネルおよびマイクロキャビティ内の速度スリップと非線形圧力勾配を捉えている。LBMはクヌーゼン数に依存するスリップ速度（$V_s \propto K_n^2$ のスケーリング）と、線形性からの非単調な圧力偏差をモデル化でき、解析的モデルおよび実験結果と一致することを検証した。

ABSTRACT

The lattice Boltzmann model is a simplified kinetic method based on the particle distribution function. We use this method to simulate problems in MEMS, in which the velocity slip near the wall plays an important role. It is demonstrated that the lattice Boltzmann method can capture the fundamental behavior in micro-channel flow, including velocity slip and nonlinear pressure drop along the channel. The Knudsen number dependence of the position of the vortex center and the pressure contour in micro-cavity flows is also demonstrated.

研究の動機と目的

• 連続体仮説が成立しないマイクロエレクトロメカニカルシステム（MEMS）における希薄ガス流れのシミュレーションに、格子ボルツマン法（LBM）の適用可能性を調査すること。
• MEMSデバイスにおける質量移動および熱移動に顕著な影響を与える、マイクロチャネルおよびマイクロキャビティ内の壁面における速度スリップをモデル化すること。
• マイクロスケールの幾何形状における流れ構造（たとえば渦の位置や圧力等高線）がクヌーゼン数（$K_n$）にどのように依存するかを検討すること。
• マイクロチャネルにおける質量流量および圧力勾配について、LBMの結果を解析的解および実験データと照合して検証すること。
• 非線形圧力プロファイルおよびスリップ速度のスケーリングを通じて、連続体から希薄流への挙動の遷移をシミュレーションで探ること。

提案手法

• 規則的なグリッド上で離散ボルツマン方程式を解くために、2次元9速度（D2Q9）の格子ボルツマンモデルとBGK衝突近似を用いる。
• 圧縮性に類似したマイクロフローにおける可変粘性を考慮するため、密度に依存する緩和時間 $\tau' = \frac{1}{2} + \frac{1}{\rho}(\tau - \frac{1}{2})$ を導入する。
• 固体壁ではバウンスバック境界条件を、流入・流出端では圧力境界条件を適用してマイクロチャネル流れをシミュレートする。
• チャップマン・エンスコグ解析を用いて、低$K_n$極限における連続体流体力学とLBMを結びつける、スリップ境界条件付きのナビエ＝ストークス方程式を導出する。
• 最小二乗フィッティングを用いて速度プロファイルからスリップ速度 $V_s$ を抽出し、解析的モデルと比較する。
• 上部壁が動くマイクロキャビティの流れをシミュレートし、クヌーゼン数の変化に伴う渦の形成および圧力構造を研究する。

実験結果

リサーチクエスチョン

• RQ1 格子ボルツマン法はマイクロチャネル流れにおける速度スリップをどの程度正確に捉えることができ、スリップ速度はクヌーゼン数にどのように依存するか？
• RQ2 マイクロチャネルにおける圧力勾配の性質は何か？高クヌーゼン数において線形挙動から逸脱するか？
• RQ3 マイクロキャビティ流れにおいてクヌーゼン数は、渦中心の位置および圧力等高線構造にどのように影響を与えるか？
• RQ4 LBMは希薄流れで観測された線形からの非単調な圧力偏差を再現できるか？
• RQ5 マイクロチャネルにおける質量流量は、LBMシミュレーション、解析的モデル、実験測定と比較してどの程度一致するか？

主な発見

• 格子ボルツマン法はマイクロチャネルにおける速度スリップを正確に捉えており、スリップ速度は $V_s = 8.7K_n^2$ とスケーリングされ、理論的予測と一致する。
• $K_n \geq 0.2$ の場合、圧力勾配は線形挙動から負に逸脱するという非単調効果を示し、これは以前の実験では観測されておらず、LBMモデルによって予測されたものである。
• マイクロチャネルにおける正規化された質量流量は、$\eta = 2$ の場合、解析的モデル $M_f = 1 + 24.1K_n^2$ とよく一致しており、LBMの正確性が裏付けられる。
• マイクロキャビティ流れでは、クヌーゼン数が増加するにつれて渦中心が上方に移動し、壁面スリップに起因する運動量移動の減少により質量流量が減少する。
• 圧力等高線は連続体極限では円形であるのに対し、高クヌーゼン数ではほぼ直線に近づき、連続体挙動からの顕著な逸脱を示している。
• $K_n = 0.165$ において、LBMシミュレーションは質量流量について実験データおよび解析的モデルと非常に良好に一致しており、特に圧力比 $\eta \geq 1.8$ の場合に顕著である。