傳質(zhì)微分方程及擴散傳質(zhì)

傳質(zhì)微分方程及擴散傳質(zhì)第1頁，共28頁，2023年，2月20日，星期三2§3.1傳質(zhì)微分方程§3.1.1傳質(zhì)微分方程推導(dǎo)：A通過微分體積表面x處表面的量，在處故在方向上A的凈流出速度：

A在微分單元內(nèi)的累積速度：

設(shè)為單位體積單位時間內(nèi)由于化學(xué)反應(yīng)產(chǎn)生的A量，則A產(chǎn)生速度為：。第2頁，共28頁，2023年，2月20日，星期三質(zhì)量守恒原理[體系內(nèi)的積累]=[通過體系邊界的凈流入量]

+[體系內(nèi)的凈生成量]根據(jù)物理化學(xué)基本原理推導(dǎo)出來的所有模型，都是以上式表示的原理為基礎(chǔ)。上式對質(zhì)量、動量、能量都適用。3第3頁，共28頁，2023年，2月20日，星期三4根據(jù)質(zhì)量守恒原理，有：兩邊同時除以，則有：

即同樣可推導(dǎo)在三維方向上：簡寫成稱為傳質(zhì)微分方程。同理，對于以摩爾通量表示的形式為：第4頁，共28頁，2023年，2月20日，星期三

§3.1.2菲克第二定律

擴散介質(zhì)中小體積單元如圖所示截面1進入體積單元的擴散流截面2流出的擴散流第5頁，共28頁，2023年，2月20日，星期三6假設(shè)：擴散過程無化學(xué)反應(yīng)，那么擴散進入體積單元的量減去流出體積單元的量等于體積單元內(nèi)物質(zhì)的積累量。以代表體積單元的截面積，則（3–1）因為擴散流隨變化，故

（3–2）將（3–2）代入（3–1），得（3–3）根據(jù)菲克第一定律代入（3–3）得（3–4）第6頁，共28頁，2023年，2月20日，星期三7

當(dāng)D為常數(shù)，即不隨擴散距離、濃度變化時，有（3–5）即為菲克第二定律。（3–5）可簡寫為當(dāng)組分向三維空間擴散時，則有（直角坐標(biāo)系）（圓柱坐標(biāo)系）（球坐標(biāo)體系）求解三維擴散方程非常復(fù)雜，所以一般在制定實驗方案時，近似地安排成一維擴散，在特定邊界條件下解（3–5）式一元二階微分方程。第7頁，共28頁，2023年，2月20日，星期三8

菲克第二定律成立的條件

①無擴散引起的對流傳質(zhì)；②擴散體系內(nèi)無化學(xué)反應(yīng)；③為常數(shù)。適用于固體或靜止流體中的擴散。穩(wěn)態(tài)擴散

當(dāng)達到穩(wěn)態(tài)時，故，即第8頁，共28頁，2023年，2月20日，星期三9菲克第二定律的應(yīng)用菲克第二定律：已知條件：1）初始條件（時間上的已知值）：

t=0時，

t=0時，

t=∞時，

2）邊界條件（某空間上的已知值）：①規(guī)定某界面的濃度；②規(guī)定某界面的化學(xué)反應(yīng)速度；③將對流傳質(zhì)作為邊界條件。第9頁，共28頁，2023年，2月20日，星期三10例1

在1273K時，用混合氣體對低碳鋼（）進行滲碳，設(shè)鋼板內(nèi)部擴散為過程的控制步驟。鋼板表面碳平衡濃度。求滲碳6小時后鋼鐵表面下處的碳濃度。

已知：。解：這是固體內(nèi)部的擴散過程，故適用菲克第二定律，取擴散方向為軸。

起始條件：擴散開始前，體系內(nèi)濃度完全均勻，而為C0(初始濃度)。邊界條件：①擴散開始后，界面的濃度立即為Cs，并且在擴散過程中保持不變，即②當(dāng)，，。將上述已知條件代入，解方程得：第10頁，共28頁，2023年，2月20日，星期三11例1

在1273K時，用混合氣體對低碳鋼（）進行滲碳，設(shè)鋼板內(nèi)部擴散為過程的控制步驟。鋼板表面碳平衡濃度。求滲碳6小時后鋼鐵表面下處的碳濃度。

已知：。解：這是固體內(nèi)部的擴散過程，故適用菲克第二定律，取擴散方向為軸。

起始條件：擴散開始前，體系內(nèi)濃度完全均勻，而為C0(初始濃度)。邊界條件：①擴散開始后，界面的濃度立即為Cs，并且在擴散過程中保持不變，即②當(dāng)，，。將上述已知條件代入，解方程得：第11頁，共28頁，2023年，2月20日，星期三12

erf(x)稱為高斯誤差函數(shù)，1-erf(x)稱為補余誤差函數(shù)。誤差函數(shù)的性質(zhì)：

erf(-x)=-erf(x)erf(0)=0erf(1)=11-erf(x)=erfc(x)erfc()=0,erfc(0)=1erfc(x)為補余誤差函數(shù)。用作縱軸，為橫軸作圖

第12頁，共28頁，2023年，2月20日，星期三13

