try ai
科普
编辑
分享
反馈
  • 高超声速流

高超声速流

SciencePedia玻尔百科
主要结论
  • 高超声速飞行的特点是存在薄激波层、粘性相互作用和高温真实气体效应等在较低速度下不会出现的极端现象。
  • 像牛顿撞击理论这样的简单模型提供了基础性的见解,而相似律则允许将实验数据按比例缩放用于飛行器设计。
  • 在极端温度下,空气会发生离解和电离,这要求我们理解化学非平衡及其对热传递和流体性质的影响。
  • 应对灾难性的热量是设计的核心挑战,可通过热防护系统(TPS)来解决,该系统利用了烧蚀和非催化表面等原理。

引言

速度超过五倍声速的飞行不仅仅是传统飞行的加速版,它是一场进入不同物理领域的旅程。这就是高超声速流的世界,一个对于从太空返回的大气再入飞行器到下一代高速全球运输等技术都至关重要的领域。在此范畴内,空气本身会发生变化,我们所熟悉的空气动力学定律被推向极限,需要一套全新的原理来理解和预测飞行器所承受的力和极端热量。低速飞行与高超声速前沿之间的这一知识鸿沟,正是现代航空航天工程面临的最大挑战所在。

本文将带领读者探索高超声速流的复杂景象。我们将首先探讨其核心的 ​​原理与机制​​,从简单的粒子模型到高温化学与边界层的复杂相互作用,逐层剖析其中的物理学。随后,我们将考察其 ​​应用与跨学科联系​​,揭示如何利用这些原理来设计和保护飞行器,并将流体力学与热力学、材料科学以及先进计算联系起来。

原理与机制

进入高超声速流的世界,就如同进入一个我们熟悉的飞行规则开始弯曲和失效的领域。这不仅是一个速度惊人的世界,更是一个能量巨大的世界,在这里,空气本身不再表现为我们日常经验中那种温和、连续的流体。在这里,我们将逐层揭开这个迷人领域的面紗,从一个极其简单的图像开始,逐步构建起支配着太空边缘飞行的复杂乃至奇异的物理学。

粒子图像:牛顿之梦

我们能想象到的物体以惊人速度在空气中飞行的最简单方式是什么?远在现代流体力学变得复杂之前,Isaac Newton 提出了一个极具直观性的优美模型。想象一下,空气不是一种连续的流体,而是大量微小、独立的粒子集合,就像一场宇宙冰雹。当一个物体冲破这场“冰雹”时,粒子撞击其表面,在一次完全非弹性碰撞中,沿着表面滑开。

在这个图像中,一个粒子将其全部的法向(垂直于表面)动量传递给表面,但其切向动量则完全不传递。我们在表面上感受到的压力就是这种持续轰击所产生的力。从这个简单假设中,我们得出了一个非常强大的结果:压力系数,一个衡量表面压力的无量纲数,由著名的​​正弦平方律​​给出:

Cp=2sin⁡2θC_p = 2 \sin^2\thetaCp​=2sin2θ

其中,θ\thetaθ 是表面相对于来流的局部夹角。如果流动平行于表面(θ=0\theta=0θ=0),则没有法向撞击,因此压力不增加。如果流动正面撞擊表面(θ=90∘\theta=90^\circθ=90∘),压力达到最大值。这个​​牛顿撞击理论​​为我们估算高超声速飞行器迎风面的巨大压力提供了一个出奇好的初步猜测。

你可能会认为这只是一个粗糙的类比。但在物理学中,简单的模型往往蕴含着深刻的真理。使用斜激波的气体动力学理论更为复杂,给出了一套复杂得多的方程。然而,如果我们把这些方程推向其物理极限——让马赫数 M∞→∞M_\infty \to \inftyM∞​→∞,并且关键地,让比热比 γ→1\gamma \to 1γ→1——它们会神奇地简化为牛顿的正弦平方律。γ→1\gamma \to 1γ→1 的极限代表一种无限可压缩的气体,可以被压缩成紧贴物体的无限薄层。在这个极限下,气体动力学图像就变成了牛頓的粒子图像。激波紧紧地贴着物体,以至于它们之间的空间消失了,流体的行为就像一股粒子流在表面上放弃了它们的动量。这两个看似不同的模型之间的统一是物理学中一个反复出现的主题,揭示了自然界内在的连贯性。

