基于金屬氫化物高溫蓄熱的氫熱耦合傳遞機理

劉洋，Iqra AYUB，楊福勝，吳震，張早校

（西安交通大學(xué)化學(xué)工程與技術(shù)學(xué)院，動力工程多相流國家重點實驗室，陜西西安 710049）

引言

由于化石能源儲量有限并且污染嚴(yán)重，太陽能等可再生能源的規(guī)模化應(yīng)用成為發(fā)展趨勢。國際能源署（IEA）預(yù)測，到2050年，全球大約20%的電力將來源于太陽能[1]。其中，太陽能熱發(fā)電（CSP）具有適應(yīng)性強、適于大規(guī)模集中式發(fā)電等特點，越來越受歡迎，并且已經(jīng)有商業(yè)化的電廠投入運行[2]。

由于太陽能具有不穩(wěn)定和不連續(xù)的特點，為保障太陽能熱發(fā)電過程能夠具有持續(xù)性并具備一定的調(diào)峰能力，需要配置一個蓄熱系統(tǒng)，將白天多余的太陽能以熱能的形式儲存起來，在太陽能供應(yīng)不足時（陰天或晚上），由蓄熱系統(tǒng)釋放熱能進(jìn)行發(fā)電。此外，蓄熱技術(shù)開發(fā)還可以推動冷熱聯(lián)供、智能電網(wǎng)、余熱回收等領(lǐng)域的發(fā)展[3]。

根據(jù)蓄熱原理的不同，高溫蓄熱技術(shù)主要分為顯熱蓄熱[4-5]、潛熱蓄熱[6-7]以及熱化學(xué)能蓄熱。熱化學(xué)能蓄熱如金屬氫化物、氫氧化物、碳酸鹽等[8-9]，其相對前兩者蓄熱密度大[10]，但技術(shù)復(fù)雜，目前處于實驗室研究階段[11]。近年來，金屬氫化物蓄熱逐漸被重視起來，其屬于熱化學(xué)蓄熱的一種，具有高的蓄熱密度且安全無腐蝕，是一種極具潛力的蓄熱技術(shù)[12]。

金屬氫化物的吸放熱原理：

式中,M 為金屬或合金；H2為氫氣；MHx為金屬氫化物；ΔH為反應(yīng)焓，kJ/mol。從式（1）中可以看到，金屬氫化物的吸放熱過程伴隨著氫氣的釋放或消耗，因此，金屬氫化物蓄熱是熱質(zhì)傳遞、氫氣流動以及化學(xué)反應(yīng)的復(fù)雜耦合過程，探究其復(fù)雜氫熱耦合傳遞機理以及明確其性能控制因素對金屬氫化物高溫蓄熱應(yīng)用是非常必要的。

目前對金屬氫化物高溫蓄熱技術(shù)進(jìn)行研究主要集中在材料改性、反應(yīng)器優(yōu)化和系統(tǒng)集成優(yōu)化三方面。Mg 基金屬氫化物是最受關(guān)注的一種蓄熱材料[13]，但是純MgH2具有熱導(dǎo)率低、反應(yīng)動力學(xué)差的問題。相比之下，Mg2FeH6具有更好的運行溫度和平衡壓力，但是Urbanczyk 等[14]發(fā)現(xiàn)在循環(huán)期間，會有少量的Fe和MgH2殘留，這說明Mg2FeH6合成較為困難。Poupin 等[15]將TiMn1.5Hx用作儲氫材料，Mg2FeH6用作高溫蓄熱材料，研究發(fā)現(xiàn)當(dāng)Mg2FeH6在350～500℃之間運行時，可以抵消TiMn1.5在氫氣吸收和釋放之間顯示的壓力滯后。Humphries 等[16]對Na2Mg2NiH6進(jìn)行了研究，發(fā)現(xiàn)在10 個循環(huán)內(nèi)，沒有發(fā)生氫損失，并且在395℃下可以實現(xiàn)快速反應(yīng)。對于材料改性，多集中于摻雜其余元素來提高M(jìn)g基金屬氫化物的性能。

