力学 · 从零学起 第 3 / 5 节

第 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$ 处化为热。惯性区(介于二者)的能谱:

$$E(k) \propto \varepsilon^{2/3}k^{-5/3}$$

推导只依赖三个假设:局部均匀各向同性、级串速率 $\varepsilon$、尺度间普适性。湍流中少数被公认的定量结论——实验谱与它的高吻合度是 20 世纪物理的漂亮战绩。

原文要点:第 3 章强调湍流的两个能量来源(平均流的剪切供给)与两个去处(分子粘性耗散),湍流粘性 $\eta_t \sim \rho v_t l_t$ 作为雷诺应力的量纲估计被引入;书中指出湍流粘性不是材料常数而是流动属性——这是与分子粘性的本质区别。

老师的提醒

  • 记住湍流的"三段尺":积分尺度 $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$,一步解出对数律:

$$\bar u / u_* = \frac{1}{\kappa}\ln\frac{yu_*}{\nu} + B, \qquad \kappa \approx 0.41$$

对数律从两个假设里直接掉出来——它统治壁面湍流 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$,一步解出对数律:

$$\bar u / u_* = \frac{1}{\kappa}\ln\frac{yu_*}{\nu} + B, \qquad \kappa \approx 0.41$$

对数律从两个假设里直接掉出来——它统治壁面湍流 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 的含义。

解答

$$Re = \frac{10\times1000}{1.5\times10^{-5}} \approx 6.7\times10^8$$
$$\eta = \frac{1000}{(6.7\times10^8)^{3/4}} = \frac{1000}{1.35\times10^6} \approx 7.4\times10^{-4}\ \mathrm{m} \approx 0.7\ \mathrm{mm}$$

DNS 网格数 $\sim (L/\eta)^3 \approx (1.35\times10^6)^3 \approx 2.5\times10^{18}$ 个点——彻底不可计算。

  • 这就是"湍流模拟必须建模"的定量理由:尺度谱跨 6 个数量级,全解小涡的算力需求随 $Re^3$ 爆炸。

本页是《流体力学(理论物理教程第六卷)》第 3 节的学习笔记:为便于初学,讲解、例题与要点概括均为本站原创整理,公式与章节定位供对照原书使用。

📑 本书目录 经典研读