중심 모멘트 격자볼츠만법의 수치 특성 종합 평가

초록

본 논문은 중심 모멘트 기반 격자볼츠만법(Cascaded LBM)의 정확도·수렴성·안정성을 체계적으로 검증한다. Poiseuille 흐름, 감쇠 Taylor‑Green 와류, 뚜껑 구동 캐비티 등 세 가지 표준 벤치마크에 대해 2차 정확도와 미분량(변형률 텐서)까지의 수렴을 확인하고, 기존 MRT·BGK 모델 대비 중간 정도의 연산 비용 증가만으로 현저한 안정성 향상을 입증한다.

상세 분석

이 연구는 중심 모멘트(Lattice Central Moments)를 이용한 다중 이완 시간(MRT) 충돌 연산자를 ‘계단식(cascaded)’ 구조로 구현한 LBM을 정량적으로 평가한다. 먼저, 확산 스케일링(diffusive scaling) 하에서 격자 간격 Δx와 시간 간격 Δt를 동시에 축소했을 때 속도와 압력 구배, 그리고 변형률 텐서 성분이 모두 2차 수렴을 보임을 수치 실험을 통해 입증한다. 이는 원시 모멘트 기반 MRT와 비교했을 때, 중심 모멘트가 Galilean 불변성을 자연스럽게 보존함으로써 비선형 대류항에 대한 수치 오차가 감소함을 의미한다.

Poiseuille 흐름에서는 이론적 포아송 해와 거의 일치하는 속도 프로파일을 얻었으며, 전단 응력(∂u/∂y)까지 2차 정확도를 달성했다. 감쇠 Taylor‑Green 와류에서는 시간에 따라 진동이 감쇠되는 정확한 해를 재현함과 동시에, 와류의 에너지 스펙트럼이 이론적 -2 지수 법칙을 따르는지를 확인했다. 특히, 고레졸루션에서 발생할 수 있는 비물리적 진동을 억제하는 데 중심 모멘트 MRT가 효과적임을 보여준다.

뎁스가 높은 뚜껑 구동 캐비티 흐름에서는 레이놀즈 수(Re=1000)까지 안정적인 수렴을 확보했으며, 기존 BGK·SRT 모델이 발산하거나 비정상적인 회전 코어를 형성하는 경우에도 본 방법은 수렴성을 유지한다. 이는 충돌 연산자에서 각 차수별 모멘트를 독립적으로 이완시켜 고주파 모드(특히 비대칭 모드)의 과도 진동을 억제하기 때문이다.

연산 비용 측면에서는, 중심 모멘트 변환 및 역변환 과정이 추가되지만, 전체 시뮬레이션 시간은 기존 MRT 대비 10~15 % 정도만 증가한다. 반면, 안정성 향상으로 인해 동일한 정확도를 얻기 위해 필요한 격자 수가 30 % 이상 감소할 수 있어, 실질적인 효율성은 오히려 개선된다.

결론적으로, 중심 모멘트 기반 Cascaded LBM은 (1) 2차 정확도와 미분량까지의 고수준 수렴, (2) Galilean 불변성 보장으로 인한 비선형 흐름에서의 정확도 향상, (3) 기존 충돌 모델 대비 넓은 안정성 마진을 제공한다는 세 가지 핵심 장점을 가진다. 이러한 특성은 고레졸루션·고레인즈 수 흐름 시뮬레이션, 다중 물리 현상 결합 모델링 등에 유용하게 적용될 수 있다.