Temperature transformation recovering the compressible law of the wall for turbulent channel flow

本文针对可压缩湍流通道流中缺乏成熟温度变换的问题，基于动量与能量平衡方程提出了新的范德里斯特型和半局部型温度变换，并通过数值模拟验证了半局部型变换在粘性底层和缓冲层具有更优的数据坍缩效果及高精度，从而有效恢复了可压缩壁面律。

要点

这篇论文主要解决了一个流体力学中的“翻译”难题：如何把高速、可压缩（比如飞机在高空飞行时）的流体温度分布规律，翻译成大家熟悉的低速、不可压缩（比如水流过水管）的规律。

为了让你更容易理解，我们可以把这篇论文的研究内容想象成**“给不同语言的人翻译同一首诗歌”**。

1. 背景：为什么需要“翻译”？

想象一下，你有一本关于低速水流（不可压缩流）的“经典诗歌集”（这就是著名的壁面律 Law of the Wall）。在这本书里，水流的速度和温度分布非常有规律，就像诗歌的韵律一样，大家都能背下来，工程师们也能用它来设计管道或预测天气。

但是，当流体速度变得很快（比如超音速飞机），空气会被压缩、变热，就像**“方言”**一样，原本那本“经典诗歌集”里的规律就不灵了。速度和温度的关系变得乱糟糟的，工程师们很难直接套用旧公式。

以前的做法：
科学家们尝试过几种“翻译方法”（变换公式），试图把高速流的数据强行塞进低速流的框架里。

• 问题：有些翻译在“近墙”的地方（比如紧贴着墙壁的一层）翻译得很准，但到了“中间”或“外层”就乱了；有些翻译在处理“绝热墙”（不传热）和“等温墙”（恒温）混合的情况时，会出现奇怪的“尖刺”或错误。

2. 这篇论文做了什么？（新的翻译器）

作者（Youjie Xu 等人）发明了一套全新的“翻译器”（温度变换公式），包含两种版本：

1. VD 型（Van Driest 型）：像是一个老派的翻译，简单直接，但在某些复杂情况下不够精准。
2. SL 型（半局部型）：这是一个升级版的高级翻译。它不仅能处理简单的情况，还能在靠近墙壁的“粘性底层”和“缓冲层”里，把高速流的数据完美地“折叠”回低速流的规律上。

核心创新点（三个“修正系数”）：
为了让翻译更准确，作者引入了三个关键的“修正系数”（论文里叫 $\psi_1, \psi_2, \psi_3$），我们可以把它们想象成翻译时的**“微调旋钮”**：

• 旋钮 1（混合长度模型 $\psi_1$）：
  • 比喻：就像翻译诗歌时，要考虑到不同方言的“语调起伏”。以前大家假设语调是直线的，但作者发现，在管道中心，语调其实是抛物线形状的。调整这个旋钮，能让翻译在管道中心更自然。
• 旋钮 2（外力做功 $\psi_2$）：
  • 比喻：想象你在推一辆车。如果是平地推（等温墙），和上坡推（有外力驱动），车的速度分布不一样。这个旋钮负责计算“推力”对温度分布的影响，防止翻译时出现能量不平衡导致的“乱码”（尖刺）。
• 旋钮 3（湍流动能通量 $\psi_3$）：
  • 比喻：这是最难的“高深词汇”。湍流就像一锅沸腾的汤，里面有很多小漩涡在乱窜，带走热量。以前的翻译经常忽略这些“小漩涡”。作者发现，在等温墙的情况下，必须把这些“小漩涡”的热量搬运算进去，否则翻译出来的温度就会偏高或偏低。

3. 实验结果：翻译得有多好？

作者用超级计算机模拟了各种复杂情况（就像用各种方言念同一首诗），然后测试他们的“翻译器”：

• 效果惊人：在使用了**SL 型（半局部型）**变换后，原本杂乱无章的高速流温度数据，竟然完美地重合到了低速流的“经典诗歌”上！
• 误差极小：在大多数情况下，翻译的误差只有 2% 左右。这就像把一篇几千字的古文翻译成现代文，只有几个字意思稍微有点偏差，整体读起来非常通顺。
• 解决了“尖刺”问题：以前的公式在管道中间会出现奇怪的数值跳动（尖刺），新的公式通过引入“外力做功”和“湍流动能”的修正，把这些尖刺抚平了。

