[論文レビュー] Wet granular bed eroded by a dry granular flow

Lama Braysh, Patrick Mutabaruka+2 more|Archimer (Ifremer)|Feb 23, 2026

Granular flow and fluidized beds被引用数 0

ひとこと要約

この論文は3D DEMシミュレーションを用いて、乾燥粒子流による湿潤結合ゲル床の侵食を研究し、遅い侵食と速い侵食の2つの領域を同定し、慣性数と結合指標に基づく侵食速度のスケーリング則を提案する。

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.

研究の動機と目的

• 傾斜面上の乾燥粒子の重ね覆いが作用する湿潤結合床の侵食ダイナミクスを調査する。
• 界面での結合強度と流れの強さが質量取り込みにどう影響するかを理解する。
• 侵食の領域を特徴づけ、速い領域の表面侵食率の予測的スケーリングを開発する。
• 慣性力と毛細力を侵食閾値と速度へ結びつける枠組みを開発する。

提案手法

• ペンダント型液橋モデルによる毛細力による結合床のモデリング。
• 速度-Verlet積分スキームを用いた3次元Discrete Element Method (DEM)シミュレーション。
• 破断距離d_ruptを用いたバイナリの距離依存結合として毛細力を実装し、不可逆な結合破壊を含む。
• 慣性と結合の競合に焦点を当て、粘性せん断応力を無視する条件で運用。
• 侵食領域を特徴づける無次元量I（慣性数）とξ（結合指標）を定義・計算。
• 速い領域の侵食速度K_eのスケーリング則を導出・検証し、I/ξに対するマスタカーブを提示。

実験結果

リサーチクエスチョン

• RQ1 結合界面での遅い侵食と速い侵食の転換は何によって支配されるのか。
• RQ2 速い領域において慣性数Iと結合指標ξは表面侵食率K_eをどう支配するのか。
• RQ3 競合する慣性力と毛細力に基づく質量取り込みを予測するスケーリング則を構築できるか。
• RQ4 傾斜角と液面張力は侵食ダイナミクスと閾値にどのように影響するか。

主な発見

• 2つの侵食領域を同定：運動エネルギーのゆらぎによる遅い侵食と、界面の全体的な機械的不安定性による速い侵食。
• 速い侵食領域では(I sqrt(xi)/dot_gamma) K_eをI/ξでプロットしたときデータが崩れ、I/ξが転換と速度を支配することを示唆。
• 速い領域のマスターカーブは(I sqrt(xi)/dot_gamma) K_e = A[(I/xi) − (I/xi)^*]^κのべき法則形に従い、κ ≈ 1.97、A ≈ 6、(I/xi)^* ≈ 0.02。
• 観測された閾値(I/ξ)^*と安定性解析から予測されるa/αとの理論的一致を見出し、パラメータ推定 a ≈ 0.095 および α ≈ 5 を用いて(I/ξ)^* ≈ 0.019を得る。