法夫尔平均

玻尔百科

核心要点

法夫尔平均是一种密度加权技术，它简化了变密度湍流的控制方程。
它巧妙地解决了传统雷诺平均中使问题复杂化的未封闭湍流质量通量项。
简化方程的代价是脉动（或涨落）的定义更为复杂，因为法夫尔脉动的简单时间平均值不为零。
该框架对于模拟燃烧、火箭推进和高超声速空气动力学等复杂现象至关重要。

引言

湍流，即流体混沌且不可预测的运动，是经典物理学中最后几个尚未解决的重大问题之一。为了设计从喷气发动机到天气模型等各种系统，我们无法追踪每一个旋转的涡旋。相反，我们依赖统计平均来理解流动的平均、可预测的行为。标准方法——雷诺平均，对于像水这样的恒密度流体效果极佳。然而，当密度发生波动时——如在火焰、激波或发动机燃烧室中——这种方法会引入一连串新的、复杂的项，从而掩盖了其底层的物理原理。

本文探讨了解决此问题的优雅方案：法夫尔平均。这种强大的、密度加权的视角，为变密度流动的控制方程恢复了简洁性和物理直觉。我们将首先深入探讨原理与机制，比较法夫尔平均和雷诺平均，以揭示视角的简单改变如何驾驭可压缩湍流的复杂性。随后，关于应用与跨学科联系的章节将展示这一理论框架如何成为模拟燃烧、设计高超声速飞行器以及进行先进计算流体力学模拟不可或缺的工具。

原理与机制

要理解流体力学的世界，我们常常必须应对湍流——那种在从翻腾的河流到蜡烛的烟羽等一切事物中都能看到的液体和气体的混沌、旋转之舞。为了理解这种混沌，物理学家和工程师们并不试图追踪每一丝卷流和涡旋，而是采用平均的方法。他们寻求理解流动的平均行为，将其与狂乱的脉动部分分离开来。但这种看似简单的平均行为，背后隐藏着一种微妙而精巧的视角选择，尤其是在流体密度并非恒定时。

两种平均方法的故事：与可变密度的斗争

想象一股稳定的水流过一根管道。水的密度是恒定的。如果我们想描述湍流，最自然的方法就是所谓的雷诺平均。我们在一个固定点长时间测量速度。速度会围绕某个平均值摆动。我们可以将瞬时速度 $u$ 写成其时间平均值 $\overline{u}$ 和一个围绕该平均值波动的脉动值 $u'$ 的和。根据定义，脉动值的平均值为零： $\overline{u'} = 0$ 。这是一种清晰、直观的划分。当我们将这种平均方法应用于流体运动的基本方程（纳维-斯托克斯方程）时，我们发现平均流的方程与原始方程惊人地相似，只是增加了一个称为雷诺应力的新项。该项涉及像 $\overline{u'u'}$ 这样的相关项，代表了由湍流涡旋引起的动量输运。虽然对其建模是湍流理论的核心挑战——即著名的“封闭问题”——但该框架本身是优雅的。

但当流体密度也发生波动时会发生什么呢？想象一下热路面上方闪烁的空气，喷气发动机中燃料和空气的剧烈混合，或者爆炸中热气的迅速膨胀。在这里，密度 $\rho$ 和速度 $\mathbf{u}$ 都在瞬间剧烈波动。让我们看看如果我们固执地坚持简单的雷诺平均会发生什么。

我们从最基本的定律之一：质量守恒定律开始。其守恒形式写作：

\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{u}) = 0

这个方程仅仅说明，一个小体积内密度的变化率等于流过其边界的净质量通量。让我们将我们的平均算子——上划线——应用于这个方程。假设平均与求导可交换（这是一个标准且合理的假设），我们得到：

\frac{\partial \overline{\rho}}{\partial t} + \nabla \cdot (\overline{\rho \mathbf{u}}) = 0

到目前为止，一切顺利。但现在我们必须处理 $\overline{\rho \mathbf{u}}$ 这一项，即平均质量通量。使用雷诺分解， $\rho = \overline{\rho} + \rho'$ 和 $\mathbf{u} = \overline{\mathbf{u}} + \mathbf{u}'$ ，该项变为：

\overline{\rho \mathbf{u}} = \overline{(\overline{\rho} + \rho')(\overline{\mathbf{u}} + \mathbf{u}')} = \overline{\rho}\,\overline{\mathbf{u}} + \overline{\rho'\mathbf{u}'}

将此代回，我们优美简洁的守恒定律长出了一个丑陋的新项：

