try ai
科普
编辑
分享
反馈
  • 流体动力学中的膨胀耗散

流体动力学中的膨胀耗散

SciencePedia玻尔百科
核心要点
  • 流体运动可以分解为保持体积的涡旋(螺线运动)和改变体积的挤压(膨胀运动),每种运动都对应一种粘性耗散模式。
  • 膨胀耗散是由于流体压缩和膨胀的摩擦而导致的动能不可逆损失,由体积粘度决定。
  • 膨胀耗散的重要性与湍流马赫数的平方(Mt2M_t^2Mt2​)成正比,使其成为高速可压缩流中的一个关键因素。
  • 标准湍流模型需要明确的可压缩性修正来考虑膨胀耗散,这对于准确预测航空航天应用中的阻力、热量和混合至关重要。

引言

在高速流体动力学的世界里,空气可以像弹簧一样被压缩,一种微妙而强大的现象主导着能量的流动:膨胀耗散。尽管在低速场景中常被忽视,但理解这一过程对于准确预测超音速射流、再入飞行器乃至宇宙现象的行为至关重要。许多为不可压缩流体开发的传统工程模型未能考虑这种能量损失机制,可能导致危险的设计缺陷。本文通过将膨胀耗散分解为其核心组成部分来揭开其神秘面纱。第一部分“原理与机制”将深入探讨其基本物理学,区分流体的“挤压”运动与“涡旋”运动,并解释每种运动如何导致能量损失。随后,“应用与跨学科联系”部分将展示这一概念的深远影响,从改进航空航天工程的计算模型到解释天体物理激波的动力学。

原理与机制

要真正掌握高速飞行的本质,我们必须深入了解运动中流体的内部运作。想象一种既能像星系一样涡旋,又能像弹簧一样压缩的物质。这种双重特性是关键所在。膨胀耗散的故事,就是流体如何为其被压缩付出代价的故事,这个代价以损失的能量和产生的热量来衡量。

流体运动的两种灵魂:涡旋与挤压

让我们从一幅简单的画面开始。流体的任何复杂运动都可以被看作是两种基本运动类型的组合。第一种是涡旋、翻滚、旋转的运动——想象一下将奶油搅入咖啡时形成的美丽涡流。这是一种保持体积的运动;流体单元被剪切和旋转,但没有被挤压。用物理学的语言来说,我们称之为流动的​​螺线​​分量,因为它在数学上是无散的。

第二种运动是挤压和拉伸,膨胀和压缩——想象一下推动装有空气的注射器活塞。这种运动改变了流体单元的体积。我们称之为​​膨胀​​分量,因为它与密度的变化相关。该分量在数学上是无旋的,意味着它缺乏螺线部分那种局部旋转的特性。

著名的​​亥姆霍兹分解​​(Helmholtz decomposition)定理告诉我们,任何速度场都可以严格且唯一地分解为这两个部分:一个螺线(涡旋)的灵魂和一个膨胀(挤压)的灵魂。对于我们视为不可压缩的流体,比如厨房水槽里的水,情况很简单。流体就像一袋不可挤压的弹珠;它只能涡旋。它的膨胀灵魂处于休眠状态。但对于可压缩流体,比如流过超音速飞机机翼的空气,两种灵魂都完全苏醒,并以复杂的舞蹈形式相互作用。

运动的代价:粘性耗散

运动不是没有代价的。当你搅拌咖啡时,涡旋不会永远持续下去;它会逐渐平息。运动的能量,即动能,被转化为热量,使咖啡略微变暖。这个过程被称为​​粘性耗散​​,它是流体内部摩擦的结果。

如果流体运动有两个灵魂,那么它必须为这两种运动付出两种不同的代价。事实证明,粘性耗散也可以完美地分为两部分,每种运动各占其一。单位体积的总能量耗散率,我们称之为 Φ\PhiΦ,是这两个贡献的总和:

Φ=2μSijSij+ζθ2\Phi = 2\mu S_{ij}S_{ij} + \zeta \theta^2Φ=2μSij​Sij​+ζθ2

我们不必被这些符号吓倒。这个方程讲述了一个优美的物理故事。

