反應器中的混合對反應的影響

反應器中的混合對反應的影響第1頁，課件共96頁，創(chuàng)作于2023年2月概述返混定義理想置換反應器（PFR）理想混合反應器（CSTR）兩種流體流動極限模型反應器物料在反應器中的混合在空間上的混合在時間上的混合這種混合過程返混簡單混合返混第2頁，課件共96頁，創(chuàng)作于2023年2月理想反應器的流動模式----平推流和全混流。間歇釜全混釜u=const平推流

理想的平推流和間歇釜停留時間均一，無返混。全混釜反應器的返混無窮大，出口物料停留時間分布與釜內物料的停留時間相同。第3頁，課件共96頁，創(chuàng)作于2023年2月第4頁，課件共96頁，創(chuàng)作于2023年2月實際反應器流動形式的復雜性u溝流回流存在速度分布存在死區(qū)和短路現(xiàn)象存在溝流和回流

偏離理想流動模式，反應結果與理想反應器的計算值具有較大的差異。DeadzoneShortcircuiting第5頁，課件共96頁，創(chuàng)作于2023年2月

實際流動與PF的偏離(軸向返混)

1.并行流存在速度分布2.存在分子熱運動DistributionH2空氣TT0熱反饋QEddy渦流DispersionVLaminar第6頁，課件共96頁，創(chuàng)作于2023年2月

實際流動與CSTR的偏離短路活塞流滯流離子反應

非CSTR胃曝氣池

生化反應CSTR第7頁，課件共96頁，創(chuàng)作于2023年2月影響反應結果的三大要素：停留時間分布（residencetimedistribution,RTD）凝集態(tài)（stateofaggregation）早混或遲混（earlinessandlatenessofmixing）RTD對反應的影響注入檢測第8頁，課件共96頁，創(chuàng)作于2023年2月平均停留時間為實際停留時間ti不盡相同，轉化率x1,x2,…,x5亦不相同。出口轉化率應為5個質點轉化率的平均值，即InjectionDetection第9頁，課件共96頁，創(chuàng)作于2023年2月聚集態(tài)的影響理想反應器假定混合為分子尺度，實際工程難以達到，如結團彌散兩種體系的反應程度顯然應該是不同的。工程中，盡量改善體系的分散尺度，以達到最有效的混合，從而改善反應效果。鼓泡氣體液體噴霧第10頁，課件共96頁，創(chuàng)作于2023年2月混合遲早度的影響早混晚混即使兩反應體系的空時相同，由于反應混合的遲早不同，反應結果也不相同。第11頁，課件共96頁，創(chuàng)作于2023年2月第一節(jié)連續(xù)反應器中物料混合狀態(tài)分析一、混合現(xiàn)象的分類按混合對象的年齡可以把混合分成兩種：（1）同齡混合：相同年齡物料之間的混合

（2）返混（backmixing）：不同年齡物料之間的混合

混合尺度的大小混合也可分為兩種類型：

宏觀(macro-)混合：設備尺度上的混合

微觀(micro-)混合：物料微團尺度上的混合第12頁，課件共96頁，創(chuàng)作于2023年2月混合的機理總體流動：攪拌器旋轉時使釜內液體產(chǎn)生一定途徑的循環(huán)流動

設備尺度上的宏觀均勻

高速旋轉的旋渦與液體微團產(chǎn)生相對運動和剪切力

更小尺度上的均勻

分子擴散

微團最終消失

微觀均勻均相第13頁，課件共96頁，創(chuàng)作于2023年2月提高混合效果的措施：

消除打旋現(xiàn)象：螺旋漿式渦輪式第14頁，課件共96頁，創(chuàng)作于2023年2月提高混合效果的措施：加設導流筒螺旋式渦輪式第15頁，課件共96頁，創(chuàng)作于2023年2月攪拌器的型式高轉速攪拌器

1、螺旋槳式攪拌器

螺旋槳式攪拌器的總體循環(huán)流動推進式三葉片式第16頁，課件共96頁，創(chuàng)作于2023年2月2、渦輪式攪拌器渦輪式攪拌器的總體循環(huán)流動a-直葉圓盤渦輪b-彎葉圓盤渦輪c-直葉渦輪d-折葉渦輪e-彎葉渦輪第17頁，課件共96頁，創(chuàng)作于2023年2月大葉片低轉速攪拌器