在反應(yīng)器優(yōu)化方面，大部分研究集中在通過強化床層和換熱管的傳熱效果來改善反應(yīng)器性能[17]。Keshari 等[18]將銅針翅片和換熱管結(jié)合，通過對粉末床反應(yīng)器內(nèi)瞬態(tài)行為的研究，驗證了該結(jié)構(gòu)的優(yōu)越性。Ayub 等[19]設(shè)計了錐形翅片式粉末床反應(yīng)器，發(fā)現(xiàn)翅片結(jié)構(gòu)能顯著縮短換熱時間。Tong等[20]設(shè)計了兩種新型粉末床反應(yīng)器，包括螺旋換熱管和直換熱管的組合以及同心雙螺旋換熱管。結(jié)果表明，同心雙螺旋換熱管結(jié)構(gòu)能有效地提高反應(yīng)器的效率。Bao等[21]提出了一種以石墨含量為梯度的改良壓塊層結(jié)構(gòu)。模擬結(jié)果表明，改良壓塊層結(jié)構(gòu)可以提高反應(yīng)器的性能。相比于粉末床，壓塊床具有良好的導(dǎo)熱性和穩(wěn)定性，但其床層結(jié)構(gòu)限制了換熱管的形式。

對于金屬氫化物蓄熱系統(tǒng)的研究多集中于研究其動態(tài)特性以及驗證其可行性。Nyamsi 等[22]以LaNi5為儲氫材料，Mg、Mg2Ni 和Mg2FeH6分別作為蓄熱材料，對蓄熱系統(tǒng)性能進(jìn)行了研究。結(jié)果表明，在適宜的操作條件下，Mg 基金屬氫化物與LaNi5具有良好的耦合效應(yīng)?；?0 MJ 的蓄熱容量，Jenne等[23]建立了MgH2-MnNi4.6Al0.4H6的耦合系統(tǒng)模型，并對其動態(tài)特性進(jìn)行了研究。結(jié)果表明，蓄熱與儲氫反應(yīng)器之間的氫氣流動會被限制，采用壓縮機可以有效解決這一問題。d'Entremont 等[24]研究了以NaMgH2F-TiCr1.6Mn0.2為基礎(chǔ)的蓄熱系統(tǒng)的可行性，結(jié)果表明，所提出的金屬氫化物系統(tǒng)可以與高溫蒸汽發(fā)電裝置相結(jié)合。

目前對于金屬氫化物氫熱耦合機理的研究報道大多是圍繞中低溫材料展開的。Hahne 等[25]通過測量中低溫金屬氫化物材料的有效熱導(dǎo)率，發(fā)現(xiàn)氫氣壓力對床層有效熱導(dǎo)率的影響規(guī)律。Sun 等[26]基于均一分布模型建立了金屬氫化物床層有效熱導(dǎo)率的理論模型，并與實驗結(jié)果吻合較好。Matsushita等[27]將吸放氫過程中的膨脹收縮現(xiàn)象的影響引入接觸面積模型，修正后的模型與實驗結(jié)果吻合較好。Ishido 等[28]測量了MgH2粉末床層有效熱導(dǎo)率隨著孔隙率、氫氣壓力和溫度的變化情況，并指出可以忽略顆粒之間的輻射換熱。趙長穎等[29]研究表明熱量的傳遞是影響系統(tǒng)蓄放熱效率的關(guān)鍵因素之一。Shen 等[30]采用了多孔金屬來增加床層熱導(dǎo)率，并對其放熱過程進(jìn)行了熱分析。綜上所述，對于金屬氫化物蓄熱技術(shù)的研究大多集中于材料改性、反應(yīng)器優(yōu)化和系統(tǒng)分析三個方面，但是，當(dāng)前對于金屬氫化物在高溫下的氫熱耦合機理研究仍不夠充分。

為了明確金屬氫化物高溫蓄熱反應(yīng)器中復(fù)雜的氫熱耦合傳遞機理，本文建立了鎂-氫化鎂高溫蓄熱反應(yīng)器的多物理場耦合模型，以期為金屬氫化物高溫蓄熱性能進(jìn)一步優(yōu)化研究提供明確方向。

1 蓄熱反應(yīng)器數(shù)學(xué)模型

1.1 幾何模型

