《热学热力学与统计物理》课程教学资源（教材讲义，上）第八章 非平衡热力学（输运现象）非平衡态相变

第八章非平衡热力学（输运现象）非平衡态相变$8.1输运现象的经验规律我们前面处理的热力学体系均处于热力学平衡态，所经历的过程均是可逆过程（即准静态过程），对不可递过程仅涉及用热力学函数来判断过程进行的方向：我们也计算了一些特殊情况下的不可逆过程的变。但并未涉及不可逆过程的核心问题，即过程进行的速率。本章主要介绍近平衡的物理现象（输运现象）。所谓近平衡是指离平衡态不远，“力”和“流”成线性关系，如由于温度梯度（称“力”）引起的热流（称“流”）之间的关系是线性的，即温度梯度比较小的情况。当温度梯度比较大时，称远平衡，“力”和“流”呈非线性关系。本章对远平衡的现象（非平衡态相变）仅作简单介绍。下面先介绍近平衡时的几个经验规律。当在物体中存在温度梯度时，就会产生一个热流（能量的输运)。如一根金属棒，长度为1，两端温度分别为T和T,，T>T（图8.1)，则Q7图8.1，热传导单位时间内通过截面积为A的热量为：O=KAL-T(8.1)10用j。代表单位时间内流过单位面积的热量，称热流密度，则上式可写成矢量形式：AJ。 =-x(T)VT,(8.2)在一维情况下可写成：dlJgx = -x(T)(8.3)dx上两式称傅里叶定律，x为热导率。若在导体两端有一个电势梯度√β，则在导体中产生一个电流（电量的输运），电流密度为:J.=oE=-OVp,(8.4)

第八章 非平衡热力学（输运现象） 非平衡态相变 §8.1 输运现象的经验规律 我们前面处理的热力学体系均处于热力学平衡态，所经历的过程均是可逆过程（即准静 态过程），对不可逆过程仅涉及用热力学函数来判断过程进行的方向；我们也计算了一些特 殊情况下的不可逆过程的熵变。但并未涉及不可逆过程的核心问题，即过程进行的速率。本 章主要介绍近平衡的物理现象(输运现象)。所谓近平衡是指离平衡态不远，“力”和“流”成 线性关系，如由于温度梯度（称“力”）引起的热流（称“流”）之间的关系是线性的，即温 度梯度比较小的情况。当温度梯度比较大时，称远平衡，“力”和“流”呈非线性关系。本 章对远平衡的现象(非平衡态相变)仅作简单介绍。 下面先介绍近平衡时的几个经验规律。当在物体中存在温度梯度时，就会产生一个热流 (能量的输运)。如一根金属棒，长度为 l ，两端温度分别为 T1 和 T2 ，T2 > T1 （图 8.1），则 图 8.1，热传导 单位时间内通过截面积为 A 的热量为： l T T Q A 2 − 1  =  ， (8.1) 用 A Q J q   = 代表单位时间内流过单位面积的热量，称热流密度，则上式可写成矢量形式： J q = −(T )T  ， (8.2) 在一维情况下可写成： ( ) dx dT J qx = − T (8.3) 上两式称傅里叶定律，  为热导率。 若在导体两端有一个电势梯度  ，则在导体中产生一个电流(电量的输运)，电流密度 为： Je =E = −   ， (8.4)

此即欧姆定律，为电导率。dvy如果在一个流体中，流体的速度在x方向，并在方向存在一个速度梯度则两层dz流体间有相互作用，上层流体对下层流体施加的力，大小由牛顿定律给出（图8.2）：dvAS,f=n(8.5)dz>Vxdvxdzz1图8.2黏滞流速度分布其中n为黏滞系数，△S为作用的面积。流体的黏滞运动是动量的传输，因此上式可写成：△Kdyf.△tPax=(8.6)=ndzASAtASAt其中P.表示动量流密度，即单位时间内在单位面积上传输的动量，动量的方向在x方向，动量传输的方向在2方向。若在混合物中存在浓度n的梯度，将引起物质流，此现象称扩散现象（物质的输运），遵守斐克（Fick）定律：dndM=-dSdt,(dz)20j,=dM=-DVn,(8.7)dSdt式中J，为粒子流密度（混合物中某组元物质在单位时间内流过单位面积的粒子数），Vn为浓度梯度，D为扩散系数。在化学反应中，化学势梯度也产生物质流。非平衡态热力学，也称不可逆过程热力学是把以上的经验规律提高到热力学的高度找出一般规律。$8.2基本假设为了讨论近平衡时的非平衡态现象，我们还要以平衡态热力学为基础，再引进局部平衡

