文丘里引射反應器的噴嘴優化及壓降分析

王曉靜，唐 磊，江 增

(天津大學化工學院，天津 300072)

文丘里引射反應器是利用噴射技術和射流紊動擴散作用研制出來的流體機械［1］。由于本身沒有運動部件，因而其主要優點是結構簡單，制造容易，工作可靠，安裝維護方便，密封性好，便于綜合利用，在許多工藝流程中如化學工業、石油開發等領域應用具有明顯的優越性和不可替代性，特別適合于水下和存在放射、易燃、易爆等可能性的危險場合。

黃磷是眾多含磷化合物的基本原料。通常，黃磷脫砷工藝利用機械攪拌器使黃磷與氧化性介質混合反應，進而達到凈化效果［2-3］。然而，使用機械攪拌存在很多不足之處，主要是黃磷易于結塊而導致脫砷反應不充分，同時需要耗費大量的人力、物力，伴隨著對攪拌器的檢修和維護工作。為了解決這個問題，研究以文丘里兩相射流混合技術，用雙氧水連續氧化脫除工業黃磷 (砷質量分數在400×10-6以上)中的砷從而制取低砷磷 (砷質量分數低于10×10-6)的一種全新工藝為背景，著重研究文丘里反應器的混合分散機理以及其流體力學行為，進一步研究其對黃磷脫砷反應過程的影響，以求得到適合這一反應過程的文丘里雙相射流反應器的優化結構參數和工藝設計與計算方法。目前已完成了小試可視化實驗裝置和中試裝置的制備，并進行了黃磷脫砷的研究和計算流體力學模擬計算，為該工藝的工業化放大提供了有利的參考依據。

目前，針對文丘里的尺寸優化［4-6］、流場模擬［7-10，15-16］和實驗研究［11-13］等成果較多，但通常是從理論分析上進行相關研究，而將其用于具體的化工過程---黃磷脫砷過程的研究并不多見，鑒于黃磷脫砷工藝在整個黃磷生產中扮演著重要的角色，因此有必要對黃磷脫砷工藝中的文丘里引射反應器的噴嘴出口直徑及相對應的壓降變化進行研究和分析，以獲得可供參考的數據和結論，指導實際應用。

在文丘里引射反應器的研究中，主要有理論計算、試驗研究和數值模擬3種方法。由于黃磷的燃點較低 (約40℃)，因此為了防止其自燃必須保證黃磷不接觸空氣。傳統以實驗為基礎的研究方法需要有相應的安全保護措施和大量的人力物力，研究周期較長。而數值模擬以其高效、低成本、能適應多種可變因素等優勢，在工程實踐中得到了廣泛的應用。本研究應用FLUENT 6.3軟件，分別計算了不同噴嘴直徑下反應器的性能及對應的壓降變化情況，以效率最高為原則，在指定流量比要求下，找出一種最適合黃磷脫砷工藝的尺寸參數。

1 計算模型

1.1 工作原理和結構尺寸參數

文丘里引射反應器主要由吸入管、吸入室、噴嘴、喉管(混合室)、擴散管、引流管等部分組成，其基本結構如圖1所示。它的工作原理是:一定壓力的流體(雙氧水)通過噴嘴后以一定速度射出，在射流紊動擴散作用及卷吸作用下，吸入低壓流體(液態黃磷)，兩股不同壓力的流體在喉管(混合室)內產生動量、質量和能量交換，工作流體速度及壓力減小，被吸入流體速度增加，在喉管出口處兩股流體的流速漸趨于均勻，混合流體經擴散管將部分動能轉化為壓能，增壓后通過引流管輸出。

圖1 文丘里引射反應器結構圖Fig.1 Structure of Ventu ri ejector reactor

在文丘里引射反應器的結構參數中，面積比受噴嘴出口直徑變化的直接影響，過小的噴嘴出口直徑能夠產生極大的射流速度和梯度，可以獲得較高的分散效果，但與此同時也會增加內部流場的局部阻力損失，進而改變反應器的性能。因此，針對反應器指定的流量比要求，應綜合考慮反應器操作的穩定性、高效性等因素來決定噴嘴的直徑大小。

文丘里引射反應器的結構尺寸根據黃磷脫砷工藝的流程處理需要，將雙氧水與液態黃磷的流量之比控制在10∶1到5∶1的范圍時脫砷效果較為良好，在此基礎上參照經驗公式［1］確定如表1所示。分別計算了經驗公式內的6、8、10和12 mm等4個噴嘴出口直徑下(其對應的喉管與噴嘴出口截面積之比分別為 9.00、5.06、3.24和 2.25)文丘里引射反應器的內部流場及性能，并進行了對比分析，在不考慮發生汽蝕的條件下，通過改變不同的噴嘴出口直徑和操作條件，確定對應的最適宜工作參數，為實際黃磷脫砷工藝的節能和高效提出可參考的數據范圍。