本文以管式金屬氫化物高溫蓄熱反應(yīng)器為研究對象，如圖1（a）所示。為了簡化計算，基于對稱性，選取反應(yīng)器床層的一部分為研究對象，建立二維軸對稱模型對反應(yīng)器進(jìn)行研究，如圖1（b）所示。金屬氫化物床層內(nèi)部填充了球磨的MgH2粉末，用作高溫蓄熱材料。換熱管位于床層中央，內(nèi)外徑分別為4 mm 和6.4 mm，由316L 不銹鋼制成。為了研究反應(yīng)器高溫時的放熱性能，這里采用過熱蒸汽作為換熱流體。

圖1 金屬氫化物高溫蓄熱反應(yīng)器幾何模型Fig.1 Geometric model of metal hydride high temperature thermal storage reactor

1.2 數(shù)值模型

在金屬氫化物反應(yīng)器建模過程中，采用的主要假設(shè)如下：

（1）在反應(yīng)過程中床層內(nèi)部處于局部熱平衡狀態(tài)[31]；

（2）在反應(yīng)過程中床層材料的物性保持恒定[31]；

（3）忽略黏性耗散；

（4）忽略反應(yīng)器對環(huán)境的熱耗散[31]；

（5）忽略輻射傳熱。

將金屬氫化物床層看成多孔介質(zhì)，其能量方程為：

式中，ρe為床層有效密度，kg/m3；cp,e為床層有效比熱容，J/(kg·K)；Tbed為床層溫度，℃或K；λe為床層有效熱導(dǎo)率，W/(m·K)；k為單位體積摩爾反應(yīng)速率，mol/(m3·s)；t為時間，s。

其中：

式中，ρg、ρs分別為氫氣和固體的密度，kg/m3；cp,g、cp,s分別為氫氣和固體的比熱容，J/(kg·K)；ε為孔隙率；csat為單位體積金屬氫化物中H2的飽和濃度，mol/m3；X為反應(yīng)分?jǐn)?shù)；wt為金屬氫化物中H 的飽和質(zhì)量含量；MH2為氫氣的摩爾質(zhì)量，kg/mol。

Mg的吸氫動力學(xué)方程[32]如下：

式中，Ca為吸氫速率常數(shù)，s-1；pea為吸氫平衡壓力，MPa；Ea為反應(yīng)活化能，kJ/mol；Rg為通用氣體常數(shù)，J/(mol·K)。

平衡壓力由Van't Hoff 方程確定：

式中，ΔS為反應(yīng)熵，J/(mol·K)；參考壓力pref=0.1 MPa。

氫氣被認(rèn)為是可壓縮氣體，氫氣在多孔床層內(nèi)的流動采用Darcy定律描述，換熱流體管壁采用固體導(dǎo)熱方程，采用k-?模型描述換熱流體的流動情況。

1.3 模型驗證

關(guān)于MgH2高溫蓄熱的實驗研究，考慮到Chaise等[32]的研究所采用的反應(yīng)器和研究工況與本文最為接近，本文的模型驗證采用其實驗數(shù)據(jù)作為參照，實驗工況為反應(yīng)器初始溫度300℃，氫氣入口壓力0.77 MPa。數(shù)值模型采用COMSOL Multiphysics v5.1進(jìn)行求解。對比結(jié)果如圖2所示。從模擬結(jié)果與實驗結(jié)果的對比來看，床層中心溫度和平均反應(yīng)分?jǐn)?shù)的變化曲線均有較好吻合，說明本文所建立的MgH2反應(yīng)器多物理場耦合模型是合理且有效的。其中，模擬和實驗的平均反應(yīng)分?jǐn)?shù)吻合性較好，但是模擬和實驗的床層中心溫度之間存在一定的誤差，這是因為在實驗過程中反應(yīng)器會對外界環(huán)境散熱，越靠近反應(yīng)器壁面，散熱量越高，然而在模擬工作中忽略了環(huán)境熱耗散。這就導(dǎo)致隨著反應(yīng)的進(jìn)行，模擬溫度和實驗溫度之間存在一定的誤差。同時所選模型經(jīng)過了網(wǎng)格無關(guān)性驗證。

圖2 金屬氫化物高溫蓄熱反應(yīng)模型驗證Fig.2 Model verification

2 氫熱耦合機理分析與討論