此即欧姆定律，  为电导率。 如果在一个流体中，流体的速度在 x 方向，并在 z 方向存在一个速度梯度 dz dvx ，则两层 流体间有相互作用，上层流体对下层流体施加的力，大小由牛顿定律给出（图 8.2）： S dz dv f x =  ， (8.5) 图 8.2 黏滞流速度分布 其中  为黏滞系数， S 为作用的面积。流体的黏滞运动是动量的传输，因此上式可写成： dz dv S t f t S t K P x zx =    =    = ， (8.6) 其中 Pzx 表示动量流密度，即单位时间内在单位面积上传输的动量，动量的方向在 x 方向， 动量传输的方向在 z 方向。 若在混合物中存在浓度 n 的梯度，将引起物质流，此现象称扩散现象(物质的输运)，遵守 斐克（Fick）定律： dSdt dz dn dM D z0       = − ， D n dSdt dM J n = = −   ， (8.7) 式中 n J  为粒子流密度（混合物中某组元物质在单位时间内流过单位面积的粒子数），n 为 浓度梯度， D 为扩散系数。在化学反应中，化学势梯度也产生物质流。 非平衡态热力学，也称不可逆过程热力学是把以上的经验规律提高到热力学的高度找出 一般规律。 §8.2 基本假设 为了讨论近平衡时的非平衡态现象，我们还要以平衡态热力学为基础，再引进局部平衡

