填料塔中亞硝酸甲酯再生過程模擬

李文龍，王瑋涵，李振花，馬新賓

?

填料塔中亞硝酸甲酯再生過程模擬

李文龍，王瑋涵，李振花，馬新賓

（天津大學(xué)化工學(xué)院，綠色合成與轉(zhuǎn)化教育部重點實驗室，天津 300072）

在氣相偶聯(lián)合成草酸酯工藝過程中，亞硝酸甲酯（MN）再生過程與偶聯(lián)反應(yīng)有效匹配是實現(xiàn)整個系統(tǒng)綠色無污染的關(guān)鍵，也是影響草酸酯合成工業(yè)放大過程的關(guān)鍵一步。對于MN再生過程的工藝優(yōu)化和宏觀動力學(xué)方面進行了較多研究和報道，而對MN再生過程的數(shù)學(xué)模擬研究較少。通對MN再生過程中發(fā)生的主要化學(xué)反應(yīng)進行分析簡化，并結(jié)合MN再生反應(yīng)動力學(xué)參數(shù)，在填料塔中建立了MN再生反應(yīng)數(shù)學(xué)模型，并計算了主要工藝參數(shù)對MN收率的影響。模擬結(jié)果表明，隨溫度升高、氣體流速增大、N2體積分?jǐn)?shù)增大，MN的收率減小，液體流速增大和NO/O2摩爾比增大，MN收率增大。通過實驗驗證了該模型的準(zhǔn)確性，MN收率的模擬結(jié)果與實驗數(shù)據(jù)的最大相對偏差為-4.39%，為MN再生反應(yīng)過程的設(shè)計和操作提供了理論依據(jù)；最后，對于兩種不同的NO和O2混合方式進行了比較，通過模擬計算比較了這兩種不同操作方式對MN收率的影響。

亞硝酸甲酯；填料塔；數(shù)學(xué)模擬；再生；吸收

引 言

草酸酯是重要的有機化工原料，廣泛應(yīng)用于精細(xì)化工制備各種染料、醫(yī)藥、重要的溶劑、萃取劑以及各種中間體[1]。目前，草酸酯的合成工藝有很多，其中CO氣相催化偶聯(lián)反應(yīng)因其反應(yīng)條件溫和、催化劑的損失小、固定投資和動力消耗低等優(yōu)點，成為國內(nèi)一碳化學(xué)及有機化工領(lǐng)域中的研究熱點。CO氣相催化偶聯(lián)反應(yīng)制取草酸二甲酯（dimethyl oxalate，DMO）包括CO偶聯(lián)反應(yīng)和亞硝酸甲酯（methyl nitrite， MN）再生反應(yīng)，反應(yīng)方程式如下：

偶聯(lián)反應(yīng)

再生反應(yīng)

總反應(yīng)

偶聯(lián)反應(yīng)以再生反應(yīng)生成物MN為載氧體，生成物NO作為再生反應(yīng)的反應(yīng)物循環(huán)使用，因此，再生反應(yīng)和偶聯(lián)反應(yīng)的有效匹配是實現(xiàn)整個系統(tǒng)綠色無污染的關(guān)鍵[2]。國內(nèi)外的學(xué)者[3-5]對CO偶聯(lián)反應(yīng)過程已做了較多的研究和報道，而針對MN再生反應(yīng)所做的研究和報道較少，且一般都局限在工藝條件和宏觀動力學(xué)方面[6-8]，對MN再生反應(yīng)的數(shù)學(xué)模型的研究較少。