以充氫放熱過程為例，研究金屬氫化物高溫蓄熱反應(yīng)器的氫熱耦合傳遞機理。圖3為金屬氫化物高溫蓄熱反應(yīng)器的瞬態(tài)輸出溫度和床層平均反應(yīng)分?jǐn)?shù)隨時間的變化。在本文的研究中，將整個熱輸出分為三個階段：初始階段、平臺階段和衰退階段。此外，以峰值溫度為界，將初始階段分為初始上升階段和初始下降階段。

圖3 床層平均反應(yīng)分?jǐn)?shù)和換熱流體輸出溫度隨時間的變化Fig.3 Variation curves of average bed reaction fraction and heat transfer fluid output temperature with time

2.1 初始上升階段分析

根據(jù)圖3 可知，初始上升階段為0～36 s，選取5個時間點進(jìn)行研究，取床層壓力和床層溫度為研究對象，如圖4所示。

從圖4(a)中能夠看到，氫氣并不是一開始就充滿整個床層，而是逐步充滿的。氫氣是反應(yīng)的必備條件，因此反應(yīng)床層溫度也會逐步自上而下升高[圖4(b)]。對比任何同一時間段下的溫度和壓力圖可以發(fā)現(xiàn)，高溫區(qū)的幾何長度要比高壓區(qū)長。這是因為，反應(yīng)放出的熱量不僅會傳遞給換熱流體，還會有一部分在床層內(nèi)部傳遞。總之，初始上升階段是氫氣充滿整個床層的過程，即反應(yīng)器床層內(nèi)部材料自上而下逐步反應(yīng)的過程；對于換熱流體而言，沿著流體流向，有效換熱長度迅速升高，導(dǎo)致出口溫度迅速升高。由于氫壓的變化會造成床層充滿氫氣時間的變化，進(jìn)而影響床層反應(yīng)，因此初始上升階段的控制因素是氫壓。

圖4 初始階段不同時刻下的床層參數(shù)分布（從左至右時刻依次為7、14、21、28、36 s）Fig.4 Distribution of bed parameters at different time of initial stage

2.2 初始下降階段分析

根據(jù)圖3 可知，初始下降階段為36～200 s，選取5 個時間點進(jìn)行研究，取床層反應(yīng)分?jǐn)?shù)為研究對象，如圖5所示。

圖5 初始下降階段不同時刻下的床層反應(yīng)分?jǐn)?shù)分布（從左至右時刻依次為36、80、120、160、200 s）Fig.5 Distribution of bed reaction fraction at different time in initial descending stage

床層反應(yīng)分?jǐn)?shù)代表了床層的反應(yīng)完成程度，反應(yīng)分?jǐn)?shù)越高表示反應(yīng)完成程度越高。由圖5 可見，整個初始下降階段，反應(yīng)都發(fā)生在床層和換熱流體的接觸面附近。這是因為經(jīng)過初始上升階段，氫氣雖然充滿整個床層，床層內(nèi)部均發(fā)生反應(yīng)，但是，由于金屬氫化物的反應(yīng)是可逆的，過高的床層溫度會限制反應(yīng)的進(jìn)一步進(jìn)行，只有將充氫放熱過程產(chǎn)生的熱量及時導(dǎo)出，才會促進(jìn)充氫放熱反應(yīng)的進(jìn)行。因此，對于整個床層來說，靠近換熱流體管束的區(qū)域會率先完全反應(yīng)。同時，觀察圖5可以發(fā)現(xiàn)：隨著時間的推移，換熱管壁面附近的床層反應(yīng)是沿著流體流向逐步完成的。這是因為沿著流體流向，換熱流體的溫度逐漸升高，傳熱溫差逐漸變小，導(dǎo)致床層反應(yīng)速率逐步變小。此外，同時觀察到：在200 s時，換熱管壁面附近均完全反應(yīng)，并且在已經(jīng)反應(yīng)完全和未反應(yīng)完全區(qū)域之間形成一條明顯的分界線，即反應(yīng)鋒面，代表了床層內(nèi)部的反應(yīng)程度。因此，初始下降階段是反應(yīng)鋒面逐漸脫離換熱管壁的過程，換熱流體和床層的直接接觸面積減小，導(dǎo)致了出口溫度迅速下降。因此初始下降階段的控制因素是接觸熱阻。

2.3 平臺階段分析

