混合长度模型

玻尔百科

核心要点

混合长度模型将湍流涡旋比作气体分子，提出流体微团在一个称为“混合长度”的特征距离内输运其动量。
其最大的成功是推导出了对数形式的“壁面律”，它精确地描述了固体表面附近湍流的速度剖面。
该模型的主要局限性在于其局部性，无法解释湍流的输运或历史效应，因此对于非平衡流动是不准确的。
Prandtl 的核心概念在现代计算流体动力学（CFD）中仍然存在，并构成了大涡模拟（LES）中亚格子尺度模型的基础。

引言

湍流是经典物理学中最艰巨的挑战之一，它是一种难以预测的混沌流体运动状态。虽然纳维-斯托克斯方程（Navier-Stokes equations）控制着这种行为，但其在湍流状态下的复杂性使得使用简化模型成为必要。本文深入探讨了其中最具开创性的模型之一：Ludwig Prandtl 的混合长度模型。该模型通过为动量输运提供一个直观的物理类比，解决了我们在量化湍流应力能力方面的关键空白。在接下来的章节中，我们将首先探讨该模型的核心 原理与机制，从其与气体动理论的类比，到其成功推导出壁面律的过程。随后，我们将考察其广泛的 应用与跨学科联系，展示这个百年历史的思想如何在从水力工程到现代气候模拟等领域中仍然是一个至关重要的概念。

原理与机制

理解湍流，就是要去解决经典物理学中最后几个重大的未解难题之一。当流体运动时，其方式有两种根本性的不同。它可以以平滑、优雅、可预测的层次流动——我们称之为层流。或者，它会陷入一种充满涡旋和涡流的混沌、旋转、不可预测的运动状态——湍流。虽然我们有优美的方程——纳维-斯托克斯方程——来描述所有流体运动，但它们在湍流情况下的完整解是如此复杂，以至于我们至今无法掌握。那么，物理学家或工程师该怎么办呢？我们建立模型。而在所有构想出的模型中，最优美、最直观且出人意料地强大的模型之一，便是 Ludwig Prandtl 的混合长度模型。

来自原子世界的类比

要领会 Prandtl 的天才之处，让我们先退一步，思考一个更简单的画面：气体动理论。想象两列长长的火车在平行轨道上以不同速度行驶。它们之间的空气由无数个随机飞驰的微小分子组成。靠近慢车的一个分子，平均而言，拥有慢车的动量。如果它随机跳到与快车一起移动的空气区域，它会把它的“慢”动量带过去。反之，一个来自快车侧的分子可能会跳到慢车侧，携带“快”动量。每一次跳跃都是一次从快流层到慢流层的微小动量传递。这种传递就像一种拖曳力，一种摩擦形式。一个分子在碰撞并分享其动量之前所经过的平均距离被称为平均自由程。这种微观的动量交换是我们所说的粘性的起源。

Prandtl 观察湍流，比如河水中表层流速快于河床附近的情况，并从中看到了一个类比。他看到的不是单个分子，而是大团的、旋转的流体块——涡旋（eddies）。他提出，我们可以将这些涡旋看作是具有内聚力的流体微团，它们被湍流的混沌特性从一个流层踢到另一个流层。

涡旋之舞

这就是混合长度模型的核心所在。Prandtl 做了一个大胆的简化假设：当一个流体微团从高度 $y_1$ 的流层侧向移动到一个新的高度为 $y_2$ 的流层时，它会瞬间保持其原始的平均速度，或者更准确地说，是其在流动方向上的线性动量。这就像慢车上的一个乘客，在跳上快车后的短暂瞬间，相对于轨道，他仍然以慢车的速度移动。

这个位置不当的流体微团，现在会产生一个速度脉动。如果一个来自慢流层（速度为 $\bar{u}(y_1)$ ）的微团向上移动到一个更快的流层（速度为 $\bar{u}(y_2) > \bar{u}(y_1)$ ），它的速度现在就比它的新邻居们小。这会产生一个负的脉动， $u' = \bar{u}(y_1) - \bar{u}(y_2) < 0$ 。如果一个来自快流层的微团向下移动，它会产生一个正的脉动， $u' > 0$ 。

Prandtl 为一个微团在被抹平并与其新环境混合之前所行进的特征距离命名为混合长度， $l_m$ 。这是分子平均自由程在湍流中的类似物。它代表了流动中该位置上起主导作用的、携带动量的涡旋的大小。

从涡旋到应力：一种“湍流粘性”

