---
title: "弹性失稳：四辊磨中的对称破缺互动实验｜Sandbox Physics"
description: "一个完全没有惯性的缓慢高分子流动，怎么会失去对称、开始振荡？靠的只是弹性应力。"
canonical_url: "https://sandboxphysics.com/zh/labs/four-roll-elastic-instability/"
language: "zh-CN"
page_type: "lab"
alternate_languages:
  en: "https://sandboxphysics.com/labs/four-roll-elastic-instability/"
last_modified: "2026-10-29"
---

FR26 · 弹性失稳

# 弹性失稳：四辊磨中的对称破缺

**一个完全没有惯性的缓慢高分子流动，怎么会失去对称、开始振荡？** 靠的只是弹性应力。中心驻点把高分子拉成一条又细又强的应力带，这股应力反过来推动流动。Weissenberg 数超过约 5 后，应力带的微小偏斜不再消失而是增长，四个涡失去对称；到 10 左右，流动持续振荡，把各涡之间的流体混在一起，全程没有惯性。

四个辊子在没有惯性的条件下搅动高分子溶液。Weissenberg 数越过临界值后，中心被拉伸的高分子打破对称，流动开始振荡。

> 交互仿真：在浏览器中打开本页，即可调整参数并运行模型。

模型假设

模型是 Thomases 与 Shelley 的二维周期 Stokes–Oldroyd-B 四辊磨：固定体力驱动四个反向旋转的涡胞，中心是一个伸长驻点；没有惯性，高分子与溶剂的黏度比为二分之一。高分子应力带有很小的人工扩散，与已发表的计算相同；它会抹平细薄的应力片，并改变长时间动力学。Lab 在 64 乘 64 网格上运行，比已发表的 256 乘 256 计算粗。它的对称破缺阈值和 Weissenberg 数为 10 时的演化与我们的 128 乘 128 计算一致，但应力带形成时拉伸的最小特征值会短暂低于零，说明这一阶段分辨率不足。弹性湍流的弯曲流道与旋转圆盘实验没有建模。所有场都是模拟，未展示实测数据。

探针处的流速

离对称有多远

第 I 象限的染料从哪来

更多数值与计算方法

- **模型时间**: $0$
- **高分子耦合系数**: $0$
- **人工应力扩散**: $0$
- **网格与时间步长**: $0$
- **中心驻点的位置**: $0$
- **最大拉伸（构象张量的迹）**: $0$
- **本次运行以来拉伸的最小特征值**: $0$
- **Stokes 能量恒等式的相对残差**: $0$
- **应力谱在最高三分之一波数中的份额**: $0$
- **对流 Courant 数**: $0$
- **染料释放时刻**: $0$
- **每秒推进的模型时间**: $0$

不对称量是应力与它在辊子驱动力各对称操作下的平均之差的均方根；对称态为零。混合指数在染料均匀分布于四个象限时为一，仍留在起始象限时为零。最小特征值为负说明那里的拉伸场不符合物理，即网格对应力带太粗。所有曲线都来自这次模拟，不是实测数据，也不是文献曲线。

物理教程

## 没有惯性也会失去的对称

**背景** 四个反向旋转的辊子以极慢的速度搅动液体，惯性完全不起作用。普通液体此时严格跟随驱动力：流动与搅动有同样的对称，并且一直不变。加入长链高分子后，同样缓慢的搅动却能让流动偏向一侧，并持续振荡。

**为什么值得理解** 四个辊子中间有一个驻点，流动在那里无限制地拉伸高分子。被拉开的高分子形成一条细薄的应力带，反过来推动流动。这是只靠弹性应力打破对称的最简单场景，也是通往弹性湍流的第一步。

先抓重点

- **核心问题**: 一个完全没有惯性的缓慢高分子流动，怎么会失去对称、开始振荡？
- **一句话直觉**: Weissenberg 数 $\mathrm{Wi}$ 比较的是高分子记住一次拉伸的时间与流动时间；流动时间指纯溶剂流在中心把流体拉长到 e 倍所需的时间。$\mathrm{Wi}$ 越过约 4.9 以后，对称的应力带不再稳定：应力带稍一倾斜，流动就跟着弯曲，并把它朝同一方向拉得更歪。

### 核心数学模型

#### 由辊子和高分子共同驱动的慢流

$$
-\Delta\mathbf u+\nabla p=\beta\,\nabla\cdot\mathbf S+\mathbf f,\qquad \nabla\cdot\mathbf u=0
$$

Stokes 流：黏性与固定的辊子驱动力 $\mathbf f$ 以及高分子应力平衡。当 $\mathbf S=\mathbf I$ 时，流动是围绕中心驻点的四个对称涡胞。

#### 高分子怎样被拉伸、怎样松弛

$$
\partial_t\mathbf S+\mathbf u\cdot\nabla\mathbf S-\nabla\mathbf u\,\mathbf S-\mathbf S\,\nabla\mathbf u^{\mathsf T}=-\frac{\mathbf S-\mathbf I}{\mathrm{Wi}}+\nu_p\Delta\mathbf S
$$

Oldroyd-B 构象张量 $\mathbf S$ 被流动携带、拉伸，又向 $\mathbf I$ 松弛。很小的扩散 $\nu_p$ 是数值手段，已发表的计算也这样做，它不是真实高分子的性质。

#### 驻点处的拉伸