MN再生反應(yīng)是NO和O2與甲醇進行的氣液反應(yīng)，屬于NO化學(xué)吸收的范疇，是復(fù)雜的吸收-反應(yīng)過程之一。目前，已見報道的關(guān)于NO吸收的數(shù)學(xué)模型大多數(shù)是用于描述NO在水[9-10]、酸液[11-12]、堿液[13-15]和氧化劑溶液[16]中的吸收過程，而針對NO在甲醇中的吸收和反應(yīng)過程的研究很少，主要原因是缺乏可靠的MN再生反應(yīng)動力學(xué)模型。柳 剛[17]在雙攪拌反應(yīng)器中對合成MN的宏觀動力學(xué)進行了研究，并根據(jù)MN動力學(xué)的研究結(jié)果，在填料塔中對合成MN進行了模擬計算。考察了反應(yīng)溫度、氣體流量、混合氣體中惰性氣體的含量、NO與NO2比例、N2O3的轉(zhuǎn)化率、填料的類型對所需填料層高度的影響，為工業(yè)放大裝置的設(shè)計、操作提供了理論指導(dǎo)。郭新連等[18]采用Aspen軟件模擬 MN 反應(yīng)塔的工藝，根據(jù)實際裝置的反應(yīng)轉(zhuǎn)化率與反應(yīng)時間，驗證并修正了MN動力學(xué)參數(shù)，同時利用Aspen軟件的靈敏度，分析了優(yōu)化塔的操作參數(shù)，模擬結(jié)果對 MN 反應(yīng)塔的工藝設(shè)計與工業(yè)化有一定指導(dǎo)意義。魯文質(zhì)[19]同樣采用Aspen模擬計算了反應(yīng)條件(溫度和壓力)和精餾條件(反應(yīng)精餾塔的精餾段和提餾段塔板數(shù)、回流量)對反應(yīng)結(jié)果的影響，并將反應(yīng)精餾塔小試實驗與模擬結(jié)果進行對比，模擬結(jié)果偏高。Li等[8]在獲得了亞硝酸甲酯再生反應(yīng)動力學(xué)模型的基礎(chǔ)上，建立了描述半連續(xù)填料鼓泡塔中MN再生反應(yīng)過程的數(shù)學(xué)模型，并很好地預(yù)測了填料鼓泡塔中的MN再生反應(yīng)結(jié)果，模擬的結(jié)果與實驗結(jié)果吻合度良好。本文采用已報道的MN再生反應(yīng)動力學(xué)參數(shù)，在填料塔中建立MN再生的數(shù)學(xué)模型，分析主要工藝參數(shù)對MN再生的影響，并驗證模型的準(zhǔn)確性，從而為MN再生過程的設(shè)計和操作提供理論依據(jù)，最后比較了兩種不同操作方式(NO與O2先混和后混)對MN收率的影響。

1 數(shù)學(xué)模型建立

1.1 氣相反應(yīng)

MN再生反應(yīng)的氣相原料為O2、NO，以N2作為稀釋氣和保護氣，NO與O2發(fā)生不可逆的氧化反應(yīng)生成NO2，氧化速率對于NO是二級反應(yīng)，對于O2是一級反應(yīng)，其反應(yīng)方程式、速率方程相關(guān)參數(shù)[20]如下

同時氣相中也進行一些可逆的反應(yīng)，如

對于可逆反應(yīng)式(5)和式(6)，反應(yīng)速率極快，能夠迅速達到平衡狀態(tài)，可以認(rèn)為氣相主體各物質(zhì)很快達到平衡狀態(tài)[21]。雖然各個物質(zhì)很快達到平衡，但是受到熱力學(xué)的限制，N2O3和N2O4在氣相中的量相對于NO和NO2的量卻很少。

1.2 液相反應(yīng)

氣相主體的氮氧化合物（NO）與液相甲醇接觸時會發(fā)生如下反應(yīng)

反應(yīng)式(7)是一個快速反應(yīng)，是生成MN的主要反應(yīng)，反應(yīng)式(8)和式(9)除了生成MN外，還生成副產(chǎn)物HNO3，反應(yīng)過程中檢測到的HNO3量相對于MN的生成量是非常小的[6,8]，因此可以忽略反應(yīng)式(8)和式(9)。由于文獻中缺乏N2O3與甲醇反應(yīng)的動力學(xué)數(shù)據(jù)，且NO和NO2的混合氣與N2O3氣體具有相似的物化性質(zhì)，所以用NO和NO2混合氣代替N2O3與甲醇反應(yīng)，式(7)可以改為