流层之间流体微团的这种持续交换构成了一种强大的动量输运。在湍流中，这种输运远比层流中引起粘性的分子扩散有效得多。我们给这种湍流动量输运一个特殊的名字：雷诺剪应力，记为 $\tau_t = -\rho \overline{u'v'}$ ，其中 $u'$ 和 $v'$ 分别是流动方向和横向的速度脉动，上划线表示时间平均。

Prandtl 的模型为我们提供了一种计算这种应力的方法，而无需追踪每一个涡旋。速度脉动 $u'$ 与平均速度在混合长度上的变化量成正比，因此 $u' \propto l_m \frac{d\bar{u}}{dy}$ 。同样，假设将微团侧向踢出的横向速度 $v'$ 与 $u'$ 具有相同的数量级也是合理的。结合这些思想，取决于 $u'$ 和 $v'$ 乘积的雷诺应力，必定与该数量的平方成正比。这便引出了著名的混合长度公式：

\tau_t = \rho l_m^2 \left| \frac{d\bar{u}}{dy} \right| \frac{d\bar{u}}{dy}

或者，更常见的是 $\tau_t = \rho l_m^2 \left( \frac{d\bar{u}}{dy} \right)^2$ ，因为在简单的边界层中，速度梯度通常为正。

这个公式通常用一个方便的虚拟概念来表达，称为涡粘性（eddy viscosity）， $\nu_t$ 。我们将其定义为使湍流应力看起来就像层流的粘性应力一样： $-\overline{u'v'} = \nu_t \frac{d\bar{u}}{dy}$ 。将其与 Prandtl 的公式进行比较，可以得出涡粘性为：

\nu_t = l_m^2 \left| \frac{d\bar{u}}{dy} \right|

在这里我们必须非常小心。流体（如水或蜂蜜）的分子粘性是该流体的一种真实属性。你可以在手册中查到它。但是，涡粘性 $\nu_t$ 不是流体的属性；它是流动的属性。它随点而变，取决于当地的速度梯度和涡旋的大小。如一个计算示例所示，大气中10米高度的涡粘性可以比空气的分子粘性大数千倍，这突显了湍流在混合方面的效率有多高。

模型的胜利：壁面律

到目前为止，这是一个优美但未经证实的思想。为了使其具有预测性，我们需要一个混合长度 $l_m$ 的模型。它应该是什么样的呢？Prandtl 推断，在固体壁面附近，涡旋的大小不能超过其到壁面本身的距离。一个巨大的涡旋不可能紧挨着地面存在。最简单的可能假设是，混合长度与到壁面的距离 $y$ 成正比：

l_m = \kappa y

这里， $\kappa$ （kappa）是一个无量纲的比例常数，后来的实验表明其值约为 0.41。它现在被称为冯·卡门常数（von Kármán constant）。

现在是见证奇迹的时刻。让我们考虑壁面附近的区域（比如大气中的地面或管道内部），这里的流动是湍流。在这个“内层”中，一个非常好的近似是总剪应力是恒定的，且等于壁面上的应力 $\tau_w$ 。让我们把我们的简单假设 $l_m = \kappa y$ 代入 Prandtl 的应力公式：

\tau_w \approx \rho (\kappa y)^2 \left( \frac{d\bar{u}}{dy} \right)^2

通过简单地重新排列这个方程，我们得到了对速度梯度的预测：

\frac{d\bar{u}}{dy} = \frac{\sqrt{\tau_w/\rho}}{\kappa y} = \frac{u_*}{\kappa y}

其中 $u_* = \sqrt{\tau_w/\rho}$ 是一个特征速度尺度，称为摩擦速度。这是一个简单的微分方程。当我们对其进行积分以求得速度剖面 $\bar{u}(y)$ 时，我们发现它必须具有以下形式：

\bar{u}(y) = \frac{u_*}{\kappa} \ln(y) + \text{constant}

这就是著名的对数壁面律。仅仅通过两个简单、物理上直观的假设——即携带动量的微团在距离 $l_m$ 内保持其特性，以及这个距离与到壁面的距离成正比——我们就推导出了整个流体力学中最基本、验证最充分的定律之一。这一定律描述了沙漠平原上的风速剖面、河流中的水流以及管道中的油流。这是物理建模力量与美感的一个绝佳例证。

天才的边界：模型失效之处

尽管取得了种种成功，混合长度模型并非一个完整的湍流理论。它的优雅源于其简单性，而这种简单性也正是其最大的弱点。该模型本质上是局部的。它假设空间中某一点的湍流应力仅由同一点的平均速度梯度决定。这隐含地假设了湍流处于局部平衡状态，即由剪切产生的湍流的速率与它耗散成热量的速率立即达到平衡。