1、槳式攪拌器

2、框式和錨式攪拌器

3、螺帶式攪拌器

a-錨式bc-框式d-螺帶式第18頁，課件共96頁，創(chuàng)作于2023年2月二、連續(xù)反應過程的考察方法在連續(xù)攪拌反應釜或管式反應器中進行反應，如果反應物料的微觀混合程度不同，則考察方法即研究方法就不同。微觀混合有兩種極限狀態(tài)，完全混合和完全不混合，研究方法完全不同。

第19頁，課件共96頁，創(chuàng)作于2023年2月1、以反應容積（或微元）為對象

當物料微觀混合為完全混合時，物料呈分子狀均勻分散，物料不存在微團。對于攪拌反應器，物料以反應器為邊界，對于管式反應器，物料以dVR為邊界，所以研究的對象分別為反應器容積VR和反應器微元容積dVR。第20頁，課件共96頁，創(chuàng)作于2023年2月2、以反應物料為對象

當物料微觀混合為完全不混合時，物料呈微團獨立運動，物料的邊界為微團的邊界，所以以微團為研究對象，結合物料的停留時間分布函數(shù)和動力學方程方程可以有定量結果。

第21頁，課件共96頁，創(chuàng)作于2023年2月

如果微觀混合處于中間狀態(tài)，則幾個微團可以組成微元，研究對象為微元，目前只有定性的認識，沒有定量結果。綜上所述，考察對象都是物料，不同的是按照微觀混合的程度劃分考察的范圍：完全混合——反應容積VR或dVR；中間狀態(tài)——微元（由微團組成）；完全不混合——微團第22頁，課件共96頁，創(chuàng)作于2023年2月第二節(jié)停留時間分布停留時間分布的數(shù)學描述停留時間分布的實驗測定幾種流型的停留時間分布函數(shù)與分布密度停留時間分布的應用第23頁，課件共96頁，創(chuàng)作于2023年2月一、流體在反應器內的停留時間分布返混程度停留時間分布規(guī)律定量確定停留時間概率函數(shù)概率密度函數(shù)數(shù)學期望方差兩個特征值兩個函數(shù)描述第24頁，課件共96頁，創(chuàng)作于2023年2月1、停留時間分布的數(shù)學描述全混流反應器：機械混合最大逆向混合最大返混程度無窮大平推流反應器：機械混合為零逆向混合為零返混程度等于零間歇反應器：機械混全最大逆向混合為零返混程度等于零反應器內的返混程度不同—停留時間不同—濃度分布不同—反應速率不同—反應結果不同—生產(chǎn)能力不同非理想流動反應器：介于兩種理想情況之間停留時間是隨機變量，因此停留時間分布是一種概率分布。第25頁，課件共96頁，創(chuàng)作于2023年2月1）、停留時間分布的定量描述借用人口統(tǒng)計學（Population）兩個統(tǒng)計參數(shù)a)社會人口的年齡分布和b)死亡年齡分布，在反應工程中假設：

在反應器出口流體微元的停留時間分布：E(t),F(t)反應器注入檢測原料出料a)b)各微元保持獨立身份(identification),即微元間不能混合c)不研究微元在反應器內的歷程,只研究它在反應器內的停留時間,即壽命。則定義：第26頁，課件共96頁，創(chuàng)作于2023年2月停留時間（壽命）的概念？？例：在連續(xù)操作的反應器內，如果在某一瞬間（t=0）極快地向入口物流中加入100個紅色粒子，同時在系統(tǒng)的出口處記下不同時間間隔流出的紅色粒子數(shù)，結果如下表。如果假定紅色粒子和主流體之間除了顏色的差別以外，其余所有性質都完全相同，那么就可以認為這100個粒子的停留時間分布就是主流體的停留時間分布。停留時間范圍t→t+△t0-22-33-44-55-66-77-88-99-1010-1111-1212-14出口流中的紅色粒子數(shù)02612182217126410分率△N/N00.020.060.120.180.220.170.120.060.040.010第27頁，課件共96頁，創(chuàng)作于2023年2月以時間t為橫坐標，出口流中紅色粒子數(shù)為縱坐標，將上表作圖：若以停留時間t為橫坐標，為縱坐標作圖，則每一個長方形的面積為即表示停留時間為t→t+△t的物料占總進料的分率。第28頁，課件共96頁，創(chuàng)作于2023年2月假如示蹤劑改用紅色流體，連續(xù)檢測出口中紅色流體的濃度，如果將觀測的時間間隔縮到非常小，得到的將是一條連續(xù)的停留時間分布曲線。E(t)tt+dt