對于式(10)，Liu等[6]和Li等[8]在氣液雙攪拌裝置中對MN再生反應(yīng)進行了詳細(xì)的研究，得到了其動力學(xué)模型參數(shù)，但是Liu等[6]是在NO:NO21:1的條件下獲得MN再生反應(yīng)的動力學(xué)參數(shù)，使用范圍受到限制，而Li等[8]研究所得的MN動力學(xué)模型不受NO與NO2濃度的限制，因此本文模擬計算采用該動力學(xué)參數(shù)，其表達式如下為

1.3 填料塔中模型的建立

在填料塔中（圖1）進行MN再生反應(yīng)模擬，假設(shè)：①氣相遵循理想氣體行為，并且以活塞流的形式從塔底進入，從塔頂流出；②液相采用無水甲醇，從塔頂經(jīng)液體分布與氣相逆向接觸反應(yīng)，從塔底流出，液相中生成的副產(chǎn)物水的量相對于甲醇的量很少，近似認(rèn)為甲醇的濃度恒定；③反應(yīng)生成的MN很快擴散到氣相中，其在液相中的溶解損失量可以忽略；④全塔處于穩(wěn)定操作狀態(tài)下，且持液量不隨塔高變化。

圖 1 填料塔

氣相活塞流模型，各個物質(zhì)經(jīng)過微元高度d后的物料衡算方程為

1.4 設(shè)計參數(shù)的估算

在模擬過程中，氣含率、有效傳質(zhì)比表面積采用以下經(jīng)驗方程式進行估算。

（1）有效傳質(zhì)比表面積 對于散堆填料的有效傳質(zhì)比表面積，Puranik等[22]給出了相關(guān)的關(guān)聯(lián)式進行估算

（2）氣含率 在填料塔中，氣體的氣含率可以通過填料的空隙率減去填料的持液量求得，即

持液量采用Billet[23]模型進行估算。

持液量

靜持液量

當(dāng)L＜5時

當(dāng)L≥5時

液體Reynolds數(shù)

液體Froude數(shù)

其中

2 實驗結(jié)果與討論

在中型填料塔中進行MN再生反應(yīng)的模擬計算，該塔主要的設(shè)計參數(shù)和操作參數(shù)為：填料塔內(nèi)徑為0.6 m；采用金屬不銹鋼鮑爾環(huán)填料25 mm，無催化作用，填料高度2 m；液體流速為4 mm·s-1；氣體流速為0.1 m·s-1，惰性氣體占總氣體流量的70%～85%，NO和O2從塔底進入，摩爾比(4:1)～(7:1)，壓力101.3 kPa。

2.1 參數(shù)的敏感性

2.1.1 溫度的影響 模擬反應(yīng)條件：壓力101.3 kPa，氣體流速G0.1 m·s-1，液體流速L4 mm·s-1，惰性氣體N2占總氣體體積含量為80%，NO與O2的摩爾比為4:1，改變溫度參數(shù)，所得氧氣轉(zhuǎn)化率和MN收率的模擬結(jié)果如圖2所示。從圖2可以看出，在選定的溫度范圍內(nèi)，溫度對于MN收率的影響不大，隨著溫度升高，MN的收率略有降低。對于反應(yīng)式(4)，其反應(yīng)速率隨著溫度升高而減小，溫度升高，NO的氧化速率降低，NO2的生成量減少，MN生成速率降低，致使MN收率隨溫度升高而降低。由反應(yīng)式(10)的動力學(xué)參數(shù)可知MN生成速率卻是隨著溫度升高而增加的，也正是由于這兩個相反的作用，致使MN再生反應(yīng)受溫度的影響變化不大。