这个假设何时会失效？在任何湍流具有“历史”或从一处输运到另一处的流动中，它都会显著失效。一个经典的例子是流经后台阶的流动，就像从卡车后部脱离的流动一样。强烈的湍流在台阶拐角处的高剪切层中产生。然后，这些湍流被扫（平流输运）到下游，进入台阶后的大回流区。在这个区域，平均速度梯度非常小。局部的混合长度模型看到小梯度，会预测湍流应力几乎为零。但实际上，由于所有来自上游的“输入”湍流，应力非常高。该模型没有记忆性；它无法解释湍流能量的输运和历史。

在某些复杂的流动中，当观察到一种称为逆梯度输运的现象时，会出现更深层次的失效。在这些奇特的状况下，动量实际上是从运动较慢的区域输运到运动较快的区域——它“逆流而上”，对抗着平均速度梯度。根据 Prandtl 的模型，应力与梯度的负值成正比，这是不可能的。为了使模型与这样的观察相匹配，混合长度 $l_m$ 必须是一个虚数，这在物理上是荒谬的。

这些局限性并没有减损 Prandtl 洞察力的美妙。它们只是界定了其适用范围。混合长度模型为理解和量化湍流效应提供了卓越的第一步。它教会了我们涡旋作为动量载体的关键概念，以及涡粘性是流动本身属性的思想。对于复杂的非平衡流动，我们现在使用更先进的方法，例如单方程和双方程模型，这些模型显式地求解湍动能及其耗散率等湍流量的输运方程。这些模型赋予了湍流“记忆”，使它们能够解决那些简单而优美的混合长度模型达到其极限的问题。

应用与跨学科联系

既然我们已经掌握了混合长度模型的内部工作原理，让我们退后一步，欣赏一下全局。理解一个原理是一回事；看到它在实际工作中塑造我们周围的世界，则是另一回事。而它所塑造的世界是多么广阔！Ludwig Prandtl 关于“混合长度”的优美而简单的思想，并非尘封在旧教科书中的古老遗物。它是一个活生生的概念，其思想的后继者处于现代工程与科学的核心。从你家管道中流动的水，到塑造我们星球气候的风，混合长度的影子无处不在。

工程学的核心：壁面束缚流

让我们从模型的诞生地开始：附着在固体表面的湍流边界层。想象一下，水在巨大的管道网络中奔流，或者空气流过飞机机翼。在这些“壁面束缚流”中，流体与静止的壁面进行着持续的斗争。混合长度模型告诉我们一些关于这场斗争的深刻道理。通过假设湍流涡旋的大小，即我们的混合长度 $l_m$ ，与离壁距离 $y$ 成正比—— $l_m = \kappa y$ ——该模型预测流速不应线性变化，而应是对数变化！这便产生了著名的“壁面律”，流体力学的一块基石。它解释了为什么湍流管道中的速度剖面比层流的平缓抛物线曲线要平坦和“饱满”得多；核心区域剧烈的混合作用就像一个强大的均衡器，有效地将动量从较快的流层输运到较慢的流层。

但现实世界是复杂的。管道并非完美光滑。那会发生什么呢？模型会失效吗？完全不会！它会适应。我们可以通过告诉模型需要考虑另一个长度尺度来让它了解粗糙度：即粗糙元的高度 $k_s$ 。在非常靠近粗糙壁面的地方，涡旋的大小不能超过它们正在穿越的凸起。模型巧妙地解释了这一点，指出真实的混合长度是离壁距离和粗糙度尺寸之间的一种折衷。这个简单的扩展使该模型成为水力工程师设计运河或机械工程师优化带纹理表面的热交换器的强大工具。

同样，该模型也可以被“教导”关于远离壁面的边界层“边缘”的情况，从而对简单的对数律进行微小但重要的修正。其美妙之处在于模型的灵活性；它提供了一个坚固的框架，可以被充实以捕捉越来越多的物理细节。

进入开放空间：自由剪切流

掌握了“附着”于壁面的流动之后，让我们来释放我们的湍流。考虑一下喷气发动机排出的湍流射流，或者河流中桥墩后方的旋转尾流。这些“自由剪切流”不受固体边界的约束，而是受其自身动量的约束。在这里，混合长度的思想同样提供了直接的洞见。湍流区域的特征尺寸是什么？它就是射流或尾流本身的宽度！通过假设混合长度与流动的局部宽度成正比， $l_m \propto b(x)$ ，该模型成功地预测了这些流动如何扩展以及其中轴线速度如何随距离衰减。它解释了为什么船后的尾流会散开并减弱，这一现象对于从船舶设计到烟囱污染物扩散等所有事情都至关重要。其基本原理是相同的，但背景已经改变，揭示了该模型的统一能力。

