Tunable shear thickening in active non-Brownian suspensions

该研究通过粒子模拟发现，在非布朗摩擦颗粒活性悬浮液中，粒子自推进产生的各向同性动力学能竞争并抑制流动诱导的摩擦接触，从而在高应力下实现粘度降低的“去剪切增稠”效应，且该可调控流变行为遵循统一的标度框架。

要点

这篇论文讲述了一个非常有趣的现象：如何让一种“越搅越稠”的液体，通过加入“自己会动”的微粒，重新变回“稀汤”。

为了让你轻松理解，我们可以把这篇论文的研究对象想象成一锅正在煮的、非常浓稠的玉米糊（或者像玉米淀粉和水混合的“非牛顿流体”）。

1. 什么是“剪切增稠”？（越搅越硬）

想象一下，你手里拿着一杯这种特殊的玉米糊。

• 慢慢搅动时：它像水一样流动，很顺滑。
• 用力快速搅动时：它突然变得像石头一样硬，你的勺子甚至搅不动了。

这就是论文里说的**“剪切增稠”（Shear Thickening）。
为什么会这样？ 论文解释说，当搅拌速度很快时，里面的颗粒（玉米粒）来不及互相滑过，它们会互相“卡”在一起，形成像脚手架一样的“力链网络”**（Force Chains）。这些颗粒互相摩擦、锁死，导致液体瞬间变硬。这在工业上很麻烦，比如管道里的泥浆突然变硬堵住了管子。

2. 怎么解决？（让颗粒“自己跑起来”）

通常，人们想解决这个问题，会往里面加肥皂水（表面活性剂）让颗粒变滑，或者用超声波去“震”它们。但这篇论文提出了一个更酷、更“主动”的方法：

给颗粒装上“小马达”！

研究人员在模拟实验中，让悬浮液里的颗粒变成了**“活性颗粒”**（Active Particles）。

• 比喻：想象这锅玉米糊里，每一粒玉米都变成了一只微小的、会自己乱跑的蚂蚁。
• 现象：当你用力搅拌（施加剪切力）时，这些“蚂蚁”因为自己也在到处乱窜（自驱动），它们就没有时间乖乖地排好队、互相卡死形成“脚手架”了。
• 结果：原本应该变硬的液体，因为颗粒的“捣乱”和“乱跑”，反而变稀了！论文把这种现象称为**“去增稠”（Dethickening）**。

3. 核心发现：如何“调音”？

研究的关键在于**“调音”**。

• 搅拌力度（外部应力）：就像你推车的力气，力气越大，颗粒越想卡死。
• 蚂蚁的活跃度（活性应力）：就像蚂蚁跑得多快。
  • 如果蚂蚁跑得太慢，它们还是会被挤在一起卡死（液体变硬）。
  • 如果蚂蚁跑得刚刚好，它们就能有效地把卡在一起的颗粒“踢开”，破坏那个“脚手架”，让液体变稀。
  • 如果蚂蚁跑得太疯（活跃度太高），它们自己撞来撞去，反而又会让液体变稠（就像一群疯跑的蚂蚁把路堵死了）。

结论：只要控制好“蚂蚁”跑动的强度，就能像调收音机一样，精准地控制这锅液体的粘度。你想让它像水一样流，就调高活跃度；想让它变硬，就调低活跃度。

4. 科学上的“通用公式”

这篇论文最厉害的地方在于，他们发现这种“让蚂蚁乱跑”的方法，和以前科学家发现的“用超声波震动”或“垂直方向摇晃”的方法，本质上是同一个道理。

• 比喻：以前人们发现，如果你一边推车，一边垂直方向猛摇车，车里的东西就不容易卡死。现在发现，让东西自己乱跑，效果是一样的。
• 统一理论：研究人员建立了一个数学公式（标度框架），把“外部推力”和“内部乱跑”统一起来。这意味着，无论你是用超声波、机械震动，还是让颗粒自己动，只要用这个公式算一下，就能预测液体是变硬还是变软。

5. 这对我们有什么意义？

想象一下未来的应用：

