[論文レビュー] A minimal model for slow dynamics: Compaction of granular media under vibration or shear

Stefan Luding, Maxime Nicolas+1 more|arXiv (Cornell University)|Mar 10, 2000

Granular flow and fluidized beds参考文献 20被引用数 8

ひとこと要約

本論文は、振動またはせん断による粒子凝集の最小モデルを提案する。このモデルは、不均一なランダムエネルギー障壁の上をランダムウォークするもので、温度に類似した駆動力がエネルギー極小値の探索を制御する。モデルは、タッピングおよび周期的せん断実験で観察される遅い対数的凝集ダイナミクスとヒステリシス行動を定性的に再現し、巨視的密度の変化を不規則なポテンシャルにおける確率的ダイナミクスと結びつける。

ABSTRACT

Based on experiments of the compaction of granular materials under periodic shear of a packing of glass beads, a minimal model for the dynamics of the packing density as a function of time is proposed. First, a random ``energy landscape'' is created by a random walk (RW) in energy. Second, an ensemble of RWs is performed for various temperatures in different temperature-time sequences. We identify the minimum (mean) of the energy landscape with the maximum (random) density. The temperature scaled by the step-size of the energy landscape determines the dynamics of the system and can be regarded as the agitation or shear amplitude. The model reproduces qualitatively both tapping and shear experiments. {\bf Key-words}: Compaction, crystallization, Sinai-diffusion, random walks, quenched disorder

研究の動機と目的

• 周期的なタッピングまたは周期的せん断の下で観察される遅い対数的凝集ダイナミクスを説明すること。
• 詳細な粒子レベルのシミュレーションを伴わずに凝集の主要な特徴を捉える最小モデルを構築すること。
• 外部からの駆動（タッピング／せん断）と熱的温度との類似性が、不規則なポテンシャルにおけるエネルギー極小値への緩和を記述するのに有効であるかを検討すること。
• 凝集におけるヒステリシス行動が、不規則なポテンシャル障壁における非平衡ダイナミクスに起因するかを調査すること。

提案手法

• 粒子凝集体の複雑なポテンシャルエネルギー表面を表すために、ランダムウォークを用いて不均一な1次元エネルギー障壁を生成する。
• 局所的なエネルギー差と温度パラメータ T に依存する遷移確率を持つ、バイアス付きランダムウォークを多数のウォークャーが行う。
• 温度 T はエネルギーステップサイズにスケーリングされ、せん断またはタッピング振幅に類似した制御パラメータとして機能する。
• 系の密度はエネルギー障壁の最小値に対応させられ、平均密度は時間とともに到達するグローバル最小値に対応する。
• シミュレーションでは、確率分布の時間発展を離散的マスター方程式で記述し、時間はモンテカルロステップで測定する。
• 周期的な温度変化（T を増加してから減少）を用いて、振幅が変化する実験的手順を模倣する。

実験結果

リサーチクエスチョン

• RQ1 不均一なランダムエネルギー障壁における単純なランダムウォークが、粒子実験で観察される遅い対数的凝集ダイナミクスを再現できるか？
• RQ2 本モデルは、タッピングまたはせん断実験で振幅を増加させてから減少させた際のヒステリシス行動を捉えられるか？
• RQ3 外部からの駆動（タッピング／せん断）と温度との類似性が、不規則系における非平衡緩和を記述するのに妥当か？
• RQ4 系の緩和時間は温度パラメータ T に対してどのようにスケーリングされ、拡散的から準拡散的行動への遷移は何かの要因によって決定されるか？
• RQ5 急激なせん断振幅の変化（それに類する温度の急変）に伴い観察される体積分率の急速で不可逆的なジャンプをモデルは説明できるか？

主な発見

• モデルは、タッピングおよび周期的せん断実験の両方で観察されるように、時間とともにゆっくりと対数的に密度が増加するダイナミクスを再現する。
• ヒステリシス行動がモデルで観察される：温度を減少させる際の密度は、増加させる際の密度より低く、実験で観察される不可逆的凝集と整合する。
• 温度の急激な変化（せん断振幅の変化に類似）により、直ちに逆方向の密度ジャンプが発生する。これは20サイクル以内に発生し、初期状態に依存しない。
• 低温ではランダムウォークが局所的極小値に捕らわれ、準拡散的行動を示すが、高温では拡散的運動が支配的になる。
• 拡散定数は温度が低下するにつれて減少し、T=0 でゼロに近づくことから、深い極小値での局在化が示唆される。
• 温度を低下させると、均一なランダムウォーク行動からシナ型拡散に移行する。大規模系では、グローバル最小値に到達するには長時間が必要となる。