相似的力量

求解完整的高超声速流方程是一项艰巨的任务。但如果我们不必这样做呢?如果我们能找到一个“窍门”,一个相似性原理,告诉我们不同的流动情况在某种根本方式上是相同的,那会怎样?这就是​​高超声速相似律​​的魔力。

考虑一个细长体以非常高的马赫数 M∞M_\inftyM∞​ 和小攻角 α\alphaα 飞行。真正支配流动行为的是什么?与其说是总速度,不如说是垂直于物体表面的速度分量,也就是引起压缩的部分。这个法向速度大约是 U∞αU_\infty \alphaU∞​α。控制其物理过程的关键参数原来是这个法向速度与声速之比,一种“法向马赫数”。这便引出了​​高超声速相似参数​​:

K=M∞αK = M_\infty \alphaK=M∞​α

其深刻含义在于,对于同一个细长体,在不同马赫数 M1M_1M1​ 和 M2M_2M2​ 下的两种流动,如果它们的相似参数相同,那么它们的压力分布将几乎完全相同。也就是说,如果我们在马赫数为10、攻角为 2∘2^\circ2∘ 的条件下进行风洞试验,我们只需将攻角减半至 1∘1^\circ1∘,就可以预测马赫数为20时飞行的压力分布,使得 M1α1=M2α2M_1 \alpha_1 = M_2 \alpha_2M1​α1​=M2​α2​。这个原理是空气动力学设计的基石,它允许工程师利用縮比模型的实验数据来设计全尺寸飞行器。这种相似性概念甚至延伸到空气本身的复杂化学反应,其中必须保持另一个参数,即自由来流压力与物体长度的乘积(p∞Lp_\infty Lp∞​L),为常数,以确保化学反应在模型和全尺寸飞行器之间能够正确地缩放。

当粘性与无粘相遇:相互作用的边界层

在经典空气动力学中,我们通常将世界分为两部分:一个靠近表面、粘性起主导作用的薄而“粘滞”的​​边界层​​,以及一个广阔的、无摩擦的“无粘”外部流场。我们分别分析它们,然后将它们拼凑在一起。在高超声速流中,这种权力的和平分立被彻底打破。这两个区域展开了一场剧烈的对话,这一现象被称为​​高超声速粘性相互作用​​。

它的作用原理如下。对于一个以高超声速运动的物体,边界层内的粘性摩擦是巨大的。这种摩擦产生大量的热量,将靠近表面的气体温度提高到数千度。根据理想气体定律,这种高温高压的气体膨胀,密度变得小得多——它会“膨胀”起来。

这形成了一个出人意料的厚边界层,并且随着流动沿物体表面发展而迅速增长。从外面冷的超声速流的角度看,这个迅速增厚的边界层就像一个实体斜坡。外部流动被迫向上偏转,远离物体。当超声速流发生偏转时,会产生一道激波。在这种情况下,会形成一道弱斜激波,它不是由物体本身引起的,而是由粘性边界层形成的“等效体”形状所诱导的。

这道激波反过来又增加了外部流动的压力。这个更高的压力随后通过边界层传递到飛行器表面。结果是惊人的:一个以零攻角飞行的薄平板,无粘理论预测其表面压力不会有任何变化,但实际上在其前缘附近却承受着非常高的压力,这一切都源于流动中粘性部分和无粘部分之间这种错综复杂的反馈回路。

激波层的微妙秘密

