湍流怎样从局部扰动，变成能够持续存在的状态？

水一直流，混乱为什么不一定一直在

在一根光滑直管里，让水稳定流动，再在入口附近给一次短促扰动。你可能看到一小段杂乱的流动向下游移动，前后仍是平滑的层流。把这个实验重复很多次，有些局部湍流走得很远，有些突然消失；再提高流速，有的会在下游产生另一团，后来出现一串被平滑区隔开的湍流。1 2

这些局部结构叫 puff，本文称“湍流团”。它们是不断有新水流入、旧水流出的运动结构，并非一包固定的水被搅乱后原样漂走。

这使题目多出一个关键层次：整条管道里始终有湍流，并不要求最初那一团永远不死。 只要它在消失前产生后代，后代继续产生后代，局部事件就可能维持一个长时存在的群体。不过，后代能否留下来还取决于它落在哪里、旁边的流动是否适合它生存。

本文研究光滑直管中这种持续状态怎样出现，再用弯管、加热管和其他剪切流检查它的边界。需要追踪的链条是：平均流如何供能，局部湍流怎样消亡和扩散，以及这些过程怎样组成宏观状态。高雷诺数下的能量级联和Navier–Stokes整体正则性属于其他问题。

三个“开始湍流”，要分别测量

管流常用雷诺数描述流速、管径与黏性的相对关系：

这里的 $U$ 是截面平均流速，$D$ 是管径，$\nu$ 是运动黏度。本文涉及的其他几何可能采用不同速度和长度约定，不能直接把其雷诺数与直管数值大小比较。

“到了某个雷诺数就会湍流”其实可能指三件不同的事。

在通常的直管转捩范围，层流可以对无穷小扰动保持线性稳定，却被合适的有限扰动触发湍流。2011年的实验装置经过控制，在雷诺数约4400时仍能保持层流；研究者再用可控扰动创建湍流团。这说明维持湍流的条件与必然从层流自发触发的条件不能合并。1 2

还有一个看似矛盾、其实测量对象不同的现象：一小块区域从平滑流动变成湍流时，局部波动可以一下变得很强；但整条管道的湍流占比，仍可随着控制参数逐渐从零增加。局部振幅的跳变与空间占比的连续转变能够同时发生。4

一团湍流靠什么活着

先看能量从哪里来。驱动管流的压差不断做功，维持平均速度分布。普通层流里，靠近中心的水流得快，靠近管壁的水流得慢。扰动如果产生横向运动，就能把不同速度的流体搬到不同位置，形成沿流向的快、慢条带；条带和涡旋的相互作用让波动继续从平均剪切中取得能量。黏性则不断耗散这些波动。1

对转捩初期的湍流团，仅说“有剪切就能活”还不够。它的上游界面面对速度剖面较尖的层流，流入界面的动量和能量输运支持局部湍动生产。湍流又反过来搅混流体，使平均速度剖面变平。在较低雷诺数下，单靠这段已被改变的局部流动，未必能维持同样的活动；离开湍流团之后，剖面需要一段距离恢复。1 3 4

可以把过程想成：湍流一面依赖周围流动，一面改变周围流动，改变后的环境又会限制它自身和邻居。 这里的“供能”始终以压差驱动和具体的速度剖面为条件，不是湍流凭空自我制造能量。

这个图像能解释为什么转捩初期会留下有限大小的湍流团：上游有利于活动的维持，下游的恢复过程限制持续扩展。定量上还要比较两个界面的传播速度；平均速度相同时，团可以保持近似大小，速度分离则会改变它的长度，不能单靠一张供能示意图算出团有多长。1 它也给出一个可检查的预期——在一个湍流团下游挤入另一个，下游团接到的流动可能更不利，反而更快消失。3

这不是“扰动越大越湍”的单调故事。相关实验和模型表明，改变平均速度剖面的方式可能同时增强局部波动、削弱后续维持波动的条件。要判断结果，必须跟到能量供应和剖面恢复这两步。1 4

寿命有限，为什么群体仍能持续

