[論文レビュー] MC-DEM: a novel simulation scheme for modeling dense granular media

Nicolas Brodu, J.A. Dijksman+1 more|arXiv (Cornell University)|Oct 23, 2014

Granular flow and fluidized bedsEngineering参考文献 16被引用数 19

ひとこと要約

本稿では、球状粒子上に存在する複数の接触点間の相互作用を明示的にモデル化する新しい離散要素法MC-DEMを提案する。線形弾性に基づく変形を組み込むことで、密な粉体系の力の予測を向上させる。この手法は、3次元実験データの力分布およびスケーリング則と優れた一致を示し、従来のDEM（接触を独立に扱う）よりも顕著に優れており、有限要素法を用いる必要がない。

ABSTRACT

This article presents a new force model for performing quantitative simulations of dense granular materials. Interactions between multiple contacts (MC) on the same grain are explicitly taken into account. Our readily applicable method retains all the advantages of discrete element method (DEM) simulations and does not require the use of costly finite element methods. The new model closely reproduces our recent experimental measurements, including contact force distributions in full 3D, at all compression levels up to the experimental maximum limit of 13\%. Comparisons with traditional non-deformable spheres approach are provided, as well as with alternative models for interactions between multiple contacts. The success of our model compared to these alternatives demonstrates that interactions between multiple contacts on each grain must be included for dense granular packings.

研究の動機と目的

• 従来のDEMシミュレーションと密な粉体系の実験データとの間で観察される定量的不一致を是正すること。
• 標準的なDEMが無視している、1つの粒子上に存在する複数の接触の相互影響をモデル化すること。
• DEMのシンプルさを保ちつつ、密で準ジャム状態の系の精度を向上させる計算効率の良い手法を開発すること。
• 全3次元実験測定による接触力およびマクロな応答と、モデルの妥当性を検証すること。
• 粒子スケールの変形および接触間力の結合が、密な粉体媒体の正確なシミュレーションに不可欠であることを示すこと。

提案手法

• 点力が作用する球体に対する解析的線形弾性解を用いて、局所的な粒子変形を計算する。
• 他の接触点kが同一粒子上に及ぼす力によって引き起こされる接触cにおける法線方向変位δk→cを補正項として導入する。
• 有効な重なり量はδc + Σk δk→cとして更新され、これによりヘルツ則F ∝ (δc + Σk δk→c)^3/2が適用される。
• 境界条件の影響を補正するための経験的係数γを導入し、実験データへの適合性を向上させる。
• すべての変形は、長距離相互作用や有限要素法（FEM）計算を伴わず、粒子単位で計算される。
• 形状-変形関係が既知の他の粒子系（例：エマルジョンや泡）へも一般化可能である。

実験結果

リサーチクエスチョン

• RQ1 同一粒子上に複数の接触がある場合の相互作用を考慮した修正DEMモデルは、密な3次元粉体系の実験的力分布をよりよく再現できるか？
• RQ2 弾性変形を介した接触間力結合を組み込むことで、密な粉体材料のマクロおよびミクロ的応答にどのような影響を与えるか？
• RQ3 標準的なDEMにおける接触独立仮定は、ジャム状態に近い密な粉体系の真の力学的挙動をどれほど誤って捉えているか？
• RQ4 多接触効果に対する単純な解析的補正は、高コストなFEMベースのシミュレーションに代わって、高い精度を維持できるか？
• RQ5 粒子の変形および力の伝播は、密度およびひずみに伴う接触力のスケーリングにどのように影響を与えるか？

主な発見

• MC-DEMは3DXP研究から得られた3次元実験データ（全3次元力分布およびマクロな力-ひずみ応答）と優れた一致を示した。
• 実験のスケーリング則F ∝ f Z φ b と高い忠実度で一致し、標準的なDEMよりも顕著に優れた性能を示した。
• MC-DEMシミュレーションでは、荷重解除段階でわずかな曲がりが観察され、標準的なDEMよりも実験データに近いが、わずかなヒステリシスは残存している。
• MC-DEMにおける最大上部力Fは、実験のRtop値とわずかなFminオフセットの範囲内で一致しており、強い定量的整合性を示した。
• 粒子変形に起因する新規接触の形成を、標準的なDEMが欠落している重要なメカニズムとして、モデルが正しく捉えた。
• 経験的係数γは境界条件の影響を効果的に補正できており、将来的なモデルで複雑な制約を粒子レベルの統計からモデル化できる可能性を示唆している。