流体力学现象

通过可视化理解 Reynolds 数与流动转捩

管内流动的转捩大约从Re 2300开始。圆柱尾流的St停在0.2附近，Re提高10倍时，DNS网格大约增到178倍左右。

2026.03.25

泊肃叶流动：解析解推导与有限差分验证

h=G=μ=1时，通道中心速度为1/2、流量为2/3。有限差分的L2从2.78e-16起停在舍入误差，格子再密也不下降。

2026.03.28

昆虫在蜂蜜中游泳 —— 雷诺数改变的流动世界

同样的水里，蓝鲸雷诺数是2.5×10^8，细菌是3×10^-5。粘性项前面的系数1/Re把两个世界隔开大约13个数量级。

2026.04.21

水柱为何会断裂成水滴？—— Plateau–Rayleigh 不稳定性与碎裂机制

半径 9.02 倍是增长最快的波长，对应 kR 约 0.697。Weber 与 Ohnesorge 再把粗液滴和雾化分成两边。

2026.04.23

漩涡的两副面孔 — 自由涡、强迫涡与亥姆霍兹

当 U=1、R=0.5 时，强迫涡核心涡度是 4（即 2ω）。自由涡外侧除中心一点外为 0，兰金接缝落在半径 r=R 上。

2026.04.28

20℃的水也会沸腾 — 空化的物理

20℃的饱和蒸气压只有2.34 kPa，60℃会升到19.9 kPa。喉部静压跌破该值时，空化数σ转为负，水便开始沸腾。

2026.04.30

水流崩散的瞬间 — 雷诺临界数与湍流边界层的多层结构

圆管的工程临界大约在Re 2300。安静的管子能把层流留到五万，贴壁流动再用y+分成粘性底层、缓冲层和对数层，第一点也跟着变。

2026.05.05

旋转的球为什么会拐弯 — 环量、马格努斯效应与库塔–儒科夫斯基定理

在U=10、R=0.5的圆柱表面上，三组环量让上表面Cp从对称的-3走到-15。到临界环量时，下表面驻点处的Cp等于1。

2026.05.07

超音速下要扩张才能加速 — 马赫数、激波与拉瓦尔喷管的悖论

面积比2.5的设计出口马赫数约为2.44，静压降到滞止压力的6.3%。背压一旦高于设计值，扩张段内部就立起激波，出口马赫随之掉下去。

2026.05.12

草履虫游过的水就是人类的蜂蜜 — 粘性的分子起源,以及 Newton、Fourier、Fick 写出同一行的原因

1个大气压的空气里，平均自由程约68 nm，平均速率约470 m/s。μ≃ρc̄ℓ/3这一行把热传导和扩散记在同一本账上。

2026.05.14

多孔介质的三张面孔 — Darcy、Forchheimer 与 Brinkman

按砾石层常数,粘性损失与Forchheimer惯性损失相抵的Darcy速度是0.1 m/s,过了这一点直线会把压降读小。

2026.05.19

两个函数相加,一个圆柱便出现了 — 复势与达朗贝尔佯谬

均匀流加上一个偶极子之后，不必先画圆，半径为R的流线自己就会闭合。环量为零时，沿圆柱表面把压力积分，阻力只剩下舍入误差。

2026.05.21

烟囱里的烟为何卷着上升而不是平躺 — 横流中的射流与反向旋转涡对

r=2、4、8时，x=10d处的平均高度分别为5.44d、8.66d、13.78d。尾流St不到同尺寸圆柱0.21的一半。

2026.05.26

流动的三幅肖像 — 流线、迹线、脉线为何分道扬镳

把收缩比降到0.25时，喉部速度变成入口的4倍，而局部项却仍是0。当V0=0.6且ω=2时，流线的y始终都停在±0.6以内。

2026.05.28

湍流同时使用三把尺子 — Kolmogorov级联与-5/3定律

Reynolds 数增加到十倍,η/L 约缩到 0.18 倍,Re=10^4 的三维 DNS 已经要大约 10^9 个单元。

2026.06.02

静止中被加热的流体 — 自然对流与瑞利数1708

一块高 0.30 m 的 60 度壁面，在 20 度空气里给出 Ra 7.03e7、Nu 54.9 和热流 197.6。

2026.06.04

粗糙的球为何飞得更远 — 阻力危机与边界层分离

光滑球与凹坑球在同一Re=1.99e5、70 m/s下分成2.049 N和0.488 N。差别是分离角82°还是120°。

2026.06.09

加热得更烫反而更凉了 — 沸腾曲线与临界热流密度

Zuber 的 CHF 是 1.259 MW/m²。Rohsenow 在 20°C 为 1.094，到 30°C 升到 3.692，已经越过天花板。

2026.06.11

想超声速就把管子加宽 — 缩放喷管与壅塞

出口面积比取 4 时设计马赫数为 2.940，出口压力约为滞止压力的 0.0298 倍。γ=1.4 的流量系数约为 0.5787。

2026.06.18

只测阻力就能算出传热 — 雷诺比拟

在Re_x=10^5的平板上，Cf/2保持为1.050e-3，当Pr为0.7、1和7时St·Pr^(2/3)都与之重合。

2026.06.23

汽油发动机为何止步于压缩比12 — 奥托与狄塞尔循环

相同压缩比下奥托在r=8为56.5%、r=12为63.0%，狄塞尔(rc=2)为49.0%与56.7%。爆震把汽油机挡在12附近。

2026.06.25

理想气体开始说谎的地方 — 压缩因子 Z 与 van der Waals 状态方程

100 bar、200 K下的氮气，压缩因子Z为0.7905。理想气体密度168.4比实际的213.0大约轻了26%左右。

2026.06.30

厨房水槽里每天都有的激波 — 弗劳德数与水跃

