[论文解读] 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

一句话总结

本文证明了格子Boltzmann方法（LBM）能够准确模拟微机电系统（MEMS）中的微尺度流体流动，捕捉到微通道和微腔体中的速度滑移及非线性压力降。结果表明，LBM可模拟与克努森数（Knudsen number）相关的滑移速度（其标度关系为 $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.

研究动机与目标

• 研究格子Boltzmann方法（LBM）在微机电系统（MEMS）中模拟稀薄气体流动的适用性，此时连续介质假设不再成立。
• 在微通道和微腔体中模拟壁面滑移速度，该效应显著影响MEMS器件中的质量与热量传递。
• 研究流动结构（如涡心位置和压力等值线）在微尺度几何结构中随克努森数（$K_n$）的变化规律。
• 通过微通道中的质量流量和压力降数据，将LBM模拟结果与解析解和实验数据进行验证。
• 通过模拟非线性压力分布和滑移速度标度关系，探索从连续流到稀薄流行为的过渡过程。

提出的方法

• 采用二维九速度（D2Q9）格子Boltzmann模型，结合BGK碰撞近似，在规则网格上求解离散Boltzmann方程。
• 引入依赖于密度的松弛时间 $\tau' = \frac{1}{2} + \frac{1}{\rho}(\tau - \frac{1}{2})$，以考虑可压缩类微流中粘度的变化。
• 在固体壁面采用反弹边界条件，在进出口采用压力边界条件，以模拟微通道中的流动。
• 利用Chapman-Enskog分析推导出带有滑移边界条件的Navier-Stokes方程，将LBM与低$K_n$极限下的连续介质流体力学联系起来。
• 采用最小二乘拟合方法从速度分布中提取滑移速度 $V_s$，并与解析模型进行比较。
• 对带有移动顶壁的微腔体流动进行模拟，研究不同克努森数下涡旋形成与压力结构的变化。

实验结果

研究问题

• RQ1 格子Boltzmann方法在多大程度上能准确捕捉微通道流动中的速度滑移？滑移速度如何随克努森数变化？
• RQ2 微通道中的压力降具有何种特性？在较高克努森数下是否偏离线性行为？
• RQ3 克努森数如何影响微腔体流动中涡心位置和压力等值线结构？
• RQ4 LBM能否再现稀薄流中观测到的压力偏离线性的非单调行为？
• RQ5 微通道中的质量流量在LBM模拟、解析模型与实验测量之间有何比较？

主要发现

• 格子Boltzmann方法成功捕捉了微通道中的速度滑移，滑移速度标度关系为 $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$）条件下。