[논문 리뷰] Contact of Viscoelastic Spheres

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

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}의 관계를 따르며, 실험 데이터와 수치 시뮬레이션과 일치한다.
• ε(g)에 대한 분석적 해는 ε = 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)를 도출할 수 있다.