第 8 章 · 非平衡热力学(输运现象)· 第 5 份 / 共 5 份

8.6 非平衡态相变

走出近平衡:对流花样、化学振荡与耗散结构
目录(点击展开,也可直接点条目跳转)
    1 / 10
    第零步 · 先看地图

    1. 温差大到一定程度,液体会自己排出花样

    0:00 / 0:00

    教师式重点讲解;切换部分时会自动暂停。

    先看一件事,它不需要任何前置知识:

    取一个平底容器,装一薄层液体(厚度 $d$,几毫米),底下均匀加热。

    温差小的时候:液体宏观上一动不动,热老老实实从底下传到上面,满足傅里叶定律。

    把底下的温度一点点往上调:到某一个温度差,液体突然动起来——而且不是乱动,液面上浮现出一片规则的六角形格子,边界清清楚楚,像蜂巢。再等下去,这个花样也不会散,它就这么稳稳地待着。

    请注意这里最反直觉的一点:我们只是把「推得更狠」,系统却变得更「有条理」了。加热是一件彻底不带任何空间信息的事——底部处处一样热,容器是圆的或方的,没有谁规定六角形该长在哪里。可是花样出现了,而且尺寸、形状都可以算出来。

    这一节讲的就是这类现象。它在这一章里位置特殊:它是前四份成立的那个前提被打破之后的世界。

    前四份一直站在哪块地基上

    回头数一下,第 1 份到第 4 份的每一个结论,都靠三条前提撑着:

    前提在哪里立的
    ①局部平衡假设:每一小块虽然整体不平衡,但小块内部可以当作平衡态,$T$、$p$、$\mu$ 处处有定义第 1 份第 7 段
    ②线性唯象方程:流与力成正比,$J_j=\sum_i L_{ji}X_i$第 1 份第 9 段
    ③昂萨格倒易关系:$L_{ij}=L_{ji}$第 1 份第 10 段、第 2 份

    这三条合起来,就是「近平衡」这个词的全部含义。第 3 份的电动效应、第 4 份的热电效应,全部是在这块地基上盖起来的。

    现在把力加大——温度梯度不再是「小」的。②首先失效。而一旦②失效,③也就无从谈起(倒易关系说的是唯象系数之间的对称,没有唯象系数就没有它)。

    这里要说准一句:失效的不是每一点上的傅里叶定律或黏滞定律——它们照旧成立。失效的是「整个系统的流正比于整个系统的力」这件事,根源是流体的运动方程里本来就藏着非线性项,液体一旦动起来,它就露头了。

    失效之后会发生什么,是这一节唯一的问题。

    朴素的猜测是「更乱」——推得越狠越无序。实际情况恰恰相反:在近平衡区,系统只会老老实实地按比例响应,绝不会自己长出结构;而远离平衡时,系统反而可以从无序变得有序。

    自然界里这类例子非常多:生物由几个化学元素组成蛋白质、再组成有机体;天上的云会排成鱼鳞状、条状的规则结构。本节以液体中的对流图案为例,讲清这件事。

    这一份和前四份的写法不一样

    必须提前说明:非线性问题没有「昂萨格关系」那样的一锤定音的定理。这一节的定量结论(临界值、花样的尺寸)来自对流体力学方程组做线性稳定性分析和数值计算,数学工具远超本书范围。

    所以本份的做法是:结论直接给,但每一个数都会交代它从哪来、凭什么是这个数、怎么和实验对上。你不需要会解那些方程,但读完应该能回答「1707.762 这个数是什么意思」。

    第 2–3 段第 4–5 段第 6–7 段第 8 段第 9 段
    现象:贝纳德对流是什么,长什么样机制:为什么会流起来(两种驱动力)判据:瑞利数 $R$ 与临界值 $R_c$换个舞台:化学振荡统一解释:耗散结构与熵账
    8.1–8.2 · 已讲四条经验定律 → 力与流的框架
    8.3 · 已讲熵产生率:力和流已定准
    8.4 · 已讲电动效应:四个压成两个
    8.5 · 已讲热电效应:三个系数只需一个 ε
    8.6 · 本份走出线性区:非平衡态相变
    回顾一分钟温习第 1 份的两个结论
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    力与流。把「不平衡的程度」写成力 $X_i$(温度不均、浓度不均、电势不均……),把「因此发生的输运」写成流 $J_j$(热流、粒子流、电流……)。

    线性唯象方程。只要偏离平衡不远,每一个流都正比于所有的力:

    $$J_j=\sum_i L_{ji}X_i$$

    对角项 $L_{jj}$ 就是老定律(傅里叶、欧姆、菲克……);非对角项 $L_{ji}$ 是交叉效应——第 3、4 份整两份都在讲它。

    昂萨格倒易关系。$L_{ij}=L_{ji}$。它是白送的一半信息:只要测出一个方向的交叉效应,反方向的不用测也知道。

    关键在「只要偏离平衡不远」这半句。它是一句限制,不是一句客套。本节做的事,就是把这半句去掉,看会发生什么。
    易错崩掉的是②,不是①
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    起点 · 主干刚列出的三条前提
    $$\text{①局部平衡}\qquad \text{②}J_j=\sum_i L_{ji}X_i\qquad \text{③}L_{ij}=L_{ji}$$

    不是。在贝纳德对流里,①仍然成立得好好的。

    液体里每一小块都还有确定的温度、压强、密度;傅里叶定律在每一点上依然成立(热流仍然正比于当地的温度梯度),黏滞定律也一样。液体离「分子层面的平衡」远得很吗?并不——常温常压下的水,分子碰撞把局部平衡维持得非常牢。

    那到底是什么变成非线性的?是方程本身,不是本构关系。流体的运动方程里有一项

    $$(\mathbf{v}\cdot\nabla)\mathbf{v},\qquad (\mathbf{v}\cdot\nabla)T$$

    叫对流项:液体一动,就把自己和自己的热量搬着走。速度乘速度、速度乘温度——这是两个未知量相乘,货真价实的非线性。液体不动($\mathbf{v}=0$)时它是零,所以近平衡时看不见它;一旦动起来,它就是主角。

    所以本节离开的准确说法是:「线性响应」,不是「局部平衡」。

    第 1–5 份里,「流正比于力」是把整个系统看成一个输入输出盒子时成立的;到了这一节,盒子内部自己组织起来了,整体的流和整体的力不再成比例——但盒子里每一点上,老定律还在用。
    第一步 · 现象

    2. 贝纳德的实验:从静止到六角格子

    0:00 / 0:00

    教师式重点讲解;切换部分时会自动暂停。

    把第 1 段那个装置说细一点。

    容器里装一薄层液体,厚度为 $d$,从容器底部均匀加热。记底部温度 $T_0$、液体表面温度 $T_1$,且 $T_1<T_0$。

    温差小:什么都不发生

    当 $T_0-T_1$ 之差较小时,稳恒状态下的热流满足傅里叶定律——热量靠传导一路向上,液体宏观上静止。

    但请留意随之而来的一件事:由于存在热膨胀,与温度梯度相伴随,液体里有一个线性的密度梯度。下面热、所以轻;上面冷、所以重。

    这是一个「头重脚轻」的排布。

    重的在上、轻的在下——从能量上说,只要把上下调换一下,重心就会降低,系统会更「舒服」。那它为什么不立刻翻过来?

    因为有东西按着它。这个「按着它的东西」是什么、什么时候按不住了,就是第 5、6 段的全部内容。阈值的存在,从这里就已经埋下了。

    温差到某个值:突然对流,而且很规则

    1900 年,法国物理学家贝纳德(Bénard)用熔化的鲸蜡作为液体,以 $100\ ^\circ\mathrm{C}$ 的水蒸气在底部加热,发现:

    当上下温差大到一定程度时,液体开始出现对流,且图案非常规则——

    · 在水平面上呈现六角图案的花样;
    · 六角的线性尺寸由液体层的厚度 $d$ 决定(课本写作「约等于 $d$」——第 7 段会把这句话算准:一次翻转的水平距离确实 ≈ $d$,而整个六角的直径是 $d$ 的两三倍);
    · 液体从每个六角的中心向上升,沿半径向外侧流动,再沿六角的边向下流动。

    这个现象就称为贝纳德对流。

    剖面 自由液面 T₁(冷) 底部均匀加热 T₀(热) d 上升与下沉交替排列;间距完全由 d 决定(第 7 段会算出来) 俯视 中心(●)上升 → 沿半径向外流 → 沿六角的边下沉
    图 1 上:剖面。液体分成一个个反向旋转的「胞」,一胞上升、邻胞下降。下:俯视。液面被划成规则的六角格子,中心上升、边缘下沉。

    三件事值得单独记住:

    观察到的为什么值得惊讶
    有阈值:温差小时什么都没有,过了某个值才发生不是「温差越大流得越多」这种连续的响应——是开关式的
    有尺度:胞的大小 ≈ $d$,不是随机的加热是均匀的,容器也没规定尺寸,尺度是系统自己选的
    有对称性破缺:原本水平方向处处等价,现在分出了「中心」和「边」方程和边界条件都是水平各向同性的,解却不是

    第三条是最深的一条。加热、容器、液体本身,在水平方向上都没有任何偏好;可是结果挑出了一套六角格子,把「处处一样」变成了「有中心、有边」。这在物理里有个专门的名字:自发对称性破缺——和铁磁体在居里点以下自己挑一个磁化方向是同一件事。

    旁通同一个机制,换了个场合
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    • 熬粥、煮味噌汤。锅底加热,汤面上会出现一格一格缓慢翻动的花纹,界线相当清楚——这就是贝纳德胞,只是形状被锅的形状和杂质扰乱了。
    • 天上的「鱼鳞云」「条状云」。大气底层被地面加热,上层冷,是一个巨大的贝纳德系统。规则的云街和蜂窝状云(open / closed cells,卫星云图上非常清楚),尺度几公里到几十公里,机制完全相同。
    • 油漆和涂层的「橘皮纹」。溶剂从表面蒸发带走热量、也改变表面张力,涂层里就长出规则的胞,干了以后固定下来——这是工业上真要花力气去避免的贝纳德对流。
    • 太阳表面的「米粒组织」。太阳对流层的顶端就是一层贝纳德胞,每个胞直径约 1000 km,中心亮(热气上升)、边缘暗(冷气下沉)——和图 1 下半张是同一幅画。
    这些例子的尺度从毫米到千公里,跨了九个数量级,却是同一套物理。这正是第 6 段那个无量纲数的价值:它把「多大、多黏、多热」全部折算成一个数,尺度就不重要了。
    易错「不稳定」和「会发生」是两件事
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    如果液体是理想流体(没有黏性、没有热传导),那么只要下面比上面轻,哪怕差一丁点,静止状态就立刻失稳——这叫瑞利-泰勒不稳定性,没有阈值。

    真实液体不是这样,因为有两件事在拖后腿:

    • 黏滞。液元想动,周围的液体就磨它,把动能变成热。动得越快磨得越狠。
    • 热扩散。液元之所以轻,是因为它比周围热。可它一边上升,热量一边往周围漏;温差没了,浮力也就没了。

    所以真正的问题不是「上下重不重」,而是:一个微小的扰动,是会长大还是会衰减?

    浮力让扰动长大,黏滞和热扩散让扰动衰减。两边打成平手的那一点,就是临界点。把这句话写成公式,就是第 6 段的瑞利数。
    第一步 · 现象

    3. 六角还是卷?看两个面是什么面

    0:00 / 0:00

    教师式重点讲解;切换部分时会自动暂停。

    贝纳德那锅鲸蜡有一个容易被忽略的细节:它的上表面是敞着的。

    如果液体表面暴露在空气中,称此界面为自由界面;底部和容器接触,称此界面为刚性界面。两者的差别是实打实的力学差别:

    刚性界面自由界面
    液体在界面上能不能滑动不能。紧贴固体的那层液体速度为零(无滑移)能。液面可以整体被拉着走
    切向的力固体可以给液体任意大的切向阻力只有空气,几乎给不了阻力;但液面自己的表面张力可以拉它
    后果更「刹得住」,需要更大的温差才动更容易动,而且多了一种全新的驱动力

    贝纳德的原始实验是「刚性底 + 自由面」。他看到的六角花样,用硅油拍摄的照片和理论计算的一个六角胞,就是课本图 8.10 的 (a) 和 (b)。

    后来的实验:换成两面都刚性,花样就变了

    后来又做了很多对流实验,绝大部分用上下面均与刚性界面接触,无自由表面(也就是上面盖一块板)。结果是:

    温度差大到一定程度以后会突然出现对流「卷」(roll)——液体卷成一根根平行的滚筒,一根顺时针、相邻一根逆时针。

    · 假如容器是长方形,则对流卷是直的(课本图 8.11(a) 的侧面照片);
    · 如果容器是圆形,则形成环型卷(图 8.11(b) 的同心圆)。

    课本图 8.12 则是理论计算的自由-自由界面的对流卷图案,放在一起以作比较。
    不要把「六角」和「卷」当成两种不同的现象。

    它们是同一个失稳在不同边界条件下选出的不同花样:都是过了阈值以后静止态失稳,都是水平对称性破缺,只是破缺的方式不同——卷只破缺一个方向(沿卷的方向仍然平移不变),六角破缺两个方向。

    粗略地说:有自由表面(表面张力参战)容易出六角,上下都刚性(纯浮力)容易出卷。

    为什么把它叫做「相变」

    课本给了一句很硬的判断:这些图案一经形成就相当稳定,是一个相变。

    「相变」这个词用在这里不是打比方。把它和你熟悉的相变逐条对照:

    相变的特征在贝纳德对流里
    有一个临界值,越过它性质突变有。临界温差 / 临界瑞利数 $R_c$(第 6、7 段)
    对称性破缺:新相的对称性低于旧相有。静止态水平方向连续平移对称;对流态只剩离散的格子对称
    序参量随驱动长大有。序参量是对流速度的振幅 $A$。上下边界对称的卷:$A\propto\sqrt{R/R_c-1}$,从零连续长大(像二级相变);六角:阈值处跳变、带滞后(像一级相变)
    新相稳定,扰动它会自己回来有。这正是「一经形成就相当稳定」
    但和水沸腾那种相变有一个根本差别,必须说清:

    水和水蒸气两相共存,靠的是化学势相等——那是一个平衡态条件,撤掉一切外界作用,冰还是冰。

    而对流花样没有任何平衡态条件可言。你一断掉底部的加热,花样几分钟内就消失得干干净净。它必须靠持续的能量流「养」着才存在。

    所以它叫非平衡态相变——「相变」取的是「阈值+对称破缺+序参量」这一层含义,不是「两相共存」那一层。这也正是本节标题的由来。
    深入一个两行的模型,就能给出这条平方根
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    起点 · 主干表格里那一行
    $$A\propto\sqrt{\dfrac{R}{R_c}-1}$$

    设对流的强弱可以用一个数 $A$(振幅)刻画。它随时间怎么变?先写最粗的一项:

    $$\frac{\mathrm{d}A}{\mathrm{d}t}=\varepsilon A,\qquad \varepsilon\equiv\frac{R}{R_c}-1$$

    $R<R_c$ 时 $\varepsilon<0$,扰动指数衰减;$R>R_c$ 时 $\varepsilon>0$,扰动指数长大——这就是失稳。

    但指数长大不可能一直持续下去,对流强到一定程度会自己限制自己(流动越强,它搬运热量越有效,把驱动它的温度梯度削平了)。补上最低阶的自限制项:

    $$\frac{\mathrm{d}A}{\mathrm{d}t}=\varepsilon A-gA^3,\qquad g>0$$

    为什么是 $A^3$ 而不是 $A^2$:如果把整个花样反过来(上升变下降)也是一个同样合法的解,那么方程必须在 $A\to-A$ 下不变,偶次项被对称性禁掉。

    令 $\mathrm{d}A/\mathrm{d}t=0$:

    $$A\left(\varepsilon-gA^2\right)=0\ \Longrightarrow\ A=0\quad\text{或}\quad A=\sqrt{\varepsilon/g}$$
    $R<R_c$:只有 $A=0$(静止)这一个解,且稳定。
    $R>R_c$:$A=0$ 变得不稳定,新解 $A=\sqrt{\varepsilon/g}\propto\sqrt{R/R_c-1}$ 出现并接管。

    这叫超临界叉式分岔。它和朗道二级相变理论里 $m\propto\sqrt{T_c-T}$ 是同一个数学结构——这就是「相变」这个词用得理直气壮的原因。

    但那个「如果」是有条件的,条件不满足时结论要改。上下边界同种、且液体的物性可以当作不随温度变($A\to-A$ 才真是一个对称操作)时,得到的是卷,上面那套成立。

    而六角不满足:中心上升和中心下沉不等价——这正是六角能分出「中心」和「边」的原因。对称性一破,$A^2$ 项就被允许写进方程:

    $$\frac{\mathrm{d}A}{\mathrm{d}t}=\varepsilon A+\beta A^{2}-gA^{3}$$
    有了 $\beta A^{2}$,分岔就不再是那条平滑的平方根曲线:振幅在阈值处跳变,而且降温回去时走另一条路(滞后)——更像一级相变。

    所以确切的说法是:卷像二级相变,六角像一级相变。两者都是非平衡态相变,只是分岔的类型不同。
    易错照片和计算混在一起,容易读串
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    图界面是什么看到什么
    8.10(a)刚性-自由照片(硅油)一大片六角格子
    8.10(b)刚性-自由理论计算单个六角胞里的流线
    8.11(a)刚性-刚性照片(侧面)方容器里的直卷
    8.11(b)刚性-刚性照片(俯视)圆容器里的环形卷
    8.12自由-自由理论计算一排平行的卷
    留意自由-自由这一档只有计算、没有照片。原因很实在:上下两面都自由的液层,在地面上根本做不出来——总得有个东西托着液体。

    它之所以还被反复计算,是因为它的边界条件最干净,是唯一能手算出解析解的一档。第 7 段那个漂亮的 657.511 就出在这里。
    第二步 · 为什么会流起来

    4. 机制一:把液面拉开的是表面张力

    0:00 / 0:00

    教师式重点讲解;切换部分时会自动暂停。

    理论上做了大量的工作,来理解对流产生的原因。下面给出定性上的解释——两种界面组合,对应两套完全不同的驱动力。先讲贝纳德本人那一种。

    对于底部是刚性界面、而上面是自由界面的情况,形成对流的驱动力是表面张力。

    链条一步一步是这样的:

    第几步发生了什么
    1由于底部温度高,分子运动的速度较大,同时存在涨落
    2某个速度大的液体元偶然到达上层,引起自由表面上一小块面积元的局部加热
    3而液体的表面张力随温度的升高而减小,所以自由面上这块被加热的小面积元,表面张力显著变小
    4于是周围(较冷、张力较大)的表面元把液体从热处向外拉
    5此处液面被拉走、压强减小,于是引起液体从下面向上升来补充
    6上来的液体又是热的 → 回到第 3 步。正反馈闭合,对流就开始了。
    表面张力大的一侧赢,把液面往两边拉走 σ 大 σ 大 底部加热 热点:σ 小 热液体被吸上来
    图 2 表面张力驱动的对流(上边界为与空气接触的自由表面)。热点处 $\sigma$ 小,拉不过周围,液面被向两侧拉开;该处压强降低,下方的热液体补上来,正反馈闭合。
    请注意第 3 步是全链条的枢纽:$\sigma$ 随 $T$ 升高而减小。

    如果表面张力随温度增大,第 4 步的箭头就会反过来——周围的液面被拉向热点,热点处压强升高,液体下沉,扰动被压回去,永远不会有对流。

    整套机制的正负号,全押在 $\mathrm{d}\sigma/\mathrm{d}T<0$ 这一条物性上。

    怎么证明驱动力真的是表面张力,而不是重力?

    课本给了一个很硬的证据。当然这里还有重力的影响(液体在地面上,浮力一直都在),要判断到底谁说了算,最干脆的办法是把重力去掉:

    1972 年在太空船上成功地做了类似的实验,看到六角对称的图案。

    太空船上的重力加速度仅为 $10^{-6}g$——浮力实际上等于零。由此证明驱动力是表面张力,与重力无关。

    这是一个教科书式的控制变量实验:把嫌疑人之一(重力)拿掉,现象照旧,于是它被排除。

    严格一点说,这个实验证明的是「对这样的薄液层,表面张力足以独立驱动对流,不需要重力帮忙」。它并不意味着浮力在哪里都不重要——液层一厚,主角就换人了(第 5 段讲的正是那一支)。

    回顾$\sigma$ 的两种读法,和它的温度趋势
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    两种等价的读法。

    • 力的读法:液面上任画一条线,线两边的液面互相拉扯,单位长度上的拉力就是 $\sigma$,单位 $\mathrm{N/m}$。
    • 能量的读法:要把液面扩大一个单位面积,需要做的功就是 $\sigma$,单位 $\mathrm{J/m^2}$。两个单位是同一个量纲。

    本段用的是力的读法:热点周围的液面「拉得更用力」,所以把热点处的液体拽走。

    为什么随温度升高而减小?表面张力的来源是「表面上的分子比内部少了一半邻居,能量偏高」。温度升高时,分子热运动加剧、液体密度下降、分子间平均距离拉大,表面分子与内部分子的差别被抹平,这份额外能量就变小。

    极端情形最好记:到了临界温度,液相和气相不再有分别,表面根本不存在了,于是 $\sigma\to0$。既然 $\sigma$ 在临界点必须降到零,它一路随温度下降就毫不奇怪了。

    量级:水在室温附近 $\sigma\approx0.072\ \mathrm{N/m}$,$\mathrm{d}\sigma/\mathrm{d}T\approx-1.5\times10^{-4}\ \mathrm{N/(m\cdot K)}$;一般的油约 $-(5\sim8)\times10^{-5}$。
    深入$Ma$ 与 $R$:谁说了算,由液层厚度决定
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    表面张力驱动这一支,判据的写法和第 6 段的瑞利数完全同构——推动量除以两个刹车量:

    $$Ma=\frac{\left(-\dfrac{\mathrm{d}\sigma}{\mathrm{d}T}\right)\Delta T\, d}{\rho\,\nu\,D_T}$$

    称马兰戈尼数,临界值约 $Ma_c\approx80$。它和 $R$ 的差别只在两处:推动量换成了表面张力随温度的变化率,而厚度只出现一次方(因为驱动力只作用在表面这一层),不是 $d^3$。

    两支的强弱之比,把 $\Delta T$、$\nu$、$D_T$ 全约掉了:

    $$\frac{Ma}{R}=\frac{-\mathrm{d}\sigma/\mathrm{d}T}{\rho\,g\,\alpha\,d^{2}}\ \propto\ \frac{1}{d^{2}}$$
    结论只和厚度有关:液层越薄,表面张力越占优。

    代入一般油类的量级($-\mathrm{d}\sigma/\mathrm{d}T\approx6\times10^{-5}\ \mathrm{N/(m\cdot K)}$,$\rho\approx900$,$\alpha\approx8\times10^{-4}$),令 $Ma/Ma_c$ 与 $R/R_c$ 相等,解出交叉厚度约 1 cm:

    · $d\lesssim1\ \mathrm{cm}$:表面张力主导(贝纳德那层约 1 mm 的鲸蜡,正落在这一档,而且优势约 100 倍);
    · $d\gtrsim1\ \mathrm{cm}$:浮力主导(大气、地幔、太阳对流层,全都远在这一侧)。

    历史上这是一桩公案:贝纳德 1900 年的实验,长期被当作浮力对流的典范;直到 1956 年布洛克(Block)和 1958 年皮尔逊(Pearson)才指出,那层薄薄的液体里真正说了算的是表面张力。所以今天的术语把两者分开:瑞利-贝纳德对流指浮力那一支,贝纳德-马兰戈尼对流指表面张力这一支。

    第二步 · 为什么会流起来

    5. 机制二:浮力,以及两个把它按住的东西

    0:00 / 0:00

    教师式重点讲解;切换部分时会自动暂停。

    换一种边界条件:对于上下两面均为刚性界面,对流是由重力引起的。上面盖了板,没有自由表面,表面张力这条路被堵死了,剩下的只有浮力。

    课本的链条同样是三步:

    第几步发生了什么
    1当容器底部加热时,较下面的液体层受热膨胀,产生一个浮力,企图向上升
    2但液体的黏滞力将阻止它上升
    3当底部的温度进一步升高时,浮力将克服黏滞力,而上面较冷的液体将向下运动,从而形成宏观的对流

    「浮力克服黏滞力」这句话是全节的转折点:它把一个定性的故事,变成了一个可以比大小的不等式。既然是比大小,就一定能写成一个比值——这个比值就是下一段的主角。

    刹车其实有两个,不是一个

    课本第 2 步只提了黏滞力。但真正让阈值出现的,是两件事一起在拖后腿:

    底部加热 → 下层液体膨胀、变轻 热液元 浮力(推动) 黏滞:把动能磨掉 热扩散:把它和周围的 温差抹平 → 浮力也没了
    图 3 浮力驱动的三方角力。浮力推它上去;黏滞把它的动能磨掉;热扩散把它和周围的温差抹平,连浮力本身都一起抹掉。
    它做的事它拖慢的是
    浮力热液元比周围轻,被往上推——(这是推动方)
    黏滞周围液体磨它,动能变成热耗散掉它动的能力
    热扩散热量从液元漏向周围,温差被抹平它轻的能力(浮力的来源)
    第二个刹车比第一个更隐蔽,但同样致命。

    液元之所以有浮力,唯一的理由是它比周围热。它上升需要时间;在这段时间里,热量一直在往周围漏。如果漏得比升得快,等它升到一半,它已经和周围一样冷了——浮力自己消失了。

    所以这里比的是两个速率:热扩散抹平温差的快慢,和浮力驱动液元上升的快慢。这也解释了为什么下一段的判据里,热扩散系数会和黏滞系数并排出现在分母上。
    易错问题不是「谁轻谁重」,是「扰动能不能长大」
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    「热的轻所以往上浮」这句话,作为对流的方向说明是对的,但它解释不了本节要解释的两件事:

    • 为什么有阈值。按这句话,只要下面比上面热一点点,就该马上流起来。可实验说不会。
    • 为什么有确定的尺寸。这句话里根本没有长度,凭什么胞的宽度是 $d$ 而不是别的?

    把问法换掉,两个问题一起解决:不问「静止状态合不合理」,而问——

    给静止状态加一个任意小的扰动,它是长大还是衰减?

    这个问法自动带来两样东西:
    ① 长大和衰减之间有一条分界线,那就是阈值;
    ② 扰动可以有各种波长,不同波长长大的难易不同,最容易长大的那一个波长胜出——这就是尺寸的来源。

    这套问法叫线性稳定性分析:它不解非线性方程,只问「零附近的小扰动往哪个方向走」。第 7 段那三个精确到小数点后三位的临界值,就是这么算出来的。

    顺带回答一个常见困惑:为什么会有「向下运动的液体」?

    液体不可压缩,有上升就必须有等量的下沉。这就是为什么对流一定成「胞」——不可能整层液体一起往上跑。图 1 剖面里一升一降相邻排列,根源就在这条。
    第三步 · 判据

    6. 把「浮力克服黏滞力」写成一个数

    0:00 / 0:00

    教师式重点讲解;切换部分时会自动暂停。

    第 5 段结束在一句话上:浮力要克服黏滞力。要把它变成能算的东西,就得把「推动」和「阻碍」各自量化,然后相除。对此瑞利建立了一个模型,出现对流的条件用瑞利数 $R$ 来衡量,定义为

    $$R=\frac{g\,\alpha\,\Delta T}{\eta\,D_T}\cdot d^{3}$$
    (8.6.1)

    逐项读一遍。课本把它分成两组:

    组符号是什么为什么在这一边
    使对流发生
    (分子)
    $g$重力加速度没有重力就没有浮力(第 4 段的太空船实验)
    $\alpha$热膨胀系数热了要能变轻才有浮力;$\alpha=0$ 的液体绝不对流
    $\Delta T$上下底之间的温度差温差越大,上下密度差越大
    阻止对流发生
    (分母)
    $\eta$黏滞系数磨掉液元的动能
    $D_T$热扩散系数抹平液元与周围的温差,把浮力本身抹掉
    几何$d^{3}$液体的厚度见下面

    分子那三个正是第 5 段的「推动」,分母那两个正是那两个刹车。$R$ 就是它们的比值:$R$ 大 = 推得过;$R$ 小 = 按得住。

    (8.6.1) 里的 $\eta$,必须理解成「运动黏度」$\nu=\eta_{\text{动力}}/\rho$,单位 $\mathrm{m^2/s}$。

    这一点课本没有点明,但绕不过去:只有取运动黏度,$R$ 才是无量纲的;若按中学和普物里常用的动力黏度($\mathrm{Pa\cdot s}$)代入,算出来的 $R$ 带着 $\mathrm{m^3/kg}$ 的量纲,和 1707 这种纯数根本没法比大小。

    两者只差一个密度:$\nu=\eta_{\text{动力}}/\rho$,而水在室温下这两个数相差整整 1000 倍。用数值算题时若发现结果差了三四个数量级,十有八九就是栽在这里。

    它是一个纯数,这一点比它的形式更重要

    $R$ 没有单位。这意味着比较不同的系统时,只需要比这一个数:一锅毫米厚的油、一杯水、几公里厚的大气、几千公里厚的地幔,尺度跨了九个数量级,只要 $R$ 相同,行为就相同。这也是为什么实验室里几毫米厚的一层硅油,能拿来研究地幔里的对流。

    于是判据可以写得非常干脆。当 $R<R_c$(临界瑞利数)时,液体中发生的是线性的热传导过程;当 $R>R_c$ 时,出现对流,并按不同的边界条件出现不同的规则图案。

    $R<R_c$$R>R_c$
    液体的状态宏观静止规则对流
    热怎么走纯导热,傅里叶定律导热 + 对流搬运
    属于本章哪一半第 1–4 份的线性区本节的非线性区

    $d^{3}$:为什么厚一点差别就那么大

    厚度进入的是三次方,这是 (8.6.1) 里最有威力的一项:层厚翻一倍,同样的温差算出的 $R$ 变成 8 倍;反过来说,临界温差降到原来的 $1/8$。

    所以一杯稍微深一点的水,一丁点温差就已经在悄悄翻腾;而上下封死的薄液层想靠浮力对流,难得很——1 mm 厚的水要差 120 K。

    但别把这条用到贝纳德那锅蜡上。

    $d^{3}$ 说的是浮力这一支。贝纳德的液层上面是自由面,说了算的是表面张力(第 4 段),那一支里厚度只是一次方——1 mm 的液层大约差一度就动了,一点也不难。

    他之所以要用 $100\ ^\circ\mathrm{C}$ 的蒸汽,是因为鲸蜡在四五十摄氏度才熔化,得先让它化开并保持液态,不是为了凑出上百度的温差。
    易错三行量纲验算,一次记牢
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    起点 · 主干的式 (8.6.1),原样照抄
    $$R=\frac{g\,\alpha\,\Delta T}{\eta\,D_T}\cdot d^{3}$$

    分子的量纲:

    $$[\,g\,\alpha\,\Delta T\,d^{3}]=\frac{\mathrm{m}}{\mathrm{s^2}}\cdot\frac{1}{\mathrm{K}}\cdot \mathrm{K}\cdot \mathrm{m^3}=\frac{\mathrm{m^4}}{\mathrm{s^2}}$$

    要让整个 $R$ 无量纲,分母必须也是 $\mathrm{m^4/s^2}$。热扩散系数是死的:$[D_T]=\mathrm{m^2/s}$。所以

    $$[\eta]=\frac{\mathrm{m^4/s^2}}{\mathrm{m^2/s}}=\frac{\mathrm{m^2}}{\mathrm{s}}$$

    $\mathrm{m^2/s}$ 正是运动黏度 $\nu$ 的单位。

    若误用动力黏度($\mathrm{Pa\cdot s}=\mathrm{kg/(m\cdot s)}$):

    $$\frac{\mathrm{m^4/s^2}}{\dfrac{\mathrm{kg}}{\mathrm{m\cdot s}}\cdot\dfrac{\mathrm{m^2}}{\mathrm{s}}}=\frac{\mathrm{m^3}}{\mathrm{kg}}\ \neq\ 1$$
    两者只差一个密度:$\nu=\eta_{\text{动力}}/\rho$。水在 $20\ ^\circ\mathrm{C}$:$\eta_{\text{动力}}\approx1.0\times10^{-3}\ \mathrm{Pa\cdot s}$,$\rho\approx10^{3}$,故 $\nu\approx1.0\times10^{-6}\ \mathrm{m^2/s}$——整整差了 1000 倍,算出来的 $R$ 也就差 1000 倍。

    记法:分母上两个量必须是同一种东西——$\nu$ 是「动量的扩散系数」,$D_T$ 是「热量的扩散系数」,它们的比 $\nu/D_T$ 叫普朗特数,本来就是个纯数。能和 $D_T$ 并排站的,只能是 $\nu$。
    实例把水的数代进 (8.6.1),看 $d^{3}$ 有多凶
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    起点 · 主干的式 (8.6.1) 与判据 $R>R_c$
    $$R=\frac{g\,\alpha\,\Delta T}{\nu\,D_T}\,d^{3}\ >\ R_c$$

    水在 $20\ ^\circ\mathrm{C}$:$\alpha=2.07\times10^{-4}\ \mathrm{K^{-1}}$,$\nu=1.00\times10^{-6}\ \mathrm{m^2/s}$,$D_T=1.43\times10^{-7}\ \mathrm{m^2/s}$,取 $g=9.8\ \mathrm{m/s^2}$。上下都用刚性面,$R_c=1707.8$(第 7 段)。

    先算「每 1 K 温差贡献多少 $R$」:

    $$\frac{R}{\Delta T}=\frac{9.8\times2.07\times10^{-4}}{1.00\times10^{-6}\times1.43\times10^{-7}}\,d^{3}=1.42\times10^{10}\,d^{3}\ \ \mathrm{K^{-1}}\quad(d\ \text{以米计})$$
    层厚 $d$$R/\Delta T$(每 K)临界温差 $\Delta T_c$
    1 mm14.2120 K(实际上做不到,水早沸了)
    2 mm11315 K
    5 mm17730.96 K
    1 cm141860.12 K
    厚度只涨了 10 倍,临界温差降了 1000 倍——这就是 $d^{3}$。

    最后一行值得回味:一杯 1 cm 深的水,上下差十分之一度就已经在对流了。所以「静止的水靠传导散热」这件事,在日常尺度上几乎从不成立;反过来,贝纳德必须把液层做到毫米级、并且烧到 $100\ ^\circ\mathrm{C}$,才能在可控的条件下把阈值卡出来。
    列传瑞利数 $R$
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    起点 · 主干第一次写下它的地方,原样照抄
    $$R=\frac{g\,\alpha\,\Delta T}{\eta\,D_T}\cdot d^{3}$$
    幕一 · 身世:没有它,哪句话说不出来

    「浮力克服黏滞力」——这句话没法用。浮力和黏滞力的量纲相同,但它们各自依赖液体厚度、温差、材料,把它们分别算出来再比大小,换一种液体、换一个容器就得重算一遍,而且永远得不到一个能写进教科书的判据。

    更糟的是:实验上换了液体、换了厚度,临界温差从零点几度跳到上百度,看起来毫无规律。没有 $R$,这些实验就是一堆互相对不上的数据。有了它,全部塌缩成一句话:过 1708 就对流。

    幕二 · 名片
    读法
    瑞利数,Rayleigh number;英文读 /ˈreɪli/(「瑞-利」,不是「瑞雷」)
    量纲
    无量纲,纯数,没有单位
    记号
    课本写 $R$;文献里更常见 $Ra$,避免和电阻、气体常数撞车
    亲戚
    雷诺数 $Re$、马兰戈尼数 $Ma$、普朗特数 $Pr$、努塞尔数 $Nu$——流体力学里一大家子无量纲数,个个都是「推动 ÷ 阻碍」
    幕三 · 脾气:它平时怎么表现
    • 只有正值有意义。$\Delta T$ 指下热上冷;反过来(上热下冷)$R<0$,那是稳定分层,永远不对流。
    • 它不是「有多少对流」,而是「够不够格对流」。$R=1706$ 和 $R=1$ 一样是静止的;跨过 1708 才有质变。它是一个阈值型的量。
    • 对厚度极其敏感($\propto d^{3}$),对温差只是一次方。想让系统对流,加厚比加热有效得多。
    • 它没有上限。地幔 $R\sim10^{7}$,太阳对流层 $R\sim10^{20}$。$R$ 远超 $R_c$ 时规则花样会失稳、破碎,最终变成湍流——本节讲的是刚过阈值那一小段。
    幕四 · 去向:以后还会在哪儿碰见它

    本份第 7 段会算出它的三个临界值并和实验对照,第 9 段用它来说明「花样必须靠持续供能养着」。走出本书,它是气象、海洋、地球物理和恒星内部对流的标准判据;工程上换热器、晶体生长、建筑通风的设计里,它是必查的第一个数。

    一句话认人:$R$ 是「想翻腾的劲」除以「按得住的劲」,比值过了一千七百零八,液体就不肯再老实待着了。
    第三步 · 判据

    7. 三个临界值,和它们与实验的对照

    0:00 / 0:00

    教师式重点讲解;切换部分时会自动暂停。

    $R_c$ 不是量出来的,是算出来的。理论上做了大量的计算工作,主要是分析一组平衡方程:

    方程管的是
    有局域质量密度平衡方程液体不会凭空增减(不可压缩)
    包括速度的动量平衡方程浮力、压强、黏滞力怎么推动液体
    内能的平衡方程热怎么被传导、又怎么被流动带着走

    根据不同的边界条件,计算出何时发生出现对流图案的相变,以及对流图案的花样。前面的图 8.10(b) 就是计算得到的刚性-自由界面的一个六角的对流花样,而图 8.12 是对自由-自由界面计算得到的对流卷图案。

    三个数

    理论上计算的瑞利数为

    $$\begin{aligned} R_c&=657.511 &&\text{(对自由-自由界面)}\\ R_c&=1100.650 &&\text{(对自由-刚性界面)}\\ R_c&=1707.762 &&\text{(对刚性-刚性界面)} \end{aligned}$$

    实验上得到的瑞利数 $R_c=1700\pm51$,可见符合相当好。

    这个对照要看两层意思。

    ① 数值对得上:$1707.762$ 落在 $1700\pm51$ 的区间里(偏离仅 $0.15$ 个标准差)。

    ② 挑对了那一档才对得上:实验值该和三个数里的哪一个比?和 $1707.762$——因为第 3 段说过,绝大部分对流实验是上下面均为刚性界面。理论和实验能对上,前提是边界条件也对上了。

    三个数的大小顺序本身也有物理含义:$657<1101<1708$,刚性面越多,需要的推动越大。因为刚性面上液体不能滑动(速度必须为零),液体要在贴壁处从零加速起来,速度梯度大,黏滞耗散就大——每多一个刚性面,就多一份刹车。

    657.511 可以手算出来

    三档里,自由-自由那一档的边界条件最干净,能得到解析解。理论给出的中性稳定曲线(即「扰动既不长大也不衰减」所对应的 $R$)是

    $$R(a)=\frac{(\pi^{2}+a^{2})^{3}}{a^{2}}$$

    其中 $a$ 是水平方向扰动的无量纲波数(以液层厚度 $d$ 为长度单位)。这条曲线是一条 U 形曲线:

    扰动波数 a(无量纲) R 最低点 657.5 a꜀=2.22 曲线以上:扰动长大 → 对流 曲线以下:扰动衰减 → 纯导热
    图 4 自由-自由界面的中性稳定曲线。曲线以上的扰动会长大,以下的会衰减。最先失稳的是曲线最低点对应的那个波数,最低点的高度就是 $R_c$。

    为什么是最低点?因为把 $R$ 从小往大调,第一个碰到曲线的地方就是最低点——那个波数的扰动最先长大,它挑走了整个花样的尺寸。所以:

    $$R_c=\min_a R(a),\qquad a_c=\arg\min_a R(a)$$

    令 $x=a^{2}$,求极小:

    $$\frac{\mathrm{d}}{\mathrm{d}x}\frac{(\pi^{2}+x)^{3}}{x}=\frac{(\pi^{2}+x)^{2}\left[3x-(\pi^{2}+x)\right]}{x^{2}}=0\ \Longrightarrow\ 2x=\pi^{2}$$

    即 $a_c=\pi/\sqrt2\approx2.221$。代回去:

    $$R_c=\frac{\left(\tfrac{3}{2}\pi^{2}\right)^{3}}{\tfrac{1}{2}\pi^{2}}=\frac{27}{4}\pi^{4}=657.511\cdots$$
    课本上那个看起来像是数值计算出来的 $657.511$,其实是一个闭式:$\dfrac{27}{4}\pi^{4}$。

    另外两档没有这么好的运气——刚性边界让本征值问题变成超越方程,$1100.650$ 和 $1707.762$ 只能数值求解。这也解释了课本为什么把三个数都写到小数点后三位:它们是本征值,不是测量值,想要多少位就有多少位。

    顺带就把「花样有多大」算了出来

    三档的临界波数分别是 $a_c=2.221$、$2.682$、$3.117$(都是浮力这一支的值)。波数对应的水平波长(以 $d$ 为单位)是 $\lambda=2\pi/a$,于是 $\lambda=2.83\,d$、$2.34\,d$、$2.02\,d$——都是两三倍的层厚。

    从波长换算到「一个胞有多大」,卷和六角要分开算:

    • 卷:一上一下才是一个完整周期,所以一根卷只占半个波长,宽度 $=\lambda/2$。
    • 六角:它是三个夹 $120^\circ$ 的波叠出来的,上升中心排成一个三角格子,相邻两个中心的距离不是 $\lambda$,而是 $\dfrac{2}{\sqrt3}\lambda\approx1.15\,\lambda$。
    界面$\lambda$卷宽 $\lambda/2$六角直径 $2\lambda/\sqrt3$
    自由-自由$2.83\,d$$1.41\,d$$3.27\,d$
    自由-刚性$2.34\,d$$1.17\,d$$2.70\,d$
    刚性-刚性$2.02\,d$$1.01\,d$$2.33\,d$
    回头把第 2 段那句「六角的线性尺寸约等于 $d$」读准。

    · 如果「线性尺寸」指整个六角的直径(相邻两个上升中心的距离),那是 $2\sim3$ 倍的 $d$,不是 $d$;
    · 如果指液体一次翻转所走的水平距离(从中心的上升流到边上的下沉流),那是上面数字的一半,$1.2\sim1.6\,d$,这才是 ≈ $d$。

    对卷来说两种读法合一,都是 $\approx d$(刚性-刚性 $1.01\,d$,一根卷的截面几乎是正方形)。所以课本那句话按「一次翻转的距离」理解才准确。

    还有一处要交代清楚:上面三个 $a_c$ 属于浮力那一支。贝纳德本人那个实验是表面张力驱动的(第 4 段),它有自己的临界波数 $a_c\approx1.99$,给出的波长 $3.16\,d$、六角直径 $3.65\,d$——比浮力那一支更大,也和实验上看到的、明显宽于层厚的胞径相符。

    但无论按哪种读法、哪一支机制,结论不变:花样的尺度只由一个长度决定,就是液层厚度 $d$。

    这才是本节最值得带走的一件事。

    加热是均匀的,容器没有规定任何尺寸,可是花样的大小被算了出来,而且只依赖一个长度——液层厚度 $d$。

    系统自己选择了一个尺度。这就是「自组织」这个词的确切含义:结构的尺寸不是外界强加的,是系统在众多可能的波长里,挑了最容易长大的那一个。
    深入不解方程,只问「小扰动往哪走」
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    起点 · 主干的中性曲线
    $$R(a)=\frac{(\pi^{2}+a^{2})^{3}}{a^{2}}$$

    完整推导要用到流体力学,这里只交代它是怎么被逼出来的,三步:

    第一步:把解写成「静止解 + 小扰动」。静止的纯导热态是方程组的一个精确解(速度处处为零、温度线性分布)。设速度、温度各偏离它一点点,代入方程,把所有二次以上的小量丢掉——包括那个惹事的对流项 $(\mathbf v\cdot\nabla)$,因为它是两个小量相乘。剩下的是一组线性方程。

    第二步:按波数分解。线性方程的好处是不同波长的扰动互不干扰,可以一个一个算。于是设扰动形如

    $$\propto \sin(a x)\,\sin(n\pi z)\,e^{st}$$

    水平方向按波数 $a$ 分解,竖直方向要满足上下边界条件(自由-自由时正好是正弦,这就是这一档能手算的原因),$s$ 是增长率。

    第三步:令 $s=0$。$s>0$ 扰动长大、$s<0$ 衰减,分界线是 $s=0$。把 $s=0$ 代进去,方程给出 $R$ 与 $a$ 的一个关系——就是那条 $R(a)$。竖直方向取最低的模式 $n=1$(能量最省),得到主干那个式子。

    注意这套方法的边界:它只能告诉你「什么时候开始动」,不能告诉你「动起来以后有多强」。

    因为第一步就把非线性项扔了。$R$ 刚过 $R_c$ 时扰动确实很小,答案可信;$R\gg R_c$ 时对流已经很强,被扔掉的项重新变成主角,这套分析就失效了。

    「动起来以后有多强」要靠第 3 段深入框里那个 $A^{3}$ 项——把非线性捡回最低的一阶,才算得出振幅。
    旁通一个数顶一整套实验
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    $R$ 不是孤例,它有一大家子亲戚,构造方式完全一样:

    数推动 ÷ 阻碍阈值意味着
    雷诺数 $Re=\dfrac{vL}{\nu}$惯性 ÷ 黏滞层流 → 湍流
    瑞利数 $R$浮力 ÷(黏滞 × 热扩散)导热 → 对流
    马兰戈尼数 $Ma$表面张力梯度 ÷(黏滞 × 热扩散)静止 → 表面驱动对流
    马赫数 $M=\dfrac{v}{c}$流速 ÷ 声速亚声速 → 超声速

    为什么非要凑成无量纲?因为物理规律不该依赖你用米还是用厘米。一个只含无量纲组合的判据,换单位不变、换尺度也不变。

    这带来一个极实用的后果:相似原理。想知道一架飞机的空气动力性能,不必造出真飞机——做个小模型,把 $Re$ 调成一样,风洞里的流动和真机就是同一个流动。

    同理,实验室里几毫米厚的硅油层,只要 $R$ 相同,就能拿来研究几千公里厚的地幔。这就是第 2 段那些跨九个数量级的例子能被放在一起讲的根据。
    第四步 · 不只是流体

    8. 化学反应也会「排出花样」——只不过排在时间上

    0:00 / 0:00

    教师式重点讲解;切换部分时会自动暂停。

    到这里为止,所有例子都是流体。如果非平衡态相变只发生在液体里,它就只是流体力学的一个专题。另一个例子是化学反应,它把这件事的普遍性一下子推开了。

    一些化学系统在很好定义的条件下,若干反应物的浓度呈现很复杂的随时间的振荡。

    而且这个振荡:可以重复;与初始条件无关;对小的外界扰动是稳定的。

    请对照第 3 段的清单看这三条——它们和「对流花样一经形成就相当稳定」是同一件事,只不过贝纳德的结构长在空间里,化学振荡的结构长在时间里。

    BZ 反应

    研究最多的是 Belousov-Zhabotinskii 反应(贝洛乌索夫-扎鲍廷斯基反应,简称 BZ 反应)。配方是:

    成分角色
    溴酸钾 $\mathrm{KBrO_3}$氧化剂,反应的「燃料」之一
    丙二酸 $\mathrm{CH_2(COOH)_2}$被氧化的底物,真正被消耗掉的那个
    硫酸铈催化剂,$\mathrm{Ce^{3+}}\leftrightarrow\mathrm{Ce^{4+}}$ 来回倒
    硫酸 $\mathrm{H_2SO_4}$溶剂兼强酸环境

    把它们溶解于硫酸中,搅拌溶液,颜色周期性地来回变化(跟着 $\mathrm{Ce^{3+}}$、$\mathrm{Ce^{4+}}$ 的比例走),周期从几十秒到几分钟,肉眼可见,可以持续几十个周期。课本描述的「红蓝之间振荡」是演示时最常见的样子——不过那抹鲜艳的红与蓝其实来自另加的指示剂,纯铈体系本身只是「无色 ↔ 淡黄」。

    lg[Br⁻] lg([Ce⁴⁺]/[Ce³⁺]) 时间 t 诱导期(不振荡) 一个周期
    图 5 BZ 反应的示意曲线(对应课本图 8.13)。两条曲线反相:$\mathrm{Br^-}$ 缓慢下降,一旦跌破阈值,$\mathrm{Ce^{4+}}$ 猛增;随后 $\mathrm{Br^-}$ 被重新造出来,把反应压回去。注意不是正弦波——是「慢慢积累、突然放电」的弛豫振荡。
    课本图 8.13 有两个细节值得注意,不看清会读错:

    ① 纵轴是 $\lg[\mathrm{Br^-}]$ 和 $\lg\left([\mathrm{Ce^{4+}}]/[\mathrm{Ce^{3+}}]\right)$,都取了对数。$\mathrm{Br^-}$ 的浓度在一个周期里变化三个数量级以上,不取对数根本画不进一张图。

    ② 图的左段(约前 600 s)没有振荡,之后才开始。这段叫诱导期:反应刚开始时体系还没被推到足够远离平衡的地方,等某些中间产物积累够了,振荡才启动。「远离平衡」是有门槛的,这张图把门槛画出来了。

    什么样的系统才振荡得起来

    课本把条件列得很明确:产生化学振荡的条件很苛刻,要满足:远离平衡,对物质源及能量源要开放,反应机制必须复杂(没有一个简单系统有振荡),速率方程必须包含非线性项。

    条件为什么非要不可
    远离平衡平衡态附近一切都单调地趋向平衡;细致平衡原理直接禁止了持续的循环
    对物质源、能量源开放振荡要靠消耗丙二酸维持。封住系统,燃料烧完,振荡必然停下来——BZ 在烧杯里也只能振几十个周期
    机制复杂要有自催化(产物促进自己生成)和反馈抑制两个环节,缺一不可
    速率方程含非线性项这是数学上的硬要求:只有非线性系统才可能有「振幅由自己决定、被扰动后能自己回来」的持续振荡

    把这四条和贝纳德那一套对照,会发现是同一张清单:远离平衡($R>R_c$)、持续供能(底部一直加热)、有正反馈(第 4 段那条闭合链)、方程非线性(对流项)。两个毫不相干的系统,长出结构所需要的条件完全一样。这正是下一段要收的网。

    易错红蓝来自指示剂,铈本身是「无色 ↔ 淡黄」
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    纯铈体系的颜色对比很弱:$\mathrm{Ce^{3+}}$ 无色,$\mathrm{Ce^{4+}}$ 淡黄。肉眼能看出振荡,但绝谈不上「红蓝分明」。

    课堂演示和网上的视频里那种鲜艳的红 ↔ 蓝,来自额外加入的氧化还原指示剂「邻二氮菲亚铁」(ferroin):还原态是红色,氧化态是蓝色。它跟着铈的价态一起翻转,把一个微弱的信号放大成了强烈的视觉效果。

    这不影响课本任何一句话的物理——振荡的主角仍然是 $\mathrm{Ce^{3+}}/\mathrm{Ce^{4+}}$ 和 $\mathrm{Br^-}$,图 8.13 画的也确实是铈的价态比。

    指示剂只是一块表盘。但如果你照着课本的配方去做实验,只放铈盐而不加 ferroin,然后疑惑「说好的红蓝呢」,那就白折腾了。
    深入会振荡,但振不出「振幅由自己定」的那一种
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    起点 · 主干那四个条件里的最后一条
    $$\text{速率方程必须包含非线性项}$$

    先承认一半:线性系统当然可以振荡——弹簧振子就是。但线性系统的振荡只有三种命运:

    • 有阻尼 → 振幅指数衰减,最后停下;
    • 负阻尼 → 振幅指数长大,最后炸掉;
    • 零阻尼 → 一直振,但振幅完全由初始条件决定,推它一下,新振幅就换了,永远回不来。

    而课本要求的振荡是:与初始条件无关、对小扰动稳定。三种命运没有一种做得到。

    能做到的只有「极限环」:相平面上一条孤立的闭合轨道,附近所有轨道都被它吸过去。

    被推离了会自己回来,起点不同也殊途同归——这正是「振幅由系统自己定」的几何图像。而极限环只能出现在非线性系统中:自治线性系统的闭轨要么一个没有,要么是一整族同心闭轨(没有哪一条是「孤立」的、能吸引邻居的)。

    (「自治」二字不能省:外加一个周期驱动的线性系统,比如受迫阻尼振子,确实有唯一稳定的周期解。但那个振幅和频率是外界给的,不是系统自己定的——BZ 反应没有人在旁边给它打拍子。)

    BZ 反应里的非线性项在哪?核心是溴的自催化步骤:亚溴酸 $\mathrm{HBrO_2}$ 参与反应,结果生成两份 $\mathrm{HBrO_2}$——

    $$\mathrm{BrO_3^-}+\mathrm{HBrO_2}+2\mathrm{Ce^{3+}}+3\mathrm{H^+}\ \longrightarrow\ 2\,\mathrm{HBrO_2}+2\mathrm{Ce^{4+}}+\mathrm{H_2O}$$

    它对 $[\mathrm{HBrO_2}]$ 的贡献是 $+k[\mathrm{HBrO_2}]$——越多长得越快,指数爆炸。刹住它的是自身的歧化 $-2k'[\mathrm{HBrO_2}]^{2}$,这是一个平方项,货真价实的非线性。

    而 $\mathrm{Br^-}$ 是抑制剂:它在的时候会把 $\mathrm{HBrO_2}$ 吃掉,自催化点不着火。一个周期是这样走的:$\mathrm{Br^-}$ 被慢慢消耗 → 跌破阈值 → 自催化点火、$\mathrm{Ce^{4+}}$ 猛增(颜色翻转)→ $\mathrm{Ce^{4+}}$ 氧化丙二酸,重新放出 $\mathrm{Br^-}$ → 抑制剂回来,火被扑灭 → 从头再来。

    「自催化 + 延迟的负反馈」——这就是全部秘密。对照第 4 段:热点让表面张力变小(正反馈),液体流走又把热带走(负反馈)。两个系统,同一个骨架。
    溯源为什么当年的化学家不肯相信
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    1951 年,苏联生物化学家别洛乌索夫(Б. П. Белоусов)在研究一个模拟三羧酸循环的体系时,看到溶液颜色周期性地来回变化。他投出论文,被拒,理由是「明显违反热力学第二定律」;改投另一家,再次被拒。他最终只在一份不起眼的会议摘要集上发表了它,此后离开了这个课题。

    直到 1961 年前后,扎鲍廷斯基重做并系统研究了这个反应,它才被接受。1972 年 Field、Kőrös、Noyes 给出了详细的反应机理(FKN 机理),把那个自催化步骤揪了出来。

    当年那句拒稿理由,错在哪?

    「化学反应只能单调地趋向平衡」这条判断,对封闭的、接近平衡的体系是对的。审稿人的错误在于把它当成了无条件的定律。

    BZ 反应整体上确实一路朝平衡走——丙二酸在被不可逆地消耗,自由能在单调下降,熵产生率始终为正。振荡发生在「走向平衡的路上」,而不是绕过了平衡。燃料一旦耗尽,振荡立刻停止,体系安安静静躺在平衡态。

    第二定律要求的是终点,从来没要求过路径必须单调。这一课的分量,不比任何一个公式轻。
    第五步 · 收网

    9. 花样是被「养」出来的:耗散结构

    0:00 / 0:00

    教师式重点讲解;切换部分时会自动暂停。

    现在回答那个从第 1 段就悬着的问题:系统自发地从无序变得有序,熵不是减少了吗?这不违反第二定律?

    把熵账拆成两笔

    对任何一个和外界有交换的系统,熵的变化都可以拆成两项:

    $$\mathrm{d}S=\mathrm{d}_eS+\mathrm{d}_iS,\qquad \mathrm{d}_iS\geqslant0$$
    项是什么符号能不能为负
    $\mathrm{d}_eS$与外界交换带进带出的熵可以为负——往外排熵完全合法
    $\mathrm{d}_iS$系统内部不可逆过程产生的熵永不为负。这才是第二定律
    第二定律管的是 $\mathrm{d}_iS\geqslant0$,不是 $\mathrm{d}S\geqslant0$。

    「熵只增不减」那句话,成立的条件是孤立系统($\mathrm{d}_eS=0$)。贝纳德那锅液体不是孤立系统——它下面接着热源、上面接着冷源,两头都在和外界交换。

    稳态的账:进来的熵少,出去的熵多

    花样形成之后是一个稳定状态:宏观上不再变化,所以 $\mathrm{d}S/\mathrm{d}t=0$。于是

    $$\frac{\mathrm{d}_eS}{\mathrm{d}t}=-\frac{\mathrm{d}_iS}{\mathrm{d}t}\ <\ 0$$

    系统必须持续地向外净排出熵,排出的速率恰好等于内部产生的速率。这笔账可以算得很具体。设底部热源温度 $T_h$、顶部冷源温度 $T_c$,稳态下从底部进来多少热量,就得从顶部出去多少(记为 $\dot Q$):

    $$\frac{\mathrm{d}_eS}{\mathrm{d}t}=\underbrace{\frac{\dot Q}{T_h}}_{\text{随热流进来的熵}}-\underbrace{\frac{\dot Q}{T_c}}_{\text{随热流出去的熵}}\ <\ 0\quad(\because T_h>T_c)$$

    因此

    $$\frac{\mathrm{d}_iS}{\mathrm{d}t}=\dot Q\left(\frac{1}{T_c}-\frac{1}{T_h}\right)>0$$

    这个式子你在第 2 份见过——它就是两个热源之间热传导的熵产生率。绕了一整章,回到了同一个表达式。

    关键在于:能量是「等量进出」的,熵却不是。

    同样一份热量 $\dot Q$,在高温处携带的熵 $\dot Q/T_h$ 少,在低温处携带的熵 $\dot Q/T_c$ 多。

    所以系统做的事是:收进「高质量」的能量,吐出「低质量」的能量,差额那部分熵被排到了外界。它靠这个差额,替自己内部维持了一个低熵的、有结构的状态。

    代个数:$\dot Q=1\ \mathrm{W}$、$T_h=310\ \mathrm{K}$、$T_c=300\ \mathrm{K}$,则每秒排出 $1.1\times10^{-4}\ \mathrm{J/K}$ 的熵——一小时约 0.39 J/K。花样就是拿这个换来的。

    这类结构有个专门的名字

    靠持续耗散能量来维持的有序结构,称为耗散结构(dissipative structure)。这个概念由普里高津(I. Prigogine)在 1960 年代末提出,他因此获得 1977 年诺贝尔化学奖。

    它和你熟悉的那种有序,是两回事:

    平衡态的有序(晶体、雪花)非平衡的有序(对流胞、化学振荡)
    靠什么维持什么都不靠。孤立起来照样存在靠持续的能量/物质流。断掉就消失
    熵产生率零(已经在平衡态了)大于零,而且持续
    怎么形成降温,让系统趋近平衡加大驱动,让系统远离平衡
    尺度由什么定分子间相互作用(原子尺度)宏观参量(液层厚度 $d$)
    撤掉外界作用纹丝不动几分钟内消失

    最后一行是最好的判别法,也最容易在家里验证:关掉炉子,花样就没了。而冰箱里的冰晶不会因为你关掉冰箱就立刻消失——它只是慢慢化掉,那是另一回事(趋向新的平衡)。

    易错最小熵产生定理有它的适用范围,本节正好在范围之外
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    起点 · 主干刚算出的稳态熵产生率
    $$\frac{\mathrm{d}_iS}{\mathrm{d}t}=\dot Q\left(\frac{1}{T_c}-\frac{1}{T_h}\right)$$

    近平衡区确实有这么一条定理(普里高津的最小熵产生定理):只要①流与力线性、②昂萨格关系成立、③唯象系数可视为常数,那么在给定的外界约束下,稳态就是熵产生率取极小的那个状态,而且这个稳态唯一、稳定。

    这条定理其实解释了第 1 段那句「近平衡不可能从无序到有序」:稳态唯一,就没有第二个状态可选;力变一点,稳态跟着变一点,永远不会突然跳到另一副样子。

    但把它用到本节就错了。看上面这个式子:$T_h$、$T_c$ 固定时,$\mathrm{d}_iS/\mathrm{d}t$ 正比于 $\dot Q$。而对流一旦发生,同样的温差下热流 $\dot Q$ 反而变大了(对流搬热比导热有效,这正是它出现的后果)。

    也就是说:系统选的是熵产生率更大的那个状态,不是更小的。

    这不构成矛盾——定理的三个前提在这里全都不成立($R>R_c$ 时早就不是线性区了)。它只是不适用。

    请连同这条一起记住:近平衡热力学里的漂亮定理,一条都不能不加检查地搬到远离平衡的地方去。这一章前四份和第五份的分界线,就划在这里。

    旁通能推得很远,但有一条必须守住的界线
    0:00 / 0:00

    补充内容独立讲解;收起面板时会自动暂停。

    先说能推的部分。第 1 段提到「生物由各个化学元素组成蛋白质,再组成有机体」,这句话在耗散结构的语言里可以说得很确切:

    • 生命是典型的耗散结构。一个人每天摄入约 $10^{7}\ \mathrm{J}$ 的食物能量,几乎等量地以低温热的形式排给环境。能量收支大致平衡,熵收支不平衡——正是靠净排熵,才维持住体内那套精密的秩序。停止进食与散热,结构在几天到几周内瓦解。
    • 大气与洋流:吸收高温的太阳辐射($\sim6000\ \mathrm{K}$),向宇宙辐射低温红外($\sim255\ \mathrm{K}$),中间那份熵差养出了环流、气旋、洋流——全是巨型贝纳德胞。
    • 台风是最教科书式的一个:暖海面供能,高空排热,结构(眼、眼壁、螺旋雨带)尺度确定、相当稳定,一旦登陆离开暖水面就迅速瓦解。
    再说必须守住的界线。「耗散结构」是一个必要条件的框架,不是万能解释:

    · 它能告诉你为什么有序不违反第二定律、什么条件下可能出现结构;
    · 它不能告诉你会出现哪一种结构,更不能推出「生命必然产生」「社会必然如何演化」。

    上世纪八九十年代,这套词汇被大量借到经济学、社会学乃至管理学里,多数属于类比而非推导。作为比喻可以,作为论证不行——判据很简单:你能写出那个系统的 $R$ 吗?能算出它的 $R_c$ 吗?写不出来,就只是借了个词。
    收尾

    10. 这一节,和这一整章

    0:00 / 0:00

    教师式重点讲解;切换部分时会自动暂停。

    把 8.6 的链条串一遍

    步骤结论
    前提被打破温度梯度增大 → 线性唯象关系失效(失效的是「整体流正比于整体力」,不是局部平衡)
    现象贝纳德对流:过阈值后突然出现规则花样,尺寸 ≈ $d$,中心上升、边缘下沉
    花样取决于边界刚性-自由 → 六角;刚性-刚性 → 卷(方容器直卷、圆容器环卷)
    机制(一)自由面:$\sigma$ 随 $T$ 减小 → 热点被拉开 → 下方补流。太空船实验($10^{-6}g$)证明与重力无关
    机制(二)两面刚性:浮力 vs 黏滞 + 热扩散
    判据$R=\dfrac{g\alpha\Delta T}{\eta D_T}d^{3}$;$R<R_c$ 导热,$R>R_c$ 对流
    临界值657.511 / 1100.650 / 1707.762;实验 $1700\pm51$ 与刚性-刚性一档相符。$657.511=\tfrac{27}{4}\pi^{4}$,可手算
    另一个舞台BZ 反应:浓度的时间振荡,四个必要条件与对流完全同构
    统一解释耗散结构:$\mathrm{d}S=\mathrm{d}_eS+\mathrm{d}_iS$,稳态靠 $\mathrm{d}_eS<0$ 持续排熵养着,断供即消失

    这一整章走过的路

    五份读完,第 8 章的结构其实非常清晰——它是一部「把不可逆过程装进热力学」的连续剧:

    份做的事留下的东西
    第 1 份
    8.1–8.2
    四条经验定律(傅里叶、欧姆、牛顿、菲克)长得一模一样 → 抽象成「力 × 流」的框架;提出局部平衡假设、线性唯象方程、昂萨格倒易关系一套语言
    第 2 份
    8.3
    从熵的平衡方程算出熵产生率 $\dfrac{\mathrm{d}_iS}{\mathrm{d}t}=\sum_j J_jX_j$,把「力和流该怎么配对」这件事定死一把尺
    第 3 份
    8.4
    电动效应:四个实验效应被昂萨格关系压成两个独立系数,得到萨克斯关系第一次兑现
    第 4 份
    8.5
    热电效应:三个系数全部由一个 $\varepsilon$ 写出,开尔文两条关系式从猜想变成定理第二次兑现,更漂亮
    第 5 份
    8.6
    把「近平衡」这个前提拆掉:远离平衡时,系统自己长出结构一扇门
    整章的分水岭,就在最后这一份。

    第 1–4 份属于线性非平衡态热力学:有完整的框架、有普遍的定理(昂萨格关系)、有可以一路推到底的演绎链。它成熟、封闭、漂亮。

    第 5 份进入的是非线性区:没有那样的普遍定理,每一类系统要单独分析,答案往往依赖数值计算。它不封闭,但它才是真实世界的常态——大气、海洋、地幔、恒星、以及一切活着的东西,全都在这一边。

    本份小结

    这一份的逻辑链条是:把第 1–4 份的第②条前提(线性唯象关系)拆掉 → 贝纳德实验事实:过阈值突然出现规则花样,六角尺寸 ≈ $d$,中心上升沿边下沉,一经形成相当稳定,是一个非平衡态相变 → 花样由界面决定:刚性-自由出六角、刚性-刚性出卷 → 两种驱动机制:自由面靠表面张力($\mathrm{d}\sigma/\mathrm{d}T<0$,$10^{-6}g$ 的太空实验排除了重力),两面刚性靠浮力(对手是黏滞与热扩散两个刹车)→ 把「浮力克服黏滞力」写成无量纲判据 $R=\frac{g\alpha\Delta T}{\eta D_T}d^{3}$($\eta$ 须取运动黏度),$R<R_c$ 导热、$R>R_c$ 对流 → 三个临界值 657.511 / 1100.650 / 1707.762 由中性曲线 $R(a)=\frac{(\pi^2+a^2)^3}{a^2}$ 的极小给出,实验 $1700\pm51$ 与刚性-刚性相符,且由 $a_c$ 反推出花样尺度只由 $d$ 定(卷宽 $\approx d$,六角直径 $2\sim3\,d$) → BZ 化学振荡给出同一张条件清单(远离平衡、开放、机制复杂、非线性)→ 统一解释:耗散结构靠 $\mathrm{d}_eS=-\mathrm{d}_iS<0$ 持续排熵维持,断掉供能即消失。

    第 8 章到此结束。从四条互不相干的经验定律出发,一路走到熵产生率、昂萨格关系、电动与热电效应,最后走出线性区,看到系统自己长出结构——这一章真正教会的,是怎么把「不可逆」这件事请进热力学的框架里,以及这个框架在哪里为止。