[논문 리뷰] Wet granular bed eroded by a dry granular flow

Lama Braysh, Patrick Mutabaruka외 2명|Archimer (Ifremer)|2026. 02. 23.

Granular flow and fluidized beds인용 수 0

한 줄 요약

본 논문은 3D DEM 시뮬레이션을 이용해 응집성 젖은 입상층이 위에 흐르는 마른 입상 흐름에 의해 침식되는 현상을 연구하고, 느린 침식과 빠른 침식 체계를 식별하며 관성수와 응집 지수를 기반으로 한 침식 속도 스케일링 법칙을 제안한다.

ABSTRACT

Using three-dimensional Discrete Element Method (DEM) simulations, we investigate the erosion dynamics of a cohesive bed composed of wet spherical particles subjected to the shear flow of an overlying non-cohesive granular layer. Cohesion is modeled through a capillary attraction law, where the erosion process is governed by the irreversible rupture of liquid bridges at the interface. By systematically varying the liquid-vapor surface tension and the inclination angle of the bed, we analyze the influence of cohesive strength and flow intensity on the mass entrainment rate. Our results identify two distinct erosion regimes: a slow, stochastic regime driven by granular temperature fluctuations, and a fast, collective regime characterized by a global mechanical instability of the interface. We propose a robust scaling law for the surface erosion rate in the fast erosion regime, based on the interplay between two dimensionless parameters: the inertial number ($I$) and the cohesion index ($ξ$). This framework reveals that the threshold for the fast erosion regime is determined by the ratio of the geometric mean of the driving stresses (kinetic and pressure) to the cohesive resistance. These findings provide a comprehensive description of the coupling between inertial and capillary forces, offering a predictive tool for the stability of cohesive interfaces in sheared granular flows.

연구 동기 및 목표

• 경사면 위의 흐르는 더러운 입자 층에 노출된 응집성 젖은 입상층의 침식 역학을 조사한다.
• 응집 강도와 흐름의 강도가 표면에서의 질량 포섭에 어떤 영향을 미치는지 이해한다.
• 침식 체계의 특성을 규명하고 빠른 체계에서의 표면 침식 속도에 대한 예측 가능한 스케일링을 개발한다.
• 관성력과 모세관 힘을 침식 임계치 및 속도와 연결하는 체계를 개발한다.

제안 방법

• 모듈러 공간에서의 캡럴리 결합으로 응집성 층을 모델링하고 pendular-state 액체 다리 모델을 이용해 응집을 나타낸다.
• 속도-베레트(velocity-Verlet) 적분 스킴을 쓰는 3차원 Discrete Element Method (DEM) 시뮬레이션을 사용한다.
• 끊김 거리 d_rupt와 되돌릴 수 없는 결합 파손으로 이진, 거리 의존 결합으로 모세관 힘을 구현한다.
• 점성 스트레스가 무시 가능한 상태에서 작동하여 관성-응집 간 경쟁에 집중한다.
• 침식 체계를 특징짓기 위한 무차원 수 I(관성 수)와 xi(응집 지수)를 정의하고 계산한다.
• 빠른 체계에서의 침식 속도 K_e를 위한 스케일링 법칙을 도출하고 검증하며 I/xi에 대한 마스터 커브를 제시한다.

실험 결과

연구 질문

• RQ1 응집 인터페이스에서 느린 침식과 빠른 침식 체계 간의 전이에 무엇이 작용하는가?
• RQ2 I(관성 수)와 xi(응집 지수)가 빠른 체계에서의 표면 침식 속도 K_e를 어떻게 제어하는가?
• RQ3 경합하는 관성력과 모세관력에 기초해 질량 포섭을 예측하는 스케일링 법칙을 제시할 수 있는가?
• RQ4 기울기 각도와 액체 표면 장력이 침식 역학과 임계값에 어떤 영향을 미치는가?

주요 결과

• 두 가지 침식 체계가 확인된다: 키네틱한 변동성에 의해 좌우되는 느린 침식과 인터페이스의 글로벌 기계적 불안정에 의해 좌우되는 빠른 침식.
• 빠른 침식 체계는 I sqrt(xi)/dot_gamma를 K_e에 대해 나타낼 때 I/xi를 인자로 취하면 데이터가 분석에 의해 한 축으로 모이는 것을 보인다, 이는 I/xi가 전이와 속도를 결정한다는 것을 나타낸다.
• 빠른 체계의 마스터 커브는 (I sqrt(xi)/dot_gamma) K_e = A[(I/xi) − (I/xi)^*]^κ의 거듭제곱 형태를 따르며, κ ≈ 1.97, A ≈ 6, (I/xi)^* ≈ 0.02이다.
• 관찰된 임계값 (I/xi)^*와 안정성 분석에서 예측된 a/α 사이의 이론적 일치가 발견되며, 매개변수 추정치 a ≈ 0.095와 α ≈ 5를 사용하면 (I/xi)^* ≈ 0.019이 된다.