[論文レビュー] Contact of Viscoelastic Spheres

Thomas Schwager, Thorsten Poeschel|arXiv (Cornell University)|Mar 22, 2002

Adhesion, Friction, and Surface Interactions被引用数 6

ひとこと要約

本稿では、粘弾性球体衝突における速度に依存する法線復元係数を、粘弾性接触方程式の摂動級数解法を用いて導出する。解析的に 1−ϵ ∼ g^{1/5} が示され、衝突速度の増加に伴う復元係数の低下が実験的に観察される現象を説明し、粗粒ガス力学における定数復元係数モデルの不備を解消する。

ABSTRACT

In a recent paper an implicit equation for contacting viscoelastic spheres was derived (Brilliantov et al. PRE 53, 5382 (1996)). Integrating this equation it can be shown that the coefficient of normal restitution depends on the impact velocity as 1-ε\sim g^{1/5}.

研究の動機と目的

• 粗粒ガスのモデル化においてエネルギー損失によってクラスタリング不安定性を示す、定数復元係数の不一致を解消すること。
• 粘弾性接触力学に基づき、粘弾性球体の物理的に正確な速度に依存する法線復元係数を導出すること。
• 復元係数および衝突時間の解析的解を提供し、数値的・実験的結果と一致させること。
• 古典的ヘルツ接触モデルに、減衰パラメータの摂動級数展開を用いて粘弾性減衰効果を拡張すること。

提案手法

• 弾性および粘性項を含む、粘弾性球体接触の運動方程式を定式化し、圧縮量 ξ および有効質量・半径を用いる。
• ρ および衝突速度 g を用いて時間、速度、圧縮量をスケーリングし、減衰パラメータ α を持つ次元なし方程式 (7) を導出する。
• η(Θ) に対して √Θ の級数展開を用い、式 (11) を解く。x(Θ) = v^{4/5}Θ(1+η(Θ)) かつ η(0)=0 を仮定する。
• 摂動理論を α² 項まで適用し、最大圧縮量および衝突時間を計算する。境界条件の一致には逆衝突法を用いる。
• 中間衝突で前向きおよび逆向きの衝突解を一致させ、v′ を v および α の関数として導出し、復元係数 ε = v′/v を得る。
• ε および衝突時間 t_c の漸近的表現を導出し、g^{1/5} への明示的依存性を示す。

実験結果

リサーチクエスチョン

• RQ1 粘弾性球体衝突において、法線復元係数 ε は衝突速度 g にどのように依存するか？
• RQ2 定数 ε の仮定が粗粒ガスの挙動を記述できないのはなぜか？ その代替モデルは何か？
• RQ3 粘弾性接触方程式から導出された復元係数 ε(g) の解析的形は何か？
• RQ4 粘弾性減衰が存在する場合、衝突時間 t_c は衝突速度 g に対してどのようにスケーリングするか？
• RQ5 数値積分を伴わずに、摂動級数解法が速度に依存する復元係数の挙動を正確に捉えられるか？

主な発見

• 法線復元係数は衝突速度の増加に伴い減少し、1−ϵ ∼ g^{1/5} の関係に従う。これは実験データおよび数値シミュレーションと一致する。
• 解析的解として ε = 1 − (3/2)C₁Aρ^{2/5}g^{1/5} + O(g^{2/5}) が得られ、材料定数および衝突速度への明示的依存性が示された。
• 衝突時間 t_c は t_c ∼ g^{−1/5} に比例し、摂動展開において g^{1/5} の補正項が存在する。
• 最大圧縮は、減衰がある（粘弾性）衝突では、ない（弾性）衝突よりも早く発生し、衝突時間も長くなる。これは図1のスケーリングされた力学的挙動で確認された。
• α² 項までの摂動級数解法は、数値結果を正確に再現し、粘弾性衝突におけるエネルギー散逸の本質的物理を捉えている。
• 逆衝突法により、中間衝突で前向きおよび後向きの解を一致させることができ、Θ_c⁰ を超える直接積分を伴わずに ε(v) を導出可能である。