t圖中曲線下微小面積E(t)dt表示停留時間在t和t+dt之間的物料占t=0時進料的分率。第29頁，課件共96頁，創(chuàng)作于2023年2月(1)分布密度函數(shù)E(t)E(t)定義為在t=0時刻進入反應器的流體微元，在t時離開反應器的概率根據(jù)定義，E(t)應具有歸一性，即因為當時間無限長時，t=0時刻加入的流體質點都會流出反應器第30頁，課件共96頁，創(chuàng)作于2023年2月(2)分布函數(shù)F(t)F(t)函數(shù)定義為在t時刻離開反應器流體質點數(shù)占總示蹤流體質點數(shù)的分率E(t)與F(t)的關系：

第31頁，課件共96頁，創(chuàng)作于2023年2月

二、停留時間分布規(guī)律的實驗測定入口處加“激勵”出口處研究“響應”示蹤法停留時間分布規(guī)律確定研究階躍示蹤法脈沖示蹤法從原物料到示蹤劑的階躍變化用示蹤劑代替原物料的瞬間變化第32頁，課件共96頁，創(chuàng)作于2023年2月停留時間分布的測定一般采用示蹤技術，示蹤劑選用易檢測其濃度的物質，根據(jù)其光學、電學、化學及放射等特性，采用比色、電導、放射檢測等測定濃度。選擇示蹤劑要求：

1)與主流體物性相近；2)高低濃度均易檢測，以減少示蹤劑的用量；3)不產(chǎn)生相變或相轉移；4)示蹤劑濃度易轉變?yōu)楣庑盘柣螂娦盘?，以便于計算機數(shù)據(jù)采集和處理。停留時間的測定方法根據(jù)示蹤劑的加入方式分為脈沖法、階躍法和周期輸入法，前兩者應用較廣。第33頁，課件共96頁，創(chuàng)作于2023年2月1、階躍示蹤法從某一時刻起，將原物料全部切換為示蹤物，示蹤物濃度階躍突變。01.0(a)“激勵”曲線01.0(b)“響應”曲線圖4-3階躍法測定停留時間分布函數(shù)第34頁，課件共96頁，創(chuàng)作于2023年2月由圖可知，在t=0時，C=0；t，CC0

時間為t時，出口物料中示蹤劑濃度為C(t)，物料流量為V，所以示蹤劑流出量為VC(t)，又因為在時間為t時流出的示蹤劑，也就是反應器中停留時間小于t的示蹤劑，按定義，物料中停留時間小于t的粒子所占的分率為F(t)，因此，當示蹤劑入口流量為VC0時，出口流量VC0F(t)，所以有：因此，用此法可直接方便地測定實際反應器的留時間分布函數(shù)。第35頁，課件共96頁，創(chuàng)作于2023年2月2、脈沖示蹤法在盡可能短的時間內，用示蹤物瞬間代替不含示蹤物的物料，然后又立刻恢復為原來的物料。即給物料一個脈沖訊號，與之同時測定出口流的響應曲線圖4-4脈沖法測定停留時間分布密度函數(shù)01.0(b)“響應”曲線0(a)“激勵”曲線第36頁，課件共96頁，創(chuàng)作于2023年2月設Δt0時間內注入示蹤劑的總量為M（mol），出口處濃度隨時間變化為C(t)，在示蹤劑注入后tt+dt時間間隔內，出口處流出的示蹤劑量占總示蹤劑量的分率：若在注入示蹤劑的同時，流入反應器的物料量為N，在注入示蹤劑后的t

t+dt時間間隔內，流出物料量為dN，則在此時間間隔內，流出的物料占進料的分率為：第37頁，課件共96頁，創(chuàng)作于2023年2月示蹤劑的停留時間分布就是物料質點的停留時間分布，即：因此：有：只要測得V，M和C（t），即可得物料質點的分布密度。第38頁，課件共96頁，創(chuàng)作于2023年2月由于M=VC0Δt0，C0及Δt0難以準確測量，故示蹤劑的總量可用出口所有物料的加和表示：因此，利用脈沖法可以很方便的測出停留時間分布密度。第39頁，課件共96頁，創(chuàng)作于2023年2月

