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

8.1–8.2 输运现象与「力·流」框架

从四条经验定律,到一套能装下它们的统一语言
目录(点击展开,也可直接点条目跳转)
    1 / 11
    第零步 · 先看地图

    1. 前面七章漏掉了一件事:快慢

    0:00 / 0:00

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

    先说一件容易被忽略、但决定了这一整章存在意义的事。

    前面七章的热力学,处理的全都是平衡态,走的过程全都是可逆过程(也就是准静态过程)。它非常擅长回答两类问题:

    • 这件事能不能发生?——比如两个物体接触,热会不会自己从冷的流向热的(不会,熵减了)。
    • 如果能,最后停在哪?——比如化学反应进行到什么程度会停下来(停在 $\sum_i \mu_i\nu_i = 0$ 那一点)。

    但它完全不回答第三类问题:

    这件事发生得多快?一杯热水放在桌上会变凉,热力学能告诉你它一定会凉、最后凉到室温;但它一个字也说不出「要凉多久」。

    「多快」这个问题,热力学不是不想答,而是结构上答不了——它的全部工具($U$、$S$、$G$、各种偏导数关系)里根本就没有出现过时间 $t$。

    根子在一条更基本的限制上,这条限制这一整章都要反复用到,请先记牢:

    $T$、$S$、$p$、$\mu$ 这些热力学量,本来就只对平衡态才有定义。

    一杯左边 80 ℃、右边 20 ℃ 的水,你说它的温度是多少?没有答案——此刻「这杯水的温度」这个说法本身就不成立。

    为了让这些量在整个过程中一直有定义,热力学只好要求过程每一瞬间都是平衡态,也就是走得无限慢(准静态)。而问一个「无限慢」的理论「这要多久」,它当然只能回答「无限久」。

    所以这一章要做的第一件事,是把这个前提松开一点:允许系统偏离平衡,但只偏离一点点,再看看热力学还能站得住多少。

    深入时间是怎么被赶出热力学的
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干里那句「工具里没有时间」
    $$dU = T\,dS - p\,dV + \sum_i \mu_i\,dN_i$$

    看看这个式子里有什么:$U$、$S$、$V$、$N_i$、$T$、$p$、$\mu_i$。一个 $t$ 都没有。这不是作者忘了写,是被前提排除掉的。

    前提是这样一环扣一环的:

    • 热力学量($T$、$p$、$S$)只对平衡态有定义。一杯左边 80 ℃、右边 20 ℃ 的水,你说它的温度是多少?没有答案——「温度」这个词此刻不成立。
    • 所以要让这些量在整个过程中一直有定义,过程就必须每一瞬间都是平衡态——这就是准静态过程。
    • 而「每一瞬间都平衡」意味着要走得无限慢。一旦走得有限快,系统内部就会出现温度差、浓度差,就不是平衡态了。

    于是绕了一圈回到原点:热力学为了让自己的量有定义,把过程的速率设成了「无限慢」。问一个「无限慢」的理论「这要多久」,它当然只能回答「无限久」。

    这一章要做的,正是把这个前提松开一点点——允许系统偏离平衡,但只偏离一点点,然后看看热力学还能不能站得住、能多说出些什么。

    松开多少?——近平衡与远平衡

    把前提松开,可以松一点,也可以松很多,效果完全不同。这一章把它分成两档:

    档位什么意思会发生什么本章篇幅
    近平衡 偏离平衡不远,比如温度梯度比较小 「原因」和「结果」成正比(线性关系)。温度差大一倍,热流就大一倍 8.1–8.5,绝大部分
    远平衡 偏离平衡很远,比如温度梯度很大 正比关系失效,成非线性关系。系统反而会自发长出规则的结构 8.6,只做简介

    这一章的名字叫非平衡热力学,也叫不可逆过程热力学——两个名字指同一件事,看到哪个都别慌。

    回顾三个最容易混的前置词
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    平衡态:系统放着不管,各处的温度、压强、成分都不再变化了,也没有任何宏观的流动。这是一个状态。

    准静态过程:过程走得足够慢,慢到系统在每一瞬间都来得及回到平衡态。这是一个过程。它是一串平衡态连成的线。

    可逆过程:过程能够原路返回,并且让系统和外界都回到原样,不留任何痕迹。准静态是可逆的必要条件,但还不够——比如有摩擦的准静态过程就不可逆,因为摩擦生的热收不回来。

    本章的关键对照:凡是真实发生的、有确定速率的过程(热传导、扩散、通电),全都是不可逆的。所以「不可逆过程热力学」这个名字,说的其实就是「真实过程的热力学」。
    落点 · 主干正在把它松开的就是这个前提
    准静态过程 $\;\Longleftrightarrow\;$ 过程速率 $\to 0$

    这一章长什么样,这一份在哪

    第 8 章的内容分成了五份讲义,下面是整章的地图。你现在在第 1 份。

    第 1–7 章平衡态热力学:能不能发生、往哪走
    8.1–8.2 · 本份四条经验定律 → 抽象成「力」与「流」
    8.3熵产生率:把力和流定准
    8.4电动效应
    8.5热电效应
    8.6远离平衡:非平衡态相变

    这一份要走的路只有两步:先看四条实验总结出来的经验定律(8.1),再把它们抽象成一套统一的语言(8.2)。四条定律你多半见过其中一两条,但很可能没意识到它们其实是同一件事。

    第一步 · 四条经验定律(一)

    2. 热往哪流、流多快:傅里叶定律

    0:00 / 0:00

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

    先看最熟悉的一个:一根金属棒,长度 $l$,两端温度分别是 $T_1$ 和 $T_2$,并且 $T_2 > T_1$(左热右冷)。棒的横截面积是 $A$。

    热流:每秒 Q 焦耳 A T₂ (热) T₁ (冷) l
    图 1 热传导。热量从高温端 $T_2$ 流向低温端 $T_1$,每秒穿过端面 $A$ 的热量记作 $\dot Q$。

    实验发现,单位时间内通过截面 $A$ 的热量是:

    $$\dot{Q} = \kappa A\,\frac{T_2 - T_1}{l}$$
    (8.1.1)

    这里 $\kappa$ 叫热导率,是材料自己的性质(铜大、木头小),$\dot Q$ 上面那个点表示「每秒多少」。

    先把这个式子读懂:$\dfrac{T_2-T_1}{l}$ 是「每走一米,温度差多少度」。式子说的是——温度沿途掉得越陡,热流越大;截面越粗,热流越大;材料越导热,热流越大。三件事都符合直觉。

    从「一根棒」升级到「每一点」

    (8.1.1) 有个毛病:它只能描述整根棒。换一根粗细不同、长短不同的棒,就得重新算。我们想要的是一个描述空间中每一点的式子。

    为此定义热流密度:

    $$J_q = \frac{\dot{Q}}{A}$$

    它表示单位时间内流过单位面积的热量(单位 $\mathrm{W/m^2}$)。有了它,(8.1.1) 可以写成不依赖具体几何形状的矢量式:

    $$\boldsymbol{J}_q = -\kappa(T)\,\nabla T$$
    (8.1.2)

    如果只沿 $x$ 一个方向,就退化成:

    $$J_{qx} = -\kappa(T)\,\frac{\mathrm{d}T}{\mathrm{d}x}$$
    (8.1.3)

    这两个式子就是傅里叶定律。

    回顾梯度 ∇:把「往哪个方向变得最快」写成一个矢量
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    温度 $T$ 在空间中每一点都有一个值,写成 $T(x,y,z)$,这叫一个标量场。梯度就是把这个场「哪个方向升得最快、有多快」提取出来:

    $$\nabla T = \Big(\frac{\partial T}{\partial x},\ \frac{\partial T}{\partial y},\ \frac{\partial T}{\partial z}\Big)$$

    它是个矢量,两个信息合在一起:

    • 方向:指向温度上升最快的那一边。(注意是上升,不是下降——负号就是从这来的。)
    • 大小:沿那个方向,每走一米温度升高多少度(单位 $\mathrm{K/m}$)。

    一个最好用的图像:把温度场想成一片山地,高处热、低处冷。那么 $\nabla T$ 就是站在每一点上朝正上坡方向指的一根箭头,坡越陡箭头越长。

    在一维情形,梯度就退化成一个普通导数 $\dfrac{\mathrm{d}T}{\mathrm{d}x}$——(8.1.3) 就是这么来的。

    常见的口误:「温度梯度」不是「温度差」。温度差的单位是 K,梯度的单位是 K/m。两端差 10 K 的棒,长 1 cm 和长 1 m,梯度差 100 倍。
    落点 · 主干正在用它的地方
    $$\boldsymbol{J}_q = -\kappa(T)\,\nabla T$$
    桥从「整根棒」到「每一点」的三步
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 就是上面主干里的 (8.1.1)
    $$\dot{Q} = \kappa A\,\frac{T_2 - T_1}{l}$$
    第 1 步 · 两边除以面积 $A$,摆脱「棒有多粗」

    左边 $\dot Q/A$ 正好就是我们定义的热流密度 $J_q$:

    $$J_q = \kappa\,\frac{T_2 - T_1}{l}$$

    这一步之后,式子里再也没有 $A$ 了——同样的材料,粗棒细棒的 $J_q$ 一样。

    第 2 步 · 把「两端之差」换成「当地的斜率」,摆脱「棒有多长」

    关键是把坐标轴摆好。让 $x$ 轴顺着热流方向,即从热端指向冷端。那么在这个坐标下:

    • 起点($x$ 小的那头)是热端,温度 $T_2$;
    • 终点($x$ 大的那头)是冷端,温度 $T_1$。

    所以沿 $x$ 走,温度是下降的,当地斜率为负:

    $$\frac{\mathrm{d}T}{\mathrm{d}x} = \frac{T_1 - T_2}{l} = -\,\frac{T_2 - T_1}{l} \;<\; 0$$

    注意这里出现的是 $T_1 - T_2$(终点减起点),而 (8.1.1) 里写的是 $T_2 - T_1$。两者刚好差一个负号——负号就是在这一步产生的。

    第 3 步 · 代回去

    把 $\dfrac{T_2-T_1}{l} = -\dfrac{\mathrm{d}T}{\mathrm{d}x}$ 代入第 1 步的结果:

    $$J_q = \kappa \cdot \Big(-\frac{\mathrm{d}T}{\mathrm{d}x}\Big) = -\kappa\,\frac{\mathrm{d}T}{\mathrm{d}x}$$
    终点 · 回到主干的式 (8.1.3)
    $$J_{qx} = -\kappa(T)\,\frac{\mathrm{d}T}{\mathrm{d}x}$$

    再把一维的 $\dfrac{\mathrm{d}T}{\mathrm{d}x}$ 换成三维的 $\nabla T$,就得到主干的 (8.1.2) 了。

    负号的物理意思,一句话:$\nabla T$ 指向热的那边,而热要往冷的那边流。两个方向正好相反,所以中间必须有个负号。这个负号不是数学习惯,它就是热力学第二定律——热自发地从高温流向低温。等号后面的负号一去掉,热就自己从冷流向热了。
    列传流密度 J
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干第一次写下它的地方,原样照抄
    $$J_q = \frac{\dot{Q}}{A}$$
    幕一 · 身世:没有它,哪句话说不出来

    在 $J$ 出生之前,描述传热只能说「这根棒每秒传 5 焦耳」。这句话有个致命毛病:它是整根棒的属性。换一根粗一倍的棒,5 就得改成 10;棒粗细不均匀,这句话干脆没法说;棒里某一处热得快、某一处热得慢,它一个字也表达不了。

    物理学要写的是微分方程,而微分方程说的永远是「空间中这一点此刻怎么样」。所以必须有一个量,它属于点而不属于物体。$\boldsymbol J$ 就是为这件事造出来的:把「总共流过多少」除以面积,「多快」这件事就从物体身上剥离,交给了空间的每一点。

    幕二 · 名片
    读法
    就念字母 $J$。中文叫流密度,也叫通量密度;英文 flux density,或直接叫 flux
    下标
    下标标明「流的是什么东西」:$\boldsymbol J_q$ 热流、$\boldsymbol J_e$ 电流、$\boldsymbol J_n$ 粒子流。$q$ 取自 heat 的传统记号 $Q$,$e$ 取自 electric,$n$ 取自粒子数密度 $n$
    量纲
    统一是「被输运的东西 ÷ 面积 ÷ 时间」。热流密度 $\mathrm{W/m^2}$(即 $\mathrm{J\cdot m^{-2}s^{-1}}$);电流密度 $\mathrm{A/m^2}$;粒子流密度 $\mathrm{m^{-2}s^{-1}}$
    别名
    在电磁学里,同一个模式叫坡印廷矢量(能流密度);在流体力学里叫质量流密度
    幕三 · 脾气:它平时怎么表现
    • 它是矢量,不是数。「每秒流过多少」还不够,还得说清楚往哪个方向流。这是它和一般密度(如 $\rho$、$n$)最大的区别。
    • 它是场量,每一点都有自己的值。和 $T$、$p$ 一样,$\boldsymbol J$ 写出来是 $\boldsymbol J(x,y,z,t)$。
    • 乘回面积和时间,就还原成总量。$\boldsymbol J \cdot \Delta \boldsymbol A \cdot \Delta t$ 就是这段时间内穿过这块面的总量。这条是它和实验测量之间的桥——实验上能测的永远是总量。
    • 它和「守恒」天生绑在一起。如果某处的东西既不凭空产生也不消失,那么某个小盒子里存量的减少,就必须等于从盒子表面流出去的量。这句话写成公式,就是散度定理——第 2 份一开头就要用它。
    幕四 · 去向:以后还会在哪儿碰见它

    本份之内:第 3 段是电流密度 $\boldsymbol J_e$,第 4 段是动量流密度 $P_{zx}$,第 5 段是物质流密度 $\boldsymbol J_n$——四条经验定律,左边全是它。到第 6 段你会看到,这不是巧合。第 8 段起,它会被进一步抽象成一个统称:流。

    本份之外:第 2 份里的 $-\dfrac{\partial(\rho u)}{\partial t} = \mathrm{div}\,\boldsymbol J_q$ 会把它和能量守恒正式绑在一起;第 4 份的热电效应里,$\boldsymbol J_q$ 和 $\boldsymbol J_e$ 会同时出现、互相牵扯。走出本书,它还是电磁学、流体力学、乃至概率论(概率流)的通用零件——这里只点名。

    一句话认人:$\boldsymbol J$ 回答的是「此时此地,每秒每平方米,有多少东西朝哪个方向穿过去」。它把「一共流了多少」这种关于物体的话,翻译成了关于空间每一点的话——物理学能写微分方程,全靠这一步翻译。
    第一步 · 四条经验定律(二)

    3. 欧姆定律,其实是「电量的输运」

    0:00 / 0:00

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

    换一个场景:不加温度差,改成在导体两端加一个电势差。导体里会产生电流——也就是电量在流动。

    实验规律是:

    $$\boldsymbol{J}_e = \sigma \boldsymbol{E} = -\sigma \nabla \varphi$$
    (8.1.4)

    这就是欧姆定律,$\sigma$ 叫电导率。$\boldsymbol J_e$ 是电流密度,$\varphi$ 是电势,$\boldsymbol E$ 是电场强度。

    请把它和上一段的傅里叶定律并排放:
    $\boldsymbol{J}_q = -\kappa\,\nabla T\qquad$ 热流 $=-$ 系数 $\times$ 温度梯度
    $\boldsymbol{J}_e = -\sigma\,\nabla \varphi\qquad$ 电流 $=-$ 系数 $\times$ 电势梯度

    把 $T$ 换成 $\varphi$、$\kappa$ 换成 $\sigma$,两个式子一模一样。这不是巧合,第 6 段会把这件事说透。

    负号的含义也完全平行:$\nabla\varphi$ 指向电势升高的方向,而正电荷要往电势低的方向流,所以差一个负号。

    回顾电场就是电势的「下坡方向」
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    电势 $\varphi$ 是个标量场(每点一个数),电场 $\boldsymbol E$ 是个矢量场(每点一个箭头)。两者的关系正是上一段刚回顾过的梯度:

    $$\boldsymbol{E} = -\nabla \varphi$$

    沿用「山地」的比喻:$\varphi$ 是地形高度,那么 $\boldsymbol E$ 就是每一点上朝正下坡方向指的箭头,坡越陡箭头越长。正电荷放上去就顺着 $\boldsymbol E$ 往下滚——这正是「电势能高的地方往低处走」。

    所以 (8.1.4) 里那个负号,跟热传导里的负号是同一件事在两个场景里的表现:梯度指上坡,东西往下坡走。

    落点 · 主干正在用它的地方
    $$\boldsymbol{J}_e = \sigma\boldsymbol{E} = -\sigma\nabla\varphi$$
    旁通通向电路学:同一条定律的两种写法
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干里的这个式子
    $$\boldsymbol{J}_e = \sigma \boldsymbol{E}$$

    是同一个定律,(8.1.4) 是它的局部形式(说每一点),$U=IR$ 是它的整体形式(说整根导线)。当年做实验先发现的是后者,这里的写法是把它翻译到了每一点上——和第 2 段把 (8.1.1) 翻译成 (8.1.2) 是完全一样的操作。

    验算一下。取一根长 $l$、截面积 $A$、电导率 $\sigma$ 均匀的导线:

    • 两端电压 $U$ 均匀分摊,所以场强 $E = U/l$;
    • 电流是电流密度乘面积,$I = J_e A$。

    代进 $J_e = \sigma E$:

    $$\frac{I}{A} = \sigma\,\frac{U}{l} \;\Longrightarrow\; U = \frac{l}{\sigma A}\,I$$

    括号里那一坨就是电阻:

    $$R = \frac{l}{\sigma A} = \frac{\rho\, l}{A},\qquad \rho = \frac{1}{\sigma}$$

    $\rho = 1/\sigma$ 就是电阻率。于是 $U = IR$ 完整地掉了出来。

    顺带看清一件事:$R$ 里含着 $l$ 和 $A$,所以电阻是这根导线的属性,换根导线就变;而 $\sigma$ 里不含几何量,它是材料的属性。这正是「流密度」列传里说的那件事——除以面积,是为了把物体的属性还原成材料和空间的属性。
    第一步 · 四条经验定律(三)

    4. 最反直觉的一条:流动的是「动量」

    0:00 / 0:00

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

    第三个场景在流体里。设流体沿 $x$ 方向流动,速度是 $v_x$;但不同高度流得不一样快——在 $z$ 方向存在一个速度梯度 $\dfrac{\mathrm{d}v_x}{\mathrm{d}z}$。

    v x 越往上 流得越快 z x
    图 2 黏滞流的速度分布。速度方向都是 $x$,但大小沿 $z$ 变化。

    这时相邻两层流体之间会互相拉扯:快的一层拖着慢的一层加速,慢的一层拖着快的一层减速。上层对下层施加的这个力,大小由牛顿黏滞定律给出:

    $$f = \eta\,\frac{\mathrm{d}v_x}{\mathrm{d}z}\,\Delta S$$
    (8.1.5)

    其中 $\eta$ 叫黏滞系数(水小、蜂蜜大),$\Delta S$ 是两层之间接触的面积。

    关键的转换:把「力」看成「动量的传递」

    到这里为止,(8.1.5) 看起来跟前两条完全不像——前两条讲的是「什么东西在流」,这一条讲的却是「有多大的力」。

    但换个角度看就一样了。力作用一段时间,效果就是传过去一份动量。快的那层把动量交给慢的那层,慢的就变快了。所以黏滞现象本质上是动量的输运。照着 $J = \dfrac{\text{东西}}{\text{面积}\times\text{时间}}$ 的老办法定义动量流密度:

    $$P_{zx} = \frac{\Delta K}{\Delta S\,\Delta t} = \frac{f\cdot \Delta t}{\Delta S\,\Delta t} = \eta\,\frac{\mathrm{d}v_x}{\mathrm{d}z}$$
    (8.1.6)

    $P_{zx}$ 表示:单位时间内、单位面积上传输的动量。两个下标各有分工——

    • 后一个下标 $x$:动量本身的方向是 $x$(流体是沿 $x$ 流的);
    • 前一个下标 $z$:动量传输的方向是 $z$(动量是从上层传到下层的)。
    为什么这条最难接受:前两条里,流动的东西(热、电荷)本身是标量,一个数就说清了。而这里流动的东西是动量,它自己就是矢量。于是「什么在流」和「往哪流」变成了两个独立的方向,非得两个下标才描述得了。理解到这一步就够了,不必更深究。
    这里有一个必须当场说清的细节:(8.1.6) 少了一个负号。

    前两条定律都写成「流 $=-$ 系数 $\times$ 梯度」,这一条却是正的。原因是 (8.1.5) 讨论的是「上层对下层的作用力」,方向已经在文字里交代过,公式只负责给大小——这是流体力学里的常见写法。

    但物理上,动量的流向和前两条完全一致:从 $v_x$ 大的地方流向 $v_x$ 小的地方。严格按矢量约定写出来应当是
    $$P_{zx} = -\,\eta\,\frac{\mathrm{d}v_x}{\mathrm{d}z}$$
    请记住这个补回来的负号,第 6 段要靠它,四条定律才能真正对齐。
    桥从「力」到「动量流密度」
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 就是上面主干里的 (8.1.5)
    $$f = \eta\,\frac{\mathrm{d}v_x}{\mathrm{d}z}\,\Delta S$$
    第 1 步 · 请出牛顿第二定律的另一种写法

    中学最熟的写法是 $f = ma$。但牛顿本人写的其实是动量形式:力等于动量的变化率。

    $$f = \frac{\mathrm{d}K}{\mathrm{d}t}\qquad(K \text{ 是动量})$$

    两边乘上一小段时间 $\Delta t$:

    $$\Delta K = f\cdot \Delta t$$

    这就是动量定理:力乘上作用时间,等于传过去的动量。所以「上层对下层施加力 $f$,持续 $\Delta t$」这句话,可以原地翻译成「上层向下层交付了 $f\Delta t$ 这么多动量」。这两句话说的是同一件事,只是换了个说法。

    第 2 步 · 套用流密度的通用定义

    流密度的定义永远是同一个模子:

    $$\text{流密度} = \frac{\text{传过去的东西}}{\text{面积}\times\text{时间}}$$

    这次「传过去的东西」是动量 $\Delta K$,接触面积是 $\Delta S$,时间是 $\Delta t$:

    $$P_{zx} = \frac{\Delta K}{\Delta S\,\Delta t}$$
    第 3 步 · 把第 1 步的 $\Delta K = f\Delta t$ 代进去,$\Delta t$ 自动约掉
    $$P_{zx} = \frac{f\cdot\Delta t}{\Delta S\,\Delta t} = \frac{f}{\Delta S}$$

    再把起点的 (8.1.5) 除以 $\Delta S$ 代入:

    终点 · 回到主干的式 (8.1.6)
    $$P_{zx} = \eta\,\frac{\mathrm{d}v_x}{\mathrm{d}z}$$

    顺带一个有意思的观察:中间那步 $P_{zx} = f/\Delta S$——动量流密度就是单位面积上的力,也就是压强的量纲($\mathrm{Pa}$)。这不是巧合。压强本来就是气体分子把动量砸到器壁上的结果,它本身就是一种动量流密度。

    易错(8.1.6) 缺的那个负号
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干里的这个式子
    $$P_{zx} = \eta\,\frac{\mathrm{d}v_x}{\mathrm{d}z}$$

    这是个很值得警惕的地方,因为它看起来像是打破了规律,其实没有。

    正确的理解:(8.1.6) 写的是大小,负号被省掉了。物理上动量的流向和前两条完全一致——

    $$\text{动量从 } v_x \text{ 大的地方} \;\longrightarrow\; v_x \text{ 小的地方}$$

    也就是顺着速度梯度往下坡走。跟热从高温流向低温、电荷从高电势流向低电势,是同一个模式。严格按矢量约定写出来,应当是:

    $$P_{zx} = -\,\eta\,\frac{\mathrm{d}v_x}{\mathrm{d}z}$$

    错误的理解(要避免的):「黏滞和前两条不一样,它没有方向性上的限制。」这是错的。真要没有负号,就意味着动量会自发从慢的一层流向快的一层——快的越来越快、慢的越来越慢,两层流体自己分道扬镳。这直接违反热力学第二定律,现实中永远不会发生。

    为什么课本这里不写负号:因为 (8.1.5) 讨论的是「上层对下层的作用力」,方向已经在文字里交代过了,公式只负责给大小。这在流体力学里是常见写法。但你心里要清楚它该有负号——到第 6 段总结四条定律的共同形式时,正是靠这个负号,四条才能真正对齐。
    实例同样的动作,水和甘油差 1400 倍
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干的式 (8.1.5),原样照抄
    $$f = \eta\,\frac{\mathrm{d}v_x}{\mathrm{d}z}\,\Delta S$$

    设一个具体装置:两块面积 $\Delta S = 0.01\ \mathrm{m^2}$(10 cm × 10 cm)的平板,间距 $1\ \mathrm{mm}$,中间灌满液体。让上板以 $0.1\ \mathrm{m/s}$ 平移,下板固定。

    速度梯度:

    $$\frac{\mathrm{d}v_x}{\mathrm{d}z} = \frac{0.1\ \mathrm{m/s}}{0.001\ \mathrm{m}} = 100\ \mathrm{s^{-1}}$$

    20 ℃ 下的黏滞系数大致是:水 $\eta \approx 1.0\times10^{-3}\ \mathrm{Pa\cdot s}$,甘油 $\eta \approx 1.4\ \mathrm{Pa\cdot s}$。代进去:

    液体$\eta\ (\mathrm{Pa\cdot s})$需要的力 $f$
    水$1.0\times10^{-3}$$1.0\times10^{-3}\times100\times0.01 = 1\ \mathrm{mN}$
    甘油$1.4$$1.4\times100\times0.01 = 1.4\ \mathrm{N}$

    1 毫牛大约是一只蚊子的重量,1.4 牛大约是一瓶 140 克矿泉水的重量——相差 1400 倍。这就是「在水里划手很轻松、在甘油里划手很费劲」的定量版本。

    顺带验证了主干的一件事:$\eta$ 的单位是 $\mathrm{Pa\cdot s}$,也就是「压强 × 时间」。这正好印证桥里那个观察——$P_{zx} = \eta\,\dfrac{\mathrm{d}v_x}{\mathrm{d}z}$ 的左边是压强量纲,右边 $\dfrac{\mathrm{d}v_x}{\mathrm{d}z}$ 的单位是 $\mathrm{s^{-1}}$,所以 $\eta$ 必须是 $\mathrm{Pa\cdot s}$。量纲对上了,说明式子没写错。
    第一步 · 四条经验定律(四)

    5. 一滴墨水化开:菲克定律

    0:00 / 0:00

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

    最后一个场景:一杯清水,滴一滴墨水进去,不搅拌。墨水会自己慢慢散开,直到浓度处处相同。

    这里存在的是浓度梯度。设 $n$ 是某种组元的粒子数密度(单位体积里有多少个粒子),实验规律是:

    $$\mathrm{d}M = -D\Big(\frac{\mathrm{d}n}{\mathrm{d}z}\Big)\,\mathrm{d}S\,\mathrm{d}t$$

    写成流密度的形式,就是菲克(Fick)定律:

    $$\boldsymbol{J}_n = \frac{\mathrm{d}M}{\mathrm{d}S\,\mathrm{d}t} = -D\,\nabla n$$
    (8.1.7)

    $\boldsymbol J_n$ 是粒子流密度(单位时间内流过单位面积的粒子数),$D$ 叫扩散系数。这个现象就叫扩散。

    负号还是老意思:$\nabla n$ 指向浓度高的地方,而粒子要往浓度低的地方跑。墨水从浓处流向淡处,不会反过来自己聚回一滴。

    这里要预先埋一个伏笔,因为它是后面几份的关键:

    浓度梯度并不是扩散的唯一原因。在化学反应中,化学势 $\mu$ 的梯度同样会产生物质流。事实上到第 2 份你会看到,$\nabla \mu$ 才是更本质的那个原因——浓度梯度只是它在特殊情形下的表现。这里先记住有这回事。
    回顾三个都叫「多少」的量,分清楚
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    量意思单位是不是场
    $N$ 粒子数
    $N_i$ 摩尔数
    整个系统里一共有多少个 / mol否,属于整个系统
    $n$ 粒子数密度某一点附近单位体积里有多少$\mathrm{m^{-3}}$是,每点一个值
    $x_i$ 摩尔分数第 $i$ 种占总量的比例纯数,无量纲可以是

    为什么 (8.1.7) 里必须用 $n$ 而不是 $N$:因为要写梯度。梯度是「每走一米变化多少」,只有逐点有定义的量才谈得上梯度。$N$ 是整个系统的一个数,它没有「在这一点的值」,也就没有梯度。

    这跟第 2 段用 $J_q$ 而不用 $\dot Q$、第 3 段用 $\sigma$ 而不用 $R$,是同一个动作的第三次出现:凡是要写成场方程的,都得先把「整体量」换成「逐点量」。

    落点 · 主干正在用它的地方
    $$\boldsymbol{J}_n = -D\,\nabla n$$
    实例算一下,你会明白为什么泡茶要搅
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干的式 (8.1.7),原样照抄
    $$\boldsymbol{J}_n = -D\,\nabla n$$

    要估时间,需要借用一个结论:扩散在时间 $t$ 内大致走过的距离是

    $$L \sim \sqrt{2Dt}$$
    这个式子本身要到统计物理(随机行走)才推得出来,这里只借它来看数量级,不必深究。注意它的形状:距离和时间的平方根成正比,不是正比——这是扩散最要命的特点。

    常温水中,小分子(墨水、糖、盐)的扩散系数大致 $D \approx 1\times10^{-9}\ \mathrm{m^2/s}$。反过来解出走 $L$ 需要的时间 $t = L^2/(2D)$:

    要走的距离 $L$所需时间 $t = L^2/2D$体感
    $1\ \mathrm{mm}$约 $500\ \mathrm{s}$8 分钟
    $1\ \mathrm{cm}$约 $5\times10^{4}\ \mathrm{s}$14 小时
    $10\ \mathrm{cm}$(一杯水)约 $5\times10^{6}\ \mathrm{s}$近两个月

    距离涨 10 倍,时间涨 100 倍。所以纯靠扩散把一杯茶泡匀,要等两个月;而搅一下只要两秒——因为搅拌靠的是对流(整块液体搬运),不是扩散。

    这也解释了为什么细胞必须很小(几十微米量级):细胞内部的物质运输全靠扩散,尺寸一大,营养还没扩散到中心就先饿死了。

    第二步 · 看出它们是同一件事

    6. 四条定律,其实是一条

    0:00 / 0:00

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

    把前四段的结果并排放在一起(黏滞那条,用第 4 段末尾补回负号之后的写法):

    定律什么在流什么梯度在推式子系数
    傅里叶热(能量)温度 $T$$\boldsymbol{J}_q = -\kappa\nabla T$热导率 $\kappa$
    欧姆电荷电势 $\varphi$$\boldsymbol{J}_e = -\sigma\nabla\varphi$电导率 $\sigma$
    牛顿黏滞动量速度 $v_x$$P_{zx} = -\eta\,\dfrac{\mathrm{d}v_x}{\mathrm{d}z}$黏滞系数 $\eta$
    菲克物质(粒子)浓度 $n$$\boldsymbol{J}_n = -D\nabla n$扩散系数 $D$

    四条定律来自四个毫不相干的实验领域——热学、电学、流体力学、化学。但它们的骨架完全相同:

    $$\underbrace{\text{流}}_{\boldsymbol J} \;=\; -\,\underbrace{\text{系数}}_{\kappa,\ \sigma,\ \eta,\ D} \;\times\; \underbrace{\text{某个量的梯度}}_{\nabla T,\ \nabla\varphi,\ \nabla v_x,\ \nabla n}$$

    三个共同点,每一个都值得单独说一句:

    • 左边永远是一个「流」。虽然流的东西五花八门(能量、电荷、动量、粒子),但都被写成了同一种形式:单位时间、单位面积上过去多少。
    • 右边永远是一个「梯度」。没有梯度就没有流——温度处处相同,热就不流;浓度处处相同,墨水就不动。梯度是「原因」,流是「结果」。
    • 永远带负号,而且系数永远为正。东西总是从高处流向低处,从来不反过来。这不是四个巧合,这是热力学第二定律在四个场景里的四张面孔。

    还有第四个共同点,它不像上面三个那么显眼,但后面几段全靠它:四条定律都是正比关系——梯度大一倍,流就大一倍。

    这个「正比」并不是天经地义的,它其实是「偏离平衡不远」这个前提的直接后果。道理很短:流 $J$ 是梯度 $X$ 的某个未知函数 $J(X)$,但我们知道 $X=0$ 时必有 $J=0$(没有梯度就没有流)。把它在 $X=0$ 附近展开:

    $$\begin{aligned}J(X) &= \underbrace{J(0)}_{=\,0} + J'(0)\,X + O(X^2)\\[3pt] &\approx L\,X \qquad (X \text{ 很小时})\end{aligned}$$

    头一项因「无梯度无流」而消失,剩下的第一项就是线性项。所以:

    • 近平衡(梯度小):$O(X^2)$ 可以扔掉,线性关系成立——8.1 到 8.5 全在这一档;
    • 远平衡(梯度大):高次项不能扔了,线性失效——这正是 8.6 的领域。

    这也顺带解释了 (8.1.2) 里为什么把热导率写成 $\kappa(T)$ 而不是常数:线性说的是「$J$ 正比于 $X$」,并没有说系数必须是常数,切线斜率本身可以随当地温度变化。

    看出这一点,就看出了这一整章的立足点:既然四条经验定律长得一模一样,那它们背后一定有一条更普遍的规律,只是被四个学科各自发现了一次、各自起了个名字。8.2 节要做的,就是把这条共同的规律直接写出来——用一套不区分热、电、动量、物质的抽象语言。
    溯源统一了以后,能问出原来问不出的问题
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干刚看出来的共同骨架
    $$\text{流} = -\,\text{系数}\times\text{梯度}$$

    如果只是想算热传导,那确实用傅里叶定律就够了,统一没有任何好处。统一的价值在别处:它让一类原本提不出来的问题变得可以提。

    注意四条定律都是「一个梯度引起一个流」的形式。可现实中,好几种梯度常常同时存在。比如一根金属棒,两端既有温差、又通着电——这时会怎样?

    • 温度梯度产生热流,电势梯度产生电流。这是四条定律各自管的部分。
    • 但温度梯度会不会也产生电流?电势梯度会不会也产生热流?

    会。这正是第 4 份要讲的热电效应——你的手机充电宝可能就用它测温度,热电偶温度计、热电制冷片都是它。四条定律各自为政的时候,这种「交叉」的问题根本没有位置可以提出来:傅里叶定律里没有电,欧姆定律里没有热。

    把四条写成同一种形式之后,「交叉项」就成了一个自然而然的问题:既然流和梯度可以配对,那能不能让每个流对每个梯度都有响应?8.2 节的线性唯象方程就是这么写出来的。

    一句话总结这个动机:统一不是为了好看,是为了把「交叉效应」这件事变得可以言说。这也是整章从 8.2 到 8.5 的主线。
    深入线性区其实宽得惊人
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干里那个「$X$ 很小时」的条件
    $$J(X) \approx L\,X \qquad (X \text{ 很小})$$

    主干说「梯度小的时候线性成立」,但「小」是相对于什么小?梯度的单位是 K/m,一个有单位的量本身无所谓大小,必须和某个同单位的参照量比。

    参照量藏在微观尺度上。真正的判据是:走过一个分子平均自由程 $\lambda$ 的距离,温度变化了多少,相比温度本身算不算大:

    $$\varepsilon = \frac{\lambda\,|\nabla T|}{T} \;\ll\; 1\ ?$$

    为什么是自由程?因为分子每走一个自由程就碰一次、重新「随大流」一次。只要这段路上环境几乎没变,它感受到的就是一个局部平衡的世界。

    代具体数字(常温常压空气):$\lambda \approx 68\ \mathrm{nm} = 6.8\times10^{-8}\ \mathrm{m}$,$T \approx 300\ \mathrm{K}$。取一个相当陡的温度梯度——比如暖气片旁边,$|\nabla T| = 100\ \mathrm{K/m}$:

    $$\varepsilon = \frac{6.8\times10^{-8}\times 100}{300} \approx 2\times10^{-8}$$

    两千万分之一。也就是说,日常生活里所有的传热现象,在分子看来都紧贴着平衡态,线性关系好得不能再好。这解释了一件本来挺费解的事:为什么傅里叶定律这种「一次近似」,实际用起来精度却高得像是精确定律。

    要把 $\varepsilon$ 推到 1 附近,需要在 68 纳米内造出 300 K 的温差——那是激波内部、飞秒激光烧蚀这一类极端场合。

    那 8.6 节的「远离平衡」是怎么回事?——这里必须说准,否则会误解:

    8.6 里的贝纳德对流并不是傅里叶定律失效了。失效的是整个系统的响应:温度梯度大到一定程度,流体的浮力压过了黏滞阻力,系统整体跳到了另一种运动方式(对流卷)。非线性来自各个过程之间的耦合与反馈,不是来自某一条唯象定律本身垮掉。

    换句话说,「远离平衡」说的是系统层面离平衡远,而不是每个小区域内部离平衡远——局部平衡假设在贝纳德对流里依然成立。

    8.1 小结

    这一半的逻辑链条是:热力学不管快慢 → 先看实验总结的四条经验定律(傅里叶、欧姆、牛顿黏滞、菲克)→ 每条都把「整体量」换成「逐点的流密度 $\boldsymbol J$」→ 发现四条的骨架完全相同:流 $=-$ 系数 $\times$ 梯度 → 负号统一来自热力学第二定律,线性来自「偏离平衡不远」。接下来的 8.2 节,要把这个共同骨架抽象成一套通用语言:力与流。

    第三步 · 搭框架(一)· 先解决资格问题

    7. 非平衡态里,凭什么还能说「温度」?

    0:00 / 0:00

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

    8.2 节要做的事,第 6 段末尾已经交代了:把四条定律的共同骨架抽象成通用语言。但在动手之前,必须先解决一个资格问题——而且这个问题相当致命。

    致命的地方在这:第 1 段那条限制说得很清楚——$T$、$S$、$\mu$ 这些量只对平衡态才有定义。可这一章讨论的偏偏全是非平衡态:一根两端温度不同的棒,整体根本不处于平衡。

    那么 (8.1.2) 里的 $\nabla T$ 是什么意思?「温度的梯度」这句话,得先承认「每一点都有温度」才说得出口。而这恰恰是被禁止的。

    这一章绕开它的办法,叫局部平衡假设。它是整个非平衡热力学的地基,只有一句话:

    局部平衡假设:把整个系统切成许多宏观小、微观大的小区域。整体虽然不平衡,但每一个小区域可以近似看成处于热力学平衡态,因而可以用热力学的那一套量去描述它、可以定义它的热力学函数。

    「宏观小、微观大」这六个字是关键,它要求小区域的尺寸卡在一个夹缝里:

    • 宏观小——小到在这个区域内部,温度、压强等等基本没有变化,可以当成一个均匀的东西;
    • 微观大——又要大到里面装着足够多的分子(比如 $10^{10}$ 个以上)。因为温度、压强本质上是大量分子的统计平均,只有三五个分子的地方谈不上「温度」。

    在这个假设下,每一个小区域里,吉布斯方程仍然成立:

    $$T\,\mathrm{d}S = \mathrm{d}U + p\,\mathrm{d}V - \sum_i \mu_i\,\mathrm{d}N_i$$

    于是 $T$、$p$、$N_i$、$\mu_i$ 以及所有热力学函数,都从「整个系统的一个数」变成了空间和时间的函数:$T(x,y,z,t)$、$\mu_i(x,y,z,t)$……这正是我们要的。

    再利用熵的可加性(整体的熵等于各部分之和),整个系统的熵就可以写成积分:

    $$S = \int s(x,y,z,t)\,\mathrm{d}x\,\mathrm{d}y\,\mathrm{d}z$$

    这里 $s$ 是熵密度——单位体积(或单位质量)的熵。请记住这个 $s$,第 2 份从头到尾都在算它。

    回顾热力学基本方程,逐项拆开
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    它通常写成 $\mathrm{d}U = T\mathrm{d}S - p\mathrm{d}V + \sum_i\mu_i\mathrm{d}N_i$,移项就是主干那个形式。这个式子说的是:系统的内能可以从三条渠道改变。

    项物理意思怎么记
    $T\,\mathrm{d}S$以热的方式吸收的能量可逆过程中 $\delta Q = T\mathrm{d}S$
    $-p\,\mathrm{d}V$以做功的方式吸收的能量体积被压缩($\mathrm{d}V<0$)时内能升高
    $\sum_i\mu_i\,\mathrm{d}N_i$因为物质进出而带来的能量每加进一点第 $i$ 种物质,内能涨 $\mu_i$

    化学势 $\mu_i$ 是这里最需要盯住的一个量,因为第 2 份起它就是主角。它的定义是:

    $$\mu_i = \Big(\frac{\partial G}{\partial N_i}\Big)_{T,\,p}$$

    意思是:在温度、压强不变的条件下,往系统里再加一点点第 $i$ 种物质,吉布斯函数 $G$ 增加多少。单位是 $\mathrm{J/mol}$。

    最好用的直觉是把它类比成温度和压强:

    • 温度不等 $\Rightarrow$ 热会流动,从高温流向低温;
    • 压强不等 $\Rightarrow$ 体积会变化,从高压推向低压;
    • 化学势不等 $\Rightarrow$ 物质会流动,从高 $\mu$ 流向低 $\mu$。

    三者是同一张桌子上的三个位置。第 2 份会把最后一条严格证明出来(现在只是类比)。

    落点 · 主干正在用它的地方
    $$T\,\mathrm{d}S = \mathrm{d}U + p\,\mathrm{d}V - \sum_i \mu_i\,\mathrm{d}N_i$$
    易错矛盾是真的,但它是可控的
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干里的这个假设
    整体非平衡态 $\;+\;$ 每个小区域平衡态

    先承认:这确实是一对矛盾的概念,课本自己也这么说。不要试图说服自己它其实不矛盾——那样反而学不懂。正确的态度是搞清楚它在什么条件下矛盾得足够轻微,轻微到可以忽略。

    正确的理解:小区域并不是真的平衡。它和邻居时时刻刻在交换能量和物质,所以严格说它也在变化。假设要求的是——

    $$\frac{\text{小区域与邻居交换的速率}}{\text{小区域内部自己弛豫到平衡的速率}} \;\ll\; 1$$

    分母那个「内部弛豫速率」由分子碰撞决定,快得惊人:常温常压气体里分子平均每 $10^{-10}$ 秒碰一次,小区域内部在纳秒量级就重新平衡了。而分子那个宏观过程(热传导、扩散)通常在秒的量级。两者相差十个数量级以上,所以这个比值确实远小于 1。

    换句话说:小区域内部「跑得飞快」,一直追得上外界的变化,所以在任何一个瞬间看它,它都已经稳定下来了。这就是课本那句「变化过程要很缓慢,仍可看成是准静态过程」的实际含义——不是宏观过程真的慢,而是它相对于分子碰撞慢。

    错误的理解(要避免的):「既然小块是平衡的,那把它们拼起来整体也该是平衡的。」错在平衡态不只要求每块自己稳定,还要求块与块之间没有差别。一根左端 80 ℃、右端 20 ℃ 的棒,每一小段自己都有确定的温度(局部平衡),但相邻段温度不同,所以整体不平衡——正是这个「不同」造成了热流,也正是本章要研究的东西。

    那它什么时候会失效?当宏观变化快到能和分子碰撞比拼的时候。例如激波(一层几个分子自由程厚的突变面)、超快激光加热(飞秒量级)、极稀薄气体(分子几乎不碰撞)。这些场合下本章的整套理论都不适用,需要另起炉灶。知道一个理论的边界在哪,和知道它说什么同样重要。
    第三步 · 搭框架(二)· 给两类角色起名

    8. 两个新名字:力 X 和流 J

    0:00 / 0:00

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

    资格问题解决了,现在可以正式抽象。

    系统处于非平衡态,必然意味着某些量在空间上不均匀——也就是必然出现梯度:温度梯度、速度梯度、电势梯度、浓度梯度。而一旦有梯度,体系中就会出现相应的流:热流、动量流、电流、物质流。

    非平衡热力学干脆给这两类东西各起一个统称:

    记号叫什么指的是在第 6 段那张表里是
    $X$力(广义力、热力学力)温度梯度、速度梯度、电势梯度、浓度梯度……右边那一栏
    $J$流(广义流、热力学流)热流、动量流、电流、物质流……左边那一栏

    于是四条经验定律的共同骨架,可以用一行字写完:

    $$\text{力 } X \;\longrightarrow\; \text{流 } J$$
    为什么这一步值得单独占一段:因为它把讨论的对象换了。原来我们在讨论「热怎么传、电怎么导」这些具体现象;现在我们讨论的是「力和流之间是什么关系」这个抽象问题。抽象问题的答案,会自动对所有具体现象都成立。

    而且更要紧的是——只有站在这个高度上,「一种力能不能引起另一种流」这个问题才有地方安放。四条定律各管一摊的时候,这个问题连提都提不出来:傅里叶定律里没有电,欧姆定律里没有热。

    关于「力」这个词,课本有一句提醒值得抄下来:「这样做和力学吻合,但比力学上的概念要广泛得多。」「吻合」在于两边都是梯度驱动——力学里最典型的力也都是势的梯度 $\boldsymbol F = -\nabla V$。而「广泛得多」体现在两处:

    • 热力学力不必是矢量。比如化学反应的「力」是化学亲合势 $A=-\sum_i\mu_i\nu_i$(第 7 章出现过),那是个标量。力学里没有标量的力。
    • 热力学力不必有统一的量纲。温度梯度是 $\mathrm{K/m}$,电势梯度是 $\mathrm{V/m}$,量纲各不相同,却都叫力。它们唯一的共同点是扮演同一个角色:造成流。
    宽松是有代价的,这里先把它挑明,第 11 段要专门处理:既然「力」只按角色定义、量纲又不统一,那么同一个物理情形可以写出好几套不同的「力」。

    就拿热传导来说,温度不均匀这件事,可以写成 $\nabla T$,也可以写成 $\nabla\!\left(\dfrac{1}{T}\right)$——两者都在温度均匀时为零,都是合格的候选。那到底该选哪一个?会不会随便选都行?
    溯源「力」这个词是怎么被借过来的
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干里的这个命名
    温度梯度、浓度梯度…… $\;\equiv\;$ 力 $X$

    先承认你的疑虑是对的:$\nabla T$ 的单位是 $\mathrm{K/m}$,跟牛顿 $\mathrm{N}$ 一点关系都没有。它不是力学意义上的力。

    这个词是按「角色」借的,不是按「量纲」借的。看看力学里,力扮演什么角色:

    力学非平衡热力学
    没有它时物体静止不动系统处于平衡,什么都不流
    有了它物体开始运动系统出现流
    它越大运动越剧烈流越大
    它来自势能的梯度 $\boldsymbol F=-\nabla V$某个热力学量的梯度

    最后一行是这个类比真正的落点。力学里最典型的力(重力、弹力、静电力)都是势的梯度:$\boldsymbol F = -\nabla V$。而热力学力也全是梯度。两边在「梯度驱动」这件事上结构完全一样,所以借这个词并不牵强。

    还有一层更实际的理由,跟历史有关。这套语言是昂萨格、普利高津他们在 1930–1950 年代建立的,而当时已经有一个现成的范本摆在那里:力学里「力做功」的整套形式体系。借用「力」这个词,等于宣布「我们要照着那套模式来处理不可逆过程」。

    照着的具体是哪一点?是配对相乘得到一个有物理意义的量这个模式:

    领域配对相乘得到
    力学力 $\times$ 位移功(能量)
    热力学$T\,\mathrm{d}S$、$-p\,\mathrm{d}V$、$\mu_i\,\mathrm{d}N_i$能量
    非平衡热力学力 $X\;\times$ 流 $J$熵产生率

    最后一行就是这套语言真正的落点。第 11 段和整个第 2 份,都在围着它转——正是「$X$ 与 $J$ 相乘应当给出熵产生率」这个要求,最终决定了力必须怎么写。

    旁通通向分析力学与热力学的共同模式
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干里的力与流
    $$X \;\longrightarrow\; J$$

    见过,而且你早就用过了——就在吉布斯方程里。回头看第 7 段那个式子:

    $$\mathrm{d}U = T\,\mathrm{d}S - p\,\mathrm{d}V + \sum_i \mu_i\,\mathrm{d}N_i$$

    右边每一项都是「某个强度量 × 某个广延量的微分」:

    项强度量(广义力)广延量(广义位移)
    $T\,\mathrm{d}S$温度 $T$熵 $S$
    $-p\,\mathrm{d}V$压强 $p$体积 $V$
    $\mu_i\,\mathrm{d}N_i$化学势 $\mu_i$物质的量 $N_i$
    力学里的 $\boldsymbol F\cdot \mathrm{d}\boldsymbol r$力 $\boldsymbol F$位移 $\boldsymbol r$

    这种成对出现的量叫共轭对。每一对相乘都得到能量,而每一对里总是一个强度量、一个广延量。物理学里这个模式随处可见(电学的 $\varphi\,\mathrm{d}q$、磁学的 $\boldsymbol H\cdot\mathrm{d}\boldsymbol B$、表面的 $\gamma\,\mathrm{d}A$……)。

    所以那边解决了什么问题,这边也想解决同一个问题:共轭配对的价值在于,只要认出谁和谁配对,就知道该往式子里写什么,不用每个场景重新想一遍。

    非平衡热力学的 $X$ 与 $J$ 正是这个模式的速率版:那边是 $\sum(\text{力})\times(\text{位移})=\text{能量}$,这边则是 $\sum X_j J_j = \text{熵产生率}$。第 11 段和第 2 份整份,都在围着这个式子转。
    第三步 · 搭框架(三)· 写下总方程

    9. 让每个流对每个力都有反应

    0:00 / 0:00

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

    有了 $X$ 和 $J$ 的名字,就可以写出这一章的总方程了。

    近平衡情况下(第 6 段末尾已经论证过为什么是线性的),假如系统里同时存在一组力 $X_1, X_2, \dots, X_k$,那么第 $j$ 个流是:

    $$J_j = \sum_{i=1}^{k} L_{ji}\,X_i$$
    (8.2.1)

    这就是线性唯象方程,$L_{ji}$ 叫唯象系数(也叫动理系数)。

    请注意这一步真正的动作:不是「第 $j$ 个流由第 $j$ 个力决定」,而是每一个流都对所有的力有反应。这就是第 8 段末尾说的那件事——四条定律各自为政时提不出来的「交叉」问题,现在被直接写进了方程。

    拆开看:两个力的情形

    抽象的求和号不好体会,取 $k=2$ 写开:

    $$\begin{aligned}J_1 &= L_{11}X_1 + L_{12}X_2\\[2pt] J_2 &= L_{21}X_1 + L_{22}X_2\end{aligned}$$

    如果这两个力具体取为 $X_1 = \nabla T$(温度梯度)、$X_2 = \nabla n$(浓度梯度),那么:

    $$J_1 = L_{11}\nabla T + L_{12}\nabla n,\qquad J_1 \text{ 为热流}$$
    (8.2.2)
    $$J_2 = L_{21}\nabla T + L_{22}\nabla n,\qquad J_2 \text{ 为物质流}$$
    (8.2.3)

    四个系数,两个是「老熟人」,两个是全新的:

    项物理意义身份
    $L_{11}\nabla T$温度梯度产生的热流就是傅里叶定律($L_{11}=-\kappa$)
    $L_{22}\nabla n$浓度梯度产生的物质流就是菲克定律($L_{22}=-D$)
    $L_{12}\nabla n$浓度梯度产生的物质流携带内能,从而引起热流交叉项:Dufour 效应(1872 年发现)
    $L_{21}\nabla T$温度梯度引起的热扩散,从而引起物质流交叉项:Soret 效应(1893 年发现)
    这就是这套理论的第一份「利润」:对角项 $L_{11}$、$L_{22}$ 只是把已知的两条定律重写了一遍,没有新东西。真正的新内容全在非对角项 $L_{12}$、$L_{21}$ 上——它们对应的两个效应,是这套框架搭起来之后才被系统地认识的现象。「这样从热力学角度推广了一些定律」,指的就是这件事。
    实例温差把混合物「分」开
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干的式 (8.2.3),原样照抄
    $$J_2 = L_{21}\nabla T + L_{22}\nabla n$$

    看一个具体的实验安排:一根装着均匀盐溶液的管子,两端保持不同温度,其余什么都不做。

    按照第 5 段的菲克定律,浓度已经均匀了,$\nabla n = 0$,物质流应该为零、什么都不该发生。但 (8.2.3) 说还有一项 $L_{21}\nabla T$ —— 只要有温差,物质就会流。

    于是溶质开始朝一头聚集,浓度变得不均匀,$\nabla n$ 慢慢长出来。等到两项恰好抵消:

    $$L_{21}\nabla T + L_{22}\nabla n = 0$$

    物质流停止,管子里就稳定地维持着一个浓度差——温差把原本均匀的混合物分开了一点。

    有多大?习惯上用索雷特系数 $S_T$ 描述稳定后的效果:

    $$\frac{\Delta c}{c} \approx S_T\,\Delta T$$

    多数水溶液的 $S_T$ 在 $10^{-3}\sim10^{-2}\ \mathrm{K^{-1}}$ 量级。取 $S_T = 5\times10^{-3}\ \mathrm{K^{-1}}$、$\Delta T = 20\ \mathrm{K}$:

    $$\frac{\Delta c}{c} \approx 5\times10^{-3}\times 20 = 0.1 = 10\%$$

    百分之十的浓度差,这不是可以忽略的小数。

    它有过一个很硬的用途:克鲁修斯–迪克尔热扩散塔——一根竖管,中心一根热丝,管壁冷却。靠 Soret 效应加上对流反复富集,可以分离同位素(比如铀 235 和铀 238)。同位素的化学性质几乎完全相同,常规化学方法分不开,而这个交叉效应能分开。

    另外,$L_{12}$、$L_{21}$ 这两个交叉系数不是各自独立的——下一段的昂萨格关系会把它们锁死在一起。
    深入因为「载体」是同一批分子
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干里的交叉项
    $$J_1 = L_{11}X_1 + \underbrace{L_{12}X_2}_{\text{交叉项}}$$

    交叉项之所以存在,根子在一件很朴素的事上:搬运不同东西的,是同一批分子。

    拿 Dufour 效应($L_{12}$)来说。浓度梯度让分子从浓处往淡处跑,这是物质流。但分子不是空手走的——每个分子身上带着自己的动能和势能,也就是内能。分子跑过去,内能也就跟着搬了过去。而「内能的搬运」正是热流的定义。于是一个纯粹由浓度梯度驱动的过程,顺带产生了热流。

    反过来的 Soret 效应($L_{21}$)也一样:温度高的地方分子运动剧烈、跑得更远,两种分子「跑得更远」的程度又不一样(轻的比重的更敏感),于是净效果就是某一种分子朝一头偏移——温度梯度产生了物质流。

    由此可以推出一条很好用的判断法则:两个过程之间有没有交叉效应,取决于它们是不是共用同一批载体。

    金属里的导电和导热都由自由电子承担 $\Rightarrow$ 有强烈的交叉效应(这就是第 4 份的热电效应,也是第 2 份第 4 个情形要算的那个)。
    而热传导(电子)和黏滞(原子实)载体不同 $\Rightarrow$ 基本没有交叉。

    还有一层更深的限制,这里只点名:居里原理——在各向同性的介质中,张量性质(阶数)不同的力和流之间不能耦合。所以标量的化学反应速率,不能和矢量的热流直接耦合。这解释了为什么并非所有 $L_{ij}$ 都非零。要展开这一条需要张量语言,超出本章范围。

    第三步 · 搭框架(四)· 这一章最漂亮的结论

    10. 昂萨格倒易关系:$L_{12} = L_{21}$

    0:00 / 0:00

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

    上一段留下四个系数 $L_{11}, L_{12}, L_{21}, L_{22}$。看起来它们应该各是各的、互不相干。

    但它们不是。在 $J_1$ 和 $J_2$ 的系数之间存在一个关系,叫昂萨格(Onsager)倒易关系:

    $$L_{21} = L_{12}$$

    写成一般形式:

    $$L_{ij} = L_{ji}$$
    (8.2.4)

    在有磁场 $\boldsymbol B$ 的情况下,要改写成:

    $$L_{ij}(\boldsymbol{B}) = L_{ji}(-\boldsymbol{B})$$
    (8.2.5)
    先把 (8.2.4) 用大白话说一遍:
    「温度梯度引起物质流的本事」,恰好等于「浓度梯度引起热流的本事」。

    这两件事听上去毫不相干——一个是 1893 年发现的,另一个是 1872 年发现的,隔了二十年,而且是在不同的实验里。但它们的系数严格相等。这不是近似,是等式。

    它值多少钱

    用一句话概括昂萨格关系的作用:它把独立的实验测量减少了一半。

    • $k$ 个力、$k$ 个流,系数矩阵 $L_{ij}$ 有 $k^2$ 个元素;
    • 有了 $L_{ij}=L_{ji}$,矩阵是对称的,独立元素只剩 $\dfrac{k(k+1)}{2}$ 个。

    $k=2$ 时,4 个变 3 个;$k=3$ 时,9 个变 6 个。省下来的不是纸面上的符号,是实打实的实验——每一个系数都要靠一套独立的实验去测。

    但真正的威力还不在省事,而在「预言」。第 3 份(8.4 电动效应)你会看到一个具体案例:那里有四个看似独立的不可逆效应,用上昂萨格关系之后,只剩两个是独立的——另外两个可以从前两个算出来。

    换句话说,你在实验室里测了效应 A,就不必再做实验,直接可以预言效应 B 的数值。课本称这是「不可逆过程热力学得到的典型结果」。
    深入它的根在微观,不在宏观
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干里的 (8.2.4)
    $$L_{ij} = L_{ji}$$

    先说清楚它的地位:在本章的框架里,昂萨格关系是一条外来的输入,不是从 (8.2.1) 推出来的。课本明说了:「Onsager 关系要在统计物理中才能证明。」

    它的根子是微观可逆性,也叫时间反演对称性:把牛顿方程(或薛定谔方程)里的 $t$ 换成 $-t$,方程形式不变。粒子运动的录像倒着放,看上去同样合乎力学规律。

    昂萨格 1931 年的论证,桥梁是涨落回归假设:

    • 平衡态并不是死水一潭,各处的密度、温度时刻在做微小的自发涨落;
    • 涨落出现之后会自己衰减回去,而它衰减的规律,和宏观上外界强加一个同样大小的偏离后的弛豫规律,是同一条;
    • 微观可逆性要求这些涨落的关联函数对称;
    • 这份对称性传递到宏观,就成了 $L_{ij}=L_{ji}$。

    这项工作让昂萨格拿了 1968 年的诺贝尔化学奖——一条式子拿一个诺奖,可见分量。

    顺便解释 (8.2.5) 里那个奇怪的 $-\boldsymbol B$:磁场下洛伦兹力 $q\boldsymbol v\times\boldsymbol B$ 里含着速度 $\boldsymbol v$,时间倒放时 $\boldsymbol v \to -\boldsymbol v$,力就反向了——时间反演对称被磁场破坏了。要修复它,必须同时把 $\boldsymbol B$ 也反向,两个负号抵消。所以对称关系变成了 $L_{ij}(\boldsymbol B) = L_{ji}(-\boldsymbol B)$。

    这里有一条对本章更要紧的附加条件,第 11 段就要用:课本紧接着说——「必须选择适当的力和流才满足此关系。」

    也就是说,$L_{ij}=L_{ji}$ 不是对任意一种写法都成立的。力和流选错了,这条漂亮的关系就会失效。这是下一段的全部内容。
    实例数一数要做几次实验
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干的 (8.2.2)(8.2.3) 那组方程
    $$\begin{aligned}J_1 &= L_{11}X_1 + L_{12}X_2\\[2pt] J_2 &= L_{21}X_1 + L_{22}X_2\end{aligned}$$

    假设你要完整刻画某种溶液的热扩散行为。没有昂萨格关系时,四个系数得各测各的:

    要测怎么测
    $L_{11}$令 $X_2=0$(浓度均匀),加温差,测热流
    $L_{21}$令 $X_2=0$,加温差,测物质流 ← 难,量很小
    $L_{22}$令 $X_1=0$(温度均匀),造浓度差,测物质流
    $L_{12}$令 $X_1=0$,造浓度差,测热流 ← 最难,要在有扩散的同时精确测微小热量

    四次实验里,两个交叉项那两次又难做又不准——交叉效应本来就弱,信号淹在噪声里。

    有了 $L_{12}=L_{21}$ 之后:挑四次里最容易做的那一次去测交叉系数,另一个直接抄过来。三次实验就够了,而且可以避开最难的那一次。

    还有一个额外的好处,实验物理里非常实用:如果你把 $L_{12}$ 和 $L_{21}$ 都测了,理论说它们必须相等——那么两次测量的结果对不对得上,就成了一次对整套实验装置的独立校验。对不上,说明装置有系统误差。一条理论关系,反过来成了实验的标尺。
    第四步 · 收尾 · 并交给第 2 份一个问题

    11. 一个陷阱:$\nabla T$ 还是 $\nabla\!\left(\frac{1}{T}\right)$?

    0:00 / 0:00

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

    框架搭完了,但有一颗地雷必须在这里引爆,否则第 2 份会看不懂。

    回到第 8 段末尾挑明的那个隐患:热力学力的定义很宽松,同一个物理情形可以写出好几套不同的「力」。热传导的力,写成 $\nabla T$ 可以,写成 $\nabla\!\left(\dfrac{1}{T}\right)$ 也可以——两者都是「温度不均匀」的度量,都在温度均匀时为零。

    看起来随便选一个都行。但课本在这里给了一句很重的话:

    「上面我们写出的两个方程 (8.2.2) 和 (8.2.3) 中的力和流的形式,它们的系数不满足 Onsager 关系。」

    也就是说——我们上一段刚刚欣赏过的 $L_{12}=L_{21}$,对第 9 段写下的那两个方程并不成立。力选错了,昂萨格关系就垮了。

    这是个相当难受的处境:一套框架的核心结论,居然依赖于「你把力写成什么样」这种看起来很随意的事。

    那到底该按什么标准选?

    课本给出的答案是一个明确的判据:

    正确的力和流的形式,要从计算「熵产生率」得到。

    具体做法是:由熵密度产生率的表达式

    $$\frac{\mathrm{d}_i S}{\mathrm{d}t} = \sum_{j=1}^{k} J_j X_j$$
    (8.2.6)

    来定义与流 $J_j$ 对应的力 $X_j$。在这样得到的 $J$ 和 $X$ 之间,线性关系的系数才满足昂萨格关系。

    先认一下 $\mathrm{d}_i S$ 这个下标 $i$,它不是求和指标。系统的熵变有两个来源,非平衡热力学必须把它们分开:
     $\mathrm{d}_e S$ —— 从外界流进来的熵(下标 $e$ 取自 external)。可正可负。
     $\mathrm{d}_i S$ —— 系统内部自己产生的熵(下标 $i$ 取自 internal)。永远 $\ge 0$,这就是热力学第二定律。

    总熵变是两者之和:$\mathrm{d}S = \mathrm{d}_e S + \mathrm{d}_i S$。所以「熵产生率」$\dfrac{\mathrm{d}_i S}{\mathrm{d}t}$ 指的只是内部产生的那一份,不含流进流出。不可逆性有多强,全写在这个量里——它为零,过程就是可逆的。

    换句话说,力和流不是可以随手指定的,它们必须成对地从同一个式子里读出来:先老老实实算出这个过程的熵产生率,把它整理成「若干项相加、每项都是两个因子相乘」的形式,然后——

    $$\frac{\mathrm{d}_i S}{\mathrm{d}t} \;=\; \underbrace{J_1}_{\text{流}}\cdot\underbrace{X_1}_{\text{力}} \;+\; \underbrace{J_2}_{\text{流}}\cdot\underbrace{X_2}_{\text{力}} \;+\;\cdots$$

    ——每一项里的两个因子,就是一对配对的流和力。这就是共轭配对:力学里力与位移配对相乘得到功,热力学里 $T$ 与 $\mathrm{d}S$、$p$ 与 $\mathrm{d}V$ 配对相乘得到能量,而这里 $X$ 与 $J$ 配对相乘得到的是熵产生率。

    易错问题不在那个因子,在于「只换了力」
    0:00 / 0:00

    补充内容独立播放;关闭补充框会自动暂停。

    起点 · 主干里的两个候选
    $$X = \nabla T \qquad\text{还是}\qquad X = \nabla\!\Big(\frac{1}{T}\Big)\ ?$$

    先把两者的关系算出来。用链式法则求 $1/T$ 的梯度:

    $$\nabla\!\Big(\frac{1}{T}\Big) = \frac{\mathrm{d}}{\mathrm{d}T}\Big(\frac{1}{T}\Big)\cdot\nabla T = -\frac{1}{T^2}\,\nabla T$$

    错误的理解:「两者只差一个因子 $-1/T^2$,把它吸收进系数 $L$ 里不就完了。」

    为什么错——但也不是因为「这个因子随位置变」。真正的判据只有一条:力和流相乘,必须还原出熵产生率。你可以随便换写法,但必须成对地换:只要同时改流和力、让乘积保持不变(哪怕改动的因子随位置变),系数矩阵做的是 $L' = A^{\mathsf{T}}LA$ 这一类变换,对称性原样保留。

    出事的是「只换力、不换流」。把力从 $\nabla\!\left(\frac{1}{T}\right)$ 换成 $\nabla T$,而流仍旧用物理上的热流 $\boldsymbol{J}_q$——这一对相乘就不再等于熵产生率,它根本不是一对共轭的流与力,昂萨格关系对它不作任何保证。

    具体算一下。从正确的一套出发(系数 $L$,满足 $L_{12}=L_{21}$),只把第 1 个力换掉、流不动,由 $\nabla\!\left(\frac{1}{T}\right) = -\frac{1}{T^2}\nabla T$ 得

    $$L'_{12} = L_{12},\qquad L'_{21} = -\frac{L_{21}}{T^{2}}$$

    换掉第一个力,动的是系数矩阵的第一列($L_{11}$ 与 $L_{21}$);而 $L_{12}$ 是「力 2 对流 1 的响应」,跟力 1 怎么写毫无关系,纹丝不动。原本相等的两个,现在差了一个 $-1/T^2$。

    一句话记住这个陷阱:昂萨格关系说的是「配对正确的那一套系数是对称的」。而配对正确的唯一标准,就是两者相乘等于熵产生率。拆散了配对再去换写法,对称性当然就没了。所以必须有一个标准,昂萨格关系才有确定的含义。

    (8.2.6) 就是那个标准。它不是四个候选里比较好看的一个,它是唯一一个能让昂萨格关系成立的选法。

    这一份到此为止,问题交给下一份

    于是本份的旅程停在一个明确的未完成状态——框架已经搭好,但它有一个空缺没填:

    已经有了还差什么
    力 $X$ 与流 $J$ 的概念(第 8 段)具体到某个现象,$X$ 和 $J$ 各自到底是什么
    线性唯象方程 $J_j=\sum_i L_{ji}X_i$(第 9 段)方程里该代入哪一套 $X$
    昂萨格关系 $L_{ij}=L_{ji}$(第 10 段)它成立的前提——正确的力和流
    判据 $\dfrac{\mathrm{d}_iS}{\mathrm{d}t}=\sum_j J_jX_j$(本段)熵产生率本身怎么算

    最后那一行,就是第 2 份(8.3 节)的全部内容:老老实实地、一个情形一个情形地把熵产生率算出来,然后从结果里读出真正的力和流。

    先剧透一个结果,你会看到它有多反直觉:算完之后会发现,产生热流的力并不是以 $\nabla T$ 的形式出现的,而是

    $$\boldsymbol{X}_q = \nabla\!\Big(\frac{1}{T}\Big)$$

    ——正是这一段一直在纠结的那个「另一个候选」。它赢了。而且赢得有理有据:它是从熵产生率的表达式里自己掉出来的,不是谁挑的。

    本份小结

    整份的逻辑链条是:热力学不管快慢,因为它的前提是「无限慢」→ 松开这个前提,先看四条经验定律(傅里叶、欧姆、牛顿黏滞、菲克)→ 四条骨架相同:流 $=-$ 系数 $\times$ 梯度 → 要在非平衡态里谈 $T$、$\mu$,先靠局部平衡假设拿到资格 → 把梯度统称为力 $X$、把流统称为流 $J$ → 写出线性唯象方程 $J_j=\sum_i L_{ji}X_i$,其中交叉项 $L_{12}$、$L_{21}$ 是全新的内容(Dufour / Soret 效应)→ 昂萨格关系 $L_{ij}=L_{ji}$ 把独立系数减半,还能预言新效应 → 但它只对正确选取的力和流成立,而正确的选法必须由熵产生率 $\dfrac{\mathrm{d}_iS}{\mathrm{d}t}=\sum_j J_jX_j$ 来定。

    下一份(8.3 节)就去把这个熵产生率算出来。