玻尔百科从管道中油气的翻腾，到碳酸饮料的气泡升腾，涉及多相——气、液、固——的流动无处不在。当这些流动变为湍流时，它们便转化为流体动力学中最复杂、最具挑战性的现象之一。单相流动的简单、可预测的规律，让位于界面、涡旋和能量交换的混沌之舞。这种复杂性带来了一个巨大的知识鸿沟：我们如何能准确预测和控制一个看似难以简单描述，却对无数工业过程和自然现象至关重要的系统？本文旨在应对这一挑战，带领读者开启一场进入湍流多相流世界的概念之旅。
在第一部分“原理与机制”中，我们将剖析其基础物理学，探索不同流型如何出现，各相如何与湍流相互作用，以及我们用来模拟这种复杂性的策略。随后，在“应用与跨学科联系”中，我们将见证这些原理的实际应用，揭示它们在生物过程、材料科学乃至天体物理学等不同领域中的关键作用。我们的探索将从支配这种复杂行为的基础元素开始。
想象一下，你正在观察一个大型水平管道的横截面，管道中输送着来自海上油井的原始混合物。流速缓慢，近乎慵懒。你会看到什么？你会发现大自然简单而优雅的组织法则在起作用：重力。这种由天然气、原油和盐水组成的混合物，会沉降成整齐的水平层。在最底部，你会发现密度最大的流体——盐水。其上是一层较轻的原油，而天然气则漂浮在最顶层。这种稳定的、按重力分层的排列，就是物理学家所说的分层流。它是所有流型中最简单的一种，是一种基准状态，其中一切都处于其能量最有利的位置。
但是现在，让我们加大马力，提高流速。这个平静分层的世界被打破了。比液体移动快得多的气体，开始将油的表面搅起波浪。这些波浪不断增长，很快，其中一个变得足够大，触及管道顶部，完全桥接了整个横截面。这个巨大的波浪，一个由油和气泡组成的翻滚、泡沫状的混合物，像攻城槌一样被推向管道下游。这就是段塞流，它与平静的分层流判若云泥。观察管道上的一个固定点，你会看到液位随着这些液体段塞与大气穴的交替出现而急剧升降。流动不再是时间上恒定的，所以我们称之为非定常流。在任何一个瞬间，沿管道长度方向的景象都各不相同，所以我们称之为非均匀流。某些部分，如段塞突兀的前缘，在很短的距离内发生变化（急变流），而其他部分，如段塞之间液膜的平缓斜坡，则变化缓慢（缓变流）。
这种从宁静分层到剧烈脉动野兽的戏剧性转变，揭示了我们主题的核心。静水力学的简单规则让位于复杂的力量之舞，而所有的动作、所有的戏剧性，都源于界面——那个气液相遇、不断变化的边界。
界面不仅仅是一个被动的边界；它是一个动量交换的活跃战场。快速移动的气体拖拽着液体，而缓慢移动的液体则阻碍着气体。正是这场拉锯战创造了我们观察到的波浪、段塞和所有其他复杂的模式。理解这种相互作用就是理解一切。
但是，我们该如何开始描述这样一个混乱、波动的系统的几何形状呢？想象一下，一种简单的单相流体在非圆形管道中流动。为了使用我们熟悉的圆形管道摩擦力公式，我们发明了一个名为水力直径（）的概念。这是一种为任何形状定义一个“有效”直径的巧妙方法。它直接源于力的基本平衡：推动流体前进的压力与壁面拉回它的剪切应力相平衡。对于任何管道，这种平衡都给出了一个非常普适的定义：，其中 是流动的横截面积， 是湿周长——流体接触的管壁长度。
现在，让我们回到我们的两相流。这里的“周长”是什么？绝妙的是，答案是：这取决于你在问什么！
假设我们对液相所受的曳力感兴趣。液体被管壁的摩擦力和流过其表面的气体剪切力所阻碍。对于液体的动量而言，界面就像另一面“墙”。因此，要为液体定义一个动量等效直径，我们必须在计算中同时包含湿壁周长（）和界面周长（）：。
但如果我们感兴趣的是液体从热管壁吸收了多少热量呢？气液界面并不从管壁传递热量。只有与液体直接接触的那部分管道才起作用。所以，对于传热而言，界面是无关紧要的，正确的等效直径是 。这是一个优美而深刻的观点：问题的几何形状本身会根据我们所考虑的物理过程而改变。不存在一个“真正”的等效直径；只存在适用于正确问题的正确直径。界面对于某些现象是边界，而对其他现象则是不可见的。
这些流动中的运动不仅是非定常的；它还是湍动的——一个由各种尺寸的旋转涡旋组成的混沌漩涡，层层嵌套。试图追踪湍流中每一个流体质点的运动是一项徒劳的任务，就像试图绘制飓风中每一滴雨的路径一样。因此，我们采用一种称为雷诺平均的统计方法。我们对流动特性进行时间平均，以获得平均流动的图像。
但是，当我们对运动方程进行平均时，会发生一件有趣的事情。描述流体如何承载自身动量的动量输运非线性项 ，在平均后会留下一个残余项。这个残余项，一个被称为雷诺应力张量的新项，代表了由湍流脉动本身引起的动量净输运。它是我们平均掉的涡旋的幽灵。这个项是未知的；我们已经平均掉了计算它所需的信息。这就是著名的湍流封闭问题。整个湍流建模领域，在本质上，就是一场寻找智能“封闭关系”的探索，用已知的平均量来表达这些未知的湍流应力。
在湍流多相流中，这个故事又增添了引人入胜的新篇章。分散相——气泡、液滴或颗粒——并不仅仅是随波逐流。它主动参与湍流，与连续相流体进行动态的相互作用。
想象一下一群气泡在液体中上升。每个气泡在穿过流体时，都会在身后留下一串扰动、旋转的水流。这些尾迹是新湍流的来源。气泡不断地搅动液体，向湍流级联注入能量。因此，气泡流本身就比单独的液体更具湍动性。为了在我们的模型中捕捉这一点，我们必须在方程中加入一个特殊的气泡诱导湍流生成项。该项告诉我们气泡正在产生多少湍动能 。这些能量从何而来？它来自于气泡滑过液体时曳力所做的功。这种能量注入的速率与滑移速度的立方（）成正比，并且对于固定的气体体积，与气泡直径（）成反比。这意味着一群微小气泡比几个大气泡是更有效的湍流生成器！
现在，考虑相反的情景：一种载有重质固体颗粒的流体，比如水中的沙子。带有湍流涡旋的流体试图搅动这些重颗粒。但这些颗粒具有惯性“惰性”——它们抗拒被加速。流体必须对颗粒做功才能拖动它们，而这个功直接从湍流涡旋中消耗能量。颗粒充当了一个能量汇，抑制了湍流，使流动变得不那么混乱。因此我们得到了一个美丽的对称性：轻的、快速移动的气泡产生湍流，而重的、响应慢的颗粒破坏湍流。第二相从来不只是一个乘客；它是流动基本特性的主动调节者。
因此，我们面对的是一个由界面处的复杂物理驱动、并由自身被各相修饰的湍流所搅动的狂野流型景观。我们如何才能期望预测这样一个系统的行为呢？这就是建模的艺术与科学发挥作用的地方，它为我们在计算机内部构建一条“虚拟河流”提供了不同的蓝图。不同的方法反映了物理保真度与实践简单性之间的基本权衡。
一种方法是务实的工程捷径，典型的例子是经典的 Lockhart-Martinelli 模型。该方法巧妙地回避了繁琐的细节。它不试图直接计算界面曳力和湍流相互作用，而是提出了一个更简单的问题：这种两相流的压降与气体或液体单独流动时的压降有何关系？然后，它使用一个巧妙的经验关联式，一个“摩擦乘子”，将所有复杂的两相物理打包成一个单一的修正因子。这是一个非常有用的工具，但它并没有真正“解释”物理；它只是将其关联起来。
在另一个极端是雄心勃勃的双流体模型，通常称为欧拉-欧拉模型。在这里，我们卷起袖子，为每个相分别写下质量、动量和能量守恒的基本定律。我们将气体和液体视为两种相互渗透的流体，每种流体都有自己的速度和温度。但现在这两组方程是相互脱离的。为了让它们相互“对话”，我们必须提供描述所有界面交换的封闭定律。我们需要一个曳力模型，一个相间传热模型，一个蒸发或冷凝速率模型，当然，还需要各相如何调节湍流的模型。我们刚才讨论的所有物理——气泡诱导的生成，颗粒诱导的衰减——都成为这些方程中的显式项。这种方法功能强大且物理意义深刻，但其预测的准确性完全取决于我们输入的封闭定律的好坏。
最后，还有第三种哲学，即流体体积（VOF）模型。这就像一个高分辨率数码相机。它不进行任何平均。它通过标记哪些计算单元格充满了液体、哪些充满了气体、哪些位于边界上，来精确追踪界面的位置。然后，它为组合的“混合流体”求解一组单一的动量方程。VOF 可以生成流动极其精细的图像，捕捉波浪和液滴的精确形状。但这种保真度是有代价的：它需要非常精细的计算网格和巨大的计算能力，这使得它对于非常大规模的系统来说不切实际。
归根结底，没有一个单一的“最佳”模型。蓝图的选择取决于我们希望回答的问题。我们是在为一条长管道寻找压降的快速粗略估算吗？像 Lockhart-Martinelli 这样的经验模型可能就非常完美。我们需要理解滑移速度如何影响湍流和相分布的详细机制吗？双流体模型就是我们的工具。我们想看到单个液滴撞击表面时那美丽而复杂的飞溅吗？那么我们转向 VOF。从观察简单的分层流到构建这些复杂的计算工具的历程，证明了我们不懈地努力去理解和预测这个复杂、湍动而又美丽的多相流世界。
现在我们已经探讨了不同物质形态间这场美丽而混沌之舞的基本原理，你可能会好奇：我们实际上在哪里能看到这些现象？如果这一切都只是学术操练，那么理解气泡、液滴和颗粒之间复杂的华尔兹又有什么用呢？答案是，我们几乎无处不见——这证明了大自然深刻的统一性。从一杯苏打水中的气泡到恒星的形成，湍流多相流的相同基本原理都在发挥作用。
在本章中，我们将踏上一段旅程，离开纯粹的理想化方程世界，去探索那个混乱、复杂而迷人的真实世界。我们将看到，掌握湍流多相流不仅是工程师和科学家的工具，更是一个变革性的视角，通过它我们可以理解和操纵微观和宇宙尺度的世界。我们的旅程将从工业机械的心脏地带延伸到天体物理学的前沿，在每一站都揭示我们刚刚学到的概念的力量与优雅。
现代工程的大部分故事都是关于控制的故事——驯服自然的狂野力量以满足人类的需求。在湍流多相流领域，这一点尤为真实，因为成功往往取决于我们精确塑造混乱的能力。
想象你是一个微生物，比如一种细菌或酵母细胞，在一个巨大的工厂里被用来生产救命的抗生素或生物燃料。像任何生物一样，你需要呼吸。对于许多这样的微生物“主力军”来说，这意味着要获得稳定的氧气供应。这里存在一个巨大的挑战：氧气出了名地不爱溶于水。我们如何确保装在巨大钢罐里的数万亿个细胞不会窒息？
答案在于生物反应器的设计，这是一个多相工程的奇迹。其策略在概念上很简单：向罐中鼓入空气。但要使其奏效则是一门艺术。氧气从气泡进入液体的速率取决于两个关键因素：所有气泡的总表面积，以及跨越每个气泡表面薄液膜的输运效率。工程师们将这些效应汇总成一个至关重要的参数，即体积传质系数，或 。为了让我们的微生物大军保持愉快和高产，我们需要使这个 尽可能大。
这是如何做到的呢？通过明智地应用湍流。通过用叶轮剧烈搅拌培养液，我们实现了两个目标。首先，剧烈的搅动运动将大的、懒散的气泡粉碎成一团细小、闪亮的微泡云，从而极大地增加了总界面面积 。其次，湍流不断地冲刷每个气泡的表面，使停滞的液膜变薄，从而提高了传质系数 。
但并非所有湍流都是平等的。如果你只是在一个高大的圆柱形罐子里放一个螺旋桨而没有任何其他装置，你会得到一个糟糕的设计。液体只会像水流入下水道一样打旋。气泡被吸入涡流中，然后迅速冲出罐子，几乎没有与液体接触。顶部和底部的细胞会因停滞而“挨饿”。这正是真正的工程巧思发挥作用的地方。一个设计精良的生物反应器使用挡板——沿着罐壁的垂直板——来打破这种懒散的涡旋，并将旋转运动转化为复杂、混沌且高效的混合模式。此外，它不使用一个大叶轮，而是在轴上分级安装多个较小的叶轮。这使得能量输入更均匀地分布，确保整个体积都处于一个粉碎气泡、翻搅液体的湍流漩涡中。正是通过以这种方式塑造流动，从低效的涡旋转变为均匀、各向同性的湍流状态，工程师才能为大规模生命创造一个理想的环境。
让我们从孕育生命转向一个看似相反的过程：液体的剧烈分解。无论是为了冷却滚烫的计算机芯片，将燃料喷入发动机气缸以实现高效燃烧，还是将牛奶喷雾干燥成粉末，我们常常需要将大块液体变成细小的液滴薄雾。这个过程称为雾化，是多相流动力学的纯粹应用。
实现这一魔术的设备是喷嘴，它们的设计千奇百怪，每一种都利用不同的物理机制来撕裂液体。考虑三种典型设计：
每种设计都是一个量身定制的解决方案，通过操纵惯性力、粘性力和表面张力的相互作用来实现所需的喷雾特性，将界面不稳定性的原理转化为强大的技术工具。
到目前为止，我们讨论的都是希望增强相间相互作用的情况。但当这种相互作用导致不良后果时会发生什么呢？一个典型的例子是污垢，即不希望的物质——如矿物水垢、铁锈或生物黏膜——在表面逐渐积聚。在换热器中，这种积聚物就像一层绝缘毯，严重削弱其性能。在管道中，它会收缩流动通道，需要越来越大的泵送功率。这几乎是每个行业都面临的一个价值数十亿美元的头痛问题。
污垢是一个典型的湍流多相问题，一场在流固界面上展开的战斗。在这里，湍流扮演着一个有趣而矛盾的双重角色。一方面，湍流涡旋是一种高效的“递送服务”。它们将溶解的致垢物或悬浮颗粒从主流区输送到近壁区，促进沉积。这种递送服务的强度由传质系数决定，正如我们在生物反应器中看到的，它随湍流的增强而增加。另一方面，正是相同的湍流在壁面产生剪切应力，这是一种“清洁服务”，试图冲刷掉任何新沉积的物质。净污垢速率是这两种相反效应竞争的结果。
几十年来，对抗污垢的主要策略就是简单地增加流速，希望增强“清洁服务”的效果超过“递送服务”。但现代科学提供了一个更优雅的解决方案：如果我们能重新设计表面本身呢？这就是超疏水表面的前景。这些是微结构化的材料，当浸入水中时可以捕获一层空气（即“Cassie-Baxter 状态”），从而创造一个几乎无摩擦的光滑界面。
这有什么帮助呢？通过在壁面产生滑移，该表面从根本上改变了近壁湍流的结构。对于给定的流率，光滑表面需要较小的驱动压力，这意味着壁面剪切应力——即清洁服务的力量——减小了。但更重要的是，整个近壁区的湍流引擎被抑制了。湍动能的产生减弱了。这对递送服务产生了深远的影响。不仅对流输运减少了，而且一种称为“湍流泳”的微妙机制（它将颗粒推向湍流梯度下降的方向）也减弱了。最终结果是，壁面实际上被“饿死”了致垢物。懒散的递送服务比减弱的清洁服务影响更大，从而导致污垢的大幅减少。
当然，这也是有代价的。如果压力过高或表面受损，脆弱的空气层可能会塌陷，水将填充微观缝隙（即“Wenzel 状态”）。表面会从超滑转变为超糙，这会极大地增加阻力，并为污垢附着提供了完美的避风港。这种表面科学、材料工程和流体力学之间的联系是一个充满活力的研究领域，预示着未来将出现自清洁管道和永久高效的换热器。
湍流多相流一些最引人注目的应用发生在那些运行在刀刃上的系统中，在这些系统中，一个小小的变化就可能引发突然的、往往是灾难性的转变。理解和预测这些不稳定性对于安全性和可靠性至关重要。
沸腾是传递热量最有效的方式之一，这就是为什么它在从蒸汽发电厂到高性能电子设备冷却系统的所有事物中都处于核心地位。当你提高浸入液体中的受热表面的温度时，气泡会越来越猛烈地形成，带走大量的能量。但是有一个极限。在某个点，即所谓的临界热通量（CHF），蒸汽的产生变得如此之多，以至于气泡合并成一层连续的蒸汽膜，覆盖整个表面。这层蒸汽膜是热的不良导体，起着绝缘层的作用。热量无法再散发出去，表面温度可能在几秒钟内急剧飙升，导致物理烧毁。这就是“沸腾危机”。
是什么导致了这种灾难性的转变？早期的流体动力学模型，由 Zuber 开创，将其设想为一种“蒸汽交通堵塞”。离开表面的蒸汽射流受到界面不稳定性（特别是 Kelvin-Helmholtz 波）的阻碍，这类似于导致液体射流破裂的波。当蒸汽速度过高时，这些波会增长，有效地堵塞了逃逸路径，导致蒸汽积聚成膜。但这幅图景完整吗？如果蒸汽射流本身是湍流的呢？这就引入了一个有趣的竞争。是阻碍界面的相干波的增长会胜出，还是湍流射流的混沌涡旋会在这些波增长之前将其破坏，从而可能推迟危机的到来？要回答这个问题，必须比较不稳定性增长的特征时间尺度与湍流涡旋的周转时间。这样的分析揭示，结果取决于流体性质和流动条件的微妙平衡，这显示了科学家如何通过探究多种物理机制的相互作用来不断完善他们的理解。
沸腾系统中的不稳定性并不仅限于单个表面。考虑沸腾液体流过一个受热通道，这是蒸汽发生器和核反应堆中常见的设置。在某些条件下，系统会开始自发振荡，流率、压力和蒸汽含量会在一个危险的、自我维持的循环中来回摆动。这些被称为密度波振荡（DWO）。
其机制是一个经典的反馈回路。入口流率的短暂下降意味着流体在受热通道中停留的时间更长。它吸收更多的热量，更多的液体变成蒸汽。这个低密度蒸汽的大气泡随后离开通道。由于通道的总压降部分取决于流体柱的重量，这团轻质蒸汽会导致所需驱动压力的瞬时下降，这反过来又将更多的流体吸入入口。这种过度修正导致流量激增，从而产生更少的蒸汽，循环往复。
预测这些振荡何时发生是一个关键的安全问题。最简单的模型，如均匀平衡模型（HEM），将两相混合物视为单一流体，其中蒸汽和液体完美地一起移动（），并且始终处于相同的饱和温度。在某些情况下，比如在细管中的高速流动，这是一个合理的近似。然而，在许多现实场景中，这种简单的图景会彻底失败。例如，在宽管中的较低流速下，轻的蒸汽气泡可以比重的液体上升得快得多（滑移比 ）。在其他情况下，特别是当液体进入时远低于沸点，热量需要有限的时间才能产生蒸汽，这意味着各相不处于热平衡状态。这些效应——滑移和非平衡——引入了额外的延迟并改变了系统的动态响应。像 HEM 这样忽略这些效应的模型可能会预测系统是稳定的，而实际上，这些微妙的相滞后可能恰恰提供了触发剧烈振荡所需的确切反馈。这说明了在对复杂系统建模时的一个关键教训：有时，魔鬼确实藏在相间物理的细节中。
我们的旅程已经穿越了工厂、发电厂和引擎。在我们的最后一站，让我们将目光从地球上的机器投向浩瀚星空。这似乎是一个巨大的飞跃，但支配我们饮料中气泡的相同基本定律，也塑造了我们银河系的结构。
在我们进行这一飞跃之前，有必要停下来思考一下为什么这些问题如此具有挑战性。在单相流中，动量、热量和质量输运之间存在着优美而强大的类比。如果你知道一个流动受到的摩擦力有多大，你通常可以非常准确地预测它传递热量或质量的效果如何。
然而，在湍流多相世界中，这些优雅的类比常常完全失效。考虑一个两相流通过多孔介质，比如水和油通过岩石，或者反应物通过填充式催化剂床。为什么在这里类比会失败？有两个深层次的原因。首先，热量和质量存在根本不同的输运路径。热量可以通过多孔基质的固体骨架传导，形成一个绕过停滞流体区域的“捷径”。而质量没有这样的捷径；它被限制在流体相中。其次，流动的几何结构本身就是一个由气、液、固界面构成的复杂、变化的织锦。这些界面在整个体积内充当热源或汇、质源或汇，而不仅仅是在边界上。这完全打破了单相类比所依赖的数学相似性。认识到简单模型在何处失效与知道它们在何处成功同样重要；正是在这些理解的前沿，新的科学得以诞生。
恒星之间的空间，即星际介质（ISM），并非空无一物。它是一种广阔、稀薄且极其动态的多相流体，是气体、等离子体和尘埃的混合物，被超新星爆发和强大的恒星风搅成一片湍流。在银河尺度上，能量在最大尺度上注入，并像在地球上的流动一样，逐级向下传递到越来越小的涡旋。
但这些能量去向何方？在不可压缩流体中，它最终会通过粘性在最小尺度上耗散为热量。在 ISM 的高度可压缩和稀薄的气体中，一种主要的耗散模式是通过一个由弱激波组成的网络。穿过这些激波的气体被加热，然后这些热气体将其能量辐射到太空中。这种辐射的特性由一个“冷却函数”描述，它告诉我们气体在给定温度下能多有效地冷却。
现代天体物理学中一个宏伟而优美的思想是，大部分 ISM 处于一种统计稳态，其中湍流级联提供的能量速率与辐射冷却的速率精确平衡。通过将某个尺度上的湍流特征时间尺度与由相同尺度激波产生的温度下气体的冷却时间尺度相等同，天体物理学家可以建立预测 ISM 热状态和结构的模型。他们可以估算出维持不同气体相所观测到的温度所需的能量耗散。
于是，我们的旅程回到了原点。我们用来设计地球上化学反应器的湍流能量级联和特征时间尺度的概念，同样被用来解码我们银河系的结构和演化。从发酵罐中气泡的微观之舞，到星际气云的宇宙芭蕾，湍流多相流是一个统一的主题，是物理学基本定律普适力量和广度的明证。