圖2 溫度對MN再生反應(yīng)結(jié)果的影響

(G0.1 m·s-1,L4 mm·s-1, NO/O2molar ratio 4:1,80%)

2.1.2 進料中NO與O2摩爾比的影響 模擬反應(yīng)條件：壓力101.3 kPa，反應(yīng)溫度35℃，氣體流速G0.1 m·s-1，液體流速L4 mm·s-1，惰性氣體N2占總氣體體積含量為80%，改變NO與O2的摩爾比，所得氧氣轉(zhuǎn)化率和MN收率的模擬結(jié)果如圖3所示。從圖3可以看出，MN的收率隨著NO/O2的摩爾比增加而增大。較高的NO/O2摩爾比，提高了O2的轉(zhuǎn)化率和NO2的生成速率，進而提高了MN的收率。在合成草酸酯的工藝中，MN再生反應(yīng)中未完成轉(zhuǎn)化的O2會造成后續(xù)偶聯(lián)反應(yīng)催化劑的失活，為了使O2盡可能完全轉(zhuǎn)化，一般采用較高的NO/O2摩爾比。

圖3 NO與O2摩爾比對MN再生反應(yīng)結(jié)果的影響

(G0.1 m·s-1,L4 mm·s-1,80%,35℃)

2.1.3 氣體流速的影響 模擬反應(yīng)條件：壓力101.3kPa，反應(yīng)溫度35℃，液體流速L4 mm·s-1，惰性氣體N2占總氣體體積含量為80%，NO與O2的摩爾比為4:1，改變氣體流速，所得氧氣轉(zhuǎn)化率和MN收率的模擬結(jié)果如圖4所示。從圖4可以看出，隨空塔氣速的增加，O2的轉(zhuǎn)化率略有減小，MN的收率減小。這是由于氣速增加，氣體在填料塔中的停留時間變短，氣液接觸時間變短。NO的氧化反應(yīng)[式(4)]是在氣相中發(fā)生的，受此影響較小，而MN再生反應(yīng)是在氣液兩相界面處發(fā)生的，停留時間變短，氣相中的反應(yīng)物與液體甲醇接觸反應(yīng)時間變短，MN的生成速率減小，使MN的收率降低。

圖4 空塔氣速對MN再生反應(yīng)結(jié)果的影響

(L4 mm·s-1,80%, NO/O2molar ratio 4:1,35℃)

2.1.4 液體流速的影響 模擬反應(yīng)條件：壓力101.3 kPa，反應(yīng)溫度35℃，氣體流速G0.1 m·s-1，惰性氣體N2占總氣體體積含量為80%，NO與O2的摩爾比為4:1，改變液體流速，所得氧氣轉(zhuǎn)化率和MN收率的模擬結(jié)果如圖5所示。從圖5中可以看出，液體流速增加，O2的轉(zhuǎn)化率降低，MN的收率增大。液體流速增加，液體的噴淋密度增大，填料層中液體的持液量增加，填料間空隙減小，NO的氧化反應(yīng)空間減小，氧化速率降低，O2轉(zhuǎn)化率減??；對于MN生成反應(yīng)，填料中液體持液量的增加，增大了氣液接觸的有效傳質(zhì)面積，促進了MN的生成，因此MN收率增加。

圖5 液體流速對MN再生反應(yīng)結(jié)果的影響

(G0.1 mm·s-1,80%, NO/O2molar ratio 4:1,35℃)

2.1.5 進料氣中惰性氣體的百分含量的影響 模擬反應(yīng)條件：壓力101.3 kPa，反應(yīng)溫度35℃，氣體流速G0.1 m·s-1，液體流速L4 mm·s-1， NO與O2的摩爾比為4:1，改變N2在總氣體流量中的體積分?jǐn)?shù)，所得氧氣轉(zhuǎn)化率和MN收率的模擬結(jié)果如圖6所示。由圖6可以看出，隨著N2含量減少，O2轉(zhuǎn)化率和MN收率均增大。在保持氣相總流量不變的前提下，降低N2含量可以增大O2和NO在氣相中的濃度，從而增大NO氧化反應(yīng)速率，O2的轉(zhuǎn)化率增大，生成的NO2在氣相中的濃度也相應(yīng)增大，使得氣相中NO和NO2的濃度相應(yīng)增大，可有效促進其在液相中的擴散，進而促進MN再生反應(yīng)的進行，因此MN收率增加。

