Skip to main content
QUICK REVIEW

[论文解读] Hard discs under steady shear: comparison of Brownian dynamics simulations and mode coupling theory

Oliver Henrich, Fabian Weysser|arXiv (Cornell University)|Jan 19, 2010
Material Dynamics and Properties参考文献 21被引用 47
一句话总结

本研究首次对二维双分散硬盘在稳态剪切下的布朗运动模拟与模式耦合理论(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是否能准确预测剪切作用下剪切诱发玻璃态的稳态微观结构与流变响应。
  • 探究理论与模拟之间差异的根源,特别是关于结构各向异性和应力高估的问题。
  • 在无流体动力学相互作用的均匀、无序剪切熔融玻璃态背景下,验证MCT-ITT框架的有效性。
  • 评估MCT-ITT在玻璃态中预测的非线性响应及无线性响应区间的有效性。

提出的方法

  • 采用Lees-Edwards边界条件和恒温器,对二维双分散硬盘系统进行布朗运动模拟,以维持均匀剪切流。
  • 使用半事件驱动算法模拟粒子碰撞,并通过碰撞过程中的动量传递计算时间平均应力。
  • 利用角度平均的对径相关函数,计算等时结构因子Sq(˙γ)及其剪切诱导畸变(Sq(˙γ) − Sq)/Sq。
  • 应用MCT-ITT框架求解单分散硬盘在剪切作用下的稳态应力与结构因子响应。
  • 在相同的剪切速率和密度下,通过一致的归一化与平均程序,将模拟结果与MCT-ITT预测值直接比较。
  • 表征MCT-ITT计算中至关重要的瞬态密度关联函数,尤其是在中间时间区间。

实验结果

研究问题

  • RQ1MCT-ITT能否定量预测剪切诱发玻璃态中的剪切应力与法向应力?
  • RQ2MCT-ITT在描述二维硬盘玻璃微观结构中剪切诱导的结构各向异性方面有多准确?
  • RQ3为何MCT-ITT相比模拟高估了结构畸变与应力的幅值?
  • RQ4如MCT-ITT所预测,在玻璃态中不存在线性响应区间的结论,在模拟中是否成立?
  • RQ5MCT-ITT中的近似(如平均场解耦与忽略流体动力学相互作用)在多大程度上导致了差异?

主要发现

  • MCT-ITT定性地再现了稳态性质的非解析行为以及剪切熔融玻璃态中无线性响应区间的特征。
  • 该理论将剪切诱导的结构各向异性高估了高达10倍,尤其在压缩轴与拉伸轴方向更为显著。
  • 在ϕ = 0.79与Pe₀ = 2×10⁻⁴条件下,相对畸变(Sq(˙γ) − Sq)/Sq在玻璃态中表现出与平衡态20%的偏离,与理论预测的非线性响应一致。
  • 模拟显示在主峰低q端翼区呈现各向同性偏离,而MCT-ITT则预测在压缩轴方向(qx = −qy)出现更强的各向异性抑制。
  • MCT-ITT的高估主要源于瞬态相关函数过慢,导致在结构松弛发生前已积累过多应变。
  • 双分散模拟与单分散理论之间平衡结构因子Sq的差异可能对偏差有所贡献,但并非误差的主要来源。

更好的研究,从现在开始

从阅读论文到最终审阅,大幅缩短您的研究时间。

无需绑定信用卡

本解读由 AI 生成,并经人工编辑审核。