第一项 2μSijSij2\mu S_{ij}S_{ij}2μSij​Sij​ 是​​螺线耗散​​。它代表了因涡旋和剪切运动的摩擦而损失的能量。量 SijS_{ij}Sij​ 是应变率张量,它衡量流体单元变形的速度,而 μ\muμ 是我们熟悉的​​动力粘度​​——这个属性衡量流体对剪切的阻力,就像蜂蜜和水的区别一样。这是不可压缩流体中存在的唯一耗散形式。

第二项 ζθ2\zeta \theta^2ζθ2 是我们今天的主角:​​膨胀耗散​​。这是由于体积变化的摩擦而损失的能量。量 θ\thetaθ 是膨胀率,代表流体在某一点上膨胀或压缩的速率。而 ζ\zetaζ (zeta) 是一个不那么出名但同样重要的流体属性,称为​​体积粘度​​。它衡量流体对被压缩或膨胀的内部阻力。动力粘度 μ\muμ 关乎对形状变化的阻力,而体积粘度 ζ\zetaζ 关乎对尺寸变化的阻力。除非流体正在被主动压缩或膨胀(θ≠0\theta \neq 0θ=0),否则这一项为零。

可压缩交响曲:湍流中的能量交换

当流动变为湍流时,它是一个由各种尺寸的涡流组成的混沌漩涡。在可压缩流中,这种混沌不仅包括涡旋,还包括波动的压缩和膨胀区域,就像一个狂乱、无序版本的注射器活塞。在这种剧烈环境中,我们必须区分涉及膨胀的两种截然不同的过程。

首先,我们有湍流版本的膨胀耗散,通常表示为 εd\varepsilon_dεd​。这是一个真正的耗散性、不可逆的过程。它是粘性摩擦中专门针对压缩运动的分量,将其动能直接并永久地转化为热量。像所有形式的摩擦一样,这是一条由热力学第二定律决定的单行道。它始终是动能的一个汇。

其次,我们有一个更微妙、更有趣的过程,称为​​压力-膨胀相关性​​,Πd=p′θ′‾\Pi_d = \overline{p' \theta'}Πd​=p′θ′​。这一项代表了脉动压力场(p′p'p′)对脉动膨胀(θ′\theta'θ′)所做的功。与耗散不同,这个过程原则上是可逆的。想象一个高压区压缩一个流体单元;这将动能转化为内能。然而,该流体块随后可能移动到一个低压区并膨胀,获得一个“反冲”,将内能转换回动能。

打个比方,膨胀耗散就像一个漏水的桶——水(能量)永远地流失了。压力-膨胀就像水在两个相连的桶之间来回晃动,一个代表动能,另一个代表内能。从长远来看,取决于压力和膨胀脉动的相关方式,它可能成为湍流动能的净源或净汇。在许多具有激波类结构的高速流中,高压与强压缩相关,使得压力-膨胀成为一个显著的动能净汇。

马赫数的指挥棒:引导能量流

是什么决定了涡旋运动和挤压运动之间的平衡?是什么决定了有多少能量被引导到膨胀模式中进行耗散?这场交响乐的指挥家是一个称为​​湍流马赫数​​(MtM_tMt​)的无量纲量。它定义为湍流脉动的特征速度与当地声速之比:Mt=2k/aM_t = \sqrt{2k}/aMt​=2k​/a,其中 kkk 是湍流动能,aaa 是声速。

不要将此与飞机的飞行马赫数混淆。湍流马赫数描述的是湍流本身的可压缩性。

在低 MtM_tMt​ 时,湍流表现得好像是不可压缩的。涡旋占主导地位,挤压可以忽略不计。但随着 MtM_tMt​ 的增加,一件非凡的事情发生了。流动的非线性动力学,主要通过压力场,开始将能量从强大的螺线模式转移到膨胀模式。混沌的涡旋实际上开始产生声波和压缩脉动。

结果,驻留在膨胀分量中的动能比例 kdk^dkd 增长。因此,膨胀耗散 εd\varepsilon_dεd​ 成为湍流能量越来越重要的一个汇。严谨的理论分析和高保真度计算机模拟揭示了一个优美简洁的标度律:膨胀耗散率相对于螺线耗散率,与湍流马赫数的平方成正比。

εd∝εsMt2\varepsilon_d \propto \varepsilon_s M_t^2εd​∝εs​Mt2​

这种二次方标度律是现代湍流理论的基石。它告诉我们,当湍流的内部马赫数加倍时,膨胀耗散的重要性将增加四倍。

