流体力学 · 从对流到混沌

Rayleigh–Bénard Convection

从微小扰动开始积分三个对流模态。旋转查看冷热交换与示踪轨迹，比较温度剖面、热输运和相空间轨道，看稳态怎样出现，又怎样失去稳定。

模型时间

\tau

$0$

Rayleigh

\mathrm{Ra}

$3287.6$

环流模态

X

$0.1$

板面热流

\mathrm{Nu}_p

$1$

体积热输运

\mathrm{Nu}_v

$1$

初生增长率

\lambda_+

$2.882$

邻近状态距离

\delta

$10^{-6}$

线性稳定性稳态分支稳定

热对流工作台

让扰动自己演化

示踪粒子与轨迹沿同一动态速度场积分；颜色表示所在位置的温度。温度切片同一二维场的三个切片，无展向动力学。同步定量图表热流历史、平均温度剖面与模态相图。相对浮力驱动

r=\mathrm{Ra}/\mathrm{Ra}_c

5.0

动量与热扩散比

\mathrm{Pr}=\nu/\kappa

10.0

初始模态扰动

X_0=Y_0

0.10

每秒播放推进的模型时间

d\tau/dt_{\mathrm{play}}

1.00

播放速率只改变时钟，不改变方程。物理控制会重新播种实验并保留暂停状态。

模型来源：E. N. Lorenz，Deterministic Nonperiodic Flow（1963），方程 23–27。

Lorenz 三模态模型采用自由滑移等温边界， $a^2=1/2,\quad \mathrm{Ra}_c=27\pi^4/4\approx657.5$ 。空间视图为二维场挤出，不是刚性壁面或三维湍流求解。

先抓重点

核心问题: 为什么有的扰动会消失，有的形成稳定环流，还有的让环流不断换向？
一句话直觉: 驱动比 $r$ 决定导热态是否稳定，但失稳不等于扰动凭空出现。“无扰动”仍保持精确导热解；非零扰动才沿不稳定方向增长。

核心数学模型

浮力、阻尼与几何

\mathrm{Ra}=\frac{g\alpha\Delta T H^3}{\nu\kappa},\quad r=\frac{\mathrm{Ra}}{\mathrm{Ra}_c},\quad \mathrm{Pr}=\frac{\nu}{\kappa},\quad \mathrm{Ra}_c=\frac{27\pi^4}{4}

所选水平波数给出自由滑移等温边界的最低初生阈值。刚性壁面不是同一个阈值。

三个动态模态

\frac{dX}{d\tau}=\mathrm{Pr}(Y-X),\quad \frac{dY}{d\tau}=rX-Y-XZ,\quad \frac{dZ}{d\tau}=XY-bZ,\quad b=\frac{8}{3}

环流、冷热对比和平均温度剖面形变一起演化。RK4 使用固定模型时间步长；邻近轨道仅在初始环流上增加极小偏移。

重建不可压缩场

\psi=\frac{\sqrt{2}\kappa(1+a^2)}{a}X\sin\!\left(\frac{\pi ax}{H}\right)\sin\!\left(\frac{\pi z}{H}\right),\quad u=-\partial_z\psi,\quad w=\partial_x\psi,\quad a^2=\frac12

流函数保证散度为零，上下板无穿透；水平方向周期。展向挤出不增加第三维动力学。

不截断温度重建

\frac{T-T_c}{\Delta T}=1-\frac{z}{H}+\frac{\sqrt{2}Y\cos(\pi ax/H)\sin(\pi z/H)-Z\sin(2\pi z/H)}{\pi r}

上下板温度固定。仅颜色映射在板温范围外饱和；低阶模型可能出现温度超调，强驱动下尤其明显。

两种热输运诊断

\mathrm{Nu}_p=1+\frac{2Z}{r},\quad \mathrm{Nu}_v=1+\frac{2XY}{br},\quad \mathrm{Nu}_v-\mathrm{Nu}_p=\frac{2}{br}\frac{dZ}{d\tau}

瞬态储热时，板面温度梯度与体积平均总热流不同；稳态时一致。两者都由当前模态计算，不是贴在动画上的经验曲线。

稳态基准与时间尺度

X_*=Y_*=\pm\sqrt{b(r-1)},\quad Z_*=r-1,\quad \mathrm{Nu}_*=3-\frac{2}{r},\quad \tau=\pi^2(1+a^2)\frac{\kappa t}{H^2}

非零分支在初生之上存在。显示时钟是无量纲模型时间，不是物理秒；强驱动 Lorenz 动力学不是实验室湍流的定量预测。

常见难点

临界值取决于边界

容易误解所有上下加热的流体都有同一个初生阈值。

正确理解这里用自由滑移、等温界面和固定周期；不能把其阈值与刚性无滑移壁面混用。

相图不是流体粒子轨迹

容易误解双翼形模态轨道就是两个真实涡旋。

正确理解相图记录整个流场的模态状态；装置中的示踪粒子才是位置轨迹。它们由同一状态驱动，但处在不同空间。

失稳不保证真实湍流

容易误解稳态分支失稳后，图像就是解析充分的三维湍流。

正确理解高驱动下仅保留三个模态会失去真实流体的高波数结构。混沌换向属于该截断系统，不能用作真实湍流热输运定律。

动手实验

01
让扰动衰减
选择“扰动衰减”，推进几个模型时间单位。
应该观察到： 环流趋近零，平均温度回到线性导热剖面，两种热流趋近一。
02
形成稳定滚轴
选择“稳定滚轴”，比较早期与后期的热流曲线。
应该观察到： 瞬态中两种热流不同；收敛后接近虚线稳态基准，模态轨道停在某个稳态分支。
03
观察确定性换向
选择“Lorenz 换向”，继续播放或提高播放速率。
应该观察到： 环流反复变号，轨道绕两个分支运动，微小初值差可能逐渐放大。播放速率不会修改方程。
04
没有种子会怎样
选择“无扰动”并推进时间，然后给初始扰动一个非零值。
应该观察到： 精确的零解不会自行生出环流；设置种子后，实验从零时刻重新演化。