• 智能护具：现在的防弹衣或护膝，遇到撞击会变硬保护你。如果加入这种“活性颗粒”，平时穿起来像普通布料一样柔软（因为颗粒在乱跑，不卡死），只有当你受到剧烈冲击且没有额外“活性”干扰时，它才会瞬间变硬。
• 工业管道：在输送高浓度泥浆的管道里，如果管道要堵了，我们可以瞬间激活管道里的颗粒让它们“跑起来”，把堵塞疏通，而不需要停机清理。

总结

这篇论文就像是在说：“别只想着怎么把颗粒粘在一起，不如给它们一点‘自由’，让它们自己动起来，这样就能把原本死板、堵塞的液体，变成听话、流动的液体。”

这是一种通过**“微观的混乱（活性）”来对抗“宏观的秩序（剪切导致的堵塞）”**的巧妙策略。

技术摘要

这是一份关于论文《Tunable shear thickening in active non-Brownian suspensions》（活性非布朗悬浮液中的可调剪切增稠）的详细技术总结。

1. 研究背景与问题 (Problem)

• 剪切增稠现象 (Shear Thickening)： 在高浓度悬浮液（如浆料）中，随着剪切速率或施加应力的增加，相对粘度会显著上升。这种现象通常源于颗粒间摩擦接触的形成，当应力克服颗粒间的排斥力时，颗粒从润滑状态转变为摩擦状态，甚至导致剪切致塞（Shear Jamming）。
• 工业挑战： 剪切增稠会导致管道堵塞、流体流动极度迟缓，给食品加工等工业过程带来困难。
• 现有调控手段的局限： 传统的调控方法通常涉及改变配方（如添加表面活性剂），这无法在流动过程中动态调整。虽然已有研究通过正交叠加剪切（OSP）、振动或声扰动（AP）在流动中引入机械微扰来破坏力链网络，从而降低粘度（即“去增稠”），但利用**活性（Active）**机制（即颗粒自身的自驱动）来动态调控剪切增稠的研究尚不充分。
• 核心问题： 活性颗粒的自推进运动（Self-propulsion）能否作为一种有效的“开关”，在不改变材料成分的情况下，动态地抑制摩擦接触，从而实现对剪切增稠行为的调控？

2. 研究方法 (Methodology)

• 模拟方法： 采用基于粒子的离散元方法（Discrete Element Method, DEM）变体进行数值模拟。
• 系统模型：
  • 颗粒： $N=10^3$ 个非布朗（Non-Brownian）颗粒，具有库仑摩擦系数 $\mu_p=1$。颗粒半径有两种（$a$ 和 $1.4a$），在周期性立方盒子中混合。
  • 流体环境： 悬浮在密度匹配、粘度为 $\eta_f$ 的流体中。
  • 相互作用力：
    • 流体阻力与力矩： 考虑斯托克斯拖曳力和扭矩。
    • 短程流体动力学相互作用： 模拟润滑力（Lubrication forces），防止颗粒直接穿透。
    • 接触力： 当颗粒间距 $\delta_{ij} < 0$ 时，模拟历史依赖的胡克接触力（法向和切向），引入摩擦约束。
    • 排斥力： 颗粒间存在库仑排斥力，防止在低应力下过早接触。
    • 活性力（关键）： 每个颗粒受到随机方向的自推进力 $f_a$。颗粒遵循“运行 - 翻滚”（Run-and-tumble）模型，具有持久时间 $\tau_p$。
• 控制参数：
  • 无量纲施加应力： $\sigma^* = \sigma_{xy} a^2 / f_r$（表征流体应力与排斥力的比值）。
  • 无量纲活性应力： $\sigma^*_a = f_a / (6\pi\eta_f \dot{\gamma} a^2)$（表征活性力与流体剪切力的比值）。
  • 持久时间： 设定 $\tau_p \dot{\gamma} \sim 1$，以模拟与正交剪切频率相当的动态效应。

3. 主要发现与结果 (Key Results)

A. 粘度调控行为

• 剪切增稠的增强与逆转：
  • 在固定活性应力（$\sigma^*_a = 0$）下，系统表现出典型的剪切增稠行为：随着 $\sigma^*$ 增加，粘度平滑上升，随后在临界应力处发生不连续跳跃（DST）。
  • 去增稠（Dethickening）： 当固定高剪切应力（$\sigma^* = 7.29$，此时被动系统处于高粘度摩擦态）并增加活性应力 $\sigma^*_a$ 时，粘度显著下降。
  • 非单调性： 粘度随 $\sigma^*_a$ 的变化呈现非单调性：先急剧下降至最小值，随后由于活性诱导的扩散效应，粘度再次上升。这种“去增稠”效应可以通过调节 $\sigma^*_a$ 进行精确调控。

