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

读曲线前，先找到模型

打开 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}

重新计算的无量纲平均场涡量平方谱，小波数斜率趋近正一，大波数斜率趋近负二。 — 按 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 秒时间间隔。

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

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

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

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

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

读曲线前，先找到模型

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

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

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

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

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

试着自己解释

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

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

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

可以放在旁边读的来源