일관된 다중완화시간 격자볼츠만법을 이용한 부피 평균 나비에 스토크스 방정식 구현

초록

본 논문은 부피 평균 나비에-스토크스 방정식(VANSE)을 정확히 복원하는 다중완화시간(MRT) 격자볼츠만(LB) 모델을 제안한다. 임시 상태 방정식을 도입해 공극률과 밀도를 분리하고, 모멘트 공간에 페널티 소스항을 추가해 비갤러킨 불변성을 보장한다. Chapman‑Enskog 분석과 다양한 수치 검증을 통해 2차 정확도와 큰 공극률 구배에서도 안정적인 계산이 가능함을 확인하였다.

상세 분석

이 연구는 기존의 밀도 기반 BGK(LB) 모델이 공극률이 급격히 변하는 다상 흐름에서 비물리적인 압력 진동과 스퓨리어스 속도를 야기한다는 근본적인 문제를 정확히 짚어낸다. 저자들은 먼저 압력‑밀도 결합을 풀어내기 위해 임시 상태 방정식 (p=\kappa\rho) 을 도입하고, 이를 통해 공극률 (\phi) 와 밀도 (\rho) 를 독립적으로 취급한다. 이 과정에서 3차 차수의 헤르미트 전개를 이용한 새로운 평형 분포 (f_i^{eq}) 를 정의했으며, 이는 0차·1차 모멘트는 기존 모델과 동일하게 유지하면서 2차 모멘트의 압력 항을 올바르게 수정한다. 그러나 3차 모멘트에서 발생하는 비갤러킨 불변성 위반을 해결하기 위해, 다중완화시간(MRT) 충돌 연산자를 채택하고, 모멘트 공간에 페널티 소스 (\mathbf{C}) 를 삽입한다. 이 소스항은 비대칭 점성 텐서 오류를 정량적으로 보정하여, 연산자 고유값 (s_e, s_\epsilon, s_q) 을 적절히 선택함으로써 수치 안정성을 크게 향상시킨다. Chapman‑Enskog 전개를 통해, 연속 방정식과 모멘텀 방정식이 정확히 VANSE 형태로 복원됨을 증명했으며, 특히 압력 구배 항에 대한 추가 강제항 (\mathbf{F}_p = -\nabla(\kappa\rho) + \nabla(p)) 이 공극률 구배와 무관하게 작용함을 확인했다.

수치 실험에서는 (1) 제조해법을 이용한 정밀 수렴 테스트, (2) 급격한 공극률 변화를 갖는 다공성 채널 흐름에서의 스퓨리어스 속도 억제, (3) 시간·공간적으로 변동하는 공극률 필드에 대한 2차 정확도 검증을 수행했다. 모든 경우에서 MRT‑LB 모델은 기존 SR‑TB 기반 압력 모델보다 넓은 레이놀즈 수와 점성 범위에서 수렴성을 유지했으며, 특히 (\phi<0.5) 영역에서도 안정적으로 계산이 가능했다.

결과적으로, 본 논문은 부피 평균 Navier‑Stokes 방정식에 대한 일관된 LB 스키마를 제공함으로써, 다상 흐름, 다공성 매질, 입자‑유체 결합 시뮬레이션 등 다양한 공학·과학 분야에서 기존 방법이 갖는 한계를 극복할 수 있는 강력한 도구를 제시한다.