圖6 N2的體積含量對MN再生反應(yīng)結(jié)果的影響

(G0.1 mm·s-1，L4 mm·s-1，NO/O24:1，35℃)

2.2 模擬與實驗結(jié)果的對比

受實驗室條件的限制，本文采用文獻報道的MN再生反應(yīng)的實驗條件進行模擬，并將模擬結(jié)果與文獻中的實驗結(jié)果進行對比。具體的實驗條件和設(shè)備參數(shù)如下[19]：塔徑50 mm，填料采用3 mm×3 mm不銹鋼絲網(wǎng)θ環(huán)，填料高度2 m；原料氣總流量2.789 L·min-1（標(biāo)況），體積組成NO 9.68%、O22.47%、N287.7%，甲醇流量3.6 g·min-1，壓力101.3 kPa。

從圖7 (a)可以看出，O2轉(zhuǎn)化率的模擬結(jié)果比實驗結(jié)果偏低，而MN收率的模擬結(jié)果與實驗結(jié)果能較好地吻合[圖7 (b)]。這是由于模擬計算過程中MN的生成是以NO和NO2為反應(yīng)物，所以MN的生成會消耗部分NO，降低NO在氣相的濃度，使得NO的氧化速率減小，O2的轉(zhuǎn)化率降低。另外，實際反應(yīng)過程中會生成部分副產(chǎn)物HNO3，增加了O2的消耗，而這部分消耗在模擬計算中被忽略了，因此使得模擬計算所得O2的轉(zhuǎn)化率比實際反應(yīng)結(jié)果所得O2的轉(zhuǎn)化率較低，模擬結(jié)果與實驗結(jié)果偏差見表1。從表1可以看出，MN收率和O2轉(zhuǎn)化率的相對偏差都在允許的偏差范圍內(nèi)。

圖7 計算值與實驗值的對比

表1 不同溫度下模擬結(jié)果與實驗結(jié)果偏差

3 先混操作

前文所描述的是將NO和O2分別從塔底引進，氧化反應(yīng)和MN再生反應(yīng)同時在塔中進行，這種操作叫后混方式[圖8 (a)]，與之對應(yīng)的一種操作是在填料塔之前增加一個氧化容器，使得NO與O2預(yù)先充分混合，然后再引入填料塔中進行反應(yīng)，這種操作叫做先混方式[圖8 (b)]，其優(yōu)點是可以降低反應(yīng)器的反應(yīng)高度，減少填料用量。根據(jù)動力學(xué)方程和氣液反應(yīng)機理可以定性得出如下結(jié)論：在使用同樣操作條件和設(shè)備尺寸下，與后混方式相比，采取先混方式操作時MN的收率較高，反應(yīng)速率較快。然而鑒于文獻中缺乏兩種操作情況下的定量比較結(jié)果，本文對于先混方式的反應(yīng)結(jié)果也進行了模擬計算，并在不同反應(yīng)溫度下與后混反應(yīng)結(jié)果進行了比較。

圖8 后混和先混操作方式

為了計算先混方式下的反應(yīng)結(jié)果，需要建立相應(yīng)的模型并進行模擬計算。在建立先混方式下的數(shù)學(xué)模型之前，同樣需要做一些假設(shè)條件，除了后混反應(yīng)所做的假設(shè)之外，增加一條假設(shè)，即假設(shè)NO與O2在進入填料塔之前已經(jīng)反應(yīng)完全，O2被完全消耗。因此，先混反應(yīng)各物質(zhì)在填料塔中的物料衡算微分方程為

