教师式重点讲解;切换部分时会自动暂停。
先看一件事,它不需要任何前置知识:
请注意这里最反直觉的一点:我们只是把「推得更狠」,系统却变得更「有条理」了。加热是一件彻底不带任何空间信息的事——底部处处一样热,容器是圆的或方的,没有谁规定六角形该长在哪里。可是花样出现了,而且尺寸、形状都可以算出来。
这一节讲的就是这类现象。它在这一章里位置特殊:它是前四份成立的那个前提被打破之后的世界。
回头数一下,第 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$ | 换个舞台:化学振荡 | 统一解释:耗散结构与熵账 |
补充内容独立讲解;收起面板时会自动暂停。
力与流。把「不平衡的程度」写成力 $X_i$(温度不均、浓度不均、电势不均……),把「因此发生的输运」写成流 $J_j$(热流、粒子流、电流……)。
线性唯象方程。只要偏离平衡不远,每一个流都正比于所有的力:
对角项 $L_{jj}$ 就是老定律(傅里叶、欧姆、菲克……);非对角项 $L_{ji}$ 是交叉效应——第 3、4 份整两份都在讲它。
昂萨格倒易关系。$L_{ij}=L_{ji}$。它是白送的一半信息:只要测出一个方向的交叉效应,反方向的不用测也知道。
补充内容独立讲解;收起面板时会自动暂停。
不是。在贝纳德对流里,①仍然成立得好好的。
液体里每一小块都还有确定的温度、压强、密度;傅里叶定律在每一点上依然成立(热流仍然正比于当地的温度梯度),黏滞定律也一样。液体离「分子层面的平衡」远得很吗?并不——常温常压下的水,分子碰撞把局部平衡维持得非常牢。
那到底是什么变成非线性的?是方程本身,不是本构关系。流体的运动方程里有一项
叫对流项:液体一动,就把自己和自己的热量搬着走。速度乘速度、速度乘温度——这是两个未知量相乘,货真价实的非线性。液体不动($\mathbf{v}=0$)时它是零,所以近平衡时看不见它;一旦动起来,它就是主角。
教师式重点讲解;切换部分时会自动暂停。
把第 1 段那个装置说细一点。
容器里装一薄层液体,厚度为 $d$,从容器底部均匀加热。记底部温度 $T_0$、液体表面温度 $T_1$,且 $T_1<T_0$。
当 $T_0-T_1$ 之差较小时,稳恒状态下的热流满足傅里叶定律——热量靠传导一路向上,液体宏观上静止。
但请留意随之而来的一件事:由于存在热膨胀,与温度梯度相伴随,液体里有一个线性的密度梯度。下面热、所以轻;上面冷、所以重。
1900 年,法国物理学家贝纳德(Bénard)用熔化的鲸蜡作为液体,以 $100\ ^\circ\mathrm{C}$ 的水蒸气在底部加热,发现:
这个现象就称为贝纳德对流。
三件事值得单独记住:
| 观察到的 | 为什么值得惊讶 |
|---|---|
| 有阈值:温差小时什么都没有,过了某个值才发生 | 不是「温差越大流得越多」这种连续的响应——是开关式的 |
| 有尺度:胞的大小 ≈ $d$,不是随机的 | 加热是均匀的,容器也没规定尺寸,尺度是系统自己选的 |
| 有对称性破缺:原本水平方向处处等价,现在分出了「中心」和「边」 | 方程和边界条件都是水平各向同性的,解却不是 |
第三条是最深的一条。加热、容器、液体本身,在水平方向上都没有任何偏好;可是结果挑出了一套六角格子,把「处处一样」变成了「有中心、有边」。这在物理里有个专门的名字:自发对称性破缺——和铁磁体在居里点以下自己挑一个磁化方向是同一件事。
补充内容独立讲解;收起面板时会自动暂停。
补充内容独立讲解;收起面板时会自动暂停。
如果液体是理想流体(没有黏性、没有热传导),那么只要下面比上面轻,哪怕差一丁点,静止状态就立刻失稳——这叫瑞利-泰勒不稳定性,没有阈值。
真实液体不是这样,因为有两件事在拖后腿:
所以真正的问题不是「上下重不重」,而是:一个微小的扰动,是会长大还是会衰减?
教师式重点讲解;切换部分时会自动暂停。
贝纳德那锅鲸蜡有一个容易被忽略的细节:它的上表面是敞着的。
如果液体表面暴露在空气中,称此界面为自由界面;底部和容器接触,称此界面为刚性界面。两者的差别是实打实的力学差别:
| 刚性界面 | 自由界面 | |
|---|---|---|
| 液体在界面上能不能滑动 | 不能。紧贴固体的那层液体速度为零(无滑移) | 能。液面可以整体被拉着走 |
| 切向的力 | 固体可以给液体任意大的切向阻力 | 只有空气,几乎给不了阻力;但液面自己的表面张力可以拉它 |
| 后果 | 更「刹得住」,需要更大的温差才动 | 更容易动,而且多了一种全新的驱动力 |
贝纳德的原始实验是「刚性底 + 自由面」。他看到的六角花样,用硅油拍摄的照片和理论计算的一个六角胞,就是课本图 8.10 的 (a) 和 (b)。
后来又做了很多对流实验,绝大部分用上下面均与刚性界面接触,无自由表面(也就是上面盖一块板)。结果是:
课本给了一句很硬的判断:这些图案一经形成就相当稳定,是一个相变。
「相变」这个词用在这里不是打比方。把它和你熟悉的相变逐条对照:
| 相变的特征 | 在贝纳德对流里 |
|---|---|
| 有一个临界值,越过它性质突变 | 有。临界温差 / 临界瑞利数 $R_c$(第 6、7 段) |
| 对称性破缺:新相的对称性低于旧相 | 有。静止态水平方向连续平移对称;对流态只剩离散的格子对称 |
| 序参量随驱动长大 | 有。序参量是对流速度的振幅 $A$。上下边界对称的卷:$A\propto\sqrt{R/R_c-1}$,从零连续长大(像二级相变);六角:阈值处跳变、带滞后(像一级相变) |
| 新相稳定,扰动它会自己回来 | 有。这正是「一经形成就相当稳定」 |
补充内容独立讲解;收起面板时会自动暂停。
设对流的强弱可以用一个数 $A$(振幅)刻画。它随时间怎么变?先写最粗的一项:
$R<R_c$ 时 $\varepsilon<0$,扰动指数衰减;$R>R_c$ 时 $\varepsilon>0$,扰动指数长大——这就是失稳。
但指数长大不可能一直持续下去,对流强到一定程度会自己限制自己(流动越强,它搬运热量越有效,把驱动它的温度梯度削平了)。补上最低阶的自限制项:
为什么是 $A^3$ 而不是 $A^2$:如果把整个花样反过来(上升变下降)也是一个同样合法的解,那么方程必须在 $A\to-A$ 下不变,偶次项被对称性禁掉。
令 $\mathrm{d}A/\mathrm{d}t=0$:
但那个「如果」是有条件的,条件不满足时结论要改。上下边界同种、且液体的物性可以当作不随温度变($A\to-A$ 才真是一个对称操作)时,得到的是卷,上面那套成立。
而六角不满足:中心上升和中心下沉不等价——这正是六角能分出「中心」和「边」的原因。对称性一破,$A^2$ 项就被允许写进方程:
补充内容独立讲解;收起面板时会自动暂停。
| 图 | 界面 | 是什么 | 看到什么 |
|---|---|---|---|
| 8.10(a) | 刚性-自由 | 照片(硅油) | 一大片六角格子 |
| 8.10(b) | 刚性-自由 | 理论计算 | 单个六角胞里的流线 |
| 8.11(a) | 刚性-刚性 | 照片(侧面) | 方容器里的直卷 |
| 8.11(b) | 刚性-刚性 | 照片(俯视) | 圆容器里的环形卷 |
| 8.12 | 自由-自由 | 理论计算 | 一排平行的卷 |
教师式重点讲解;切换部分时会自动暂停。
理论上做了大量的工作,来理解对流产生的原因。下面给出定性上的解释——两种界面组合,对应两套完全不同的驱动力。先讲贝纳德本人那一种。
对于底部是刚性界面、而上面是自由界面的情况,形成对流的驱动力是表面张力。
链条一步一步是这样的:
| 第几步 | 发生了什么 |
|---|---|
| 1 | 由于底部温度高,分子运动的速度较大,同时存在涨落 |
| 2 | 某个速度大的液体元偶然到达上层,引起自由表面上一小块面积元的局部加热 |
| 3 | 而液体的表面张力随温度的升高而减小,所以自由面上这块被加热的小面积元,表面张力显著变小 |
| 4 | 于是周围(较冷、张力较大)的表面元把液体从热处向外拉 |
| 5 | 此处液面被拉走、压强减小,于是引起液体从下面向上升来补充 |
| 6 | 上来的液体又是热的 → 回到第 3 步。正反馈闭合,对流就开始了。 |
课本给了一个很硬的证据。当然这里还有重力的影响(液体在地面上,浮力一直都在),要判断到底谁说了算,最干脆的办法是把重力去掉:
这是一个教科书式的控制变量实验:把嫌疑人之一(重力)拿掉,现象照旧,于是它被排除。
严格一点说,这个实验证明的是「对这样的薄液层,表面张力足以独立驱动对流,不需要重力帮忙」。它并不意味着浮力在哪里都不重要——液层一厚,主角就换人了(第 5 段讲的正是那一支)。
补充内容独立讲解;收起面板时会自动暂停。
两种等价的读法。
本段用的是力的读法:热点周围的液面「拉得更用力」,所以把热点处的液体拽走。
为什么随温度升高而减小?表面张力的来源是「表面上的分子比内部少了一半邻居,能量偏高」。温度升高时,分子热运动加剧、液体密度下降、分子间平均距离拉大,表面分子与内部分子的差别被抹平,这份额外能量就变小。
补充内容独立讲解;收起面板时会自动暂停。
表面张力驱动这一支,判据的写法和第 6 段的瑞利数完全同构——推动量除以两个刹车量:
称马兰戈尼数,临界值约 $Ma_c\approx80$。它和 $R$ 的差别只在两处:推动量换成了表面张力随温度的变化率,而厚度只出现一次方(因为驱动力只作用在表面这一层),不是 $d^3$。
两支的强弱之比,把 $\Delta T$、$\nu$、$D_T$ 全约掉了:
历史上这是一桩公案:贝纳德 1900 年的实验,长期被当作浮力对流的典范;直到 1956 年布洛克(Block)和 1958 年皮尔逊(Pearson)才指出,那层薄薄的液体里真正说了算的是表面张力。所以今天的术语把两者分开:瑞利-贝纳德对流指浮力那一支,贝纳德-马兰戈尼对流指表面张力这一支。
教师式重点讲解;切换部分时会自动暂停。
换一种边界条件:对于上下两面均为刚性界面,对流是由重力引起的。上面盖了板,没有自由表面,表面张力这条路被堵死了,剩下的只有浮力。
课本的链条同样是三步:
| 第几步 | 发生了什么 |
|---|---|
| 1 | 当容器底部加热时,较下面的液体层受热膨胀,产生一个浮力,企图向上升 |
| 2 | 但液体的黏滞力将阻止它上升 |
| 3 | 当底部的温度进一步升高时,浮力将克服黏滞力,而上面较冷的液体将向下运动,从而形成宏观的对流 |
「浮力克服黏滞力」这句话是全节的转折点:它把一个定性的故事,变成了一个可以比大小的不等式。既然是比大小,就一定能写成一个比值——这个比值就是下一段的主角。
课本第 2 步只提了黏滞力。但真正让阈值出现的,是两件事一起在拖后腿:
| 它做的事 | 它拖慢的是 | |
|---|---|---|
| 浮力 | 热液元比周围轻,被往上推 | ——(这是推动方) |
| 黏滞 | 周围液体磨它,动能变成热耗散掉 | 它动的能力 |
| 热扩散 | 热量从液元漏向周围,温差被抹平 | 它轻的能力(浮力的来源) |
补充内容独立讲解;收起面板时会自动暂停。
「热的轻所以往上浮」这句话,作为对流的方向说明是对的,但它解释不了本节要解释的两件事:
把问法换掉,两个问题一起解决:不问「静止状态合不合理」,而问——
这个问法自动带来两样东西:
① 长大和衰减之间有一条分界线,那就是阈值;
② 扰动可以有各种波长,不同波长长大的难易不同,最容易长大的那一个波长胜出——这就是尺寸的来源。
这套问法叫线性稳定性分析:它不解非线性方程,只问「零附近的小扰动往哪个方向走」。第 7 段那三个精确到小数点后三位的临界值,就是这么算出来的。
教师式重点讲解;切换部分时会自动暂停。
第 5 段结束在一句话上:浮力要克服黏滞力。要把它变成能算的东西,就得把「推动」和「阻碍」各自量化,然后相除。对此瑞利建立了一个模型,出现对流的条件用瑞利数 $R$ 来衡量,定义为
逐项读一遍。课本把它分成两组:
| 组 | 符号 | 是什么 | 为什么在这一边 |
|---|---|---|---|
| 使对流发生 (分子) | $g$ | 重力加速度 | 没有重力就没有浮力(第 4 段的太空船实验) |
| $\alpha$ | 热膨胀系数 | 热了要能变轻才有浮力;$\alpha=0$ 的液体绝不对流 | |
| $\Delta T$ | 上下底之间的温度差 | 温差越大,上下密度差越大 | |
| 阻止对流发生 (分母) | $\eta$ | 黏滞系数 | 磨掉液元的动能 |
| $D_T$ | 热扩散系数 | 抹平液元与周围的温差,把浮力本身抹掉 | |
| 几何 | $d^{3}$ | 液体的厚度 | 见下面 |
分子那三个正是第 5 段的「推动」,分母那两个正是那两个刹车。$R$ 就是它们的比值:$R$ 大 = 推得过;$R$ 小 = 按得住。
$R$ 没有单位。这意味着比较不同的系统时,只需要比这一个数:一锅毫米厚的油、一杯水、几公里厚的大气、几千公里厚的地幔,尺度跨了九个数量级,只要 $R$ 相同,行为就相同。这也是为什么实验室里几毫米厚的一层硅油,能拿来研究地幔里的对流。
于是判据可以写得非常干脆。当 $R<R_c$(临界瑞利数)时,液体中发生的是线性的热传导过程;当 $R>R_c$ 时,出现对流,并按不同的边界条件出现不同的规则图案。
| $R<R_c$ | $R>R_c$ | |
|---|---|---|
| 液体的状态 | 宏观静止 | 规则对流 |
| 热怎么走 | 纯导热,傅里叶定律 | 导热 + 对流搬运 |
| 属于本章哪一半 | 第 1–4 份的线性区 | 本节的非线性区 |
厚度进入的是三次方,这是 (8.6.1) 里最有威力的一项:层厚翻一倍,同样的温差算出的 $R$ 变成 8 倍;反过来说,临界温差降到原来的 $1/8$。
所以一杯稍微深一点的水,一丁点温差就已经在悄悄翻腾;而上下封死的薄液层想靠浮力对流,难得很——1 mm 厚的水要差 120 K。
补充内容独立讲解;收起面板时会自动暂停。
分子的量纲:
要让整个 $R$ 无量纲,分母必须也是 $\mathrm{m^4/s^2}$。热扩散系数是死的:$[D_T]=\mathrm{m^2/s}$。所以
$\mathrm{m^2/s}$ 正是运动黏度 $\nu$ 的单位。
若误用动力黏度($\mathrm{Pa\cdot s}=\mathrm{kg/(m\cdot s)}$):
补充内容独立讲解;收起面板时会自动暂停。
水在 $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$」:
| 层厚 $d$ | $R/\Delta T$(每 K) | 临界温差 $\Delta T_c$ |
|---|---|---|
| 1 mm | 14.2 | 120 K(实际上做不到,水早沸了) |
| 2 mm | 113 | 15 K |
| 5 mm | 1773 | 0.96 K |
| 1 cm | 14186 | 0.12 K |
补充内容独立讲解;收起面板时会自动暂停。
「浮力克服黏滞力」——这句话没法用。浮力和黏滞力的量纲相同,但它们各自依赖液体厚度、温差、材料,把它们分别算出来再比大小,换一种液体、换一个容器就得重算一遍,而且永远得不到一个能写进教科书的判据。
更糟的是:实验上换了液体、换了厚度,临界温差从零点几度跳到上百度,看起来毫无规律。没有 $R$,这些实验就是一堆互相对不上的数据。有了它,全部塌缩成一句话:过 1708 就对流。
本份第 7 段会算出它的三个临界值并和实验对照,第 9 段用它来说明「花样必须靠持续供能养着」。走出本书,它是气象、海洋、地球物理和恒星内部对流的标准判据;工程上换热器、晶体生长、建筑通风的设计里,它是必查的第一个数。
教师式重点讲解;切换部分时会自动暂停。
$R_c$ 不是量出来的,是算出来的。理论上做了大量的计算工作,主要是分析一组平衡方程:
| 方程 | 管的是 |
|---|---|
| 有局域质量密度平衡方程 | 液体不会凭空增减(不可压缩) |
| 包括速度的动量平衡方程 | 浮力、压强、黏滞力怎么推动液体 |
| 内能的平衡方程 | 热怎么被传导、又怎么被流动带着走 |
根据不同的边界条件,计算出何时发生出现对流图案的相变,以及对流图案的花样。前面的图 8.10(b) 就是计算得到的刚性-自由界面的一个六角的对流花样,而图 8.12 是对自由-自由界面计算得到的对流卷图案。
理论上计算的瑞利数为
实验上得到的瑞利数 $R_c=1700\pm51$,可见符合相当好。
三个数的大小顺序本身也有物理含义:$657<1101<1708$,刚性面越多,需要的推动越大。因为刚性面上液体不能滑动(速度必须为零),液体要在贴壁处从零加速起来,速度梯度大,黏滞耗散就大——每多一个刚性面,就多一份刹车。
三档里,自由-自由那一档的边界条件最干净,能得到解析解。理论给出的中性稳定曲线(即「扰动既不长大也不衰减」所对应的 $R$)是
其中 $a$ 是水平方向扰动的无量纲波数(以液层厚度 $d$ 为长度单位)。这条曲线是一条 U 形曲线:
为什么是最低点?因为把 $R$ 从小往大调,第一个碰到曲线的地方就是最低点——那个波数的扰动最先长大,它挑走了整个花样的尺寸。所以:
令 $x=a^{2}$,求极小:
即 $a_c=\pi/\sqrt2\approx2.221$。代回去:
三档的临界波数分别是 $a_c=2.221$、$2.682$、$3.117$(都是浮力这一支的值)。波数对应的水平波长(以 $d$ 为单位)是 $\lambda=2\pi/a$,于是 $\lambda=2.83\,d$、$2.34\,d$、$2.02\,d$——都是两三倍的层厚。
从波长换算到「一个胞有多大」,卷和六角要分开算:
| 界面 | $\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$ |
还有一处要交代清楚:上面三个 $a_c$ 属于浮力那一支。贝纳德本人那个实验是表面张力驱动的(第 4 段),它有自己的临界波数 $a_c\approx1.99$,给出的波长 $3.16\,d$、六角直径 $3.65\,d$——比浮力那一支更大,也和实验上看到的、明显宽于层厚的胞径相符。
但无论按哪种读法、哪一支机制,结论不变:花样的尺度只由一个长度决定,就是液层厚度 $d$。
补充内容独立讲解;收起面板时会自动暂停。
完整推导要用到流体力学,这里只交代它是怎么被逼出来的,三步:
第一步:把解写成「静止解 + 小扰动」。静止的纯导热态是方程组的一个精确解(速度处处为零、温度线性分布)。设速度、温度各偏离它一点点,代入方程,把所有二次以上的小量丢掉——包括那个惹事的对流项 $(\mathbf v\cdot\nabla)$,因为它是两个小量相乘。剩下的是一组线性方程。
第二步:按波数分解。线性方程的好处是不同波长的扰动互不干扰,可以一个一个算。于是设扰动形如
水平方向按波数 $a$ 分解,竖直方向要满足上下边界条件(自由-自由时正好是正弦,这就是这一档能手算的原因),$s$ 是增长率。
第三步:令 $s=0$。$s>0$ 扰动长大、$s<0$ 衰减,分界线是 $s=0$。把 $s=0$ 代进去,方程给出 $R$ 与 $a$ 的一个关系——就是那条 $R(a)$。竖直方向取最低的模式 $n=1$(能量最省),得到主干那个式子。
补充内容独立讲解;收起面板时会自动暂停。
$R$ 不是孤例,它有一大家子亲戚,构造方式完全一样:
| 数 | 推动 ÷ 阻碍 | 阈值意味着 |
|---|---|---|
| 雷诺数 $Re=\dfrac{vL}{\nu}$ | 惯性 ÷ 黏滞 | 层流 → 湍流 |
| 瑞利数 $R$ | 浮力 ÷(黏滞 × 热扩散) | 导热 → 对流 |
| 马兰戈尼数 $Ma$ | 表面张力梯度 ÷(黏滞 × 热扩散) | 静止 → 表面驱动对流 |
| 马赫数 $M=\dfrac{v}{c}$ | 流速 ÷ 声速 | 亚声速 → 超声速 |
为什么非要凑成无量纲?因为物理规律不该依赖你用米还是用厘米。一个只含无量纲组合的判据,换单位不变、换尺度也不变。
教师式重点讲解;切换部分时会自动暂停。
到这里为止,所有例子都是流体。如果非平衡态相变只发生在液体里,它就只是流体力学的一个专题。另一个例子是化学反应,它把这件事的普遍性一下子推开了。
请对照第 3 段的清单看这三条——它们和「对流花样一经形成就相当稳定」是同一件事,只不过贝纳德的结构长在空间里,化学振荡的结构长在时间里。
研究最多的是 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+}}$ 的比例走),周期从几十秒到几分钟,肉眼可见,可以持续几十个周期。课本描述的「红蓝之间振荡」是演示时最常见的样子——不过那抹鲜艳的红与蓝其实来自另加的指示剂,纯铈体系本身只是「无色 ↔ 淡黄」。
课本把条件列得很明确:产生化学振荡的条件很苛刻,要满足:远离平衡,对物质源及能量源要开放,反应机制必须复杂(没有一个简单系统有振荡),速率方程必须包含非线性项。
| 条件 | 为什么非要不可 |
|---|---|
| 远离平衡 | 平衡态附近一切都单调地趋向平衡;细致平衡原理直接禁止了持续的循环 |
| 对物质源、能量源开放 | 振荡要靠消耗丙二酸维持。封住系统,燃料烧完,振荡必然停下来——BZ 在烧杯里也只能振几十个周期 |
| 机制复杂 | 要有自催化(产物促进自己生成)和反馈抑制两个环节,缺一不可 |
| 速率方程含非线性项 | 这是数学上的硬要求:只有非线性系统才可能有「振幅由自己决定、被扰动后能自己回来」的持续振荡 |
把这四条和贝纳德那一套对照,会发现是同一张清单:远离平衡($R>R_c$)、持续供能(底部一直加热)、有正反馈(第 4 段那条闭合链)、方程非线性(对流项)。两个毫不相干的系统,长出结构所需要的条件完全一样。这正是下一段要收的网。
补充内容独立讲解;收起面板时会自动暂停。
纯铈体系的颜色对比很弱:$\mathrm{Ce^{3+}}$ 无色,$\mathrm{Ce^{4+}}$ 淡黄。肉眼能看出振荡,但绝谈不上「红蓝分明」。
课堂演示和网上的视频里那种鲜艳的红 ↔ 蓝,来自额外加入的氧化还原指示剂「邻二氮菲亚铁」(ferroin):还原态是红色,氧化态是蓝色。它跟着铈的价态一起翻转,把一个微弱的信号放大成了强烈的视觉效果。
补充内容独立讲解;收起面板时会自动暂停。
先承认一半:线性系统当然可以振荡——弹簧振子就是。但线性系统的振荡只有三种命运:
而课本要求的振荡是:与初始条件无关、对小扰动稳定。三种命运没有一种做得到。
BZ 反应里的非线性项在哪?核心是溴的自催化步骤:亚溴酸 $\mathrm{HBrO_2}$ 参与反应,结果生成两份 $\mathrm{HBrO_2}$——
它对 $[\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^-}$ → 抑制剂回来,火被扑灭 → 从头再来。
补充内容独立讲解;收起面板时会自动暂停。
1951 年,苏联生物化学家别洛乌索夫(Б. П. Белоусов)在研究一个模拟三羧酸循环的体系时,看到溶液颜色周期性地来回变化。他投出论文,被拒,理由是「明显违反热力学第二定律」;改投另一家,再次被拒。他最终只在一份不起眼的会议摘要集上发表了它,此后离开了这个课题。
直到 1961 年前后,扎鲍廷斯基重做并系统研究了这个反应,它才被接受。1972 年 Field、Kőrös、Noyes 给出了详细的反应机理(FKN 机理),把那个自催化步骤揪了出来。
教师式重点讲解;切换部分时会自动暂停。
现在回答那个从第 1 段就悬着的问题:系统自发地从无序变得有序,熵不是减少了吗?这不违反第二定律?
对任何一个和外界有交换的系统,熵的变化都可以拆成两项:
| 项 | 是什么 | 符号能不能为负 |
|---|---|---|
| $\mathrm{d}_eS$ | 与外界交换带进带出的熵 | 可以为负——往外排熵完全合法 |
| $\mathrm{d}_iS$ | 系统内部不可逆过程产生的熵 | 永不为负。这才是第二定律 |
花样形成之后是一个稳定状态:宏观上不再变化,所以 $\mathrm{d}S/\mathrm{d}t=0$。于是
系统必须持续地向外净排出熵,排出的速率恰好等于内部产生的速率。这笔账可以算得很具体。设底部热源温度 $T_h$、顶部冷源温度 $T_c$,稳态下从底部进来多少热量,就得从顶部出去多少(记为 $\dot Q$):
因此
这个式子你在第 2 份见过——它就是两个热源之间热传导的熵产生率。绕了一整章,回到了同一个表达式。
靠持续耗散能量来维持的有序结构,称为耗散结构(dissipative structure)。这个概念由普里高津(I. Prigogine)在 1960 年代末提出,他因此获得 1977 年诺贝尔化学奖。
它和你熟悉的那种有序,是两回事:
| 平衡态的有序(晶体、雪花) | 非平衡的有序(对流胞、化学振荡) | |
|---|---|---|
| 靠什么维持 | 什么都不靠。孤立起来照样存在 | 靠持续的能量/物质流。断掉就消失 |
| 熵产生率 | 零(已经在平衡态了) | 大于零,而且持续 |
| 怎么形成 | 降温,让系统趋近平衡 | 加大驱动,让系统远离平衡 |
| 尺度由什么定 | 分子间相互作用(原子尺度) | 宏观参量(液层厚度 $d$) |
| 撤掉外界作用 | 纹丝不动 | 几分钟内消失 |
最后一行是最好的判别法,也最容易在家里验证:关掉炉子,花样就没了。而冰箱里的冰晶不会因为你关掉冰箱就立刻消失——它只是慢慢化掉,那是另一回事(趋向新的平衡)。
补充内容独立讲解;收起面板时会自动暂停。
近平衡区确实有这么一条定理(普里高津的最小熵产生定理):只要①流与力线性、②昂萨格关系成立、③唯象系数可视为常数,那么在给定的外界约束下,稳态就是熵产生率取极小的那个状态,而且这个稳态唯一、稳定。
这条定理其实解释了第 1 段那句「近平衡不可能从无序到有序」:稳态唯一,就没有第二个状态可选;力变一点,稳态跟着变一点,永远不会突然跳到另一副样子。
请连同这条一起记住:近平衡热力学里的漂亮定理,一条都不能不加检查地搬到远离平衡的地方去。这一章前四份和第五份的分界线,就划在这里。
补充内容独立讲解;收起面板时会自动暂停。
先说能推的部分。第 1 段提到「生物由各个化学元素组成蛋白质,再组成有机体」,这句话在耗散结构的语言里可以说得很确切:
教师式重点讲解;切换部分时会自动暂停。
| 步骤 | 结论 |
|---|---|
| 前提被打破 | 温度梯度增大 → 线性唯象关系失效(失效的是「整体流正比于整体力」,不是局部平衡) |
| 现象 | 贝纳德对流:过阈值后突然出现规则花样,尺寸 ≈ $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 份的第②条前提(线性唯象关系)拆掉 → 贝纳德实验事实:过阈值突然出现规则花样,六角尺寸 ≈ $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 章到此结束。从四条互不相干的经验定律出发,一路走到熵产生率、昂萨格关系、电动与热电效应,最后走出线性区,看到系统自己长出结构——这一章真正教会的,是怎么把「不可逆」这件事请进热力学的框架里,以及这个框架在哪里为止。