[論文レビュー] Geometry of Frictionless and Frictional Sphere Packings

Leo Silbert, Deniz Ertaş+3 more|OpenSIUC (Southern Illinois University Carbondale)|Nov 8, 2001

Pickering emulsions and particle stabilization被引用数 39

ひとこと要約

本研究では分子動力学シミュレーションを用いて摩擦なしおよび摩擦ありの球体パッケージングを調査し、摩擦なしの硬い球体パッケージングは普遍的に静定的（z=6）である一方、摩擦ありパッケージングは摩擦が増加するにつれて協調数が6から4に滑らかに減少する過静定的挙動を示すことが明らかになった。これは静定性仮説に反する。結果は摩擦係数、初期密度、および系の歴史に強く依存しており、摩擦ありパッケージングにおける力のネットワークが幾何学的要因のみによって一意に決定されないことを示している。

ABSTRACT

We study static packings of frictionless and frictional spheres in three dimensions, obtained via molecular dynamics simulations, in which we vary particle hardness, friction coefficient, and coefficient of restitution. Although frictionless packings of hard-spheres are always isostatic (with six contacts) regardless of construction history and restitution coefficient, frictional packings achieve a multitude of hyperstatic packings that depend on system parameters and construction history. Instead of immediately dropping to four, the coordination number reduces smoothly from $z=6$ as the friction coefficient $μ$ between two particles is increased.

研究の動機と目的

• 摩擦なしおよび摩擦ありの球体パッケージングが一般的なパッケージング条件下で静定的協調数に達するかどうかを調査すること。
• 粒子の硬さ、摩擦係数、回復係数、および初期条件がパッケージングの最終的な幾何学的形状および機械的安定性に与える影響を特定すること。
• 沈殿によって形成された三次元パッケージングにおける摩擦あり球体の静定性仮説（z=4）の妥当性を検証すること。
• 構築歴およびエネルギー散逸が最終的なパッケージング構造および力のネットワークに与える影響を検討すること。
• 弾性で有限剛性を有するシミュレーションから外挿することで、硬球極限（k_n → ∞）を意味的に到達可能とすることができるかどうかを明確にすること。

提案手法

• 3次元周期的ボックス内にN=20,000個の単分散球体を用いた分子動力学シミュレーション。底面は粗い境界、上面は開放されている。
• 球体は法線方向および接線方向に線形ばね-ダンパー接触則に従い、パラメータとして法線剛性k_n、摩擦係数μ、および回復係数εを用いる。
• 初期状態は初期パッケージング密度φ^iを0.02から0.3の範囲で変化させた上で、重力による沈殿を経て静的パッケージングを形成する。
• 時間積分に速度-ベルレ法を用い、運動エネルギーがしきい値を下回るまで系を進化させ、静的平衡状態に達する。
• パrameter空間全体で協調数z、最終パッケージング密度φ^f、および接触力比F_t/F_nを測定・分析する。
• 接線力と法線力の比ζ = F_t / (μ F_n)を用い、接触が摩擦限界以下にあるかどうかを評価する。

実験結果

リサーチクエスチョン

• RQ1 摩擦あり球体パッケージングにおける協調数zが摩擦係数μが増加するにつれて静定値4に近づくか？
• RQ2 最終パッケージング密度φ^fおよび協調数zは初期パッケージング密度φ^iおよび回復係数εにどのように依存するか？
• RQ3 構築歴および系パラメータ（k_n、μ、ε）が摩擦ありパッケージングの最終的幾何学的形状および力のネットワークに及ぼす影響の程度はどの程度か？
• RQ4 過静定的摩擦ありパッケージングにおける力のネットワークは、粒子の幾何学的形状と外部荷重によって一意に決定されるのか、それとも歴史依存性を示すのか？
• RQ5 有限剛性シミュレーションからの外挿によって、硬球極限（k_n → ∞）を意味的に到達可能とすることができるか？

主な発見

• 摩擦なしの硬い球体パッケージングは、構築歴や回復係数に依存せず、常に静定的（z=6）である。
• 摩擦ありパッケージングは過静定的であり、摩擦係数μが増加するにつれて協調数zが6から4に滑らかに減少するが、これは静定性仮説に反する。
• μ ≥ 1の場合、接触力比F_t / F_nの分布がζ = F_t / (μ F_n)に対して重ね合わさるため、有効に無限摩擦挙動を示している。
• 最終パッケージング密度φ^fは初期パッケージング密度φ^iが小さくなるにつれて増加し、経験的フィット式φ^f = 0.5778 + 0.0567 exp(−4.3φ^i)に従う。
• φ^i → φ^fの極限において協調数zがz=4に外挿されないため、過静定パッケージングが静定パッケージングと統計的に同等ではないことが示された。
• 摩擦ありパッケージングにおける力のネットワークは、幾何学的形状と荷重によって一意に決定されないため、過静定領域では問題が不適切に定式化されている。歴史および準備経路が結果に顕著に影響を及ぼす。