工程师的困境:模拟不可见之物

所有这些讨论不仅仅是学术性的;它对工程学具有深远的实际影响。为了设计超音速飞机或再入舱,工程师们依赖于​​计算流体动力学(CFD)​​,它使用湍流模型来预测阻力和热量等问题。

湍流的主力模型,如标准的 ​​k-ε 模型​​,是为不可压缩流开发和校准的。它们的“DNA”中只包含有关螺线耗散的信息。它们在结构上对膨胀耗散和压力-膨胀的存在是盲目的。

当这些模型被用于具有显著 MtM_tMt​ 的高速流动时,它们会感觉到能量消失得比预期的要快。由于不知道膨胀耗散这个独立的通道,它们犯了一个关键错误:它们归咎于它们所知的唯一机制——螺线耗散,并人为地夸大其值。这就像一个修理工通过过度加速发动机来修理汽车的爆胎。其结果是对湍流结构的根本性错误预测,通常导致对湍流混合和耗散的过高预测,这可能转化为不准确的设计。

解决方案源于我们刚刚讨论的物理学,即引入一个​​可压缩性修正​​。工程师们明确地在他们的模型中添加一个新项,以解释膨胀耗散。其中最著名的修正之一由 Sarkar 提出,其形式为:

εtotal=εs+C⋅εs⋅Mt2\varepsilon_{\text{total}} = \varepsilon_s + C \cdot \varepsilon_s \cdot M_t^2εtotal​=εs​+C⋅εs​⋅Mt2​

这里,εs\varepsilon_sεs​ 是基准模型的标准螺线耗散,第二项是 εd\varepsilon_dεd​ 的模型,完美地捕捉了我们发现的 Mt2M_t^2Mt2​ 标度律。这是一个将基础物理学直接转化为工程工具的美丽范例,帮助我们建造更安全、更高效的高速飞行器。

当挤压无关紧要时:Morkovin 假设

为了完成我们的旅程,我们必须问最后一个关键问题:只要飞行速度高,可压缩性对湍流结构就总是重要的吗?令人惊讶的答案是否定的。

我们必须再次区分飞行马赫数 M∞M_{\infty}M∞​ 和湍流马赫数 MtM_tMt​。一架飞机可能以 5 马赫的速度飞行,但在紧贴其机翼的边界层深处,湍流涡旋的脉动速度可能远低于当地声速,导致 MtM_tMt​ 很小。

这一见解由 Mark Morkovin 在其著名的 ​​Morkovin 假设​​中正式提出。该假设指出,如果湍流马赫数 MtM_tMt​ 很小(通常小于约 0.3),那么可压缩性对湍流结构的直接影响可以忽略不计。在这种状态下,膨胀耗散和压力-膨胀项的量级为 Mt2M_t^2Mt2​,可以安全地忽略。湍流在其核心上表现得好像是不可压缩的。

那么,在这些高速流动中,可压缩性的主要影响是什么呢?是由于气动加热导致的平均流体性质——尤其是密度和粘度——的巨大变化。为了处理这个问题,建模者使用一种称为 ​​Favre(或密度加权)平均​​ 的技术,它优雅地简化了方程。Morkovin 假设的惊人结果是,一旦我们使用这种平均技术,我们就可以将我们信任的不可压缩湍流模型应用于广泛的高超音速流动,例如预测到飞行器表面的热传递。

这使我们的故事形成了一个完整的闭环。我们从将运动分为涡旋和挤压开始。我们了解到每一种运动都有其代价,一种耗散形式。我们看到了湍流马赫数如何指挥它们之间的交响乐,以及工程师们必须如何教他们的模型去聆听这首音乐。最后,我们了解到,有时,即使在高超音速风的咆哮中,湍流本身也以不可压缩的旋律低语,而挤压则淡入沉寂。

应用与跨学科联系

在探讨了膨胀耗散的原理和机制之后,我们现在踏上一段旅程,去看看这个迷人的概念在实践中的应用。就像一把万能钥匙,它解锁了我们对跨越惊人尺度范围的现象的理解,从尖啸着掠过高超音速机翼的稀薄空气,到在星系间涟漪般传播的寂静而巨大的激波。我们将看到,这不仅仅是专家们使用的一个深奥的修正项,而是一条美丽而统一的线索,编织在流体动力学的结构之中。