表1 結構尺寸參數Table 1 Geom etric param eters of VER

1.2 網格劃分和湍流模型

由于文丘里引射反應器的內部流場屬于不規則區域的有限空間射流，因此在整個計算域內，網格的劃分不宜采用單一的網格，而應根據速度梯度的大小來適當調整網格的疏密，以保證流場各部分的節點距離相對穩定，確保計算結果的可靠性。針對反應器結構來說，在噴嘴出口至喉管這一范圍內，兩種流體存在較大的梯度差異，二者在接觸混合的同時伴隨著劇烈的湍動和傳遞現象，必須考慮其紊流剪切作用和壓力脈動作用，因此在這一區域需要適當加大網格密度。同時，在近壁區域，由于邊界層的存在也應當適當加密網格。網格劃分形式如圖2所示。

圖2 網格結構示意圖Fig.2 G rid d ivision m odel of Venturi ejector reactor

求解器設置選擇分離式隱式求解器，Realizable k-ε雙方程紊流模型，壓力方程采用標準離散格式，多相流采用歐拉兩相流模型(流體的相關物性數據見表2)，而壓力-速度耦合方程采用 Phase Coupled Simp le格式，其余方程采用一階迎風離散格式［14］。在邊界條件的設定中，工作流體和引射流體進口均設為壓力進口條件，混合液出口設為壓力出口條件。引射流體的進口壓力考慮液態黃磷重力造成的靜壓，設置為8 kPa，出口壓力根據工程應用需要，保證混合液壓能可以克服沿程阻力到達指定儲罐，設置為18 kPa。壁面邊界采用計算效率高、實用性強的壁面函數法。在此基礎上對文丘里引射反應器全流場進行了數值模擬。

表2 工作介質物性數據Table 2 Properties of work ing m edium

計算過程中監視計算域的速度、壓力、湍流動能等變量，當監視參數區域穩定同時殘差達到預定值10-5，并且引射流體的質量流量基本保持不變，總的質量不平衡率小于總的質量流量的0.5%時，認為計算收斂。計算中沒有考慮溫度的影響。

2 性能參數

由于文丘里引射反應器的設計原理及工作過程與射流泵相似，故可以使用射流泵的相關性能參數作為衡量和評價該反應器的部分性能指標。

射流泵的流量、壓力及主要幾何參數用以下無因次函數表示為

式(1)～(4)中，Qs為吸入口流量(液態黃磷);Q0為工作流量(雙氧水);p0為工作壓力;ps為吸入口壓力;pc為出口壓力;ρ0為工作流體密度;ρs為引射流體密度;ρc為工作流體密度;v0為工作速度;vs為吸入口速度;vc為出口速度;z為位置高程;Ath為喉管截面積;A為噴嘴收縮處截面積。

3 噴嘴直徑對性能指標的影響

3.1 噴嘴直徑對q-h曲線和q-η曲線的影響

在中試實驗裝置中，調節雙氧水的入口壓力分別為0.2、0.3和0.4 MPa，實驗中測得雙氧水工作流量分別為1.13、1.41和1.65 m3/h，而在相同的操作條件下，數值模擬所得到的雙氧水工作流量分別為1.19、1.35和1.79 m3/h，對應的相對誤差分別為5.3%、4.3%和7.8%，相對誤差滿足要求，說明數值模擬的結果是可靠的。

根據工程需要，在引射壓力為8 kPa，出口背壓為18 kPa的條件下，改變工作壓力，分別對噴嘴直徑為 6、8、10和 12 mm(對應的面積比為 9.00、5.06、3.24和2.25，下同)進行數值計算，繪制出文丘里引射反應器的工作特性曲線如圖3和圖4。

圖3 噴嘴直徑對q-h曲線的影響Fig.3 The im pact of different nozzle on the curve of q-h

圖4 噴嘴直徑對q-η曲線的影響Fig.4 The im pact of different nozzle on the curve of q-η

由圖3和圖4可以看出，隨著噴嘴直徑的減小，對應面積比的增大，q-h曲線和 q-η曲線均向下偏移;而在相同流量比的條件下，小直徑的噴嘴由于局部阻力較大，導致壓力損失過多而不能獲得較高的效率和較大的壓力比，說明在其他尺寸參數相同的情況下，使用大直徑的噴嘴有助于提升反應器的工作效率和壓力比，但是，由于噴嘴直徑的增大，導致了噴嘴出口速度的降低，由此可能造成高速射流效果不佳，進而無法獲得充分分散的混合反應結果。