根據(jù)圖3 可知，平臺階段為200～9000 s，選取5 個時間點進(jìn)行研究，同樣取床層反應(yīng)分?jǐn)?shù)為研究對象，如圖6所示。

圖6 平臺階段不同時刻下的床層反應(yīng)分?jǐn)?shù)分布(從左至右時刻依次為500、2000、4000、6000、9000 s)Fig.6 Distribution of bed reaction fraction at different time in platform stage

從圖6可以看出：在整個平臺階段，反應(yīng)鋒面逐漸向床層外側(cè)移動，直到換熱流體進(jìn)口側(cè)的床層完全反應(yīng)。結(jié)合圖3中的換熱流體出口溫度的變化趨勢，可以看出在平臺階段，換熱流體的輸出溫度較為平穩(wěn)，是反應(yīng)器的有效熱輸出階段。同時，在9000 s時，床層的進(jìn)口側(cè)完全反應(yīng)，但是沿著流體流向，由于換熱流體的傳熱溫差逐漸減小，未反應(yīng)區(qū)域逐漸增加，出現(xiàn)非均勻反應(yīng)。可以看到，非均勻反應(yīng)導(dǎo)致?lián)Q熱流體輸出溫度的提前下降，是導(dǎo)致?lián)Q熱流體有效熱輸出時間縮短的主要原因。因此，平臺階段反應(yīng)鋒面在床層內(nèi)部移動，床層熱阻成為影響平臺階段溫度輸出的主要因素。

2.4 衰退階段分析

根據(jù)圖3 可知，衰退階段為9000～14000 s，但是在10500 s 時床層內(nèi)部均已反應(yīng)完全，因此以10500 s為時間節(jié)點，將衰退階段分為兩個階段。在9000～10500 s 內(nèi)選取5 個時間點進(jìn)行研究，取床層反應(yīng)分?jǐn)?shù)為研究對象（圖7）；在10500～14000 s 內(nèi)選取5 個時間點進(jìn)行研究，取床層溫度為研究對象（圖8）。

從平臺階段的分析可知，在9000 s 時床層內(nèi)部出現(xiàn)非均勻反應(yīng)。由圖7 可見，在9000～10500 s 內(nèi)，床層的非均勻反應(yīng)部分逐漸完全反應(yīng)，到10500 s時床層反應(yīng)完全，因此在衰退階段的前期，床層的非均勻反應(yīng)部分是換熱流體溫升的主要原因。從圖8可以看到10500 s時床層完全反應(yīng)，但是床層內(nèi)部還有余熱存在。這是因為金屬氫化物床層的導(dǎo)熱性較差，產(chǎn)生的熱量并不能夠及時傳遞到換熱流體，尤其是床層外側(cè)，會造成床層內(nèi)部存在反應(yīng)余熱。由圖8 可見，10500～14000 s 是床層余熱逐漸被換熱流體吸收的過程。綜上，在衰退階段的前期，床層的非均勻反應(yīng)部分是換熱流體溫升的主要原因，衰退階段后期，床層余熱是換熱流體溫升的主要原因。因此，衰退階段的控制因素是床層熱阻（熱導(dǎo)率）。

圖7 衰退階段不同時刻下的床層反應(yīng)分?jǐn)?shù)分布(從左至右時刻依次為9000、9375、9750、10125、10500 s)Fig.7 Distribution of bed reaction fraction at different time in decline stage

圖8 衰退階段不同時刻下的床層溫度分布(從左至右時刻依次為10500、11000、12000、13000、14000 s)Fig.8 Distribution of bed temperature at different time in decline stage

2.5 不同控制因素對反應(yīng)器性能的影響

由圖9可見：隨著進(jìn)口氫壓的增大，換熱流體輸出溫度整體有小幅度升高，熱輸出時間縮短。這是因為進(jìn)口氫壓增加，意味著床層內(nèi)部材料的反應(yīng)平衡溫度升高，反應(yīng)速率也會提高，所以換熱流體輸出溫度會升高，但是床層總的放熱量不變，所以熱輸出時間會縮短。同時，由于床層熱導(dǎo)率不高，換熱能力有限，限制了氫壓增大帶來的反應(yīng)速率變化，所以從整體上看換熱流體輸出溫度增幅不大。進(jìn)一步，由圖10可見，隨著進(jìn)口氫氣壓力的增加，峰值點對應(yīng)溫度升高，對應(yīng)時間縮短，說明氫壓的變化可以影響初始上升階段的反應(yīng)器性能。