由，得知欲求時間t的擴散濃度，須先求出的值，再由圖得值。查圖得，。滲碳6小時，鋼鐵表面處的碳濃度

可以由這樣的方法求固體中的擴散系數(shù)。

第13頁，共28頁，2023年，2月20日，星期三14§3.2伴有非均相反應(yīng)的擴散整個過程包括擴散：反應(yīng)物向界面擴散化學(xué)反應(yīng)：在界面上發(fā)生化學(xué)反應(yīng)，此時擴散仍遵守菲克第二定律?；瘜W(xué)反應(yīng)提供重要邊界條件。例2

設(shè)有一球狀的碳粒與氧反應(yīng)生成一氧化碳試計算碳的燃燒速度。假定氣膜層中的氧和一氧化碳不發(fā)生反應(yīng)。第14頁，共28頁，2023年，2月20日，星期三15解：計算碳的燃燒速度，計算氧氣消耗速度即可。取球形坐標(biāo)，在擴散過程中沒有化學(xué)反應(yīng)發(fā)生，在穩(wěn)態(tài)條件下,則傳質(zhì)微分方程（取球坐標(biāo)）：簡化成：

（3-89）

或（3-90）

邊界條件：，，R為碳粒半徑；，，xO2為空氣中的氧分壓。有上述兩個邊界條件，可以求出。第15頁，共28頁，2023年，2月20日，星期三16

而我們要求的是,因為根據(jù)菲克第一定律

對于碳和氧化學(xué)反應(yīng)生成一氧化碳，有

（每擴散進來1摩爾氧氣，就有2摩爾一氧化碳在相反方向上擴散出去）（3-95）

（3-98）在整個擴散進程中WO2保持為常數(shù)（這一點從3-89式可以得出，或者，因為是串聯(lián)過程，擴散過程中通過不同球面的氧氣質(zhì)量是相等的）即：對（3-98）整理得：第16頁，共28頁，2023年，2月20日，星期三17兩邊分別求積分：為積分常數(shù)。當(dāng)時，，故;

當(dāng)時,，故(3–99)

若知道氧的量即可換算成碳的燃燒速度。1）假定化學(xué)反應(yīng)速度比擴散快得多，則,2）假定化學(xué)反應(yīng)速度和擴散速度相差不大，則非均相反應(yīng)中化學(xué)反應(yīng)速度能提供一個重要的邊界條件：即時,

假設(shè)為一級基元反應(yīng)，第17頁，共28頁，2023年，2月20日，星期三18將代入(3–99)，得當(dāng)大時，上式中對數(shù)可按泰勒級數(shù)展開，方程簡化為：

(3–104)關(guān)于泰勒級數(shù):上式叫做函數(shù)f(x)在點x0處的泰勒級數(shù)。該式叫做函數(shù)f(x)在點x0處的泰勒展開式。第18頁，共28頁，2023年，2月20日，星期三19§3.3伴有均相化學(xué)反應(yīng)的擴散傳質(zhì)微分方程

第19頁，共28頁，2023年，2月20日，星期三20兩種類型化學(xué)反應(yīng);1傳質(zhì)相中均相反應(yīng)2傳質(zhì)相中的非均相反應(yīng)。反應(yīng)速度作為邊界條件（因為控制體積內(nèi)無化學(xué)反應(yīng)發(fā)生）本節(jié)只討論第一種情況。例3

過程：氣體A被液體B吸收，同時在液體中發(fā)生一級基元反應(yīng)而被消耗。試求液面上氣體A的的摩爾通量。分析：液面，，第20頁，共28頁，2023年，2月20日，星期三21假設(shè)B無限大，擴散達到穩(wěn)態(tài)，A在液體中流動性很小，液體中流動很小，則有：由傳質(zhì)微分方程又故有：第21頁，共28頁，2023年，2月20日，星期三22邊界條件：，，解出：雙曲正弦函數(shù)：雙曲余弦函數(shù)：故，液面上的摩爾通量（現(xiàn)對于靜止坐標(biāo)）

雙曲正切函數(shù)第22頁，共28頁，2023年，2月20日，星期三23表示了化學(xué)反應(yīng)的影響，它是一個無因次量，稱為哈特數(shù)（Hattanumber）。當(dāng)很大時，,方程簡化為：

(3–115)

上述傳質(zhì)又稱為液膜傳質(zhì)，與傳質(zhì)方程式比較：

液膜傳質(zhì)系數(shù)與擴散系數(shù)的1/2次方成正比。第23頁，共28頁，2023年，2月20日，星期三24例：熔融金屬在一窄管中的揮發(fā)，金屬蒸汽在熔體上方和氣相中的組分發(fā)生反應(yīng)。如工業(yè)實踐中粗硒的提純。a控制A的揮發(fā)速度很小XA→0；

b用惰性氣體稀釋XB→0c窄管中平行流(非渦流)NAy+NBy→0

列方程：傳質(zhì)微分方程：A：

(1)(穩(wěn)態(tài))例4

窄管中金屬A揮發(fā)并與氣相中B反應(yīng)——發(fā)生在管內(nèi)

B：(2)

第24頁，共28頁，2023年，2月20日，星期三25化學(xué)反應(yīng)速度：設(shè)為不可逆二級反應(yīng)(3)

菲克第一定律：(5)