2011年的关键进展，是同时比较湍流团消亡与产生后代的时间尺度，而不是只寻找一个单团寿命突然发散的流速。2

在排除初始形成阶段后，所研究管流中单类消亡或分裂事件的等待统计可以近似呈指数形式。这意味着已经活了一段时间，并不会让一个处于相同统计条件的湍流团自动变得更“老”或更“年轻”。用 $\tau_d$ 表示衰亡特征时间、$\tau_s$ 表示分裂特征时间，相应速率为：

提高雷诺数时，测得的平均寿命迅速变长，而分裂等待迅速缩短。2011年研究据两种过程的平衡，把持续湍流的临界点估计在 $\mathrm{Re}\approx2040\pm10$。论文已经明确提醒：空间相互作用会使真正的临界点偏离孤立团两速率相等的位置；两条曲线对雷诺数极其敏感，才使二者在当时实验精度下很接近。2

下面重绘2024年补充材料中孤立团极限的拟合函数。每个“特征时间”都针对相应事件，不是两种事件竞争时“下一件事”发生的共同等待时间。

图上交点与约2040不同，并不构成我们发现了实验矛盾。这里用的是2024年公布的拟合参数；原作者也讨论了不同装置的绝对校准差异，以及参数微小变化如何明显影响交点。更基本的一点是：孤立团速率交点、模型临界点、实验临界点，是三种不同层次的量。 能把模拟的临界位置写出很多小数位，不意味着实验也有同样的绝对精度。3

真正困难的地方是时间。靠近转捩，两类事件的特征时间能达到千万个 $D/U$ 的量级；整个群体达到统计稳定还可能需要更久。有限长管道里，湍流团早已从末端流出，研究者却还没等到足够多次生灭。把“很久没有消失”当作永远持续，正是在这里容易出错。2 3

邻居并不只是帮助传播

2024年的研究把湍流团之间的作用拆成几项，在不同间距下直接测量。3

首先是位置：上游团会改变下游团的传播速度，二者靠近时，下游团有被“推远”的趋势。这是经过流场传递的有效相互作用，不是两团固体碰撞。

其次是寿命：一个湍流团太靠近上游邻居，会更容易消亡。再次是分裂：上游团若在下游方向已经有近邻，新后代的形成会受抑制，平均分裂等待变长。传播并非只要有一次产生后代的事件，就一定增加一个可长期留下来的团。

左图中，把上游间距从141个管径缩短到29个管径，寿命点估计由约 $2.85\times10^6$ 降到 $3.09\times10^3$ 个 $D/U$，相差接近三个数量级。但远距估计的区间很宽，不能把这当作精确的“922倍规律”。右图中，36个管径间距的分裂等待点估计约为134个管径时的3.2倍；这也只是特定雷诺数和间距对的统计比较。3

本库从补充材料提取了48个寿命、44个分裂时间汇总条目，并核对表格页面；图中选取的两组用于展示邻居作用方向。我们没有取得原始逐事件记录，因而这里是公开统计量重绘与计算，不是重新拟合实验寿命分布。

这些作用组合后，群体会出现看似奇怪的组织方式：局部湍流内部依旧混乱，湍流团之间却可能形成较规则的间隔。论文把这种短程排列称为拥堵或类晶体状态。它不意味着水结晶，也不意味着整条管道所有湍流团都锁定不动；临界附近的稀疏状态与更远处的局部有序状态还要分开。3

一个十格小例子：只数出生与死亡为什么不够

暂时把流体细节放下，保留一个可以算到底的示意系统。十个格子围成环，每格只有“活跃”或“空”两种状态。活跃格平均活一个时间单位，并以速率2尝试让另一个格子活跃；目标已经活跃时，这次尝试不增加任何新格子。没有外部重新点火，全部变空以后就无法自行恢复。

我们比较两种传播方式：一种只能随机选左右邻居，另一种可以随机选其余九格中的任何一格。两者每个活跃格的总尝试率相同，差别只有尝试能送到哪里。它们都不是管流的定量模型，格子也没有被指定成多少个管径。

