[논문 리뷰] Residence time of inertial particles in 3D thermal convection: implications for magma reservoirs

Vojtěch Patočka, Nicola Tosi외 1명|arXiv (Cornell University)|2022. 02. 25.

Particle Dynamics in Fluid Flows참고 문헌 57인용 수 14

한 줄 요약

이 연구는 레일라이히-베나르 대류에서 아인ertia 입자 체류 시간에 대한 3차원 수치 모델을 제시하며, 레이놀즈 수가 10⁷ 이상일 경우 입자 침강이 스토크스 속도와 유동 속도 비율에 따라 꾸준히 대류 지배에서 자유 낙하 역학으로 전이됨을 보여준다. 모델은 2차원 시뮬레이션에서 입자 유인 현상으로 인해 체류 시간을 과대평가함을 드러내며, 고체화된 마그마 침전층의 결정 크기 분포는 과거 마그마 대류의 활발함을 추론하는 데 활용될 수 있다.

ABSTRACT

The dynamic behavior of crystals in convecting fluids determines how magma bodies solidify. In particular, it is often important to estimate how long crystals stay in suspension in the host liquid before being deposited at its bottom (or top, for light particles). We perform a systematic 3D numerical study of particle-laden Rayleigh-Benard convection, and derive a robust model for the particle residence time. For Rayleigh numbers higher than 10^7, inertial particles' trajectories exhibit a monotonic transition from fluid tracer-like to free-fall dynamics, the control parameter being the ratio between particle Stokes velocity and the fluid velocity. The average settling rate is proportional to the particle Stokes velocity in both the end-member regimes, but the distribution of the residence times differs markedly from one to the other. For lower Rayleigh numbers (<10^7), an interaction between large-scale circulation and particle motion emerges, increasing the settling rates on average. Nevertheless, the mean residence time does not exceed the terminal time, i.e. the settling time from a quiescent fluid, by a factor larger than four. An exception are simulations with only a slightly super-critical Rayleigh number (~ 10^4), for which stationary convection develops and some particles become trapped indefinitely. 2D simulations of the same problem overestimate the flow-particle interaction - and hence the residence time - for both high and low Rayleigh numbers, which stresses the importance of using 3D geometries for simulating particle-laden flows. We outline how our model can be used to explain depth changes of crystal size distribution in sedimentary layers of magmatic intrusions that are thought to have formed via settling of a crystal cargo, and discuss how the micro-structural observations of solidified intrusions can be used to infer the past convective velocity of magma.

연구 동기 및 목표

• 3차원 열대류에서 아인ertia 입자(결정)가 침강하기 전에 얼마나 오랫동안 떠 있을지 이해하는 것.
• 특히 중간 스토크스 속도를 가진 입자에 대해 대규모 유체 순환 구조가 마그마 챔버 내 입자 침강에 어떤 영향을 미치는지 규명하는 것.
• 입자 랜든 대류를 모델링할 때 2차원 시뮬레이션의 한계와 고체화된 화산 텍스처 해석에 미치는 영향을 평가하는 것.
• 고체화된 마그마 내재 구조물에서 관측 가능한 결정 크기 분포와 입자 체류 시간을 연결하는 견고한 모델을 개발하는 것.
• 고체화된 내재물의 미세 구조적 특징이 과거 마그마 저수지의 대류 속도를 어떻게 추론하는 데 사용될 수 있는지 탐색하는 것.

제안 방법

• 레이놀즈 수 범위(10⁴ ~ 10⁹)와 프란틀 수(1 ~ 100)를 변화시켜 희박한 관성 입자 분산계를 포함한 레일라이히-베나르 대류에 대한 3차원 직접 수치 시뮬레이션을 수행한다.
• 시간에 따라 변화하는 유동장 내에서 관성 입자 운동을 기술하는 맥스레이-로빈슨 방정식을 사용하여 입자 궤적을 추적한다.
• 입자가 바닥(또는 떠오르는 입자인 경우 상단)에 침착하기 전까지 떠 있는 시간으로 체류 시간을 정의하며, 경계면에서의 입자 유량으로 계산한다.
• 입자 스토크스 속도와 루트 평균 제곱 유속(|vt|/urms)의 비율을 핵심 제어 변수로 삼아, 유동 지배(먼지처럼)와 침강 지배(돌처럼)의 두 영역으로 동역학을 분류한다.
• 3차원 결과를 2차원 시뮬레이션과 비교하여 입자 유인과 체류 시간 추정에 미치는 기하학적 편향을 정량화한다.
• 다양한 입자 유형과 유동 영역에서의 침강 곡선과 체류 시간 분포를 분석하여 고레이놀즈 수 대류에 대해 일반화 가능한 모델을 유도한다.

실험 결과

연구 질문

• RQ1 3차원 열대류에서 관성 입자의 체류 시간은 레이놀즈 수와 입자 스토크스 속도에 따라 어떻게 변화하는가?
• RQ2 대규모 순환 구조는 정적인 조건에 비해 마그마 챔버 내에서 입자 침강을 얼마나 늦추거나 가속화하는가?
• RQ3 왜 2차원 시뮬레이션은 입자 랜든 대류에서 체류 시간을 체계적으로 과대평가하는가? 이는 마그마 텍스처 해석에 어떤 영향을 미치는가?
• RQ4 고체화된 내재물의 침전층에서의 결정 크기 분포는 과거 마그마 저수지의 대류 활발함(즉, 유속)을 추론하는 데 사용될 수 있는가?
• RQ5 3차원 대류에서 대류 지배와 침강 지배의 입자 운동 간 전이의 본질은 무엇인가?

주요 결과

• 레이놀즈 수가 10⁷ 이상일 경우, 입자 체류 시간은 스토크스 속도와 유속 비율에 따라 꾸준히 대류 지배에서 자유 낙하 역학으로 전이된다.
• 고레이놀즈 수 영역에서는 평균 침강 속도가 입자 스토크스 속도에 비례하지만, 먼지처럼(유동 지배)과 돌처럼(침강 지배)의 한계에서 체류 시간 분포는 상당히 다르다.
• 레이놀즈 수가 10⁷ 이하일 경우, 대규모 순환은 평균적으로 침강 속도를 높이며, 그러나 평균 체류 시간은 정적인 유체에서의 최종 침강 시간을 넘지 않으며 최대 4배 이내이다.
• 매우 낮은 레이놀즈 수(~10⁴)에서는 정적인 대류가 발생하고 일부 입자는 무한정으로 갇혀 비정상적으로 긴 체류 시간을 경험한다.
• 2차원 시뮬레이션은 고레이놀즈 수와 저레이놀즈 수 모두에서 입자 유인 현상으로 인해 체류 시간을 과대평가한다. 이는 입자 랜든 대류 모델링에서 3차원 기하학의 중요성을 강조한다.
• 마그마 침전층의 결정 크기 분포는 고체화 과정 중에 지배적인 역학(유동 지배 대비 침강 지배)을 반영하므로, 과거 마그마 대류 속도를 추론하는 데 유용한 지표가 될 수 있다.