[論文レビュー] Internal state of granular assemblies near random close packing

Jean-Noël Roux|arXiv (Cornell University)|May 17, 2004

Granular flow and fluidized beds参考文献 4被引用数 5

ひとこと要約

本研究では、分子動力学シミュレーションを用いて、ランダム密 Packing (RCP) に近い系が密度、配位数、織構造において独立に変化することを示した。これは、RCP が一意の内部状態を意味するという考えに挑戦するものである。Φ ≈ 0.637 の同一密度でも、潤滑処理と振動による圧密プロトコルでは異なる構造が得られる。潤滑処理試料は高い配位数 (z* ≈ 5.79) と高い弾性率を示すが、振動処理試料は低い配位数 (z* ≈ 4.56) であり、実験的な音速とよく一致する。この結果、弾性率が隠れた内部状態変数を記録していることが示された。

ABSTRACT

The structure of random sphere packings in mechanical equilibrium in prescribed stress states, as studied by molecular dynamics simulations, strongly depends on the assembling procedure. Frictionless packings in the limit of low pressure are devoid of dilatancy, and consequently share the same random close packing density, but exhibit fabric anisotropy related to stress anisotropy. Efficient compaction methods can be viewed as routes to circumvent the influence of friction. Simulations designed to resemble two such procedures, lubrication and vibration (or ``tapping'') show that the resulting granular structures differ, the less dense one having, remarkably, the larger coordination number. Density, coordination number and fabric can thus vary independently. Calculations of elastic moduli and comparisons with experimental results suggest that measurable elastic properties provide information on those important internal state variables.

研究の動機と目的

• 摩擦のない系において、ランダム密 Packing (RCP) が一意の内部状態に対応すると仮定する考えに挑戦すること。
• 潤滑処理と振動処理という異なる圧密プロトコルが、RCP に近い摩擦のない球状粒子パッキングの内部構造に与える影響を調査すること。
• 弾性率が、配位数や織異方性といった隠れた内部状態変数を推定するためのプローブとして機能できるかどうかを検証すること。
• シミュレートされた弾性特性を、ガラスビーズパッキングにおける実験的音速測定と比較すること。

提案手法

• ヘーツ力の法則を用いた分子動力学シミュレーションと、接線方向力に簡略化されたカタネオ-ミンドリン-デレシエヴィッチ関係を用いる。
• 歪みパラメータ εα を基準ボックス寸法に対して定義し、ポテンシャルエネルギー W = -P∑εα の制約付き最小化により等方的圧力 P を制御する。
• 周期的境界条件を用い、4000個の球体を含む5つの独立なサンプルの平均をとることで統計的妥当性を確保する。
• 2つの圧密経路のシミュレーション：(A) 理想的潤滑（摩擦なし、低圧力圧縮）、および (C) 振動タッピング（タッピング誘発の緩和）。
• 標準的な関係式を用いて体積弾性率 (B) とせん断弾性率 (G) から音速を計算する：VP = √[(B + ⁴⁄₃G)/ρ], VS = √[G/ρ]。
• シミュレートされた音速 (VP, VS) を、Domenico (1977) および Jia & Mills (2001) の実験データと、10 kPa から 100 MPa の圧力範囲で比較する。

実験結果

リサーチクエスチョン

• RQ1 異なる圧密歴史が同じ密度で異なる内部構造を生じる場合、ランダム密 Packing は一意に定義可能だろうか？
• RQ2 RCP に近い領域において、摩擦なしの潤滑処理と振動処理の両方の圧密プロトコルで、配位数と織異方性はどのように変化するか？
• RQ3 高密度の粒状パッキングにおいて、音速 (VP, VS) が固体分率ではなく、配位数とどの程度相関するか？
• RQ4 圧力下での実験的音速測定を用いて、配位数や織といった内部状態変数を推定できるか？

主な発見

• 潤滑処理 (A) により、Φ = 0.637 ± 0.002、z* = 5.79 ± 0.007、f₀ = 1.7% のラッタラーを有するランダム密 Packing が得られ、理想 RCP と整合的である。
• 振動処理 (C) では、密度はわずかに高い (Φ = 0.635 ± 0.001) が、顕著に低い配位数 (z* = 4.56 ± 0.03) と高いラッタラー分率 (f₀ = 13.3% ± 0.6%) を示す。
• 同一密度でも、潤滑処理と振動処理の試料では配位数と織構造に顕著な差が生じており、密度、配位数、織構造が独立に変化可能であることを実証した。
• 振動処理試料 (C) の音速 (VP と VS) は、Domenico (1977) および Jia & Mills (2001) の実験データと、特に log(V) 対 log(P) の傾きにおいて、潤滑処理試料 (A) よりも良好に一致する。C 試料では VS の傾きが約 0.25、VP が約 0.22 に近く、実験と一致する。
• 潤滑処理試料の音速の圧力依存性は、有効媒体理論に一致する指数 ~1/6 のべき乗則に従うが、振動処理試料はより急で、普遍的ではないスケーリングを示す。
• 潤滑処理試料の高い弾性率は、密度ではなく高い配位数に起因しており、配位数が弾性応答の主要要因であることが示された。