밀도 분산과 마하수의 새로운 관계 아음속 초음속 등온 난류

초록

이 연구는 등온 조건에서 아음속(M≈0.1)부터 초음속(M≈15)까지 다양한 마하수 구간의 난류를 1024³ 격자 해상도의 수치 실험으로 조사한다. 난류 구동은 완전 발산성(솔레노이달)과 완전 회전성(컴프레시브) 두 가지 극단적인 강제 방식을 적용하였다. 결과는 밀도 분포의 표준편차(σ_ρ)와 전체 마하수(M) 사이에 선형 관계가 존재함을 확인하고, 비례 상수는 구동 방식에 따라 달라진다. 또한, 압축성 속도 성분(M_c)와 σ_ρ 사이에도 강제에 거의 의존하지 않는 새로운 선형 관계가 제시된다.

상세 분석

본 논문은 등온 가스의 난류에서 밀도 PDF가 어떻게 변하는지를 정량적으로 규명하고자, 고해상도(최대 1024³) 유체역학 격자 시뮬레이션을 수행하였다. 마하수 M을 0.1에서 15까지 연속적으로 변화시키면서, 두 가지 극단적인 구동 메커니즘을 적용하였다. 첫 번째는 발산이 없는 솔레노이달(∇·F=0) 강제로, 이는 회전 운동을 주로 유발해 압축성 성분을 최소화한다. 두 번째는 회전이 없는 컴프레시브(∇×F=0) 강제로, 순수히 압축 파동을 생성해 밀도 변동을 극대화한다. 이러한 구동 방식은 기존 연구에서 제시된 σ_ρ–M 관계의 비례 상수 b가 1/3(솔레노이달)에서 1(컴프레시브)까지 크게 달라지는 원인으로 알려져 있다.

시뮬레이션은 무한히 큰 정적 격자와 주기적 경계조건을 사용했으며, 초기 조건은 균일 밀도와 무작위 위상만을 갖는 작은 진동으로 설정하였다. 시간에 따라 난류는 통계적 평형에 도달했으며, 이후 10~20개의 교차시간에 걸쳐 밀도와 속도 필드를 샘플링했다. PDF는 로그정규 형태를 보였으며, 표준편차 σ_ρ는 로그밀도 s=ln(ρ/ρ₀)의 표준편차와 직접 연관된다.

핵심 결과는 σ_ρ와 전체 마하수 M 사이의 선형 관계 σ_ρ = b M이 광범위한 M 구간(0.1 ≤ M ≤ 15)에서 유지된다는 점이다. 여기서 b는 구동 방식에 따라 달라지며, 솔레노이달 경우 b≈0.38, 컴프레시브 경우 b≈1.05로 측정되었다. 이는 기존 이론적 예측(b≈1/3~1)과 일치하면서도, 특히 초음속 영역에서 압축성 파동이 지배적인 경우 b가 1에 가까워짐을 확인한다.

가장 혁신적인 발견은 압축성 마하수 M_c, 즉 속도장의 압축성 성분만을 고려한 표준편차와 σ_ρ 사이의 관계이다. 압축성 성분은 Helmholtz 분해를 통해 얻어졌으며, M_c = ⟨v_c²⟩¹ᐟ² / c_s 로 정의된다. 저자들은 σ_ρ와 M_c 사이에 σ_ρ = b_c M_c 라는 새로운 선형 관계를 제시했으며, b_c≈0.9 정도로 거의 구동 독립적인 값을 보였다. 이는 압축성 모드가 밀도 변동을 직접적으로 주도한다는 물리적 직관과 부합한다. 또한, M_c가 작을 때(아음속)와 클 때(초음속) 모두 동일한 비례 상수를 유지함으로써, 기존의 전체 마하수 기반 관계보다 보편성이 높다.

이러한 결과는 별 형성 이론, 은하 내 가스 흐름, 그리고 수치적 서브그리드 모델링 등에 중요한 함의를 가진다. 특히, 밀도 분산을 직접적으로 예측할 수 있는 간단한 관계식이 제공되므로, 복잡한 난류 시뮬레이션 없이도 평균 마하수 혹은 압축성 마하수만으로 밀도 통계량을 추정할 수 있다. 또한, 구동 방식에 따른 b값 차이를 압축성 모드로 보정함으로써, 다양한 물리적 상황(예: 초신성 충격, 대규모 은하 충돌)에서도 일관된 모델링이 가능해진다.