[논문 리뷰] Intermittent distribution of tracers advected by a compressible random flow

Jérémie Bec, Krzysztof Gawędzki외 1명|ArXiv.org|2003. 10. 13.

Fluid Dynamics and Turbulent Flows참고 문헌 1인용 수 3

한 줄 요약

이 논문은 압축 가능한 랜덤 유동에서 트레이서 밀도의 다중분포 스케일링과 유동의 유한시간 리아프노프 지수(편차율)의 극단값 간의 정량적 연관성을 수립한다. 압축 가능한 크라이히난 모델을 사용하여 질량 모멘트의 스케일링 지수 ξₙ에 대한 해석적 표현을 유도하며, 간헐성은 희귀한 강한 편차를 보이는 유동 실현에서 기인하며, 수치 시뮬레이션을 통해 예측된 스케일링과 임계 순서 n_cr에서의 포화가 확인된다.

ABSTRACT

Multifractal properties of a tracer density passively advected by a compressible random velocity field are characterized. A relationship is established between the statistical properties of mass on the dynamical fractal attractor towards which the trajectories converge and large deviations of the stretching rates of the flow. In the framework of the compressible Kraichnan model, this result is illustrated by analytical calculations and confirmed by numerical simulations.

연구 동기 및 목표

• 압축 가능한 무작위 유동에서 트레이서 밀도의 다중분포 성질을 기술하기 위해.
• 트레이서 질량 모멘트의 스케일링 지수 ξₙ를 흐름의 편차율(유한시간 리아프노프 지수)의 극단값과 연관짓기 위해.
• 안정화된 헨슈텔-프로카치아 스펙트럼 HP_an(n)과 동적 아트레터의 통계적 성질 간의 연결 고리를 수립하기 위해.
• 장시간 평균을 고려한 크라이히난 모델의 수치 시뮬레이션을 통해 이론적 예측을 검증하기 위해.

제안 방법

• 유동장이 가우시안이고 시간에 대해 백색임을 가정하는 압축 가능한 크라이히난 모델을 사용하여 피동 트레이서의 이송을 모델링한다.
• 라그랑주 경로가 수렴하는 동적 아트레터 상의 불변 측도로 트레이서 밀도 ρ*를 정의한다.
• 소규모 질량 분포를 기술하기 위해 국소 질량 m_r(x)와 그 스케일링 지수 h_r(x)를 도입한다.
• 스케일링 ⟨m_r^n⟩ ∼ r^ξₙ에 기반하여 안정화된 헨슈텔-프로카치아 스펙트럼 HP_an(n)을 유도하며, 이는 비율 함수 S(h)의 레전드르 변환과 관련된다.
• ξₙ를 변형 텐서의 고유값의 공동 분포의 생성 함수를 통해 편차율 σ_i의 극단값과 연관지운다.
• 최대 10⁵개의 트레이서와 10⁸개의 회전 시간을 사용한 수치 시뮬레이션을 수행하여 모멘트를 계산하고 스케일링 지수를 검증한다.

실험 결과

연구 질문

• RQ1 압축 가능한 유동에서 트레이서 질량 모멘트의 스케일링 지수 ξₙ는 어떻게 유한시간 리아프노프 지수의 극단값과 관련이 있는가?
• RQ2 ξₙ가 포화되는 임계 순서 n_cr는 무엇에 의해 결정되며, 이는 압축성에 따라 어떻게 달라지는가?
• RQ3 트레이서 측도의 다중분포 스펙트럼 f(h)는 편차율 통계로부터 유도될 수 있는가?
• RQ4 모멘트 ⟨m_r^n⟩는 r에 대해 어떻게 스케일링되는가? 이는 트레이서 분포의 간헐성과 비정규성에 어떤 의미를 갖는가?

주요 결과

• 트레이서 분포 질량 모멘트의 스케일링 지수 ξₙ는 흐름의 편차율의 극단값에 의해 결정되며, 정밀한 해석적 관계가 도출되었다.
• 압축 가능한 크라이히난 모델에서 지수 ξₙ는 n에 대해 비선형적 의존성을 보이며, 이는 m_r의 분포에서 강한 간헐성과 비정규 꼬리의 존재를 시사한다.
• ξₙ가 포화되는 임계 순서 n_cr는 동축 극한에서 포아송 연산자의 영모드 해의 적분 가능성에 의해 결정되며, ℘ < 1/6일 경우 n_cr = 1/(6℘), ℘ > 1/6일 경우 n_cr = 1/(6(1−℘))이다.
• 최대 10⁸개의 회전 시간을 사용한 수치 시뮬레이션은 이론적 예측을 확인하였으며, 해석적 스케일링 지수와의 양호한 일치 및 n_cr에서의 포화를 보였다.
• 정수 n에 대해 모멘트 ⟨m_r^n⟩는 트레이서 밀도의 정적 상관 함수 F_{n+1}과 연결되며, 이는 두 번째 순서 미분 연산자 M_{n+1}^†의 영모드이고, 그 동차도가 ξₙ를 제공한다.
• 고정된(fixed) 및 안정화된 헨슈텔-프로카치아 스펙트럼은 다르며, HP_qu(n) ≥ HP_an(n)이며, 안정화된 스펙트럼은 ξₙ의 스케일링에 의해 완전히 결정된다.