\frac{\partial \overline{\rho}}{\partial t} + \nabla \cdot (\overline{\rho}\,\overline{\mathbf{u}}) = -\nabla \cdot (\overline{\rho'\mathbf{u}'})

这个新项 $\overline{\rho'\mathbf{u}'}$ 是密度脉动和速度脉动之间的相关项。它是一个湍流质量通量，是平均质量流的一个额外贡献，纯粹由湍流的关联混沌引起。如果我们看动量方程，情况会更糟。对流动量项 $\overline{\rho u_i u_j}$ 会分解成一大堆新的相关项，不仅包括雷诺应力，还包括像 $\overline{\rho'u_i'u_j'}$ 这样的三阶相关项。我们优雅的方程变成了一团糟。在物理学中，这通常表明我们没有从正确的角度看待问题。

视角的转变：按质量平均

在变密度流中，雷诺平均的问题在于它是在空间中的一个固定点上对时间进行平均。但如果我们改变视角呢？与其考虑一个点的平均速度，或许我们应该考虑通过该点的质量的平均速度。这是一个微妙但深刻的转变。我们不再仅仅观察空间中的一个点，而是在追随质量本身。

这引出了密度加权平均的思想，也称为法夫尔平均，以法国物理学家 Alexandre Favre 的名字命名。对于任何量 $\phi$ ，其法夫尔平均（用波浪号表示）定义为：

\tilde{\phi} \equiv \frac{\overline{\rho \phi}}{\overline{\rho}}

你可以将 $\overline{\rho \phi}$ 看作是“ $\rho \phi$ ”属性的总通量，将 $\overline{\rho}$ 看作是质量的总通量。因此， $\tilde{\phi}$ 是每单位质量的 $\phi$ 的平均量。它是流动中一个典型分子所经历的属性的平均值。

现在，让我们重新审视我们的质量守恒方程：

\frac{\partial \overline{\rho}}{\partial t} + \nabla \cdot (\overline{\rho \mathbf{u}}) = 0

根据法夫尔平均速度 $\tilde{\mathbf{u}}$ 的定义，我们有 $\overline{\rho \mathbf{u}} = \overline{\rho} \tilde{\mathbf{u}}$ 。代入后，方程奇迹般地简化为：

\frac{\partial \overline{\rho}}{\partial t} + \nabla \cdot (\overline{\rho} \tilde{\mathbf{u}}) = 0

看！平均量的方程与原始瞬时方程的形式完全相同。那个恼人的湍流质量通量项消失了。当然，它并没有真正消失——它被吸收到平均速度 $\tilde{\mathbf{u}}$ 的定义中，而 $\tilde{\mathbf{u}}$ 现在代表平均质量通量的速度，这在可压缩流中是一个更具物理意义的量。

这种优雅也延伸到了动量方程。当我们对对流项 $\rho u_i u_j$ 进行法夫尔平均时，我们之前看到的混乱的相关项整洁地折叠成两项：一个平均分量 $\overline{\rho}\tilde{u}_i\tilde{u}_j$ 和一个单一的未封闭项 $\overline{\rho u_i''u_j''}$ ，其中 $u_i''$ 是新的“法夫尔脉动”（ $u_i'' = u_i - \tilde{u}_i$ ）。这个未封闭项是法夫尔平均雷诺应力，它代表了由质量加权速度脉动引起的动量输运。封闭问题依然存在，但现在它被清晰地表述在一个单一、明确定义的张量中，类似于不可压缩流中的雷诺应力。我们恢复了守恒定律内在的美和统一性。

天下没有免费的午餐：法夫尔脉动的本质

这种美妙的简化伴随着一个虽小但重要的概念代价。我们必须重新审视我们对“脉动”的概念。对于雷诺平均，脉动 $u'$ 是直观的：其简单平均值 $\overline{u'}$ 为零。那么法夫尔脉动 $u''$ 呢？

让我们对法夫尔脉动进行简单的雷诺平均，即 $\overline{\phi''}$ 。根据定义 $\phi'' = \phi - \tilde{\phi}$ ，我们有 $\overline{\phi''} = \overline{\phi - \tilde{\phi}} = \overline{\phi} - \overline{\tilde{\phi}}$ 。由于 $\tilde{\phi}$ 已经是平均量，其平均值就是它本身，所以 $\overline{\phi''} = \overline{\phi} - \tilde{\phi}$ 。通常，在变密度流中，一个量的雷诺平均和法夫尔平均是不相等的。实际上，通过一点代数运算可以证明，它们的差异与我们之前试图摆脱的湍流质量通量直接相关：

\overline{\phi''} = \overline{\phi} - \tilde{\phi} = - \frac{\overline{\rho'\phi'}}{\overline{\rho}}