机器之心:铸就飞行的未来

膨胀耗散最直接、最深刻的应用或许是在高速飞行领域。当飞机以超音速冲破空气时,低速流动的简单、平滑模式让位于由激波、膨胀和剧烈湍流组成的漩涡。为了设计能够承受这些条件的飞行器,我们必须能够预测它们,而我们的故事就从这里开始。

想象一个超音速混合层,这是一个湍流的战场,两股不同速度的气流在这里发生冲突。在不可压缩的低速流中,这个混合层通过贪婪地卷入周围的空气而增长,这是一个由湍流涡驱动的过程。但随着马赫数的增加,一件非凡的事情发生了:混合层的增长变得慢得多。湍流似乎变得迟缓,效率降低。为什么呢?一个关键原因就是膨胀耗散。湍流涡内部的快速压缩和膨胀开始从湍流中消耗能量,起到强大的制动作用。在我们的计算机模型中,包含一个膨胀耗散项使我们能够捕捉到这种对混合的真实物理抑制,这对于设计从超音速燃烧室到喷气发动机降噪系统的一切都至关重要。

当我们考虑整个飞机上的流动时,故事变得更加关键。想象一下流动从表面分离,就像空气在陡峭倾斜的襟翼或后向台阶后脱离——这是航空航天设计中的一个常见特征。这个流动在下游重新附着到表面的点决定了压力、力和结构载荷。早期的湍流模型,对可压缩性的制动效应视而不见,通常会预测这种再附着发生得太早。通过考虑膨胀耗散,我们的模拟给出了一个更真实、更长的再附着长度,更好地匹配我们在实验中观察到的情况,从而实现更安全、更高效的设计。

现在,让我们将风险提高到大气再入的残酷环境中。返回地球的太空舱或高超音速飞行器会经受极端的加热。预测这种气动热载荷是生死攸关的问题。在这里,膨胀耗散的作用变得至关重要。在这些流动中,会形成强大的激波,与它们相互作用的湍流被强烈压缩。如果没有一个合适的膨胀耗散模型,模拟可能会极大地过高预测激波下游的湍流水平。这种“过剩”的湍流反过来又会导致对传递到飞行器表面的热量的灾难性过高预测。通过正确模拟压缩如何从湍流中消耗能量,膨胀耗散项有效地“驯服”了模拟中的湍流,使壁面摩擦和热通量的预测与现实更加吻合。它帮助我们避免了模拟中可怕的(且不符合物理的)“热峰值”,并设计出有效的热防护系统。

建模者的艺术:将物理学构建于机器之中

那么,我们如何在一台计算机模拟中捕捉这种复杂的物理现象呢?在这里,我们窥探计算物理学家的思想。这是一种由严谨原则指导的艺术形式。

第一个原则是所有原则中最基本的:能量守恒。当膨胀耗散从湍流速度脉动中移除能量时,那部分能量并不会凭空消失。它必须去往某个地方。去哪里?它被直接转化为热能,加热流体。因此,我们添加到湍流动能方程 kkk 中以模拟这种耗散的任何“汇”项,都必须在流体内能方程中显示为一个大小相等、符号相反的“源”项。这种优美而简单的平衡确保了我们的模拟尊重热力学第一定律,将湍流统计的抽象世界与温度这一具体属性联系起来。

但物理学的内容更加丰富。湍流不仅仅是一团随机的能量;它具有结构和方向性,这一特性被称为各向异性。想象一下被激波挤压的湍流。它不仅仅是失去能量;它还会改变形状。脉动可能在一个方向上被压扁,而在另一个方向上被拉长。为了捕捉这一点,物理学家们使用了更先进的工具,如雷诺应力模型,该模型追踪湍流应力张量 RijR_{ij}Rij​ 每个分量的演化。在这个复杂的框架内,可压缩性以两种不同的方式表现出来。首先是​​膨胀耗散​​,这是一个各向同性的过程,从所有分量中消耗能量,代表湍流“死亡”为热量。其次是​​压力-膨胀​​项,它描述了流体的平均压缩或膨胀如何直接改变湍流涡的形状,在不同分量之间传递能量。区分这两种效应对于准确预测像激波-边界层相互作用这样复杂的现象至关重要,在这些现象中,整个湍流结构被剧烈地重组。

