删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

流化床膜反应器甘油强化重整的数值模拟

本站小编 哈尔滨工业大学/2020-12-05

流化床膜反应器甘油强化重整的数值模拟

王帅1,2,王会宁3,唐宇翔2,杨学松2

(1. 广东省新能源和可再生能源研究开发与应用重点实验室,广州 510640; 2.哈尔滨工业大学,能源科学与工程学院,哈尔滨 150001; 3. 北京航天动力研究所,北京 100076)



摘要:

为实现生物柴油副产物甘油在制氢行业的发展,以流化床甘油重整制氢为研究对象,基于双流体模型和颗粒动力学理论,结合甘油重整反应动力学模型,并嵌入二氧化碳吸附动力学模型和氢气膜分离模型来描述两种强化重整方法的作用. 对流化床反应器生物甘油强化重整制氢过程开展了数值模拟,对反应器内颗粒浓度、组分浓度、温度进行预测,探究重整过程中气固两相流动与反应特性,分析氢气膜分离和二氧化碳吸附两种强化重整方法的相互作用规律,评价操作参数对重整性能的影响. 结果表明:二氧化碳吸附可以抑制浓度极化阻力,提高氢气渗透速率;吸附剂与催化剂比例为1∶1时,与没有吸附剂相比,氢气相对产量提高了5%;氢气分离膜厚的减少会进一步提高二氧化碳吸附速率,当膜厚从300 μm减少到30 μm时,吸附速率提高1.4%;催化-吸附双功能颗粒的使用可以加强二氧化碳的吸附水平,同时促进氢气分离过程,相较于无吸附强化,氢气渗透量提高了近20%.

关键词:  甘油重整  氢气分离  二氧化碳吸附  流化床  膜分离

DOI:10.11918/201910087

分类号:TK91

文献标识码:A

基金项目:广东省新能源和可再生能源研究开发与应用重点实验室资助(Y909kq1001)



Numerical simulation of enhanced reforming of glycerol in a membrane-assisted fluidized bed

WANG Shuai1,2,WANG Huining3,TANG Yuxiang2,YANG Xuesong2

(1.Guangdong Provincial Key Laboratory of New and Renewable Energy Research and Development, Guangzhou 510640, China; 2. School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001,China; 3. Beijing Aerospace Propulsion Institute, Beijing 100076, China)

Abstract:

To achieve the development of glycerol as biodiesel byproduct in the industry of hydrogen production, the glycerol reforming process for hydrogen production in a membrane-assisted fluidized bed reactor is numerically simulated on the basis of the two-fluid model and the kinetic theory of granular flow coupled with the glycerol reforming kinetic model, where the CO2 sorption kinetic model and hydrogen separation model are implemented. The gas component and particle concentrations as well as temperature are predicted, and the multiphase flow behaviors and reaction characteristics during the reforming process are evaluated. The mutual interaction mechanism of the two enhancing methods including membrane hydrogen separation and carbon dioxide sorption is discussed, and the impact of operating parameters on reforming performance is examined. The result reveals that the concentration polarization resistance will be restricted with the rising hydrogen permeation. At the sorbent to catalyst ratio of 1∶1, the relative hydrogen yield is improved by 5% compared to the reforming process without sorbents. When the membrane thickness is reduced from 300 μm to 30 μm, the CO2 sorption rate can be increased by 1.4%. The utilization of catalyst-sorbent bi-functional particles can enhance CO2 sorption and hydrogen separation. The hydrogen permeation is improved by almost 20%.

Key words:  glycerol reforming  hydrogen separation  carbon dioxide sorption  fluidized bed  membrane separation


王帅, 王会宁, 唐宇翔, 杨学松. 流化床膜反应器甘油强化重整的数值模拟[J]. 哈尔滨工业大学学报, 2020, 52(7): 89-95. DOI: 10.11918/201910087.
WANG Shuai, WANG Huining, TANG Yuxiang, YANG Xuesong. Numerical simulation of enhanced reforming of glycerol in a membrane-assisted fluidized bed[J]. Journal of Harbin Institute of Technology, 2020, 52(7): 89-95. DOI: 10.11918/201910087.
基金项目 广东省新能源和可再生能源研究开发与应用重点实验室资助(Y909kq1001) 作者简介 王帅(1985—),男,副教授,博士生导师 通信作者 王帅,shuaiwang@hit.edu.cn 文章历史 收稿日期: 2019-10-14



Abstract            Full text            Figures/Tables            PDF


流化床膜反应器甘油强化重整的数值模拟
王帅1,2, 王会宁3, 唐宇翔2, 杨学松2    
1. 广东省新能源和可再生能源研究开发与应用重点实验室, 广州 510640;
2. 哈尔滨工业大学, 能源科学与工程学院, 哈尔滨 150001;
3. 北京航天动力研究所, 北京 100076

收稿日期: 2019-10-14
基金项目: 广东省新能源和可再生能源研究开发与应用重点实验室资助(Y909kq1001)
作者简介: 王帅(1985—),男,副教授,博士生导师
通信作者: 王帅,shuaiwang@hit.edu.cn


