[논문 리뷰] Caustics and Intermittency in Turbulent Suspensions of Heavy Particles

Jérémie Bec, Luca Biferale외 3명|ArXiv.org|2009. 05. 08.

Particle Dynamics in Fluid Flows참고 문헌 1인용 수 5

한 줄 요약

이 연구는 관성 입자가 농축된 난류에서 속도 통계를 조사하여, 점성 하위영역에서 캐우스틱스에 의해 유도된 준이중프랙탈 스케일링과 관성 영역에서 척도에 의존하는 허더 지수를 규명한다. 고해상도 시뮬레이션은 캐우스틱스가 고차 구조 함수를 지배하며, 입자 운동은 척도에 의존하는 스토크스 수에 의해 결정됨을 보여주며, 이는 구름 내 조기 빗방울 형성과 관련이 있다.

ABSTRACT

The statistics of velocity differences between very heavy inertial particles suspended in an incompressible turbulent flow is found to be extremely intermittent. When particles are separated by distances within the viscous subrange, the competition between quiet regular regions and multi-valued caustics leads to a quasi bi-fractal behavior of the particle velocity structure functions, with high-order moments bringing the statistical signature of caustics. Contrastingly, for particles separated by inertial-range distances, the velocity-difference statistics is characterized in terms of a local Hölder exponent, which is a function of the scale-dependent particle Stokes number only. Results are supported by high-resolution direct numerical simulations. It is argued that these findings might have implications in the early stage of rain droplets formation in warm clouds.

연구 동기 및 목표

• 점성 하위영역에서 관성 입자 간 속도 차이 통계에서 캐우스틱스와 선호적 농축의 역할을 정량화하는 것.
• 비균일성과 다중스케일 역학이 따뜻한 구름 물리학에 관련된 희박한 분산에서 충돌 빈도에 미치는 영향을 규명하는 것.
• 캐우스틱스와 입자 관성의 상호작용으로 인해 속도 구조 함수가 이중프랙탈 또는 다중스케일 행동을 보이는지 조사하는 것.
• 척도에 의존하는 스토크스 수와 관성 영역 내 국소 스케일 지수 간의 연결 고리를 설정하는 것.
• 이러한 결과가 입자 충돌이 핵심적인 조기 단계의 빗방울 형성에 어떻게 관련되는지 평가하는 것.

제안 방법

• Rλ ≈ 400 및 Rλ ≈ 185에서 비압축성 나비에-스토크스 난류의 직접 수치 시뮬레이션(DNS)을 통해 전체 스케일 범위를 해상한다.
• 입자 궤적은 τ d²X/dt² = −dX/dt + u(X,t) 운동 방정식을 통해 계산되며, τ는 스토크스 시간이고, St = τ/τη는 스토크스 수이다.
• 거리 r로 분리된 입자 쌍에 대해 속도 구조 함수 Sₚ(r;St) = ⟨|Ṙ|ᵖ | R = r⟩를 계산하며, |Ṙ|은 종방향 속도 차이를 나타낸다.
• 충돌 커널을 계산하기 위해 고스트 입자 방법을 사용하며, κ(r;St) ∝ r^{D₂−1} S₁(r;St)로 표현되어 구조 함수와 충돌 빈도를 연결한다.
• 국소 허더 지수 ξₚ(r;τ) = d(ln Sₚ)/d(ln r)는 척도에 의존하는 스토크스 수 St(r) = ε^{1/3}τ / r^{2/3}의 함수로 분석된다.
• 다양한 r 및 St에서의 속도 차이 PDF 통계 분석은 거듭제곱 법칙 핵심부와 늘어진 지수 尾부를 보이며, 캐우스틱스 기여를 시사한다.

실험 결과

연구 질문

• RQ1 점성 하위영역에서 캐우스틱스는 관성 입자 간 속도 차이 통계에 어떻게 영향을 미치는가?
• RQ2 비균일성 있는 캐우스틱스로 인해 속도 구조 함수 Sₚ(r;St)가 어느 정도 이중프랙탈 또는 다중스케일 행동을 보이는가?
• RQ3 국소 스케일 지수 ξₚ(r;τ)는 자가유사 유동 모델이 예측한 것처럼 척도에 의존하는 스토크스 수 St(r)에만 의존하는가?
• RQ4 선호적 농축과 속도 차이의 상호작용은 희박한 분산에서 충돌 커널에 어떻게 영향을 미치는가?
• RQ5 이러한 결과는 입자 충돌이 핵심적인 따뜻한 구름의 조기 단계 빗방울 형성에 어떻게 관련되는가?

주요 결과

• 점성 하위영역(r ≪ η)에서 속도 구조 함수 Sₚ(r;St)는 정상 영역과 다중값 캐우스틱스 간의 경쟁으로 인해 준이중프랙탈 행동를 보이며, 고차 모멘트는 캐우스틱스 기여에 의해 지배된다.
• 비율화된 속도 차이 σ = Ṙ/R의 확률 밀도 함수(PDF)는 |σ| ≲ σ*(r)에서 거듭제곱 법칙 행동 ∝ |σ|^{-(1+ξ∞)}을 보이며, p > ξ∞일 경우 발산 모멘트를 나타낸다.
• |σ| ≳ σ*(r)에서의 늘어진 지수 꼬리 부분은 캐우스틱스 크기 통계와 관련되며 다중스케일을 유도할 수 있으나, 구조 함수에 대한 기여는 아직 확정되지 않았다.
• 관성 영역(η ≪ r ≪ L)에서 첫 번째 국소 허더 지수 ξ₁(r;τ)는 척도에 의존하는 스토크스 수 St(r)에 대해 그려졌을 때 수렴되며, 이는 St(r)가 r과 τ보다 지배적임을 확인한다.
• 작은 St(r)에서는 ξ₁(St(r))가 K41 예측값(ξ₁ = 1/3)에 수렴하여 트레이서 유사 행동을 나타내지만, 더 큰 St(r)에서는 크게 이격되어 비균일성이 증가한 것을 반영한다.
• 결과적으로, 거의 단일 크기의 수증기 방울로 이루어진 조기 단계 따뜻한 구름에서는 중력 침강이 지배적이기 전에도 캐우스틱스에 의해 입자 충돌 빈도가 크게 증가할 수 있음을 시사한다.