这些模型项的创建本身就是一段发现之旅。通常,我们没有一个从第一性原理出发的完美理论。相反,我们运用物理推理。例如,为了模拟可压缩性如何影响湍流尺度方程(即 ε\varepsilonε 方程),我们可以通过施加约束来推导出修正项的形式:它必须具有正确的物理量纲,它必须在不可压缩极限下消失,并且其效应在极端压缩下应达到饱和。这种推理路线导出了优雅的数学形式,这些形式可以通过与实验和更详细的模拟进行对比来测试和校准,从而使我们能够系统地构建更准确、更稳健的模型。

这种对物理的深刻理解甚至可以指导计算过程本身。在一种称为目标导向的网格自适应技术中,我们可以告诉计算机将其计算资源集中在最重要物理现象发生的区域。如果我们想准确计算一个流动中的总膨胀耗散,还有什么比膨胀耗散本身更好的指导呢?我们可以使用其大小作为指标,告诉模拟在强压缩区域自动放置更多的网格点,确保我们以高保真度捕捉这些关键事件。从某种意义上说,是物理学告诉计算机如何看待它自己。

回声与低语:天体物理学与声学

膨胀耗散的影响远远超出了我们的地球机器,延伸到宇宙的寂静深处和声音的微妙传播中。

考虑一个穿过一片湍流空气的声波。这个波是压缩与膨胀的舞蹈。当它穿过湍流时,湍流涡被反复地挤压和拉伸。这个过程激活了膨胀耗散,导致涡流从声波中提取能量,并将其转化为热量。结果是声波被阻尼;它的能量被湍流实实在在地“吃掉”了。从宏观角度看,湍流赋予了流体一种有效的“体积粘度”,这是一个量化其对体积变化阻力的属性。这一现象是气动声学的基石,解释了湍流如何在从喷气发动机排气到大气本身的各种事物中帮助减弱声音。

现在,让我们把目光从可闻之声转向天文现象。星际介质并非宁静的虚空;它是一个充满暴力的地方,被来自超新星爆发和恒星风的巨大激波纵横交错。激波前缘是一个极其薄的区域,在这里,气体的宏观动能被剧烈且不可逆地转化为热能,将气体从数千度加热到数百万度。这种转变的引擎是粘性。在激波内部,流体以惊人的速率被压缩。这个巨大的速度散度 ∇⋅u\nabla \cdot \boldsymbol{u}∇⋅u 驱动了强大的粘性耗散,成为主要的加热机制。这个过程,在本质上,是膨胀耗散的一次戏剧性而宏伟的展示,它提供了使超新星遗迹发光数千年的热能。

一个基本问题:斯托克斯的幽灵

在这次宏大的巡游之后,很自然地会问:这个已被证明如此多才多艺的效应,其最终来源是什么?要找到答案,我们必须回到流体力学的基础——流体中应力与应变的关系。

对于一个简单的牛顿流体,我们将粘性应力与流体变形率联系起来。这种关系包含两个常数:控制对形状变化阻力的剪切粘度 μ\muμ,以及控制对体积变化阻力的体积粘度 ζ\zetaζ。我们一直在讨论的膨胀耗散与这个体积粘度 ζ\zetaζ 直接成正比。

然而,在19世纪,乔治·斯托克斯爵士提出了一个方便的命题。基于对单原子气体的论证,他假设这两种粘度是相关的,即 ζ=0\zeta=0ζ=0。一个多世纪以来,​​斯托克斯假设​​一直是流体动力学的基石,被内置于教给每个学生的标准纳维-斯托克斯方程中。根据这个假设,对于大量问题,膨胀耗散从一开始就被定义为不存在!

那么我们为什么花了这一章来讨论它呢?因为在湍流和高能等离子体的复杂世界里,斯托克斯假设失效了。一个湍流流动,从宏观上看,其行为并不像一个简单的气体。它会产生一个有效的体积粘度,即使其底层的流体分子遵守斯托克斯定律。无数涡旋与压缩应变相互作用的集体行为,产生了平均流所感受到的对体积变化的阻力。激波内部的等离子体也是如此。

因此,膨胀耗散是机器中的幽灵。它是一个方便、简化的假设失效的物理体现。它提醒我们,宇宙的丰富性常常在我们最简单的模型失效的地方被发现,并且通过探索这些前沿,我们揭示了对世界更深刻、更统一的理解。