[論文レビュー] Sound wave velocities in dry and lubricated granular packings packings: numerical simulations and experiments

Ivana Agnolin, Jean-Noël Roux+3 more|arXiv (Cornell University)|Jan 20, 2009

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

ひとこと要約

本研究では、数値シミュレーションと実験を組み合わせて、密度が低いにもかかわらず潤滑された粒子パッケージングは、連絡数が高いために乾燥パッケージングよりも顕著に高い音波速度を示すことを示した。シミュレーションでは、準備方法、特に潤滑処理が密度よりも微視的構造に与える影響が顕著であり、連絡数が剛性および波動伝播に影響を与える主因であることが明らかになった。これは、密度のみが密度の高い粒子系材料の弾性特性を決定すると仮定する考え方に挑戦するものである。

ABSTRACT

Numerical simulations are used to investigate the origins of the different wave velocities measured in dense granular samples assembled with different methods. Glass bead packings are prepared in the lab either by pouring and vibrating the dry material in a container, or by mixing with a very small amount of a viscous lubricant. Lubricated samples, although less dense, exhibit significantly higher wave velocities for confining pressures in the 100 kPa range. Numerical predictions for elastic moduli agree much better with experimental results when the computational preparation of the samples mimics the laboratory one, albeit in a simplified manner. A plausible explanation to the laboratory observations is that the coordination number, which influences the material stiffness more than its density, is notably higher in lubricated packings.

研究の動機と目的

• 潤滑された粒子パッケージングが密度が低いにもかかわらず、乾燥パッケージングよりも高い弾性波速度を示す理由を理解すること。
• 乾燥パッキングと潤滑混合という異なるサンプル準備法が、連絡数や剛性といった微視的構造的性質に与える影響を調査すること。
• 数値シミュレーションが、粒子パッケージングにおける実験的波速度測定値を正確に再現できるかどうかを評価すること。
• 密度と連絡数のどちらが、密な粒子系材料における弾性波速度に支配的要因であるかを特定すること。
• 簡略化された数値モデルが、実験室での準備プロトコルが粒子の微視的構造および機械的応答に与える影響を的確に捉えられるかどうかを評価すること。

提案手法

• 実験室での準備方法の簡略化表現を用いた離散粒子動力学による数値シミュレーションを用いて、『振動させた』乾燥（C型）および『潤滑された』（A型およびB型）アセンブリをモデル化した。
• 定圧下での力学的平衡状態をシミュレートした後、接触ネットワークから剛性行列を計算し、弾性モジュラスを導出する。
• 剛性行列から有効体積率（B）および剪断率（G）を計算するために、Voigt-Reuss-Hill平均化法を用いた。
• 標準的な弾性理論式を用いて体積弾性率（VP）および横波速度（VS）を導出する：VP = √[(B + 4G/3)/ρ]、VS = √[G/ρ]。
• 乾燥（E1）および潤滑（E2）されたガラスビードパッケージングの超音波測定による実験データと、異なる定圧（30–1000 kPa）下でのシミュレート波速度を比較した。
• シミュレートされたサンプルの連絡数および接触ネットワーク構造を分析し、微視的構造とマクロな弾性応答の相関を検証した。

実験結果

リサーチクエスチョン

• RQ1 潤滑された粒子パッケージングは、密度が低いにもかかわらず、なぜ乾燥パッケージングよりも高い音波速度を示すのか？
• RQ2 特に潤滑処理を含むサンプル準備法が、密な粒子パッケージングにおける連絡数および弾性モジュラスにどの程度の影響を与えるのか？
• RQ3 実験室の準備手順を簡略化した数値シミュレーションが、粒子系材料における実験的波速度トレンドを再現できるか？
• RQ4 連絡数と密度のどちらが、粒子パッケージングにおける弾性波速度により大きな影響を与えるのか？
• RQ5 異なる定圧下で、数値的にシミュレートされたパッケージングの機械的応答は、実験結果とどの程度一致するのか？

主な発見

• 潤滑パッケージング（E2）は、固体分率が0.62（乾燥パッケージングの約0.64）であるにもかかわらず、縦波および横波速度が10–20%高い。
• 数値シミュレーションでは、C型（振動乾燥パッケージングを模倣）のサンプルが、乾燥パッケージング（E1）の実験的波速度をA型（完全潤滑）やB型（不完全潤滑）のシミュレーションよりもはるかによく再現した。
• 潤滑パッケージングの高い波速度は、剛性を密度よりも顕著に高める高い連絡数に起因する。
• シミュレーションでは、密な粒子パッケージングにおいて、連絡数がパッケージング密度よりも弾性波速度に顕著に大きな影響を与えることが明らかになった。
• 本研究では、数値モデルが最終的な密度だけでなく、準備プロトコルを正確に再現する必要があることが示された。なぜなら、連絡数は微視的構造の独立変数であるためである。
• 低定圧（約100 kPa）では、高・低連絡数パッケージングの差が重要であるが、高圧（>1 MPa）では、シミュレーションタイプ間の差が小さくなり、実験トレンドと一致するようになった。