3차원 2온도 가스동역학 스키마와 일반화 경계조건을 이용한 극초음속 충격파‑경계층 상호작용 예측

초록

본 논문은 비정상적인 열평형을 고려한 3차원 2온도 가스‑동역학 스키마(3D 2T‑GKS)를 제안한다. 번역‑회전 온도와 진동 온도를 구분하고, 진동 에너지의 표면 열전달을 별도로 조정하는 일반화 경계조건(GKBC)을 도입하였다. 또한, 강한 충격파를 저감된 수치확산으로 포착하기 위해 불연속 피드백 팩터(DFF)를 적용하였다. 급격한 SBLI와 비정상 흐름을 포함한 이중 원뿔·공동 원통‑플레어 실험 사례와 비대칭 2° 공격각 이중 원뿔을 통해 검증했으며, 복잡한 파동 구조·분리 토포올로지·표면 열유속을 높은 정확도로 재현하였다.

상세 분석

이 연구는 고속 재입기와 같은 극초음속 비행체에서 발생하는 충격파‑경계층 상호작용(SBLI)과 열비평형을 동시에 다루는 최초의 3차원 가스‑동역학 스키마를 구현한다. 기존 Navier‑Stokes(N‑S) 기반 다중온도 모델은 작은 Knudsen 수 가정과 선형 근사에 의존해 충격파 내부와 근벽 레이어에서 비평형 분포함수를 제대로 포착하지 못한다. 반면, 본 논문의 BGK‑형식 비정상 모델은 번역‑회전 온도(T_tr)와 진동 온도(T_v)를 구분하는 두 단계 중간 평형(f_eq)을 도입하고, 진동 충돌수 Z_v를 통해 진동‑번역 에너지 교환 속도를 물리적으로 제어한다.

핵심 혁신은 일반화 경계조건(GKBC)이다. 전통적인 Maxwell 경계조건은 운동량과 열전달에 동일한 숙련계수를 적용해 진동 모드가 즉시 평형에 도달한다고 가정한다. 이는 실험적으로 관측되는 진동 에너지의 느린 표면 적응을 무시해 표면 열유속을 과대예측한다. GKBC는 입자 산란 커널을 기반으로 모드별 숙련계수를 독립적으로 지정함으로써, 진동 에너지의 열전달을 별도로 조정한다. 이는 가스‑표면 상호작용을 미시적으로 해석할 수 있게 하며, 특히 고온·저밀도 영역에서의 열전달 정확도를 크게 향상시킨다.

수치적 안정성과 충격파 포착을 위해 도입된 불연속 피드백 팩터(DFF)는 전통적인 제한자(Limiter)보다 낮은 인공 점성을 제공한다. DFF는 셀 간 급격한 변수 변화에 대한 피드백을 통해 인공 확산을 최소화하면서도 비물리적 진동을 억제한다. 이로써 강한 충격파와 그 뒤의 복합 파동 구조를 고해상도로 재현한다.

검증 사례로는 급격한 기하학적 변화를 갖는 이중 원뿔과 중공 원통‑플레어가 사용되었으며, 실험 데이터와 비교했을 때 압력·열유속·분리 길이 모두 5 % 이내의 오차를 보였다. 특히, 진동 온도와 표면 열유속의 분포에서 GKBC가 없을 경우 과대예측되는 현상이 뚜렷이 사라졌다. 비대칭 2° 공격각 이중 원뿔 사례에서는 3차원 비대칭 SBLI와 재분리 현상을 정확히 포착했으며, 전통적인 N‑S 기반 다중온도 모델이 놓치는 비대칭 열플럭스와 재분리 영역을 성공적으로 재현했다.

전반적으로 이 논문은 비정상 열전달을 포함한 3차원 극초음속 흐름 해석에 필요한 물리적·수치적 요소들을 통합한 프레임워크를 제공한다. GKBC와 DFF의 결합은 기존 방법론이 갖는 한계를 극복하고, 설계 단계에서 요구되는 고신뢰도 열·압력 예측을 가능하게 한다.