后混和先混操作的模擬采用相同的操作條件：填料塔內(nèi)徑為0.6 m；采用金屬不銹鋼鮑爾環(huán)填料25 mm，無催化作用，填料高度2 m；壓力101.3 kPa，氣體流速G0.1 m·s-1，液體流速L4 mm·s-1，惰性氣體N2占總氣體體積含量為80%，NO與O2的摩爾比為4:1。則不同反應(yīng)溫度下兩種不同操作方式的模擬結(jié)果如圖9所示。

圖9 不同溫度下先混與后混操作時MN收率對比

(G0.1 m·s-1,L4 mm·s-1, NO/O2molar ratio 4:1,80%)

從圖9可以看出，不同溫度下，先混操作所得MN的收率要比后混操作的高8.82%～13.36%。這是由于采用先混操作，反應(yīng)器底部入口和反應(yīng)床層氣相中的NO2濃度比后混操作的NO2濃度大，先混操作的MN生成速率比后混的大。同時，MN的生成反應(yīng)與NO氧化反應(yīng)屬于連串反應(yīng)，MN生成速率在一定程度上受到NO氧化速率的影響，造成后混操作比先混操作MN的收率低。另外，從圖9還可以看出，隨著反應(yīng)溫度升高，MN收率間的差距有逐漸增大的趨勢。這是因為在先混操作中，填料塔中只進行MN再生反應(yīng)，而MN再生反應(yīng)的生成速率隨溫度升高而增加，因此溫度升高，MN的收率增加。對于后混操作，NO氧化生成氮氧化合物與液相中氮氧化合物與甲醇反應(yīng)生成MN的反應(yīng)具有連串關(guān)系，溫度升高降低了NO的氧化速率，使NO2的生成濃度減小，MN的生成量減小，因此，MN收率隨反應(yīng)溫度升高而減少。

采用兩種不同操作方式時，MN收率隨反應(yīng)床層的變化如圖10所示，這表明采用兩種不同的進料方式，達到同樣的MN收率所使用的填料高度不同。表2為上述反應(yīng)條件下采用兩種操作方式，MN收率為75%時，所需的填料層高度。從表2中可以看出，在同一溫度下，后混操作所需的填料高度比先混操作所需要的大很多，并且隨著反應(yīng)溫度升高，差距增大，45℃時，后混操作所需填料高度約為先混操作的2倍。因此，采用先混操作有利于減少生產(chǎn)設(shè)備尺寸，減少填料量。

圖10 先混與后混操作中MN收率隨床層高度變化

(G0.1 m·s-1,L4 mm·s-1, NO/O2molar ratio 4:1,80%,35℃)

表2 不同溫度下兩種操作所需的填料高度

4 結(jié) 論

（1）在填料塔中建立了MN再生的數(shù)學(xué)模型，模擬計算結(jié)果表明：隨著溫度升高、氣體流速增大、N2體積分?jǐn)?shù)增大，MN的收率減??；隨著液體流速增大和NO/O2摩爾比增大，MN收率增大。

（2）將模擬結(jié)果與實驗數(shù)據(jù)進行對比，MN收率的模擬值與實驗數(shù)據(jù)的最大相對偏差為-4.39%；O2轉(zhuǎn)化率的模擬值與實驗數(shù)據(jù)的最大相對偏差 為-8.15%。表明所建立的數(shù)學(xué)模型對于MN再生反應(yīng)器放大設(shè)計有一定的指導(dǎo)意義。

（3）通過比較兩種不同操作方式對MN收率和填料層高度的影響發(fā)現(xiàn)，與后混操作方式相比，先混操作可獲得較高的MN收率；MN收率一定時，先混操作所用的填料高度遠(yuǎn)低于后混操作。

符 號 說 明

