다차원 이산 볼츠만 프레임워크를 통한 압축성 유동 시뮬레이션

초록

본 논문은 1차원 이산 볼츠만 모델에 추가 자유도와 고대칭 이산 속도 집합을 도입해 비압축성 비정상 효과를 최소화하고, 연산자 분할 기법으로 1D, 2D, 3D 흐름을 하나의 프레임워크에서 일관되게 시뮬레이션할 수 있음을 보여준다. Sod·Lax 충격파, 평행 이동, 음파 전파 등 네 가지 벤치마크에서 높은 정확도와 안정성을 검증하였다.

상세 분석

이 연구는 기존 DBM(Discrete Boltzmann Method)의 차원 의존성을 근본적으로 재구성한다는 점에서 혁신적이다. 먼저, 1차원 D1V5 모델을 선택하고, 자유도 I=4(회전·진동) 를 추가함으로써 비열비 γ=(D+I+2)/(D+I) 를 자유롭게 조절할 수 있다. 이는 고정된 격자와 속도 집합을 유지하면서도 다양한 물리적 기체를 모델링할 수 있게 한다.
속도 집합은 (0, ±1, ±5) 로 구성되며, 대칭성이 뛰어나 Galilean 불변성을 보장한다. 특히, 이산 속도와 연속 속도 사이의 변환이 선형이므로, 흐름이 전체 좌표계에서 이동해도 수치적 편향이 발생하지 않는다.
연산자 분할은 Godunov 1차 분할을 채택했으며, 3차원 볼츠만 방정식을 x‑, y‑, z‑축 순서로 독립적인 1차원 흐름으로 분해한다. 각 단계는 (i) 경계조건 적용, (ii) 로컬 평형분포 f_eq 계산, (iii) BGK 충돌‑전파 연산, (iv) 매크로 변수 복원 순으로 진행된다. 이 구조는 코드 재사용성을 극대화하고, 차원 확대 시 추가적인 격자·속도 설계 없이 기존 1D 모듈만으로 2D·3D 시뮬레이션을 수행할 수 있게 한다.
수치 검증에서는 고해상도(Δx=2×10⁻⁴, Δt=5×10⁻⁶)와 5000 격자를 사용해 Sod 및 Lax 충격파 문제를 풀었으며, 밀도·압력·속도·온도 프로파일이 정확한 라인어 솔루션과 거의 일치한다. 이는 충격 전파와 희박파(rarefaction) 모두를 정확히 포착함을 의미한다.
평행 이동 테스트에서는 대각선 속도 (0.5, 0.5) 로 전체 유체 블록을 이동시켰고, 시간 t=0.4 에서 이론적 이동 거리 Lx와 일치하는 결과를 얻어 Galilean 불변성을 실증했다.
음파 전파 실험에서는 1D, 2D, 3D 각각에서 작은 압력 교란을 초기화하고, 전파 속도 v_s=√(γT) 와 일치하는 전파 위치를 확인했다. 특히 2D·3D에서 원형·구형 파동이 정확히 방사형으로 전파되는 점은 연산자 분할이 차원 간 상호작용을 적절히 처리함을 보여준다.
한계점으로는 2D·3D 시뮬레이션에서 비정상(NE) 효과가 충분히 반영되지 않아 복잡한 난류·충격 구조에 대한 정밀도는 떨어질 수 있다. 향후 고차 순간(3차, 4차)까지 보존하는 확장형 DBM이나 다중 속도 집합을 도입하면 이 문제를 보완할 수 있을 것이다.