脈沖法和階躍法的比較脈沖法階躍法示蹤劑注入方法在原有的流股中加入示蹤劑，不改變原流股流量將原有流股換成流量與其相同的示蹤劑流股E(t)可直接測得F(t)可直接測得第40頁，課件共96頁，創(chuàng)作于2023年2月研究不同流型的停留時間分布，通常是比較它們的統(tǒng)計特征值。常用的特征值有兩個：數(shù)學期望—平均值方差—離散程度平均停留時間它是指整個物料在設備內的停留時間，而不是個別質點的停留時間。不管設備型式和個別質點的停留時間，只要反應體積與物料體積流量比值相同，平均停留時間就相同。三、停留時間分布的數(shù)字特征第41頁，課件共96頁，創(chuàng)作于2023年2月E(t)dt=dF(t)F(t)：所有停留時間為0-t的質點所占的分率F(t+dt)：所有停留時間為0-t+dt的質點所占的分率dF(t)=F(t+dt)-F(t)dF(t)：所有停留時間為t-t+dt的質點所占的分率數(shù)學期望：所有質點停留時間的“加權平均值”（4-17）,（4-18）第42頁，課件共96頁，創(chuàng)作于2023年2月在等時間間隔取樣時：tt1t2t3…….E(t)E(t1)E(t2)E(t3)……對于離散型測定值，可以用加和代替積分值（4-19）第43頁，課件共96頁，創(chuàng)作于2023年2月方差是停留時間分布離散程度的量度方差越小，越接近平推流對平推流，各物料質點的停留時間相等，故方差為零。方差：各個物料質點停留時間t與平均停時間差的平方的加權平均值。（4-20）第44頁，課件共96頁，創(chuàng)作于2023年2月如果是離散型數(shù)據(jù)，將積分改為加和：取樣為等時間間隔時：（4-20）第45頁，課件共96頁，創(chuàng)作于2023年2月平均對比時間：停留時間為t時，，因此，θ和t一一對應，且有：F（θ）=F（t），此時：歸一性：對比時間（無因次時間）：（4-22）第46頁，課件共96頁，創(chuàng)作于2023年2月用θ表示的方差：（4-24）第47頁，課件共96頁，創(chuàng)作于2023年2月描述停留時間分布的兩個函數(shù)：（4-24）第48頁，課件共96頁，創(chuàng)作于2023年2月四、理想流型反應器停留時間分布平推流模型全混流模型第49頁，課件共96頁，創(chuàng)作于2023年2月1、平推流：“激勵”曲線“響應”曲線圖4-5理想置換反應器“激勵”與“響應”曲線兩條曲線完全一樣，只是“響應”曲線比“激勵”曲線平移了一段時間。第50頁，課件共96頁，創(chuàng)作于2023年2月所有物料質點的停留時間都相同，且等于整個物料的平均停留時間tm，停留時間分布函數(shù)與分布密度為：由方差定義，第51頁，課件共96頁，創(chuàng)作于2023年2月設進行階躍注入實驗，反應器的容積為VR，物料的體積流量為V，達到穩(wěn)態(tài)后，從t=0開始，將進料切換為含示蹤劑濃度為C0的物料，在切換后某dt時間內，對全釜作物料衡算：進入的示蹤劑量=流出的示蹤劑量+示蹤劑的積累量全混流第52頁，課件共96頁，創(chuàng)作于2023年2月給全混流反應器一個階越“激勵”“激勵”曲線“響應”曲線

全混流反應器“激勵”與“響應”曲線第53頁，課件共96頁，創(chuàng)作于2023年2月對反應器做示蹤物的物料衡算：第54頁，課件共96頁，創(chuàng)作于2023年2月例：現(xiàn)有一個停留時間相同的全混釜和平推流反應器當二者以不同的方式進行串聯(lián)時，其各自的最終轉化率為多少？已知：停留時間為1min，反應速率常數(shù)為1.0m3/kmol.min,液相反應物濃度為1kmol/m解：全混釜在前，平推流在后全混釜停留1min后反應物的濃度：第55頁，課件共96頁，創(chuàng)作于2023年2月反應物濃度為cA1后進入平推流在平推流停留1min后反應物的濃度：平推流在前停留1min后反應物的濃度：第56頁，課件共96頁，創(chuàng)作于2023年2月在全混釜中再停留1min后反應物的濃度：第57頁，課件共96頁，創(chuàng)作于2023年2月

