第 3 章 · 湍流:朗道对未解难题的陈述
湍流是经典物理最后的大难题。朗道这一章没有给出解——他给出的是问题的正确提法:雷诺数足够大时流动自发失稳、时间上变得不可预测,能量从大涡逐级传向小涡。这一章的每个概念至今仍是湍流研究的通用语言。
失稳与湍流的发生
层流在 $Re$ 超过临界值后失稳:扰动被流动放大(而非被粘性抹平),级级失稳后流场在时间上呈现不可预测的脉动——湍流的操作性定义:速度场的混乱时间依赖 + 空间上的多尺度涡旋。朗道-霍普夫的原始图像(逐次分岔叠加新频率)后来被费根鲍姆的倍周期与 Ruelle-Takens 的混沌路线修正,但"失稳序列"的框架不变。
统计描述与雷诺分解
把速度分成平均与脉动 $\mathbf{v} = \bar{\mathbf{v}} + \mathbf{v}'$,N-S 方程平均后多出雷诺应力 $-\rho\overline{v_i'v_j'}$——平均运算产生新未知量(闭合问题):湍流模型(混合长、$k$-$\varepsilon$)全是对雷诺应力的模型化尝试。混合长(普朗特):$v'l \sim l\,\partial\bar v/\partial y$,给出湍流边界层的对数速度剖面。
柯尔莫哥洛夫 $-5/3$ 律:本章的皇冠
能量级串的统计理论:大涡注入能量、逐级传给小涡、在耗散尺度 $\eta$ 处化为热。惯性区(介于二者)的能谱:
推导只依赖三个假设:局部均匀各向同性、级串速率 $\varepsilon$、尺度间普适性。湍流中少数被公认的定量结论——实验谱与它的高吻合度是 20 世纪物理的漂亮战绩。
老师的提醒
- 记住湍流的"三段尺":积分尺度 $L$(含能涡)、惯性区、耗散尺度 $\eta \sim L\,Re^{-3/4}$——$Re$ 越大尺度谱越宽,DNS 代价随 $Re^3$ 爆炸。
- 工程上的诚实态度:湍流模型是插值工具不是物理定律,重要设计必须有实验或大涡模拟背书。
自检
- 雷诺应力从哪来、闭合问题是什么?
- 推导 $-5/3$ 律的量纲论证。
- 为什么 DNS 不能覆盖高雷诺数工程问题?
更深一层:湍流的"能量经济学"
湍流的柯尔莫哥洛夫图像是一套能量经济学:大涡是"央行"(从平均流吸储能量)、各级涡是"商业银行"(无耗散地传递)、耗散尺度 $\eta$ 的涡是"消费者"(把动能烧成热)。级串的恒定流量假设(惯性区无生产无耗散)直接产出 $-5/3$ 谱:$\varepsilon$(单位质量耗散率)是唯一货币单位,量纲分析一步到位 $E(k) \propto \varepsilon^{2/3}k^{-5/3}$。实验验证的成功(海洋、大气、风洞谱的 $-5/3$ 段)与间歇性修正的失败(精细统计偏离)并存——1941 K41 是"半步真理":均值谱对,高阶统计需要间歇性修正(K62、多分形)。
湍流模拟的三档价位也因此清晰:RANS(只算均值,模型全部脉动——最便宜、信息最少)、LES(解析大涡、模化小涡——贵 10-100 倍、瞬态结构可信)、DNS(全部解析——科研特权,$Re^3$ 成本)。选档依据:你需要多"瞬态"的信息。
推导细读:混合长理论的"代数裸奔"
普朗特混合长理论的推导是"量纲裸奔"的教科书:湍流脉动速度 $v' \sim l\,|\partial\bar u/\partial y|$(涡携带流体移动混合长 $l$ 所造成的速度差),雷诺应力 $-\rho\overline{u'v'} \sim \rho l^2(\partial\bar u/\partial y)^2$。对壁面湍流设 $l = \kappa y$(与壁距离成正比)+ 近壁剪切恒定 $\tau \approx \tau_w$,一步解出对数律:
对数律从两个假设里直接掉出来——它统治壁面湍流 70 年(从管道到大气边界层),$\kappa$ 与 $B$ 是全书最常用的两个经验常数。混合长的"量纲裸奔"精神(用尺度分析硬算、常数留给实验)贯穿整个湍流工程。
在真实世界里:湍流模型的"信任等级"
工程湍流模型的信任分级:对数律区($\pm 5\%$,可信赖)、分离预测(RANS 模型分歧大,SST 相对可靠)、转捩位置(标准模型不预测,需转捩模型或实验)、曲率/旋流(标准模型系统性偏差,需曲率修正)。"用对模型的场合"比"用好模型"更重要——这是湍流工程的第一诫。
专题深入:湍流的"能量经济学"完整版
湍流的柯尔莫哥洛夫图像是一套能量经济学:大涡是"央行"(从平均流吸储能量)、各级涡是"商业银行"(无耗散地传递)、耗散尺度 $\eta$ 的涡是"消费者"(把动能烧成热)。级串的恒定流量假设(惯性区无生产无耗散)直接产出 $-5/3$ 谱:$\varepsilon$(单位质量耗散率)是唯一货币单位,量纲分析一步到位 $E(k) \propto \varepsilon^{2/3}k^{-5/3}$。实验验证的成功(海洋、大气、风洞谱的 $-5/3$ 段)与间歇性修正的失败(精细统计偏离)并存——1941 K41 是"半步真理":均值谱对,高阶统计需要间歇性修正(K62、多分形)。
湍流模拟的三档价位也因此清晰:RANS(只算均值,模型全部脉动——最便宜、信息最少)、LES(解析大涡、模化小涡——贵 10-100 倍、瞬态结构可信)、DNS(全部解析——科研特权,$Re^3$ 成本)。选档依据:你需要多"瞬态"的信息。
数学细节展开:混合长理论的"代数裸奔"完整版
普朗特混合长理论的推导是"量纲裸奔"的教科书:湍流脉动速度 $v' \sim l\,|\partial\bar u/\partial y|$(涡携带流体移动混合长 $l$ 所造成的速度差),雷诺应力 $-\rho\overline{u'v'} \sim \rho l^2(\partial\bar u/\partial y)^2$。对壁面湍流设 $l = \kappa y$(与壁距离成正比)+ 近壁剪切恒定 $\tau \approx \tau_w$,一步解出对数律:
对数律从两个假设里直接掉出来——它统治壁面湍流 70 年(从管道到大气边界层),$\kappa$ 与 $B$ 是全书最常用的两个经验常数。混合长的"量纲裸奔"精神(用尺度分析硬算、常数留给实验)贯穿整个湍流工程。
工程案例细节:湍流模型的"信任等级"工程版
工程湍流模型的信任分级:对数律区($\pm 5\%$,可信赖)、分离预测(RANS 模型分歧大,SST 相对可靠)、转捩位置(标准模型不预测,需转捩模型或实验)、曲率/旋流(标准模型系统性偏差,需曲率修正)。"用对模型的场合"比"用好模型"更重要——这是湍流工程的第一诫。
例题
大气边界层:$U = 10\ \mathrm{m/s}$、$L = 1000\ \mathrm{m}$、$\nu = 1.5\times10^{-5}$(空气)。估算柯尔莫哥洛夫耗散尺度 $\eta = L\,Re^{-3/4}$,并说明对 DNS 的含义。
解答
DNS 网格数 $\sim (L/\eta)^3 \approx (1.35\times10^6)^3 \approx 2.5\times10^{18}$ 个点——彻底不可计算。
- 这就是"湍流模拟必须建模"的定量理由:尺度谱跨 6 个数量级,全解小涡的算力需求随 $Re^3$ 爆炸。