让我们更仔细地看看激波层本身——即弓形激波与物体之间的热压缩气体区域。在一个钝头飞行器上,激波是弯曲的:在头部正前方是一道非常强的正激波,随着它绕过肩部弯曲,激波逐漸变弱,变得更斜。

热力学的一条基本定律告诉我们,激波越强,产生的熵就越多。这意味着激波层中的流体带有它穿过激波哪一部分的记忆。一个穿过头部强激波的流体粒子具有非常高的熵,而一个穿过远处较弱激波的粒子则具有低得多的熵。这就形成了一个分层的流动,熵存在着强烈的梯度。这个紧贴物体的高熵流体区域被称为​​熵层​​。

但故事并未就此结束。伟大的物理学家 Aurel Stodola 以及后来的 Luigi Crocco 发现了熵和涡量之间的深刻联系。Crocco定理实质上告诉我们,在流动中存在熵梯度的地方,必然会产生涡量(衡量局部涡旋或旋转的物理量)。因此,我们的熵层也是一个​​涡层​​。

随着粘性边界层沿物体表面增长,它可以“吞噬”或吸入这种高熵、高涡量的流体。当这种​​熵层吞噬​​发生时,边界层的游戏规则就改变了。它的外缘现在接触到的流体比正常情况下更热、密度更低、涡旋更强。这极大地改变了边界层内的温度分布,并且对飞行器设计至关重要的是,它可以显著增加到表面的热流率。这是一个由激波几何形状产生的微妙效应,但对飞行器的热防护系统而言,却可能带来生死攸关的后果。

是气体,但非完美:真实气体效应

到目前为止,我们大多假设空气是一种“完全气体”。但在高超声速激波后惊人的温度下——通常比太阳表面还要热——空气远非完全气体。构成空气的分子本身开始表现出奇怪的行为。这就是​​真实气体效应​​的领域。

像氮气(N2\text{N}_2N2​)或氧气(O2\text{O}_2O2​)这样的分子以多种方式储存能量:通过移动(平动)、翻滚(转动)以及原子像由弹簧连接一样来回振动(振动)。当一定体积的气体被激波剧烈压缩和加热时,能量并不是均匀分布的。平动和转动模式几乎立刻被激发——在几次分子碰撞内。分子立即开始更快地移动和翻滚。然而,振动模式更“僵硬”,需要更长的时间来吸收它们的能量份额。

这就产生了一种奇异的​​热力学非平衡​​状态。气体不能再用单一的温度来描述。相反,我们必须使用一个​​双温度模型​​:

  • 一个​​平动-转动温度(TtT_tTt​)​​,反映分子移动和翻滚的动能。这个温度在穿过激波时几乎瞬间跃升到一个非常高的值。
  • 一个​​振动温度(TvT_vTv​)​​,表征储存在振动中的能量。这个温度滞后,从低值开始,然后慢慢“弛豫”上升到 TtT_tTt​。

这带来了深远的影响。压力是分子与表面碰撞的结果,这是一个由其平动运动支配的过程。因此,压力由平动-转动温度(TtT_tTt​)决定。然而,气体的总内能,决定其热容和其他性质,是所有能量模式的总和,因此同时取决于 TtT_tTt​ 和 TvT_vTv​。在能量更高的情况下,振动会变得如此剧烈,以至于分子本身会分解(​​离解​​),最终,原子可能会被剥离电子(​​电离​​),使空气变成发光的等离子体。这些真实气体效应从根本上改变了物理过程,改变了激波脱体距离、压力分布,尤其是热传递。

超越连续介质:当空气变得“粗糙”

我们整个讨论都建立在一个最后的、巨大的假设之上:即空气是一种​​连续介质​​,一种光滑、连续的物质。这个假设是著名的​​Navier-Stokes方程​​的基础,后者是流体力学的支柱。但这个假设只有在分子两次碰撞之间行进的距离——即​​平均自由程​​(λ\lambdaλ)——远小于我们问题的特征长度尺度(LLL)时才有效。这两个长度的比值就是至关重要的​​Knudsen数​​,Kn=λ/LKn = \lambda/LKn=λ/L。