摘要: 为实现生物柴油副产物甘油在制氢行业的发展, 以流化床甘油重整制氢为研究对象, 基于双流体模型和颗粒动力学理论, 结合甘油重整反应动力学模型, 并嵌入二氧化碳吸附动力学模型和氢气膜分离模型来描述两种强化重整方法的作用.对流化床反应器生物甘油强化重整制氢过程开展了数值模拟, 对反应器内颗粒浓度、组分浓度、温度进行预测, 探究重整过程中气固两相流动与反应特性, 分析氢气膜分离和二氧化碳吸附两种强化重整方法的相互作用规律, 评价操作参数对重整性能的影响.结果表明:二氧化碳吸附可以抑制浓度极化阻力, 提高氢气渗透速率; 吸附剂与催化剂比例为1 : 1时, 与没有吸附剂相比, 氢气相对产量提高了5%;氢气分离膜厚的减少会进一步提高二氧化碳吸附速率, 当膜厚从300 μm减少到30 μm时, 吸附速率提高1.4%;催化-吸附双功能颗粒的使用可以加强二氧化碳的吸附水平, 同时促进氢气分离过程, 相较于无吸附强化, 氢气渗透量提高了近20%.
关键词: 甘油重整    氢气分离    二氧化碳吸附    流化床    膜分离    
Numerical simulation of enhanced reforming of glycerol in a membrane-assisted fluidized bed
WANG Shuai1,2, WANG Huining3, TANG Yuxiang2, YANG Xuesong2    
1. Guangdong Provincial Key Laboratory of New and Renewable Energy Research and Development, Guangzhou 510640, China;
2. School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, China;
3. Beijing Aerospace Propulsion Institute, Beijing 100076, China



Abstract: To achieve the development of glycerol as biodiesel byproduct in the industry of hydrogen production, the glycerol reforming process for hydrogen production in a membrane-assisted fluidized bed reactor is numerically simulated on the basis of the two-fluid model and the kinetic theory of granular flow coupled with the glycerol reforming kinetic model, where the CO2 sorption kinetic model and hydrogen separation model are implemented. The gas component and particle concentrations as well as temperature are predicted, and the multiphase flow behaviors and reaction characteristics during the reforming process are evaluated. The mutual interaction mechanism of the two enhancing methods including membrane hydrogen separation and carbon dioxide sorption is discussed, and the impact of operating parameters on reforming performance is examined. The result reveals that the concentration polarization resistance will be restricted with the rising hydrogen permeation. At the sorbent to catalyst ratio of 1 : 1, the relative hydrogen yield is improved by 5% compared to the reforming process without sorbents. When the membrane thickness is reduced from 300 μm to 30 μm, the CO2 sorption rate can be increased by 1.4%. The utilization of catalyst-sorbent bi-functional particles can enhance CO2 sorption and hydrogen separation. The hydrogen permeation is improved by almost 20%.
Keywords: glycerol reforming    hydrogen separation    carbon dioxide sorption    fluidized bed    membrane separation    
氢能被誉为21世纪的清洁高效能源,具有高热值、燃烧性能好、清洁无污染等优点,在世界各国加速能源转型的战略背景下,具有巨大的发展潜力[1].催化重整是使用最广泛的制氢方式之一,然而, 传统的固定床反应器因受热不均,颗粒表面易形成积碳,进而导致燃料转化率和氢气产量的下降.流化床反应器具有良好的传质、传热性能和连续运行潜力,在制氢工艺中逐渐受到了研究者们的关注[2].

为了进一步实现催化重整过程的强化,打破原有热力学平衡的限制,利用渗透膜将氢气从反应体系中分离,可大幅度提高氢气的产率,同时可以实现高纯度氢气的生产[3]. Rakib等[4]对流化床膜反应器内丙烷蒸汽重整制氢进行了实验研究,结果表明,与传统反应器相比,流化床膜反应器在较低的温度下可以获得相同的产氢率,有效地抑制了催化剂烧结和高温焦炭的形成. Mahecha-Botero等[5]利用实验室规模的流化床膜反应器对天然气蒸汽重整过程进行了研究,结果表明,减小渗透侧氢气分压力、增大膜面积和减小膜厚都有利于工艺整体效率的提升. Chen等[6]通过实验手段评价了反应侧压力对鼓泡流化床膜反应器内甲烷蒸汽重整和自热重整过程的反应性能的影响,认为较高的压力并不总能提高产氢率,因为较高的氢渗透率被较低的平衡转化率所抵消.二氧化碳吸附,作为另一种强化重整的有效手段,通过添加固体吸附剂捕获反应体系的二氧化碳,突破热力学限制,从而提高了重整过程中氢气产量和浓度,同时吸附反应是放热反应,为吸热的重整反应减小了外界能量输入[7]. Lindborg等[8]通过模拟手段,比较了有无二氧化碳吸附的鼓泡流化床反应器中甲烷重整反应性能,结果表明,二氧化碳吸附可以显著提高氢气浓度. Hildenbrand等[9]利用实验手段研究流化床反应器氧化钙吸附强化甲烷蒸汽重整过程.发现吸附剂多次循环再生后,氢氧化钙的形成会改变实际的水碳比,导致甲烷转化率和氢气产率降低.

上述研究表明,膜分离和二氧化碳吸附两种强化重整手段通过打破热力学平衡,可以进一步提高燃料转化率和氢气产率. Silva等[10]通过热力学分析,研究了氢气分离与二氧化碳吸附一体化对甘油蒸汽重整制氢过程的影响,发现二者结合不仅大大提高了氢气产率,而且吸附反应释放的热量可以使重整系统达到自热状态.然而,将两种强化技术同时应用于反应器重整制氢过程的研究较少,开展膜分离与二氧化碳吸附协同强化重整技术,探究两者之间的相互作用规律,对反应器的设计和过程的调控具有重要的意义.