所以，法夫尔脉动的简单平均值不为零！这起初可能感觉很奇怪。然而，根据定义，其密度加权平均值为零：

\overline{\rho \phi''} = 0

这就是权衡。我们接受一个更复杂的脉动定义，以换取一套更简洁的平均流方程。这就像一个聪明的会计师重组一本混乱的账簿。总余额保持不变，但条目被以一种使整体结构清晰易懂的方式重新分组。

更大的图景：一个适用于复杂世界的统一框架

法夫尔平均的威力在于它为分析密度变化的湍流提供了一个一致且优雅的框架，这种情况在航空航天工程、燃烧科学和大气模拟中很常见。它清晰地将质量、动量和能量的平均输运与其湍流输运分离开来。此外，在恒密度流的极限下， $\rho'$ 变为零，法夫尔平均变得与雷诺平均相同（ $\tilde{\phi} = \overline{\phi}$ ），整个框架无缝地简化为我们熟悉的不可压缩湍流理论世界。这种一致性是稳健物理理论的标志。

当然，法夫尔平均并不能解决我们所有的问题。它优雅地构建了封闭问题，但并未解决它。并且，可压缩流包含其不可压缩对应物中所没有的物理现象。例如，在湍动能的收支方程中出现了一个新项，称为压力-膨胀相关项， $\overline{p' \nabla \cdot \mathbf{u}'}$ 。该项描述了湍流的动能与流体的内能（热能）之间的可逆交换，就像一个由微观弹簧组成的湍流场被压缩和膨胀。这是一个与粘性耗散完全不同的过程，后者是运动到热量的不可逆、单向转换。

理解和模拟这些可压缩性特有的效应是现代湍流研究的前沿。但得益于法夫尔平均提供的思想清晰性，我们有了一个坚实、统一的基础，可以从此出发进行这些探索，使我们能够解决我们这个时代一些最复杂和最重要的流体力学问题。

应用与跨学科联系

在了解了法夫尔平均的原理之后，我们现在来到了探索中最激动人心的部分：见证这一强大思想的实际应用。正是在这里，在轰鸣的火箭发动机、超声速飞机和火焰复杂舞动的真实世界中，一个物理概念的真正价值才得以显现。法夫尔平均不仅仅是一种数学上的便利；它是一把钥匙，解锁了我们理解、预测和设计一些科学技术中最复杂、最重要现象的能力。它让我们能够在变密度湍流令人眼花缭乱的混沌中，感知到一种隐藏的简洁性。

火箭的轰鸣与火焰的低语

想象一下试图描述蜡烛火焰中的“平均”温度。在任何一个小体积内，你可能会发现冷的、致密的蜡蒸汽和空气，也可能会发现极热的、稀薄的二氧化碳和水蒸气。对该体积进行简单的、非加权的平均会得出一个误导性的数字，一个几乎不存在的“温吞”的温度。然而，真正的物理学关注的是质量和能量的守恒。真正重要的是每单位质量的平均属性。

这正是法夫尔平均所提供的洞见。通过用局部密度 $\rho$ 对我们的平均值进行加权，我们本质上是在追随质量，而不仅仅是体积。这种视角转变的显著结果是，质量、动量和能量的平均方程呈现出一种优美简洁且我们熟悉的形式。平均连续性方程变为 $\frac{\partial \overline{\rho}}{\partial t} + \nabla \cdot (\overline{\rho} \tilde{\mathbf{u}}) = 0$ ，任何输运方程中的对流项都简化为一个平均输运部分和一个单一湍流通量项的和，例如 $\overline{\rho u_i \phi} = \overline{\rho} \tilde{u}_i \tilde{\phi} + \overline{\rho u_i'' \phi''}$ 。那些涉及密度脉动的、困扰传统雷诺平均的相关项的杂乱集合，奇迹般地消失了，被吸收到法夫尔平均量的定义之中。这种数学上的优雅使得模拟湍流燃烧这项艰巨的任务变得易于处理。

这种质量加权观点的威力延伸到了燃烧模拟中使用的复杂统计工具。例如，在“小火焰”模型中，我们常常需要知道在给定的燃料-空气混合比下的平均反应速率或温度。这是一个条件平均。在这里，法夫尔的哲学同样指导着我们。我们不使用基于体积的条件平均，而是用法夫尔条件平均，它问的是：“在这个特定混合比下，一个质量微团的平均属性是什么？”这个以质量为中心的问题是物理上正确的问题，对于准确预测湍流中的化学反应至关重要。

