파라메트릭 격자볼츠만법을 통한 고효율 이소열·열압축 흐름 시뮬레이션

초록

본 논문은 라그랑주 보간 다항식 계수를 이용해 이산 평형 분포를 구성하고, 이들을 매개변수화함으로써 기존 LBGK 모델보다 넓은 속도·온도 범위와 향상된 안정성을 제공하는 새로운 파라메트릭 격자볼츠만법(PLBM)을 제안한다. 3·4·5개의 온-격자 속도를 사용해 이소열 및 열압축 Navier‑Stokes 방정식을 정확히 회복하고, 충격파 튜브와 켈빈‑헬름홀츠 불안정성 사례에서 뛰어난 성능을 보인다.

상세 분석

이 논문은 격자볼츠만법(LBM)의 핵심인 이산 평형 분포 f_eq 를 새로운 수학적 틀로 재구성한다. 저자는 이산 속도 v_i 와 크로네커 델타 δ_ik 를 통과하는 라그랑주 보간 다항식의 계수 c_ij 를 이용해 r_i (재분배 규칙)를 정의하고, 이를 맥스웰‑볼츠만 분포의 모멘트 μ_n 과 일치시키는 q 개의 속도 집합을 선택한다. 이때 q 개의 속도는 n* = q‑1 이라는 관계를 만족하도록 설계되어, n* = 3이면 이소열 압축 Navier‑Stokes 방정식을, n* = 4이면 열압축 Navier‑Stokes 방정식을 정확히 회복한다.

특히, 3개의 속도(0, ±√ζ θ₀) 를 사용한 경우 ζ=3이면 기존 LBGK와 동일하고, ζ=4이면 속도 u 에 대한 양의 r_i 조건이 |u|≤√3 까지 확대된다. 이는 LBGK의 |u|≤√2 보다 넓은 안정 영역을 의미한다. 저자는 이론적 분석과 함께 충격파 튜브 실험을 수행해, ζ=4 모델이 LBGK와 엔트로피 모델보다 진동 없이 정확한 전단 속도와 압력 프로파일을 제공함을 입증한다. 또한, 4개의 속도를 이용한 온-격자 모델은 기존 오프‑격자 Gauss‑Hermite 방식보다 적은 5 개의 온-격자 속도로 열압축 흐름을 시뮬레이션할 수 있다. 이는 b=2a 와 같은 비율 조정을 통해 온-격자 속도를 격자점 간 이동이 가능한 “온‑격자” 형태로 만들 수 있기 때문이다.

2차원 켈빈‑헬름홀츠 불안정성 실험에서는 ζ=4 파라메트릭 모델을 텐서곱으로 확장해 D2Q9와 비교했을 때, 더 긴 시간 동안 안정적으로 소용돌이를 유지하고, 격자 해상도에 대한 수렴성 검증에서도 LBGK보다 낮은 오차를 보였다. 이는 r_i 의 양의성 확보와 모멘트 정확도가 향상된 결과로, 고속·고온 흐름에서의 수치 확산과 비선형 왜곡을 최소화한다는 장점을 시사한다.

결론적으로, PLBM은 매개변수 ζ 또는 a, b 등을 조정함으로써 원하는 정확도와 안정성을 설계 단계에서 선택할 수 있게 하며, 기존 LBGK·엔트로피·Gauss‑Hermite 기반 모델을 포괄하는 일반화된 프레임워크를 제공한다. 이는 복잡한 열·압축 현상을 저비용·고효율로 구현하려는 CFD 연구에 중요한 전진을 의미한다.