고속 압축성 흐름을 위한 3차원 라티스 볼츠만 모델

초록

본 논문은 Kataoka와 Tsutahara가 제시한 3차원 라티스 볼츠만(LB) 모델을 고속 압축성 흐름에 적용하도록 개선하였다. 비진동(NND) 스킴으로 대류항을 이산화하고, 추가 점성·확산 항을 도입해 수치적 안정성을 높였다. von Neumann 안정성 기준을 만족하도록 파라미터를 선정하고, 라크스 충격관 및 고마하 충격관, 기포·액적에 대한 충격파 상호작용 등 다양한 벤치마크 실험을 통해 서브소닉부터 초음속까지의 흐름을 정확히 재현함을 확인하였다.

상세 분석

이 연구는 기존 2차원·저속 전용 LB 모델의 한계를 극복하고, 3차원 고속 압축성 유동을 직접 시뮬레이션할 수 있는 프레임워크를 제시한다. 핵심은 두 가지 수치적 개선에 있다. 첫째, 대류항을 비진동(NND) 스킴으로 이산화함으로써 충격파와 같은 급격한 불연속면에서 발생하는 Gibbs 현상을 효과적으로 억제한다. NND는 자유 파라미터가 없으며, 고차 정확도와 비진동 특성을 동시에 만족한다는 점에서 기존 Lax‑Wendroff·MUSCL 등과 차별화된다. 둘째, 점성·열전도와 같은 물리적 확산 효과를 모델링하는 동시에, 인위적인 인공 점성 항을 추가하여 수치적 안정성을 강화한다. 이 인공 점성은 실제 물리적 점성계수와는 별개로, 고마하 수치에서 발생할 수 있는 비물리적 진동을 완화하도록 설계되었다.

안정성 분석에서는 von Neumann 선형 안정성 이론을 적용해 파라미터(격자 간격 Δx, 시간 간격 Δt, 충돌 연산자 τ 등)의 허용 범위를 도출하였다. 특히, 충돌 연산자의 이완 시간 τ를 1보다 크게 설정함으로써 점성 항이 충분히 작용하도록 하였으며, 이는 고마하 흐름에서의 비선형 안정성에도 긍정적인 영향을 미친다.

검증 사례로는 (i) 라크스 충격관 문제와 고마하(>5) 충격관 문제, (ii) 충격파가 액적·기포에 충돌하는 복합 현상 등을 선택하였다. 라크스 충격관에서는 압력·밀도·속도 프로파일이 정확히 해석해와 일치했으며, 고마하 충격관에서는 충격 전후의 급격한 변수 변화에도 불구하고 비진동 스킴과 인공 점성 덕분에 수치 발산 없이 안정적인 결과를 얻었다. 액적·기포 실험에서는 충격파에 의해 발생하는 변형, 파동 전파, 그리고 기포 붕괴 현상이 물리적 기대와 일치했으며, 기존 LB 모델이 보였던 비대칭 및 과도한 진동이 현저히 감소하였다.

결과적으로, 제안된 3차원 고속 압축성 LB 모델은 물리적 정확도와 수치적 안정성을 동시에 만족시키며, 기존 LB 방법이 다루기 어려웠던 초음속·극초음속 흐름을 효율적으로 시뮬레이션할 수 있음을 입증한다.