a，ae——填料有效傳質(zhì)比表面積，a的單位cm-1，ae的單位m-1 ap——填料幾何比表面積，m-1 CE——甲醇的濃度，mol·L-1 Ch——常數(shù)，與填料結(jié)構(gòu)有關(guān) Di——物質(zhì)i在液相中擴散系數(shù)，cm2·s-1 FrL——液體Froude數(shù) g——重力加速，m·s-2 Hi——物質(zhì)i的亨利系數(shù)，kPa·cm3·mmol-1 hL——填料層持液量 hLS——靜持液量 ko——氧化速率常數(shù)，kPa-2·s-1 kr——MN生成速率常數(shù)，ml2·mmol-2·s-1 Pi——氣相主體中物質(zhì)i的壓力，kPa R——摩爾氣體常數(shù)，8.314 J·mol-1·K-1 Ro——NO氧化速率，kmol·m-3·s-1 Rr——MN生成速率，mmol·ml-2·s-1 ReL——液體Reynolds數(shù) T——溫度，K uf——泛點氣速，m·s-1 uG——空塔氣速，cm·s-1 uL——液體流速，m·s-1 z——填料層高度，cm εG——空隙率 μL——液體的黏度，Pa·s ρL——液體密度，kg·m-3 σc——填料材質(zhì)表面張力，N·m-1 σL——液體表面張力，N·m-1

References

[1] Fang Jingang (房金剛), Xu Genhui (許根慧). Research progress in the synthesis of oxalate [J].(化學(xué)推進劑與高分子材料), 2004, 2 (2): 18-21

[2] Wang Weihan (王瑋涵), Li Zhenhua (李振花), Ma Xinbin (馬新賓). Regeneration process of methyl nitrite in a bubble reactor [J].(化學(xué)反應(yīng)工程與工藝), 2011, 27 (5): 450-453

[3] Meng Fandong, Xu Genhui, Guo Rongqun, Yan Hongfei, Chen Manqiao. Kinetic study of carbon monoxide coupling reaction over supported palladium catalyst [J].:, 2004, 43 (6): 785-790

[4] Wang Baowei, Ma Xinbin, Xu Genhui, Mao Lanfang. Optimization of cyclical process for CO coupling regeneration reactions [J]., 2000, 24 (2-7): 1337-1341

[5] Gao X, Zhao Y, Wang S, Yin Y, Wang B, Ma X, Ji Y, Liu G, Li W, Xiao W D. A Pd-Fe/α-Al2O3/cordierite monolithic catalyst for CO coupling to oxalate [J]., 2011, 66 (15): 3513-3522

[6] Liu Gang, Ji Yang, Li Wei, Berski S, Latajka Z, Gordon A J. Kinetic study on methyl nitrite synthesis from methanol and dinitrogen trioxide [J]., 2010, 157 (2/3): 483-488

[7] Wang Huajun, Li Guangxing. Kinetic study on the synthesis of ethyl nitrite by the reaction of C2H5OH, O2, and NO in a trickle bed reactor [J]., 2010, 163 (3): 422-428

[8] Li Zhenhua, Wang Weihan, Lü Jing, Ma Xinbin. Modeling of a packed bubble column for methyl nitrite regeneration based on reaction kinetics and mass transfer [J]., 2013, 52 (8): 2814-2823

[9] Ramanand S B, Phaneswara Rao D, De Paiva J L, Kachan G C. Modelling and simulation of NOabsorption into water in a countercurrent flow packed column [J]., 1996, 20 (8): 1059-1063

[10] Yasuda M, Tsugita N, Ito K, Yamauchi S, Glomm W R, Tsuji I, Asano H. High-efficiency NOabsorption in water using equipment packed with a glass fiber filter [J]., 2011, 45 (5): 1840-1846

[11] Loutet K G, Mahecha-Botero A, Boyd T, Buchi S, Reid D, Brereton C M H. Experimental measurements and mass transfer/reaction modeling for an industrial NOabsorption process [J]., 2011, 50 (4): 2192-2203

