• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    數(shù)值模擬和優(yōu)化變壓吸附流程研究進(jìn)展

    2014-04-04 13:09:09銀醇彪張東輝魯東東張正旺
    化工進(jìn)展 2014年3期
    關(guān)鍵詞:變壓數(shù)學(xué)模型穩(wěn)態(tài)

    銀醇彪,張東輝,魯東東,張正旺

    (天津大學(xué)化工學(xué)院化學(xué)工程聯(lián)合國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津 300072)

    變壓吸附技術(shù),是一種基于不同壓力下氣相吸附質(zhì)與固相吸附劑間作用力大小不同的新型氣體分離與凈化技術(shù),它可以通過周期性地變換壓力實(shí)現(xiàn)氣體混合物的分離或提純[1]。1960年 Skarstormsh循環(huán)專利的申請(qǐng),極大地促進(jìn)了變壓吸附技術(shù)在氣體分離領(lǐng)域的研究、應(yīng)用和發(fā)展[2]。目前,變壓吸附技術(shù)已經(jīng)發(fā)展成為一種先進(jìn)的氣體分離技術(shù),其廣泛應(yīng)用于氣體干燥、有機(jī)溶劑回收、空氣分離、氫氣提純、沼氣升級(jí)、煤層氣富甲烷化、二氧化碳捕集、一氧化碳/氫氣分離、乙醇脫水等領(lǐng)域[3]。在變壓吸附技術(shù)應(yīng)用不斷擴(kuò)展的同時(shí),其裝置規(guī)模也不斷擴(kuò)大。目前,全球建成的最大變壓吸附制氫裝置已經(jīng)達(dá)到340 000 m3/h。因此,隨著變壓吸附技術(shù)的不斷推廣和應(yīng)用,變壓吸附裝置規(guī)模的不斷擴(kuò)大,數(shù)值模擬變壓吸附流程顯得尤為重要,這是由于數(shù)值仿真能夠突破實(shí)驗(yàn)裝置的規(guī)模限制,同時(shí)可以節(jié)約?;瘜?shí)驗(yàn)的成本。數(shù)值優(yōu)化變壓吸附流程,是指對(duì)流程的各種工藝參數(shù)進(jìn)行最優(yōu)化設(shè)計(jì),如吸附和解吸壓力、各步驟運(yùn)行時(shí)間、吹掃和置換量、進(jìn)出口流量、閥門開度、吸附塔高徑比等。數(shù)值優(yōu)化可以提升變壓吸附技術(shù)的性能表現(xiàn),如提高目標(biāo)產(chǎn)物的純度和回收率,提升吸附劑的單位產(chǎn)量,降低過程的單位能耗。目前,國(guó)內(nèi)外對(duì)數(shù)值模擬和優(yōu)化變壓吸附流程展開了廣泛的研究,并已取得了顯著的研究成果[4-14]。

    1 變壓吸附流程的數(shù)學(xué)模型

    變壓吸附流程涉及的模型主要有吸附塔、閥門、緩沖罐、真空泵、壓縮機(jī),其中吸附塔模型是整個(gè)流程模型的核心。吸附塔數(shù)學(xué)模型包括質(zhì)量衡算、動(dòng)量衡算、能量衡算、吸附特性方程、質(zhì)量傳遞模型、進(jìn)出口邊界條件、氣體狀態(tài)方程,這些偏微分方程和代數(shù)方程組共同構(gòu)成了吸附塔的數(shù)學(xué)模型。

    動(dòng)量衡算方程有 4種形式:①Darcy’s law equation,將吸附塔壓降描述成與氣體速度成線性常數(shù)正比的關(guān)系;②Karman-Kozeny equation[15],適合于描述層流條件下的吸附塔壓降;③Burke-Plummer equation[15],適合于湍流條件下的吸附塔壓降;④Ergun equation[15],是前兩種形式的結(jié)合,適合于描述層流和湍流混合條件下的吸附塔壓降。Ergun equation是目前描述吸附塔壓降最普遍的形式。能量衡算模型可分為等溫和非等溫模型,數(shù)值模擬和優(yōu)化較小規(guī)模且較小吸附熱的變壓吸附分離過程,為簡(jiǎn)化過程計(jì)算量均采用等溫模型,如小規(guī)??辗诌^程[16-17];數(shù)值模擬和優(yōu)化任意規(guī)模較大吸附熱、吸附熱小但規(guī)模大的變壓吸附分離過程,為保障計(jì)算結(jié)果的準(zhǔn)確性,能量衡算模型必須采用非等溫模型,如變壓吸附捕集二氧化碳[11]、大規(guī)模變壓吸附分離空氣[18]等。此外,對(duì)于忽略吸附塔壁與外部環(huán)境傳熱行為的非等溫模型又可以稱為絕熱非等溫模型[19]。吸附特性方程用于描述多組分氣體競(jìng)爭(zhēng)吸附下各組分的平衡吸附量,它主要受床層溫度、氣體分壓、氣體相互作用等因素影響,其主要包含以下幾種形式:①Henry equation,只有在氣體分壓很低的情況下適用,其作為吸附特性方程實(shí)用性不強(qiáng);②Langmuir equation,是目前應(yīng)用最為廣泛的吸附特性方程,其可以細(xì)分為 single-site Langmuir、dual-site Langmuir[20]、mult-site Langmuir[21]3 種形式,這 3種形式的準(zhǔn)確度和復(fù)雜度依次遞增;③Langmuir-Freundlich equation[22],它通常是在Langmuir equation描述失效的情況下使用,該方程的形式較復(fù)雜,使用時(shí)會(huì)增加數(shù)值模擬和優(yōu)化變壓吸附流程的難度和計(jì)算量;④Exten-Virial equation,是由Taqvi等[23]提出并運(yùn)用于多組分氣體分離,該方程形式非常復(fù)雜,具體的實(shí)際應(yīng)用并不多。質(zhì)量傳遞分為3種模型:①平衡模型,該模型忽略掉氣固兩相間的傳質(zhì)阻力,認(rèn)為氣固相的傳質(zhì)速率很大,但在實(shí)際應(yīng)用中,不管是基于平衡還是基于動(dòng)力學(xué)的吸附劑,氣固相的傳質(zhì)阻力均不能忽略;②線性推動(dòng)力模型,該模型考慮了氣固兩相間的傳質(zhì)阻力,且假設(shè)吸附劑孔隙內(nèi)外的吸附質(zhì)組成一致,是目前應(yīng)用廣泛的質(zhì)量傳遞模型。Sicar等[24]比較研究了線性推動(dòng)力模型和費(fèi)克擴(kuò)散模型在描述吸附塔內(nèi)傳質(zhì)行為的應(yīng)用,得出前者能夠有效地描述吸附劑和吸附質(zhì)間的傳質(zhì)行為;③孔擴(kuò)散模型,包含大孔(大孔和中孔)擴(kuò)散和微孔擴(kuò)散,其考慮了吸附劑空隙內(nèi)外吸附質(zhì)組成相異的情況,又被稱為雙線性推動(dòng)力模型。Rodrigues等[25-27]對(duì)孔擴(kuò)散模型進(jìn)行了完善,使得該模型在描述氣固兩相間的傳質(zhì)行為更準(zhǔn)確。氣體狀態(tài)方程可分為理想和非理想氣體狀態(tài)方程,如果變壓吸附流程的操作壓力處于中低壓范圍,采用理想狀態(tài)方程足以描述氣體吸附質(zhì)的狀態(tài);如果壓吸附流程的操作壓力處于高壓范圍,則采用非理想狀態(tài)方程描述氣體吸附質(zhì)的狀態(tài),如維里方程[28]。吸附塔進(jìn)出口邊界條件是構(gòu)成封閉變壓吸附流程數(shù)學(xué)模型的必要條件。Delgado等[29]比較和分析了數(shù)值模擬變壓吸附流程中各步驟的邊界條件,并重點(diǎn)研究了均壓步驟的邊界條件。

    質(zhì)量、能量衡算模型是吸附塔數(shù)學(xué)模型中最為復(fù)雜的兩項(xiàng),是一組對(duì)于時(shí)間和空間的偏微分方程組,其中空間上可以分為一維、二維、三維這3個(gè)層次。這些方程組的復(fù)雜程度直接取決于吸附塔建模時(shí)所考慮的空間維數(shù)。國(guó)內(nèi)外學(xué)者對(duì)一維、二維的吸附塔數(shù)學(xué)模型都進(jìn)行了研究報(bào)道,一維模型是數(shù)值模擬和優(yōu)化變壓吸附流程普遍采用的數(shù)學(xué)模型。三維模型由于形式過于復(fù)雜,目前鮮有研究報(bào)道。

    1.1 一維數(shù)學(xué)模型

    一維數(shù)學(xué)模型是指建立固定床吸附塔模型時(shí)假定吸附塔內(nèi)流體屬于活塞流,只考慮吸附塔的軸向濃度梯度、壓力梯度和溫度梯度,而忽略吸附塔徑向和吸附劑顆粒內(nèi)外的濃度、壓力、溫度差異。一維模型簡(jiǎn)化了變壓吸附流程的數(shù)學(xué)模型,降低了數(shù)值模擬和優(yōu)化該流程的計(jì)算量、計(jì)算時(shí)間和計(jì)算難度,這對(duì)本已復(fù)雜程度過高、計(jì)算難度過大的固定床吸附塔數(shù)學(xué)模型而言是必要的,因而它是目前國(guó)內(nèi)外學(xué)者在模擬和優(yōu)化變壓吸附流程中普遍采用的數(shù)學(xué)模型。

    1.2 二維數(shù)學(xué)模型

    二維數(shù)學(xué)模型是指在建立吸附塔數(shù)學(xué)模型時(shí)同時(shí)考慮了濃度、壓力和溫度在固定床吸附塔軸向和徑向上的梯度。該模型是對(duì)一維數(shù)學(xué)模型的改進(jìn),提高了描述變壓吸附流程的準(zhǔn)確性,但同時(shí)也較大程度地增加了數(shù)值模擬和優(yōu)化的難度、計(jì)算量和計(jì)算時(shí)間。Kim等[30]研究比較了一維和二維吸附塔模型在數(shù)值模擬變溫吸附中的應(yīng)用,得出了應(yīng)用兩種不同數(shù)學(xué)模型所產(chǎn)生的結(jié)果差別較小,且一維模型計(jì)算結(jié)果也已滿足要求的結(jié)論。Mohamadinejad等[31]比較了一維和二維模型在數(shù)值模擬變壓吸附吸附/解吸二氧化碳中的應(yīng)用,指出通過考慮平均速率和平均孔隙率的影響,能夠提高一維模型的準(zhǔn)確性。Zheng等[32]通過建立吸附塔的二維模型,數(shù)值模擬了簡(jiǎn)單變壓吸附空分流程,給出了固定床吸附塔內(nèi)的濃度場(chǎng)、溫度場(chǎng)和速度場(chǎng)結(jié)果,指出運(yùn)用二維模型模擬變壓吸附流程能夠指導(dǎo)設(shè)計(jì)吸附塔結(jié)構(gòu)。

    1.3 三維模型

    三維模型進(jìn)一步提高了固定床吸附塔模型的準(zhǔn)確性和復(fù)雜性,充分考慮了吸附塔軸向和徑向、吸附劑顆粒徑向這3個(gè)方向的濃度、壓力、溫度變化,從而使吸附塔數(shù)學(xué)模型變成了在三維空間和時(shí)間上的偏微分和代數(shù)方程組。該模型理論上能夠更精確地描述吸附塔內(nèi)的傳質(zhì)、傳熱和傳動(dòng)現(xiàn)象,適宜用于指導(dǎo)設(shè)計(jì)吸附塔結(jié)構(gòu)。然而,由于該模型復(fù)雜程度過高、方程組之間的耦合過大,如果采用它描述固定床吸附塔的數(shù)學(xué)模型,則數(shù)值模擬和優(yōu)化變壓吸附時(shí)對(duì)計(jì)算機(jī)的技術(shù)和性能、求解方法的精度和收斂性等均有著更高的要求。目前,國(guó)內(nèi)外鮮有關(guān)于應(yīng)用三維數(shù)學(xué)模型描述固定床吸附塔的研究報(bào)道。

    2 變壓吸附流程的數(shù)值模擬

    變壓吸附流程的數(shù)學(xué)模型是由關(guān)于空間、時(shí)間的復(fù)雜偏微分方程組和代數(shù)方程組構(gòu)成的混合方程組。目前,求解復(fù)雜解偏微分方程組時(shí)難以給出其解析解,而只能給出它們的數(shù)值解。數(shù)值求解偏微分方程組主要涉及兩方面的選擇,即偏微分方程的離散形式和離散方法。此外,變壓吸附流程屬于周期循環(huán)過程,周期中每個(gè)步驟都有相應(yīng)的進(jìn)出口邊界條件用以構(gòu)成封閉的變壓吸附流程數(shù)學(xué)模型。因此,數(shù)值模擬整個(gè)周期的變壓吸附流程,需要在相鄰兩步驟變換時(shí)切換吸附塔進(jìn)出口邊界條件。

    2.1 偏微分方程的離散形式

    偏微分方程離散的形式分為兩種,即全離散和部分離散。全離散是指對(duì)偏微分方程在空間和時(shí)間上同時(shí)進(jìn)行離散,從而將復(fù)雜的偏微分方程轉(zhuǎn)化為簡(jiǎn)單的代數(shù)方程,進(jìn)而采用代數(shù)方程的求解方法(如牛頓迭代法)求解。部分離散是指只對(duì)偏微分方程在空間上進(jìn)行離散,而時(shí)間上仍保持連續(xù),從而將復(fù)雜的偏微分方程轉(zhuǎn)化為常微分方程,進(jìn)而采用常微分方程的求解方法(如龍格庫(kù)塔積分法)求解。部分離散形式又被稱為線性法,它是目前數(shù)值求解變壓吸附流程模型的常用離散形式。Nilchan等[33]采用全離散形式,數(shù)值模擬了簡(jiǎn)單結(jié)構(gòu)的快速變壓吸附空分流程,取得了較好的模擬結(jié)果。Ko等[34]分別采用全離散和部分離散形式,數(shù)值模擬了復(fù)雜周期結(jié)構(gòu)的變壓吸附捕集二氧化碳流程,結(jié)果表明全離散形式雖然縮短了數(shù)值模擬所需的運(yùn)算時(shí)間,但是得到的數(shù)值結(jié)果不準(zhǔn)確。Biegler等[35]評(píng)價(jià)了全離散和部分離散在數(shù)值模擬變壓吸附流程中的應(yīng)用,指出采用全離散形式應(yīng)用于數(shù)值模擬簡(jiǎn)單結(jié)構(gòu)的變壓吸附流程具有優(yōu)勢(shì),但應(yīng)用于復(fù)雜結(jié)構(gòu)的多塔變壓吸附流程時(shí)對(duì)設(shè)置初值要求高,且運(yùn)算不易收斂和求解結(jié)果可能失真。

    2.2 偏微分方程的離散方法

    離散方法是指采用全離散形式或者部分離散形式,將偏微分方程離散為代數(shù)方程或常微分方程的過程而采取的方法。目前,常用的離散方法包括有限差分法、有限體積法、有限元法和自適應(yīng)多分辨率方法。Santos等[36]采用有限中心差分離散法和部分離散形式,數(shù)值模擬了一種基于新循環(huán)結(jié)構(gòu)的變壓吸附流程,該新循環(huán)的采用顯著地提高了吸附劑的使用效率。周圓圓等[37]采用有限向上差分離散法和部分離散形式,對(duì)三塔真空變壓吸附富氧工藝進(jìn)行了數(shù)值模擬,考察了產(chǎn)品氣流量、產(chǎn)品氣升壓量、吸附塔高徑比等因素對(duì)富氧效果的影響。張湜等[38]采用有限中心差分離散法和部分離散形式,數(shù)值模擬了五塔真空變壓吸附凈化沼氣流程,分析了緩沖罐中雜質(zhì)濃度對(duì)吸附步驟出口濃度的影響。Grande等[26]采用正交有限元離散法和部分離散形式,對(duì)變壓吸附分離丙烷/丙烯流程進(jìn)行了數(shù)值模擬,得到了高純度、高回收率的丙烯產(chǎn)品。Nikolic等[39]基于開放編程軟件 gPROMS建立了能模擬多塔變壓吸附流程的通用平臺(tái),其采用正交有限元離散法和部分離散形式,考察了多塔變壓吸附流程的設(shè)備成本、操作成本、分離性能之間的關(guān)系,并闡述了正交有限元離散法的應(yīng)用優(yōu)勢(shì)。

    Liu等[40]采用正交有限元離散法和部分離散形式,數(shù)值模擬了數(shù)組復(fù)雜周期的多塔真空變壓吸附捕集二氧化碳流程,該研究對(duì)于工業(yè)真空變壓吸附捕集煙道氣二氧化碳有一定的指導(dǎo)意義。Cruz等[41]采用自適應(yīng)多分辨率離散法和部分離散形式,數(shù)值模擬了兩塔的變壓吸附和真空變壓吸附空分流程,求解結(jié)果表明變壓吸附制氧更有優(yōu)勢(shì),指出該離散方法相比差分法和有限元法具有更好的求解穩(wěn)定性和精確度。Jiang等[42]采用有限體積離散法和線性法,數(shù)值模擬了單塔和雙塔真空變壓吸附空氣分離,指出變壓吸附模型中涉及的衡算偏微分方程屬于雙曲線形式的偏微分方程,有限體積離散法相比有限差分離散法能夠更有效地處理這類方程,且有限體積離散法本身也是一種守衡算法。戴先知等[18]采用有限體積離散法和全離散形式,對(duì)變壓吸附空分制氧非等溫過程進(jìn)行數(shù)值模擬,結(jié)合實(shí)驗(yàn)得出了小型裝置中吸附熱對(duì)產(chǎn)品氣的影響不大。

    2.3 吸附塔進(jìn)出口邊界條件的切換

    邊界條件的切換方式有兩種,即機(jī)械指定各步驟進(jìn)出口邊界條件和根據(jù)進(jìn)出口速度方向自主判定。Ko等[43]基于開放編程軟件gPROMS數(shù)值模擬單塔不同循環(huán)變壓吸附/真空變壓吸附捕集二氧化碳流程,機(jī)械指定了原料氣升壓、吸附、順放、逆放、吹掃等步驟對(duì)應(yīng)的吸附塔進(jìn)出口邊界條件。田春剛[44]基于開放編程軟件 MATLAB數(shù)值模擬了兩塔變壓吸附制氮流程,采用機(jī)械指定切換升壓、吸附、降壓、反吹4個(gè)步驟對(duì)應(yīng)的吸附塔進(jìn)出口邊界條件。Nikolic等[39]基于開放編程軟件gPROMS建立了能夠應(yīng)用于多塔變壓吸附數(shù)值模擬的通用結(jié)構(gòu),采用自主判定方法切換吸附周期中各步驟對(duì)應(yīng)的吸附塔進(jìn)出口邊界條件。何東榮等[9]基于專業(yè)的變壓吸附模擬軟件ASPEN-ADSIM,模擬了虛擬三組分的變壓吸附分離過程,其采用的邊界條件切換方式即為自主判定。機(jī)械指定和自主判定這兩種邊界條件的切換方式,在具體的應(yīng)用中各有優(yōu)劣。前者適用于包含空白步驟的變壓吸附循環(huán)結(jié)構(gòu),這是由于空白步驟中吸附塔進(jìn)出口速度方向會(huì)一直發(fā)生微小的變化,這對(duì)于基于進(jìn)出口速度方向而切換邊界條件的自主判定法會(huì)比較困難,而機(jī)械指定則能夠有效地應(yīng)對(duì)這種微小變化。然而,對(duì)于相鄰步驟吸附塔進(jìn)出口速度方向發(fā)生驟變的情況,如先順放后抽真空,機(jī)械指定進(jìn)出口邊界條件則可能會(huì)造成一定的計(jì)算誤差,而自行判定則可以有效地應(yīng)對(duì)。因此,如果將機(jī)械指定和自行判定兩種切換方式結(jié)合使用,則能夠更有效地切換數(shù)值模擬變壓吸附流程中相鄰步驟的進(jìn)出口邊界條件。

    3 變壓吸附流程的數(shù)值優(yōu)化

    數(shù)值優(yōu)化變壓吸附流程能夠提高該分離過程的性能表現(xiàn),如提高產(chǎn)品氣的回收率、純度和單位產(chǎn)率,降低生產(chǎn)過程的單位能耗。目前,國(guó)外關(guān)于變壓吸附流程的數(shù)值優(yōu)化已經(jīng)有了廣泛的研究報(bào)道,且這些研究已應(yīng)用于變壓吸附制氧、制氫和二氧化碳捕集等方面,取得了顯著的成效。Jiang等[42]在對(duì)單塔六步驟變壓吸附制氧流程進(jìn)行了數(shù)值優(yōu)化,在產(chǎn)品氣中氧氣的純度達(dá)到了設(shè)計(jì)值(95%),其制氧的單位能耗下降了 4.87%~7.95%。在隨后的研究中,Jiang等[4]對(duì)五塔十一步驟變壓吸附分離氫氣/氮?dú)?二氧化碳/甲烷/一氧化碳進(jìn)行數(shù)值優(yōu)化,在滿足產(chǎn)品氣中一氧化碳濃度不高于 10 mg/L的條件下,氫氣的回收率相比設(shè)計(jì)值提高了 2.22%~3.32%。Agarwal等[45]基于超結(jié)構(gòu)最優(yōu)化思想,數(shù)值優(yōu)化了變壓吸附捕集煙道氣中二氧化碳的流程,得到了最優(yōu)的兩塔六步驟真空變壓吸附循環(huán)結(jié)構(gòu),并在滿足二氧化碳純度和回收率分別為90%、85%的設(shè)計(jì)條件下,實(shí)現(xiàn)了最優(yōu)捕集二氧化碳的單位能耗為464.76 kW·h/t,遠(yuǎn)低于目前研究報(bào)道的捕集二氧化碳的單位能耗值(765~950 kW·h/t)。

    3.1 循環(huán)穩(wěn)態(tài)定義方法

    變壓吸附流程,其本質(zhì)是一個(gè)包含多個(gè)步驟的周期性動(dòng)態(tài)過程,因而它不具有真正意義的穩(wěn)定狀態(tài)。對(duì)于任意兩個(gè)相鄰的運(yùn)行周期,如果吸附塔內(nèi)的溫度、濃度、壓力等變量的變化相同,則認(rèn)為該變壓吸附流程運(yùn)行達(dá)到了循環(huán)穩(wěn)態(tài)。循環(huán)穩(wěn)態(tài)對(duì)于任何一個(gè)變壓吸附過程均至關(guān)重要,能反映這個(gè)變壓吸附流程運(yùn)行的真正狀態(tài)。對(duì)變壓吸附流程進(jìn)行數(shù)值優(yōu)化,其最優(yōu)化結(jié)果只有出現(xiàn)在循環(huán)穩(wěn)態(tài)或者近似循環(huán)穩(wěn)態(tài)下才有實(shí)際意義。因此,數(shù)值優(yōu)化變壓吸附流程時(shí)需要對(duì)循環(huán)穩(wěn)態(tài)進(jìn)行定義,使最優(yōu)化求解結(jié)果和循環(huán)穩(wěn)態(tài)同時(shí)達(dá)到。

    目前,國(guó)內(nèi)外研究報(bào)道的循環(huán)穩(wěn)態(tài)定義方法主要有4類,即循環(huán)迭代法、直接決定法、全離散法、單塔結(jié)構(gòu)法。循環(huán)迭代法是一種傳統(tǒng)的循環(huán)穩(wěn)態(tài)定義方法,是指在數(shù)值模擬和優(yōu)化變壓吸附流程時(shí),以前個(gè)周期運(yùn)算的結(jié)束值作為后個(gè)周期運(yùn)算的初始值,如此依次進(jìn)行迭代運(yùn)算直至滿足循環(huán)穩(wěn)態(tài)條件。Cruz等[41]在對(duì)變壓吸附流程進(jìn)行數(shù)值優(yōu)化時(shí),采用循環(huán)迭代法對(duì)循環(huán)穩(wěn)態(tài)進(jìn)行限定,限定條件是產(chǎn)品氣的純度和回收率在任意相鄰兩個(gè)周期內(nèi)變化不超過0.01%。Huang等[46]在數(shù)值優(yōu)化變壓吸附富集等離子體反應(yīng)氣中的氫氣時(shí),其采用循環(huán)迭代法對(duì)循環(huán)穩(wěn)態(tài)進(jìn)行限定,限定條件是氣相組成和固相濃度的數(shù)值在任意相鄰兩個(gè)周期內(nèi)變化不超過 0.001。Khajuria等[14]在數(shù)值優(yōu)化變壓吸附分離氫氣和甲烷時(shí),分別以產(chǎn)品氣純度和回收率變化、氣固相濃度變化、溫度變化為限定條件,通過大量的數(shù)值模擬實(shí)驗(yàn)來(lái)考察采用循環(huán)迭代法達(dá)到近似循環(huán)穩(wěn)態(tài)所需要的周期數(shù)目。

    直接決定法是一種基于 Newton算法和精確靈敏度分析的循環(huán)穩(wěn)態(tài)定義方法,它是針對(duì)應(yīng)用循環(huán)迭代法達(dá)到循環(huán)穩(wěn)態(tài)所需計(jì)算時(shí)間長(zhǎng)、運(yùn)算周期多而提出的一種改進(jìn)方法。通過 Newton算法和精確靈敏度分析,直接決定法能夠不斷升級(jí)和優(yōu)化每個(gè)周期運(yùn)算的初始條件,這有別于循環(huán)迭代法中繼承式的初始條件,從而有效地減少了達(dá)到循環(huán)穩(wěn)態(tài)所需的運(yùn)算時(shí)間和周期。 Croft等[47]運(yùn)用直接決定法來(lái)模擬周期變壓吸附的循環(huán)穩(wěn)態(tài),并對(duì)這種方法的穩(wěn)定性進(jìn)行了研究,指出該法相比傳統(tǒng)的循環(huán)迭代法計(jì)算效率更高。Ding等[48]采用直接決定法對(duì)循環(huán)穩(wěn)態(tài)進(jìn)行限定,通過對(duì)周期變壓吸附流程進(jìn)行數(shù)值優(yōu)化,得出了直接決定法能夠有效地定義循環(huán)穩(wěn)態(tài)的結(jié)論。Jiang等[4,42]采用直接決定法對(duì)循環(huán)穩(wěn)態(tài)進(jìn)行限定,數(shù)值優(yōu)化了單塔和多塔變壓吸附流程,計(jì)算結(jié)果表明該法能有效地減少運(yùn)算時(shí)間和周期。Glover等[49]對(duì)影響吸附塔穿透行為的相關(guān)參數(shù)進(jìn)行靈敏度分析時(shí),采用了直接決定法對(duì)循環(huán)穩(wěn)態(tài)進(jìn)行限定。

    全離散法是指對(duì)變壓吸附數(shù)學(xué)模型在空間和時(shí)間上同時(shí)離散,從而將偏微分方程轉(zhuǎn)化為代數(shù)方程進(jìn)行求解,與此同時(shí),循環(huán)穩(wěn)態(tài)的限定條件也被轉(zhuǎn)化為代數(shù)方程。因此,全離散法不是嚴(yán)格的循環(huán)穩(wěn)態(tài)定義方法,因?yàn)樗鼘⒆儔何搅鞒痰膭?dòng)態(tài)模型轉(zhuǎn)化為穩(wěn)態(tài)模型,以致求解只需一次迭代便可達(dá)到循環(huán)穩(wěn)態(tài)。Nilchan等[33]采用全離散法,結(jié)合非線性規(guī)劃對(duì)單塔變壓吸附制氧進(jìn)行了數(shù)值優(yōu)化。Agarwal等[45,50]采用全離散法和精確靈敏度分析,對(duì)變壓吸附捕集二氧化碳的流程進(jìn)行了數(shù)值優(yōu)化。單塔結(jié)構(gòu)法是指假定變壓吸附流程中所有吸附塔運(yùn)行的周期結(jié)構(gòu)完全相同,通過數(shù)據(jù)儲(chǔ)存單塔運(yùn)算中每個(gè)步驟的計(jì)算結(jié)果,待下個(gè)周期運(yùn)算開始時(shí),依次將上個(gè)周期對(duì)應(yīng)步驟的計(jì)算結(jié)果作為初值賦給該周期的對(duì)應(yīng)步驟進(jìn)行迭代運(yùn)算。因此,單塔結(jié)構(gòu)法也不是嚴(yán)格的循環(huán)穩(wěn)態(tài)定義方法,它只是由于通過假設(shè)縮短了運(yùn)算達(dá)到循環(huán)穩(wěn)態(tài)所需的時(shí)間和周期,但還需另外對(duì)循環(huán)穩(wěn)態(tài)進(jìn)行定義。Kumar等[51]首次描述了單塔結(jié)構(gòu)法,指出該法相比傳統(tǒng)的多塔方法具有計(jì)算時(shí)間短、計(jì)算周期少。Jiang等[4]在數(shù)值模擬和優(yōu)化多塔變壓吸附流程時(shí),對(duì)比了真實(shí)的多塔結(jié)構(gòu)和假設(shè)的單塔結(jié)構(gòu),認(rèn)為后者雖然能預(yù)測(cè)最終的循環(huán)穩(wěn)態(tài),但由于其假設(shè)所有塔運(yùn)行的周期結(jié)構(gòu)完全相同,導(dǎo)致該模型對(duì)變壓吸附流程的中間過程描述不準(zhǔn)。

    直接決定法、全離散法和單塔結(jié)構(gòu)法,它們相對(duì)于傳統(tǒng)的循環(huán)迭代法都屬于加速循環(huán)穩(wěn)態(tài)運(yùn)算的方法,因而一定程度地克服了循環(huán)迭代法計(jì)算量大、運(yùn)算周期多的缺點(diǎn)。然而,在這3種改進(jìn)方法中,直接決定法雖也能預(yù)測(cè)達(dá)到循環(huán)穩(wěn)態(tài)前的部分中間狀態(tài)變化,但是由于其不斷升級(jí)周期運(yùn)算的初始條件,導(dǎo)致其不能完整地預(yù)測(cè)整個(gè)變壓吸附流程中間狀態(tài)的變化;全離散法則完全不能預(yù)測(cè)中間狀態(tài)變化,且其預(yù)測(cè)地循環(huán)穩(wěn)態(tài)不一定準(zhǔn)確;單塔結(jié)構(gòu)法雖能預(yù)測(cè)循環(huán)穩(wěn)態(tài),但周期結(jié)構(gòu)一致的假設(shè)導(dǎo)致它對(duì)中間狀態(tài)的預(yù)測(cè)失真。此外,循環(huán)迭代法相比這3種改進(jìn)方法具有更好的運(yùn)算穩(wěn)定性、收斂性。

    3.2 數(shù)值優(yōu)化方法

    最優(yōu)化方法是數(shù)值優(yōu)化變壓吸附流程中的另一重要因素。相比傳統(tǒng)的穩(wěn)態(tài)分離手段,如精餾分離和吸收分離,變壓吸附是一個(gè)周期性動(dòng)態(tài)分離過程,因而對(duì)變壓吸附流程進(jìn)行數(shù)值優(yōu)化時(shí)需要采取更為復(fù)雜的動(dòng)態(tài)優(yōu)化方法。Smith等[52-53]在數(shù)值優(yōu)化變壓吸附流程領(lǐng)域作了開創(chuàng)性的工作,先后應(yīng)用混合整數(shù)線性規(guī)劃(MILP)和混合整數(shù)非線性規(guī)劃(MINLP)對(duì)變壓吸附流程的操作成本進(jìn)行數(shù)值優(yōu)化。Nilchan等[33]在對(duì)偏微分方程全離散成非線性代數(shù)方程的條件下,采用非線性規(guī)劃(NLP)對(duì)變壓吸附分離空氣的平均能耗進(jìn)行數(shù)值優(yōu)化。Knaebel等[54]基于開放編程軟件gPROMS,結(jié)合循環(huán)迭代方法對(duì)循環(huán)穩(wěn)態(tài)進(jìn)行定義,采用簡(jiǎn)化空間的序貫二次規(guī)劃法(r-SQP),在確保產(chǎn)品氣中氫氣純度不低于99.5%的前提下,數(shù)值優(yōu)化變壓吸附分離氫氣/甲烷中氫氣的回收率。Jiang等[42]提出一種基于精確靈敏度分析和Newton算法的r-SQP優(yōu)化算法,同時(shí)結(jié)合置信區(qū)間和比例縮放因子,在指定能耗的前提下對(duì)單塔變壓吸附分離空氣的氧氣回收率進(jìn)行數(shù)值優(yōu)化。在隨后的研究中,他們[4]對(duì)計(jì)算靈敏度的順序進(jìn)行了改進(jìn),使之前的滯后靈敏度分析改進(jìn)為與模型方程并行求解的靈敏度分析,大幅度地降低了數(shù)值優(yōu)化變壓吸附流程中所需的迭代計(jì)算時(shí)間。

    Cruz等[41]基于自適應(yīng)多分辨率離散方法,利用序貫二次規(guī)劃(SQP)對(duì)變壓吸附制氧過程進(jìn)行數(shù)值優(yōu)化。Nikolic等[55]基于開放編程軟件 gPROMS建立了能夠數(shù)值優(yōu)化多塔變壓吸附流程的通用平臺(tái),結(jié)合非精確靈敏度分析和單塔結(jié)構(gòu)法,運(yùn)用r-SQP實(shí)現(xiàn)了對(duì)變壓吸附流程的塔數(shù)目進(jìn)行數(shù)值優(yōu)化,這是一種先進(jìn)的優(yōu)化思想。Ayoub等[56]提出一種混合動(dòng)態(tài)優(yōu)化思想,并比較了滯后算靈敏度計(jì)算和并行靈敏度計(jì)算的應(yīng)用效果,結(jié)果表明基于并行靈敏度計(jì)算的方法能夠獲得更優(yōu)的計(jì)算結(jié)果和更快的計(jì)算速度。Agarwal等[45,50]采用全離散形式,結(jié)合靈敏度分析和 r-SQP,并對(duì)代數(shù)方程組的求解算法進(jìn)行改進(jìn)以避免采用全離散導(dǎo)致失真的求解結(jié)果,提出了一種基于超結(jié)構(gòu)變壓吸附流程的數(shù)值優(yōu)化思想,并將這種思想應(yīng)用于變壓吸附捕集二氧化碳流程的數(shù)值優(yōu)化。與前述數(shù)值優(yōu)化不同,基于超結(jié)構(gòu)思想的數(shù)值優(yōu)化不需要設(shè)定變壓吸附流程的循環(huán)結(jié)構(gòu),而是通過最優(yōu)化結(jié)果來(lái)確定。Dowling等[13]基于超結(jié)構(gòu)思想,采用更為合理的部分離散形式,提出一種適用于大規(guī)模變壓吸附流程數(shù)值優(yōu)化的優(yōu)化策略。Khajuria等[14]基于開放編程軟件gPROMS,利用非精確靈敏度分析和 r-SQP,結(jié)合循環(huán)迭代法定義循環(huán)穩(wěn)態(tài),對(duì)原料氣組成變動(dòng)的變壓吸附流程的設(shè)計(jì)和控制同時(shí)進(jìn)行數(shù)值優(yōu)化,并比較了模型預(yù)測(cè)控制(MPC)和微分積分控制(PI)的應(yīng)用效果。此外,Sankararao等[6]采用多目標(biāo)優(yōu)化方法對(duì)兩塔變壓吸附制氧流程進(jìn)行了數(shù)值優(yōu)化,他們共執(zhí)行了兩組雙目標(biāo)優(yōu)化和一組四目標(biāo)優(yōu)化。Haghpanah等[57]采用多目標(biāo)優(yōu)化方法對(duì)單塔變壓吸附捕集二氧化碳流程進(jìn)行數(shù)值優(yōu)化,求得了產(chǎn)品純度和回收率、單位能耗和單位產(chǎn)率的帕累托解。

    3.3 變壓吸附優(yōu)化的商業(yè)軟件簡(jiǎn)介

    變壓吸附以其模型方程耦合程度高、模型方程數(shù)量大、模型變量數(shù)量級(jí)差別大等特點(diǎn),對(duì)它進(jìn)行數(shù)值優(yōu)化仍存在較大的困難。目前,據(jù)報(bào)道已經(jīng)實(shí)現(xiàn)商業(yè)的變壓吸附優(yōu)化軟件主要有英國(guó)PSE公司的gPROMS軟件和美國(guó)Aspen公司的Aspen Adsim軟件。Khajuria[14]、Nikolic[55]等采用gPROMS研究了多塔變壓吸附分離多組分氣體流程的數(shù)值優(yōu)化,提高了該變壓吸附分離過程中產(chǎn)品的回收率、純度等表現(xiàn)變量。美國(guó)Aspen公司在用戶使用說明書中簡(jiǎn)要介紹了Aspen Adsim數(shù)值優(yōu)化變壓吸附的功能,然而鮮有文獻(xiàn)報(bào)道該軟件在變壓吸附優(yōu)化方面的研究實(shí)例。此外,美國(guó)卡耐基梅隆大學(xué)Bigler教授課題組[13,42,45,50]在其公開發(fā)表的關(guān)于變壓吸附優(yōu)化的研究論文中,對(duì)他們所采用的優(yōu)化程序代碼進(jìn)行了簡(jiǎn)單介紹,但并未說明其具體的商業(yè)應(yīng)用情況。

    4 數(shù)值模擬和優(yōu)化面臨的問題

    然而,數(shù)值模擬和優(yōu)化變壓吸附流程的研究在取得顯著進(jìn)展的同時(shí),也存在著以下3點(diǎn)可能限制其發(fā)展和推廣的問題。①變壓吸附的理論研究落后于它的實(shí)際應(yīng)用。變壓吸附技術(shù)已經(jīng)成為一項(xiàng)主要的氣體分離技術(shù),雖然國(guó)內(nèi)外關(guān)于變壓吸附吸附單組分氣體或分離雙組分氣體的實(shí)驗(yàn)研究較多,但是對(duì)多組分氣體分離的實(shí)驗(yàn)研究相當(dāng)缺乏,而多組分氣體分離的相關(guān)數(shù)據(jù)又難以通過目前已有的吸附理論模型獲取,這導(dǎo)致設(shè)計(jì)分離多組分氣體的變壓吸附流程時(shí)缺乏重要的基礎(chǔ)數(shù)據(jù),如吸附平衡數(shù)據(jù)、動(dòng)力學(xué)數(shù)據(jù)、吸附熱等。重要基礎(chǔ)數(shù)據(jù)的缺失,會(huì)降低數(shù)值模擬和優(yōu)化結(jié)果的可信度。②變壓吸附流程數(shù)學(xué)模型的準(zhǔn)確度。相比精餾或吸收等傳統(tǒng)分離手段,雖然變壓吸附流程的設(shè)備更簡(jiǎn)單、操作更容易實(shí)現(xiàn)自動(dòng)化,但是它的數(shù)學(xué)模型遠(yuǎn)遠(yuǎn)較前二者復(fù)雜繁瑣。因此,為了有效地降低數(shù)值模擬和優(yōu)化變壓吸附流程的計(jì)算量和計(jì)算時(shí)間,國(guó)內(nèi)外研究普遍采用一維活塞流模型來(lái)描述變壓吸附流程,并作了一些其他方面的假設(shè),然而這些簡(jiǎn)化都會(huì)降低模型預(yù)測(cè)結(jié)果的準(zhǔn)確性。③計(jì)算機(jī)性能和技術(shù)、求解算法的限制。建立準(zhǔn)確度高的二維或三維數(shù)學(xué)模型來(lái)描述變壓吸附流程,雖然能夠提高數(shù)值模擬和優(yōu)化結(jié)果的可信度,但同時(shí)也大幅度提高了待求解數(shù)學(xué)模型的復(fù)雜程度,使得數(shù)值模擬和優(yōu)化的求解難度大幅度提高。求解高維數(shù)的變壓吸附數(shù)學(xué)模型,這對(duì)計(jì)算機(jī)技術(shù)和性能,求解算法的收斂性、穩(wěn)定性和精準(zhǔn)度均有著更高的要求。

    5 結(jié) 語(yǔ)

    數(shù)值模擬和優(yōu)化變壓吸附流程具有三方面的用途,即節(jié)約模化實(shí)驗(yàn)的成本、指導(dǎo)工藝流程的設(shè)計(jì)、提升工藝流程的性能表現(xiàn)。隨著吸附理論研究的深入開展、變壓吸附分離多組分氣體實(shí)驗(yàn)研究的擴(kuò)展、變壓吸附流程數(shù)學(xué)模型的完善、計(jì)算機(jī)性能和技術(shù)的提高、數(shù)值求解算法的改進(jìn),數(shù)值模擬和優(yōu)化變壓吸附流程的研究必將更為深入,進(jìn)而使模擬和優(yōu)化的數(shù)值結(jié)果能夠更好地應(yīng)用于指導(dǎo)和提升變壓吸附工藝的設(shè)計(jì),促進(jìn)變壓吸附技術(shù)更廣泛的應(yīng)用。

    [1]Ruthven D M,F(xiàn)arooq S,Knaebel K S.Pressure Swing Adsorption[M].New York:VCH Publishers Press,1994.

    [2]Skarstrom C W.Method and apparatus for fractionating gaseous mixtures by adsorption:US,2944627[P].1960-07-12.

    [3]Sircar S.Pressure swing adsorption[J].Industrial and Engineering Chemistry Research,2002,41:1389-1392.

    [4]Jiang L,F(xiàn)ox V G,Biegler L T.Simulation and optimal design of multiple-bed pressure swing adsorption systems[J].AIChE Journal,2004,50: 2904-2917.

    [5]Kostroski K P,Wankat P C.High recovery cycles for gas separations by pressure-swing adsorption[J].Industrial and Engineering Chemistry Research,2006,45:8117-8133.

    [6]Sankararao S,Gupta S K.Multi-objective optimization of pressure swing adsorption for air separation[J].Industrial and Engineering Chemistry Research,2007,46:3751-3765.

    [7]Reynolds S P,Mehrotra A,Ebner A D,et al.Heavy reflux PSA cycles for CO2recovery from flue gas:Part Ⅰ.Performance evaluation[J].Adsorption,2008,14:399-413.

    [8]Agarwal A,Biegler L T,Zitney S E.Simulation and optimization of pressure swing adsorption systems using reduced-order modeling[J].Industrial and Engineering Chemistry Research,2009,48:2327-2343.

    [9]何東榮,周向輝,張東輝.利用ASPEN-ADSIM模擬變壓吸附分離過程[J].天然氣化工,2009,34(3):11-15.

    [10]Santos M P S,Grande C A,Rodrigues A E.Pressure swing adsorption for biogas upgrading.Effect of recycling streams in pressure swing adsorption design[J].Adsorption,2011,50:974-985.

    [11]沈春枝,孫玉柱,李平,等.真空變壓吸附過程捕獲煙道氣中CO2的數(shù)值模擬[J].華東理工大學(xué)學(xué)報(bào):自然科學(xué)版,2011,37(5):524-531.

    [12]周圓圓,楊華偉,張東輝.甲烷/氮?dú)庾儔何椒蛛x的實(shí)驗(yàn)與模擬[J].天然氣化工,2011,36(5):21-27.

    [13]Dowling A W,Vetukuri S R R,Biegler L T.Large-scale optimization strategies for pressure swing adsorption cycles synthesis[J].AIChE Journal,2012,58(12):3777-3791.

    [14]Khajuria H,Pistikopoulos E N.Optimization and control of pressure swing adsorption process under uncertainly[J].AIChE Journal,2013,59(1):120-131.

    [15]Bird R B,Stewart W E,Lightfoot E N.Transport Phenomena[M].New York:John Wiley & Sons Press,1960.

    [16]Sivakumar S V,Rao D P.Modified duplex psa.1.sharp separation and process intensification for CO2-N2-13X zeolite system[J].Industrial and Engineering Chemistry Research,2011,50:3426-3436.

    [17]Sivakumar S V,Rao D P.Modified duplex PSA.2.Sharp separation and process intensification for CO2-N2-5A zeolite system[J].Industrial and Engineering Chemistry Research,2011,50:3437-3445.

    [18]戴先知,劉應(yīng)書.變壓吸附空分制氧過程非等溫模擬[J].低溫工程,2006(4):47-53.

    [19]王嘯,馬正飛,姚虎卿.真空變壓吸附空分制氧等溫與非等溫過程模擬比較[J].南京工業(yè)大學(xué)學(xué)報(bào),2003,25(4):49-54.

    [20]Yang R T.Gas Separation by Adsorption Process[M].London:Imperial College Press,1997.

    [21]Nitta T,Shigetomi T,Kuro-Oka M,et al.An adsorption isotherm of multi-site occupancy model for homogeneous surface[J].Journal of Chemical Engineering of Japan,1984,17:994-1003.

    [22]Yang R T.Gas Separation by Adsorption Process[M].Boston:Butter-worths Press,1987.

    [23]Taqvi S M,LeVan M D.Virial description of two-component adsorption on homogeneous and heterogeneous surfaces[J].Industrial and Engineering Chemistry Research,1997,36:2197-2206.

    [24]Sircar S,Hufton J R.Why does the linear driving force model for adsorption kinetics work?[J].Adsorption,2000,6:137-147.

    [25]Azevedo D C S,Rodrigues A E.Bilinear driving force approximation in the modeling of a simulated moving bed using bidisperse adsorbents[J].Industrial and Engineering Chemistry Research,1999,38:3519-3529.

    [26]Grande C A,Rodrigues A E.Propane/propylene separation by pressure swing adsorption using zeolite 4A[J].Industrial and Engineering Chemistry Research,2005,44:8815-8829.

    [27]Gavenati S,Grande C A,Rodrigues A E.Removal of carbon dioxide from natural gas by vacuum pressure swing adsorption[J].Energy &Fuels,2006,20:2648-2659.

    [28]Grande C A,Blom R,M?ller A,et al.High-pressure separation of CH4/CO2using activated carbon[J].Chemical Engineering Science,2013,89:10-20.

    [29]Delgado J A,Rodrigues A E.Analysis of the boundary conditions for the simulation of the pressure equalization step in PSA cycles[J].Chemical Engineering Science,2008,63:4452-4463.

    [30]Kim Min-Bae,Moon Jong-Ho,Lee Chang-Ha,et al.Effects of heats transfer on the transient dynamics of temperature swing adsorption process[J].Korean Journal of Chemical Engineering,2004,21(3):703-711.

    [31]Mohamadinejad H,Knox J C,Smith J E.Experimental and numerical investigation of two-dimensional CO2adsorption/desorption in packed sorption beds under non-ideal flows[J].Separation Science and Technology,2003,38(16):3875-3904.

    [32]Zheng Xingang,Liu Yingshu,Liu Wenhai.Two-dimensional modeling of the transport phenomena in the adsorber during pressure swing adsorption process[J].Industrial and Engineering Chemistry Research,2010,49:11814-11824.

    [33]Nilchan S,Pantelides C C.On the optimization of periodic adsorption process[J].Adsorption,1998,4:113-147.

    [34]Ko D,Siriwardane R,Biegler L T.Optimization of a pressure-swing adsorption process using zeolite 13X for CO2sequestration[J].Industrial and Engineering Chemistry Research,2003,42:339-348.

    [35]Biegler L T,Jiang L,F(xiàn)ox V G.Recent advances in simulation and optimal design of pressure swing adsorption systems[J].Separations and Purification Reviews,2004,33(1):1-39.

    [36]Santos M S,Grande C A,Rodrigues A E.New cycle configuration to enhance performance of kinetic PSA processes[J].Chemical Engineering Science,2011,66:1590-1599.

    [37]周圓圓,韋向攀,張東輝.三塔真空變壓吸附富氧工藝過程模擬[J].化工進(jìn)展,2011,30(S1):263-267.

    [38]張湜,陳文亮,李暉,等.真空變壓吸附沼氣凈化過程的仿真研究[J].化工學(xué)報(bào),2013,64(4):1300-1305.

    [39]Nikolic D,Giovanoglou A,Georgiadis M C,et al.Generic modeling framework for gas separations using multibed pressure swing adsorption process[J].Industrial and Engineering Chemistry Research,2008,47:3156-3169.

    [40]Liu Z,Grande C A,Li P,et al.Multi-bed vacuum pressure swing adsorption for carbon dioxide capture from flue gas[J].Separation and Purification Technology,2011,81:307-317.

    [41]Cruz P,Santos J C,Mendes A,et al.Cyclic adsorption separation processes:Analysis strategy and optimization procedure[J].Chemical Engineering Science,2003,58:3143-3158.

    [42]Jiang L,Biegler L T.Simulation and optimization of pressure-swing adsorption systems for air separation[J].AIChE Journal,2003,49(5):1140-1157.

    [43]Ko D,Siriwardane R,Biegler L T.Optimization of pressure swing adsorption and fractionated vacuum pressure swing adsorption process for CO2capture[J].Industrial and Engineering Chemistry Research,2005,44:8084-8094.

    [44]田春剛.變壓吸附制氮?dú)獾臄?shù)值模擬[D].武漢:華中科技大學(xué),2007.

    [45]Agarwal A,Biegler L T.A superstructure-based optimal synthesis of psa cycles for post-combustion CO2capture[J].AIChE Journal,2010,56(7):1813-1828.

    [46]Huang Q,Malekian A,Eic M.Optimization of PSA process for producing enriched hydrogen from plasma reactor gas[J].Separation and Purification Technology,2008,62:22-32.

    [47]Croft D T,Le van M D.Periodic states of adsorption cycles-Ⅰ.Direct determination and stability[J].Chemical Engineering Science,1994,49(11):1821-1829.

    [48]Ding Y,Croft D T,LeVan M D.Periodic states of adsorption cyclesⅣ.Direct optimization[J].Chemical Engineering Science,2002,57(21):4521-4531.

    [49]Glover T G,LeVan M D.Sensitivity analysis of adsorption bed behavior:Examination of pulse inputs and layered-bed optimization[J].Chemical Engineering Science,2008,63:2086-2098.

    [50]Agarwal A,Biegler L T,Zitney S E.Superstructure-based optimal synthesis of pressure swing adsorption cycles for precombustion CO2capture[J].Industrial and Engineering Chemistry Research,2010,49:5066-5079.

    [51]Kumar R,F(xiàn)ox V G,Hartzog D G,et al.A versatile process simulator for adsorptive separations[J].Chemical Engineering Science,1994,49:3115-3125.

    [52]Smith O J,Westerberg A W.Mixed-integer programming for pressure swing adsorption cycle scheduling[J].Chemical Engineering Science,1990,45(9):2833-2842.

    [53]Smith O J,Westerberg A W.The optimal design of pressure swing adsorption systems[J].Chemical Engineering Science,1991,46(12):2967-2976.

    [54]Knaebel S P,Ko D,Biegler L T.Simulation and optimization of a pressure swing adsorption system:Recovering hydrogen from methane[J].Adsorption,2005,11:615-620.

    [55]Nikolic D,Kikkinides E S,Georgiadis M C.Optimization of multibed pressure swing adsorption process[J].Industrial and Engineering Chemistry Research,2009,48:5388-5398.

    [56]Ayoub S,Tondeur D,Latifi M A.New formulation of optimizationbased simulation of a pressure swing adsorption process:Hybrid dynamic optimization[J].AIChE Journal,2011,57(5):1367-1372.

    [57]Haghpanah R,Majumder A,Nilam R,et al.Multiobjective optimization of a four-step adsorption process for postcombustion CO2captureviafinite volume simulation[J].Industrial and Engineering Chemistry Research,2013,52(11):4249-4265.

    猜你喜歡
    變壓數(shù)學(xué)模型穩(wěn)態(tài)
    基于變壓吸附分離技術(shù)在氣體中的應(yīng)用研究
    可變速抽水蓄能機(jī)組穩(wěn)態(tài)運(yùn)行特性研究
    AHP法短跑數(shù)學(xué)模型分析
    活用數(shù)學(xué)模型,理解排列組合
    一種變頻變壓的雙頻注入絕緣監(jiān)測(cè)方法
    碳化硅復(fù)合包殼穩(wěn)態(tài)應(yīng)力與失效概率分析
    電廠熱力系統(tǒng)穩(wěn)態(tài)仿真軟件開發(fā)
    煤氣與熱力(2021年4期)2021-06-09 06:16:54
    元中期歷史劇對(duì)社會(huì)穩(wěn)態(tài)的皈依與維護(hù)
    中華戲曲(2020年1期)2020-02-12 02:28:18
    對(duì)一個(gè)數(shù)學(xué)模型的思考
    航空用24脈波自耦變壓整流裝置的研究
    АⅤ资源中文在线天堂| 午夜免费男女啪啪视频观看 | 久久香蕉精品热| 国产精品无大码| 久久久久久伊人网av| 国产伦人伦偷精品视频| 亚洲一区高清亚洲精品| 中出人妻视频一区二区| 欧美3d第一页| 成人性生交大片免费视频hd| 国产人妻一区二区三区在| 国产精品美女特级片免费视频播放器| 国产视频内射| 在线国产一区二区在线| 最好的美女福利视频网| 我要搜黄色片| 亚洲精品456在线播放app | 亚洲黑人精品在线| 成熟少妇高潮喷水视频| 亚洲av二区三区四区| 男人舔奶头视频| 色哟哟·www| 春色校园在线视频观看| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩瑟瑟在线播放| 国内揄拍国产精品人妻在线| 久久天躁狠狠躁夜夜2o2o| 色av中文字幕| 色综合色国产| 亚洲美女视频黄频| 午夜福利在线观看吧| 亚洲av电影不卡..在线观看| 亚洲av五月六月丁香网| 99久久久亚洲精品蜜臀av| 22中文网久久字幕| 精品久久久久久久久亚洲 | 性色avwww在线观看| 久久精品国产亚洲av香蕉五月| 亚洲精品色激情综合| 亚洲五月天丁香| 内地一区二区视频在线| 又紧又爽又黄一区二区| 草草在线视频免费看| av视频在线观看入口| 成人二区视频| 欧美人与善性xxx| 99久久久亚洲精品蜜臀av| 国产精品免费一区二区三区在线| 在现免费观看毛片| 真人做人爱边吃奶动态| 精品人妻视频免费看| 中文字幕熟女人妻在线| 99久国产av精品| 不卡视频在线观看欧美| 久久精品国产亚洲av香蕉五月| 毛片一级片免费看久久久久 | 午夜福利在线在线| 亚洲七黄色美女视频| 麻豆av噜噜一区二区三区| 联通29元200g的流量卡| 老熟妇仑乱视频hdxx| 亚洲18禁久久av| 国内精品宾馆在线| 免费av观看视频| 亚洲成人久久性| 九色成人免费人妻av| 麻豆成人午夜福利视频| 在线播放无遮挡| 日本与韩国留学比较| 免费看美女性在线毛片视频| 久久久国产成人精品二区| 久久久久久久久久成人| 国产精品av视频在线免费观看| 久久久久久久精品吃奶| 又黄又爽又刺激的免费视频.| 久久精品91蜜桃| 五月玫瑰六月丁香| 国产视频内射| 一个人看视频在线观看www免费| 午夜福利18| 在线观看美女被高潮喷水网站| 国产69精品久久久久777片| 国产午夜精品论理片| 自拍偷自拍亚洲精品老妇| 亚洲av一区综合| 91精品国产九色| 亚洲欧美日韩东京热| 少妇人妻一区二区三区视频| 亚洲五月天丁香| 最近最新中文字幕大全电影3| 在线天堂最新版资源| 免费不卡的大黄色大毛片视频在线观看 | 成人特级av手机在线观看| 中文字幕精品亚洲无线码一区| 在线国产一区二区在线| 成人三级黄色视频| 69av精品久久久久久| 免费在线观看成人毛片| 精品久久久久久成人av| 色综合色国产| 亚洲精品色激情综合| 别揉我奶头~嗯~啊~动态视频| a级一级毛片免费在线观看| 人人妻,人人澡人人爽秒播| 一级毛片久久久久久久久女| 久久久久精品国产欧美久久久| 国产亚洲精品久久久久久毛片| 久久久久久久久久黄片| 久久精品国产亚洲网站| 国产白丝娇喘喷水9色精品| 国产精品久久久久久精品电影| 午夜老司机福利剧场| 又黄又爽又刺激的免费视频.| 不卡一级毛片| 男女边吃奶边做爰视频| 精品日产1卡2卡| 日韩在线高清观看一区二区三区 | 国产主播在线观看一区二区| 我的女老师完整版在线观看| 一夜夜www| bbb黄色大片| 欧美日本视频| 免费看a级黄色片| 天堂√8在线中文| 欧美激情在线99| 国产麻豆成人av免费视频| 亚洲精华国产精华液的使用体验 | 22中文网久久字幕| 99国产精品一区二区蜜桃av| 成人一区二区视频在线观看| 91午夜精品亚洲一区二区三区 | 亚洲成人免费电影在线观看| 精品乱码久久久久久99久播| 日本a在线网址| 国产高清有码在线观看视频| a级毛片a级免费在线| 亚洲国产精品成人综合色| 亚洲国产色片| 精品人妻一区二区三区麻豆 | 亚洲精品一区av在线观看| 免费观看人在逋| 亚洲男人的天堂狠狠| 一区二区三区激情视频| 亚洲性夜色夜夜综合| 老司机午夜福利在线观看视频| 网址你懂的国产日韩在线| 久久久久久久久久黄片| 国产69精品久久久久777片| 国产白丝娇喘喷水9色精品| 免费不卡的大黄色大毛片视频在线观看 | 三级毛片av免费| 老司机深夜福利视频在线观看| 国内揄拍国产精品人妻在线| 午夜激情欧美在线| 国产老妇女一区| 九色成人免费人妻av| 成年免费大片在线观看| 999久久久精品免费观看国产| 久久久久久久午夜电影| 成熟少妇高潮喷水视频| 亚洲中文字幕一区二区三区有码在线看| aaaaa片日本免费| 91在线观看av| av中文乱码字幕在线| 国产精品日韩av在线免费观看| 日日干狠狠操夜夜爽| 一本一本综合久久| 别揉我奶头 嗯啊视频| 亚洲无线观看免费| 在线观看av片永久免费下载| 日本精品一区二区三区蜜桃| 国产精品美女特级片免费视频播放器| 亚洲经典国产精华液单| 久久久久久大精品| 精品久久久久久久久久久久久| 久久热精品热| 午夜福利在线观看吧| 国产精品野战在线观看| 久久精品国产自在天天线| 毛片女人毛片| 欧美日本视频| 亚洲精品色激情综合| 桃色一区二区三区在线观看| 国产极品精品免费视频能看的| 亚洲内射少妇av| 国产成人aa在线观看| 天堂av国产一区二区熟女人妻| 国产免费一级a男人的天堂| 我的女老师完整版在线观看| 久久久精品欧美日韩精品| 亚洲av免费高清在线观看| 精品久久久久久久久久久久久| 亚洲色图av天堂| 黄色女人牲交| 亚洲一区高清亚洲精品| 身体一侧抽搐| 久久午夜福利片| 赤兔流量卡办理| 日韩精品有码人妻一区| x7x7x7水蜜桃| 999久久久精品免费观看国产| 偷拍熟女少妇极品色| 午夜久久久久精精品| 亚洲av中文av极速乱 | 我的老师免费观看完整版| 制服丝袜大香蕉在线| 观看免费一级毛片| 黄片wwwwww| av女优亚洲男人天堂| 少妇熟女aⅴ在线视频| 亚洲精品色激情综合| 久久精品国产99精品国产亚洲性色| .国产精品久久| 露出奶头的视频| 成年女人毛片免费观看观看9| 一个人看视频在线观看www免费| 日韩中文字幕欧美一区二区| 欧美一区二区亚洲| 中文字幕久久专区| 久久99热这里只有精品18| 精品一区二区三区视频在线| 啦啦啦韩国在线观看视频| 不卡一级毛片| 男女之事视频高清在线观看| 99久久精品一区二区三区| 一a级毛片在线观看| 国产真实乱freesex| 看十八女毛片水多多多| 国产伦人伦偷精品视频| 日日摸夜夜添夜夜添小说| 久久这里只有精品中国| 性欧美人与动物交配| 三级国产精品欧美在线观看| 啪啪无遮挡十八禁网站| 亚洲,欧美,日韩| 深夜精品福利| 精品人妻一区二区三区麻豆 | 国产精品久久久久久精品电影| 久久久久性生活片| 12—13女人毛片做爰片一| 毛片一级片免费看久久久久 | 亚洲18禁久久av| 精品日产1卡2卡| 男女啪啪激烈高潮av片| 男插女下体视频免费在线播放| 一个人看视频在线观看www免费| 国产国拍精品亚洲av在线观看| 老司机福利观看| 97碰自拍视频| 久久久国产成人免费| 色综合色国产| 一个人免费在线观看电影| 人妻丰满熟妇av一区二区三区| 一级a爱片免费观看的视频| 午夜精品久久久久久毛片777| 国产精品98久久久久久宅男小说| 成人午夜高清在线视频| 91狼人影院| 九色成人免费人妻av| 高清在线国产一区| 日本色播在线视频| 色视频www国产| 我的老师免费观看完整版| 69人妻影院| 尤物成人国产欧美一区二区三区| 亚洲无线观看免费| 日日啪夜夜撸| 极品教师在线免费播放| 成人性生交大片免费视频hd| 赤兔流量卡办理| 网址你懂的国产日韩在线| 免费在线观看影片大全网站| av在线天堂中文字幕| 老司机福利观看| 欧美日韩综合久久久久久 | 中国美女看黄片| av专区在线播放| 亚洲精品在线观看二区| 欧美黑人巨大hd| 91麻豆av在线| 亚洲性久久影院| 成人三级黄色视频| 日本 欧美在线| 国产色爽女视频免费观看| 国产午夜福利久久久久久| 亚洲精品粉嫩美女一区| 成人鲁丝片一二三区免费| 国产精品一区二区三区四区免费观看 | 亚洲电影在线观看av| 久久久久久国产a免费观看| 非洲黑人性xxxx精品又粗又长| 在线观看66精品国产| 深爱激情五月婷婷| 欧美高清性xxxxhd video| 最近中文字幕高清免费大全6 | 精品国内亚洲2022精品成人| 欧美日韩精品成人综合77777| 成人av在线播放网站| 亚洲,欧美,日韩| 久久精品影院6| 一级av片app| 国产免费一级a男人的天堂| 啪啪无遮挡十八禁网站| 国产亚洲精品综合一区在线观看| 久久中文看片网| 欧美中文日本在线观看视频| 麻豆成人av在线观看| 亚洲精品国产成人久久av| 国产 一区 欧美 日韩| 熟女电影av网| 久久精品夜夜夜夜夜久久蜜豆| 丝袜美腿在线中文| 麻豆精品久久久久久蜜桃| 一级黄色大片毛片| 久久亚洲精品不卡| 国内毛片毛片毛片毛片毛片| 一卡2卡三卡四卡精品乱码亚洲| 不卡一级毛片| 丰满的人妻完整版| 久久国内精品自在自线图片| 国产精品亚洲美女久久久| 日本 av在线| 国产高清有码在线观看视频| 91麻豆av在线| 99久久久亚洲精品蜜臀av| 我要看日韩黄色一级片| 中文字幕免费在线视频6| av福利片在线观看| 亚洲精品日韩av片在线观看| 九九热线精品视视频播放| 97超级碰碰碰精品色视频在线观看| 国产 一区 欧美 日韩| 三级男女做爰猛烈吃奶摸视频| 久久亚洲精品不卡| 日本免费a在线| 欧美色欧美亚洲另类二区| 在线观看舔阴道视频| 日韩欧美国产一区二区入口| 亚洲精华国产精华精| 欧美国产日韩亚洲一区| 亚洲专区中文字幕在线| 欧美一区二区国产精品久久精品| 在线看三级毛片| 黄色欧美视频在线观看| 国产三级中文精品| 禁无遮挡网站| 精品久久久久久久久久久久久| 少妇高潮的动态图| 99热只有精品国产| 无人区码免费观看不卡| 中国美白少妇内射xxxbb| 成人三级黄色视频| 国产精品嫩草影院av在线观看 | 国产伦精品一区二区三区视频9| 国产中年淑女户外野战色| 在线播放无遮挡| 亚洲av美国av| 亚洲avbb在线观看| av天堂在线播放| 久久久久久国产a免费观看| 高清在线国产一区| 亚洲成人精品中文字幕电影| 日本 欧美在线| 一本精品99久久精品77| 国产高潮美女av| 18禁在线播放成人免费| 国产av一区在线观看免费| 国产精品福利在线免费观看| 黄色视频,在线免费观看| 亚洲在线自拍视频| 精品一区二区三区视频在线观看免费| 99久久中文字幕三级久久日本| x7x7x7水蜜桃| 亚洲成人免费电影在线观看| 国产精品精品国产色婷婷| 精品久久久久久久久久久久久| 国产女主播在线喷水免费视频网站 | 日本成人三级电影网站| 亚洲一级一片aⅴ在线观看| 国产v大片淫在线免费观看| 免费在线观看影片大全网站| 免费人成视频x8x8入口观看| 成年女人毛片免费观看观看9| 精品一区二区三区人妻视频| 男女边吃奶边做爰视频| 亚洲人与动物交配视频| 亚洲性久久影院| 男女做爰动态图高潮gif福利片| 夜夜夜夜夜久久久久| 美女被艹到高潮喷水动态| 久久欧美精品欧美久久欧美| 午夜福利视频1000在线观看| 看黄色毛片网站| 国产 一区精品| 乱系列少妇在线播放| 成人高潮视频无遮挡免费网站| 国产精品免费一区二区三区在线| 亚洲欧美精品综合久久99| 此物有八面人人有两片| 美女cb高潮喷水在线观看| 亚洲色图av天堂| 深夜精品福利| x7x7x7水蜜桃| 国产亚洲精品久久久久久毛片| 女生性感内裤真人,穿戴方法视频| 国产大屁股一区二区在线视频| 久久精品人妻少妇| 亚洲精品影视一区二区三区av| 99久久精品一区二区三区| 99热只有精品国产| 国产v大片淫在线免费观看| 嫩草影视91久久| 午夜福利在线观看吧| 亚洲国产日韩欧美精品在线观看| 99热只有精品国产| 九九爱精品视频在线观看| 日韩人妻高清精品专区| 日本爱情动作片www.在线观看 | 国语自产精品视频在线第100页| 天堂√8在线中文| 欧美国产日韩亚洲一区| 精品国产三级普通话版| 国产精品美女特级片免费视频播放器| 欧美日本视频| 亚洲一级一片aⅴ在线观看| 人妻少妇偷人精品九色| 国产伦精品一区二区三区四那| 日韩欧美国产一区二区入口| 黄色视频,在线免费观看| 小蜜桃在线观看免费完整版高清| 我的女老师完整版在线观看| 日本黄色片子视频| 日韩人妻高清精品专区| 免费搜索国产男女视频| 午夜日韩欧美国产| 日本精品一区二区三区蜜桃| 精品福利观看| 99九九线精品视频在线观看视频| 啦啦啦观看免费观看视频高清| 国产精品久久久久久精品电影| 夜夜夜夜夜久久久久| 久久久久久九九精品二区国产| 色播亚洲综合网| 久久久久久久久久成人| 午夜亚洲福利在线播放| 变态另类成人亚洲欧美熟女| 麻豆成人午夜福利视频| 亚洲自拍偷在线| 亚洲国产精品久久男人天堂| 五月玫瑰六月丁香| 日日摸夜夜添夜夜添av毛片 | 日韩中文字幕欧美一区二区| 国产美女午夜福利| 亚洲国产精品sss在线观看| 校园春色视频在线观看| 国产精品久久视频播放| 亚洲人成网站在线播放欧美日韩| 乱系列少妇在线播放| 婷婷精品国产亚洲av在线| 国产v大片淫在线免费观看| 久久久成人免费电影| 国产色婷婷99| 日韩欧美免费精品| 又紧又爽又黄一区二区| 国产真实乱freesex| 午夜福利18| 欧美国产日韩亚洲一区| 精品久久国产蜜桃| 亚洲国产精品合色在线| 成年人黄色毛片网站| 嫩草影视91久久| 女的被弄到高潮叫床怎么办 | 午夜日韩欧美国产| 免费搜索国产男女视频| 1000部很黄的大片| 久久草成人影院| 精品一区二区免费观看| 看片在线看免费视频| а√天堂www在线а√下载| 熟女电影av网| 久久久久性生活片| 日韩欧美一区二区三区在线观看| 免费一级毛片在线播放高清视频| а√天堂www在线а√下载| 国产熟女欧美一区二区| 51国产日韩欧美| eeuss影院久久| 一本精品99久久精品77| 狂野欧美激情性xxxx在线观看| 黄色一级大片看看| 久久久国产成人免费| 精品久久久久久,| 中文字幕人妻熟人妻熟丝袜美| 成年女人永久免费观看视频| 日韩精品中文字幕看吧| 在现免费观看毛片| 美女大奶头视频| 老女人水多毛片| 亚洲狠狠婷婷综合久久图片| av福利片在线观看| 在线免费观看不下载黄p国产 | 中文字幕久久专区| 校园人妻丝袜中文字幕| 日本-黄色视频高清免费观看| 在线看三级毛片| 久久久成人免费电影| 真人做人爱边吃奶动态| 99久久无色码亚洲精品果冻| 国产高清视频在线播放一区| 午夜福利成人在线免费观看| 哪里可以看免费的av片| 黄色配什么色好看| 精品一区二区三区视频在线| 五月伊人婷婷丁香| 大又大粗又爽又黄少妇毛片口| 午夜亚洲福利在线播放| 国产一级毛片七仙女欲春2| 99热这里只有是精品在线观看| www.色视频.com| 在线国产一区二区在线| videossex国产| 亚洲精品亚洲一区二区| 亚洲精品乱码久久久v下载方式| 亚洲精品粉嫩美女一区| 白带黄色成豆腐渣| 亚洲人与动物交配视频| 九九久久精品国产亚洲av麻豆| 看片在线看免费视频| 啪啪无遮挡十八禁网站| 男女啪啪激烈高潮av片| 欧美国产日韩亚洲一区| 狂野欧美激情性xxxx在线观看| 波野结衣二区三区在线| 午夜激情福利司机影院| 他把我摸到了高潮在线观看| netflix在线观看网站| 午夜福利在线观看吧| 国产精品久久久久久亚洲av鲁大| 国产久久久一区二区三区| 国产一区二区亚洲精品在线观看| 老熟妇乱子伦视频在线观看| 一进一出抽搐动态| 黄色女人牲交| 国产精品人妻久久久久久| 色精品久久人妻99蜜桃| 婷婷精品国产亚洲av在线| 天堂影院成人在线观看| 免费看美女性在线毛片视频| 永久网站在线| 人妻久久中文字幕网| 国产成人福利小说| 亚洲精品一卡2卡三卡4卡5卡| 美女 人体艺术 gogo| 他把我摸到了高潮在线观看| 搡老熟女国产l中国老女人| 国产精华一区二区三区| 久久国产乱子免费精品| 久久久精品大字幕| 岛国在线免费视频观看| 在线观看av片永久免费下载| 亚洲图色成人| 欧美性感艳星| 国产一级毛片七仙女欲春2| 婷婷丁香在线五月| 琪琪午夜伦伦电影理论片6080| av福利片在线观看| 麻豆av噜噜一区二区三区| 国产精品人妻久久久影院| 亚洲va日本ⅴa欧美va伊人久久| 在线免费十八禁| 国产伦精品一区二区三区视频9| 亚洲精华国产精华液的使用体验 | 日韩av在线大香蕉| 日韩欧美国产一区二区入口| 国产aⅴ精品一区二区三区波| 日本免费a在线| 91久久精品电影网| 在线播放国产精品三级| 超碰av人人做人人爽久久| av在线蜜桃| 亚洲美女搞黄在线观看 | 干丝袜人妻中文字幕| 亚洲熟妇中文字幕五十中出| 老熟妇乱子伦视频在线观看| 一级黄片播放器| 深夜精品福利| 成年版毛片免费区| 在线看三级毛片| 国产精品一区二区三区四区久久| 老熟妇乱子伦视频在线观看| 亚洲欧美激情综合另类| 日本精品一区二区三区蜜桃| 白带黄色成豆腐渣| xxxwww97欧美| 天堂√8在线中文| av天堂中文字幕网| 久久久国产成人精品二区| 精品久久久久久久久av| 国产高清不卡午夜福利| 亚洲一区高清亚洲精品| 精品久久国产蜜桃| 国产视频一区二区在线看| 中出人妻视频一区二区| 九九热线精品视视频播放| 亚洲电影在线观看av| 中文字幕免费在线视频6| 亚洲狠狠婷婷综合久久图片| 国内精品美女久久久久久| 久久精品国产亚洲av天美| 日韩欧美免费精品| 国产爱豆传媒在线观看| 国产 一区 欧美 日韩| 久久久久久久精品吃奶| av中文乱码字幕在线| 免费观看精品视频网站|