第二章 流体力学基础:流体的定义与分类
各位同学,欢迎来到流体力学基础这一章。说实话,很多刚接触CFD的朋友,一上来就急着学软件操作,结果算出来的结果连自己都不敢信。为什么?因为基础没打牢。今天我就带大家把流体的“脾气秉性”摸清楚。
核心知识点速览
- 流体的定义:能流动的物质,没有固定形状
- 流体分类:牛顿流体 vs 非牛顿流体
- 关键属性:密度、粘度、压力
- 流动状态:层流与湍流
- 边界层:贴近壁面的“特殊区域”
2.1 流体的定义——它和固体到底差在哪?
流体,说白了就是“会流动的东西”。你想想看,你把一块铁放在桌上,它不动;你把一杯水倒出来,它立刻散开。这就是区别——流体在微小剪切力作用下就会持续变形。
我个人习惯把流体想象成“一群不团结的分子”。固体分子之间抱得很紧,你拉我扯;流体分子呢?各玩各的,稍微给点力就散伙。嗯,这个比喻虽然糙了点,但道理没错。
我在项目中遇到过一位机械工程师,他坚持认为空气“不算流体”,因为看不见摸不着。结果做风冷散热仿真时,边界条件设错了,温度高了20度。所以记住:气体和液体都是流体,只是压缩性不同。
2.2 牛顿流体与非牛顿流体——你的牙膏是流体吗?
这个问题很有意思。我们先把流体分成两大类:
| 类型 | 特点 | 常见例子 |
|---|---|---|
| 牛顿流体 | 粘度不随剪切速率变化 | 水、空气、甘油、机油 |
| 非牛顿流体 | 粘度随剪切速率变化 | 牙膏、血液、番茄酱、泥浆 |
牛顿流体最听话。你给它多大的力,它就流多快,线性关系。我刚开始做CFD时,默认所有流体都是牛顿流体,结果算牙膏挤出过程,完全不对。后来才意识到——非牛顿流体的粘度是“看人下菜碟”的。
非牛顿流体又分几种:
- 剪切稀化:越搅越稀(番茄酱、油漆)——你拍瓶子底部的动作就是在剪切它
- 剪切稠化:越搅越稠(玉米淀粉溶液)——你跑上去它就硬,你慢慢走它就软
- 宾汉流体:有“屈服应力”,力不够就不动(牙膏、牙膏)——挤牙膏时第一下最费劲
避坑指南:我曾经在仿真血液流动时,默认用了牛顿模型。结果导师说:“你血管里的血细胞是活的,不是水分子。”后来改用非牛顿模型,结果才合理。所以做仿真前,先搞清楚你的流体“是什么性格”。
2.3 流体三大属性——密度、粘度、压力
这三个参数,是CFD仿真的“三驾马车”。少一个,你的模型就瘸腿。
2.3.1 密度(ρ)
密度就是单位体积的质量。水的密度是1000 kg/m³,空气是1.225 kg/m³(标准状态)。
我建议你记住一个关键点:密度变化意味着可压缩性。气体密度随压力变化大,液体几乎不变。做高速空气动力学时,必须考虑可压缩性;做水管流动时,忽略它问题不大。
2.3.2 粘度(μ)
粘度是流体的“粘稠程度”,说白了就是流体内部的摩擦力。蜂蜜粘度大,水粘度小。
这里有个容易混淆的概念:动力粘度μ和运动粘度ν。ν = μ/ρ。我刚开始学的时候也搞混过,后来记住一句话:运动粘度是“去掉密度影响后的粘度”,用来比较不同流体的流动性更公平。
常见流体粘度参考(25°C)
- 水:0.001 Pa·s
- 空气:1.8×10⁻⁵ Pa·s
- 机油:0.1-0.5 Pa·s
- 蜂蜜:2-10 Pa·s
2.3.3 压力(p)
压力是流体对壁面的“推力”。在CFD中,我们通常处理的是静压和总压。
静压是流体本身的压力,总压是静压加上动压(1/2 ρv²)。伯努利方程说的就是这三者的关系。我在做管道仿真时,经常用总压来判断系统有没有能量损失——如果总压下降了,说明有摩擦或涡流。
2.4 流动状态——层流与湍流
这个问题,我敢说每个CFD工程师都纠结过。为什么?因为层流和湍流的处理方法完全不同。
判断标准是雷诺数Re:
Re = ρvL / μ
其中v是流速,L是特征长度(比如管径),μ是粘度。
- Re < 2300:层流。流体像叠在一起的薄片,各层之间不混合。我做过微流控芯片仿真,里面全是层流,两种液体可以并排流动而不混合——很神奇吧?
- Re > 4000:湍流。流体混乱、随机、有涡旋。你想想看,河流里的漩涡、烟囱冒出的烟,都是湍流。
- 2300 < Re < 4000:过渡区。这个区域最麻烦,层流和湍流交替出现,仿真时很难收敛。
注意:我曾经接手一个项目,同事把管道流动全当成层流算,结果Re=5000。算出来的压降比实测低了30%。后来改成湍流模型,结果才对。所以一定要先算雷诺数,别凭感觉猜。
湍流仿真有几种常用方法:
- RANS(雷诺平均):工程中最常用,计算量适中
- LES(大涡模拟):精度更高,但计算量大
- DNS(直接数值模拟):最精确,但只有超级计算机才跑得动
2.5 边界层概念——贴近壁面的“秘密世界”
这个知识点,我当年学的时候觉得“不就是靠近壁面的一层吗?”后来做仿真才发现,边界层决定了整个流动的阻力。
边界层是流体在壁面附近速度从0(无滑移条件)变化到主流速度的区域。它分两种:
- 层流边界层:速度分布平滑,摩擦阻力小,但容易分离
- 湍流边界层:速度分布更饱满,摩擦阻力大,但不容易分离
我在做机翼仿真时,发现一个有趣的现象:层流边界层虽然摩擦阻力小,但容易在逆压梯度下分离,导致升力骤降。而湍流边界层虽然摩擦大,但能“粘”在壁面上更久。所以很多飞机机翼故意加扰流器,让边界层提前转捩成湍流——这就是工程中的权衡。
仿真建议:做CFD时,边界层网格必须足够密。我一般会在壁面附近布置10-15层棱柱层网格,第一层高度根据y+值计算。y+是边界层网格的无量纲距离,层流要求y+≈1,湍流要求y+≈30-300(取决于壁面函数)。这个后面会详细讲。
2.6 本章知识体系总览
下面这张图,是我自己总结的流体力学基础框架。你把它记在脑子里,后面学CFD会轻松很多。
好了,这一章的内容就到这里。流体力学基础是CFD的“内功心法”,你花时间把它吃透了,后面学软件操作、调参数、分析结果,都会顺手很多。下一章我们开始讲控制方程——纳维-斯托克斯方程,那才是CFD真正的“发动机”。
公众号:蓝海资料掘金营,微信deep3321