在高超声速飞行器通常飞行所在的极高海拔处,空气极其稀薄,因此平均自由程 λ\lambdaλ 很大。这增加了Knudsen数,空气开始看起来不像一种流体,而更像是一堆独立的粒子。连续介质假设开始失效。

但高超声速飞行还设下了另一个更隐蔽的陷阱。即使飞行器在总Knudsen数很小(例如,Kn≪1Kn \ll 1Kn≪1)的高度飞行,流动本身也可能产生​​连续介质假设失效​​的局部区域。激波和边界层非常薄,有时只有几毫米厚。在这些区域内,温度和速度的梯度是巨大的。局部特征长度尺度(LgL_gLg​,即物理性质发生显著变化的距离)变得与平均自由程相当。局部Knudsen数,Kng=λ/LgKn_g = \lambda/L_gKng​=λ/Lg​,可能会变得很大,就在流动的核心区域,Navier-Stokes方程不再有效。

那时会发生什么?我们已经到达了流体力学的前沿。物理学家和工程师们已经发展出更高阶的理论,如​​Burnett方程​​和​​正则化13矩(R13)方程​​,这些理论试图通过包含从动理学理论推导出的更复杂的项来修补连续介质模型。对于更稀薄的流动,他们完全放弃了连续介质模型,回归到基于粒子的视角,使用强大的计算方法,如​​直接模拟蒙特卡洛(DSMC)​​。在一个美妙的循环中,我们发现自己再次追踪单个粒子,正如牛顿所想象的那样,但现在我们掌握了几个世纪以来关于分子在天空边缘复杂舞蹈的知识。

应用与跨学科联系

在经历了高超声速流基本原理的旅程之后,我们现在来到了故事中激动人心的部分:见证这些原理的实际应用。高超声速飞行的研究并非某种抽象的学术活动;它是关于极端的科学,是建造能够勇敢面对大气进入的狂暴或在几分钟内跨越大陆的机器的艺术。这是一个我们那些整洁的低速理论被撕碎的领域,迫使我们锻造新的概念工具。在这里,我们将探讨这些工具如何被应用,揭示一个流体力学、热力学、化学、材料科学乃至计算理论必须联手的景象。

近似的艺术:驾驭力

远在超级计算机时代之前,航空航天工程的先驱们面临着一项艰巨的任务:如何预测一个以难以想象的速度运动的物体所承受的巨大作用力和压力?他们无法求解完整而复杂的方程。相反,他们发展出了优美简洁、直观的模型。其中最著名的就是​​牛顿撞击理论​​。

想象一下,迎面而来的空气不是连续的流体,而是一阵微小的、独立的粒子雨,就像雨点击打挡风玻璃一样。在这个模型的最简单版本中,当这些粒子撞击飞行器表面时,它们将全部的法向动量传递给它,就好像它们在一次完全非弹性碰撞中粘在了表面上一样。切向动量则被假定不受影响。从这个近乎卡通化的简单图像中,人们可以推导出对物体压力分布的一个非常有用的初步估计。任何一点的压力都简单地与当地表面与流动方向夹角的正弦值的平方成正比。这解释了为什么再入飞行器要使用钝头体:它们向流动展示了一个大的垂直表面积,产生一个强大但脱体的弓形激波,该激波与物体保持一定距离,有助于控制热负荷。我们甚至可以改进这个模型,想象粒子以其初始法向速度的一部分反弹回来,类似于一个具有特定“弹性”或恢复系数的球。这增加了一层真实感,显示了压力如何取决于气体-表面相互作用的性质。