五、利用RTD診斷反應器內流動狀況第58頁，課件共96頁，創(chuàng)作于2023年2月第59頁，課件共96頁，創(chuàng)作于2023年2月實際反應器中的流動狀況總是偏離理想流動很難建立其真實方程可以先建立一種非理想流動模型，用它來描述實際反應器中的流動情況再通過對模型參數(shù)估值來確定偏離理想流動的具體程度常用模型有：層流模型；軸向擴散模型；多級串聯(lián)全混流模型；組合模型非理想流動是指物料在反應器內有部分返混發(fā)生，目前對非理想流動一般借助于模型加以描述。第三節(jié)非理想流動模型第60頁，課件共96頁，創(chuàng)作于2023年2月

CSTR(0維)PFR(一維)PD模型(一維)徑向擴散PFR(準二維)CoCC0CC0CC0CC1徑向射流(圓柱形三維)四角射流(異形三維)渦流V(x,y,z)C(x,y,z)T(x,y,z)圓柱座標直角座標1、反應器流動模式第61頁，課件共96頁，創(chuàng)作于2023年2月2、停留時間分布的應用確定模型參數(shù)m或Pe用多級串聯(lián)全混流模型或軸向擴散模型模擬實際反應器中的流動狀況，關鍵是確定串聯(lián)的級數(shù)m或Pe，m或Pe又與方差有關，因此，可以通過實驗確定停留時間分布，進而計算方差，m或Pe，然后求得轉化率。定性分析流動狀況活塞流全混流定量分析流動狀況實際反應器中可能存在短路與死角，使實際的平均停留時間不等于VR/V，因此可以得用停留時間分布來定量估算死角與短路的程度。第62頁，課件共96頁，創(chuàng)作于2023年2月例：某全混流反應器體積為100L，物料流率為1L/s，試求在反應器中停留時間為（1）90～110s，（2）0～100s，（3）>100s的物料占總進料的比率。解：出口物料的份額用F(t)表示，所求比率90～110s：F(110)－F(90)=0.074=7.4%小于平均停留時間的物料占63.2%大于的物料占總物料的36.8%第63頁，課件共96頁，創(chuàng)作于2023年2月

實際反應器內的流動是很復雜的，僅用理想化的平推流或全混流進行計算是不夠的，非常必要對實際的流型進行逼近模擬。常用的方法有如下幾種。第64頁，課件共96頁，創(chuàng)作于2023年2月填料塔層流模型可對平推流模型進行修正，但比較粗糙。主要模型假設為：1)反應混合物在管內作層流流動。2)反應混合物在軸徑向均不存在返混。u第65頁，課件共96頁，創(chuàng)作于2023年2月軸向混合模型在平推流模型的基礎上再迭加一個軸向混合軸向擴散模型的基本假設：第66頁，課件共96頁，創(chuàng)作于2023年2月