4. 简化版：让普通人也能用的“速查表”

虽然上面的公式很精确，但计算太复杂（需要知道很多高阶数据，比如那些“小漩涡”的具体数据）。
作者还做了一个**“简化版速查表”**（简化变换公式）：

• 他们发现，虽然“小漩涡”很重要，但可以用一个简单的数学模型（像是一个平滑的曲线函数）来近似代替。
• 这样，即使没有超级计算机提供的详细数据，工程师也能用这个简化版公式，得到非常接近的结果。

5. 这项研究有什么用？（未来的应用）

• 设计更高效的飞机和发动机：工程师可以用这套新公式，更准确地预测高速飞行时的温度分布，从而设计出更耐热的材料或更高效的冷却系统。
• 改进计算机模拟（LES）：在做大规模流体模拟时，计算机很难算到墙壁那么细。这个公式可以作为“墙壁模型”，告诉计算机“墙壁附近大概是什么样”，从而大大节省计算时间。
• 逆向工程：甚至可以反过来用。如果我们知道低速流的规律，想预测高速流会是什么样，这个公式也能帮上忙。

总结

这篇论文就像是一位精通多门方言的文学大师，他不仅发现以前翻译高速流温度规律的方法有“语病”（能量不平衡、忽略小漩涡），还发明了一套带有三个“微调旋钮”的新翻译法。

这套新法子能把高速、高温、高压的复杂流体数据，完美地“降维”成大家熟悉的低速流体规律，误差极小，而且还能简化成“速查表”供日常使用。这对于未来设计超音速飞行器、火箭发动机以及优化能源系统，都是一块非常重要的基石。

技术摘要

这是一篇关于可压缩湍流通道流中温度壁面律（Law of the Wall, LoW）重构的学术论文。文章提出了一种新的温度变换方法，旨在解决可压缩流动中温度分布难以像速度分布那样通过变换统一为不可压缩形式的问题。

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

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

• 背景：在高雷诺数不可压缩湍流边界层中，平均速度遵循著名的“壁面律”（对数律）。根据雷诺类比，温度分布也应遵循类似的规律。然而，在可压缩流动中，由于气动加热效应和速度 - 温度场的耦合，现有的温度分布难以直接应用不可压缩的壁面律。
• 现有挑战：
  • 现有的温度 - 速度关系（TV-relation）通常需要边界层外缘的速度和温度作为输入，这在内部流动（如通道流、管道流）中难以预先确定，且对复杂热边界条件（如混合等温/绝热壁）适应性差。
  • 现有的温度变换（如 Van Driest 型或半局部型）在缓冲层（buffer layer）和粘性底层（viscous sublayer）的数据坍缩效果不佳，或者在绝热壁条件下存在奇点问题。
  • 之前的变换往往忽略了混合长度模型的非线性、体积力的做功以及湍流动能（TKE）通量的影响，导致能量不平衡，进而引起变换后的温度剖面出现非物理的尖峰或偏差。

2. 方法论 (Methodology)

作者基于动量平衡和能量平衡方程，在重叠层（overlap layer）进行了详细推导，提出了新的Van Driest 型 (VD-type) 和 半局部型 (SL-type) 温度变换。

• 核心推导：
  • 从控制体积的动量和能量方程出发，考虑了分子扩散、湍流扩散、体积力做功以及湍流动能（TKE）通量（$q_t^k$）。
  • 引入了三个关键参数来修正变换核：
    1. $\psi_1$：与混合长度模型（Mixing Length Model, $l_m$）相关，用于修正粘性底层和重叠层的尺度。
    2. $\psi_2$：与体积力（Body Force）及其对流体做功相关，修正了总剪切应力分布。
    3. $\psi_3$：与湍流动能通量（TKE flux）相关，用于解决能量不平衡问题。
