---
title: "用一条谱线和一段轨迹读活性湍流"
description: "重算 Giomi 的平均场涡量谱，再读 Serra 论文的示踪拉伸图，分清理论曲线、数值场与实验数据。"
canonical_url: "https://sandboxphysics.com/zh/stories/active-matter/reading-active-turbulence/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/active-matter/reading-active-turbulence/"
last_modified: "2026-10-28"
---

[活性物质 · 边学边做实验](https://sandboxphysics.com/zh/stories/active-matter/)

L4 · 第 8 章

# 用一条谱线和一段轨迹读活性湍流

空间谱告诉我们，旋转运动主要分布在哪些尺度；示踪映射告诉我们，一小团流体经过一段时间怎样被拉开。下面固定阅读 Giomi 的 arXiv:1409.1555v2 与 Serra 等人 2023 年论文，并重算前者的一条理论谱线。

Sandbox Physics · 2026-10-28

这一章

[查看八章课程](https://sandboxphysics.com/zh/stories/active-matter/)

开始前

先读取向驱动流动和锚定两章，认识傅里叶波数、速度梯度，以及对空间或多次样本求平均。

## 读曲线前，先找到模型

打开 Giomi 2015 年 8 月 3 日的第二版。方程（1）写出耦合的向列与流体模型，第 I.2 节用弹性应力和活性应力的平衡定义活性长度。平方根关系估计的是一个尺度，不能直接当成生物样品的普适系数。

图 5C–D 画出三种活性下的数值涡量平方谱和能谱，后面的平均场计算用来解释它们的尺度规律。平均场模型把涡看成彼此不相关的圆形对象，涡面积按指数分布。这里重算方程（9）和（14）的归一化理论涡量平方谱，不从图片提取模拟点。

$$
\ell_{\rm a}\sim\sqrt{K/|\alpha|}
$$

## 用两条计算路线得到同一条谱

涡量平方衡量局部旋转的强度，谱则把它分到不同波数上。用特征涡半径缩放波数，再除去整体归一化常数，曲线便成为无量纲形式。它既能用修正贝塞尔函数计算，也能对涡半径做独立数值积分。

下图按这个理论表达式重新计算，独立积分用来检验特殊函数求值。双对数图上，小波数的斜率趋近正一，大波数趋近负二。横轴没有长度单位，纵轴也未按实验标定；读图时应当保留这个无量纲含义。

$$
\begin{aligned}\frac{\Omega(\kappa)}{C}&=\kappa e^{-\kappa^2/2}\!\left[I_0\!\left(\frac{\kappa^2}{2}\right)-I_1\!\left(\frac{\kappa^2}{2}\right)\right]\\&=\frac4\kappa\int_0^\infty t^3e^{-t^2}J_1(\kappa t)^2\,dt.\end{aligned}
$$

![重新计算的无量纲平均场涡量平方谱，小波数斜率趋近正一，大波数斜率趋近负二。](https://sandboxphysics.com/figures/active-matter/giomi-mean-field-spectrum.zh.svg)

按 Giomi 方程（9）和（14）重算的归一化平均场曲线，并以独立积分检验。这是图 5 相关的理论谱线，没有复制或数字化论文模拟点。

完整例题

## 先用一个简单流动检查示踪分析

流映射给每个起点指定一段时间之后的终点。对它求空间导数，就能知道起初很小的一片流体怎样变形；按正确顺序把导数矩阵与转置相乘，得到右 Cauchy–Green 张量。最大特征值的平方根，就是最强的局部拉伸倍数。

取一个平面伸展流：水平方向指数伸长，竖直方向以相同速率缩短，面积不变。算出的最大有限时间 Lyapunov 指数等于伸展速率的绝对值。这个例子没有缺陷，也没有活性应力，专门检验“给定速度之后，示踪分析算得对不对”。

$$
\begin{gathered}\mathbf u=(a x,-a y),\qquad F_T(x_0,y_0)=(e^{aT}x_0,e^{-aT}y_0),\\C=(\nabla F_T)^{\mathsf T}\nabla F_T,\qquad\Lambda_T=\frac{\ln\sqrt{\lambda_{\max}(C)}}{|T|}=|a|.\end{gathered}
$$

## 按物质的来路与去路读 2023 年论文

打开 Serra 等人的论文，DOI 为 10.1038/s41567-023-02062-y，发表于 2023 年 5 月 25 日。图 1 先区分初始位置与终点，再画拉伸和弯折；方程（1）至（3）定义流映射与有限时间拉伸。图 2 把它们用于微管—驱动蛋白实验测得的速度，图示前向和后向映射采用 100 秒时间间隔。

前向图放在初始位置上，追问起初靠近的粒子怎样分离；后向图放在最终位置上，追问这段时间里哪些路径汇聚。先认清标签，再比较脊线。本站浏览器提供合成取向场、速度和示踪轨迹，没有载入这份实验数据，也未计算论文的完整变形分析。

## 下一次比较，先挑能复算的一步

检验理论谱线，就保留涡模型假设，比较两种求值；检验连续场求解器，就用可解的正弦模式并加密网格；研究缺陷控制，就比较有无锚定时配对的轨迹与指定通量。每项检验各自回答一个明确问题。

下次读新论文，先记版本、图号、坐标单位、数据来自实验还是模拟，以及哪些参数已经给出。选一个能独立计算的比较。保存种子、初始场、网格和观察时间，再解释新出现的图案。

带着问题动手

## 先预测，再改变，最后解释

1. 固定弹性，活性大小增强四倍，活性长度估计会变多少？
2. 把重算的理论谱线与两端斜率比较，再在流动 Lab 里同时跟踪缺陷和附近示踪粒子。
3. 分别指出哪些结论来自理论曲线，哪些来自合成轨迹，哪些属于论文的实验观察。

[同时追踪缺陷与流体](https://sandboxphysics.com/zh/labs/active-nematic-flow/) [从取向场找半整数缺陷](https://sandboxphysics.com/zh/labs/nematic-defect-microscope/)

这些 Lab 生成粒子或连续场的合成记录，页面模型说明列出实际动力学和单位。

## 试着自己解释

### 1. 固定弹性，活性大小增至四倍，活性长度估计如何变化？

展开答案与解释

缩短为一半。带摩擦的模型还存在屏蔽长度，不一定能只用一个长度归纳全部行为。

### 2. 归一化理论谱吻合，能说明向列求解器复现了 2015 年整套湍流模拟吗？

展开答案与解释

不能。这一步检验的是特定涡假设下的平均场表达式。动力学复现还要匹配方程、系数、区域、网格、采样方法和论文的数值观测。

### 3. 平面刚体旋转的有限时间拉伸指数是多少？

展开答案与解释

零。流映射的导数是正交旋转矩阵，Cauchy–Green 张量为单位阵，所有拉伸倍数都是一。

## 可以放在旁边读的来源

1. [Giomi, Geometry and topology of turbulence in active nematics, arXiv:1409.1555v2 (2015), Physical Review X 5, 031003](https://arxiv.org/abs/1409.1555v2)

   方程（1）、第 I.2 节和图 2 提供连续场模型与活性长度背景。本站用人为指定的可解场检验流动计算，与论文的湍流模拟分开。

2. [Serra, Lemma, Giomi, Dogic & Mahadevan, Defect-mediated dynamics of coherent structures in active nematics, Nature Physics 19, 1355–1361 (2023)](https://doi.org/10.1038/s41567-023-02062-y)

   图 1 和方程（1）至（3）介绍流映射与有限时间拉伸；图 2 使用微管与驱动蛋白实验的速度场。这里作阅读引导，不复制图片，也未导入实验数据。

[上一章 让细棒偏爱一个方向，能把缺陷引开吗](https://sandboxphysics.com/zh/stories/active-matter/anchoring-and-confinement/)