[論文レビュー] Asymptotic Energy Decay in Inelastic Fluids

J.A.G. Orza, Ricardo Brito+1 more|arXiv (Cornell University)|Feb 24, 2000

Granular flow and fluidized beds被引用数 3

ひとこと要約

本研究では、N=50,000粒子のシミュレーションを用いて、自由に進化する非弾性硬球系におけるエネルギー減衰の2つの理論的予測を検証する。モード結合理論（BritoとErnstによるもの）は、E∼Aτ⁻ⁱ/²というτに依存するべきべき乗則の減衰を予測するが、Ben-Naimらの提唱する、非弾性度に依存しない普遍的なt⁻ⁱ/²の減衰は、特に中程度および高密度領域では明確に支持されない。この領域では、長時間経過しても非弾性度の効果が依然として顕著に現れる。

ABSTRACT

Replacement of cond-mat/0002383. Re-analysis of computer simulations, new figures and modification of conclusions.

研究の動機と目的

• 非弾性硬球系における長時間エネルギー減衰の2つの理論的予測の妥当性を検証すること。
• Ben-Naim らが提唱するように、漸近的に大きな時間においてエネルギー減衰が回復係数に依存しなくなるかどうかを特定すること。
• BritoとErnstが提唱するモード結合理論の適用可能性を評価すること。この理論は、Aがrおよびφに依存するE∼Aτ⁻ⁱ/²を予測する。
• 系のサイズおよび境界条件が、有限系における真の漸近的挙動を隠す役割を果たすメカニズムを調査すること。
• 短時間におけるHaffの法則と、高密度非弾性流体における長時間エネルギー減衰ダイナミクスとの間に生じる矛盾を明確にすること。

提案手法

• 2次元非弾性硬ディスク系における大規模分子動力学的シミュレーション（N=50,000）を実施する。
• モード結合理論に従い、τ(t)を非線形時間変数として用い、データの収束を図る。
• log E(τ) 対 τ のプロットを作成し、d=2における予測されたべき乗則の減衰 E∼Aτ⁻ⁱ/² を検証する。
• 逆エネルギー 1/E(t) 対 t のプロットを作成し、スティッキー・ガスの予想する漸近的非依存性を検証する。
• 異なる充填率 φ=0.05, 0.11, 0.245, 0.4 における結果を比較し、密度依存性を評価する。
• 系サイズ効果を分析し、特に L∼Lv(τ)∝√τ に達した際に周期的境界条件のアービファクトが支配的になる遷移を特定する。

実験結果

リサーチクエスチョン

• RQ1 非弾性流体におけるエネルギー減衰は、漸近的に大きな時間において、回復係数に依存しない普遍的なべき乗則 t⁻ⁱ/² に従うか？
• RQ2 BritoとErnstのモード結合理論は、τ(t) を用いた長時間エネルギー減衰をどの程度正確に予測できるか？
• RQ3 密度ゆらぎおよびクラスタリングは、非弾性流体における線形応答理論の有効性にどのように影響を与えるか？
• RQ4 有限サイズ効果がどの時間スケールで支配的になり、真の漸近的挙動を隠すようになるか？
• RQ5 系サイズは漸近的減衰則に影響を与えるか？その場合、どのような影響を及ぼすか？

主な発見

• BritoとErnstのモード結合理論は、φ≈0.05, 0.11, 0.245, 0.4 に対応する r₀≈0.6, 0.75, 0.8, 0.85 を中心とする幅Δr≈0.1 のrウィンドウ内で、E(τ) に対して定量的に正確な予測を提供する。
• 低密度（φ=0.05）では、極端な非弾性度において理論からのずれが生じ、相対的密度ゆらぎが著しく大きいことが原因で理論の破綻を示す。
• 逆エネルギー 1/E(t) 対 t のプロットにおいて、φ=0.4 であっても、異なるr値間で共通の線形挙動への収束が観察されず、長時間においても非弾性度依存性が持続していることが示唆される。
• 高密度（φ=0.4）では、1/E(t) の曲線が後期にかけて発散する傾向を示し、rに依存しない漸近的独立性が明確に示されない。
• 小系では、t→∞ においてエネルギーが E(t)∼C(N,r)/t² と減衰するが、CはNおよびrに依存する。これはHaffの法則に反し、有限サイズ効果を示している。
• 周期的境界条件の支配的になる遷移は、系サイズLがパターンサイズ Lv(τ)∝√τ と同等に達した際に発生し、この時点で小系では真の漸近的挙動に到達できなくなる。