当然,这个模型是一种 caricature(夸张的简化)。但它的力量在于其简洁性和它所提供的物理洞察力。其他巧妙的近似方法也随之而来。例如,​​切锥法​​是一个绝妙的技巧。它提出,绕一个尖锥体的复杂三维流动可以在任何一点上被一个更简单的二维流动所近似,即绕一个楔角的流动,该楔角与当地锥角相匹配。这使得工程师可以使用现成的楔形体流动解来估算锥体上的压力,而锥体是一种具有巨大实际重要性的形状。

随着我们理解的加深,我们开始窥探激波和物体之间的薄而炽热的区域——激波层。​​薄激波层理论​​提供了一个更精细的图像。其最优雅的预测之一是由离心力引起的修正。当被压缩的高密度流体尖啸着沿着钝体的曲面流动时,就像一列火车绕着弯道行驶。流体感受到一股将它拉离表面的离心力,这使得物体上的压力比简单的牛顿估计值略有降低 [@problem-id:469495]。这是一个美丽的例子,说明我们如何建立科学理解:我们从一个简单的模型开始,然后在其上叠加额外的物理效应以提高其准确性。

物理的统一:爆炸波与边界层

有时在物理学中,我们会发现看似毫不相干的现象之间存在惊人的联系。​​爆炸波相似理论​​是流体力学中最深刻的例子之一。由杰出的苏联科学家G. G. Chernyi提出,它将绕细长高超声速物体的定常流与强大爆炸产生的不 steady、膨胀的火球联系起来。

想象一下引爆一条长长的圆柱形炸药。一个圆柱形激波向外猛冲,其后的流动随时间演变。该相似理论指出,附着于高超声速细长体上的激波后的流场空间演化,与圆柱形爆炸波后的流动时间演化直接类似 [@problem_se:637582]。飞行器问题中的下游距离xxx扮演了爆炸问题中时间ttt的角色。供给爆炸波的能量类似于高超声速物体所受的阻力。这不仅仅是定性上的相似;这是一种深刻的数学等价关系,它使我们能够利用人们熟知的爆炸物理学来预测高速射弹周围的压力和激波形状。这是物理定律统一力量的明证。

另一个迷人的相互作用发生在飞行器的边缘。在低速流中,我们认为边界层是一层极薄的“粘性”空气膜。但在低密度、高空的高超声速飞行中,这一层会变得如此之厚,以至于从根本上改变了流动。这就是​​高超声速强相互作用​​的范畴。厚厚的边界层有效地“填充”了飞行器,向迎面而来的外部流动呈现出一个新的、更钝的形状。这个修正后的形状产生了自己的激波,进而增加了表面的压力。这个更高的压力又压缩了边界层,改变了其增长。这个错综复杂的反馈回路,即边界层决定压力场,压力场又决定边界层,支配着尖锐机翼和鰭片前缘的作用力和热流。

炽热的现实:热、化学与材料响应

高超声速飞行的最大挑战不是力,而是火。在马赫数10时,空气的驻点温度可以达到数千开尔文,比许多恒星的表面还要热。管理这种末日般的加热是飞行器设计的核心问题。

我们对热传递的低速直觉在这里完全失效。我们不能再认为热传递是由简单的温差驱动的。在高超聲速下,边界层内部的摩擦——一种称为​​粘性耗散​​的现象——成为热能的主要来源。能量方程必须重新构建。关键的通货不是温度,而是​​焓​​——气体的总能量含量,包括其热能和化学能。热传递的驱动势变成了边界层边缘的恢复焓与壁面气体焓之间的差值。专门的理论,如用于驻点加热的​​Fay–Riddell理论​​,就是专门为处理这种高焓环境而发展的,考虑了这个关键区域强烈的压力梯度和真实气体行为。

这种剧烈的加热会产生戏剧性的后果。其中最引人注目的是​​烧蚀​​,即表面通过燃烧消耗来散热的过程,就像流星划过夜空一样。本来会渗入飞行器结构的热量,反而被用作汽化热防护罩材料的潜熱。汽化产物随后被注入边界层,这个过程称为“吹除”,它进一步加厚了边界层,并阻挡了传入的对流热。这种牺牲性防御是几乎所有大气再入舱所使用的热防护系统(TPS)背后的原理。

