玻尔百科在模拟复杂流体运动的探索中，计算流体动力学（CFD）依赖数值方法将物理学的控制定律转化为预测。其挑战在于创造出的方法不仅要准确，还要足够稳健，以处理流体可能展现的各种现象——从微风拂面到剧烈的激波。在最成功的方法中，平流上游分裂格式（Advection Upstream Splitting Method, AUSM）族格式因其坚实的物理直觉基础而备受赞誉。AUSM并非将控制方程视为一个纯粹的数学问题，而是优雅地将其解构为基本的物理分量。
本文深入探讨AUSM+格式的原理和威力。我们将探究它所解决的知识空白：需要一种通用的求解器，既能高保真地捕捉激波和接触面等尖锐特征，又能在极大的速度范围内保持稳定和准确。通过理解其设计，您将看到物理洞察力如何带来卓越的数值性能。
首先，在“原理与机制”一章中，我们将剖析通量分裂的核心思想，揭示流体中的信息流如何被分离为对流分量和压力驱动分量。随后，“应用与跨学科联系”一章将展示该方法非凡的多功能性，演示同样的基本逻辑如何让我们能够模拟从飞机机翼、火箭发动机到遥远黑洞的天体物理射流等一切事物。
要真正领会像平流上游分裂格式（AUSM）这样一种方法的精妙之处，我们不能只看最终的方程。我们必须追溯其发现之旅，从它试图描述的物理现象本身开始。想象我们正在观察一条河流。我们同时看到两件事在发生：水本身在向下游流动，携带着树叶和嫩枝；同时，涟漪在水面上传播，携带着关于扰动的信息。AUSM的核心思想是认识到流体运动定律包含这两种截然不同的现象，并构建一个能分别处理这两种现象的数值方法。
无粘、可压缩流体（即摩擦可忽略不计的气体或液体）的运动由一套优美的原理——欧拉方程所支配。在其守恒形式下，欧拉方程指出，一个体积内某个量的变化率等于该量穿过该体积边界的净流量，即通量。对于一维流动，这可以写为 。
在这里， 是流体的状态，一个包含密度（）、动量（）和总能（）的向量。向量 是通量，描述了这些量的输运。它看起来是这样的：
乍一看，这似乎是一堆杂乱的项。但借助一些物理直觉，我们可以将其理清。AUSM方法的精妙之处在于，它看到这个通量实际上是两个不同物理过程的总和。
首先，是流体运动并携带其自身属性的简单行为。这就是对流。如果流体以速度 运动，它就携带着自身的质量、动量和能量。我们可以将这个对流通量写成速度 乘以一个包含密度、动量和总焓 的向量。总焓是随流体自然对流的量：
其次，是压力 的效应。压力是一种内力；它是流体微团相互推挤的方式。这个压力通量只包含压力力，仅影响动量方程：
如果将这两部分相加，，就能完美地恢复原始的欧拉通量 。这种分解是AUSM在概念上的核心。这不仅仅是一个数学技巧，而是基于物理现象的分离。 描述了流体的整体输运，而 则代表了压力力。
为什么这种物理分裂如此强大？因为它完美地反映了信息在流体中传播的两种方式。欧拉方程是双曲型的，这意味着信息以有限的速度作为波进行传播。当我们分析其数学性质时，会发现恰好存在两种类型的波。
平流模式：这种波以当地流体速度 传播。它携带温度和密度（在恒定压力下）等属性的变化，但不携带压力信号。可以把它想象成风中的一缕烟；它只是随气流漂移。该模式直接对应于我们的对流通量 。
声学模式：这些是声波。它们相对于流体以声速 传播。从静止的观察点看，它们的移动速度为 （下游）和 （上游）。这些波是压力变化在流体中传播的机制。在某一点的突然推动会产生一个向外传播的压力波。该模式直接对应于我们的压力通量 。
通过将通量分裂为对流和压力部分，我们已将数值问题分解为两个更简单、物理上截然不同的部分。现在我们可以设计策略来独立处理流体属性的“漂移”和压力波的“传播”。
在数值模拟中，世界被分解成离散的单元。为了计算流体如何变化，我们需要确定两个可能具有不同状态（左侧L和右侧R）的单元之间的界面通量。决定使用哪一侧的信息的过程称为迎风。
AUSM格式使用马赫数 作为主控制器——一根指挥棒——来协调这个过程。界面通量通过构建界面质量通量 和界面压力 来建立。
对于对流部分，逻辑很简单。质量通量应由“迎风”状态决定。AUSM族格式使用马赫数分裂函数 和 来构建它。对于从左到右的流动，我们期望界面马赫数 由左侧状态决定，反之亦然。组合后的界面马赫数为 。
对于压力部分，逻辑遵循声波。
这种分离处理正是该方法如此稳健的原因。物理学中仅涉及漂移的部分通过简单的质量迎风处理，而涉及复杂波相互作用的部分则通过一种复杂的、感知马赫数的压力混合来处理。
最初的AUSM是一个绝妙的想法，但用于亚声速混合的特定数学函数被发现存在微小缺陷。在非常具有挑战性的情况下，如近停滞流或非常强的激波处，它们可能会导致微小的不准确性或振荡。
AUSM+ 的发展是一个精炼的过程，其指导思想是追求数学上的优雅和物理上的完美。设计者们问道：这些混合函数必须具备的理想属性是什么？他们推断，从亚声速到超声速的过渡应该是完全平滑的。这意味着不仅函数的值，而且其斜率在音速点（）也必须匹配。这被称为确保 连续性。
当这些简单而优雅的平滑性约束被施加于所需最小阶数的多项式上时，一个唯一的解便应运而生。例如，对于马赫数分裂，亚声速区域（）中 分量的多项式为 。这个公式并非任意的；它是要求一个行为良好的格式所带来的数学结果。
此外，AUSM+引入了一个关键的新特性：目标压力耗散。激波是一个剧烈的、极薄的区域，流体属性在此发生跳跃。在数值上，这可能引起振荡。AUSM+增加了一个微小的、“智能的”耗散项，该项与跨越单元界面的压力跳跃成正比。这个项就像一个微型减震器，仅在激波处激活，以保持其稳定和清晰。关键是，这种耗散在其他任何地方都会消失，因此它不会模糊或“涂抹”流动的其他精细特征。这种智能设计优雅地解决了诸如“奇偶解耦”等棘手的数值问题，这是一种可能困扰激波模拟的棋盘状不稳定性。
这种对物理原理的细致关注带来了显著的实际好处，使AUSM族格式在众多方法中脱颖而出。
其最受称赞的成功之一是能够以极高的清晰度捕捉接触间断。接触间断是一个界面，在该界面上压力和速度是恒定的，但密度和温度发生跳跃（想象一层氦气坐在一层空气上）。因为AUSM的压力通量是由压力差驱动的，所以在接触面上它看不到任何作用对象，不会引入人为的扩散。与此同时，对流通量只是以正确的速度输运密度跳跃。结果是一个完美的、不被涂抹的清晰接触面，这是许多其他格式难以实现的。
另一个深远的优势是其“全速”能力。许多可压缩流求解器在低马赫数（）时会变得非常不准确。它们固有的数值耗散通常以声速 来衡量，而声速可能比流体速度 大数千倍。这就像试图用一台为卡车设计的秤来称量一根羽毛的重量——测量值淹没在噪声中。AUSM族格式如AUSM+及其后继者（AUSM+-up, SLAU2）巧妙地使其耗散与马赫数本身成比例。随着流动减慢，数值耗散自动减小，使得该格式能够对从通风房间中空气的轻微流动到航天器的高超声速再入等各种情况都保持其准确性。
最终，即使是最复杂的格式也必须遵守一个基本的速度限制以保持稳定。要运行一个模拟，所选的时间步长 必须足够小，以使移动最快的波不会在一步之内跳过整个网格单元。这就是Courant–Friedrichs–Lewy（CFL）条件。对于欧拉方程，最快的波以 的速度传播。因此，时间步长必须满足 。
总之，AUSM族格式是一个由物理洞察力驱动的算法设计的优美典范。通过将流体运动的复杂舞蹈解构为其基本步骤——对流和压力功——并通过马赫数来协调它们的数值处理，它在极其广泛的物理体系中实现了准确性、稳健性和效率的罕见结合。它证明了这样一个理念：最强大的工具往往是那些最仔细倾听其所描述物理现象的工具。
当我们学习物理理论的原理时，我们正在学习游戏的规则。质量、动量和能量守恒——这些是流体运动的基石规则。但知道规则是一回事，玩转游戏则是另一回事。宇宙以惊人的复杂性进行着这场游戏，创造出从火焰的精巧舞姿到恒星爆炸的巨大狂怒等一切事物。要在我们的计算机中捕捉这种复杂性，我们需要的不仅仅是对规则的字面翻译；我们需要一个具有物理直觉的工具，一个尊重游戏特性和精妙之处的工具。
平流上游分裂格式（AUSM）就是这样一种工具。我们已经看到了它背后的原理——将信息流优雅地分裂为对流部分（仅仅被携带的东西）和压力部分（宣告变化的声波）的思想。现在，让我们看看这种物理洞察力让我们能做什么。让我们探索这所开启的广阔应用领域，从设计安静的飞机到研究以接近光速运动的宇宙射流。
在画出一幅杰作之前，你必须能够画出一条干净的线条。在流体动力学的世界里，最基本的“线条”之一就是接触间断。想象一下一团热空气和一团冷空气之间的边界，两者一起运动。没有压力差，没有风切变，只有一个温度和密度的急剧变化。这听起来很简单，但对计算机来说，这可能是一场噩梦。许多数值格式要么会将这条清晰的线模糊成一团，要么更糟的是，凭空制造出虚假的压力摆动，就像在没有投石的池塘里泛起的涟漪。
然而，一个真正稳健的方法必须看清其本质：简单的平流。物质只是在移动。AUSM+格式，由于其本身的设计，就是为了出色地通过这个测试而设计的。它分离压力和对流信息的方式，使其能够在理想条件下以手术般的精度输运这种“接触间断”，保持边界完美清晰，压力场平稳安静。这不仅仅是一个学术上的好奇心；它是信任的基础。如果一个格式能正确处理这种简单、干净的情况，当我们让它去处理现实世界问题的绚丽混乱时，我们就能更加自信。
真实世界不是一个简单的一维管道；它充满了复杂形状的物体。为了模拟空气流过飞机机翼或汽车，我们必须建立该几何形状的虚拟表示。我们不能使用简单的矩形网格；我们必须使用由单元——三角形、四边形、四面体——组成的灵活网格，这种网格可以贴合任何任意形状。这被称为非结构网格。
像AUSM+这样的格式的美妙之处在于其逻辑是局域的。它在分隔两个单元的每一个面上进行操作。它只需要知道两侧流体的状态和面的方向。这使得它具有极强的适应性。我们可以将相同的基本分裂逻辑应用于我们复杂网格中的任何一个面，从而使我们能够将方法从简单的一维问题扩展到对最复杂工程设计的完全三维模拟。
但是，模拟不仅仅是中间的物体；它还关乎周围的世界。在我们的计算域的边缘会发生什么？如果我们在模拟一个喷气发动机，我们必须指定进气口和排气口的条件。这就是边界条件的艺术。在这里，AUSM+的物理性质再次提供了明确的指导。流体动力学中的特征线理论精确地告诉我们哪些信息流入一个域，哪些信息流出。例如，对于亚声速的喷气排气，外部世界可以施加其压力，但排出气体的速度和温度是由内部决定的。AUSM+中的分裂自然地反映了这种物理现实。通量的分裂部分对应于这些传入和传出的波，精确地告诉我们如何“插入”我们的边界数据而不违反物理学，确保我们的虚拟世界与真实世界正确连接。
流体的世界跨越了令人难以置信的速度范围。同样的方程既支配着通风系统中空气的温和流动，也支配着重返大气层的航天器上的高超声速流动。一个真正强大的工具必须是全速域的大师。
最具挑战性的区域之一是跨声速飞行，此时飞机的速度非常接近声速。在这里，机翼上的部分流动可能是亚声速的，而其他部分则变为超声速，产生激波。这种混合环境是出了名的不稳定。控制方程的数学特性发生变化，没有经过精心设计的格式可能会彻底失败。AUSM族格式固有的稳健性，能够保持一致水平的物理耗散，使其成为在这些险恶的跨声速水域中航行的宝贵工具，从而能够设计出高效安全的飞机。
也许更令人惊讶的是，另一个极端——非常低的速度——也带来了其独特的挑战。许多为高速、可压缩流设计的方法在马赫数接近零时变得不准确和低效。就好像在只需要低语时它们却在大喊大叫。像AUSM+-up这样的特殊变体就是为了解决这个问题而开发的。通过引入仅在低速时才激活的经过仔细校准的项，该格式能够优雅地适应，以与捕捉激波相同的保真度捕捉近乎不可压缩流的物理现象。这种“全速”能力意味着同样的代码可以模拟火箭发射，以及发射后飘散的袅袅青烟。
一个基本思想的真正力量，取决于它能达到的广度。通量分裂的概念并不仅限于纯粹的空气动力学；它提供了一种语言来描述相互作用的物理过程的交响乐。
燃烧：在喷气发动机或火箭内部，我们拥有的不仅仅是流动的空气。我们拥有燃料、氧化剂和产物的湍流、炽热的混合物，所有这些都在经历快速的化学反应。为了模拟这一点，我们不仅必须跟踪流的动量和能量，还必须跟踪每种化学物质的浓度。AUSM+的平流核心为此提供了一种自然的方式。就像我们输运动量一样，我们可以输运氢、氧和水的质量分数，使我们能够模拟我们称之为火焰的化学与流体动力学的复杂舞蹈。
多相与多组分流：自然界和工业界的许多过程都涉及不同流体的混合物——液体中的气泡、空气中的燃料液滴，甚至具有不同性质的气体混合物。在这些组分之间的界面上，对数值格式的简单应用可能会再次产生人为的压力振荡。通过调整AUSM+框架，例如使用每个组分的体积分数信息，我们可以设计出正确处理这些界面的格式，使我们能够以更高的精度模拟从燃料喷射器到气泡反应器的一切事物。
激波物理与天体物理学：当一个强激波撞击两种不同流体之间的界面时，一种被称为Richtmyer-Meshkov不稳定性的美丽而复杂的不稳定性就会出现。这个过程在超新星的爆炸中以及在试图利用核聚变的尝试中至关重要。模拟此类现象将数值方法推向了极限。不同的格式，如AUSM+、Roe或HLLC，各有其优缺点，产生不同的数值误差。在这些极端情景下对它们进行比较，有助于科学家了解其模拟的可靠性，并为工作选择最佳工具，而这项工作可能就是理解宇宙中重元素的诞生。
相对论流体动力学：现在来到最宏伟的舞台。以守恒原理表达的流体动力学定律是如此基本，以至于它们甚至延伸到了Einstein的狭义相对论领域。在宇宙中，有从黑洞附近喷射出的等离子体射流，以接近光速的速度运动。为了描述这些，我们需要相对论流体动力学。令人惊奇的是，AUSM+的核心思想——基于信息流方向分裂通量——可以被转换到这个奇异的领域。通过定义一个“相对论马赫数”并确保任何信号的传播速度都不允许超过光速，我们可以构建一个相对论AUSM格式。帮助我们设计飞机机翼的相同思维，也帮助我们模拟物质在宇宙所能提供的最极端环境中的行为。
从一条干净的线到燃烧的火焰，从一架飞机到一个黑洞，融入AUSM+格式的物理推理线索让我们能够探索一系列惊人的现象。它证明了这样一个理念：一个建立在深刻物理洞察力之上的工具，不仅是一个更好的工具，更是一把解锁更广阔理解宇宙的钥匙。