本文基于双流体模型,结合化学反应动力学理论,对流化床反应器内膜分离与二氧化碳吸附协同强化甘油重整制氢过程开展了模拟研究,获得了流化床重整制氢系统中多相流动和反应行为,探究了二氧化碳吸附和氢气膜分离过程的相互作用规律,比较了不同条件下吸附速率与膜渗透率之间的相互影响.

1 数学模型 1.1 控制方程计算模拟采取欧拉-欧拉双流体模型,该模型基于连续性假设,物质被认定为在其所占据的空间连续分布.计算时采用如下的基本守恒方程.

质量守恒方程为

$\begin{array}{l}\frac{\partial }{{\partial t}}({\varepsilon _{\rm g}}{\rho _{\rm g}}) + \nabla {\rm{\cdot(}}{\varepsilon _{\rm g}}{\rho _{\rm g}}{u_{\rm g}}) = {S_{\rm gm}} + {S_{\rm{mass, {H_2}}}}, \\\frac{\partial }{{\partial t}}({\varepsilon _{\rm m}}{\rho _{\rm m}}) + \nabla {\rm{\cdot(}}{\varepsilon _{\rm m}}{\rho _{\rm m}}{u_{\rm m}}) = {S_{\rm mg}}.\end{array}$ (1)

式中:SgmSmg为由反应引起的质量源项,Smass, H2为氢气分离引起的质量源项,u为速度,ρε分别为密度和体积分数,下角标g表示气相,m表示固相的第m相.

动量守恒方程:

$\begin{array}{l}[\frac{\partial }{{\partial t}}({\varepsilon _{\rm g}}{\rho _{\rm g}}{u_{\rm g}}) + \nabla {\rm{\cdot(}}{\varepsilon _{\rm g}}{\rho _{\rm g}}{u_{\rm g}}{u_{\rm g}})] = - {\varepsilon _{\rm g}}\nabla {p_{\rm g}} + {\varepsilon _{\rm g}}\nabla {\rm{\cdot}}{\tau _{\rm g}} + {\varepsilon _{\rm g}}{\rho _{\rm g}}g - {\beta _{\rm gm}}({u_{\rm g}} - {u_{\rm m}}) + {S_{\rm gm}}{u_{\rm g}} + {S_{\rm{mom, {H_2}}}}, \\[\frac{\partial }{{\partial t}}({\varepsilon _{\rm m}}{\rho _{\rm m}}{u_{\rm m}}) + \nabla {\rm{\cdot(}}{\varepsilon _{\rm m}}{\rho _{\rm m}}{u_{\rm m}}{u_{\rm m}})] = - {\varepsilon _{\rm m}}\nabla {p_{\rm g}} - \nabla {p_{\rm m}} + {\varepsilon _{\rm m}}\nabla {\rm{\cdot}}{\tau _{\rm m}} + {\varepsilon _{\rm m}}{\rho _{\rm m}}g + {\beta _{\rm gm}}({u_{\rm g}} - {u_{\rm m}}) + {\beta _{\rm mn}}({u_{\rm m}} - {u_n}) + {S_{\rm mg}}{u_{\rm m}}.\end{array}$ (2)

式中:p为压力,β为相间曳力系数,τ为剪切应力项,Smom, H2为氢气分离引起的动量源项.

能量守恒方程为

$\begin{array}{l}\frac{\partial }{{\partial t}}({\varepsilon _{\rm g}}{\rho _{\rm g}}{c_{\rm pg}}{T_{\rm g}}) + \nabla {\rm{\cdot(}}{\varepsilon _{\rm g}}{\rho _{\rm g}}{c_{\rm pg}}{T_{\rm g}}{u_{\rm g}}) = \nabla [{\varepsilon _{\rm g}}({\lambda _{\rm g}} + {c_{\rm pg}}\frac{{{\mu _{\rm gt}}}}{{{{\Pr }_t}}})\nabla {T_{\rm g}}] + \sum {{h_{\rm gm}}} ({T_{\rm m}} - {T_{\rm g}}) + \sum {{S_{\rm {g, i}}}{c_{pg, i}}} {T_{\rm g}} + {S_{en, {H_2}}}, \\\frac{\partial }{{\partial t}}({\varepsilon _{\rm m}}{\rho _{\rm m}}{c_{\rm pm}}{T_{\rm m}}) + \nabla {\rm{\cdot}}({\varepsilon _{\rm m}}{\rho _{\rm m}}{c_{\rm pm}}{T_{\rm m}}{u_{\rm m}}) = \nabla [{\varepsilon _{\rm m}}({\lambda _{\rm m}} + {c_{\rm pm}}\frac{{{\mu _{\rm st}}}}{{{{\Pr }_t}}})\nabla {T_{\rm g}}] + {h_{\rm gm}}({T_{\rm g}} - {T_{\rm m}}) + \sum {{S_{\rm {m, i}}}{c_{pm, i}}{T_{\rm m}}.} \end{array}$ (3)

式中:cp为比定压热容,λ为导热系数,μ为黏度,hgm为气固相间换热系数,Sen, H2为氢气分离引起的能量源项.

组分输运方程为