[12] Li Y, Liu Y Z, Zhang L Y, Su Q, Jin G L. Absorption of NOinto nitric acid solution in rotating packedbed [J]., 2010, 18 (2): 244-248

[13] Patwardhan J A, Pradhan M P, Joshi J B. Simulation of NOgas absorption under adiabatic condition and comparison with plant data [J]., 2002, 57 (22/23): 4831-4844

[14] Pradhan M P, Joshi J B. Absorption of NOgases in plate column: selective manufacture of sodium nitrite [J]., 2000, 55 (7): 1269-1282

[15] Pradhan M P, Joshi J B. Absorption of NOgases in aqueous NaOH solutions: selectivity and optimization [J]., 1999, 45 (1): 38-50

[16] Thomas D, Vanderschuren J. The absorption-oxidation of NOwith hydrogen peroxide for the treatment of tail gases [J]., 1996, 51 (11): 2649-2654

[17] Liu Gang (柳剛). Studies on the synthesizing methyl nitrite [D]. Shanghai: East China University of Science and Technology, 2010

[18] Guo Xinlian (郭新連), Zhang Jingya (張靜雅), Ye Liping (葉麗萍), Fu Zhonghe (付中赫). The simulate calculation and optimization for the reactive distillation process of producing methyl nitrite [J].(計算機與應(yīng)用化學(xué)), 2010, 27 (11): 1485-1488

[19] Lu Wenzhi (魯文質(zhì)). Simulation on ethylene glycol synthesis by super DMO-based technology [D]. Shanghai: Shanghai Jiao Tong University, 2006

[20] Suchak N J, Joshi J B. Simulation and optimization of NOabsorption system in nitric acid manufacture [J]., 1994, 40 (6): 944-956

[21] Pradhan M P, Suchak N J, Walse P R, Joshi J B. Multicomponent gas absorption with multiple reactions: modelling and simulation of NOabsorption in nitric acid manufacture [J]., 1997, 52 (24): 4569-4591

[22] Puranik S S, Vogclpohl A. Effective interfacial area in irrigated packed column [J]...., 1974, 29 (2): 501-507

[23] Billet Reinhard. Packed Towers (填料塔) [M]. Beijing: Chemical Industry Press, 1988: 98-100

Simulation of methyl nitrite regeneration in packed column

LI Wenlong, WANG Weihan, LI Zhenhua, MA Xinbin

Key Laboratory for Green Chemical TechnologySchool of Chemical Engineering and TechnologyTianjin UniversityTianjinChina

The process of methyl nitrite (MN) regeneration is a key step in scale-up of synthesis of alkyl esters. To maximize the real atom economy of the process, CO coupling reaction and MN regeneration need to be optimized and proceed at the same reaction rate. Studies on the process and kinetics of MN regeneration are reported. But few mathematical models for MN regeneration in a packed column are available. By analyzing the important reactions of MN regeneration, a mathematic model of MN regeneration in a packed column was developed with the reported kinetic parameters. The effect of main reaction parameters on MN yield was studied. The simulation results showed that MN yield increased with increasing liquid velocity and NO/O2mole ratio but decreased with increasing temperature, gas velocity and N2volume fraction. The simulation results were verified by experiments under different conditions, which proved that this model could be used for the design and operation of MN regeneration. Furthermore, the simulation results at two different mixing modes of NO and O2were presented.

methyl nitrite; packed column; mathematical modeling; regeneration; absorption

2014-08-26.

Prof.MA Xinbin, xbma@tju.edu.cn

10.11949/j.issn.0438-1157.20141300

TQ 018

A

0438—1157（2015）03—0979—08

國家“十一五”科技攻關(guān)項目（2006BAE02B03）。

2014-08-26收到初稿，2014-12-15收到修改稿。

聯(lián)系人：馬新賓。第一作者：李文龍（1989—），男，碩士研究生。

supported by the National Key Project for the 11th Five Year Plan (2006BAE02B03).