另一方面，從圖3和圖4中可以看出，高的壓力比能夠獲得同樣高的效率，這也驗證了Ning和Satofuka［15］的研究理論成果。在圖3中，每一個面積比都有一條對應的效率曲線，且都有一個最適宜值，但是，不同面積比的最佳效率值是不同的。設定從最高效率下降5%為標準(相對值)來確定每一面積比下較為適宜的流量比范圍，同時列出最高效率點處對應的最適宜流量比，結果如表3所示。

表3 不同面積比下的最適宜流量比Table 3 The best flow rate ratio w ith different area ratio

針對不同的流量比需求，可按工程需要，確定工作流體和引射流體的混合反應比，進而根據表3的最適宜流量比范圍選擇最適宜的面積比，在不考慮高速噴射造成的分散混合效應前提下，可以保證文丘里引射反應器在高效率區進行工作。

3.2 噴嘴直徑對壓降的影響

圖5是不同面積比下工作流體通過文丘里引射反應器噴嘴的壓力降，其中工作流體的入口壓力設置為0.1 MPa。

圖5 噴嘴直徑對壓力降的影響Fig.5 The im pact of d ifferent nozzle on p ressu re d rop

由圖5可以看出，面積比越大(即噴嘴直徑越小)，則壓力降越大，而壓力降可以用來評價混合效果的好壞，因此，壓力降越大，高速射流造成的卷吸作用和紊動擴散作用越強烈，進入喉管的引射流體的動量越大，速度梯度和壓力梯度越大，因此工作流體對引射流體的分散效果越好，混合效果也越好。因此，對于接觸反應的分散效果要求嚴格的工程問題，應該在綜合評價驅動力、能量供給等能耗問題的條件下，盡可能選擇高面積比的反應器，以獲得滿意的分散反應效果。

3.3 噴嘴直徑對軸向速度和壓力的影響

圖6和圖7是不同面積比下工作流體通過文丘里引射反應器后軸向速度和壓力的變化趨勢圖，其中工作流體的入口壓力設置為0.1 MPa，截面位置Z=0表示噴嘴出口截面。可以看出，在喉管和擴散管的部分，4個面積比對應的曲線差別最大，速度差異和壓力差異最為明顯;在喉管段，速度衰減較為突出，說明在文丘里引射反應器的喉管內，工作流體和引射流體發生劇烈的傳質傳能現象，二者進行混合反應;在擴散管段，隨著軸向距離的增加，工作流體速度降低，引射流體速度增加，在喉管內二者進一步混合，流速逐漸趨于均勻，混合液的動能被轉化為壓能傳輸出去。

圖6 噴嘴直徑對軸向速度的影響Fig.6 The im pact of different nozzle on axial velocity

圖7 噴嘴直徑對軸向壓力的影響Fig.7 The im pact of d ifferent nozzle on axial p ressure

另一方面，在相同工作流體入口壓能的驅動下，不同的面積比對應的反應器喉管及擴散管內部的速度和壓力梯度變化也不盡相同。圖6顯示面積比大的條件下對應的速度梯度最大，并且速度衰減也最為強烈。原因是小直徑噴嘴的局部能量損失太大，導致工作流體的壓能在高速射流后，將能量傳遞給引射流體的同時迅速衰減。由于文丘里引射反應器內的高速射流卷吸效應屬于有限空間內的射流，與自由剪切湍流不同，因此反應器內部流場受到噴嘴直徑的影響較為顯著。

3.4 中試實驗與模擬結果的比較

選取噴嘴出口直徑為10 mm的文丘里引射反應器為例，根據貴州某廠在中試實驗中得到的相關數據，將實際中雙氧水的工作流量與相同操作條件下模擬得到的雙氧水流量進行對比，得到的結果如表4所示。

表4 實驗和模擬的雙氧水流量對比Table 4 The contrast of Q0 between test and sim u lation

對比實驗和模擬數據可以得出，在相同的參數設置條件下，得到的雙氧水流量大小基本吻合，相對誤差基本保持在5%以下，說明模擬的結果是可靠的。

4 結論

應用FLUENT軟件，采用 Phase Coupled Simple算法和Realizable k-ε湍流模型，模擬了在其他尺寸參數相同的條件下不同噴嘴直徑表現出的文丘里引射反應器工作性能和內部流場變化情況。主要結論如下。

1)噴嘴直徑直接影響到文丘里引射反應器的工作性能和內部流場。在特定的結構尺寸條件下，噴嘴直徑越大(即面積比越小)，反應器所形成的壓力比越大，驅動能力越強，同時工作效率也越高。

2)在相同的工作壓力條件下，噴嘴直徑越小，產生的壓力降越大，局部能量損失也越大。另一方面，對應的接觸混合分散效果也越好。對于分散效果要求嚴格的應用場合，應該綜合考慮反應器的設計和選型，盡量保證反應器在獲得滿意分散效果的前提下，在高效工作區工作。

