传质微分方程及扩散传质

传质微分方程及扩散传质第一页，共二十八页，编辑于2023年，星期五2§3.1传质微分方程§3.1.1传质微分方程推导：A通过微分体积表面x处表面的量，在处故在方向上A的净流出速度：

A在微分单元内的累积速度：

设为单位体积单位时间内由于化学反应产生的A量，则A产生速度为：。第二页，共二十八页，编辑于2023年，星期五质量守恒原理[体系内的积累]=[通过体系边界的净流入量]

+[体系内的净生成量]根据物理化学基本原理推导出来的所有模型，都是以上式表示的原理为基础。上式对质量、动量、能量都适用。3第三页，共二十八页，编辑于2023年，星期五4根据质量守恒原理，有：两边同时除以，则有：

即同样可推导在三维方向上：简写成称为传质微分方程。同理，对于以摩尔通量表示的形式为：第四页，共二十八页，编辑于2023年，星期五

§3.1.2菲克第二定律

扩散介质中小体积单元如图所示截面1进入体积单元的扩散流截面2流出的扩散流第五页，共二十八页，编辑于2023年，星期五6假设：扩散过程无化学反应，那么扩散进入体积单元的量减去流出体积单元的量等于体积单元内物质的积累量。以代表体积单元的截面积，则（3–1）因为扩散流随变化，故

（3–2）将（3–2）代入（3–1），得（3–3）根据菲克第一定律代入（3–3）得（3–4）第六页，共二十八页，编辑于2023年，星期五7

当D为常数，即不随扩散距离、浓度变化时，有（3–5）即为菲克第二定律。（3–5）可简写为当组分向三维空间扩散时，则有（直角坐标系）（圆柱坐标系）（球坐标体系）求解三维扩散方程非常复杂，所以一般在制定实验方案时，近似地安排成一维扩散，在特定边界条件下解（3–5）式一元二阶微分方程。第七页，共二十八页，编辑于2023年，星期五8

菲克第二定律成立的条件

①无扩散引起的对流传质；②扩散体系内无化学反应；③为常数。适用于固体或静止流体中的扩散。稳态扩散

当达到稳态时，故，即第八页，共二十八页，编辑于2023年，星期五9菲克第二定律的应用菲克第二定律：已知条件：1）初始条件（时间上的已知值）：

t=0时，

t=0时，

t=∞时，

2）边界条件（某空间上的已知值）：①规定某界面的浓度；②规定某界面的化学反应速度；③将对流传质作为边界条件。第九页，共二十八页，编辑于2023年，星期五10例1

在1273K时，用混合气体对低碳钢（）进行渗碳，设钢板内部扩散为过程的控制步骤。钢板表面碳平衡浓度。求渗碳6小时后钢铁表面下处的碳浓度。

已知：。解：这是固体内部的扩散过程，故适用菲克第二定律，取扩散方向为轴。

起始条件：扩散开始前，体系内浓度完全均匀，而为C0(初始浓度)。边界条件：①扩散开始后，界面的浓度立即为Cs，并且在扩散过程中保持不变，即②当，，。将上述已知条件代入，解方程得：第十页，共二十八页，编辑于2023年，星期五11例1

在1273K时，用混合气体对低碳钢（）进行渗碳，设钢板内部扩散为过程的控制步骤。钢板表面碳平衡浓度。求渗碳6小时后钢铁表面下处的碳浓度。

已知：。解：这是固体内部的扩散过程，故适用菲克第二定律，取扩散方向为轴。

起始条件：扩散开始前，体系内浓度完全均匀，而为C0(初始浓度)。边界条件：①扩散开始后，界面的浓度立即为Cs，并且在扩散过程中保持不变，即②当，，。将上述已知条件代入，解方程得：第十一页，共二十八页，编辑于2023年，星期五12

erf(x)称为高斯误差函数，1-erf(x)称为补余误差函数。误差函数的性质：

erf(-x)=-erf(x)erf(0)=0erf(1)=11-erf(x)=erfc(x)erfc()=0,erfc(0)=1erfc(x)为补余误差函数。用作纵轴，为横轴作图