$\begin{array}{l}\frac{\partial }{{\partial t}}({\varepsilon _{\rm g}}{\rho _{\rm g}}{Y_{\rm {g, i}}}) + \nabla {\rm{\cdot}}({\varepsilon _{\rm g}}{\rho _{\rm g}}{Y_{\rm {g, i}}}{u_{\rm g}}) = \nabla {\rm{\cdot}}[{\varepsilon _{\rm g}}({\rho _{\rm g}}{D_{\rm {g, i}}} + \frac{{{\mu _{\rm gt}}}}{{S{C_t}}})\nabla {\rm{\cdot}}{Y_{\rm {g, i}}}] + {S_{\rm {g, i}}} + {S_{mass, {H_2}}}, \\\frac{\partial }{{\partial t}}({\varepsilon _{\rm m}}{\rho _{\rm m}}{Y_{\rm {m, i}}}) + \nabla {\rm{\cdot(}}{\varepsilon _{\rm m}}{\rho _{\rm m}}{Y_{\rm {m, i}}}{u_{\rm m}}{\rm{) = }}\nabla {\rm{\cdot}}[{\varepsilon _{\rm m}}({\rho _{\rm m}}{D_{\rm {m, i}}} + \frac{{{\mu _{\rm st}}}}{{S{C_t}}})\nabla {\rm{\cdot}}{Y_{\rm {m, i}}}] + {S_{\rm {m, i}}}.\end{array}$ (4)

式中:Y为质量分数,D为扩散系数.

氢气膜分离过程是通过在流化床边壁处的一层网格中加入源项来考虑的,对于气相的连续性方程和氢气组分守恒方程为

$\phi = - \frac{A}{V}\;\frac{{{P_{\rm e}}}}{l}\exp \left( { - \frac{{{E_a}}}{T}} \right)(P_{\rm{retentate, {H_2}}}^{0.5} - P_{\rm{permeate, {H_2}}}^{0.5}).$ (5)

式中:PeEal分别为指前因子、活化能和膜厚, 其值分别为1.12×10-5 mol·s-1·-1·Pa-0.5、29.16 kJ·mol-1和30 μm, pretentate, H2ppermeate, H2分别为反应器侧和渗透侧的氢气分压;AV为单元网格的面积和体积.

1.2 反应动力学模型二氧化碳吸附强化甘油水蒸气重整制氢过程由一系列复杂的反应组成,其气相组分主要包括C3H8O3、H2O、H2、CO2、CO和CH4.本研究操作条件(蒸汽甘油比为6,且温度为873 K)下,可忽略积碳生成[11].主要反应如下:

$\begin{array}{l}{{\rm{C}}_{\rm{3}}}{{\rm{H}}_{\rm{8}}}{{\rm{O}}_{\rm{3}}}{\rm{(g) + }}{{\rm{H}}_{\rm{2}}}{\rm{O(g)}} \to {\rm{C}}{{\rm{H}}_{\rm{4}}}{\rm{(g) + 2C}}{{\rm{O}}_{\rm{2}}}{\rm{(g) + 3}}{{\rm{H}}_{\rm{2}}}{\rm{(g), }}\\{\rm{C}}{{\rm{H}}_{\rm{4}}}{\rm{(g) + }}{{\rm{H}}_{\rm{2}}}{\rm{O(g)}} \leftrightarrow {\rm{CO(g) + 3}}{{\rm{H}}_{\rm{2}}}{\rm{(g) , }}\\{\rm{CO(g) + }}{{\rm{H}}_{\rm{2}}}{\rm{O(g)}} \leftrightarrow {\rm{C}}{{\rm{O}}_{\rm{2}}}{\rm{(g) + }}{{\rm{H}}_{\rm{2}}}{\rm{(g), }}\\{\rm{CaO(s) + C}}{{\rm{O}}_{\rm{2}}}{\rm{(g)}} \to {\rm{CaC}}{{\rm{O}}_{\rm{3}}}{\rm{(s)}}{\rm{.}}\end{array}$ (6)

采用Dou等[12]提出的层流有限元模型计算反应速率,并假设速率常数随温度的变化而变化.甘油重整和吸附反应的速率表达式如下:

$\begin{array}{l}{R_1} = 1.838 \times {10^5}\exp ( - 74\;210/T){C_{\rm {C_3}{H_8}{O_3}}}{C_{\rm {H_2}O}}{\rho _{\rm cat}},\\{R_2} = 1.198 \times {10^{17}}\exp ( - 26\;830/T){C_{\rm C{H_4}}}{C_{\rm {H_2}O}}{\rho _{\rm cat}},\\{R_3} = 0.017\;67\exp (4\;400/T){C_{\rm CO}}{C_{\rm {H_2}O}}{\rho _{\rm cat}},\\{R_4} = {C_{\rm CaO}}{k_c}{({C_{\rm C{O_2}}} - {C_{\rm C{O_2},eq}})^{0.37}}{(1 - X/{X_u})^{2.61}}.\end{array}$ (7)

式中:CXCCO2, eqkc分别为浓度、转化率、CO2平衡浓度和动力学系数,其具体表达式可见Iliuta等[13]的研究.对于双功能催化剂-吸附剂,反应速率可表示为[14]

${R_i} = {v_{{R_1},i}}{R_1} + {v_{{R_2},i}}{R_2} + {v_{{R_3},i}}{R_3} + {v_{{R_{CO}}_{_2},i}}{R_4}.$ (8)

式中:vR1, ivR2, ivR3, ivR4, i代表组分i的化学计量数(i =C3H8O3, H2O, H2, CO2, CO, CH4).

