Percolation of a cohesive fine particle in a static bed

技術的サマリー：静的な粒子層における凝集性微粒子の浸透

問題提起
微粒子（ファイン）が、より大きな粒子からなる静的な粒子層を通過する浸透現象は、製薬や化学プロセスなどの産業処理、および土砂輸送や地下水汚染などの地球物理的流れにおいて、基本的なプロセスである。非凝集系では、微粒子が通過するか捕捉されるかは、主に幾何学的制約、すなわち大粒子と微粒子の直径比（$R = d_l/d_f$）によって支配される。幾何学的捕捉閾値は $R_t \approx 6.464$ に存在し、$R > R_t$ の場合、微粒子は理論的には大粒子が形成する細孔の喉部を自由に通過すると期待される。しかし、現実のシステムでは、接触ダイナミクスを変化させる凝集力（ファンデルワールス力、静電気力、液体ブリッジ）がしばしば関与する。過去の研究では非凝集性の浸透や凝集性流れのレオロジーが扱われてきたが、幾何学的制約のみが通過を許容する「自由篩い」領域において、凝集性が浸透にどのように影響するかという点については理解のギャップが存在する。本研究は、$R=7$（幾何学的閾値以上）の条件下で、大粒子からなる静的な粒子層を通過する凝集性微粒子の浸透を調査し、微粒子が粒子層粒子とは相互作用するが互いには相互作用しない単一粒子限界に焦点を当てている。

手法
著者は、オープンソースソフトウェア LAMMPS を用いた離散要素法（DEM）シミュレーションを採用している。シミュレーションでは、大きなサイズ比に適した改良された速度ベルレ法が利用されている。

• システム構成：単分散の大粒子（$d_l = 4$ mm）からなる静的な粒子層が、充填率 $\phi = 0.6$ で充填されている。微粒子（$d_f \approx 0.57$ mm、$R=7$）は粒子層の上方から静止状態から放出され、重力下で落下する。
• 接触モデル：法線接触力はヘルツの弾性モデルに従う。非弾性性は、一定の復元係数（$e_n$）を介してモデル化される。接線力は、すべり摩擦（$\mu$）と転がり抵抗を伴うミンディルモデルを用いる。
• 凝集性：凝集相互作用は、主に接触中に一定の引力を仮定し、ヘルツの反発力を伴う Derjaguin-Muller-Toporov (DMT) モデルを用いてモデル化される。凝集性の強さは、剥離力と微粒子の重量の比として定義されるボンド数（$Bo = F_c / W$）によって定量化される。比較のために、Johnson-Kendall-Roberts (JKR) モデルを用いた限られた数のシミュレーションも実施された。
• パラメータ：本研究では、ボンド数（$Bo$）、すべり摩擦係数（$\mu$）、復元係数（$e_n$）、ヤング率（$E$）を体系的に変化させた。統計的な精度を確保するため、$10^4$ 個の独立した微粒子のダイナミクスを追跡した。

主要な結果

1. 非幾何学的捕捉：$R=7$（幾何学的に自由な通過を許容）であっても、十分に大きな凝集性と摩擦は粒子の捕捉をもたらす。$Bo$ が増加するにつれて、特徴的な浸透距離は減少し、粒子はより浅い深さで捕捉されるようになる。
2. 捕捉分布：捕捉された粒子の垂直分布は指数関数的減衰に従う。特徴的な捕捉距離（$\lambda$）は、浸透挙動の指標として機能する。逆捕捉距離（$d_l/\lambda$）は $Bo$ とともに単調に増加する。
3. 摩擦と凝集性の役割：
   • 低 $Bo$ において、すべり摩擦は捕捉に重要な役割を果たす。
   • 高 $Bo$ において、凝集性が支配的な捕捉メカニズムとなる。
   • 逆捕捉距離と摩擦の関係において「折れ曲がり」が観測される。この折れ曲がりは $Bo$ に伴ってシフトし、摩擦が支配的な領域と凝集性が支配的な領域との間の遷移に対応する。
   • 粒子の剛性（$E$）は捕捉に影響を与える：より柔らかい粒子（低い $E$）は、特に高い $Bo$ において捕捉されやすくなる。
4. 衝突メカニズム：捕捉プロセスは、2 つの連続するステップに分解される。
   • 付着条件：衝突時、粒子は跳ね返るか付着する。付着は、衝突速度が臨界速度（$v_{cr}$）以下の場合に発生し、これは $Bo$、$e_n$、および$E$ に依存する。
   • 接触後運動：粒子が付着した場合、それは滑って剥離するか、捕捉されたまま残る。非跳ね返り粒子が静止したまま残るかどうかは、力平衡に基づく静的平衡条件によって決定される。この条件は、再スケーリングされた摩擦係数 $\tilde{\mu} = \mu\sqrt{4Bo^2 - 1}$ に依存する。$\tilde{\mu} > 1$ の場合、粒子は静止したまま残る可能性が高く、$\tilde{\mu} < 1$ の場合、滑って剥離する可能性が高い。
5. 予測モデル：著者は、微視的ダイナミクスを巨視的捕捉距離に結びつける衝突モデルを開発した。逆捕捉距離は $1/\lambda = -k \ln(1 - P_t)$ によって予測される。ここで、$k$ は衝突頻度、$P_t$ は衝突ごとの捕捉確率（衝突速度と $v_{cr}$ の比によって決定される）である。このモデルは、衝突支配領域（高 $Bo$）における捕捉距離を正確に予測する。
6. 二重接触：摩擦または凝集性が低い領域では、粒子はしばしば滑り、二重接触（2 つの大粒子に同時に接触する状態）を形成する。この配置は、単一接触条件では剥離が示唆される場合でも捕捉をもたらすか、あるいは逆に、2 番目の接触からの跳ね返りによって剥離を可能にする。

意義と主張
本論文は、凝集性微粒子の浸透が幾何学的な到達可能性のみによって決定されるわけではないことを確立している。むしろ、凝集性、復元性、摩擦、および衝突速度の相互作用といった粒子スケールの相互作用ダイナミクスが、幾何学的な期待を上書きし、細孔の喉部が幾何学的に通過に十分な広さであっても捕捉を引き起こしうる。
著者は、凝集性粒子の流動混合物におけるより複雑な粒度分離の問題に取り組むための前提条件となる、静的な粒子層におけるこれらのメカニズムを理解するための基礎的枠組みを提供すると主張している。付着およびその後の運動を支配する特定の条件（臨界速度、静的平衡基準）を特定することにより、本研究は、基本的な材料特性に基づいた捕捉距離の予測能力を提供する。結果は、凝集性システムにおいて、微粒子の「運命」は、単に粒子層の静的幾何学ではなく、衝突、付着、接触後運動という一連の局所的相互作用によって支配されることを浮き彫りにしている。