압축성 난류에서 입자 쌍 분산의 간헐성 및 비보편성

초록

이 연구는 2차원 등온 압축성 난류에서 라그랑지 입자 쌍의 확산을 조사한다. 입자 거리의 두 배가 되는 ‘doubling time’과 절반이 되는 ‘halving time’의 음의 모멘트 지수(χ_D^p, χ_H^p)가 비선형적으로 p에 의존함을 발견했으며, halving time 지수는 다중프랙털 모델이 예측하는 보편적 관계를 만족하지만, doubling time 지수는 강제 방식(솔레노이달 vs. 이로테이션)과 마하 수(Ma)에 따라 달라진다. 솔레노이달 강제에서는 χ_D^p가 솔레노이달 속도 구조함수의 스케일 지수(ζ_s^p)와 연결되고, 이로테이션 강제에서는 알려진 관계가 성립하지 않는다. 이러한 결과는 천체 물리학적 가스 흐름에서 혼합과 화학 반응을 재고할 필요성을 시사한다.

상세 분석

본 논문은 압축성 난류에서 라그랑지 입자 쌍의 확산을 기존 리처드슨 법(Richardson’s law)의 일반화 관점에서 접근한다. 2차원 등온 이상기체를 대상으로, 대규모 솔레노이달(소용돌이성) 및 이로테이션(압축성) 강제력을 각각 적용한 네 가지 DNS(시뮬레이션) 케이스(S1, S2, C1, C2)를 수행하였다. 마하 수는 전이음속(Ma≈1)과 초음속(Ma>1) 두 범위로 설정했으며, 각 케이스에서 4096²개의 라그랑지 입자를 삽입해 입자 쌍의 거리 R(t)를 추적하였다.

전통적인 ⟨R²(t)⟩의 시간 의존성은 관측했지만, 관측 가능한 스케일 구간이 너무 짧아 정확한 지수 추정이 어려웠다. 따라서 ‘exit time’ 접근법을 채택하였다. 특정 거리 R에 대해 R가 3R/2(두 배) 혹은 3R/4(절반)으로 처음 도달하는 시간을 각각 τ_D(R)와 τ_H(R)라 정의하고, 그 음의 모멘트 ⟨τ_D^{-p}⟩∝R^{-χ_D^p}, ⟨τ_H^{-p}⟩∝R^{-χ_H^p}의 스케일 지수를 측정했다.

다중프랙털 모델에 따르면 χ_{D,H}^p = p – ζ^p (ζ^p는 전체 속도 차이의 구조함수 지수)이어야 한다. 실험 결과는 halving time에 대해서는 이 관계가 잘 맞았으며, χ_H^p는 강제 방식과 Ma에 무관한 보편성을 보였다. 반면 doubling time은 두 가지 중요한 비보편성을 나타냈다. 첫째, 솔레노이달 강제(S1, S2)에서는 χ_D^p = p – ζ_s^p 로, 전체 속도가 아니라 솔레노이달 성분(u_s)의 구조함수 지수와 일치했다. 이는 입자 거리 증가가 소용돌이 구조에 의해 주도됨을 의미한다. 둘째, 이로테이션 강제(C1, C2)에서는 χ_D^p가 어떠한 알려진 ζ^p와도 일치하지 않았으며, Ma에 대한 의존성도 관측되지 않았다. 즉, 압축성 성분(u_c)과 양의 발산 영역(∇·u>0)이 거리 증가에 복합적으로 작용한다는 것을 시사한다.

다중프랙털 모델의 Legendre 변환을 이용해 동적 지수 χ_{D,H}^p와 구조 차원 D_{D,H}(h)를 연결하였다. D_H(h)는 모든 케이스에서 h→0일 때 1에 수렴했으며, 이는 1차원 충격파가 halving time을 지배함을 의미한다. 반면 D_D(h)는 솔레노이달 강제에서 Ma가 커질수록 h≈0 근처에서 1에 가까워졌는데, 이는 충격파 주변에 형성되는 얇은 고소용돌이 패치가 거리 증가에 기여함을 보여준다. 이로테이션 강제에서는 D_D(h)도 전반적으로 1에 가까워, 압축성 영역 자체가 주요 역할을 함을 나타낸다.

결과적으로, 압축성 난류에서는 입자 쌍의 확산이 두 종류의 exit time에 대해 서로 다른 동적 멀티스케일링을 보이며, 기존 다중프랙털 모델이 포착하지 못하는 비보편적 현상이 존재한다. 이는 천체 물리학적 환경(예: 분자 구름, 원시 원반)에서 혼합 효율과 화학 반응 속도를 재평가해야 함을 암시한다.