1.3 几何模型及计算参数模拟采用的流化床膜反应器,其高度为1.5 m,直径为0.3 m,结构示意图见图 1.

Fig. 1
图 1 流化床膜反应器的结构示意 Fig. 1 Structure of a fluidized bed membrane reactor


本文采用两种吸附剂设计(传统两颗粒和双功能颗粒设计),以探究吸附对重整性能的强化作用.双功能颗粒的催化组分和吸附组分与传统两颗粒设计一致.在无吸附重整过程中,选择惰性颗粒代替吸附剂颗粒,保证催化剂用量和总颗粒数保持不变.此模拟计算所采用的初始条件、边界条件以及相关的主要物性参数见表 1.

表 1
表 1 主要模拟参数 Tab. 1 Main simulation parameters 反应器
高度,
直径/m 催化剂, 吸附剂,
惰性, 双功能颗粒直径/
μm 催化剂, 吸附剂,
惰性, 双功能颗粒密度/
(kg·m-3) 初始催化剂-
吸附剂
体积比 初始颗粒
堆积
高度/m 颗粒相
初始
体积分数 渗透膜
高度/
m 入口
气体速度/
(m· s-1)

1.5, 0.3 150,150,150,150 2650,3320,3320,2784 4:1 0.25 0.53 0.5 0.5



C3H8O3,
H2O, N2入口份额 水蒸汽密度/
(kg·m-3) 水蒸汽黏度/
(Pa·s) 水蒸汽比热容/
(kJ·kg-1·K-1) 出口压力/
MPa 入口温度/
K 渗透侧氢气分压力/
MPa

0.02,0.12,0.86 0.25 3.05×10-5 2.2 0.1 873 0.001



表 1 主要模拟参数 Tab. 1 Main simulation parameters


2 结果与讨论为了验证甘油重整和氢气分离模型,进行了出口气体组分和氢气渗透速率的模拟,结果见图 2,并与Dou等[12]和Helmi等[15]的实验结果进行比较.由图 2可知,模型预测的出口组分摩尔分数和氢气渗透速率与实验值之间的最大相对误差不超过10%,在可接受的范围内.另外,经过网格无关性验证,采用了32 000个元素的网格.

Fig. 2
图 2 实验与模拟数据对比 Fig. 2 Comparison of simulation and measured data


有、无吸附条件下,流化床内时均空隙率分布见图 3.由图 3可知,由于吸附剂、催化剂和惰性颗粒相似的物理性质,3种操作条件下床层空隙率分布和床层膨胀高度基本相同.同时可观察到,渗透膜附近存在局部的低空隙率区域,这是由于氢气分离和壁面摩擦抑制了颗粒流化,形成颗粒致密区域.同时发现使用双功能颗粒设计,颗粒致密区域范围扩大.

Fig. 3
图 3 时均空隙率分布云图 Fig. 3 Contour plots of time-averaged voidage


有、无吸附条件下,墙壁附近和中间区域两个监测点的空隙率随时间变化如图 4所示.

Fig. 4
图 4 壁面和中心检测点空隙率随时间的变化 Fig. 4 Variation of voidage with time at wall and center of reactor


由图 4可知,中心区域比墙壁附近的空隙率波动更加剧烈.这是由于致密区域的形成增加了墙壁附近气泡运动的阻力,使气泡运动主要集中在床层中间区域,从而导致了明显的空隙率振荡.与无二氧化碳吸附相比,二氧化碳吸附进一步增强了墙壁附近的空隙率波动程度.尤其是采用双功能吸附剂时,二氧化碳吸附的作用效果更加明显.这是由于二氧化碳吸附强化了催化重整反应,导致流化床气体量增大,同时扩大了致密区域的范围,导致更多的气泡沿着流化床中间区域上升.此外,二氧化碳的吸附也增大了墙壁附近的空隙率波动程度.这是由于重叠的催化剂和吸附剂局部高浓度区域形成导致了墙壁附近的气体产物也有所增加,从而加剧了墙壁附近的空隙率波动.

有、无二氧化碳吸附条件下,催化剂和吸附剂颗粒时均体积分数分布如图 5所示.氢气分离导致流化床膜周围形成催化剂和吸附剂(或惰性颗粒)的局部高浓度区.一方面由于墙壁处颗粒易受到摩擦;另一方面,此区氢气分离作用引起颗粒流化速度降低,两者共同导致墙壁处局部颗粒高浓度区.另外,二氧化碳吸附也明显导致该区域范围扩大,进一步验证了吸附对氢气分离的强化作用.虽然致密化区的形成不利于流化床床层和渗透膜间传质,且减小了气固相接触和气体停留时间,不利于流化床膜反应器性能.但此重叠高浓度区域内催化剂和吸附剂高度集中,促进催化剂和吸附剂间传质,一定程度上有利于吸附强化催化重整反应.

Fig. 5
图 5 颗粒时均体积分数分布 Fig. 5 Contour plots of time-averaged solid volume fraction


气相组分的时均质量分数分布如图 6所示.由于反应器底部的催化剂和吸附剂的浓度较高,重整反应剧烈,其组分变化较为明显.氢气分离通过直接提高反应物浓度和形成局部催化剂颗粒高浓度区域,促进了渗透膜周围甘油分解和产物生成,导致膜周围形成甘油局部低浓度区域、氢气和二氧化碳局部高浓度区.两种吸附设计条件下吸附作用是相似的,对重整过程生成的二氧化碳都实现了充分吸附,从而进一步强化了重整反应,扩大了局部氢气高浓度区域和局部甘油低浓度区域范围.尤其是采用双功能颗粒设计时,这一强化效果更加明显.

