Characterisation of rough-wall drag in compressible turbulent boundary layers

✨

这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

Each language version is independently generated for its own context, not a direct translation.

这篇论文主要研究了一个非常实际的问题：当飞机或高速飞行器在空气中高速飞行（甚至达到超音速）时，如果表面变得粗糙（比如因为积冰、沙尘撞击或材料剥落），会产生多大的额外阻力？

为了让你更容易理解，我们可以把这篇论文的研究过程想象成**“给高速飞行器表面粗糙度‘量体裁衣’"**的过程。

1. 核心问题：为什么“量尺”在高速下不准了？

想象一下，你有一把专门用来测量普通速度下（比如骑自行车）粗糙表面阻力的“尺子”（科学家称之为Nikuradse 公式）。

在低速（不可压缩流）时：这把尺子很准。只要你知道表面有多粗糙（比如贴了多粗的砂纸），就能算出它会增加多少阻力。
在高速（可压缩流，如超音速）时：空气变得像“弹簧”一样，会被压缩。当气流流过粗糙表面时，不仅会产生摩擦，还会产生激波（就像船头划开水产生的波浪，但这里是空气被挤压产生的冲击波）。这会产生额外的“波阻”。

现在的困境是：科学家发现，如果直接拿低速时测得的“粗糙度数值”去套用高速公式，结果就不准了。就像你拿着给自行车设计的轮胎花纹去跑 F1 赛车，完全对不上号。

2. 他们做了什么实验？

研究团队在南安普顿大学和德国慕尼黑联邦国防军大学，利用一个巨大的风洞（就像一个超级大的吹风机），进行了以下实验：

模拟粗糙表面：他们贴了两种不同粗细的砂纸（P60 和 P24，就像家里打磨家具用的砂纸，P24 比 P60 更粗糙）。
模拟不同速度：他们让气流以不同的速度流过砂纸，从亚音速（像普通飞机）到超音速（2.9 倍音速，像高超音速导弹）。
测量数据：他们用激光（PIV 技术）像拍高速摄影一样，捕捉空气流过砂纸时的流动情况，计算阻力。

3. 他们发现了什么？（用比喻解释）

他们发现，在高速下，粗糙表面产生的阻力确实比低速时大，而且这种增加跟速度（马赫数）有关。

这就好比：

低速时：粗糙表面就像在平滑的冰面上撒了一把沙子，阻力增加是固定的。
高速时：粗糙表面就像在冰面上撒了沙子，而且沙子还在制造小海啸（激波）。速度越快，海啸越大，阻力就越大。

之前的理论认为，只要把空气的“压缩性”修正一下，低速的公式就能直接用在高速上。但实验证明，这个修正还不够，因为粗糙表面产生的“激波”是之前公式没考虑到的。

4. 他们提出了什么解决方案？

为了解决这个问题，作者尝试了三种不同的“修正尺子”的方法，试图把高速数据强行拉回到低速的公式曲线上：

方法一：找“替身”
- 试图在低速实验中找到一个跟高速实验“长得像”的粗糙表面，直接借用它的阻力数据。
- 结果：不太成功。就像试图用一只成年猫的重量去估算一只老虎的重量，虽然都是猫科动物，但体型差异太大，很难直接对应。
方法二：调整“粘度”
- 考虑到高速下空气的粘性和温度变化，尝试用空气在粗糙表面顶端的粘度来重新定义粗糙度。
- 结果：对作者自己的实验数据有效，但对别人的旧数据效果一般。
方法三：给阻力“打折”（最成功的方法）
- 作者发现，如果把计算出来的阻力值，乘以一个跟温度比（空气温度 vs 墙壁温度）有关的修正系数，所有的数据（包括他们自己的和以前别人的）都能完美地落在同一条曲线上！
- 比喻：这就像是你发现高速下的阻力读数“虚高”了，于是你拿了一个“折扣券”（修正系数），把读数打折后，就发现它其实和低速时的规律是一脉相承的。

5. 这意味着什么？（总结与展望）