$$
\frac{d S_{xx}}{dt}=2\dot\varepsilon\, S_{xx}-\frac{S_{xx}-1}{\mathrm{Wi}},\qquad \dot\varepsilon\,\mathrm{Wi}>\tfrac12\ \Rightarrow\ S_{xx}\ \text{grows exponentially}
$$

$\dot\varepsilon$ 是以流动时间为单位的中心应变率：只有溶剂时为 1，高分子开始抵抗后变小，在 $\mathrm{Wi}$ 很小的极限下为 $2/3$。越过 $\dot\varepsilon\,\mathrm{Wi}=1/2$ 后拉伸快过松弛；没有扩散时，中心应力会无界增长，Thomases 与 Shelley（2007）发现这从 $\mathrm{Wi}\approx0.9$ 开始。这里让应力带保持有限的是扩散和它对流动的反馈。

#### 度量失去的对称

$$
A=\left\langle\left\|\mathbf Q-\mathcal P\mathbf Q\right\|_F^2\right\rangle^{1/2},\qquad \mathbf Q=\mathbf S-\mathbf I
$$

$\mathcal P$ 把应力在辊子驱动力的各个对称操作下求平均。对称态的 $A$ 为零；对称应力带一旦失稳，$A$ 就指数增长。

### 常见难点

#### 失稳一定要靠惯性？

**容易误解** 这么慢的流动不会失稳；失稳和湍流需要很大的 Reynolds 数。

**正确理解** 这里的 Reynolds 数按构造就是零。驱动失稳的非线性在高分子方程里：流动拉伸高分子，被拉伸的高分子又让流动弯曲。

#### 这些数值已经收敛

**容易误解** 这里显示的阈值和周期就是模型的精确答案。

**正确理解** 阈值 4.86 在我们的 64 乘 64 和 128 乘 128 网格之间一致，也与已发表的 4.8 到 4.9 一致。在默认网格上，应力带形成时 $\mathbf S$ 的最小特征值会短暂变负，说明这一阶段分辨率不足。改变 $\nu_p$ 会改变长时间的行为。

### 动手实验

1. 01

   #### 看振荡

   从“会搅混的振荡”开始，看探针曲线和彩色染料。

   **你会看到：** 明亮的应力带偏离中心来回摆动，探针速度以一长一短两个周期起伏；第 II、IV 象限的染料不断进入第 I 象限，对角的第 III 象限却很少进来。

2. 02

   #### 待在阈值以下

   选择“对称的应力细带”，再把 Weissenberg 滑块从 4 拖过 5。

   **你会看到：** 在 4 时应力带保持笔直，不对称量衰减。越过约 4.9 后它转而上升；到 6 时它在对数坐标上沿直线爬升，几百个流动时间后应力带倾倒。

3. 03

   #### 把应力抹平

   选择“应力扩散过强”，与参照的振荡对比。

   **你会看到：** 扩散加到已发表值的十倍后，对称破缺推迟到约 590 个流动时间，之后只剩微弱的缓慢摆动。

物理结论

只靠弹性应力就能打破慢流的对称：当高分子记住一次拉伸的时间超过约五个流动时间，中心的应力带就会倾倒，流动开始振荡并搅混。

## 相关实验与文章

复杂流体 · 阅读路径

### 松弛、拉伸、伸长与弹性失稳

这条阅读路径里用到本实验的文章，后面跟着它们推荐的其他实验。

1. [文章 **弯曲流线上的法向应力怎样引发弹性失稳**](https://sandboxphysics.com/zh/stories/complex-fluids/normal-stresses-and-curved-streamlines/)
2. [实验 **线团—伸展转变：剪切与伸长流中的聚合物** 读出法向应力与剪切应力之比](https://sandboxphysics.com/zh/labs/polymer-stretch-startup/)
3. [文章 **四辊流的驻点怎样失去对称，数值方法又改变了多少**](https://sandboxphysics.com/zh/stories/complex-fluids/elastic-instability-at-a-stagnation-point/)
4. [文章 **读懂弹性湍流：没有惯性的混乱流动有哪些证据，还有什么没解决**](https://sandboxphysics.com/zh/stories/complex-fluids/reading-elastic-turbulence/)

运动与动力系统 · 入门路径

## 从一次下落到一条轨道

先测重力，再分解运动，加入接触力与碰撞，最后把同样的定律带到轨道上。

[从这里开始 → **自由落体与重力测量** 采集几次下落，用相机记录的位置拟合重力加速度；再加入空气阻尼，看拟合结果怎样偏离真值。 观察 · 14 分钟](https://sandboxphysics.com/zh/labs/free-fall/)

1. [1 自由落体与重力测量](https://sandboxphysics.com/zh/labs/free-fall/)
2. [2 抛体飞行设计器](https://sandboxphysics.com/zh/labs/projectile-flight/)
3. [3 斜面与摩擦](https://sandboxphysics.com/zh/labs/inclined-plane/)
4. [4 一维碰撞台](https://sandboxphysics.com/zh/labs/collision-bench/)
5. [5 弹簧振子与周期标定](https://sandboxphysics.com/zh/labs/mass-spring/)
6. [6 轨道是怎样形成的](https://sandboxphysics.com/zh/labs/orbit/)

[查看完整路径 →](https://sandboxphysics.com/zh/topics/mechanics/#starter-path)