Fig. 6
图 6 气相组分时均质量分数分布 Fig. 6 Contour plots of time-averaged mass fraction of gas composition


二氧化碳吸附对氢气渗透的影响如图 7所示.由图 7可知,有无吸附条件下的氢气渗透速率的轴向分布趋势是相似的,都随着轴向高度升高而增大.这是由于随着重整反应进行,渗透膜表面氢气分压力沿轴向增大.图中同时比较了有无吸附条件下的总氢气渗透流量.由图 7可知,添加二氧化碳吸附剂大大提高了氢气渗透流量.这是因为二氧化碳的吸附强化了反应以及氢气的产量,进而提高了膜附近氢气的分压力,导致氢气的渗透量增加.使用双功能吸附剂时,二氧化碳吸附对氢渗透的增强作用更为明显.由于双功能颗粒时二氧化碳传质阻力较小,其二氧化碳吸附的作用效果更加明显.因此可以看到,吸附反应不仅强化了重整,而且在一定程度上促进了氢气的分离过程,当使用双功能颗粒时,相比于无吸附,氢气的渗透量提高了近20%.

Fig. 7
图 7 有无吸附条件下氢气渗透流率的分布 Fig. 7 Distribution of hydrogen permeation rate with or without sorption


有、无二氧化碳吸附条件下,不同轴向高度处的气相温度径向分布如图 8所示.由图 8可知,无二氧化碳吸附条件下,催化重整反应吸收大量热量,不同轴向高度处的气相温度明显低于入口温度.渗透膜周围的温度存在明显降低趋势,这是由于氢气分离强化了重整反应, 导致此区域的温度降低程度更加明显.而加入二氧化碳吸附后,吸附反应释放热量,导致不同轴向高度处的气相温度都有所升高.并且气相温度随着轴向高度的升高而增大,这说明吸附反应可以提供重整反应所需要的热量,进而实现了系统的自热.

Fig. 8
图 8 有无吸附条件下不同高度处气相温度径向分布 Fig. 8 Lateral distribution of gas temperature at different heights with or without sorption


渗透膜厚度对燃料转化率和二氧化碳吸附速率的影响如图 9所示.由图 9可知,甘油转化率和吸附速率随着膜厚的减小而升高.这是由于随着膜厚的减小,氢气分离对重整制氢的强化作用增大,提高了燃料转化率.同时氢气分离导致其他气体的浓度升高的程度也有所增大.两者共同导致二氧化碳浓度升高,从而引起二氧化碳吸附速率的增大.当膜厚从300 μm减少到30 μm时,膜分离性能提高的同时,二氧化碳吸附速率可以提高1.4%,提升的程度不是很大,这是因为本文采用的是外置膜流化床反应器,而当膜管采用内置膜布置在床层中时,随着膜管数的增多,膜分离与二氧化碳吸附的相互作用会更强烈.

Fig. 9
图 9 渗透膜厚度对燃料转化率和二氧化碳吸附速率的影响 Fig. 9 Effect of membrane thickness on fuel conversion and carbon dioxide sorption rate


图 10(a)给出了吸附剂比例对氢气产量的影响.随着吸附剂比例的增大,氢气相对产量增大.这是由于随着吸附剂比例的增大,二氧化碳吸附对催化重整的强化作用也更加明显;同时吸附反应释放的热量随之增加,升高了床层温度,两者共同导致氢气相对产量升高.当吸附剂与催化剂比例为1:1时,氢气相对产量与没有吸附剂相比,提高了5%.对二氧化碳吸附强化甘油重整进行热力学分析[16]时发现,二氧化碳吸附大幅度提高了氢气产量,每摩尔甘油可以生成接近7 mol的氢气.而相较于热力学平衡,流化床反应器中由于颗粒与燃料的接触时间的限制,进而影响了甘油的转化与氢气的生成,因此有必要对过程实施调控,以更好地提高燃料转化率和氢气产率.

Fig. 10
图 10 吸附剂比例对氢气产量及浓度极化系数的影响 Fig. 10 Effect of sorbent to catalyst ratio on hydrogen yield and EAC


为了描述浓度极化的传质限制,本文引入了有效平均浓度极化系数(EAC)评价浓度极化水平[17]. 图 10(b)给出了吸附剂比例对浓度极化系数的影响.由图 10(b)可知,随着吸附剂比例的增大,渗透膜反应器的浓度极化水平减小.由于随着吸附剂比例的增大,膜周围的局部吸附剂颗粒高浓度区域范围也增大,导致二氧化碳吸附对催化重整反应强化程度增大,同时二氧化碳吸附程度的增大也导致包括氢气在内的其他气体浓度升高程度有所增加,两者都导致渗透膜表面氢气浓度梯度减小,从而引起反应器的浓度极化水平不断减小.

3 结论1) 基于双流体模型和反应动力学模型,探究了流化床反应器内二氧化碳吸附和氢气分离两种强化重整方式的相互作用规律,二氧化碳吸附可以抑制浓度极化阻力,提高氢气渗透速率,吸附剂与催化剂比例为1:1时,氢气相对产量与没有吸附剂相比,提高了5%.

