[논문 리뷰] Hard discs under steady shear: comparison of Brownian dynamics simulations and mode coupling theory

Oliver Henrich, Fabian Weysser외 2명|arXiv (Cornell University)|2010. 01. 19.

Material Dynamics and Properties참고 문헌 21인용 수 47

한 줄 요약

이 연구는 정 steady shear 조건에서 이차원 이이성경 디스크에 대해 브라운 운동 시뮬레이션과 모드 상호작용 이론(MCT-ITT) 간의 처음으로 전적으로 정량적인 비교를 제시한다. 이는 MCT-ITT가 비선형 응답과 유리 상태에서 선형 응답의 부재를 정성적으로 잘 묘사하지만, 유리 상태에서의 비선형 응답을 정량적으로 과대평가함을 보여주며, 이는 주로 비틀림 하중 하에서의 구조적 회복 속도가 과소평가되기 때문이다. 이로 인해 응력과 구조적 이방성은 최대 10배 과대평가된다.

ABSTRACT

Brownian dynamics simulations of bidisperse hard discs moving in two dimensions in a given steady and homogeneous shear flow are presented close to and above the glasstransition density. The stationary structure functions and stresses of shear-melted glass are compared quantitatively to parameter-free numerical calculations of monodisperse hard discs using mode coupling theory within the integration through transients framework. Theory qualitatively explains the properties of the yielding glass but quantitatively overestimatesthe shear-driven stresses and structural anisotropies.

연구 동기 및 목표

• 비틀림된 디스크에서 모드 상호작용 이론(MCT-ITT)과 브라운 운동 시뮬레이션 간에 완전히 정량적이고 매개변수 없는 비교를 수행한다.
• MCT-ITT가 정 steady shear 조건에서 유리 상태의 정적 미세구조 및 유리성 응답을 정확히 예측하는지 테스트한다.
• 이론과 시뮬레이션 간의 격차의 근본 원인을 분석하며, 특히 구조적 이방성과 응력 과대평가에 초점을 맞춘다.
• 유체역학적 상호작용이 없는 균일하고 무질서한 비틀림 유리 상태에서 MCT-ITT 프레임워크를 검증한다.
• 유리 상태에서 MCT-ITT가 예측하는 비선형 응답과 선형 응답 부재 영역이 실제로 시뮬레이션에서 성립하는지 평가한다.

제안 방법

• Lees-Edwards 경계 조건과 열역학적 조절기를 사용하여 2차원 이이성경 디스크 시스템의 브라운 운동 시뮬레이션을 수행한다. 이는 균일한 비틀림 유량을 유지한다.
• 반사 이벤트 기반의 반정적 알고리즘을 사용하여 입자 충돌을 시뮬레이션하고, 충돌 이벤트 동안 운동량 이동을 통해 시간 평균 응력을 계산한다.
• 등시간 구조 인자 Sq(˙γ)와 그 비틀림 유도 왜곡 (Sq(˙γ) − Sq)/Sq를 각도 평균 쌍 상관관계를 통해 계산한다.
• MCT-ITT 프레임워크를 사용하여 비틀림 하에서 단일성 디스크의 정적 응력 및 구조 인자 응답을 해석한다.
• 동일한 비틀림 속도와 밀도에서 시뮬레이션 결과를 직접 MCT-ITT 예측과 비교하며, 일관된 정규화 및 평균화 절차를 사용한다.
• 특히 중간 시간 영역에서 중요한 임시 밀도 상관 함수를 특성화한다.

실험 결과

연구 질문

• RQ1 MCT-ITT는 비틀림된 유리 상태에서의 응력과 정규 응력에 대해 정량적으로 정확하게 예측하는가?
• RQ2 MCT-ITT는 2차원 디스크 유리의 미세구조에서 비틀림에 의한 구조적 이방성을 얼마나 정확히 묘사하는가?
• RQ3 왜 MCT-ITT는 시뮬레이션보다 구조적 왜곡과 응력의 크기를 과대평가하는가?
• RQ4 MCT-ITT가 예측하는 유리 상태에서 선형 응답 영역이 없는 현상은 시뮬레이션에서도 실제로 관찰되는가?
• RQ5 MCT-ITT의 근사치, 예를 들어 평균장 분리 및 유체역학적 상호작용 무시 등으로 인해 발생하는 격차는 어느 정도인가?

주요 결과

• MCT-ITT는 정적 성질의 비해석적 행동과 비틀림 유리 상태에서 선형 응답 영역의 부재를 정성적으로 잘 재현한다.
• 이론은 비틀림에 의한 구조적 이방성을 최대 10배 과대평가한다. 특히 압축축과 확장축 방향에서 두드러진다.
• ϕ = 0.79 및 Pe₀ = 2×10⁻⁴ 조건에서, 상대적 왜곡 (Sq(˙γ) − Sq)/Sq는 유리 상태에서 평형 상태에서의 20% 이상의 이격을 보이며, 이는 이론이 예측하는 비선형 응답과 일치한다.
• 시뮬레이션 결과는 주 피크의 저q 영역에서 등방성의 왜곡을 보이지만, MCT-ITT는 압축축 방향(qx = −qy)에서 더 강한 이방성 감소를 예측한다.
• MCT-ITT의 과대평가 원인은 임시 상관 함수가 너무 느리기 때문이며, 이는 구조적 회복 이전에 과도한 변형이 축적되기 때문이다.
• 이이성경 시뮬레이션과 단일성 이론 간의 평형 상태 구조 인자 Sq의 격차는 오차의 주요 원인이지만, 그 영향은 제한적이다.