圖9 不同壓力下床層平均反應(yīng)分?jǐn)?shù)和換熱流體輸出溫度隨時間的變化曲線Fig.9 Variation curves of average bed reaction fraction and heat transfer fluid output temperature with time under different pressures

圖10 不同壓力下?lián)Q熱流體輸出峰值點參數(shù)隨時間的變化規(guī)律Fig.10 Variation law of heat transfer fluid peak point parameters with time under different pressures

根據(jù)文獻(xiàn)[33]可知，床層的壁面?zhèn)鳠嵯禂?shù)可達(dá)1000～3000 W/(m2·K)。在該區(qū)間內(nèi)選取5 個壁面?zhèn)鳠嵯禂?shù)，研究床層平均反應(yīng)分?jǐn)?shù)和換熱流體輸出溫度的變化規(guī)律，如圖11 所示。表1 中列出了不同壁面?zhèn)鳠嵯禂?shù)下，換熱流體輸出溫度曲線的關(guān)鍵參數(shù)。從圖11 中可以看出，不同的壁面?zhèn)鳠嵯禂?shù)下，床層平均反應(yīng)分?jǐn)?shù)和換熱流體輸出溫度隨時間的變化曲線基本重合，說明壁面?zhèn)鳠嵯禂?shù)對反應(yīng)器性能的影響比較小。從表1 中峰值點對應(yīng)時間、對應(yīng)溫度和有效輸出溫度的變化可以看到，峰值點對應(yīng)時間基本不變，峰值點對應(yīng)溫度和有效輸出溫度隨壁面?zhèn)鳠嵯禂?shù)的增大有輕微增加。這是因為壁面?zhèn)鳠嵯禂?shù)越高，壁面接觸熱阻越小，傳熱效果會越好。增加量不大，說明壁面?zhèn)鳠嵯禂?shù)的變化，即接觸熱阻的變化對整體傳熱性能的影響較小。從表1反應(yīng)鋒面形成時間、形成溫度和初始下降斜率的變化可知，不同壁面?zhèn)鳠嵯禂?shù)對這三個參數(shù)有較大的影響。壁面?zhèn)鳠嵯禂?shù)越大，反應(yīng)鋒面形成時間越短，形成溫度越高，對應(yīng)的初始下降斜率越大。結(jié)合表1其余參數(shù)的變化規(guī)律可知，壁面?zhèn)鳠嵯禂?shù)（即接觸熱阻）對初始下降階段有較大的影響，但是對整體傳熱的影響并不大。

圖11 不同壁面?zhèn)鳠嵯禂?shù)下床層平均反應(yīng)分?jǐn)?shù)和換熱流體輸出溫度隨時間的變化Fig.11 Variation curves of average bed reaction fraction and heat transfer fluid output temperature with time under different wall heat transfer coefficients

從前文的分析可知，在平臺階段和衰退階段的控制因素是床層熱阻，即床層熱導(dǎo)率。圖12為不同床層熱導(dǎo)率下，床層平均反應(yīng)分?jǐn)?shù)和換熱流體輸出溫度隨時間的變化，其關(guān)鍵參數(shù)見表2。從表2中的有效輸出溫度、平臺末端時間和反應(yīng)完成時間的變化趨勢可見，隨著床層熱導(dǎo)率的增加，床層的平臺期溫度（有效輸出溫度）明顯升高，反應(yīng)時間明顯縮短，床層達(dá)到完全反應(yīng)狀態(tài)的時間明顯縮短。這是因為床層熱導(dǎo)率越高，床層熱阻越小，換熱效果越好，反應(yīng)速率越快，換熱流體的輸出溫度會越高。由于床層中材料的總量不變，反應(yīng)速率越快，反應(yīng)持續(xù)時間會越短，達(dá)到完全反應(yīng)的時間也就越短。同時觀察表2中余熱散熱時間和散熱段斜率可以看到，床層熱導(dǎo)率越高，散熱時間越短，散熱段斜率越高，散熱效率越高。綜合可以看出，床層熱導(dǎo)率可以顯著影響平臺階段和衰退階段。