此外,在这些温度下,空气本身也改变了它的特性。它不再是双原子氮和氧的简单混合物。分子间的剧烈碰撞使它们分解,这个过程称为​​离解​​。空气变成了一个由分子、原子、甚至离子和电子组成的反应性化学汤。这是一种​​真实气体效应​​。因为能量被用来打破这些化学键,与理想的惰性气体相比,真实的离解气体在激波后的温度会更低,密度会更高。这种化学活动不是瞬时发生的;它需要时间。化学反应发生所需的时间与流体流过飞行器所需的时间之间的竞争至关重要。这个比率由​​Damköhler数​​捕捉。如果流动太快以至于反应跟不上,我们就处于化学非平衡状态,这对作用在飛行器表面的压力和热传递有一阶效应。

飞行器表面本身也可以在这种化学作用中扮演积极的角色。有些材料是​​催化​​的,意味着它们充当化学“媒人”,促使边界层中的离解原子在表面直接重组成分子。这种重组将其原子的键能直接作为热量释放,从而大大增加了飞行器的热负荷。相比之下,非催化(或“被动”)表面则阻止这种重组,从而显著减少了加热。因此,TPS材料的选择不仅仅是关于其耐热能力,还关于其化学惰性。

数字风洞:计算与现代设计

高超声速流的极端复杂性——其相互作用的激波、粘性层、化学反应和表面烧蚀——远远超出了简单分析模型所能捕捉的范畴。现代高超声速飞行器的设计几乎完全依赖于​​计算流体力学(CFD)​​,它使用功能强大的计算机在数值网格上求解控制运动的方程。

CFD使我们能够模拟整个多物理场环境,考虑到流体动力学、气体中的化学动力学以及TPS材料的热和化学响应之间的耦合。例如,一个烧蚀热防护罩的模拟不仅必须包括空气的流动,还必须包括从壁面注入的热解气体、一个涉及数十种物质和数百个反应的完整气相化学反应机理,以及表面的催化特性。

然而,即使在这里,在纯数字领域,高超声速流的极端性质也带来了独特的挑战。其中最微妙和顽固的问题之一是“总能问题”。气体的总能量EtE_tEt​是其内能ρe\rho eρe(与温度有关)和动能12ρu2\frac{1}{2}\rho u^221​ρu2的总和。在高超声速流中,动能可能比内能大数千倍。具有有限精度的计算机计算总能量EtE_tEt​,然后试图通过减法计算出微小的内能来求得温度:ρe=Et−12ρu2\rho e = E_t - \frac{1}{2}\rho u^2ρe=Et​−21​ρu2。这就像试图通过称量整艘战舰(船长在船上和不在船上时)来称量船长的体重——微小的差异很容易在数值噪声中丢失。这种“相减抵消”可能导致温度 wildly inaccurate or even negative temperatures.

解决方案是一种数值上的巧妙设计。用于高超声速的稳健CFD代码基本上保留了两套账本。它们总是保守地更新总能量,这对于以正确的强度和速度捕捉激波至关重要。但它们也演化一个独立的、辅助的内能方程。然后,代码使用一个开关:在流动的大部分区域,它信任保守计算的总能量。但在它检测到动能压倒内能的区域,它会明智地切换到使用辅助内能方程的值来计算压力和温度。这种“双能量”方法确保了激波的物理特性得到完美保留,同时避免了在高速、较冷的流动部分出现灾难性的数值错误。这是深度计算科学如何被用来掌握物理世界最极端方面的完美例子。

从简单的粒子模型到庞大复杂的模拟,理解和驾驭高超声速流的旅程是一场持续的冒险,推动着工程、化学和计算的边界。这是一个要求整体观念的领域,物理世界的每一个组成部分都在一场激烈、炽热而最终美丽的舞蹈中交织在一起。