충격파 내부 구조를 위한 최소 이산 속도 모델

초록

본 논문은 2차원 격자 위에서 가상의 입자가 제한된 이산 속도 집합을 이용해 이동하도록 설계한 모델을 제시한다. 25개와 33개의 이산 속도로 구성된 두 가지 모델이 Navier‑Stokes 방정식 수준의 정확도를 달성함을 보이며, 특히 33속도 모델은 기존 37속도 모델과 동등한 정밀도를 유지하면서 계산 비용을 절감한다. 이를 기반으로 충격파 튜브 시뮬레이션 및 다중 성분·복잡 기하 흐름을 수행해 내부 구조 진화를 성공적으로 재현한다.

상세 분석

이 연구는 격자 볼츠만 방정식(Lattice Boltzmann Equation, LBE)의 이산 속도 집합을 최소화하는 문제에 초점을 맞춘다. 기존 문헌에서는 2차원 정방격자에서 37개의 속도가 Navier‑Stokes 수준의 열·압축성 흐름을 재현하는 최소 조건으로 제시되었지만, 저자는 대칭성을 활용해 33개의 속도로 구성된 새로운 집합을 설계하였다. 속도 집합은 (0,0), c(1,0), c(2,0), c(3,0), c(1,1), c(2,2), c(4,4), c(2,1) 및 이들의 대칭 복사본으로 이루어지며, c≈0.818381 로 정해진다. 각 속도에 대응하는 가중치 w_i는 고차 가우시안 적분식에 의해 결정되어 등방성을 만족한다.
가중치와 속도는 Hermite 다항식 전개를 통해 Maxwell‑Boltzmann 분포의 근사 형태인 이산 평형 분포 f_i^eq 를 구성한다. 4차 Hermite 전개까지 포함함으로써 압력, 온도, 속도 등 2차 및 4차 모멘트를 정확히 보존한다. 25속도 모델은 5속도 기본 모델의 텐서 곱을 이용해 구성되며, 고마흐 수치에서는 고마흐 수치가 큰 마하수 영역에서 안정성이 떨어지는 반면, 33속도 모델은 보다 넓은 마하수 범위에서 안정적인 시뮬레이션을 제공한다.
충격파 튜브 실험에서는 1000×8 격자에 초기 조건을 설정하고, 좌측과 우측의 밀도·압력을 4배 차이로 두어 급격한 압축파를 생성하였다. ω=1 로 설정한 충돌 연산자는 점성 효과를 포함시키며, 결과는 33속도와 37속도 모델 모두에서 거의 동일한 밀도 프로파일을 보여준다. 또한 다중 성분 흐름을 10개의 수직 스트립으로 나누어 교차 배치함으로써, 두 성분이 충격파 내부에서 어떻게 섞이고 재배열되는지를 시각화하였다. 복잡한 기하 형태(200×200 격자)에서도 동일한 방법으로 초기 밀도 분포를 지정하고, 시간 전진에 따라 내부 구조가 자연스럽게 전파되는 모습을 확인하였다.
이러한 결과는 이산 속도 집합을 최소화하면서도 열·압축성 효과를 정확히 포착할 수 있음을 증명한다. 특히, 가중치와 속도 선택이 등방성 조건을 만족하도록 설계된 점이 모델의 수치적 안정성과 정확도에 크게 기여한다. 향후 3차원 확장 및 고차 비등방성 흐름에 대한 적용 가능성도 기대된다.