上游Fr约为2.86、水深0.2 m的急流，水跃之后变成Fr 0.42、水深0.714 m，比能因此损失0.238 m。

2026.07.02

始于锅底的混沌 — 瑞利-贝纳对流与洛伦兹方程

刚性壁上，对流卷要到Ra=1708才点燃。洛伦兹截断里，加热参数ρ一过24.74，这卷就不再停在定点，而是直接塌成混沌。

2026.07.09

挂在天上的波浪 — 开尔文-亥姆霍兹不稳定性

在空气与水面相接的交界面上，最小风速大约是6.6 m/s，最先增长的波长约为1.69 cm，其增长率达到了每秒4.89。

2026.07.14

旋转的球为什么会拐弯 — 马格努斯效应与环量

无规则旋转让球形弹丸侧偏。环量Γ=2πa²ω决定侧向力,spin 2的侧偏6.90 m并不是spin 1的3.71 m的两倍。

2026.07.16

单个波纹比波群跑得快 — 波的色散与群速度

深水重力波里群速度正好是相速度的一半。波长约1.73cm处相速度最慢，只有0.232 m/s，重力与表面张力在那里打平。

2026.07.21

水要流多远才变成抛物线 — 入口段长度与Hagen–Poiseuille

层流圆管里，入口处几乎平坦的速度分布要走过约0.05倍雷诺数的管径，才会收成抛物线。在那一站之前套用64/Re，沿程压降会算错。

2026.07.23

湍流加速火焰的效果在哪里折弯 — Damköhler 褶皱假说与 bending

表格里把 u'/S_L 调到 12 时，大尺度的 Σ 仍停在 4.03，小尺度则在 Ka 超过 29 后掉到 2.00。

2026.07.28

四道激波胜过一道 — 超声速进气道的斜激波串

在 Mach 3 的超声速来流里，单独一道正激波留下的恢复系数是 0.3283，四块 8° 斜板之后升到 0.8235。

2026.08.06

杀死涡的不是粘性 — 牛顿的1/r与Lamb–Oseen涡核

自由涡在半径0.5处过了50倍扩散时间，速度仍停在小数第四位。剪切应力很大也一样。n=1时粘性力为0，带走环量的是壁面。

2026.08.11

让烟先下降1 m，抽力就少了7 Pa — 富兰克林壁炉与浮力预算

烟囱的顶端仍钉在7 m，下降1 m却让净抽力从44.72 Pa掉到37.72 Pa。浮力账本把整整7 Pa的缺口算清楚。

2026.08.20

科里奥利力做的功是零，叶轮却加了1800 J/kg — 旋转坐标系里的能量与转焓

在旋转参考系里，沿流线守恒的是转焓而不是总焓。后弯0°时科里奥利做功为零，总焓仍升1800 J/kg，扬程183.5 m。

2026.08.27

总压下降1,080 Pa的地方，耗散却是零 — 当伯努利常数横跨流线时

旋转的涡核内总压下降1,080 Pa，而粘性耗散恰好为零。本文用兰金涡划出边界：总压亏损在什么时候是真实损失，在什么时候只是涡度的几何。

2026.09.03

动态系数的 49.9% 为负 — Germano 恒等式与平均操作的位置

用 Germano 恒等式测量 Smagorinsky 系数，网格点约一半给出负的值。平均取在哪里，决定了计算的稳定性。

2026.09.15

层流是2.0000，湍流是1.0322 — 速度水头前面的那个系数

用断面平均速度拼出的速度水头，不是管内流动实际携带的动能。层流α精确等于2.0000，Re=10^6时湍流α为1.0322。

2026.09.17

速度提高4倍，楔角仍是19.47° — 水深把开尔文尾迹撑开到90°的地方

深水船后尾迹的半角不随船速变化，一直停在19.47°。水深弗劳德数靠近1时同一楔角撑向90°，超过1后又按马赫锥收窄。

2026.09.24

网格加密8倍，2.4674纹丝不动 — 1807年傅里叶级数留下的两个常数

1807年审稿时被退回的傅里叶级数收敛问题，在今天的网格上以8.9490%和π²/4这两个常数回来了。

2026.09.29

用同一个式子导出μ和k，Pr偏了1.98倍 — 决定123.5 K壁温的系数9/4

同一平均自由程导出μ和k时，空气Pr被钉在γ=1.400，实测是0.707。Eucken的9/4把Couette壁温差收到123.5 K。

2026.10.01

同一个 χ=50 返回了 1759 K 和 818 K — 火焰面库折叠的位置

把标量耗散率 χ 当作表格坐标轴，熄火附近的解会分成三支。本文用 61 点火焰面求解，看进度变量 C 如何把这个折叠展平。

2026.10.06

1 μm 铝颗粒的升温比连续介质公式慢 4.6 倍 — Knudsen 过渡区的 Sherman 调和插值

1500 K、1 atm下1 μm铝颗粒的Kn为0.67。连续介质Nu=2把加热时间低估为1/4.58。Sherman调和平均把它补回。

2026.10.15

用 Kp 去除逆反应速率常数，2500 K 下 AlO 竟解离了 99.6% — 平衡常数里藏着的 (p°/RT)^Δn

把 Kp 代入 kr = kf/K，Δn≠0 的反应平衡会偏差 (p°/RT)^Δn 倍。用 AlO 解离反应器复现 SI 与 CGS 两种错误

2026.10.22

混合分数稍有波动，平均温度就降 964 K —— presumed β-PDF 与 Z=0 处的无穷大

在化学计量 Z̄ 处，β-PDF 平均温度为何比 T(Z̄) 低 964 K；当 a<1 时，均匀网格积分为何偏高 115 K，用代码验证

2026.10.27