專利名稱:一種復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測方法
技術(shù)領(lǐng)域:
本發(fā)明涉及氣固流態(tài)化系統(tǒng)設(shè)計(jì)技術(shù)領(lǐng)域,具體地說,本發(fā)明涉及一種復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測方法。
背景技術(shù):
氣固流態(tài)化系統(tǒng)被廣泛應(yīng)用于石油、化工、制藥等行業(yè)中,其中一個(gè)典型的應(yīng)用實(shí)例就是煉油行業(yè)中的流化催化裂化(FCC)工藝。該工藝一般包含提升管、沉降器和再生器,但為了適應(yīng)日益嚴(yán)格的環(huán)境限制和逐漸豐富的產(chǎn)品需求,其具體組織形式逐漸呈現(xiàn)多樣化和復(fù)雜化。譬如,為了增加汽油辛烷值而降低其烯烴含量,中石化開發(fā)了一種以催化劑二次循環(huán)和提升管中段擴(kuò)徑為主要特征的多產(chǎn)異構(gòu)烷烴(MIP)的復(fù)雜FCC新工藝[許友好,張久順,龍軍.石油煉制與化工,2001,32:1-5]。該工藝中各反應(yīng)器或同一反應(yīng)器各段的構(gòu)型、尺寸以及操作條件(主要是指表觀氣速和顆粒通量)各不相同,而現(xiàn)有的商業(yè)計(jì)算流體 力學(xué)(CFD)軟件在處理這類具有復(fù)雜邊界和操作條件的系統(tǒng)時(shí)往往費(fèi)時(shí)費(fèi)力,有時(shí)甚至無能為力。實(shí)際上,雖然數(shù)值模擬能夠提供反應(yīng)器內(nèi)流體從非穩(wěn)態(tài)到穩(wěn)態(tài)的動態(tài)演化過程,但這類些工業(yè)裝置一般都希望在一個(gè)定態(tài)上操作,以實(shí)現(xiàn)裝置運(yùn)行和產(chǎn)品質(zhì)量的穩(wěn)定,因此穩(wěn)態(tài)動力學(xué)行為對此類工藝的設(shè)計(jì)和放大具有更重要的參考和指導(dǎo)意義。盡管已經(jīng)有大量實(shí)驗(yàn)數(shù)據(jù)和許多經(jīng)驗(yàn)關(guān)聯(lián)式能夠?qū)哂袌A柱狀等簡單構(gòu)型反應(yīng)器的穩(wěn)態(tài)動力學(xué)進(jìn)行預(yù)測,但目前仍然缺乏能夠合理描述復(fù)雜循環(huán)流化床反應(yīng)系統(tǒng)的整個(gè)回路、同一反應(yīng)器具有復(fù)雜變徑結(jié)構(gòu)的區(qū)段、或反應(yīng)系統(tǒng)操作在流域轉(zhuǎn)變點(diǎn)上時(shí)的定態(tài)流場的通用性理論或經(jīng)驗(yàn)預(yù)測方法。氣固流態(tài)化本身的復(fù)雜性主要表現(xiàn)在隨著表觀氣速和顆粒通量等操作條件的改變,顆粒層會呈現(xiàn)鼓泡、湍動、快速流化和稀相輸送等不同的被稱為“流域”的狀態(tài)。在多數(shù)流域,總體上會呈現(xiàn)上稀下濃的軸向S形濃度分布和中間稀邊壁濃的徑向環(huán)核結(jié)構(gòu),而在局部則呈現(xiàn)連續(xù)的稀相乳狀物包圍離散的密相顆粒團(tuán)聚物的特征。針對這種狀況,李靜海[顆粒流體兩相流多尺度作用模型和能量最小方法.博士學(xué)位論文.中國科學(xué)院過程工程研究所,北京.1987]提出了并流上行氣固流態(tài)化系統(tǒng)的能量最小多尺度(EMMS)模型。該模型提出并論證了氣固流化系統(tǒng)中的稀密兩相結(jié)構(gòu)是由顆粒懸浮輸送能最小決定的,且在考慮微觀流動機(jī)理(稀密相的單顆粒繞流和相界面上的壓降平衡)的基礎(chǔ)上,能夠準(zhǔn)確捕捉給定宏觀操作參數(shù)(表觀氣速和顆粒通量)、具有圓柱狀簡單構(gòu)型的氣固流化系統(tǒng)的顆粒濃度以及流域轉(zhuǎn)變等主要流體動力學(xué)特征。然而,目前的EMMS模型僅適用于圓柱狀的簡單構(gòu)型的氣固流化系統(tǒng),不能直接應(yīng)用于對具有復(fù)雜結(jié)構(gòu)和形狀的流化系統(tǒng)的整體動力學(xué)特征的預(yù)測。因此,當(dāng)前迫切需要一種能夠快速預(yù)測復(fù)雜氣固流化系統(tǒng)整體動力學(xué)特征的方法
發(fā)明內(nèi)容
本發(fā)明的目的是提供一種復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測方法。為實(shí)現(xiàn)上述發(fā)明目的,本發(fā)明提供了一種復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測方法,包括下列步驟I)將具有復(fù)雜幾何結(jié)構(gòu)和不同操作條件的氣固流化系統(tǒng)分解成若干基本流化單元,基本流化單元是形狀為圓柱狀、圓臺狀或倒圓臺狀,且操作條件單一的子流化單元;所述操作條件是指所述基本流化單元的表觀氣速和顆粒通量條件;2)將圓臺狀或倒圓臺狀的子流化單元進(jìn)一步近似為多個(gè)不同直徑圓柱狀片段的組合體;3)在已知?dú)夤塘骰到y(tǒng)氣體流率的基礎(chǔ)上設(shè)定氣固流化系統(tǒng)的顆粒流率; 4)根據(jù)當(dāng)前所設(shè)定的氣固流化系統(tǒng)的顆粒流率,計(jì)算各圓柱狀子流化單元和各圓臺狀或倒圓臺狀子流化單元的各圓柱狀片段的表觀氣速和顆粒通量條件;然后根據(jù)所計(jì)算的表觀氣速和顆粒通量條件,基于EMMS模型或相應(yīng)的經(jīng)驗(yàn)關(guān)聯(lián)式,計(jì)算各圓柱狀子流化單元或圓柱狀片段的顆粒濃度、壓降和顆粒存料量;5)根據(jù)上述各子單元顆粒濃度、壓降和顆粒存料量計(jì)算結(jié)果,在已知全系統(tǒng)氣體流率、總存料量以及顆粒流率控制閥門開度的情況下,判斷上述各圓柱狀子流化單元和圓柱狀片段的顆粒存料量和壓降之和是否同時(shí)滿足所述復(fù)雜流化系統(tǒng)的總存料量平衡、進(jìn)出口壓力平衡以及系統(tǒng)內(nèi)部各循環(huán)回路的壓力平衡條件;如果不滿足,則重新設(shè)定氣固流化系統(tǒng)的顆粒流率,回到步驟4);如果滿足,進(jìn)入步驟6);6)根據(jù)當(dāng)前所計(jì)算出的各圓柱狀子流化單元和各圓柱狀片段的顆粒濃度、壓降和顆粒存料量計(jì)算各子單元的顆粒濃度分布并按氣固流化系統(tǒng)的原有結(jié)構(gòu)進(jìn)行組合得出全系統(tǒng)的整體動力學(xué)特征。其中,所述步驟4)至5)中,利用基于Compute Unified DeviceArchitecture (CUDA)、OpenMP 或 Message Passing Interface (MPI)的并行計(jì)算技術(shù)加速EMMS模型的求解以及加速對是否滿足氣固流化系統(tǒng)各平衡條件的判斷。其中,所述步驟4)中,對于處于密相湍動流化狀態(tài)的圓柱狀子流化單元,基于相應(yīng)的經(jīng)驗(yàn)關(guān)聯(lián)式進(jìn)行計(jì)算;對于其余圓柱狀子流化單元,基于EMMS模型進(jìn)行計(jì)算。其中,所述步驟I)中,當(dāng)所述氣固流化系統(tǒng)具有非流化部分時(shí),還應(yīng)分解出相應(yīng)的非流化子單元;所述步驟4)中,基于相應(yīng)的經(jīng)驗(yàn)關(guān)聯(lián)式進(jìn)行計(jì)算所述非流化子單元的顆粒濃度、壓降和顆粒存料量。其中,所述步驟6)中,所述全系統(tǒng)的整體動力學(xué)特征是所述氣固流化系統(tǒng)的顆粒濃度分布和/或顆粒流率。與現(xiàn)有技術(shù)相比,本發(fā)明具有下列技術(shù)效果I、利用本發(fā)明所提供的方法能夠?qū)Χ喾N復(fù)雜流化系統(tǒng)進(jìn)行實(shí)時(shí)或準(zhǔn)實(shí)時(shí)全循環(huán)穩(wěn)態(tài)模擬,以為諸如FCC工藝等復(fù)雜氣固流化系統(tǒng)的動力學(xué)設(shè)計(jì)、放大以及在線優(yōu)化提供
定量參考。2、本發(fā)明尤其適用于對煉油行業(yè)中的流化催化裂化(FCC)系統(tǒng)催化劑循環(huán)量及全局催化劑濃度分布的預(yù)測。
以下,結(jié)合附圖來詳細(xì)說明本發(fā)明的實(shí)施例,其中圖I示出了本發(fā)明的原理示意圖;圖2示出了本發(fā)明的流程示意圖;圖3示出了本發(fā)明一個(gè)實(shí)施例中的多產(chǎn)異構(gòu)烷烴FCC工藝簡圖;圖4示出了本發(fā)明一個(gè)實(shí)施例中所計(jì)算出的多產(chǎn)異構(gòu)烷烴FCC裝置內(nèi)的催化劑濃度分布。
具體實(shí)施例方式EMMS模型能夠準(zhǔn)確捕捉給定宏觀操作參數(shù)(表觀氣速和顆粒通量)、具有圓柱狀簡單構(gòu)型的氣固流化系統(tǒng)的顆粒濃度以及流域轉(zhuǎn)變等主要流體動力學(xué)特征。然而,目前的 EMMS模型僅適用于圓柱狀簡單構(gòu)型的氣固流化系統(tǒng),不能直接應(yīng)用于對具有復(fù)雜結(jié)構(gòu)和形狀的流化系統(tǒng)的整體動力學(xué)特征的預(yù)測。本案發(fā)明人分析了 EMMS模型的優(yōu)勢和不足,在對復(fù)雜流化系統(tǒng)進(jìn)行幾何分解的基礎(chǔ)上,以EMMS模型的全局表達(dá)為主體輔以部分合理的經(jīng)驗(yàn)關(guān)聯(lián)式,最終實(shí)現(xiàn)了對復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測。參考圖1,本發(fā)明中,先將復(fù)雜流化系統(tǒng)分解為構(gòu)型簡單和條件單一的子系統(tǒng),然后根據(jù)EMMS模型和經(jīng)驗(yàn)關(guān)聯(lián)式計(jì)算各子系統(tǒng)滿足全系統(tǒng)物料和壓降平衡的存料量和壓降動力學(xué)特性,然后再將各子系統(tǒng)的動力學(xué)特性組合,得到復(fù)雜流化系統(tǒng)的整體動力學(xué)特性。參考圖2,本發(fā)明的復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測方法原理如下I、將具有復(fù)雜幾何結(jié)構(gòu)和不同操作條件的氣固流化系統(tǒng)分解成若干基本流化單元,基本流化單元是形狀為圓柱狀、圓臺狀或倒圓臺狀,且操作條件單一的子流化單元;所述操作條件是指所述基本流化單元的表觀氣速和顆粒通量條件;本步驟中,當(dāng)所述氣固流化系統(tǒng)具有非流化部分時(shí),還應(yīng)分解出相應(yīng)的非流化子單元。2、將圓臺狀或倒圓臺狀的子流化單元進(jìn)一步近似為多個(gè)不同直徑圓柱狀片段的組合體??梢愿鶕?jù)實(shí)際需求來確定本步驟的不同直徑圓柱狀片段的組合體與原始的子流化單元的近似精度,這是本領(lǐng)域技術(shù)人員易于理解的。比如,該組合體的體積與實(shí)際子流化單元的體積相對誤差不大于1%時(shí),即可以保證由此導(dǎo)致的顆粒存料量的相對計(jì)算誤差小于1%。3、在已知?dú)夤塘骰到y(tǒng)氣體流率的基礎(chǔ)上設(shè)定氣固流化系統(tǒng)的顆粒流率。4、根據(jù)當(dāng)前所設(shè)定的氣固流化系統(tǒng)的顆粒流率,計(jì)算各圓柱狀子流化單元和各圓臺狀或倒圓臺狀子流化單元的各圓柱狀片段的表觀氣速和顆粒通量條件。根據(jù)所計(jì)算的表觀氣速和顆粒通量條件,基于EMMS模型或相應(yīng)的經(jīng)驗(yàn)關(guān)聯(lián)式,計(jì)算各圓柱狀子流化單元或圓柱狀片段的顆粒濃度、壓降和顆粒存料量。本步驟中,對于處于密相湍動流化狀態(tài)的圓柱狀子流化單元,基于相應(yīng)的經(jīng)驗(yàn)關(guān)聯(lián)式進(jìn)行計(jì)算;對于其余圓柱狀子流化單元,基于EMMS模型進(jìn)行計(jì)算。基于相應(yīng)的經(jīng)驗(yàn)關(guān)聯(lián)式進(jìn)行計(jì)算所述非流化子單元的顆粒濃度、壓降和顆粒存料量。5、根據(jù)上述各子單元顆粒濃度、壓降和顆粒存料量計(jì)算結(jié)果,在已知全系統(tǒng)氣體流率、總存料量以及顆粒流率控制閥門開度的情況下,判斷上述各圓柱狀子流化單元和圓柱狀片段的顆粒存料量和壓降之和是否同時(shí)滿足所述復(fù)雜流化系統(tǒng)的總存料量平衡、進(jìn)出口壓力平衡以及系統(tǒng)內(nèi)部各循環(huán)回路的壓力平衡條件;如果不滿足,則重新設(shè)定氣固流化系統(tǒng)的顆粒流率,回到步驟4 ;如果滿足,進(jìn)入步驟6。6、根據(jù)當(dāng)前所計(jì)算出的各圓柱狀子流化單元和各圓柱狀片段的顆粒濃度、壓降和顆粒存料量得出全系統(tǒng)的整體動力學(xué)特征,例如所述氣固流化系統(tǒng)的顆粒濃度分布和/或顆粒流率。在一個(gè)實(shí)施例中,其中步驟4和5可以利用基于CUDA(或者OpenMP,或者M(jìn)PI等)并行計(jì)算技術(shù)加速EMMS模型的求解以及加速對是否滿足氣固流化系統(tǒng)各平衡條件的判斷。這樣可以進(jìn)一步提高預(yù)測結(jié)果的實(shí)時(shí)性。并行計(jì)算時(shí),硬件設(shè)備可采用(但不限于)GPU+CPU或者GPU+CPU異構(gòu)集群。下面以圖3的多產(chǎn)異構(gòu)烷烴FCC工藝的復(fù)雜流化系統(tǒng)為例,對本發(fā)明的快速預(yù)測方法進(jìn)行詳細(xì)說明。參考圖3,多產(chǎn)異構(gòu)烷烴FCC工藝主要由提升管、沉降器以及再生器所組成。其中,提升管反應(yīng)器由下至上又可以分為底部的裂解區(qū)、中間的擴(kuò)大段轉(zhuǎn)化區(qū)和上部的汽提區(qū)。 再生器也可以分為下部的燃燒區(qū)和上部的再生區(qū)。從沉降器出來的待生催化劑一部分返回到提升管轉(zhuǎn)化區(qū)的下部,其余部分被輸送到再生器進(jìn)行再生后返回到提升管裂解區(qū)的下部。本優(yōu)選實(shí)施例提供的復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測方法包括下列步驟I、根據(jù)幾何構(gòu)型和流化狀態(tài)對多產(chǎn)異構(gòu)烷烴FCC工藝進(jìn)行分解根據(jù)多產(chǎn)異構(gòu)烷烴FCC工藝的幾何結(jié)構(gòu)以及各部分流化狀態(tài)的不同,將其分解為流化子系統(tǒng)和非流化子系統(tǒng)兩大類。流化子系統(tǒng)包括裂解區(qū)I、裂解轉(zhuǎn)化過渡區(qū)2、轉(zhuǎn)化區(qū)3、轉(zhuǎn)化汽提過渡區(qū)4、汽提區(qū)5、燃燒區(qū)6、燃燒再生過渡區(qū)底部到再生催化劑循環(huán)管入口區(qū)域7、再生催化劑循環(huán)管入口到燃燒再生過渡區(qū)頂部區(qū)域8、再生區(qū)底部密相區(qū)9和再生區(qū)頂部稀相區(qū)10,共10個(gè)區(qū)域。非流化子系統(tǒng)包括沉降器11、待生催化劑循環(huán)管12、再生催化劑循環(huán)管13和MIP待生催化劑循環(huán)管14共4個(gè)區(qū)域。圖3中,15為再生催化劑循環(huán)閥門;16為MIP待生催化劑循環(huán)閥門。上述10個(gè)流化子系統(tǒng)均為基本流化單元,即形狀為圓柱狀、圓臺狀或倒圓臺狀,且操作條件單一的子流化單元;所述操作條件是指所述基本流化單元的表觀氣速和顆粒通量(本實(shí)施例中為催化劑通量)條件。分解出流化子系統(tǒng)后,記錄每個(gè)流化子系統(tǒng)的幾何參數(shù)以備后續(xù)計(jì)算。對于圓柱狀流化子系統(tǒng),記錄其高度和橫截面直徑,對于圓臺狀和倒圓臺狀流化子系統(tǒng),記錄其高度以及上表面和下表面的直徑。2、將圓臺狀或倒圓臺狀的流化子系統(tǒng)進(jìn)一步分解將上述10個(gè)流化子系統(tǒng)中圓臺狀或倒圓臺狀的流化子系統(tǒng)(即裂解轉(zhuǎn)化過渡區(qū)2、轉(zhuǎn)化汽提過渡區(qū)4、燃燒再生過渡區(qū)底部到再生催化劑循環(huán)管入口區(qū)域7、再生催化劑循環(huán)管入口到燃燒再生過渡區(qū)頂部區(qū)域8)進(jìn)一步近似為多個(gè)不同直徑圓柱狀片段的組合體,且該組合體的體積與實(shí)際流化子系統(tǒng)的體積相對誤差不大于I %,這樣可以保證由此導(dǎo)致的顆粒存料量的相對計(jì)算誤差小于I %,從而保證后續(xù)顆粒濃度和壓降計(jì)算的準(zhǔn)確性。當(dāng)然,體積相對誤差的上限取值并不限于I %,本領(lǐng)域技術(shù)人員可以根據(jù)實(shí)際需求來決定該體積相對誤差的上限取值。
在將圓臺狀或倒圓臺狀的流化子系統(tǒng)近似為多個(gè)不同直徑圓柱狀片段的組合體時(shí),記錄每個(gè)圓柱狀片段的高度和橫截面直徑,以備后續(xù)計(jì)算。3、設(shè)定多產(chǎn)異構(gòu)烷烴FCC工藝的操作條件和初始催化劑流率一般地,假定催化劑顆粒的直徑和密度在整個(gè)系統(tǒng)中保持不變,提升管流體從裂解區(qū)底部一次性引入。提升管流體的物性取轉(zhuǎn)化區(qū)平均溫度下的油氣和水蒸氣平均物性,再生器流體的物性取再生區(qū)平均溫度下的煙氣物性,如表I所示。表I
權(quán)利要求
1.一種復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測方法,包括下列步驟 1)將具有復(fù)雜幾何結(jié)構(gòu)和不同操作條件的氣固流化系統(tǒng)分解成若干基本流化單元,基本流化單元是形狀為圓柱狀、圓臺狀或倒圓臺狀,且操作條件單一的子流化單元;所述操作條件是指所述基本流化單元的表觀氣速和顆粒通量條件; 2)將圓臺狀或倒圓臺狀的子流化單元進(jìn)一步近似為多個(gè)不同直徑圓柱狀片段的組合體; 3)在已知?dú)夤塘骰到y(tǒng)氣體流率的基礎(chǔ)上設(shè)定氣固流化系統(tǒng)的顆粒流率; 4)根據(jù)當(dāng)前所設(shè)定的氣固流化系統(tǒng)的顆粒流率,計(jì)算各圓柱狀子流化單元和各圓臺狀或倒圓臺狀子流化單元的各圓柱狀片段的表觀氣速和顆粒通量條件;然后根據(jù)所計(jì)算的表 觀氣速和顆粒通量條件,基于EMMS模型或相應(yīng)的經(jīng)驗(yàn)關(guān)聯(lián)式,計(jì)算各圓柱狀子流化單元或圓柱狀片段的顆粒濃度、壓降和顆粒存料量; 5)根據(jù)上述各子單元顆粒濃度、壓降和顆粒存料量計(jì)算結(jié)果,在已知全系統(tǒng)氣體流率、總存料量以及顆粒流率控制閥門開度的情況下,判斷上述各圓柱狀子流化單元和圓柱狀片段的顆粒存料量和壓降之和是否同時(shí)滿足所述復(fù)雜流化系統(tǒng)的總存料量平衡、進(jìn)出口壓力平衡以及系統(tǒng)內(nèi)部各循環(huán)回路的壓力平衡條件;如果不滿足,則重新設(shè)定氣固流化系統(tǒng)的顆粒流率,回到步驟4);如果滿足,進(jìn)入步驟6); 6)根據(jù)當(dāng)前所計(jì)算出的各圓柱狀子流化單元和各圓柱狀片段的顆粒濃度、壓降和顆粒存料量計(jì)算各子單元的顆粒濃度分布并按所述氣固流化系統(tǒng)的原有結(jié)構(gòu)進(jìn)行組合得出全系統(tǒng)的整體動力學(xué)特征。
2.根據(jù)權(quán)利要求I所述的復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測方法,其特征在于,所述步驟4)中,對于處于密相湍動流化狀態(tài)的圓柱狀子流化單元,基于相應(yīng)的經(jīng)驗(yàn)關(guān)聯(lián)式進(jìn)行計(jì)算;對于其余圓柱狀子流化單元,基于EMMS模型進(jìn)行計(jì)算。
3.根據(jù)權(quán)利要求I所述的復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測方法,其特征在于,當(dāng)所述氣固流化系統(tǒng)具有非流化部分時(shí),所述步驟I)中,所述氣固流化系統(tǒng)分解成若干所述基本流化單元和非流化子單元;所述步驟4)還包括,基于相應(yīng)的經(jīng)驗(yàn)關(guān)聯(lián)式進(jìn)行計(jì)算所述非流化子單元的顆粒濃度、壓降和顆粒存料量。
4.根據(jù)權(quán)利要求I所述的復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測方法,其特征在于,所述步驟4)至5)中,利用基于CUDA、OpenMP或MPI的并行計(jì)算技術(shù)對EMMS模型進(jìn)行求解以及對是否滿足氣固流化系統(tǒng)各平衡條件進(jìn)行判斷。
5.根據(jù)權(quán)利要求I所述的復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測方法,其特征在 于,所述步驟6)中,所述全系統(tǒng)的整體動力學(xué)特征是所述氣固流化系統(tǒng)的顆粒濃度分布和/或顆粒流率。
全文摘要
本發(fā)明涉及一種復(fù)雜流化系統(tǒng)整體動力學(xué)特性的快速預(yù)測方法,包括1)將流化系統(tǒng)分解成若干圓柱狀、圓臺狀或倒圓臺狀且操作條件單一的子流化單元;2)將圓臺狀或倒圓臺狀的子流化單元近似為多個(gè)圓柱狀片段;3)設(shè)定氣固流化系統(tǒng)的顆粒流率;4)基于EMMS模型或相應(yīng)的經(jīng)驗(yàn)關(guān)聯(lián)式,計(jì)算各圓柱狀子流化單元或圓柱狀片段的顆粒濃度、壓降和顆粒存料量;5)判斷顆粒存料量和壓降之和是否同時(shí)滿足系統(tǒng)的平衡條件;如不滿足,則重新設(shè)定顆粒流率,回到步驟4);如滿足,進(jìn)入步驟6);6)計(jì)算各子流化單元的顆粒濃度分布并得出全系統(tǒng)的整體動力學(xué)特征。本發(fā)明能夠?qū)Χ喾N復(fù)雜流化系統(tǒng)進(jìn)行實(shí)時(shí)或準(zhǔn)實(shí)時(shí)全循環(huán)穩(wěn)態(tài)模擬。
文檔編號G01M10/00GK102778344SQ201110122298
公開日2012年11月14日 申請日期2011年5月12日 優(yōu)先權(quán)日2011年5月12日
發(fā)明者劉新華, 李靜海, 葛蔚 申請人:中國科學(xué)院過程工程研究所