这个框架是如此稳健，以至于它甚至可以扩展到最极端的环境，例如现代液体推进剂火箭发动机中的超临界燃烧。在这些巨大的压力下，液体和气体的区别变得模糊，流体表现为一种性质变化剧烈的、稠密的、非理想的“汤”。虽然焓和温度等变量之间的关系变得远为复杂，但植根于普适的质量守恒定律的法夫尔平均的基本定义保持不变且不可或缺。它提供了一个坚实的基础，我们可以在其上为这些前沿应用构建更复杂的热力学模型。

突破音障

让我们从燃烧的热量转向飞行的速度。当飞机以超声速飞行时，激波会压缩空气，使其密度急剧增加。人们可能期望这种流动中的湍流与低速风洞中熟悉的湍流完全不同。然而，法夫尔平均再次揭示了一种深刻而统一的简洁性，这被概括在一个名为 Morkovin 假说的优美洞见中。

该假说指出，如果流动中的脉动本身不是超声速的（即湍流马赫数 $M_t = \sqrt{2k}/a$ 很小），那么湍流的结构在很大程度上不受平均流高速的影响。即使飞机以 5 马赫的速度划过天空，其表面边界层中的微小旋转涡旋的行为方式也可能与不可压缩流非常接近。

法夫尔平均正是使这一隐藏真理清晰可见的透镜。通过使用质量加权平均，我们有效地“剔除”了与平均可压缩流相关的大尺度密度变化。剩下的部分——对湍流脉动的描述——看起来与不可压缩湍流的方程惊人地相似。这意味着，只要谨慎处理，我们可以将我们经过充分检验的低速湍流模型应用于高速空气动力学领域，这对工程师来说是一个巨大的实践优势。

这一原理不仅对于预测阻力至关重要，对于模拟空气动力加热也同样重要，后者是再入航天器和高超声速飞行器的一个关键问题。热量的传递由能量方程控制，与动量方程一样，应用法夫尔平均简化了平均能量方程。它使我们能够以一致的方式定义和模拟一个单一的湍流热通量项，从而可以预测在极端速度下经历的强烈热负荷。

从理论到虚拟现实

这些强大的理论思想是如何转化为工程师和科学家们日常使用的实用工具的呢？答案在于计算流体力学（CFD）领域，在这里，超级计算机被用来创建“虚拟风洞”来模拟这些复杂的流动。

构建一个实用的 CFD 模型的第一步是解决“封闭问题”。法夫尔平均方程包含新的未知项，如法夫尔平均雷诺应力 $-\overline{\rho u_i'' u_j''}$ ，它代表了由湍流引起的动量输运。一种常见的方法，即 Boussinesq 假设，是将此项建模为与流动的平均应变率成正比，从而引入一个“湍流粘度”或“涡粘度” $\mu_t$ 。这优雅地将未知的湍流应力与已知的平均速度场联系起来。

但是我们如何确定这个涡粘度呢？这就是像著名的 $k$ - $\epsilon$ 模型发挥作用的地方。我们不只是猜测 $\mu_t$ ，而是求解两个额外的输运方程，分别针对湍动能 $k$ 及其耗散率 $\epsilon$ 。为保持一致性，这些方程本身必须使用法夫尔平均量来构建。

然而，大自然有时会带来更多惊喜。虽然法夫尔平均简化了平均流方程，但我们必须持续关注其底层的物理学。在湍流本身高度可压缩（即湍流马赫数 $M_t$ 很高）的流动中，新的物理机制，如“压力-膨胀”，可能变得重要。这代表了湍流脉动与流体内能之间的能量交换。为了捕捉这一点，可压缩性修正被添加到 $k$ - $\epsilon$ 模型中，通常是作为在高 $M_t$ 时激活的新项。

最后，理论在最平凡的层面与实践相遇：设置模拟。为了模拟一股高速气体射流进入一个腔室，工程师必须在边界上指定入口湍流的属性。人们不能简单地凭空设定 $k$ 和 $\epsilon$ 的值。相反，我们使用理论将它们与更具物理直觉的参数联系起来，例如特征湍流长度尺度和湍流马赫数。这使得模拟能够进行一致且具有物理意义的初始化，从而弥合了抽象方程与具体数值结果之间的鸿沟。

最终我们看到，法夫尔平均远不止是一种数学技巧。它是一个指导原则，一种物理洞见，照亮了穿越变密度湍流复杂性的道路。它提供了一种统一的语言来描述像蜡烛火焰、高超声速飞机和火箭发动机这样多种多样的现象，将物理学的基本定律与工程的艺术和科学联系起来。它让我们能够找到潜在的秩序和统一，提醒我们支配宇宙的法则所固有的美。