为什么位置重要，不必先等计算结果。若环上只有两个相邻活跃格，向外扩展只有两条活跃—空格边界；把两个活跃格分开，通常会出现四条这样的边界。活跃数量一样，成功传播的机会已经不同。因此，群体数量本身不足以完整预测后续增长。

若忽略这些关联，设活跃比例为 $f$，常见的均匀平均场近似是：

这里 $\lambda$ 是单个活跃格的繁殖尝试率与衰亡率之比。第一项把目标空着的概率直接近似为 $1-f$；取 $\lambda=2$ 时，方程有一个 $f=1/2$ 的正稳态。关键不在求出这个数有多难，而在它把“邻居是谁、是否一起成团”从方程里删掉了。

本库另外枚举十格系统全部1024种构型，求解有限状态主方程，得到下图。初始条件都是十格全活跃；曲线是总体平均，不是挑出的单次轨迹。

只向左右传播时，全灭平均等待约9.54；能向任一其他格传播时，约18.43。到第20个时间单位，仍有活动的概率分别约8.1%和33.5%。同样的单格生灭尝试率，没有给出同样的群体寿命。

这也说明另一个边界：十格、没有外来点火的随机系统，最终都会全灭；有限时间里看起来稳定，不等于严格的无限寿命。统计物理讨论持续活跃相，需要认真处理系统大小、观测时间和初态的极限。这里的有限模型没有推翻这种相变，也没有算出它的临界点，只展示了空间关联和有限系统波动为何不能随手丢弃。

2024年的“定向渗流”结论，到底靠什么成立

定向渗流描述一类活动在空间中传播、沿时间延续，又能局部消亡的非平衡转变。“定向”首先指时间方向：此刻的活动影响未来，未来不会返回去生成过去。它不是说水渗入了某种多孔材料，也不要求所有空间传播只能向一个方向。

在理想化直管描述中，完全层流、没有足以触发湍流的外扰时，是一个吸收态；已有湍流则会波动、衰亡和播种邻近区域。这个对应关系提供了候选理论，但相似的故事还不够，必须检查集体统计是否也一致。

2024年的证据由几个层次接起来：3

1. 实验约束局部规则。 测湍流团间距如何影响速度、寿命与分裂，以及传播位置的随机变化。
2. 模型推进到实验难以达到的尺度。 把湍流团表示成相互作用的点，保留测得的规则；再用大尺度模拟和理论分析检查临界行为。
3. 比较多种集体签名。 不只比较湍流团密度，还比较层流空隙的时间与空间分布、不同条件下的时间演化缩放。
4. 回到实验核可达范围。 通过在出口检测团的排列、在入口重新创建相应排列，近似实现周期边界，观察更长时间的群体行为。

最后一种做法重新创建的是统计意义上的湍流团序列，不是同一小团水绕了一圈。它依赖“重新创建的团在适当发展后可代表原团统计”的条件，也利用了单团近似无记忆的性质。

本文采用的判断是：这套由实验约束的模型，为普通直管转捩的定向渗流解释提供了有分量的支持；最靠近临界点的完整标度，并没有由管道实验直接独立测齐。 原论文明确写道，已收敛的实验仍在拥堵区，未进入所讨论的渐近标度范围；某个更靠近临界的测量持续约一个月，仍无法可靠估计平均密度。

因此，模型与实验互相约束是这项工作的强处，而不是需要藏起来的弱点。它的限制也同样具体：微观规则如何外推、粗粒化是否保留了承重过程、跨装置的雷诺数怎样校准，都会影响结论的适用范围。论文中的多个测量和模拟属于同一研究链，并不自动构成几次独立复现。

改变流动条件后，会走另一条路吗

2026年的研究提供了一个很直接的切口：如果额外作用力改变平均速度剖面，层流区与湍流区之间的供能和共存条件会怎样变化？4

研究者比较弯管实验，以及受加热浮力、指定体积力等影响的管流数值模拟，并考察受磁场作用的导电通道流。在其研究条件下，原先宽广的层流—湍流共存区被显著压窄。作者提出的机制是：额外作用力改变两种状态的速度剖面，降低它们之间的某些剖面差异和空间能量输运，使局部湍流结构更难依赖邻近层流维持；直到更高流速下，湍流才可在更大范围维持。