圖12 不同床層熱導(dǎo)率下床層平均反應(yīng)分?jǐn)?shù)和換熱流體輸出溫度隨時間的變化Fig.12 Variation curves of average bed reaction fraction and heat transfer fluid output temperature with time under different bed thermal conductivity

同時，對比圖11可知，相比于壁面?zhèn)鳠嵯禂?shù)（接觸熱阻），床層熱導(dǎo)率（床層熱阻）顯然對反應(yīng)器性能的影響更為劇烈。此外，結(jié)合圖10、表1 和表2 中峰值點對應(yīng)時間和對應(yīng)溫度的變化，可以看到壁面?zhèn)鳠嵯禂?shù)和床層熱導(dǎo)率對峰值點的影響較小，氫壓是影響峰值點參數(shù)，也就是初始上升階段的主要因素。總之，幾何結(jié)構(gòu)不變的情況下，氫壓、接觸熱阻和床層熱導(dǎo)率是控制反應(yīng)器熱輸出能力的三個因素，相比于氫壓的難以調(diào)節(jié)和控制，接觸熱阻的影響較小，床層熱導(dǎo)率是調(diào)節(jié)反應(yīng)器性能的主要因素。

表1 不同壁面?zhèn)鳠嵯禂?shù)下峰值點參數(shù)、反應(yīng)鋒面參數(shù)和輸出溫度Table 1 Peak point，reaction front parameters and output temperature under different wall heat transfer coefficients

表2 不同床層熱導(dǎo)率下峰值點參數(shù)、反應(yīng)器輸出參數(shù)Table 2 Peak point and reactor output parameters under different bed thermal conductivity

3 結(jié)論

本文通過建立金屬氫化物高溫蓄熱反應(yīng)器的多物理場耦合模型，明確了氫熱耦合傳遞機理，為金屬氫化物高溫蓄熱技術(shù)的應(yīng)用提供了理論基礎(chǔ)，主要結(jié)論如下。

（1）金屬氫化物高溫蓄熱反應(yīng)器的熱輸出可以分為三個階段：初始階段、平臺階段和衰退階段。以峰值點為界，將初始階段分為初始上升和初始下降兩個階段。在初始上升階段，氫壓是主要控制因素，當(dāng)氫氣充滿整個床層時，輸出溫度達(dá)到最大值；在初始下降階段，床層與換熱管壁面的接觸熱阻是主要控制因素。在平臺階段和衰退階段，床層熱阻（熱導(dǎo)率）是主要的控制因素。進(jìn)一步，當(dāng)換熱流體進(jìn)口側(cè)完全反應(yīng)時，會出現(xiàn)非均勻反應(yīng)，進(jìn)入衰退階段。在衰退階段，非均勻反應(yīng)和床層余熱是導(dǎo)致?lián)Q熱流體溫升的主要原因。

（2）在初始下降階段，反應(yīng)鋒面逐步形成，在平臺階段，反應(yīng)鋒面向床層外側(cè)移動，并且隨著時間移動速度越來越慢。非均勻反應(yīng)是由換熱流體沿著流體流向傳熱溫差不同所導(dǎo)致的，同時非均勻反應(yīng)是導(dǎo)致?lián)Q熱流體有效輸出時間縮短的主要因素。

（3）平臺階段為反應(yīng)器的有效放熱階段，綜合整個反應(yīng)過程，接觸熱阻對反應(yīng)的影響可以忽略。氫壓、床層熱導(dǎo)率是影響床層熱質(zhì)傳遞的主要因素。增大氫壓和床層熱導(dǎo)率可以使床層快速和高效地吸放熱，但是相比于氫壓，床層熱導(dǎo)率是更加敏感和易于調(diào)節(jié)的因素。在金屬氫化物的高溫蓄熱領(lǐng)域，提高性能的最有效途徑是提高床層熱導(dǎo)率。

符號說明

h——傳熱系數(shù)，W/(m2·K)

p——壓力，MPa

Rg——通用氣體常數(shù)，J/(mol·K)

TC1——床層中心溫度，°C

ε——孔隙率

下角標(biāo)

bed——床層

e——有效值

ea——吸氫平衡態(tài)

g——氫氣

in——入口

out——輸出

s——固體

w——換熱管壁面