가변 프란틀 수와 비열비를 위한 이중분포 준평형 압축유동 모델

초록

**
본 논문은 이중분포(double‑distribution) 프레임워크와 준평형(quasi‑equilibrium) 접근법을 결합해, 프란틀 수와 비열비가 자유롭게 변하는 압축 유동을 정확히 재현하는 새로운 BGK‑형 모델을 제시한다. 고차 속도 격자와 정적 기준 프레임을 이용해 연속 방정식, 운동량 및 에너지 보존을 유지하면서 Navier‑Stokes‑Fourier(N‑S‑F) 방정식을 완전 복원한다. 열유동 Couette와 2차원 충격‑와소 상호작용을 포함한 검증 사례에서 높은 정확도, 수치적 안정성 및 갈릴레오 불변성을 확인하였다.

**

상세 분석

**
이 연구는 압축성 유동 시뮬레이션에서 가장 난제 중 하나인 ‘프란틀 수와 비열비를 동시에 자유롭게 조절할 수 있는 일관된 kinetic 모델’ 구축을 목표로 한다. 기존 BGK 모델은 프란틀 수가 1로 고정되고, 다원자 기체의 내부 자유도(비열비)를 반영하기 어려워 확장에 한계가 있었다. 저자들은 두 가지 핵심 아이디어를 결합하였다. 첫째, 이중분포(DDF) 구조를 도입해 첫 번째 분포 f는 질량·운동량을, 두 번째 분포 g는 내부 에너지(또는 전체 에너지)의 일부를 담당하도록 설계하였다. 이를 통해 비열비 γ를 자유롭게 설정할 수 있다. 둘째, 준평형(QE) 접근법을 적용해 ‘느린’ 고차 모멘트(N)를 별도의 quasi‑equilibrium 상태 f*·g*에 매핑하고, 두 단계(τ₁, τ₂)로 이완시키는 QE‑BGK 모델을 구축하였다. τ₁≤τ₂라는 이완 시간 계층은 H‑정리를 보장하고, 열전도와 점성계수를 각각 τ₂와 τ₁에 연결함으로써 프란틀 수 Pr=τ₁/τ₂(또는 그 역수) 를 자유롭게 조절한다.

수학적으로는 H‑함수 최소화 문제에 라그랑주 승수를 도입해 정확한 QEDF를 정의하고, 실제 구현에서는 Grad‑형 근사와 삼각 엔트로피 방법을 활용해 실용적인 표현을 얻었다. 모델의 연속성·운동량·에너지 보존은 M‑계와 N‑계의 정의에 의해 엄격히 유지된다.

이산화 단계에서는 고차 Hermite 다항식 기반의 고차 속도 격자(D3Q27, D3Q45 등)를 사용해 충분한 순간(모멘트) 재현성을 확보하였다. 정적 기준 프레임을 유지함으로써 격자 이동에 따른 복잡성을 최소화하고, 확장 차수와 사분면(Quadrature) 차수가 모델 정확도에 미치는 영향을 명확히 분석하였다. 특히, Mach ≈ 2~3 수준까지의 중간 초음속 흐름에서도 Galilean 불변성을 유지함을 수치 실험으로 검증하였다.

검증 사례로는 (i) 열 Couette 흐름에서 온도와 속도 프로파일이 이론 해와 거의 일치함을 보였으며, (ii) 2‑D 충격‑와소 상호작용에서 복잡한 비선형 파동 구조와 열전달이 정확히 재현되어, 모델이 강한 비선형성·불연속성을 포함한 상황에서도 안정적임을 입증했다. 또한, 모든 프란틀 수(Pr ∈