1.3 压缩性与热物性
1.3 压缩性与热物性
本节摘要：压缩性度量流体密度随压力变化的敏感程度，由体积模量 K 定量刻画——水约 2.2 GPa，常压气体等温过程约 0.1 MPa 量级，相差四个数量级。压缩性通过声速与马赫数影响流动分类，决定何时可以采用不可压缩近似。本节用 Python 计算声速、水击压强与可压缩近似的适用阈值。
深海一万米的马里亚纳海沟底部，压力约 110 MPa，相当于每平方厘米压着一吨多的东西。那里的水被压密了多少？答案约 5%。这就是液体的"硬"：体积模量 K 定义为产生单位相对密度变化所需的压力增量，K = ρ dp/dρ。水的 K ≈ 2.2 GPa，要把水压缩 1%，得上 22 MPa——日常工况里液体密度几乎是个常数，这正是"不可压缩假设"的底气。
一、体积模量与声速
体积模量不只描述"难压程度"，它还直接决定压力扰动的传播速度。推导物理图像很直观：声波是微弱的压缩-膨胀波，介质越"硬"（K 大）、越轻（ρ 小），波传播越快。理想气体中声速还有一个只依赖温度的漂亮表达式 c = √(γ R T)，γ 为比热比、R 为气体常数。空气 15 °C 时声速约 340 m/s，水约 1480 m/s——声音在水里跑得比在空气里快四倍多，鲸鱼的低频叫声能横跨半个洋盆。
import numpy as np # 声速对比 K_water, rho_water = 2.2e9, 998.0 c_water = np.sqrt(K_water / rho_water) gamma, R, T = 1.4, 287.0, 288.15 # 空气 15 摄氏度 c_air = np.sqrt(gamma * R * T) print(f"水中声速 {c_water:.0f} m/s") print(f"空气声速 {c_air:.0f} m/s") print(f"高空 -50 摄氏度空气声速 {np.sqrt(gamma*R*223.15):.0f} m/s")
第三个数很重要：高空更冷、声速更低，同样表速的飞机在高空马赫数更大。协和客机巡航在 18 公里高度，那里声速只有约 295 m/s——"超音速"的门槛比海平面低了不少。
二、什么时候必须考虑压缩性
工程上的经典判据是马赫数 Ma = U/c。经验规则：
- Ma < 0.3：密度变化不超过约 5%，可安全使用不可压缩模型（伯努利、泊肃叶都按此办理）；
- Ma > 0.3：气体的密度变化开始影响动量与能量收支，必须按可压缩流处理（第 5 章的主战场）；
- 液体：除非水击、液压冲击、水下爆炸这类压力剧变工况，一律按不可压处理。
注意"Ma<0.3"只对气体流动的动量方程宽松，热对流的密度浮升力效应（自然对流）是另一回事，微小密度差配上重力照样驱动整个流场，不能因为这个判据就把 Boussinesq 近似扔掉。
现场案例：水击——液体也会"发飙"的场景
快速关闭管道末端的阀门，流速动量在瞬间清零，动能转成压力波在管里来回拍。刚性管理论给出的压力升高估算为 Δp = ρ c ΔU。取水流速 2 m/s：
rho, c, dU = 998.0, 1480.0, 2.0 dp = rho * c * dU print(f"瞬时关阀水击压升 {dp/1e6:.1f} MPa，约 {dp/101325:.0f} 个大气压")
输出约 3.0 MPa、29 个大气压——一次猛关阀门就够让铸铁管接头开裂。真正的电站与输水工程靠缓慢关闭阀门（延长关闭时间让波反射抵消）、设置调压井或水击容纳气室来对付它。这个案例说明："液体不可压"是关于密度变化的陈述，不是"液体不会传递压力波"——波照样以 1480 m/s 拍过去，只是波形里密度几乎没变。
三、热物性一览
可压缩流动与传热问题还需要几个热学参数：比热容（定压 cp、定容 cv，空气 cp ≈ 1005 J/(kg·K)）、比热比 γ = cp/cv（空气 1.4，单原子气体 5/3）、热导率、热膨胀系数。它们通过状态方程 p = ρ R T 把压力、密度、温度锁在一起：定常能量方程因此不是可有可无的附录，而是与动量方程平起平坐的主方程。第 5 章推导等熵关系式时，γ 会反复出场。
💡 关键直觉：不可压缩不是流体的属性，而是流动的属性。同一份空气，低速风扇里可当不可压流算，涡轮喷气发动机里必须全套可压缩方程。"液体=不可压、气体=可压"是方便的记忆，不是物理定律。
四、等温与绝热：气体压缩的两种极限
同样把气体压力翻倍，等温过程（充分换热、温度不变）与绝热过程（来不及换热）消耗的功不同，"等效体积模量"也不同：等温时 K 等于压力 p 本身，绝热时 K 等于 γp。常压下空气两者都在 0.1 MPa 量级——比水小四到五个数量级，这就是"气体软、液体硬"的定量表达。声波传播是快到来不及换热的过程，所以声速公式里用的是 γ 而不是 1，这也是上一节声速表达式的来历。
工程判断口诀：过程快（声波、快速阀门、气动冲击）按绝热；过程慢且换热好（储气罐缓慢充放、长管低温输送）按等温；都靠不上时按多变过程处理，指数介于 1 与 γ 之间，由实验标定。
💡 记一个数量级锚点：要把常压空气压缩到密度明显变化只需零点几兆帕的压差；而水密度变化百分之几就要吉帕量级——差两万多倍。选"可压还是不可压"模型时，这道鸿沟就是底气。
变式思考
把体积模量与声速的关系倒过来念也成立：测出声速就能反推介质的等效体积模量，这正是超声物性检测与油品在线分析的原理——声速对温度、组分、含气量都敏感，一支声速探头就是一台成分分析仪。液态食品的浓度控制、液压油混入空气的报警，都在用这条链路。物性、声速、测量，三者在工程里互相兑换。
本节要点回顾
- 体积模量 K：单位相对密度变化所需压差；水约 2.2 GPa，气体小四个量级；
- 声速 c = √(K/ρ)（液体）或 c = √(γRT)（理想气体），声速随气温降低而减小；
- Ma < 0.3 判据：气体流动密度变化小于约 5%，可用不可压缩近似；超过则进入第 5 章的可压缩世界；
- 水击：Δp = ρ c ΔU 可达数十个大气压，是液体里压缩性最凶猛的出场方式；
- 热物性（cp、γ、状态方程）把可压缩流动的能量收支与动量收支耦合起来。
至此四张身份卡集齐。下一章让流体先静止下来，看它如何把力递到四面八方。