3)在相同的驅動工作壓力條件下，沿軸線方向的速度分布和壓力分布不盡相同。噴嘴直徑越小，對應的速度和壓力梯度越大，當工作流體從噴嘴高速射流出后，速度和壓力迅速衰減，在這過程中工作流體和引射流體進行強烈的傳質傳能接觸反應。

［1］陸宏圻.噴射技術理論及應用［M］.武漢:武漢大學出版社，2004 Lu Hongqi.The theory and application for spray techniques［M］.Wuhan:Wuhan University Press，2004(in Chinese)

［2］劉軍，李軍，任永勝，等.氧化增強劑對黃磷脫砷效果的影響［J］.無機鹽工業，2010，42(2):55-57 Liu Jun，Li Jun，Ren Yongsheng，et al.Effect of oxidation reinforcer on removing arsenic from yellow phosphorus［J］.Inorganic Chemicals Industry，2010，42(2):55-57(in Chinese)

［3］査坐統，蘇毅，李國斌.黃磷精制廢液脫砷的實驗研究［J］.化工科技，2009，17(4):13-16 Zha Zuotong，Su Yi，LiGuobin.The experimental study on the arsenic removal of waste liquid of refining yellow phosphorus［J］.Science and Technology in Chemical Industry，2009，17(4):13-16(in Chinese)

［4］龍新平，程茜，韓寧，等.射流泵最佳喉嘴距的數值模擬［J］.核動力工程，2008，29(1):35-38 Long Xinping，Cheng Qian，Han Ning，et al.Numerical simulation of optimal nozzle-to-throat clearance of jet pump［J］.Nuclear Power Engineering，2008，29(1):35-38(in Chinese)

［5］Long X，Han N，Chen Q.Influence of nozzle exit tip thickness on the performance and flow field of jet pump［J］.Mechanical Science and Technology，2008，22:1 959-1 965

［6］王常斌，林建忠，石興.射流泵最佳參數的確定方法［J］.流體機械，2004，32(9):21-25 Wang Changbin，Lin Jianzhong，Shi Xing.Method of optimal parameter ascertainment of jet pump［J］.Fluid Machinery，2004，32(9):21-25(in Chinese)

［7］王常斌，林建忠，石興.射流泵湍流場的數值模擬與實驗研究［J］.高校化學工程學報，2006，20(2):175-179 Wang Changbin，Lin Jianzhong，Shi Xing.Numerical simulation and experiment on the turbulent flow in the jet pump［J］.Chemical Engineering of Chinese Universities，2006，20(2):175-179(in Chinese)

［8］龍新平，關運生，韓寧，等.可調式射流泵性能的數值模擬［J］.排灌機械，2008，26(6):1-5 Long Xinping，Guan Yunsheng，Han Ning，et al.Numerical simulation on performance of self-regulating jet pump［J］.Drainage and Irrigation Machinery，2008，26(6):1-5(in Chinese)

［9］Ajmal S，Im ran R C，Mansoor H I.Experimental and numerical analysis of steam jet pump［J］.Multiphase Flow，2011，37:1 305-1 314

［10］董敬亮，王曉冬.水蒸氣噴射泵混合室內邊界層脫離現象的計算流體力學研究［J］.真空科學與技術學報，2010，30(4):455-458 Dong Jingliang，Wang Xiaodong.Computational fluid dynamic study of boundary layer separation in m ixing chamber of steam ejector［J］.Vacuum Science and Technology，2010，30(4):455-458(in Chinese)

［11］Eames IW，Worall M，Wu S.An experimental investigation into the integration of a jet-pump refrigeration cycle and a novel jet-spay thermal ice storage system［J］.Applied Thermal Engineering，2013，53(2):285-290

［12］Kanjanapon C，Satha A.An experimental investigation of a steam ejector refrigerator:The analysis of the pressure profile along the ejector［J］.Applied Thermal Engineering，2004，24:311-322

［13］David S，Zine A，Mohamed O.An experimental investigation of an ejector for validating numerical simulations［J］.Refrigeration，2011，34:1 717-1 722

［14］王福軍.計算流體動力學分析［M］.北京:清華大學出版社，2004 Wang Fujun.The analysis of computational fluid dynam ics［M］.Beijing:Tsinghua University Press，2004

［15］Ning T，Satofuka N.Nihon Kikai Gakkai Ronbunshu［J］.Transactions of the Japan Society of Mechanical Engineers:Part B，1997，63:1 243-1 250

［16］王曉靜，江增，唐磊.黃磷脫砷裝置中分相器的 CFD數值模擬研究［J］.化學工業與工程，2015，32(4):44-49 Wang Xiaojing，Jiang Zeng，Tang Lei.CFD numerical simulation on the phase separator equipment for dearsenication of yellow phosphorus［J］.Chem ical Industry and Engineering 2015，32(4):44-49(in Chinese)