4、流体力学基础：连续性方程、伯努利方程、纳维-斯托克斯方程简介、层流与湍流

各位工程师朋友，大家好。今天我们来聊聊流体力学。我知道，很多做散热设计的同行一听到“流体力学”四个字就头疼，觉得那是搞CFD的人才会碰的东西。其实不然。你想想看，我们做的风冷散热器，风扇吹出来的风怎么走？热管里的工质怎么循环？这些问题的底层逻辑，都离不开流体力学。

我个人习惯是，不把流体力学当成一门纯理论课来学。我更愿意把它看作一套“工具语言”。你不需要会推导所有公式，但你必须知道每个方程在说什么，以及它在我们仿真软件里扮演什么角色。今天，我就带大家把这几个核心概念捋一遍。

4.1 连续性方程：质量守恒的“记账本”

先说最简单的。连续性方程，说白了就是质量守恒。你从管道一端流进去多少质量的流体，另一端就得流出多少，中间不能凭空消失，也不能凭空产生。

数学上长这样：

∂ρ/∂t + ∇·(ρV) = 0

别被这个符号吓到。对于不可压缩流体（比如低速流动的空气、水），密度ρ是常数，方程就简化为：

∇·V = 0

这代表什么？代表流入一个控制体的流量，必须等于流出的流量。我在项目中遇到过这样一个坑：有一次做服务器机箱的散热仿真，模型跑出来温度总是偏高。我检查了半天，最后发现是某个进风口网格画得太粗，导致连续性方程没收敛，流量算少了。嗯，从那以后，我每次跑仿真都会先看一眼残差曲线里的连续性方程是否收敛到1e-3以下。

4.2 伯努利方程：能量守恒的“简化版”

伯努利方程，大家可能更熟悉一些。它描述的是理想流体（无粘性、不可压缩）沿流线运动时的能量关系：

P + 0.5ρV² + ρgz = 常数

这个公式告诉我们三件事：压力能、动能、势能，三者可以互相转换。你想想看，为什么风扇出口风速快的地方，静压会低？就是这个道理。

不过我要提醒一句：伯努利方程有严格的适用条件。它只适用于无粘、不可压缩、定常流动，而且只能沿同一条流线使用。我在做风道设计时，经常用它来估算进出口的压差，但从来不会用它来预测分离区或涡流区域的流动。

4.3 纳维-斯托克斯方程：流体力学的“终极BOSS”

好了，重头戏来了。纳维-斯托克斯方程（简称N-S方程），是所有CFD软件的底层核心。它长这样：

ρ(∂V/∂t + V·∇V) = -∇P + μ∇²V + ρg

别慌，我们拆开看：

• 左边： 惯性项，描述流体加速度。
• 右边第一项： 压力梯度，驱动流体运动的力。
• 右边第二项： 粘性项，描述流体内部的摩擦力。
• 右边第三项： 体积力，比如重力。

说白了，N-S方程就是牛顿第二定律在流体上的应用。质量×加速度 = 合力。就这么简单。

但为什么它这么难解？因为左边那个V·∇V是非线性项。这导致N-S方程几乎没有解析解，只能靠数值方法去逼近。我们做热仿真时，软件后台就是在用有限体积法或有限元法，把这个方程离散成几百万个代数方程，然后迭代求解。

4.4 层流与湍流：两种截然不同的流动状态

最后，我们聊聊流动状态。为什么同样一个风道，风速低的时候流动很“乖”，风速高了就变得乱七八糟？这就要说到雷诺数了。

Re = ρVD/μ

雷诺数Re，是惯性力与粘性力的比值。Re小的时候，粘性力占主导，流动是层流——流体分层流动，互不掺混。Re大的时候，惯性力占主导，流动变成湍流——流体质点做无规则运动，产生强烈的掺混和能量耗散。

工程上有个经验判据：

我在做散热器设计时，会刻意让空气侧流动处于湍流状态。为什么？因为湍流的换热系数通常是层流的3-5倍。你想想看，同样一个散热器，如果风道设计不合理导致局部出现层流，那一片的散热效率会急剧下降。

但湍流也有代价——压降大。所以这是一个权衡。我个人的经验是：对于强迫风冷散热，雷诺数控制在5000-20000之间比较理想。太低换热不够，太高风扇功耗受不了。

好了，以上就是流体力学在热设计中最核心的几个概念。你不需要成为流体力学专家，但掌握这些基础，能让你在看仿真结果时更有底气，知道哪些是物理现象，哪些是数值误差。

下面这张图，是我自己总结的本章知识框架，希望能帮你把这几块内容串起来。