B. 微观结构机制

• 摩擦接触数量的减少： 模拟显示，活性自推进引入了额外的各向同性动力学，破坏了剪切诱导的力链网络（Force chain networks）。
• 竞争机制： 外部剪切驱动倾向于克服排斥势垒，形成大量摩擦接触；而活性运动倾向于打散这些接触。两者竞争的结果是，在高活性应力下，颗粒间的摩擦接触数（$n_{fc}$）显著减少，系统回到低粘度状态。
• 力链图像： 在无活性或低活性下，颗粒形成贯穿系统的各向异性力链（高粘度态）；在高活性下，力链网络被破坏，接触变得孤立且各向同性（低粘度态）。

C. 普适标度律 (Universal Scaling)

• 标度框架验证： 研究将活性系统的粘度数据映射到最近提出的用于描述机械/声扰动下可调剪切增稠的普适标度框架中。
• 标度变量： 定义了活性系统的标度变量 $\tilde{x}_a = f(\sigma^*) C(\phi) g(\sigma^*_a) / (\phi_0 - \phi)$。
  • $f(\sigma^*)$ 和 $g(\sigma^*_a)$ 分别描述了应力和活性应力的影响。
  • $C(\phi)$ 是各向异性因子。
• 数据坍缩： 将不同体积分数 $\phi$、$\sigma^*$ 和 $\sigma^*_a$ 下的粘度数据坍缩到一条主曲线上。该曲线展示了从润滑主导的粘度发散（指数 $\approx 2$）到摩擦主导的粘度发散（指数 $\approx 1.48$）的交叉行为。这证明了活性调控与正交剪切调控在物理机制上的普适性。

D. 相图演变

• 剪切致塞与 DST 区域的缩小： 在 $\sigma^* - \phi$ 相图中，随着活性应力 $\sigma^*_a$ 的增加：
  • 剪切致塞（Shear Jamming）区域向更高的体积分数移动，意味着在活性存在下，系统更难发生剪切致塞。
  • 不连续剪切增稠（DST）区域的边界向随机密堆（$\phi_0$）移动，DST 区域显著缩小。
  • 尽管活性应力很大，DST 区域并未完全消失（这与某些正交剪切实验不同），因为过高的活性会导致粘度再次上升。

4. 主要贡献 (Key Contributions)

1. 提出活性调控新机制： 首次通过粒子模拟证明，非布朗活性颗粒的自推进运动可以作为一种动态手段，有效抑制摩擦接触，从而在流动过程中“关闭”剪切增稠效应。
2. 揭示微观竞争机制： 阐明了活性各向同性运动与剪切诱导的各向异性力链形成之间的竞争关系，解释了粘度降低的微观物理起源。
3. 建立普适标度理论： 成功将活性悬浮液的流变数据纳入现有的普适标度框架，证明了活性调控与外部机械/声扰动调控在数学描述和物理本质上的统一性。
4. 相图重构： 定量描述了活性应力如何改变剪切增稠和剪切致塞的相边界，为设计具有可控流变特性的智能悬浮液提供了理论依据。

5. 意义与展望 (Significance)

• 工业应用潜力： 该研究为开发新型智能流体材料提供了新思路。通过引入活性成分（如活性胶体或仿生微纳机器人），可以在不改变流体化学成分的前提下，动态调节流体的粘度，解决管道堵塞或加工困难问题。
• 理论统一性： 研究将“活性物质”物理与“颗粒流变学”紧密结合，表明活性可以被视为一种等效的“内应力”或“微扰”，与外部机械扰动具有相同的流变学调控效果。
• 未来方向： 作者指出，除了调节活性力的大小，调节颗粒的持久时间（Persistence time）可能进一步优化相图，甚至完全消除 DST 和剪切致塞区域。未来的工作将探索这些参数的联合调控效应。

总结： 该论文通过严谨的数值模拟，确立了活性自推进作为一种有效的“去增稠”机制，不仅丰富了剪切增稠的理论框架，也为活性物质在复杂流体工程中的应用开辟了新的途径。