第十二页，共二十八页，编辑于2023年，星期五13

由，得知欲求时间t的扩散浓度，须先求出的值，再由图得值。查图得，。渗碳6小时，钢铁表面处的碳浓度

可以由这样的方法求固体中的扩散系数。

第十三页，共二十八页，编辑于2023年，星期五14§3.2伴有非均相反应的扩散整个过程包括扩散：反应物向界面扩散化学反应：在界面上发生化学反应，此时扩散仍遵守菲克第二定律。化学反应提供重要边界条件。例2

设有一球状的碳粒与氧反应生成一氧化碳试计算碳的燃烧速度。假定气膜层中的氧和一氧化碳不发生反应。第十四页，共二十八页，编辑于2023年，星期五15解：计算碳的燃烧速度，计算氧气消耗速度即可。取球形坐标，在扩散过程中没有化学反应发生，在稳态条件下,则传质微分方程（取球坐标）：简化成：

（3-89）

或（3-90）

边界条件：，，R为碳粒半径；，，xO2为空气中的氧分压。有上述两个边界条件，可以求出。第十五页，共二十八页，编辑于2023年，星期五16

而我们要求的是,因为根据菲克第一定律

对于碳和氧化学反应生成一氧化碳，有

（每扩散进来1摩尔氧气，就有2摩尔一氧化碳在相反方向上扩散出去）（3-95）

（3-98）在整个扩散进程中WO2保持为常数（这一点从3-89式可以得出，或者，因为是串联过程，扩散过程中通过不同球面的氧气质量是相等的）即：对（3-98）整理得：第十六页，共二十八页，编辑于2023年，星期五17两边分别求积分：为积分常数。当时，，故;

当时,，故(3–99)

若知道氧的量即可换算成碳的燃烧速度。1）假定化学反应速度比扩散快得多，则,2）假定化学反应速度和扩散速度相差不大，则非均相反应中化学反应速度能提供一个重要的边界条件：即时,

假设为一级基元反应，第十七页，共二十八页，编辑于2023年，星期五18将代入(3–99)，得当大时，上式中对数可按泰勒级数展开，方程简化为：

(3–104)关于泰勒级数:上式叫做函数f(x)在点x0处的泰勒级数。该式叫做函数f(x)在点x0处的泰勒展开式。第十八页，共二十八页，编辑于2023年，星期五19§3.3伴有均相化学反应的扩散传质微分方程

第十九页，共二十八页，编辑于2023年，星期五20两种类型化学反应;1传质相中均相反应2传质相中的非均相反应。反应速度作为边界条件（因为控制体积内无化学反应发生）本节只讨论第一种情况。例3

过程：气体A被液体B吸收，同时在液体中发生一级基元反应而被消耗。试求液面上气体A的的摩尔通量。分析：液面，，第二十页，共二十八页，编辑于2023年，星期五21假设B无限大，扩散达到稳态，A在液体中流动性很小，液体中流动很小，则有：由传质微分方程又故有：第二十一页，共二十八页，编辑于2023年，星期五22边界条件：，，解出：双曲正弦函数：双曲余弦函数：故，液面上的摩尔通量（现对于静止坐标）

双曲正切函数第二十二页，共二十八页，编辑于2023年，星期五23表示了化学反应的影响，它是一个无因次量，称为哈特数（Hattanumber）。当很大时，,方程简化为：

(3–115)

上述传质又称为液膜传质，与传质方程式比较：

液膜传质系数与扩散系数的1/2次方成正比。第二十三页，共二十八页，编辑于2023年，星期五24例：熔融金属在一窄管中的挥发，金属蒸汽在熔体上方和气相中的组分发生反应。如工业实践中粗硒的提纯。a控制A的挥发速度很小XA→0；

b用惰性气体稀释XB→0c窄管中平行流(非涡流)NAy+NBy→0

列方程：传质微分方程：A：

(1)(稳态)例4

窄管中金属A挥发并与气相中B反应——发生在管内

B：(2)

第二十四页，共二十八页，编辑于2023年，星期五25化学反应速度：设为不可逆二级反应(3)

菲克第一定律：(5)假