[논문 리뷰] Lattice-Boltzmann Simulations of Fluid Flows in MEMS

Xiaobo Nie, Gary D. Doolen외 1명|arXiv (Cornell University)|1998. 06. 11.

Lattice Boltzmann Simulation Studies인용 수 23

한 줄 요약

이 논문은 연속체 가정이 실패하는 마이크로-전자기기(MEMS)에서 미세 체적 유동을 정확하게 시뮬레이션할 수 있음을 보여주며, 미세 채널 및 미세 캐비티에서의 속도 슬립과 비선형 압력 강하를 포착한다. LBM는 Knudsen 수에 따라 변하는 슬립 속도(스케일링: $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 장치의 질량 및 열전달에 미치는 영향을 분석하기 위해.
• 마이크로 스케일 기하구조에서 유동 구조(예: 소용돌이 중심 위치 및 압력 등고선)가 Knudsen 수($K_n$)에 따라 어떻게 변화하는지 조사하기 위해.
• 마이크로 채널에서 질량 유량 및 압력 강하에 대해 LBM 결과를 분석 모델 및 실험 데이터와 비교하여 검증하기 위해.
• 비선형 압력 프로파일과 슬립 속도 스케일링을 통한 연속체 유동에서 희박 유동으로의 전이 행동을 시뮬레이션을 통해 탐구하기 위해.

제안 방법

• 정규 격자에 정의된 이산 볼츠만 방정식을 해결하기 위해, BGK 충돌 근사법을 사용하는 2차원 9속도(D2Q9) 격자-볼츠만 모델을 사용한다.
• 압축성 유사 미세 유동에서 점도가 변화하는 것을 고려하기 위해 밀도에 따라 변하는 리듬 시간 $\tau' = \frac{1}{2} + \frac{1}{\rho}(\tau - \frac{1}{2})$ 를 도입한다.
• 고체 벽면에서는 반사 경계 조건을 적용하고, 입구 및 출구에는 압력 경계 조건을 적용하여 마이크로 채널 유동을 시뮬레이션한다.
• Chapman-Enskog 분석을 통해 나이퀴스트-스토크스 방정식과 슬립 경계 조건을 유도하여, LBM이 저-$K_n$ 한계에서 연속체 유체역학과 연결됨을 입증한다.
• 최소 제곱법을 사용하여 속도 프로파일에서 슬립 속도 $V_s$ 를 추출하고 분석 모델과 비교한다.
• 움직이는 상부 벽면을 가진 마이크로 캐비티 유동을 시뮬레이션하여 다양한 Knudsen 수에서 소용돌이 형성 및 압력 구조를 연구한다.

실험 결과

연구 질문

• RQ1 격자-볼츠만 방법은 마이크로 채널 유동에서 속도 슬립을 얼마나 정확하게 포착할 수 있으며, 슬립 속도는 Knudsen 수에 따라 어떻게 스케일링되는가?
• RQ2 마이크로 채널에서의 압력 강하는 어떤 성격을 가지며, 높은 Knudsen 수에서 선형성에서 벗어나는가?
• RQ3 Knudsen 수는 마이크로 캐비티 유동에서 소용돌이 중심 위치와 압력 등고선 구조에 어떤 영향을 미치는가?
• 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의 정확성을 검증한다.
• 마이크로 캐비티 유동에서는 Knudsen 수가 증가함에 따라 소용돌이 중심이 상향 이동하고 질량 유량이 감소하는데, 이는 벽면 슬립으로 인한 운동량 전달 감소 때문이다.
• 압력 등고선은 연속체 근사에서 원형에서 고 Knudsen 수에서는 거의 직선으로 변형되어 연속체 행동에서의 강한 이탈을 나타낸다.
• LBM 시뮬레이션은 $K_n = 0.165$ 에서 실험 데이터 및 분석 모델과 우수한 일치를 보이며, 특히 $\eta \geq 1.8$ 인 고압력 비율에서 뛰어난 일치를 보였다.