• 混合长度模型的选择：
  • 对比了线性模型（$l_m = \kappa y$）、抛物线模型（$l_m = \kappa y \sqrt{1-y/h}$）和增强型模型。研究发现，在粘性底层中，混合长度不应像传统涡粘模型那样被阻尼（damped），而应满足 $l_m \approx \kappa y$ 以保证线性律；在重叠层，抛物线或增强型模型能更好地捕捉雷诺应力。
• 简化模型：
  • 针对 $\psi_3$ 依赖高阶统计量（TKE 通量）的问题，作者分析了 $q_t^k/q_w$ 的多层结构，提出了基于双曲正切函数（tanh）的近似模型（$\beta_1, \beta_2, \beta_3$），从而导出了简化版温度变换，使其无需直接依赖高阶统计量即可应用。

3. 主要贡献 (Key Contributions)

1. 提出新的温度变换框架：首次系统地基于动量 - 能量平衡推导了包含混合长度、体积力做功和 TKE 通量修正的 VD 型和 SL 型温度变换公式。
2. 揭示能量不平衡机制：指出了以往变换（如 Huang et al. [8]）在分母中出现符号变化导致温度剖面出现尖峰（spikes）的根本原因是忽略了体积力做功和 TKE 通量导致的能量不平衡。新变换通过 $\psi_2$ 和 $\psi_3$ 有效消除了这些尖峰。
3. 发现 TKE 通量的多层结构：详细分析了可压缩通道流中湍流动能通量的分布，识别出其在粘性底层、缓冲层、对数层和尾流区的不同行为，并据此建立了实用的近似模型。
4. 提供简化变换方案：开发了不依赖高阶统计量的简化变换公式，使其更易于在工程应用（如 LES 壁面模型）中使用。

4. 研究结果 (Results)

研究使用了直接数值模拟（DNS）数据（Gerolymos & Vallet, Lusher & Coleman）和壁面解析大涡模拟（WRLES）数据进行了验证，涵盖了等温壁、绝热壁及混合热边界条件。

• 数据坍缩效果：
  • SL 型变换优于 VD 型：SL 型变换在粘性底层和缓冲层的数据坍缩效果显著优于 VD 型，能更好地恢复不可压缩温度剖面。
  • 对数律恢复：使用抛物线混合长度模型（$l_m^P$）或增强模型（$l_m^E$）时，变换后的温度在重叠层呈现出清晰的对数律特征。
  • 误差分析：对于等温壁，整个边界层的积分平均误差（Integral Mean Error）大多低于 2%，均方根误差约为 1.7%。在混合热边界条件下，等温侧的误差也控制在 3.5% 以内。
  • 绝热壁挑战：在绝热壁一侧，变换后的温度幅值仍略高于参考值（误差 >10%），这可能与混合长度在近壁区的过预测或参考不可压缩剖面的对称性假设有关，这是未来改进的方向。
• 简化模型性能：基于 $\beta$ 模型的简化变换在大多数情况下表现良好，特别是在高雷诺数下，能够准确预测不可压缩参考剖面。

5. 意义与应用 (Significance)

• 理论意义：该工作深化了对可压缩湍流中热 - 动量耦合机制的理解，特别是明确了 TKE 通量和体积力在能量平衡中的关键作用，为建立普适的可压缩壁面律提供了坚实的理论基础。
• 工程应用：
  • 近壁建模：提出的变换可用于改进可压缩流动的壁面模型（Wall Models），特别是基于对数律的代数壁面模型和控制理论方法。
  • 逆变换：可用于从不可压缩参考剖面重构可压缩流动的平均速度、温度、密度等关键物理量，辅助设计。
  • 通用性：该方法不仅适用于通道流，还通过附录推广到了管道流，并讨论了向零压力梯度边界层扩展的潜力。

总结：这篇论文通过严谨的推导和广泛的数值验证，提出了一套能够有效恢复可压缩湍流温度壁面律的新变换方法。它解决了以往变换中存在的能量不平衡和缓冲层表现不佳的问题，并为可压缩流动的数值模拟和工程建模提供了强有力的工具。