流體沿軸向有參數(shù)變化，徑向參數(shù)均一；流體流動主流為平推流但疊加一逆向渦流擴散；逆向渦流擴散遵循菲克定理，但整個反應器內擴散系數(shù)為一常數(shù)。第67頁，課件共96頁，創(chuàng)作于2023年2月軸向擴散模型仿照分子擴散中用擴散系數(shù)分子運動表征有效擴散系數(shù)表征一維返混用平推流反應器+渦流擴散運動非理想流動的管式反應器主要用于湍流流動的管式反應器、固定床反應器、塔式反應器。第68頁，課件共96頁，創(chuàng)作于2023年2月軸向混合模型對實際反應器，處理時在平推流的基礎上迭加一個軸向混合來進行校正。適合于不存在死角、短路和循環(huán)流、返混程度較小的非理想流動模型。模型參數(shù)是軸向混合彌散系數(shù)EZ，停留時間分布可表示為EZ的函數(shù)。Peclet準數(shù)：當返混程度相當小時；數(shù)學期望θ=1方差對一級不可逆反應，轉化率可表示為：第69頁，課件共96頁，創(chuàng)作于2023年2月第70頁，課件共96頁，創(chuàng)作于2023年2月00.51.01.52.00.51.01.52.00.050.10.5110θ第71頁，課件共96頁，創(chuàng)作于2023年2月00.51.01.52.00.20.40.60.81.00.0010.0120.0250.25θ第72頁，課件共96頁，創(chuàng)作于2023年2月基本假設：由m個體積相等的CSTR串聯(lián)組成；從一個CSTR到下一個CSTR之間的管道內物料不發(fā)生反應?！嗉壢炝鞔?lián)模型多釜串聯(lián)模型:用幾個等體積的全混流反應器串聯(lián)來模擬實際反應器中的流動狀況。第73頁，課件共96頁，創(chuàng)作于2023年2月00.51.01.50.51.0m=1m=2m=5m=10m→∞00.51.01.50.51.0m=1m=2m=5m=10m→∞圖4-11不同m值時的E(θ),F(θ)曲線第74頁，課件共96頁，創(chuàng)作于2023年2月每一級的停留時間ti=tm/m。模型參數(shù)為串聯(lián)級數(shù)m。方差m=1時，即為全混流模型m=∞時，即為平推流模型對一級不可逆反應，轉化率可表示為：假設實際反應器中的返混程度與m個等體積的全混流反應器串聯(lián)時相同，m是虛擬釜數(shù)，不一定是整數(shù)。第75頁，課件共96頁，創(chuàng)作于2023年2月

當m=1時，多級混合釜模型的停留時間分布函數(shù)與CSTR相同。當m?∞時，多級混合釜模型的停留時間分布函數(shù)與PFR相同。當1<m<∞時，多級混合釜模型的停留時間分布函數(shù)屬于非理想流動反應器。第76頁，課件共96頁，創(chuàng)作于2023年2月多級混合釜模型的應用1、多個CSTR串聯(lián)操作時，當釜數(shù)由一增加到無限多個，其停留時間分布規(guī)律將按CSTR到非理想流動反應器到PFR的分布規(guī)律變化。2、如果釜數(shù)選擇恰當，該模型的停留時間分布規(guī)律可以與任何一個非理想流動反應器的停留時間分布規(guī)律相同（即兩者的返混程度相同）。3、可以用多級全混流模型計算任何一種非理想流動反應器的轉化率。反應釜個數(shù)為：第77頁，課件共96頁，創(chuàng)作于2023年2月組合模型適用于上述兩種模型不能很好表達的情況將實際反應器的流動情況設想為平推流、全混流、死區(qū)、短路、循環(huán)流等部分組成組合模型的幾種典型例子：第78頁，課件共96頁，創(chuàng)作于2023年2月第五節(jié)非理想流動反應器的計算準確地計算實際反應器的生產(chǎn)能力和轉化率，需已知兩個前提條件：

1)本征動力學方程

2)反應器停留時間分布函數(shù)如實際反應器中諸微元具有獨立身份，則可以想象為實際反應器由不同長度管式反應器并聯(lián)組成：入口出口為各并聯(lián)反應器轉化率的積分平均。第79頁，課件共96頁，創(chuàng)作于2023年2月出口轉化率＝Σ停留時間為ti的轉化率×ti的流量分量即或無論何種反應器，只要已知停留時間分布函數(shù)，即可接上式計算出口第80頁，課件共96頁，創(chuàng)作于2023年2月對全混流對2級動力學第81頁，課件共96頁，創(chuàng)作于2023年2月對平推流對2級動力學第82頁，課件共96頁，創(chuàng)作于2023年2月例：應用脈沖示蹤法測定一容積為12L的反應裝置，進入反應器的流體速度0.8L/min，在定常態(tài)下脈沖地輸入80g示蹤劑A，并同時在反應器出口處記錄cA隨時間的變化，其實測數(shù)據(jù)列于下表:t(min)05101520253035cA(g/L)03554210若用其進行某一級液相不可逆反應(-rA)=kcA,k=0.307min-1,當平均停留時間為15min時，試求：（1）采用多級全混流模型時，其最終轉化率為多少？（2）采用軸向返混模型時，其最終轉化率為多少？第83頁，課件共96頁，創(chuàng)作于2023年2月t(min)05