目前的成就：作者找到了一种经验公式（虽然还没完全从物理原理上解释透，但在数学上很管用），可以让工程师们在设计超音速飞行器时，利用现有的低速粗糙度数据，通过简单的修正，估算出高速下的阻力。
未来的挑战：
- 目前的修正系数是“试”出来的（经验性的），不是从第一性原理推导出来的。
- 实验中没有直接测量墙壁的温度和摩擦力（因为风洞条件限制），而是靠估算。
- 未来的方向：需要开发一种专门针对“粗糙表面”的新的物理变换公式，不仅要考虑空气被压缩，还要考虑粗糙颗粒产生的“激波”效应。

一句话总结：
这篇论文告诉我们，在超音速飞行中，表面的粗糙度会产生额外的“激波阻力”，不能简单套用低速公式。作者通过大量实验，找到了一把新的“修正尺子”（基于温度比的修正系数），能更准确地预测高速飞行时的阻力，为未来设计更耐用的超音速飞行器提供了重要参考。

Each language version is independently generated for its own context, not a direct translation.

这是一篇关于可压缩湍流边界层（TBL）中粗糙壁面阻力特性表征的学术论文详细技术总结。该研究由 D. D. Wangsawijaya 等人完成，旨在解决如何将不可压缩流中建立的粗糙壁面阻力表征框架（基于 Nikuradse 和 Hama 的理论）有效推广到可压缩流（特别是超音速）中的问题。

以下是该论文的详细技术总结：

1. 研究背景与问题 (Problem)

背景：在不可压缩流中，粗糙壁面的阻力通常通过动量亏损 $\Delta U^+$ 和等效沙粒粗糙度 $k_s$ 来表征，遵循 Nikuradse (1933) 的对数律关系。在可压缩流中，通常先应用速度变换（如 Van Driest 变换）将可压缩速度场“拉伸”为不可压缩等效形式，然后计算 $\Delta U^+$ 。
核心问题：
1. 如果在不可压缩流中已知某种粗糙表面的 $k_s$ ，在什么条件下该值可以直接用于可压缩流？
2. 现有的基于光滑壁面推导的压缩性变换（Transformations）是否适用于粗糙壁面？
3. 可压缩流中，粗糙元产生的激波（波阻）如何影响阻力表征？现有的经验公式（如 Nikuradse 公式）在可压缩条件下是否仍然适用，或者需要修正？
现状局限：以往的可压缩粗糙壁面研究大多仅在单一马赫数（ $M$ ）和雷诺数（$Re $）下进行，导致$ k_s $被强制拟合到不可压缩的对数律关系中，缺乏对$ M $和$ Re$ 独立变化的系统性研究。

2. 方法论 (Methodology)

实验设施：利用德国慕尼黑联邦国防军大学的间歇式吹除式三音速风洞（TWM）。
实验对象：
- 粗糙表面：两种不同目数的砂纸（P60 和 P24），通过光学轮廓仪扫描获得了详细的表面形貌参数（如物理高度 $k$ 、均方根高度等）。
- 测试范围：马赫数 $0.3 \le M \le 2.9$ （涵盖亚音速、跨音速和超音速）；摩擦雷诺数 $7427 \le Re_\tau \le 30292$ 。
- 独立变量控制：在 $M=2$ 时，通过改变总压独立改变 $Re_\tau$ ，实现了 $M$ 和 $Re$ 的独立扫描。
测量技术：
- 使用粒子图像测速仪（PIV）测量流场速度剖面。
- 温度估算：由于设施限制未直接测量壁温，基于绝热壁假设，利用 Crocco-Busemann 关系式估算温度分布，进而推导密度和粘度分布。
- 壁面剪切应力（WSS）估算：未直接测量 WSS，采用改进的 Clauser 拟合法，结合壁面位置修正（针对粗糙元高度和 PIV 光反射）来反推摩擦速度 $U_\tau$ 和动量亏损 $\Delta U^+$ 。
数据处理：
- 应用了 5 种不同的速度变换函数（Van Driest, Howarth, Brun, Trettel & Larsson, Volpiani 等）将可压缩速度剖面转换为不可压缩等效形式。
- 对比分析了三种修正方案：
  1. 基于不可压等效的 $k_s$ 。
  2. 基于自由流粘度的粗糙度缩放 $k^* = k / \nu_\infty^+$ 。
  3. 基于壁面剪切应力修正因子 $\sqrt{1/F_c}$ 对 $\Delta U^+$ 进行修正。

