3.3 从欧拉方程到纳维-斯托克斯方程

3.3 从欧拉方程到纳维-斯托克斯方程

本节摘要：把控制体收缩到一个点，动量账变成微分方程：无粘时是欧拉方程，加上粘性项就是纳维-斯托克斯方程——空气动力学的终极宪法。本节逐项翻译方程的物理含义，明确"哪一项在什么条件下不可省略"，并给出它的少数几个精确解的共同特征。

3.2 节的账本框还是"有限大小"的。现在让它缩成一个点：体积趋于零，进出账目按泰勒展开保留一阶项，得到的微分方程描述每一点的流动。这是空气动力学从"算总账"走向"看细节"的关键一步。

欧拉方程：无粘世界的运动定律

对无粘流动，微团只受两类力：压强梯度力与体积力（重力，气动问题中通常可忽略）。逐点应用牛顿第二定律，得到欧拉方程。它的中文版只有一句：

微团的加速度 = 单位质量的压强差推力 + 单位质量的体积力

压强梯度力值得多看一眼：流体从高压区被推向低压区，推动力正比于压强的空间变化率（梯度）。伯努利方程里"流速快则压强低"的因果方向恰恰相反——是压差驱动了加速，速度是果不是因。这个因果关系的纠正，是初学者跨过第一道概念坎的标志。

纳维-斯托克斯方程：补上粘性的一项

真实流体中，微团还受邻居的摩擦——单位质量的粘性力，正比于速度的二阶空间导数（粘度乘以速度 curvature）。加上这一项，就是纳维-斯托克斯方程（下称 N-S 方程）。它的中文版：

惯性（随体加速度） = 压差力 + 粘性力 + 体积力

三项的相对强弱由什么决定？把方程无量纲化，惯性项与粘性项之比恰好是第 1 章的雷诺数。于是 N-S 方程的"难度旋钮"只有一个：

翼型巡航雷诺数百万级，落在最后一行——这正是普朗特边界层思想（第 5 章）的合法性来源：主流用无粘方程，粘性只在一层薄壳里修账。

图：方程家族的简化树——每一章都站在某个枝头

少数派俱乐部：能解析求解的情形

N-S 方程一般情形无解析解（三维湍流甚至让解的存在性与光滑性成为千禧年难题）。能徒手解出的成员屈指可数：静止流体、泊肃叶管流、库埃特平板流、驻点流……它们共同的特征是几何简单到速度梯度结构固定。第 4 章将迎回其中两位主力。工程上其余问题交给两条路：势流近似（第 4 章）与数值求解（第 8 章）。

逐项点名：把方程读成一句话

微分方程的恐惧感多半来自"看不懂符号在说话"。把 N-S 方程按项拆开，逐项贴上身份标签：随体加速度项（惯性）管"微团的速度变化有多急"；压强梯度项管"邻居的压差在推我"；粘性项管"速度不均匀时邻居蹭我的动量"；体积力项管"重力或电磁力这类远程作用"。全方程读成一句话就是：微团之所以这么加速，是被压差推、被粘性拖、被体积力拉三件事合起来的结果。

这项点名的价值在第 5、7 章兑现：边界层理论的本质是判定"粘性项只在贴壁薄层里跟惯性项同量级"；激波分析的本质是"在间断两侧保留动量与质量的通量账"。每一项的身份清楚了，"哪一项可以扔"就不再凭感觉。

雷诺数主导的三种人生

把 N-S 方程无量纲化之后，雷诺数是唯一站在惯性项与粘性项之间的旋钮。旋钮拧到三个位置，方程呈现三种完全不同的"性格"：

• 极低雷诺数（粘性全权）：惯性项可忽略，方程退化成线性的蠕动流方程，微生物游动、微流控芯片、蜂蜜流动都在这一档——线性意味着可以叠加、可以精确求解。
• 中雷诺数（两强相持）：卡门涡街、分离泡、转捩都在这一档，无解析解也无可靠近似，是理论最薄弱、数值最烧脑的地带。
• 极高雷诺数（惯性主导）：主流几乎无粘，粘性被压缩进贴壁薄层——普朗特分层救场，工程方法全部建立在这一档。

问题：纳维-斯托克斯方程"解的存在性与光滑性"千禧年难题，和工程计算有什么关系？

关系是间接但真实的。难题问的是：三维不可压 N-S 方程在任意初始条件下是否永远存在光滑解。工程实践默认"解存在且数值方法能逼近它"，千余年来所有工程计算都建立在这个默认上。数学家尚未证完，工程师却已飞上了天——这不是矛盾，而是学科分工：方程的完备性交给数学，方程的可用域交给工程。理解这条边界，你对 CFD 结果的态度会既敢用又不迷信。

问题：湍流的随机性是方程本身带来的吗？

主流观点认为是：N-S 方程完全确定，湍流是确定系统对扰动极端敏感的产物（混沌），不是方程里"藏着随机项"。数值上初条件差一个万分之一，足够长时间后解面目全非——所以湍流计算追求的是统计量吻合，而不是逐时刻吻合。这个认知决定了 8.2 节湍流模型的设计哲学。

问题：能不能干脆永远带着粘性项算，省掉“哪项可扔”的判断？

计算上可以，判断上不可以。永远带粘性就是全 N-S 加细网格，成本立刻跳档；更关键的是“该不该简化”的判断本身就是设计洞察——知道某流动粘性只在薄层起作用，才敢用势流快速迭代外形；知道激波损失主导，才把进气道做成多级斜面。简化决策是工程判断力的核心科目，方程只是考题库。

问题：体积力项什么时候真的能扔？

重力对气体几乎总可扔：一米量级的流动区域里，重力引起的压强变化约十帕，与气动压差（动压的量级，动辄上千帕）相比可忽略。但在两个场合要请回来：大气科学的全尺度模拟（公里到百公里，重力波与浮力是主角）与自然对流的室内气流（温差驱动的浮力正是唯一动力）。扔与不扔的判据还是量级比较——把每一项都换算成压强或加速度的单位，再排座次。

本节要点

• 欧拉方程：微团加速度 = 压差力 + 体积力；压差是因，加速是果。
• N-S 方程只比欧拉多一个粘性项，但多出来的正是全部真实感。
• 雷诺数是方程的难度旋钮：极小可解析、极大可无粘近似、中间最难缠。
• 精确解俱乐部成员极少，共同的门票是几何简单导致的固定梯度结构。
• 方程可以放心地"忘掉细节"，但每一项的物理身份（惯性、压差、粘性）必须随身携带。

方程写成了两种形式，下一节谈一个纯技术却影响深远的选择：账本按"守恒式"记还是按"展开式"记。