超越速度：物质的输运

到目前为止，我们只谈到了动量的混合——这个属性给了我们速度剖面。但湍流是一个不加选择的混合器。它会搅动它能“触及”的一切：热量、化学污染物、海洋中的盐分、空气中的水分。我们的模型能处理这个吗？当然！同样的逻辑也适用。我们可以为温度或化学物质等标量的输运设想一个“标量混合长度” $l_c$ 。然后，该物质的湍流通量可以与动量通量完全相同的方式进行建模。这引出了一个极其优雅的概念：湍流施密特数（Schmidt number，用于质量）或湍流普朗特数（Prandtl number，用于热量），它不过是动量混合长度与标量混合长度之比， $Sc_T = \frac{l_m}{l_c}$ 。这个简单的比率告诉我们，与热量相比，动量被湍流混合的效率是更高还是更低。这个思想是工业过程中传热建模、环境科学中污染物扩散以及河口淡水和咸水混合建模的基础。

能量学的视角：为湍流之火提供燃料

让我们更深入地探讨。从能量学的角度来看，湍流是什么？它是一个以平均流的能量为食的过程。剪切，即相邻流层之间的速度差，是涡旋混沌之舞的最终动力来源。混合长度模型为我们提供了一种直接量化这一过程的方法。能量从平均流转移到湍流脉动的速率——即湍动能的“产生项”（ $P_k$ ）——可以直接用混合长度和流动的平均应变率来表示。从本质上讲，该模型告诉我们，在平均流变形最剧烈、混合涡旋最大的地方，就是湍流之火燃烧得最旺的地方。这使模型从一个纯粹的描述性工具，转变为一个触及流动基本能量学的工具。

现代遗产：从计算尺到超级计算机

你可能会想，“这是一个世纪前的聪明模型。现在我们肯定有更复杂的工具了吧？” 的确如此，但 Prandtl 的思想是如此基础，以至于它们以各种形式，常常是伪装的形式，存在于我们最先进的计算方法中。

当工程师进行计算流体动力学（CFD）模拟时，他们通常从“零方程”或“单方程”模型开始，这些模型都是混合长度假说的直接后代。为了使模型在计算机代码中工作，需要一些巧妙的修正，例如“van Driest 阻尼”函数，它巧妙地让混合长度在非常靠近壁面、粘性效应完全主导时收缩至零。所以，下次当你看到一辆流线型的汽车或飞机设计时，请记住，其空气动力学效率很可能是用带有 Prandtl 原始洞察力 DNA 的工具打磨出来的。

故事并未就此结束。在最前沿的模拟中，即用于天气预报、气候建模和航空航天研究的大涡模拟（Large Eddy Simulations, LES）中，我们并不对所有湍流进行建模。我们解析大的、含能的涡旋，而只对小的、亚网格尺度的涡旋进行建模。我们如何对它们建模呢？通常使用 Joseph Smagorinsky 提出的一个模型，而这个模型——你猜对了——正是一个伪装的混合长度模型！在这里，混合长度与离壁距离无关，而是与计算网格本身的大小有关， $l_{\text{LES}} = C_s \Delta$ 。这是一个深刻的飞跃：混合尺度的物理概念被映射到了模拟的数值参数上。这种现代解释也必须适应现实世界的复杂性，例如，在浮力抑制湍流的稳定分层大气流中，通过减小有效混合长度来实现。先进的“动态”模型甚至允许有效混合长度在空间和时间上变化，自动适应流动的局部物理特性。

最后，我们仰望天空。近地大气层的风廓线由一个复杂的框架——莫宁-奥布霍夫相似性理论（Monin-Obukhov Similarity Theory）——所支配，该理论考虑了地表加热和冷却的关键效应。但是，在一个多云有风、热效应可以忽略不计的日子里会发生什么呢？该理论会急剧简化，出现的是纯粹的对数风廓线。而完美描述这个基本中性极限的理论是什么呢？正是 Prandtl 的混合长度理论，其中 $l_m = \kappa z$ 。从不起眼的管道到广袤的大气层，混合长度模型提供了一条共同的线索，证明了一个简单的物理思想在统一各种复杂现象方面的强大力量。