3. 关键贡献 (Key Contributions)

填补数据空白：提供了高雷诺数（ $Re_\tau$ 高达 30,000）下，独立变化马赫数和雷诺数的粗糙壁面 TBL 实验数据，这是以往文献中缺乏的。
验证变换的适用性：发现对于粗糙壁面， $\Delta U^+$ 对不同的速度变换函数（如 VD, TL, VO 等）相对不敏感，但在完全粗糙区，对数律的截距表现出明显的马赫数依赖性。
提出并验证修正因子：
- 否定了简单使用不可压 $k_s$ 直接映射到可压流的方法。
- 提出了基于自由流与壁面温度比（ $T_\infty/T_w$ ）的修正因子 $\sqrt{1/F_c}$ 。该因子源自光滑壁面压缩性修正，但在此被创新性地应用于粗糙壁面的动量亏损修正。
揭示波阻机制：通过纹影图像（Schlieren）和速度散度分析，直观展示了粗糙元前方产生的弓形激波（Bow shocks）及其延伸至自由流的现象，证实了可压缩流中额外的“波阻”是现有光滑壁面变换无法完全解释的原因。

4. 主要结果 (Results)

完全粗糙区特性：所有测试案例均处于完全粗糙区（ $\Delta U^+ > 7$ ）。
马赫数效应：在亚音速和跨音速（ $M \le 0.8$ ）下，数据能较好地坍缩到 Nikuradse 的不可压对数律上。但在超音速区（ $M=2, 2.9$ ）， $\Delta U^+$ 的对数律截距随马赫数增加而显著上移，表明存在额外的阻力机制（波阻）。
修正因子效果对比：
- 方案 A（等效 $k_s$ ）：尝试通过匹配不可压流中的几何参数来寻找等效 $k_s$ ，结果并不理想，无法统一不同马赫数下的数据。
- 方案 B（粘度缩放 $k^*$ ）：使用自由流粘度 $\nu_\infty$ 对粗糙度高度进行缩放，仅对当前实验数据有效，对历史数据（不同粗糙类型）的普适性较差。
- 方案 C（ $\sqrt{1/F_c}$ 修正）：将 $\Delta U^+$ 乘以 $\sqrt{1/F_c}$ （其中 $F_c$ 是光滑壁面摩擦系数的压缩性修正因子，依赖于 $T_\infty/T_w$ ）。这是表现最好的方案。它成功地将当前实验数据以及大量历史文献数据（涵盖砂纸、立方体、正弦波等多种粗糙类型，从超音速到高超音速）统一坍缩到 Nikuradse 的对数律关系内（误差在 $\pm 10\%$ 以内）。
物理图像：Schlieren 图像清晰显示，粗糙壁面在超音速下会产生一系列弓形激波，这些激波导致额外的波阻，而现有的基于光滑壁面的变换未能完全消除这种效应。

5. 意义与展望 (Significance & Future Work)

工程意义：该研究为高速飞行器（如高超音速飞行器）表面的热防护、烧蚀、积冰或尘埃撞击导致的表面粗糙化带来的阻力预测提供了更准确的经验修正方法。
理论启示：
- 证明了现有的光滑壁面压缩性变换在粗糙壁面应用中存在局限性，特别是忽略了粗糙元诱导的激波效应。
- 提出的 $\sqrt{1/F_c}$ 修正虽然有效，但本质上是经验性的，且存在“双重变换”的逻辑（先变换速度再修正阻力）。
未来工作：
- 需要开发专门针对粗糙壁面 TBL 的自定义变换函数，该函数应能同时考虑粗糙度属性、壁面条件（绝热/等温）以及粗糙元产生的波阻。
- 未来的实验应致力于直接测量粗糙壁面的壁面剪切应力和壁温，以减少估算带来的不确定性，特别是在非均匀、不规则粗糙表面（如烧蚀表面）上。

总结：本文通过系统的实验和数据分析，指出在可压缩粗糙壁面流动中，简单的不可压 $k_s$ 映射失效，而引入基于温度比的阻力修正因子 $\sqrt{1/F_c}$ 是目前统一不同马赫数和粗糙类型数据的最有效经验方法，同时也指出了未来建立物理机制更完善的粗糙壁面压缩性变换的必要性。