假设，即每一个小块宏观小、微观大的区域，仍可近似看成处于热力学平衡态，对此小区域仍可用热力学参量来描述，并可定义热力学函数。在此局部平衡假设下，对每一小区域，吉布斯方程仍然成立，即：Tds = dU + pdv-Zμ,dN,i但整个体系是非平衡态。所以对此体系用上了两个矛盾的概念，整个体系是非平衡态，而小区域又是平衡态，这只有在小区域与周围的能量和物质交换均是很小的情况下才成立，即变化过程要很缓慢，仍可看成是准静态过程。在整个体系中，温度T、压强p、N，、μ,及热力学函数均是x、J、二、t的函数。如整个体系的熵可写成（注意摘的可加性）：S=[s(x,y,z, tldxdydz。当整个体系处于非平衡态下，必然出现温度梯度、速度梯度、电势梯度、浓度梯度等等，相应在体系中会产生热流、动量流、电流、物质流等，因此在非平衡态热力学中要引进两个新的概念，即力-X和流-J，X是指上面的温度梯度、速度梯度、电势梯度、浓度梯度等，J指热流、动量流、电流、物质流等。这样做和力学吻合，但比力学上的概念要广泛得多。所以不可逆过程必与力和流联系。在近平衡情况下，力和流是线性的。假如有一组力X，，引起的流为：ZLIX,.J,=(8.8)i=l用一个例子加以说明，如果k=2，上式可写成：J =LX1+L12X2J = L2iX, + L22X2,如果两个力为：X，=VT（温度梯度），X，=Vn（浓度梯度），则上式写为：(8.9)J,=L,VT+L2Vn，J,为热流(8.10)J,=L2,VT+L2Vn，J2为物质流每一项的物理意义为：L,VT一由于温度梯度产生的热流（傅里叶定律）：L2Vn一由于浓度梯度产生的物质流，物质流携带内能从而引起热流；L2VT一由于温度梯度引起的热扩散，从而引起的物质流；L22Vn一由于浓度梯度产生的物质流（斐克定律）。浓度梯度引起的扩散，会在体系中产生温度差。此效应称Dufour效应（1872年发现）

假设，即每一个小块宏观小、微观大的区域，仍可近似看成处于热力学平衡态，对此小区域 仍可用热力学参量来描述，并可定义热力学函数。 在此局部平衡假设下，对每一小区域，吉布斯方程仍然成立，即： = + − i TdS dU pdV idNi ， 但整个体系是非平衡态。所以对此体系用上了两个矛盾的概念，整个体系是非平衡态，而小 区域又是平衡态，这只有在小区域与周围的能量和物质交换均是很小的情况下才成立，即变 化过程要很缓慢，仍可看成是准静态过程。在整个体系中，温度 T 、压强 p 、Ni 、 i 及热 力学函数均是 x、y、z、t 的函数。如整个体系的熵可写成（注意熵的可加性）： S s(x y z t)dxdydz  = ， ， ， 。 当整个体系处于非平衡态下，必然出现温度梯度、速度梯度、电势梯度、浓度梯度等等， 相应在体系中会产生热流、动量流、电流、物质流等，因此在非平衡态热力学中要引进两个 新的概念，即力- X 和流- J ，X 是指上面的温度梯度、速度梯度、电势梯度、浓度梯度等， J 指热流、动量流、电流、物质流等。这样做和力学吻合，但比力学上的概念要广泛得多。 所以不可逆过程必与力和流联系。 在近平衡情况下，力和流是线性的。假如有一组力 Xi ，引起的流为： = = k i j Lji Xi J 1 。 (8.8) 用一个例子加以说明，如果 k = 2 ，上式可写成： 1 L11X1 L12X2 J = + 2 L21X1 L22X2 J = + ， 如果两个力为： X1 = T （温度梯度）， X2 = n （浓度梯度），则上式写为： J1 = L11T + L12n， 1 J 为热流 （8.9） J 2 = L21T + L22n， 2 J 为物质流 （8.10） 每一项的物理意义为： L11T —由于温度梯度产生的热流（傅里叶定律）； L12n—由于浓度梯度产生的物质流，物质流携带内能从而引起热流； L21T —由于温度梯度引起的热扩散，从而引起的物质流； L22n—由于浓度梯度产生的物质流（斐克定律）。 浓度梯度引起的扩散，会在体系中产生温度差。此效应称 Dufour 效应（1872 年发现）

热扩散现象，又叫Soret效应（1893年发现）。这样从热力学角度推广了一些定律。在J,和J,的系数之间存在一个关系，称翁萨格（Onsager）倒易关系，即：L21 = L/2'写成一般形式：L, = Lj(8.11)在有磁场情况下为：L,(B)= L,(-B)。(8.12)Onsager关系要在统计物理中才能证明，而且必须选择适当的力和流才满足此关系。例如(1温度梯度可以写成VT或√而上面我们写出的两个方程（8.9）和（8.10）中的力(动)和流的形式，它们的系数不满足Onsager关系。正确的力和流的形式，要从计算产生率而得到，即由密度产生率的表达式：d,s_=ZJ,X,(8.13)dt定义与流J，对应的力X，，在这样得到的J和X之间，线性关系的系数才满足Onsager关系。下面我们计算摘产生率。d,s摘密度产生率和翁萨格关系$8.3dt（1）热传导。先计算由于存在温度梯度√T，在物体中产生热流的情况。先讨论简单的情况，如果有两个小块物体I和II，中间由导热板隔开，整个体系是封闭体系，即只能与外界交换能量，而不能交换物质。每一小块可看成是处于平衡态，两小块中的多数参量均相同，仅是温度不同，它们的温度分别为T和T"，令T!<T（见图8.3)。deg'导热板aee"/XVIII-dg"dig'图8.3，两物体之间的热传导整个体系的摘为：

热扩散现象，又叫 Soret 效应（1893 年发现）。这样从热力学角度推广了一些定律。 在 1 J 和 2 J 的系数之间存在一个关系，称翁萨格（Onsager）倒易关系，即： L21 = L12， 写成一般形式： Lij = Lji ， (8.11) 在有磁场情况下为： L (B) L ( B) ij ji   = − 。 (8.12) Onsager 关系要在统计物理中才能证明，而且必须选择适当的力和流才满足此关系。例如 温度梯度可以写成 T 或        T 1 ，而上面我们写出的两个方程（8.9）和（8.10）中的力 和流的形式，它们的系数不满足 Onsager 关系。正确的力和流的形式，要从计算熵产生率而 得到，即由熵密度产生率的表达式： i k j i i J X dt d s  = = 1 (8.13) 定义与流 i J 对应的力 Xi ，在这样得到的 J 和 X 之间，线性关系的系数才满足 Onsager 关 系。下面我们计算熵产生率。 §8.3 熵密度产生率 dt d si 和翁萨格关系 （1） 热传导。 先计算由于存在温度梯度 T ，在物体中产生热流的情况。先讨论简单的情况，如果有 两个小块物体 I 和 II，中间由导热板隔开，整个体系是封闭体系，即只能与外界交换能量， 而不能交换物质。每一小块可看成是处于平衡态，两小块中的多数参量均相同，仅是温度不 同，它们的温度分别为 I T 和 II T ，令 I T ＜ II T （见图 8.3）。 图 8.3，两物体之间的热传导 整个体系的熵为：

S=s'+s",每一小块的TdS方程为T'ds'= dU' +pdV -μdN,由于体积和粒子数不变，故T'ds' = du' =deo' +do',T"ds"= du"=deo" +dio",其中deQ'和deQ"分别为从外界流入小块I和小块ⅡI的热量，而diQ'是由于存在温度梯度从小块Ⅱ传入I的热量，diQ"是从小块I传入IⅡI的热量，且diQ'=-diQ"，即I吸收的热量等于ⅡI放出的热量。所以1deg'+_1dg'+1ds=ds'+ds"Tideo"+TirdgTTI等号后面的前两项为摘流，后两项为熵产生。可得摘产生率为d,s(1dio'ldt(8.14)dt上式中A()-(-)为力，而4Q'ht为流，从这里可看到，产生热流的力不是以VT的形式出现，而是V以上是以简单化的形式导出的，下面我们把它推广到热力学量连续变化的情况。考虑一根棒的两端分别与两个热源接触，温度分别为T和T，（T,>T），则在棒中就有热流通过。为计算摘产生，必须把棒分成微观大、宏观小的小块。我们引进体积元d’r=dxdydz，其中心在（x，y，z）。令u为内能密度（单位质量物质的内能），p为物质的密度，则体积元d"r中的内能为：pud"r。因热传导过程中，体积和粒子数不变，故体积元中的内能变化仅是热量的流入或流进，热量流进体积元，则内能增加，反之内能减少。所以单位体积中内能的减少为：_ a(pu) = divj,at整个体积内能的变化为以下积分： ,-,(8.15)

I II S = S + S ， 每一小块的 TdS 方程为： T dS dU pdV dN I I I = + −  ， 由于体积和粒子数不变，故 = = I I I T dS dU đe I Q +đi I Q ， = = II II II T dS dU đe II Q +đi II Q ， 其中 đe I Q 和 đe II Q 分别为从外界流入小块 I 和小块 II 的热量，而 đi I Q 是由于存在温度梯 度从小块 II 传入 I 的热量， đi II Q 是从小块 I 传入 II 的热量，且 đi I Q =-đi II Q ，即 I 吸收 的热量等于 II 放出的热量。所以 I I II T dS dS dS 1 = + = đe I Q + II T 1 đe II Q + I T 1 đi I Q + II T 1 đi II Q 等号后面的前两项为熵流，后两项为熵产生。可得熵产生率为：       = −I II i dt T T d S 1 1 đi I Q /dt (8.14) 上式中        = −       I II T T T 1 1 1 为力，而 đi I Q /dt 为流，从这里可看到，产生热流的力不是以 T 的形式出现，而是        T 1 。 以上是以简单化的形式导出的，下面我们把它推广到热力学量连续变化的情况。考虑一 根棒的两端分别与两个热源接触，温度分别为 T1 和 T2 （ T2 > T1 ），则在棒中就有热流通过。 为计算熵产生，必须把棒分成微观大、宏观小的小块。我们引进体积元 d r = dxdydz 3 ，其 中心在 （x，y，z）。令 u 为内能密度（单位质量物质的内能），  为物质的密度，则体积元 d r 3 中的内能为： ud r 3  。因热传导过程中，体积和粒子数不变，故体积元中的内能变化仅 是热量的流入或流进，热量流进体积元，则内能增加，反之内能减少。所以单位体积中内能 的减少为： ( ) q divJ t u  =   −  ， 整个体积内能的变化为以下积分： ( ) q V A q V u dxdydz divJ dxdydz dA J t    = =    −     ， (8.15)

j，为热流密度，jdA为在单位时间内从面积元dA中流出的热量。定义s(x，y，z，t)为单位质量的，称熵密度，则体积元中的熵为：p.s(x, y, z, t).dxdydz,把Tds=du代入，从方程(8.15)可得：Pa-ldnj,--PatTT整个体积为：(--a+/,dxdydz,(8.16)T上式的右边第一项为焰流，第二项为摘产生率，所以熵密度产生率是：ds-j.(8.17)dt（2）扩散考虑两个相等体积的容器和I，中间有一带小孔的隔板，它们的温度一样，但两边物质的摩尔数不一样，均有k个组元，若体积不变，则有：Tds = dU' + dU" -Z(u, dN, + μ" dN")ds=ao-(-"!,7T上式中aO=dU=dU'+dU"，dN!=-dN"。式中第一项为熵流，第二项为焰产生，所以焰产生率为：aNZ=-AadtdtTdtTMM因为ds(μ)S>0，从上式可见，若->0，即容器II中的化学势高于容器ITdt(T)T中的化学势，则>0，容器I中的粒子数增加，即物质流总是从化学势高的地方流向化dt学势低的地方，这和导体中的电势一样。推广到连续的情况，摘密度产生率可写成：ds-jVμ(8.18)dtT

q J  为热流密度， J q dA    为在单位时间内从面积元 dA  中流出的热量。 定义 s(x，y，z，t) 为单位质量的熵，称熵密度，则体积元中的熵为：  s(x，y，z，t) dxdydz， 把 Tds = du 代入，从方程(8.15)可得：       = =   −    − T J T J divJ t T s q q q 1  1    ， 整个体积为： ( )          = −  +    V q q V A dxdydz T J T J s dxdydz dA t  1    , (8.16) 上式的右边第一项为熵流，第二项为熵产生率，所以熵密度产生率是：       =  T J dt d s q i  1 。 (8.17) （2） 扩散 考虑两个相等体积的容器 I 和 II，中间有一带小孔的隔板，它们的温度一样，但两边物 质的摩尔数不一样，均有 k 个组元，若体积不变，则有： = + −( + ) i II i II i I i I i I II TdS dU dU  dN  dN T dS 1 = đQ ( ) I i i II i I i dN T −   −  1 ， 上式中 đQ= I II dU = dU + dU ， II i I dNi = −dN 。式中第一项为熵流，第二项为熵产生，所 以熵产生率为： dt dN dt T dN dt T T d S I i i i I i II i i I i i              = − − =  −     （ ） ， 因为  0 dt diS ，从上式可见，若 T T T I i II  i  i   = −      −  >0，即容器 II 中的化学势高于容器 I 中的化学势，则  0 dt dN I i ，容器 I 中的粒子数增加，即物质流总是从化学势高的地方流向化 学势低的地方，这和导体中的电势一样。推广到连续的情况，熵密度产生率可写成：        =  − T J dt d s ni i   ， (8.18)

Vu这里J为物质流，为力T（3）有电势情况下的扩散下面考虑带电粒子的扩散，情况和上面类同，但现在的粒子是带电的。这样除了有化学势差，还有电势差，即β≠β"+0。如果令每摩尔第i个组元的带电量为q，，则有：d=ao-(-u")N, -("-",TT其中最后一项是由于有了电场以后，内能上要加上一项q,βdN（α=I，IⅡI)，在交换粒子时，这部分能量也要迁移过去。令，=μ+q,P，称电化学势，它的作用类似化学势，上式可写成：a!ds=TTi摘产生率为：dNu(8.19)dtTTdtdt(T上式实际上是用电化学势代替了化学势，带电粒子流也总是从电化学势高的地方流向电化学势低的地方。假如两边化学势相等，则有外电场存在时，粒子从电势高的地方向电势低的地方流，而且只要电势足够高，粒子可从一方全部流向另一方。（4）同时存在温度梯度和电势梯度的情况金属中有正离子和导电电子，假定正离子不动，则导电和导热均由自由电子承担。金属中的共有化电子群可看成自由电子气体。现在，在导体上加一个电场（电势梯度），并有温度梯度，讨论电子气体的输运现象。定义u、S、n为单位体积的内能、摘、摩尔数，则在体积元dr中的上述量分别为ud'r、sd'r、nd'r。由能量守恒定律可得单位体积的内能变化为：du= du, + duz'其中du,为由热流引起的内能变化，由前面导出的公式可得：ou, =-dinj.atdu，是由于粒子流流过体积元而引起的内能变化，粒子的化学势为μ，即每个粒子携带的能

这里 ni J  为物质流，        − T  i 为力。 （3） 有电势情况下的扩散 下面考虑带电粒子的扩散，情况和上面类同，但现在的粒子是带电的。这样除了有化学 势差，还有电势差，即   0 I II   。如果令每摩尔第 i 个组元的带电量为 i q ，则有： T dS 1 = đQ ( ) ( ) I i i II i I i I i i II i I i q q dN T dN T −   −  −   −  1 1 ， 其中最后一项是由于有了电场以后，内能上要加上一项 qi  dN（i  = I，II） ，在交换粒子 时，这部分能量也要迁移过去。令 i = i + qi ~ ，称电化学势，它的作用类似化学势，上 式可写成： T dS 1 = đQ ( ) I i i II i I i dN T −   −  1 ~ ~ 熵产生率为： dt dN dt T dN dt T T d S I i i i I i II i i I i i              = − − =  −     ~ ~ ~ （ ） (8.19) 上式实际上是用电化学势代替了化学势，带电粒子流也总是从电化学势高的地方流向电化 学势低的地方。假如两边化学势相等，则有外电场存在时，粒子从电势高的地方向电势低的 地方流，而且只要电势足够高，粒子可从一方全部流向另一方。 （4） 同时存在温度梯度和电势梯度的情况 金属中有正离子和导电电子，假定正离子不动，则导电和导热均由自由电子承担。金 属中的共有化电子群可看成自由电子气体。现在，在导体上加一个电场（电势梯度），并有 温度梯度，讨论电子气体的输运现象。 定义 u 、 s 、 n 为单位体积的内能、熵、摩尔数，则在体积元 d r 3 中的上述量分别为： ud r 3 、 sd r 3 、 nd r 3 。 由能量守恒定律可得单位体积的内能变化为： du = du1 + du2 ， 其中 1 du 为由热流引起的内能变化，由前面导出的公式可得： q divJ t u  = −   1 2 du 是由于粒子流流过体积元而引起的内能变化，粒子的化学势为  ，即每个粒子携带的能

量，粒子带电荷e，导体电位为时，其能量为e，故总能量为=μ+ep，带电粒子流流过体积元而引起的内能变化为：ouz =-divij, =-pidij,-J,gradjat其中j,代表粒子流密度矢量。上式用了矢量公式：div(FA=fdivA+A·gradf。所以在单位时间内单位体积的内能变化为：u --divj, -nidivj,-J.gradpat由吉布斯方程：Tds = du - jidn可得：as_1oujionatT atT aton+dij,=0。这样又粒子数守恒定律：atas=-divj,1j. gradpat-TT.-divlI+J.gradJngradpT(T)T上面后一式用了散度公式。式中第一项为摘流，第二项为温度梯度引起热流而致的摘产生，第三项是电化学势梯度引起的的摘产生。故熵密度产生率为d=j,·gradj.gradj(8.20)dt力（graa-gradp）和流（j，j,）有以下关系：)-Lvj,=L.V(8.21)-LaVj, = L21(8.22)J,=e.J, - Le或(eV+Vμ)(8.23)统计物理可以导出上面力和流的两个关系式，并证明翁萨格关系（L2=L）

量，粒子带电荷 e ，导体电位为  时，其能量为 e ，故总能量为  =  + e ~ ，带电粒 子流流过体积元而引起的内能变化为：    2 ~ ~ ~ div J divJ J grad t u n n n    = − = − −   其中 n J  代表粒子流密度矢量。上式用了矢量公式： div(fA) = fdivA + A gradf    。 所以在单位时间内单位体积的内能变化为：   ~ ~ divJ divJ J grad t u q n n    = − − −   由吉布斯方程： Tds du dn ~ = − 可得： t n t T u t T s   −   =    ~ 1 ， 又粒子数守恒定律： + = 0   n divJ t n  。这样  1 1 ~ J grad T divJ t T s q n = − −       1 1 1 ~ J grad T T J J grad T div q q n  −        +       = −     ， 上面后一式用了散度公式。式中第一项为熵流，第二项为温度梯度引起热流而致的熵产生， 第三项是电化学势梯度引起的的熵产生。故熵密度产生率为：  1 1 ~ J grad T T J grad dt d s q n i  −       =    (8.20) 力（       T grad 1 ，  1 ~ grad T − ）和流（ q J  ， n J  ）有以下关系：  1 1 ~ 11  − 12       =  T L T J q L  (8.21)  1 1 ~ 21  − 22       =  T L T J n L  (8.22) 或 ( ) 1 21  − 22  +       =  =   e T e L T J e J L e e n   （8.23） 统计物理可以导出上面力和流的两个关系式，并证明翁萨格关系（ L12 = L21 ）

$8.4电动效应与上节（3)中的情况一样，有两个容器和ⅡI，中间有一带小孔的隔板，两个容器的温度相同，但压强不等，电势不等，并设带电粒子为单组元，这样两个容器的条件为：T'=T",'+@"p"+p",Vo1由于化学势：u=μa(T，pa，na)α=I，Il。其中na=，而NaNa'va为摩尔数，为比体积，n为粒子密度。由于两个容器的压强不相等，而化学势又是压强的函数，所以实际上就是两个容器中的粒子的化学势不等，故有：μ≠μ"，另有电势差，所以此情况与上面（3）相同，可得熵产生率为：[u'+qp!μ"+qpdN!dN;dN'd,sAuqApTTTdtTdtdtdt由于偏离平衡态不远，所以化学势差可写成=V):op而Ap=pl-p,Au = vAp ,摘产生率可写成另一种形式：dN!dN;d,s--ApAp(8.24)qdtTdtTdtdN!dN!J.为电流，=J.为物质流，这样我们可以写出下列两个方程：9dtdtAQ+L2TApJ, = L(8.25)IApAp+LJ.=LaT(8.26)22T只有像上面那样写出的流与力的两个方程的系数之间才有翁萨格关系：L2 = L21(8.27)以上谓电动效应，有四个不可逆过程：L12(A=1，流致电势：(8.28)41(4p)j,=0此时J。=0，但△β≠0，是由于压力差产生的物质流而引起的电位差

§8.4 电动效应 与上节(3)中的情况一样，有两个容器 I 和 II，中间有一带小孔的隔板，两个容器的温度 相同，但压强不等，电势不等，并设带电粒子为单组元，这样两个容器的条件为： I II T = T ， I II p  p ， I II    ， 由于化学势：  =  (T，p ，n )， = I，II     。其中   v n 1 = ，而    N V v = ，  N 为摩尔数，  v 为比体积，  n 为粒子密度。由于两个容器的压强不相等，而化学势又是压 强的函数，所以实际上就是两个容器中的粒子的化学势不等，故有： I II    ，另有电势 差，所以此情况与上面（3）相同，可得熵产生率为： dt dN T q dt dN dt T dN T q T q dt d S I i I i I i I I II II i    −    = −      + − + = −       。 由于偏离平衡态不远，所以化学势差可写成（ v p T n =           ，  ）：  = vp ， 而 I II p = p − p ， 熵产生率可写成另一种形式： dt dN q dt T dN v T p dt d S I i I i i    −  = −  ， (8.24) e I i J dt dN − q = 为电流， m I i J dt dN − v = 为物质流，这样我们可以写出下列两个方程： T p L T Je L  +  = 11 12  (8.25) T p L T J m L  +  = 21 22  (8.26) 只有像上面那样写出的流与力的两个方程的系数之间才有翁萨格关系： L12 = L21 (8.27) 以上谓电动效应，有四个不可逆过程： 1，流致电势： 11 12 0 L L p e J = −           =  （8.28） 此时 Je = 0 ，但   0 ，是由于压力差产生的物质流而引起的电位差

L212，电渗效应：(8.29)LiAp=0此时压差为零，即Ap=0，由于单位电流引起的物质流。=23，流致电流：(8.30)L22(J.) g0此时电势差为零，即△の0，但有电流，是由单位体积流引起的电流。L21(Ap)=-4，电渗压：(8.31)(Ap)=0L22此时物质流为零，即J=0，压力差是由单位电势引起。如果无翁萨格关系，则此四种效应是独立的，应用翁萨格关系，则四个不可逆过程只有两个是独立的了。用(8.27)式，可得：AJe) (器)((8.32)-1Ap这是不可逆过程热力学得到的典型结果。$8.5热电效应I2）外，还有三种热电除了上面讨论过的热传导现象（遵守傅里叶定律）和焦耳热（一a效应。下面先分别介绍实验定律，然后用不可逆过程热力学讨论三者之间的关系。1，塞贝克（Seebeck）效应由两种不同的金属A和B连接成一个回路，如图8.4所示，两个接点分别处在不同的温度（与恒温块热接触，但电绝缘），中间接数字电压表，两个开路的接点必须保持在同一温度下，可测得电压△Vre，称温差电动势。此效应称塞贝克效应。AT+ATTBBDAVre图8.4，热电动势的测量回路

2，电渗效应： 11 21 0 L L J J e p m =          = (8.29) 此时压差为零，即 p = 0 ，由于单位电流引起的物质流。 3，流致电流： 22 12 0 L L J J m e =         = (8.30) 此时电势差为零，即  = 0 ，但有电流，是由单位体积流引起的电流。 4，电渗压： 22 21 0 L p L Jm = −           =  (8.31) 此时物质流为零，即 J m = 0 ，压力差是由单位电势引起。 如果无翁萨格关系，则此四种效应是独立的，应用翁萨格关系，则四个不可逆过程只有 两个是独立的了。用(8.27)式，可得： =0  =0         = −           e p m J J J p e  ， =0  =0         = −             m e J J p J m (8.32) 这是不可逆过程热力学得到的典型结果。 §8.5 热电效应 除了上面讨论过的热传导现象（遵守傅里叶定律）和焦耳热（ 1 2 I  ）外，还有三种热电 效应。下面先分别介绍实验定律，然后用不可逆过程热力学讨论三者之间的关系。 1， 塞贝克（Seebeck）效应 由两种不同的金属 A 和 B 连接成一个回路，如图 8.4 所示，两个接点分别处在不同的温 度（与恒温块热接触，但电绝缘），中间接数字电压表，两个开路的接点必须保持在同一温 度下，可测得电压 ΔVre，称温差电动势。此效应称塞贝克效应。 图 8.4，热电动势的测量回路