很多人对流体力学的最初印象，来自初中物理那杯被大气压托住的倒扣水杯，或者是洗手池里靠近下水口时加速旋转的水涡。但真正在工程里打交道之后你会发现，流体力学不是一道难解的物理题，而是一整套描述“力与运动之间如何相互纠缠”的思维模型。它是航空航天、能源动力、化工过程、环境工程、土木结构、运动装备，甚至医疗器械和城市排水设计都会用到的底层框架。这篇作为知识体系的第二篇，我准备把流体力学内部最重要的东西，按一条从静止到运动、从理想到真实、从解析到计算的路线完整捋一遍。不管你是在校学生、转行工程师，还是自学科普爱好者，这篇文章都可以帮你在脑子里搭起一张可复用的知识地图。
1. “连续介质”假设：整个流体力学的地基到底长什么样
1.1 分子世界与宏观世界之间，我们选择了后者
流体力学研究的是液体和气体的运动。但真要严格较真起来，无论是水还是空气，都是由大量分子组成的离散系统——每个分子都在做无规则热运动，彼此碰撞，完全不是数学上光滑连续的。如果从分子尺度去看，密度、压强、速度这些概念几乎都没法定义。
工程上我们处理的是宏观问题：管道里每秒流过多少水，机翼上承受多大升力，城市风道里气压损失多少。这些问题关心的尺度从毫米到几十米，远远大于分子平均自由程。于是这里做了整个学科最关键的一个操作——连续介质假设：把流体看成是连续填充整个空间、没有缝隙的物质，并用密度、压力、速度等场变量来描述它的状态。
这个假设的价值在于，原本需要统计大量分子运动才能拿到的信息，被压缩成了一个个随时间和空间连续变化的物理场。你可以用偏微分方程来描述它们，可以用微积分去推导，可以用数值方法去离散求解。如果放弃这个假设，流体力学就退回到了气体分子动理论，工程上几乎寸步难行。
1.2 什么时候连续介质假设会失效
一个典型的失效场景是高空稀薄气体。当飞行器飞到上百公里高度，空气密度极低，分子的平均自由程可以长达几米甚至几百米，此时连续介质假设的基础就没有了。还有一种场景是微纳尺度流动，比如微流控芯片里几十微米宽的通道，如果气体压强低到一定程度，壁面滑移现象就会明显到不能忽略。
不过对于绝大多数工程问题，比如输水管道、建筑风压、车辆空气动力学、涡轮机械内部流动，连续介质假设都是成立的。所以它虽然只是个假设，却是整个流体力学生存的基石。你甚至可以说，流体力学本质上是研究“在连续介质假设之下，流体微团如何运动和受力的学问”。
2. 静力学与基本属性：流体的密度、粘性和可压缩性
2.1 密度与可压缩性：什么时候“体积不变”成了默认选项
流体的基本属性里，最直观的是密度。水的密度大约1000 kg/m³，常温空气的密度大约1.2 kg/m³，两者差了将近三个数量级。这就是为什么同一种风速下，水对结构的冲击载荷会比空气大得多。
密度之外，工程师更关心流体到底能不能被压缩。液体的可压缩性非常弱——常温下水每增加一个大气压，体积只缩小约万分之五。所以在大多数管道和明渠问题里，水直接被当作不可压缩流体处理。气体就不同了，气体本身很容易被压缩，但当流速较低时，密度变化引起的效应也很小。工程界有一条沿用已久的经验线：当流速对应的马赫数低于0.3时，气体也可以按不可压缩流体来算。马赫数是流速与当地声速之比，常温下声速约340 m/s，所以0.3马赫对应的风速大约102 m/s——一般建筑通风、汽车行驶、低速风洞都远达不到这个量级，这也是为什么大部分CFD（计算流体力学）软件里默认用不可压缩求解器就够用。
2.2 粘性：流体阻力的微观根源
粘性可能是流体力学里最容易被低估的物理量。所谓粘性，是流体抵抗剪切变形速度的一种属性。微观上看，粘性来自分子之间的动量交换：流速快的层里的分子，随机热运动到流速慢的层里，把动量带过去，宏观上表现为阻碍相对滑动的力。
水的粘度大约1.0×10⁻³ Pa·s，空气的大约1.8×10⁻⁵ Pa·s。有意思的是，液体的粘度随温度升高而下降，气体的粘度却随温度升高而上升。原因是液体中分子间距很小，粘性主要由分子间作用力主导，温度升高作用力减弱；气体中分子间距大，粘性取决于分子热运动的动量交换，温度越高交换越强。这个差异在实际工程里很坑：设计冷油管路和热油管路时，粘度变化会导致截然不同的压降，判断错一个数量级的事我都见过。
2.3 帕斯卡原理与静压分布：静止流体内的力场
流体静力学的核心，是研究流体处于静止状态时的应力分布。这里最著名的结论有两个。第一个是帕斯卡原理：封闭容器内的静止流体，其某一部分的压力变化会等值传递到流体的每个部分。液压千斤顶就是靠这个原理，用小活塞上的小力在大活塞上产生大力。
第二个结论是静压公式：在重力场中，静止流体某一点的压力等于p₀+ρgh，其中p₀是参考面压力，h是深度。这个公式能解释为什么大坝底部要建得更厚，为什么水越深压强越大，为什么潜水员下潜深度每增加10米，大约就要多承受一个大气压的压力增量。
静止流体中另一个高频考点是浮力。浸在流体中的物体，会受到一个向上、大小等于所排开流体重力的浮力。这个结论在船舶设计里叫排水量，在热气球升空里叫浮重差，在潜水艇里叫浮力调节。很多人把浮力当成一个初中就学过的简单常识，但工程应用时会发现它的边界并不简单：比如浸入部分会不会变化、自由液面是否倾斜、旋转坐标系下浮力方向如何，这些延伸问题很考验对压强分布的理解。
3. 流动描述与三大基本定律：从流线到方程
3.1 两种观察视角：欧拉观点与拉格朗日观点
进入运动学，流体力学首先给出的是描述流场的两种视角。
拉格朗日观点是跟着流体微团走，记录某一个微团在不同时刻的位置、速度和加速度。这就像你坐一艘小木筏顺河漂流，只关心自己这条轨迹上的变化。
欧拉观点则是站在固定空间位置上看，记录每个空间点在不同时刻的流速、压力和温度。这就像站在桥上看水流，关心的是某个固定点处的流速如何随时间变化。
工程上绝大部分都是采用欧拉观点，因为你在管道、风机、换热器里布置的传感器测的都是固定位置随时间的变化，而不是追踪某团流体。但两者之间的换算并不难，关键是要理解速度场的物态导数（随体导数）：某个物理量随流体微团的变化率，等于当地变化率加上对流变化率。你看到一台风扇前的风速稳定，叫当地稳定；但如果移到一个正在加热的管道里，某点的温度随时间变化来自两部分——热源加热的当地变化，以及上游高温流体被带过来的对流变化。两种变化的区别，是很多初学者一学伯努利就懵的根子。
3.2 连续性方程：质量守恒是最不受质疑的一条
连续性方程说的是，控制体内流体的质量变化率等于净流入的质量通量。对于不可压缩流动，它会退化成更漂亮的表达式：流入体积流量等于流出体积流量，即A₁v₁=A₂v₂。
这条公式看起来简单，用起来非常广。比如自来水管从粗管进入细管，流速必然会增大，因为体积流量要保持一致。河流从宽段流入窄段时水流加速，同样是这个道理。连续性方程也是后续推导边界层方程、数值模拟中压力修正算法的基础。它是一个纯粹的运动学约束，不涉及力，因此永远不会因为“力算错了”而失效。
3.3 伯努利方程：能量视角的锋利与边界
伯努利方程可能是流体力学里最出名的公式：p + ½ρv² + ρgz = 常数。它本质上是在流线方向上，忽略粘性损失和不可压缩条件下，对机械能守恒的一种表达。压力能、动能、势能三者沿流线守恒。
它好用，因为只需要知道两点之间的高度差和速度变化，就能反推压力变化，不需要解复杂的偏微分方程。比如机翼上方流速比下方快，按伯努利确实能解释上下表面存在压差。但这里要特别提醒：伯努利方程的适用条件是定常、不可压缩、无粘、沿流线。真实流动中如果出现分离、剧烈弯曲、旋转流动，直接套伯努利就很容易得出错误结论。
我见过不少初学者拿伯努利解释一切升力现象，说“上面流速快所以压力低”。这句话本身没错，但只讲了一半。机翼背后的增升机制还涉及流动偏转和环量，也就是动量定理的贡献。把因果简化成“流速差→升力”，会在解释等离子体翼、螺旋桨、地面效应时彻底翻车。
3.4 纳维-斯托克斯方程：流体世界里的牛顿第二定律
把牛顿第二定律应用到流体微团上，得到的正是纳维-斯托克斯方程（N-S方程）。它的物理含义非常直白：惯性力 = 压力梯度力 + 粘性力 + 体积力。对于不可压缩流体，它可以写成向量形式的偏微分方程，包含速度的三个分量和压力四个未知量，配合连续性方程正好形成封闭。
这个方程的难点在于非线性项，也就是对流加速度项。它让两个不同尺度的扰动可以互相耦合，让流动表现出丰富的失稳和转捩行为。这就是为什么N-S方程在数学上至今没有普遍的解析解，也是为什么一百多年来它一直是应用数学和工程研究的热点。
工程上对N-S方程不是非要解出通解，而是想办法针对具体问题做简化。边界层理论，是把方程限制在壁面附近的一个薄层里，忽略沿流向的粘性扩散；势流理论，是假设流体无粘无旋，然后用拉普拉斯方程去解速度场；润滑理论，则利用“几何上很薄”的特点把方程简化成雷诺方程。可以说，流体力学的工程智慧很大程度上体现在“什么时候可以忽略哪一项”。
3.5 雷诺数：一个数决定流动的性格
雷诺数是流体惯性力与粘性力的比值，公式为Re=ρvL/μ。它神奇之处在于，它把流速、特征尺度、密度和粘度四个量压缩成了一个无量纲数，而这个数的大小直接决定了流动是层流还是湍流。
管流中，Re<2300左右是层流，流动平滑，速度剖面呈抛物线；Re>4000基本进入湍流，速度剖面更扁平，掺混增强，阻力规律也完全不同。生活里有个非常直观的例子：水龙头开很小的时候水流是透明玻璃棒一样的光滑水柱，开大之后表面开始波动破碎，那就是雷诺数升高、层流失稳转为湍流的表现。
雷诺数还是模型实验的核心依据。飞机模型在风洞里测到的力，要换算到真实飞机，首先就需要保证两者的雷诺数相同。这一点在后面讲相似性的时候会再展开。
4. 量纲分析：把实验数据变成工程公式的脚手架
4.1 为什么流体力学绕不开量纲分析
N-S方程计算复杂，很多实际问题根本得不到解析解。工程师的办法是：做实验测数据，然后把数据整理成经验公式去用。可问题来了，管道压降同时受直径、长度、流速、密度、粘度五个变量影响，如果一组一组做实验，工作量是天文数字。
量纲分析的意义就在这里：它能在不做任何实验之前，先告诉你自变量应该如何组合成无量纲数，从而大幅压缩实验变量的个数。这就是白金汉π定理的核心思想：一个包含n个物理变量的物理方程，如果涉及m个基本量纲，就可以写成n-m个无量纲数之间的关系。
以圆管内不可压缩流动的压降问题为例。压降Δp可能依赖管径D、管长L、流速v、密度ρ和粘度μ。这里面基本量纲只有三个：质量、长度、时间。六个物理变量减去三个基本量纲，至少存在三个无量纲数。把长度和直径的比L/D、雷诺数Re、以及表示压力与动压比值的欧拉数Eu拉出来，实验结果就可以整理成Eu=f(Re, L/D)这个紧凑形式，一下子把原来六个变量的问题压缩成了两个变量的函数关系。
4.2 从量纲分析到相似准则：模型实验是怎么换算的
量纲分析再往前走一步，就是相似理论。如果两个流动系统满足几何相似、运动相似、动力相似，那么它们在同一无量纲数下就具有相同的无量纲流场。这种相似关系在工程实验里太重要了：直接给大船做大尺度实验不现实，但做一个缩比船模，只要保证关键无量纲准则一致，就可以把模型上测到的力换算到实船。
问题在于，一个流动往往受多个无量纲数控制。水面船舶绕流同时受重力和粘性影响，理论上要同时满足弗劳德数和雷诺数都相等。可是同一个缩比模型，要同时满足两个准则，对流体介质的要求会落到一个根本无法实现的物理量上。所以在船舶实验中通常的做法是，用重力相似（保证弗劳德数一致）来测量兴波阻力，用粘性修正公式单独估算摩擦阻力。这种“抓住主要矛盾，次要因素用修正量覆盖”的思路，是工程流体力学的典型思维模式。
4.3 无量纲数背后都是力的对抗
学量纲分析最忌讳死记公式，不理解每个无量纲数代表谁的比值。雷诺数是惯性力与粘性力之比；弗劳德数是惯性力与重力之比；欧拉数是压差力与惯性力之比；马赫数是惯性力与可压缩性力的表征。这些比值告诉你一个流动的“个性”由谁主导。
水波传播由重力主导，所以水面船兴波阻力看弗劳德数。含聚合物熔体的狭缝流动由粘性主导，所以看雷诺数。跨声速机翼表面出现激波，由可压缩性主导，所以看马赫数。你只要能在脑子里把“什么力在推动”“什么力在阻挡”这两件事理顺，大部分相似准则的选择都不会出错。
5. 边界层、阻力与升力：粘性如何主宰真实流动
5.1 无滑移条件与边界层的诞生
理想流体力学有个不太真实的前提：流体在壁面处可以自由滑动。但真实流体因为粘性的存在，在固壁表面速度必须等于壁面速度——这就是无滑移条件。由此带来的后果是，壁面附近必然存在一个速度从零逐渐过渡到主流速度的区域，这个薄薄的一层就是边界层。
边界层的最初概念，是普朗特在二十世纪初提出的。他的洞察在于：对于高雷诺数流动，粘性的影响被压缩在壁面附近一个很薄的层内，层外可以当作无粘流动处理。这一刀把N-S方程劈成了两个相对简单的问题：外部用势流理论解，内部用边界层方程解。整个现代空气动力学的地基，就建立在这个“一刀切”之上。
边界层厚度的估算，很多教材里给的是δ≈5x/√Reₓ。其中x是从前缘算起的流向距离，Reₓ是当地雷诺数。这个公式对层流边界层很好用。举个例子：翼展方向的典型弦长是1米，空气流速5 m/s，那么前缘下游1米处的雷诺数大约3.4×10⁵，边界层厚度大约8毫米。这个量级告诉你，机翼表面的粘性影响层其实非常薄，往往只有几毫米到几厘米，但就是这层薄薄的流体，决定了摩擦阻力、流动是否分离以及最终升力是否达到预期。
5.2 摩擦阻力与压差阻力：两种阻力的物理根源
阻力从来不是单一来源的。物体在流体中运动时，受到的阻力至少可以拆成两大类。
摩擦阻力来自壁面切应力，本质上是流体和壁面之间的动量交换。它与接触面积直接相关，边界层是层流时摩擦阻力小，湍流时掺混增强、速度剖面更饱满，壁面切应力反而更大。
压差阻力来自物体前后表面压力分布的不对称。当流动发生分离，物体背风面会形成一个低压回流区，前部和后部之间形成巨大的压力差。一个典型的例子是垂直于来流放置的圆盘，背风面完全被涡街和回流主导，压差阻力极有可能会占到总阻力的九成以上。
这两个阻力的规律完全不同。高层建筑、卡车、后视镜这些钝体，主要难点在压差阻力；而飞机机翼、潜艇壳体这些流线型体，主要难点在摩擦阻力。这也是为什么汽车外形从方盒子变成流线型时，风阻系数从0.6降到0.3以下，省下来的主要是压差阻力的部分。
5.3 阻力计算的工程套路
工程上计算阻力有一条通用路径：先判断绕流形态，再用阻力系数法估算。公式是F_D = ½ρv²AC_D。其中A是参考面积，C_D是阻力系数，它通常由实验给出，强烈依赖雷诺数和物体形状。
光滑球体在高雷诺数下C_D约0.4左右，低雷诺数层流状态下则满足斯托克斯阻力定律，理论解给出F_D=6πμRv，等效C_D=24/Re。这个概念经常用在微颗粒沉降、喷雾液滴、血液细胞流动这类问题里。对汽车来说，典型C_D在0.25到0.35之间；卡车因为外形方正，C_D往往在0.7到0.9。
一个非常有意思的估算：人骑自行车时的迎面阻力。假设人车正投影面积约0.5 m²，C_D取1.0，空气密度1.2 kg/m³，骑行速度10 m/s，那么阻力约30牛。这30牛的力看似不大，但在长时间骑行中会直接转化为体力的额外消耗，这也是为什么公路车手会追求低趴姿态、紧凑骑行服和圆头盔——每减少0.01的C_D和0.01平方米的迎风面积，在长距离骑行里都会被放大成秒级的差距。
5.4 升力是怎么来的：一次对“流速差”迷信的祛魅
升力的产生常被简化成“伯努利+上下表面流速差”，但真正做气动设计的人必须看得更准。升力更本质的来源是流体对物体施加的动量变化率：机翼把来流向下偏转，气流反过来给机翼一个向上的反作用力。这个向下偏转的强度与机翼表面形成的环量成正比，而与“上表面路程更长”没有直接关系。
奥地利的马格努斯效应就是一个很有说服力的例子：旋转的圆柱在均匀来流中会产生升力，其上表面和下表面根本没有“路程差”。升力的正确解释其实是环量叠加势流后的非对称流动，再加上边界层不分离作为前提。一旦大攻角下边界层严重分离，环量被破坏，升力也会骤然下降——这就是失速。
6. 湍流、仿真与工程判断：不可解方程的现实出路
6.1 湍流不是“乱流”：它有多尺度结构
湍流这个词在工程里名声不好，听起来像一团糟、无法预测。但现代湍流理论告诉我们，湍流并不是完全随机的噪声，它是由一系列尺度范围跨度极大的涡结构组成的。大尺度的涡从主流中抽取能量，然后不断分裂成更小尺度的涡，最终在小尺度上通过粘性耗散变成热能。这个能量从大涡到小涡，再到耗散的级串过程，就是湍流能量级串。
工程计算里，直接解析所有尺度涡的做法叫直接数值模拟DNS。这个概念听起来美好，但真实雷诺数下涡的尺度跨度非常大，要达到解析全部尺度的计算量，对大多数工程问题而言是天文数字。所以工程上普遍采用的方案是雷诺平均NS方程（RANS）：把速度和压力分解成时均量和脉动量，然后对脉动量带来的附加应力做模型化处理。
6.2 湍流模型的工程选择与CFD的实操经验
常用的RANS模型里，两方程的k-ε模型在工业界地位稳固，其中k代表湍动能，ε代表湍流耗散率。它的变体Realizable k-ε、RNG k-ε在分离流和旋流预测上有改善；k-ω SST模型则在翼型绕流和逆压梯度流动中表现更好。选择模型的逻辑不是越复杂越好，而是看你的问题里哪一种物理效应最重要。
CFD不是“画个网格点运行就行”。我做过多次实际计算之后有一条特别深的体会：先估算的永远是雷诺数和近壁面尺度，再用这些量去设计网格。近壁面第一层网格高度取决于你使用的壁面处理方式。如果你追求解析边界层，那么要求第一层网格的y⁺接近1，这意味着靠近壁面的网格要加密到毫米甚至微米量级；如果你使用壁面函数，第一层网格的y⁺可以放到30到300之间。很多人计算发散，排查三天后发现只是第一层网格设错了一个数量级。
还有一条很实际的经验：不要只盯着残差曲线到10⁻⁴就觉得计算收敛了。残差下降只能说明迭代过程趋于稳定，不代表物理量已经收敛。应该同时监视你关心的物理量，比如升力系数、阻力系数、出口流量有没有稳定下来。有一次我在做一个某风道系统改造的模拟，残差掉得很漂亮，但出口总压却持续缓慢振荡，最后发现是边界层分离区里出现了低频涡脱落，改用非定常计算才拿到合理结果——这类问题靠“残差好看”是骗不过去的。
6.3 工程判断力：流体力学给你的最后一份礼物
学习和使用流体力学的最终目的，不是把每个公式背出来，而是建立量级感和因果链。拿到一个项目，你需要迅速判断：这个流动是层流还是湍流？雷诺数在什么范围？是否涉及可压缩性？粘性主导还是惯性主导？重力是否重要？边界层会不会分离？这些问题本质上都是量纲分析加物理直觉的组合。
我个人的习惯是，在真正开算或者做实验之前，先手算几个无量纲数，再用简单公式估算一下目标量级。比如估算管道流量时，先算雷诺数；估算风载荷时，先拍一个C_D范围，套阻力公式算出量级。这套“先粗算、后精算”的动作，帮我在很多场合避免了拿着CFD结果却完全看不出对错的尴尬。计算机会告诉你一个数字，但它不会告诉你这个数字是否合理。流体力学给你最重要的能力，恰恰是判断那个数字靠不靠谱。
最后分享一个我保存多年的小做法：把常见场景的雷诺数量级做了一张卡片——人的游泳约10⁶，自行车骑行约10⁵，室内通风管道约10⁴到10⁵，微流控芯片约10⁻¹到10²。每次遇到新的流动问题，先对号入座，心里就有了谱。这个习惯让我在处理陌生领域时，至少能快速判断出该用粘性主导的理论还是惯性主导的理论，该选层流公式还是湍流公式。希望它对你也有用。