但不同部分的证据强度不同。弯管实验看到共存区随下游距离越来越窄，尚不能确认管道末端已到统计稳态，也没有在有限管长上证明共存区严格消失。 加热管的周期边界模拟能观察更久，在所测条件下，接近稳态的湍流占比从约0.8以上直接降到零。另一种指定体积力模拟还出现明显的历史依赖：从全湍流开始和从层流—湍流混合态开始，结果可以不同；当层流空隙在亚稳湍流中形成后，它会扩展并消除湍流。4

这些结果支持“空间共存条件能改变转捩方式”，却不能缩成“加热、弯曲或者多加一点力，必然产生同一种跳变”。力的形式、幅度、初始状态、可观察时间以及是否接近另一种线性失稳，都要具体核对。这里展示的是若干受控条件，不是所有受力流动的分类定理。

没有固定大小的湍流团，也不自动意味着不连续

一项值得与2026年论文对读的2024年研究，在平面Couette流数值实验里抑制特定大尺度横向流动，去掉常见的斜湍流带。剩下的活动仍可通过前沿移动扩展，并在内部形成层流空隙。6

作者进一步用单变量随机模型描述这种过程。模型中改变波动强度，便能改变层流空隙的产生，使转变呈连续或不连续；与其流体模拟相配的模型参数支持连续转变。需要保留的限定是：论文没有在那套流体直接模拟中测齐临界标度，普适类判断依赖模型与流体结果的对应。

这项研究与2026年的解释并不完全相同：前者强调前沿与空隙产生的统计，后者强调是什么流场过程控制这些统计。它也提醒我们，去掉某种固定大小的结构，只排除了一种传播组织方式；其他形式的间歇性仍可能维持连续转变。 能量输运、前沿速度和空隙成核，需要放在同一条机制链上检验，不能用其中一个标签代替其他环节。

从管流搬到通道流，连“无记忆”也要重测

2025年的通道流研究测量了完全局域化的湍流条带。条带会有较确定的长度增长，又会发生断裂；它的长度与“有效年龄”关联，消亡和繁殖概率因而不再像所研究管流中的湍流团那样近似恒定。实验与足够大计算区域中的条带模拟，在这一点上相互支持；小区域里无法充分局域化的条带，反倒可能看起来无记忆。5

这是对把一个宏观条带直接当作一个无记忆活动点的限制。它没有单独证明通道流一定不属于定向渗流普适类：微观粗粒化不合适，与更大尺度的最终标度不属于同一类，是两个命题。原作者在讨论中也保留了这一区分。

我们现在能解释到哪一步

对普通直管的转捩初期，较稳的机制链是：外部驱动维持平均流，局部湍流从中取得能量并改变速度剖面；团的消亡、产生后代和邻居相互作用，决定群体能否延续。单团长寿不能代替群体统计，孤立团出生率超过死亡率也不能自动省掉空间关联。

对2024年的管流研究，本文支持其定向渗流解释在明确条件下的价值，同时保留实验可达范围与模型外推之间的接缝。对不同几何和体积力，2024—2026年的新结果表明，改变局部结构、供能关系与空隙产生，就可能改变宏观转捩；目前不足以用一个简单判据包办所有剪切流。

进一步区分这些机制，最有用的观测也很具体：同时追踪界面附近的速度剖面与能量输运、团或条带的大小历史、生灭的条件概率，以及随观察距离和时间变化的湍流占比。若重新选择初态、延长观察、改变计算区域后结论改变，就应先查暂态和几何约束，再讨论普适性。

这篇研究最终解释的是一种持续方式：外部驱动不断供能，再生与传播的过程得以延续，承载过程的局部结构可以不断更替。 这种认识在直管里有明确的物理对象和证据边界；把它用到其他系统之前，仍须重新找到那些真实的供能和传播环节。

本库计算与重绘代码位于 research/calculations/sustaining_pipe_turbulence.py；可下载输入与结果说明、论文时间统计表、有限模型完整时间序列。