2) 氢气分离通过直接提高二氧化碳浓度和强化甘油催化重整反应,提高了二氧化碳吸附速率,当膜厚从300 μm减少到30 μm时,膜分离性能提高的同时,二氧化碳吸附速率可以提高1.4%.

3) 与传统双颗粒设计相比,催化-吸附双功能颗粒的设计,使二氧化碳传质阻力较小,其二氧化碳吸附的作用效果更加明显,当使用双功能颗粒时,相比于无吸附,氢气的渗透量提高了近20%.


参考文献
[1] NIKOLAIDIS P, POULLIKKAS A. A comparative overview of hydrogen production processes[J]. Renewable and Sustainable Energy Reviews, 2017, 67: 597. DOI:10.1016/j.rser.2016.09.044


[2] DURáN P, SANZ-MARTíNEZ A, SOLER J, et al. Pure hydrogen from biogas: intensified methane dry reforming in a two-zone fluidized bed reactor using permselective membranes[J]. Chemical Engineering Journal, 2019, 370: 772. DOI:10.1016/j.cej.2019.03.199


[3] JOKAR S M, PARVASI P, BASILE A. The evaluation of methane mixed reforming reaction in an industrial membrane reformer for hydrogen production[J]. International Journal of Hydrogen Energy, 2018, 43(32): 15321. DOI:10.1016/j.ijhydene.2018.06.142


[4] RAKIB M A, GRACE J R, LIM C J, et al. Steam reforming of propane in a fluidized bed membrane reactor for hydrogen production[J]. International Journal of Hydrogen Energy, 2010, 35(12): 6276. DOI:10.1016/j.ijhydene.2010.03.136


[5] MAHECHA-BOTERO A, BOYD T, GULAMHUSEIN A, et al. Catalytic reforming of natural gas for hydrogen production in a pilot fluidized-bed membrane reactor: mapping of operating and feed conditions[J]. International Journal of Hydrogen Energy, 2011, 36(17): 10727. DOI:10.1016/j.ijhydene.2011.05.178


[6] CHEN Zhongxiang, GRACE J R, LIM C J, et al. Experimental studies of pure hydrogen production in a commercialized fluidized-bed membrane reactor with SMR and ATR catalysts[J]. International Journal of Hydrogen Energy, 2007, 32(13): 2359. DOI:10.1016/j.ijhydene.2007.02.036


[7] TAN L, ROGHAIR I, VAN SINT ANNALAND M. Simulation study on the effect of gas permeation on the hydrodynamic characteristics of membrane-assisted micro fluidized beds[J]. Applied Mathematical Modelling, 2014, 38(17/18): 4291. DOI:10.1016/j.apm.2014.04.044


[8] LINDBORG H, JAKOBSEN H A. Sorption enhanced steam methane reforming process performance and bubbling fluidized bed reactor design analysis by use of a two-fluid model[J]. Industrial & Engineering Chemistry Research, 2009, 48(3): 1332. DOI:10.1021/ie800522p


[9] HILDENBRAND N, READMAN J, DAHL I M, et al. Sorbent enhanced steam reforming (SESR) of methane using dolomite as internal carbon dioxide absorbent: limitations due to Ca(OH)2 formation[J]. Applied Catalysis A: General, 2006, 303(1): 131. DOI:10.1016/j.apcata.2006.02.015


[10] SILVA J M, SORIA M A, MADEIRA L M. Thermodynamic analysis of glycerol steam reforming for hydrogen production with in situ hydrogen and carbon dioxide separation[J]. Journal of Power Sources, 2015, 273: 423. DOI:10.1016/j.jpowsour.2014.09.093


[11] AUTHAYANUN S, ARPORNWICHANOP A, PATCHARAVORACHOT Y, et al. Hydrogen production from glycerol steam reforming for low-and high-temperature PEMFCs[J]. International Journal of Hydrogen Energy, 2011, 36(1): 267. DOI:10.1016/j.ijhydene.2010.10.061


[12] DOU Binlin, SONG Yongchen. A CFD approach on simulation of hydrogen production from steam reforming of glycerol in a fluidized bed reactor[J]. International Journal of Hydrogen Energy, 2010, 35(19): 10271. DOI:10.1016/j.ijhydene.2010.07.165


[13] ILIUTA I, RADFARNIA H R, ILIUTA M C. Hydrogen production by sorption-enhanced steam glycerol reforming: sorption kinetics and reactor simulation[J]. AIChE Journal, 2013, 59(6): 2105. DOI:10.1002/aic.13979


[14] ALOISI I, JAND N, STENDARDO S, et al. Hydrogen by sorption enhanced methane reforming: a grain model to study the behavior of bi-functional sorbent-catalyst particles[J]. Chemical Engineering Science, 2016, 149: 22. DOI:10.1016/j.ces.2016.03.042


[15] HELMI A, VONCKEN R J W, RAIJMAKERS A J. On concentration polarization in fluidized bed membrane reactors[J]. Chemical Engineering Journal, 2018, 332: 464. DOI:10.1016/j.cej.2017.09.045


[16] 王帅, 王琦, 宋晓皎, 等. 生物甘油水蒸气重整制氢强化过程的参数评估[J]. 哈尔滨工业大学学报, 2018, 50(1): 102.
WANG Shuai, WANG Qi, SONG Xiaojiao, et al. Parameter evaluation of enhancing glycerol steam reforming for hydrogen production[J]. Journal of Harbin institute of Technology, 2018, 50(1): 102. DOI:10.11918/j.issn.0367-6234.201607119


