극초음속 경계층 두 번째 맥 모드 불안정성의 분자 가스 동역학 연구

초록

본 논문은 Mach 6, Reynolds 수 1.1 × 10⁷인 평판 경계층을 Direct Simulation Monte Carlo(DSMC) 방법으로 시뮬레이션하여 두 번째 맥 모드 불안정성의 자연 발생을 포착한다. 전력 스펙트럼 밀도(PSD) 분석으로 200–400 kHz의 우세 주파수를 확인하고, 동적 모드 분해(DMD)로 파장·위상속도가 선형 안정성 이론(LST)과 일치하는 음향 파동 패킷을 추출한다. 또한 300 kHz에서 진동하는 음향 진동 표면(AVS)을 이용해 불안정 주파수에서 파동을 증폭시키고, 500 kHz에서는 감쇠시키는 제어 효과를 입증한다. 결과는 DSMC가 고레인즈수 하이퍼소닉 전이 현상의 선형 메커니즘을 정확히 재현할 수 있음을 보여준다.

상세 분석

본 연구는 고속 항공우주 분야에서 핵심적인 두 번째 맥 모드 불안정성을 최초로 DSMC 기반 분자 가스 동역학 시뮬레이션으로 정량화한 점에서 의의가 크다. 먼저 저자들은 Mach 6, Re = 1.1 × 10⁷·m⁻¹이라는 실험실 수준의 고레인즈수를 설정하고, 변수 경도 구(VHS) 모델과 주요 빈도 충돌 스킴을 적용해 1.5 × 10⁹ 입자를 1 ns 타임스텝으로 0.6 ms 동안 진행하였다. 이러한 설정은 평균 자유행로(λ ≈ 6.1 × 10⁻⁷ m)와 격자 크기(Δx = 4 × 10⁻⁵ m)를 충분히 해상화하여, 전이 영역의 미세한 속도·온도 변동을 포착한다.

시간‑정확한 DSMC 데이터는 셀별 온도와 속도 시계열을 추출해 PSD 분석에 투입했으며, 200–400 kHz 대역에서 명확한 피크가 나타났다. 이 피크는 전통적인 LST가 예측한 불안정 구간(ω·δ≈0.3–0.6, 여기서 δ는 경계층 두께)과 일치한다. 특히, 불안정 영역은 전단층 근처에서 국한되며, DSMC가 제공하는 입자‑레벨 통계가 “자연 발생” 잡음으로 작용해 외부 강제 없이도 두 번째 맥 모드가 성장함을 확인했다.

동적 모드 분해(DMD)는 수천 개의 전압·속도 스냅샷을 이용해 전역적인 고유 모드를 추출했으며, 각 모드의 복소수 고유값으로부터 성장률과 위상속도를 계산했다. DMD 결과는 두 개의 주요 모드가 존재함을 보여준다. 첫 번째 모드는 파장 λ≈2 mm, 위상속도 cₚ≈1.5·U_∞(∞)에 해당하는 전형적인 음향 파동이며, 성장률이 양수인 구간이 LST가 제시한 불안정 구역과 정확히 겹친다. 두 번째 모드는 보다 짧은 파장과 낮은 위상속도를 가지며, 성장률이 거의 영에 가까워 비활성화된 고주파 잡음으로 해석된다. 이러한 DMD 기반 모드 구분은 기존의 선형 전역 안정성 해석을 데이터‑드리븐 방식으로 보완한다는 점에서 혁신적이다.

제어 실험에서는 x = 60 mm 지점에 1 mm 폭의 음향 진동 표면(AVS)을 배치하고, 300 kHz(불안정)와 500 kHz(안정) 두 주파수로 강제 진동을 가했다. 300 kHz에서는 AVS에서 방출된 음파가 기존 두 번째 맥 모드와 위상 동조를 이루어 파동 패킷이 하류에서 증폭되는 현상이 관찰되었으며, N‑factor 분석에서 추가적인 0.8 ≈ 1.0의 성장률 증가가 기록되었다. 반면 500 kHz에서는 파동이 경계층을 통과하면서 급격히 감쇠되어, 전반적인 PSD 레벨이 감소하고 DMD 모드의 성장률이 음수로 전환되었다. 이는 DSMC가 비평형 효과와 표면‑유체 상호작용을 자연스럽게 포함함으로써, 이론적 제어 전략(주파수 선택적 강제)의 실제 구현 가능성을 입증한다는 의미다.

전체적으로 본 논문은 (1) DSMC가 고레인즈수·고마하 수 하이퍼소닉 경계층에서 선형 불안정성을 정확히 재현할 수 있음을, (2) PSD와 DMD를 결합한 데이터 분석 파이프라인이 불안정 파동의 스펙트럼·공간 구조를 정량화하는 데 유효함을, (3) AVS와 같은 표면 진동을 이용한 주파수 선택적 제어가 DSMC 기반 시뮬레이션에서 실증적으로 검증될 수 있음을 보여준다. 이러한 결과는 향후 고속 재진입체의 전이 억제 설계와, 희소 가스·비평형 현상이 중요한 우주 비행체의 열·구조 해석에 DSMC를 활용할 근거를 제공한다.