玻尔百科流体的运动，从我们呼吸的空气到覆盖地球的海洋，都受基本物理定律的支配。几十年来，科学家们依赖一种被称为静力近似的强大简化方法来模拟这些庞大的系统。这种假设重力与压力之间存在完美平衡的方法，对于大尺度、缓慢演变的流动非常有效，构成了传统天气和气候模型的基石。然而，自然界常常远非平衡；它充满了剧烈的上升气流、破碎的海浪和湍流涡旋，在这些现象中，垂直加速度不仅存在，而且至关重要。这就产生了一个关键的知识空白：我们如何精确模拟这些动态的小尺度现象？
本文深入探讨了非静力模型的世界，这些工具旨在通过接纳流体运动的全部复杂性来填补这一空白。通过在控制方程中保留垂直加速度项，这些模型揭示了我们星球流体动力学更真实、更详细的面貌。接下来的章节将引导您了解这个复杂而迷人的主题。首先，“原理与机制”将揭示区分非静力模型与静力模型的根本物理原理，探讨使用它们的关键标准以及这一选择所带来的深远的计算后果。随后，“应用与跨学科联系”将展示这些模型的实际应用，揭示它们如何彻底改变我们理解和预测从单个雷暴、山地波到飓风强度的各种现象的能力。
要理解非静力模型的世界，我们必须首先深入探究我们如何描述运动中的流体。其核心在于，流体的运动，无论是大气中的空气还是海洋中的水，都遵循我们在初级物理课上学到的一个原理：牛顿第二定律，。对于一小团流体，该定律告诉我们其加速度是作用在其上的合力的结果。在垂直方向上，这个故事主要由两股主要力量之间的巨大拉锯战主导：向下拉拽流体团的重力，以及由流体团上方和下方压力差异产生的压力梯度力，后者可以向上或向下推动流体团。
想象一下广阔而宁静的太平洋。在数千公里的范围内，海水以缓慢而宏伟的环流形式运动。如果你跟踪一团水，你会发现它的垂直加速度小得惊人，几乎为零。在这种大尺度、缓慢演变的系统中，物理学家和海洋学家做出了一个绝妙的简化。他们问道：如果我们假设垂直加速度完全为零会怎样？如果 ，那么 。这场拉锯战变成了一种完美的静态平衡。来自压力梯度的向上推力恰好抵消了向下的重力拉力。这种美妙的平衡被称为静力平衡。
在这种静力近似下，压力变成了一个极其简单的量。任何给定深度的压力仅由其上方流体柱的重量决定。这意味着，如果你知道各处流体的密度，你只需从表面向下积分就可以计算出压力。这种简化是静力模型的基石，几十年来，这些模型一直是大规模气候和天气预测的主力。它们计算效率高，并且对于水平尺度远大于垂直尺度的现象来说，其准确性惊人 [@4072640]。
但自然界并不总是如此平静和广阔。想一想穿透大气层的高耸积雨云——即雷暴。在这里，暖湿空气团并非处于温和的平衡状态；它们正以猛烈的加速度向上飞窜。再想一想在海滩上翻腾欲碎的波浪，或一股寒冷咸涩的密集水流翻过水下山脉坠入深渊。在这些情况下，垂直加速度不仅仅是一个次要细节，它就是故事的全部。忽略它就等于完全错失了重点。
这就是非静力模型的用武之地。它们不做这种简化假设。它们完全接纳垂直动量方程的全部动态特性，将其视为一个真正的垂直速度预报方程。在这个世界里，加速度不为零，力的平衡是一场动态的、不断变化的戏剧。这种复杂性所带来的回报是能够模拟真实世界，展现其所有湍流、对流和类似波的辉煌景象 [@3802927]。
我们如何知道何时可以采用简单的静力世界，何时又必须冒险进入更复杂的非静力领域？答案，正如物理学中常见的那样，在于尺度的比较。
最直观的衡量标准是现象的纵横比：其高度除以宽度。对于一个宽达数千公里但仅有几公里深的广阔海盆来说，纵横比 非常小。它就像一张非常非常薄的纸。运动绝大多数是水平的，静力近似完美成立。但对于一个大约十公里宽、十公里高的雷暴来说，纵横比接近于一。它的高度和宽度相当。毫不奇怪，这里的垂直运动与水平运动同样重要 [@4072640]。
我们可以通过一个称为内部弗劳德数（）的无量纲数来更精确地描述。你可以将 看作一个比率：流体的特征速度（）与内重力波的自然传播速度（）之比，内重力波是在海洋或大气等分层流体内的密度层上传播的波。所以， [@4053471]。
让我们来看一个具体的例子。在深对流上升气流的核心区，特征垂直速度可达 ，深度为 。在典型的大气中，这可能对应于一个约为 的内部弗劳德数 [@4053471]。这个计算从第一性原理告诉我们，垂直加速度与驱动它的浮力大小处于同一量级。一个假设这些加速度为零的静力模型，将从根本上无视雷暴的核心物理过程。同样的原理也适用于流经陡峭水下海槛的密度流（此时流动可能变为超临界流并形成水跃），或是在海滩上浅化的表面波（此时纵横比 达到 1 的量级）[@3802927]。
放弃靜力近似最深远的影响不仅仅是我们得到了一个动态的垂直速度方程，而是压力的本质本身发生了转变。在非静力世界中，压力扮演了一个新的、至关重要的角色：它成为质量守恒的伟大执行者。
在大多数海洋和许多大气情景中，我们可以将流体视为基本上是不可压缩的。对于低速流（低马赫数）来说，这是一个非常好的近似，并在海洋的玻色尼斯克近似或大气的非弹性近似中得到形式化 [@3802879]。不可压缩性约束优雅而简单：。它表明速度场的散度处处为零。简单来说，流入一个微小盒子的流体不能比流出的多。流体不能在某一点上被凭空创造或毁灭。
在静力模型中，这个约束被用于在确定水平速度之后，诊断性地计算垂直速度 [@3815102]。但在非静力模型中，这个严格的数学约束是如何在每一点、每一瞬间都得到满足的呢？答案是压力。
想象一下模型在一个时间步长的计算。它计算了所有的力——平流、旋转、浮力——并确定了一个“临时”速度场。这个临时场是流体想要做的运动。但这个运动可能不滿足 ；它可能存在流体堆积或变薄的区域。为了纠正这一点，模型必须找到一个压力场，其梯度将提供恰到好处的推和拉，以调整速度场，使其无散度。
这就引出了计算流体动力学中最重要的方程之一：压力泊松方程。它的形式为 ，其中 代表所有非压力力 [@3813973]。这个方程是一个数学奇迹。方程右侧的项是临时速度场的散度——它衡量了“临时”流场在多大程度上违反了不可压缩性。该方程规定，压力场的曲率（）必须对这种违反做出响应。通过求解这个方程，我们找到了恢复质量守恒的唯一压力场 。
想一想这意味着什么。静力压力是一个局部属性，仅由正上方流体的重量决定。而非静力压力是一个全局属性。域中某一点的压力会瞬时受到其他所有地方流场运动的影响，因为约束条件 必须在全局范围内成立。压力就像一只无形的手，在整个流体中提供瞬时通信，以确保其协调一致。这与静力世界中压力简单的诊断性角色相比，是一个深刻的转变。
这种更真实、更优美的物理图像是有代价的。非静力模型的强大功能需要巨大的计算资源。
求解三维压力泊松方程是非静力模拟中计算成本最高的部分。它是一个巨大的线性方程组，将域中的每个网格点与所有其他网格点耦合在一起。这种“全局”特性使其成为并行计算中一个臭名昭著的瓶颈 [@3813973]。
此外，通过解析更精细尺度的物理过程，我们成了它们更快时间尺度的奴隶。显式数值格式的稳定性由Courant-Friedrichs-Lewy (CFL) 条件决定，该条件直观地指出信息在每个时间步长内传播的距离不能超过一个网格单元。非静力模型不仅必须遵守平流（流速本身）施加的 CFL 限制，还必须遵守它们所解析的最快波动的限制。这些波动包括内重力波，其速度可能相当高，迫使模型采用非常小的时间步长 [@3813973]。例如，在一个发展的对流羽流中，浮力可能导致垂直速度迅速加快，这反过来又收紧了垂直平流 CFL 限制，通常需要一个随着上升气流增强而缩小的自适应时间步长 [@3922276]。
如果模型是完全可压缩的，它还必须解析声波。由于空气或水中的声速极高（每秒数百至一千多米），一个完全显式的可压缩模型将被迫采用不切实际的微小时间步长 [@4068863]。这就是为什么许多非静力模型使用“声波过滤”近似的原因。玻色尼斯克近似假设密度处处恒定，除了在浮力项中，这对于密度变化很小的海洋来说是完美的。它将连续性方程简化为 。非弹性近似在大气科学中很常见，限制稍少，允许背景密度随高度变化，从而得到一个连续性方程 [@3802879]。两者都优雅地过滤掉了声波，同时保留了必要的非静力动力学。为了处理剩下的快波，模型开发者们设计了巧妙的数值方案——如分裂显式或半隐式方法——它们以不同于慢速平流项的方式处理快波项，从而允许一个更大、更高效的整体时间步长 [@4068863]。
我们为什么要费这么大的劲？因为回报是一个能够捕捉其静力对应模型所无法看到的现象的真实物理过程的模型。这对理解海洋和大气中能量如何被移动和转换具有深远的影响。
内波：这些波沿着流体内部的密度面传播，是传输能量的关键机制。非静力模型正确地预测了它们的色散关系，该关系决定了它们的速度如何依赖于波长。相比之下，静力模型对这种关系的描述从根本上是错误的，特别是对于那些垂直范围相对于水平波长较大的波。静力模型预测这些波以人为的高群速传播。这意味着静力模型会错误地模拟大部分内波谱的能量传播，这是理解海洋混合的一个关键缺陷 [@3925171]。
湍流：湍流是混合热量、盐分和动量的混乱三维搅动。在分层流体中，一个关键过程是通过浮力产生湍流动能（TKE）：温暖、轻的流体团上升（），而寒冷、密的流体团下沉（），从而将势能转化为动能。这个过程由浮力通量项 表示。在非静力模型中，垂直速度 是一个完全动态的、预报的变量，这种能量转换是 TKE收支的自然组成部分。然而，在靜力模型中， 仅仅是水平流动的诊断性附属物；没有垂直动能的预报方程。浮力产生垂直湍流的物理途径被切断了。因此，靜力模型在结构上无法表示地球物理湍流的这一基本方面 [@3802868]。
最终，在静力模型和非静力模型之间的选择，是关于我们希望模拟哪个宇宙的选择。静力世界是一个简化的、优雅的近似，对于描绘地球环流宏大而缓慢的舞蹈非常强大。非静力世界则是一个更丰富、更复杂、更真实的现实，一个充满破碎波浪、翱翔热气流和湍流涡旋的地方——在这个世界里，通过在其完整的垂直荣耀中拥抱 这一简单原则，我们开始能够以其最真实的形式看待自然。
在了解了区分静力与非静力世界的原理之后，我们现在来到了最激动人心的部分：看这些模型在实际中的应用。这个看似微小的差异——保留或舍弃垂直加速度项——究竟在何处真正重要？答案是，它解锁了一系列丰富的现象，从山间微风的低语到飓风眼墙的咆哮。它使我们能够为地球的流体系统构建更逼真的“数字孪生”。
要真正体会其中的差别，想象一下聆听一场管弦乐。静力模型就像戴着耳罩听音乐，所有高音——短笛、小提琴、铜钹——都被过滤掉了。你能从大提琴和低音提琴中捕捉到基本的节奏，但却错过了织体、细节和那些突然、尖锐的能量爆发。非静力模型则摘下了耳罩。突然之间，你听到了完整的交响乐。我们用来观察这一点的工具是动能谱，它告诉我们在从广阔的天气系统到单个云的不同空间尺度上存在多少能量。静力模型捕捉了大尺度流动的能量，但在几十公里的尺度上显示出陡峭且不切实际的能量衰减。相比之下，非静力模型揭示了一个在这些较小的“中尺度”上充满能量的世界，其能量谱斜率通常更接近 ，这是湍流运动的特征。这就是静力近似所无法看见的风暴、波浪和环流的能量。让我们来探索一些这些新听到的音符。
大气是浮力、惯性和压力相互作用的画布。非静力模型是我们描绘其最复杂细节的最精良画笔。
想象一下一股稳定的风流过山脉，就像水流过溪流中的石头。如果山脉非常宽阔平缓，如同海洋中的緩慢涌浪，空气有足够的时间进行垂直调整。流动平稳，且很大程度上是静力的。
但如果山脉尖锐而狭窄呢？空气被迫猛烈向上运动。它无法保持静力平衡；它会垂直加速。这种扰动产生了向下游涟漪般传播的波，更值得注意的是，这些波会垂直传播，远达平流层。这些就是“山地波”或“背风波”。要正确模拟它们，模型必须是非静力的。关键在于比较波的固有频率（由风速 和山脉宽度（与波数 相关）设定）与大气的自然垂直振荡频率，即 Brunt–Väisälä 频率 。当固有频率 并不远小于 时，非静力效应占主导地位。
这些波不仅仅是飞行员观察到的美丽的透镜状云。它们携带大量的动量。当它们向上传播时，可能会像海浪拍岸一样破碎。当它们破碎时，会释放其动量，对大气产生强大的“拖曳”力，这可以引导急流并影响全球天气模式。非静力模型可以通过模拟波及其动量输送来明确计算这种“解析的山地波拖曳”，这个过程由速度扰动的协方差 描述。对于那些无法解析每座山脉的气候和天气模型，这种效应必须通过“次网格地形拖曳参数化”来粗略估计。随着我们的模型变得越来越非静力化且分辨率提高，它们可以直接看到更大比例的这种关键拖曳力，使我们的全球预报更加准确。
当然，这也产生了一个新的技术问题：当这些垂直传播的波到达模型的顶部，即人造的“盖子”时，会发生什么？一个简单的刚性盖子会像镜子一样，将波能反射回来，污染模拟。为了防止这种情况，模型开发者在顶部附近放置了一个“海绵层”。这是一个在方程中加入了摩擦项（称为瑞利阻尼）的区域。这个阻尼项 吸收向上移动波的能量，模仿它们在真实、无界的的大气中自然地向外太空传播并消散的方式。
让我们从流经山脉的稳定气流转向更剧烈的现象：雷暴。雷暴的核心是一股强大的上升气流，一团因浮力而上升的暖湿空气羽流。最简单地说，我们可以将其建模为一个“上升的热泡”。非静力模型捕捉了其基本物理过程：由浮力驱动的初始垂直加速度，以及在热泡周围形成的复杂压力场，用以推开周围的空气 [@problemülü:4068819]。
一个成熟的雷暴不仅仅是一个上升气流。它还产生强大的、由寒冷的雨水冷却的空气组成的下沉气流。当这股稠密的空气撞击地面时，它会像泼洒的油漆一样水平散开。这个前进的冷空气锋面被称为“阵风锋”，是一个微型天气锋，你可以在暴风雨来临前感觉到一阵突然的、凉爽的阵风。这个现象是密度流的完美例子。它的速度不是由大尺度风决定的，而是由其自身的密度（或“折减重力”，）与其深度 之间的精细平衡决定的。非静力模型显示，传播速度 与 成比例，这种关系由一个量级为 1 的弗劳德数 表征。阵风锋前端翻滚的、湍流的头部会将前方的暖空气抬升，并可能引发新的风暴，这是一个典型的非静力特征。
在飓风的核心地带，非静力动力学的重要性无出其右。飓风眼墙中可怕的风是一个巨大旋转系统的一部分，但风暴的引擎是眼墙内强烈的对流。此处的上升气流可以以 或更快的速度尖嘯向上。这种运动是静力的吗？
我们可以像对山地波那样，对垂直运动提出同样的问题。我们可以定义一个垂直弗劳德数，，它比较了特征垂直速度 与重力波在垂直尺度 上的速度 。对于眼墙中一个典型的强上升气流，这个数值约为 0.25。这远未达到静力近似成立所需的 的条件。垂直加速度是显著的，是风暴动力学中至关重要的一部分。
这意味着，要准确模拟热带气旋的结构，以及至关重要的强度，模型必须是非静力的。静力模型根本无法产生构成眼墙的强烈而狭窄的上升气流。它会将它们抹平，削弱风暴的次级环流——这个巨大的“飞轮”在底部吸入湿气，在顶部排出空气——并导致对风暴巅峰风速的预报不佳。能够在公里级分辨率下运行非静力模型，是现代飓风预报的重大突破之一。
海洋与大气一样，也有其复杂的垂直运动，但它们常常隐藏在我们的视野之外。在这里，静力模型与非静力模型之间的选择同样是一个尺度问题。
考虑内波，这些波不是在海洋表面传播，而是在海洋深处不同密度层（温盐跃层）的界面上传播。当潮汐来回冲刷大陆架上的水体时，可以产生波长达数百公里的巨大内潮。这些是巨大、缓慢、温和的运动。尺度分析证实，它们的频率 远小于浮力频率 。对于这些庞然大物，静力近似是完全足够的。同样，对于风暴潮——飓风期间海平面的大范围上升——水平尺度巨大（），而垂直深度相对较小（）。同样，流动是静力的，静力模型是完成这项工作的正确工具。
但是，如果我们放大观察海底陡峭尖锐地貌（如水下山脉或大陆架边缘）上的流动，情况就变了。就像在大气中一样，流动可以产生尖锐的、非静力的背风波。这些非线性内波可能极其陡峭且充满能量。要捕捉它们的动力学，非静力海洋模型是必不可少的。
这种非静力能力带来了显著的计算成本。在静力模型中，压力通过简单地积分上方水的重量来找到。在非静力模型中，压力必须通过求解一个确保流动保持不可压缩性的三维椭圆（泊松）方程来计算。这带来了一个计算难题，尤其是在使用自适应网格加密（AMR）来放大感兴趣的特征时。压力方程的椭圆性质意味着信息必须在所有网格层级之间同时通信，这比静力模型简单的柱状物理过程要困难得多。
几十年来，全球气候和天气模型完全是静力的。它们的网格单元宽达数百公里，太粗糙以至于无法看到单个雷暴。这些未解析的风暴的集体效应必须使用“对流参数化”方案来近似，例如著名的 Kain-Fritsch 方案。这些方案充当了替身，试图根据大尺度条件（如对流有效位能 CAPE）来推断对流何时何地会发生，然后将产生的加热和增湿效应添加回粗糙的网格单元中。
今天，我们正处在一场革命的边缘：全球云解析模拟（GCRM）。得益于计算能力的巨大增长，我们现在可以运行网格间距仅为几公里（）的非静力全球模型。你可能会认为这意味着我们正在“解析”云，问题已经解决。但自然是微妙的。
模型的“有效分辨率”——即它能忠实模拟的最小特征——通常是其网格间距的 6 到 10 倍（）。所以，一个 网格的模型只能真正捕捉到尺度为 或更大的现象。然而，一个典型的深对流上升气流核心只有 宽。一个浅积云可能只有几百米宽。这意味着即使是我们最好的全球模型也运行在一个“对流灰色地带”。它们对于传统的参数化方案来说太精细而无法正常工作，但对于完全解析云的内部生命来说仍然太粗糙。虽然它们现在可以明确地表示风暴组织成中尺度系统的宏伟结构，但它们仍然在处理基本的羽流动力学方面存在困难。
这就是前沿。现在的挑战是开发“尺度感知”参数化方案，这些方案能够随着模型网格变得足够精细以至于能自己看到越来越多的对流过程而智能地调低自身作用。通过拥抱完整的非静力方程，我们并没有解决所有问题，但我们已经能够提出关于我们世界天气和气候多尺度性质的更深刻、更有趣的问题。我们摘下了耳罩，终于开始听到这个星球完整、复杂而美妙的交响乐。