A fictitious-source formulation for direct numerical simulation of turbulent boundary layers in prescribed pressure gradients

高速飞行器和喷气发动机依赖于一层附着在其表面的薄而无形的空气层，即边界层。当空气快速流过机翼或涡轮叶片的金属表面时，由于摩擦作用，空气会减速，从而产生一个速度、压力和温度不断变化的复杂区域。当空气运动速度快到足以产生可压缩性时（即其密度随着加速或减速而发生变化），这些变化会变得更加剧烈。工程师和科学家需要准确了解当这层空气受到压力变化（例如飞机爬升或俯冲时）的推动时是如何表现的，因为这些条件可能导致气流从表面分离，从而导致升力损失或阻力激增。长久以来的挑战在于，为了在计算机模拟中创造特定的压力变化，研究人员通常必须物理性地改变虚拟风洞壁面的形状，通过使其弯曲来迫使空气加速或减速。这引入了一个复杂的问题：人们无法分辨空气是对压力变化本身做出的反应，还是对产生该变化的弯曲壁面形状做出的反应。

一个研究团队开发了一种新方法，可以在无需改变虚拟空间壁面形状的情况下模拟这些条件。他们并没有物理性地重塑定义域，而是引入了一个巧妙的数学技巧：他们在控制流动的方程中加入了不可见的源项。可以将这些源项想象成在空气的特定位置注入或移除微量的质量、动量和能量，模拟空气在移动过程中经过变化面积时的扩张或收缩效应。通过精确计算这些增加和减少的量，该团队能够强制边界层边缘的空气遵循精确的速度和压力历史，就像它流经一个变化的面积一样，而此时的壁面本身仍保持完全平坦。这种方法使他们能够将压力梯度效应与壁面形状效应分离，而这种分离在以往的高速可压缩流研究中是非常难以实现的。

研究人员使用两种不同的场景测试了这一方法。首先，他们观察了低速流动的空气，在这种情况下，可压缩性的影响极小。他们将这种新的“带有隐形力的平直壁面”法与传统的模拟进行了对比，在传统模拟中，流动的上边界是通过物理性地扩张或收缩来产生压力变化的。结果显示两者非常相似。平直壁面模拟所产生的壁面摩擦、空气速度以及边界层内部湍流，与传统方法表现出的高度一致，以至于在工程师最关注的区域内，两种方法几乎无法区分。唯一的差异出现在远离壁面的外层流区域，在那里，传统方法中的物理扩张自然地产生了与隐形力不同的运动模式。这证实了新方法可以在不需要弯曲几何形状的情况下，准确地复制压力驱动的空气行为。

随后，研究团队将该方法推向极限，模拟了速度达到两倍音速的空气。在这个状态下，物理现象更为复杂，且壁面形状与气流之间的相互作用至关重要。他们将使用隐形源项来模拟特定压力历史的平直壁面模拟，与空气实际流经弯曲壁面的模拟进行了对比。尽管在边界层边缘，两种情况下的空气都遵循相同的速度历史，但研究人员发现流动的行为有所不同。弯曲壁面引入了与空气必须遵循的路径曲率相关的微妙效应，这改变了边界层的生长以及壁面上的摩擦。这是一个关键发现：它证明了虽然新方法可以完美控制压力历史，但它并未抹除“弯曲路径会以平直路径所不具备的方式影响流动”这一物理现实。该方法成功地分离了压力效应和曲率效应，表明它们是可以在研究中独立探讨的不同影响因素。

为了准确理解这两个因素是如何相互作用的，研究人员对所作用的力量进行了详细分析。他们将流动的动量分解为各个组成部分，观察压力梯度、空气速度以及壁面曲率各自如何贡献于整体行为。他们发现，弯曲壁面的影响并非作用于流动的单一、简单的力。相反，曲率改变了整个边界层的结构，进而改变了压力梯度和空气速度对流体的影响方式。曲率带来的直接几何效应与这些间接变化相比，显得异常微小。这意味着，当空气流经弯曲表面时，壁面的形状改变了空气对压力的响应方式，形成了一种复杂的、耦合的关系，如果不使用该团队开发的这种方法，很难将其理清。

这项工作的意义在于它为研究高速空气动力学提供了一个受控环境。通过使用平直壁面和隐形源项，研究人员现在可以精确规定沿表面的压力变化，而不受通常会导致该变化的几何形状的影响。这使得进行更清晰的各种流动条件对比成为可能，有助于科学家在没有几何复杂性带来的“噪声”干扰下，理解湍流和可压缩性的基本物理原理。该方法并非要取代对真实弯曲表面的研究，而是提供了一个隔离特定变量的强大工具。它提供了一种询问：“如果压力这样变化，但壁面保持平直，会发生什么？”并能得到可靠答案的方法。这种能力对于设计更高效、更稳定的高速飞行器至关重要，因为它允许工程师预测空气在特定压力条件下会如何表现，并将其与飞行器本身的形状区分开来。

研究还强调了该方法在不同速度下的鲁棒性。无论空气运动速度较慢还是处于超音速状态，该方法都能成功生成所需的压力梯度并保持边界层的完整性。研究人员指出，虽然该方法适用于气流贴附在表面的附着流，但在涉及气流分离或变得不稳定的情况下，可能需要进行调整。然而，对于绝大多数高速空气动力学应用而言，该技术提供了一种可靠且准确的方法来模拟复杂的压力环境。通过消除对复杂几何设置的需求，它简化了计算过程，并减少了运行此类模拟所需的时间和资源，使得探索更广泛的条件变得更加容易。

最终，这项研究代表了科学家处理模拟湍流边界层问题的一种转变。他们不再试图通过构建物理曲线模型来模仿自然，而是找到了通过数学操纵来模仿曲线效应的方法。这使他们能够专注于压力与湍流之间相互作用的核心物理机制。研究结果证实，压力梯度和壁面曲率是影响流动的两个不同因素，并且通过正确的工具是可以将它们分离的。随着团队继续完善该方法，它有望成为研究可压缩湍流的标准工具，为理解高速度下的空气行为提供新的见解，并为更先进的空气动力学设计铺平道路。这种能够独立于壁面几何形状控制压力历史的能力，为理解空气如何在表面运动的基本机制开辟了新的可能性，而这一知识库对于下一代高速飞行的发展至关重要。