[17] CARAVELLA A, SUN Yu. Correct evaluation of the effective concentration polarization influence in membrane-assisted devices[J]. International Journal of Hydrogen Energy, 2016, 41(27): 11653. DOI:10.1016/j.ijhydene.2015.12.068



相关话题/过程 新能源

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 热轧过程中高纯钴微观组织及织构演变
    热轧过程中高纯钴微观组织及织构演变李震1,宋克睿1,韩彦鹏1,肖柱1,贺昕2,陈志永1(1.中南大学材料科学与工程学院,长沙410083;2.北京有色金属研究总院有研亿金新材料股份有限公司,北京102200)摘要:金属钴具有同素异构转变特性。为探究热轧工艺对高纯钴的微观组织及织构演变规律的影响,对纯 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • AZ31B镁合金带材热轧过程组织均匀性及性能研究
    AZ31B镁合金带材热轧过程组织均匀性及性能研究曹东东1,2,梅瑞斌1,2,包立1,侯铮1,黄芸1(1.东北大学秦皇岛分校资源与材料学院,河北秦皇岛066004;2.东北大学材料科学与工程学院,沈阳110819)摘要:本文开展了变形温度为300、350、400℃和总压下率分别为15%、30%、45% ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 考虑高温水泄漏影响的凝水泵启动过程仿真分析
    考虑高温水泄漏影响的凝水泵启动过程仿真分析覃海波,金家善,倪何(海军工程大学动力工程学院,武汉430033)摘要:针对某船用给水机组凝水泵在备用状态和备用转换启动过程中,由于高温除氧水进入凝水泵及其吸入管路而在机组启动后发生汽蚀,导致凝水泵出口压力长时间达不到规定要求的问题,在基于差异演化算法辨识得 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 污水生物处理过程中溶解性有机氮分布和转化特征
    污水生物处理过程中溶解性有机氮分布和转化特征王小东1,陈明飞1,王子文1,王燕1,王硕1,2,3,李激1,2,3(1.江南大学环境与土木工程学院,江苏无锡214122;2.江苏省厌氧生物技术重点实验室(江南大学),江苏无锡214122;3.江苏省高校水处理技术与材料协同创新中心,江苏苏州215009 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 卷积神经网络多变量过程特征学习与故障诊断
    卷积神经网络多变量过程特征学习与故障诊断陈淑梅,余建波(同济大学机械与能源工程学院,上海201804)摘要:为提取复杂多变量过程的有效特征,提高故障诊断性能,提出一种基于卷积神经网络(convolutionalneuralnetwork,CNN)特征学习的多变量过程故障诊断模型.将高维过程信号归一化 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 典型台风登陆过程平均风时距转换系数分析
    典型台风登陆过程平均风时距转换系数分析潘晶晶1,赵林1,梁旭东2,朱乐东1,葛耀君1(1.土木工程防灾国家重点实验室(同济大学),上海200092;2.中国气象局北京城市气象研究所,北京100089)摘要:为掌握中国东南沿海台风登陆过程近地面风速变化规律及影响因素,根据2005年浙江省东海塘观测塔和 ...
    本站小编 哈尔滨工业大学 2019-10-24
  • 臭氧氧化过程中羟基自由基产率测定与分析
    臭氧氧化过程中羟基自由基产率测定与分析刘永泽,江进,马军,罗从伟,皇甫小留,郭忠凯(城市水资源与水环境国家重点实验室(哈尔滨工业大学),150090哈尔滨)摘要:为研究臭氧氧化过程中OH·的产生,利用“叔丁醇方法”测定了芳香类、胺类、烯烃类等活性物质与O3反应过程中OH·的产率,并进一步测定了不同水 ...
    本站小编 哈尔滨工业大学 2019-10-24
  • 控制受限的卫星轨道转移过程中的姿态控制
    控制受限的卫星轨道转移过程中的姿态控制胡庆雷,张爱华,姜成平,晋小伟(哈尔滨工业大学航天学院,150001哈尔滨)摘要:针对刚体卫星在轨道转移过程中存在由质心漂移和推力偏心产生干扰力矩的控制问题,考虑存在参数不确定性、外部干扰以及控制输入受限,提出了一种智能自适应滑模控制方法,这种方法继承了传统 ...
    本站小编 哈尔滨工业大学 2019-10-24
  • 水下航行体垂直发射过程气幕形态特性研究
    水下航行体垂直发射过程气幕形态特性研究刘海军1,2,王聪1,彭兴芝2,邹振祝1,王本利1,隗喜斌1(1.哈尔滨工业大学航天学院,150001哈尔滨;2.河北大学建筑工程学院,071002河北保定)摘要:采用VOF多相流混合模型和标准k-ε湍流模型对水下航行体气幕垂直发射过程中气幕形态特性进行 ...
    本站小编 哈尔滨工业大学 2019-10-24
  • 波浪模拟及其对水下航行体出水过程影响
    波浪模拟及其对水下航行体出水过程影响权晓波,孔德才,李岩北京宇航系统工程研究所摘要:为了研究波浪对水下航行体出水过程的影响,将二阶Stokes波的波面运动速度耦合到入射边界条件中,得到了与理论值吻合良好、精度较高的非定常数值波浪;研究了五级海情条件下,水下航行体模型水下垂直发射时波浪对其流场、力学特 ...
    本站小编 哈尔滨工业大学 2019-10-24