CN112657435B - Co2加氢合成甲醇膜反应器及总熵生产率最小优化方法 - Google Patents
Co2加氢合成甲醇膜反应器及总熵生产率最小优化方法 Download PDFInfo
- Publication number
- CN112657435B CN112657435B CN202011553614.4A CN202011553614A CN112657435B CN 112657435 B CN112657435 B CN 112657435B CN 202011553614 A CN202011553614 A CN 202011553614A CN 112657435 B CN112657435 B CN 112657435B
- Authority
- CN
- China
- Prior art keywords
- reaction
- tube
- reaction tube
- mixture
- methanol
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P20/00—Technologies relating to chemical industry
- Y02P20/50—Improvements relating to the production of bulk chemicals
- Y02P20/52—Improvements relating to the production of bulk chemicals using catalysts, e.g. selective catalysts
Landscapes
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
本发明提供了一种CO2加氢合成甲醇膜反应器及总熵生产率最小优化方法,该反应器包括:反应管;渗透管,其同轴设置在反应管内,渗透管的材料为羟基方钠石渗透膜,反应管与渗透管之间填充有催化剂。本发明的反应器,通过在反应管内设置渗透管,渗透管由羟基方钠石渗透膜材料制成,反应中生成的水蒸气进入渗透管,渗透管内的水蒸气由吹扫气带出渗透管,如此使得反应不断朝正方向移动,进而提高CH3OH的产率;本发明的最小优化方法,以固定的甲醇收率为约束,以总熵产生率最小为优化目标,得到最优的反应管管外热源温度,与参考反应器相比,在甲醇产率固定的条件下,该Case 1最优反应器的总熵产生率降低了7.96%。
Description
技术领域
本发明涉及CO2加氢合成甲醇技术领域,尤其涉及一种CO2加氢合成甲醇膜反应器及总熵生产率最小优化方法。
背景技术
甲醇广泛应用于化工行业,可用于合成甲醛和醋酸等化学衍生物。甲醇还广泛应用于能源领域,它的辛烷值高,燃烧效率高;它在常温下是液体,可以作为储氢的载体。此外,它还是生物柴油和甲醇燃料电池的关键组成部分。甲醇合成工艺包括天然气合成甲醇、煤制甲醇、页岩气合成甲醇、生物制合成甲醇和等。由于天然气清洁且储量丰富,因此天然气合成甲醇在化工生产中占据重要地位。
在天然气合成甲醇的工艺流程中,天然气可通过水蒸气重整转化为富含甲烷、氢气(H2)、一氧化碳、氮气和CO2的合成气,合成气在高压中温条件下合成甲醇。现有技术公开了对包含CO2回收单元的天然气合成甲醇工艺流程进行了炯分析,结果表明优化甲烷蒸汽重整单元和甲醇合成单元具有提升整体性能的潜力。提高CO2加氢合成甲醇(MethanolSynthesis via CO2Hydrogenation,MSCH) 反应器的性能有利于天然气合成甲醇工艺流程的性能。受热力学平衡限制, MSCH反应器的甲醇产率和CO2转化率依然较低。基于此,有必要对现有的反应器结构进行改进。
同时,基于经典热力学的MSCH反应器和膜反应器的研究不涉及化学反应速率的概念,而实际化学反应均是在“有限时间”内或“有限面积”上进行“有限速率”的反应。因此,迫切需要利用有限时间热力学对实际化学反应器进行研究。有限时间热力学是20世纪70年代由芝加哥学派创立的一个现代热力学分支(工程界也称为“熵产最小化”理论,通过将热力学、传热学、流体力学和化学反应动力学等传输科学的基础理论相结合,在“有限时间”或“有限尺寸”约束下,利用最优控制理论等方法求解各类能量转换装置在不同优化目标时的最优性能和最优构型。
在MSCH膜反应器的有限时间热力学研究方面,现有技术公开建立了 MSCH膜反应器的一维瞬态模型,以甲醇产率最大为优化目标,在最大反应温度的约束下,得到了反应进口温度、冷却气进口温度和膜厚度的最优值。
熵产生率表征了反应器内不可逆过程的不可逆程度。根据Gouy-Stodola理论可知,不可逆过程的内部炯损等于环境温度与过程熵产的乘积。现有技术还公开了以MSCH反应的熵产生率最小和总熵产生率最小为优化目标,利用非线性规划和最优控制方法,在甲醇收率固定的约束下,对MSCH反应器进行最优构型优化设计。
然而,目前并没有公开对MSCH膜反应器的总熵产生率最小的最优构型优化设计。
发明内容
有鉴于此,本发明提出了一种CO2加氢合成甲醇膜反应器及总熵生产率最小优化方法,以解决现有技术中的甲醇产率和CO2转化率较低的技术问题。
第一方面,本发明提供了一种CO2加氢合成甲醇膜反应器,包括:
反应管;
渗透管,其同轴设置在反应管内,所述渗透管的材料为羟基方钠石渗透膜,所述反应管与所述渗透管之间填充有催化剂。
可选的,所述的CO2加氢合成甲醇膜反应器,所述催化剂为CuO/ZnO/Al2O3催化剂。
第二方面,本发明还提供了一种所述的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,其特征在于,包括以下步骤:
向反应管中通入合成气,合成气在催化剂的作用下反应生成甲醇,反应中生成的水蒸气进入渗透管,渗透管内的水蒸气由吹扫气带出渗透管;
计算反应管管外热源与反应管内混合物之间的传热过程引起的局部熵产生率;
计算反应管内混合物和渗透管内混合物之间的传热过程引起的局部熵产率;
计算反应管内反应过程引起的局部熵产生率;
计算反应管内混合物粘性流动过程引起的局部熵产生率;
计算反应管内混合物和渗透管内混合物之间传质过程引起的局部熵产生率;
将上述各步骤中的局部熵产生率求和得到局部总熵产生率,局部总熵产生率在反应器轴向积分可得到总熵产生率;
以固定的甲醇收率为约束,以总熵产生率最小为优化目标,得到最优的反应管管外热源温度。
可选的,所述的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,反应管管外热源与反应管内混合物之间传热过程引起的局部熵产生率的计算方法采用以下公式:
σHT,r=πdriJq,r(1/TRM-1/Ta)
其中,σHT,r为反应管管外热源与反应管内混合物之间传热过程引起的局部熵产生率;dri表示反应管内径;Ta表示反应管管外热源温度;TRM表示反应管内混合物温度;Jq,r表示通过反应管内壁的热流密度,其计算公式为: Jq,r=Ur(Ta-TRM),Ur表示为反应管管外热源和反应管内混合物之间的换热系数。
可选的,所述的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,反应管内混合物和渗透管内混合物之间的传热过程引起的局部熵产生率的计算方法采用以下公式:
σHT,m=πdmoJq,m(1/TMM-1/TRM)
其中,σHT,m表示反应管内混合物和渗透管内混合物之间的传热过程引起的局部熵产生率;dmo表示渗透管外径;TMM表示渗透管内混合物温度;Jq,m表示通过渗透管的热流密度,其计算公式为:Jq,m=Um(TRM-TMM),Um表示反应管内混合物与渗透管混合物之间的换热系数。
可选的,所述的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,计算反应管内反应引起的局部熵产生率具体包括:
其中,i=1,2,为反应管内反应引起的局部熵产生率,ρc为催化剂的颗粒密度,Ac,r为反应管的截面积,εB为催化剂床层孔隙率,ηi为内扩散有效因子, r1为甲醇合成反应的反应速率,r2为逆水气变换反应的反应速率,-ΔrG1为甲醇合成反应的吉布斯自由能,-ΔrG2为逆水气变换反应的吉布斯自由能。
可选的,所述的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,反应管内混合物粘性流动引起的熵产生率的计算方法采用以下公式:
其中,cRM表示反应管内合成气的表面流速;
Rep为反应管内合成气在催化剂床层流动的雷诺数,PRM为反应管内混合物的压力,GRM为反应管内混合物的质量流率。
可选的,所述的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,反应管内混合物和渗透管内混合物之间传质过程引起的局部熵产生率的计算方法采用以下公式:
可选的,所述的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,所述吹扫气为氮气。
可选的,所述的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,以固定的甲醇收率为约束,以总熵产生率最小为优化目标,得到最优的反应管管外热源温度之后还包括:以得到的最优的反应管管外热源温度为基础,进一步以总熵产生率最小为优化目标,得到最优的合成气进入反应管的温度。
本发明的一种CO2加氢合成甲醇膜反应器及总熵生产率最小优化方法相对于现有技术具有以下有益效果:
(1)本发明的CO2加氢合成甲醇膜反应器,通过在反应管内设置渗透管,渗透管由羟基方钠石渗透膜材料制成,反应中生成的水蒸气进入渗透管,渗透管内的水蒸气由吹扫气带出渗透管,如此使得反应不断朝正方向移动,进而提高 CH3OH的产率;
(2)本发明的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,以固定的甲醇收率为约束,以总熵产生率最小为优化目标,得到最优的反应管管外热源温度,与参考反应器相比,在甲醇产率固定的条件下,该Case 1最优反应器的总熵产生率降低了7.96%;
(3)本发明的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,进一步在Case1反应器的基础上,以总熵产生率最小为优化目标,对反应管管外热源温度和合成气进入反应管的温度同时进行优化,与参考反应器相比,在甲醇产率固定的条件下,该Case2最优反应器的总熵产生率降低了10.89%。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的CO2加氢合成甲醇膜反应器的结构示意图;
图2为本发明的膜反应器的反应时通入吹扫气的示意图;
图3为本发明的反应温度和反应管进口CO2摩尔分数对甲醇收率的影响的仿真结果与实验结果对比图;
图4为本发明的参考反应器和Case 1最优反应器的管外热源温度、反应混合物温度和渗透管内混合物温度的沿程分布曲线图;
图5为本发明的参考反应器和Case 1最优反应器的传热过程引起的局部熵产生率(σHT,ref和σHT,Case 1)的沿程分布曲线图;
图7为本发明的Case 2最优反应器的管外热源温度、反应混合物温度和渗透管内侧混合物温度的沿程分布曲线图;
图9为本发明的Case 1最优反应器和Case 2最优反应器的传热过程与2个子传热过程引起的局部熵产生率(σHT,Case1、σHT1,Case1、σHT2,Case1、σHT,Case2、σHT1,Case2和σHT2,Case2)的沿程分布曲线图。
具体实施方式
下面将结合本发明实施方式,对本发明实施方式中的技术方案进行清楚、完整的描述,显然,所描述的实施方式仅仅是本发明一部分实施方式,而不是全部的实施方式。基于本发明中的实施方式,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式,都属于本发明保护的范围。
如图1所示,本发明提供了一种CO2加氢合成甲醇膜反应器,包括:
反应管1;
渗透管2,其同轴设置在反应管1内,渗透管2的材料为羟基方钠石渗透膜,反应管1与渗透管2之间填充有催化剂3。
需要说明的是,本申请实施例中,反应管1、渗透管2均为空心的圆柱体状,并且渗透管2同轴套设在反应管1内,反应管1与渗透管2之间填充有催化剂3,催化剂3可催化CO2与H2反应生成CH3OH,渗透管2由羟基方钠石渗透膜材料制成,水蒸气可透过渗透管2进入渗透管2内,反应管1一端为进口、另一端为出口。该羟基方钠石渗透膜(H-SOD膜)为现有渗透膜,本申请并没有对该羟基方钠石渗透膜进行改进,具体的,该羟基方钠石渗透膜(H-SOD膜)对水的渗透性在1×10-7~3×10-7,对水和氢气、水和一氧化碳、水和CO2的选择性大于75。
在一些实施例中,催化剂为CuO/ZnO/Al2O3催化剂,具体的该催化剂采用现有的催化剂。
本申请的CO2加氢合成甲醇膜反应器,如图2所示,在使用时,首先在反应管1的一端(进口端)且位于反应管1与渗透管2之间的空间通入合成气(包括CO2、H2、CO、CH4的混合气体)或直接通入CO2与H2,CO2与H2在催化剂的催化作用下在反应管1和渗透管2之间生成CH3OH,同时反应中生成的水蒸气经渗透管2渗透进入渗透管2内,同时在渗透管2的一端通入吹扫气,在吹扫气的作用下,带走水蒸气,如此使得反应不断朝正方向移动,进而提高CH3OH的产率。
基于同一发明构思,本发明还提供了一种CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,包括以下步骤:
S1、向反应管中通入合成气,合成气在催化剂的作用下反应生成甲醇,反应中生成的水蒸气进入渗透管,渗透管内的水蒸气由吹扫气带出反应管;
S2、计算反应管管外热源与反应管内混合物之间的传热过程引起的局部熵产生率;
S3、计算反应管内混合物和渗透管内混合物之间的传热过程引起的局部熵产生率;
S4、计算反应管内反应过程引起的局部熵产生率;
S5、计算反应管内混合物粘性流动过程引起的局部熵产生率;
S6、计算反应管内混合物和渗透管内混合物之间传质过程引起的局部熵产生率;
S7、将上述各步骤中的局部熵产生率求和得到局部总熵产生率,局部总熵产生率在反应器轴向积分可得到总熵产生率;
S8、以固定的甲醇收率为约束,以总熵产生率最小为优化目标,得到最优的反应管管外热源温度。
需要说明的是,本申请实施例中管外热源温度指的是反应器外对反应管进行加热的热源温度。
具体的,本申请实施例中,反应管管外热源与反应管内混合物之间传热过程引起的局部熵产生率σHT,r的计算方法采用以下公式:
σHT,r=πdriJq,r(1/TRM-1/Ta)
其中,dri表示反应管内径;Ta表示反应管管外热源温度;TRM表示反应管内混合物温度;Jq,r表示通过反应管内壁的热流密度,其计算公式为: Jq,r=Ur(Ta-TRM),Ur表示为反应管管外热源和反应管内混合物之间的换热系数。
具体的,本申请实施例中,反应管内混合物和渗透管内混合物之间的传热过程引起的局部熵产生率σHT,m的计算方法采用以下公式:
σHT,m=πdmoJq,m(1/TMM-1/TRM);
其中,dmo表示渗透管外径;TMM表示渗透管内混合物温度,具体的,本申请中渗透管内的混合物包括水蒸气和吹扫气的混合物;
Jq,m表示通过渗透管的热流密度,其计算公式为:
Jq,m=Um(TRM-TMM),
Um表示反应管内混合物与渗透管混合物之间的换热系数。
具体的,本申请实施例中,计算反应管内反应引起的局部熵产生率σri具体包括:
反应管内反应引起的局部熵产生率的计算公式为
其中,i=1,2,
ρc为催化剂的颗粒密度,Ac,r为反应管的截面积,εB为催化剂床层孔隙率,ηi为内扩散有效因子,r1为甲醇合成反应的反应速率,r2为逆水气变换反应的反应速率,ΔrG1为甲醇合成反应的吉布斯自由能,ΔrG2为逆水气变换反应的吉布斯自由能。
具体的,本申请实施例中,甲醇合成反应为放热的MSCH反应,逆水气变换反应为吸热反应,具体的,
κj=Ajexp(Bj/RgTRM),j=1,2,3,4,5;
Aj、Bj均为常数;
Rg表示热力学常数。
具体的,本申请实施例中,反应管内混合物粘性流动引起的熵产生率σVF的计算方法采用以下公式:
其中,cRM表示反应管内混合物的表面流速;
Rep为反应管内混合物在催化剂床层流动的雷诺数,PRM为反应管内混合物的压力,GRM为反应管内混合物的质量流率。
具体的,本申请实施例中,反应管内混合物和渗透管内混合物之间传质过程引起的局部熵产生率σMT的计算方法采用以下公式:
具体的,本申请实施例中,吹扫气为氮气。使用氮气不仅可以吹带出渗透管内的水蒸气,同时氮气不参与反应。
具体的,本申请实施例中,将上述各步骤中的局部熵产生率求和得到局部总熵产生率,局部总熵产生率在反应器轴向积分可得总熵产生率,以固定的甲醇收率为约束,以总熵产生率最小为优化目标,得到最优的反应管管外热源温度。
即总熵产生率为:
SgHT、Sgr、SgVF和SgMT分别为传热、化学反应、粘性流动和传质过程的熵产生率,σT为各局部总熵产生率产生率之和,其中,L为反应管的长度;Z 表示反应器轴向坐标(以反应器进口为原点)。
具体的,本申请中反应管中的能量守恒方程为:
其中,Fk和Cp,k分别为组分k的摩尔流率和定压摩尔热容,ΔrHi为反应i的反应焓变。具体的,k中1~7分别代表CO2、H2、CO、CH3OH、H2O、N2、CH4七种组分。
甲醇的收率ξ1和一氧化碳的收率ξ2由物量守恒方程描述:
式中,vk,i为反应i中组分k的化学计量系数。具体的,i=1对应甲醇合成反应,i=2对应逆水气变换反应,例如代表反应1即甲醇合成反应中CO2的化学计量系数,代表反应2即逆水气变换反应中CO2的化学计量系数。
水蒸气的摩尔流率为:
渗透管内的物料守恒方程为:
渗透管内的能量守恒方程为:
以下以具体实例进一步说明本申请的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法。
具体的,本申请中反应器做以下限定:反应混合物和吹扫气为理想气体;管外热源温度完全可控;膜反应器为一维拟均相活塞流反应器;反应器在径向的扩散和轴向的返混可以忽略;催化剂填充均匀,催化剂床层孔隙率为常数;羟基方钠石渗透膜对水蒸气的选择透过性为100%;内扩散有效因子ηi为1;忽略渗透管内的压降。
本申请以MSCH膜反应器的总熵产生率最小为优化目标,以完全可控的反应管管外热源温度(Ta)、反应混合物进口温度和吹扫气进口温度为优化变量。以固定的甲醇收率(ΔCH3OHfixed)为约束。本申请优化问题的目标和约束如下:
s.t.
本申请的优化问题中包含了7个微分方程组和若干辅助方程,方程数量多且相互耦合,因此无法直接得到优化问题的解析解,可以采用非线性规划进行数值求解。在非线性规划求解过程中,运用有限差分法(即微分方程的离散化处理),将反应器离散为n=1,2,...N个控制体。通过优化每一个控制体的管外热源温度(Ta(n))、和使得总熵产生率达到最小,其中,表示反应混合物的进口温度,表示渗透管内混合物的进口温度。因此,优化变量共有N+2 个。
该非线性规划问题可以利用MATLAB优化工具箱里的“fmincon”函数进行优化,在优化过程中选择了序列二次规划(SQP)算法。在优化前需要设定优化变量的变化范围,本文优化变量的变化范围为:
400≤Ta(n)≤600
优化目标可利用梯形积分进行表述如下:为了提高总熵产生率的计算精度,利用梯形积分对其进行计算,计算公式如下:
式中,Δz=L/N为步长。
约束为:
上述的优化问题旨在满足守恒方程和约束方程的条件下求解管外热源温度的最优轨迹使得总熵产生率最小,属于典型的最优控制问题。基于最优控制理论,可建立哈密顿函数:
上式中,哈密顿函数由两部分组成,右侧第一项为局部总熵产生率,右侧第二项为协态变量与守恒方程约束的乘积。哈密函数不显含坐标z,根据哈密函数的性质可判断H[x(z),λ(z),Ta(z)]为常数,该性质可用来检验优化结果的准确性。
根据Pontryagins极小值原理,最优控制问题的解需要满足状态方程:
协态方程:
极小值条件:
可得管外热源温度(Ta(z)):
针对本申请优化问题的非线性规划和最优控制的数值优化方法。其中,非线性规划方法最多包含了N+2个优化变量,无法确保得到的解是全局最优解,而最优控制利用数学工具推导出了管外热源温度的最优路径,将寻优问题转化为包含边界条件的微分方程组求解问题,能够确保得到的解是全局最优解,并且降低了运算时间。本申请的最优控制问题包含了14个一阶微分方程,属于典型的两点边界值问题,现阶段可以利用MATLAB的bvp4c函数对其进行求解,然而,bvp4c函数对初始值较为依赖。因此,需要设置合理的初始值使得微分方程组的解收敛。本文利用fmincon计算得到的非线性规划的优化结果作为最优控制的初始值。
MSCH膜反应器的操作参数、催化剂的几何参数、渗透膜的几何参数和物理参数等参考工业反应器的数据取值。本申请的MSCH膜反应器的几何参数由反应空速、各组分进口摩尔分数和长径比计算得到。表1列出了本申请的MSCH 膜反应器的参数设置。
表1-MSCH膜反应器的参数设置
具体的,表1中表示反应管进口CO2进口摩尔分数,同理其他的均表示不同组分的进口处的摩尔分数,表示反应管进口的氢气摩尔流率,dmi表示渗透管的内径,GHSV表示反应混合物空速,δm表示羟基方钠石渗透膜的厚度,表示反应管进口混合物温度,表示反应管进口混合物压力,表示渗透管进口处混合物温度,表示渗透管进口处混合物的压力,表示渗透管进口处水蒸气的摩尔流率。
理想的反应动力学模型和合理的参数设置是MSCH膜反应器优化的必要条件。为了验证MSCH反应动力学和热力学参数的准确性,本申请建立了等温等压反应器模型,并利用Klier等(Klier K,Chatikavanij V,Herman R G,Simmons G W.Catalytic synthesis ofmethanol from CO/H2 IV.The effects of carbon dioxide.J. Catal.,1982,74(2):343-360)的实验验证了模型的准确性。图3给出了反应温度和反应管进口CO2摩尔分数对甲醇收率的影响的仿真结果与实验结果对比。由图3可得,随着反应管进口CO2摩尔分数的增加,反应管出口的甲醇收率先增大后减小。这是由于,当增加时,MSCH反应速率增大,因此增大;当大于6%时,随着的增加,反应管进口一氧化碳摩尔分数降低,水汽变换反应减弱,氢气摩尔分数降低,因此减小。当反应管进口反应混合物温度等于523K,小于10%时,仿真结果与实验结果十分吻合。当等于508K时,仿真结果均与实验结果吻合。
为了得到MSCH膜反应器的最小总熵产生率与相应的管外热源温度最优构型、反应管进口反应混合物温度和反应管进口吹扫气温度,分析优化变量和其他操作参数与几何参数对MSCH膜反应器的总熵产生率的影响程度,本申请的优化步骤如下:
(1)Case 1:以管外热源温度(Ta)曲线为优化变量,优化得到Case 1最优反应器的最小总熵产生率(SgT,Case1)和管外热源温度(Ta,Case1)的最优构型。
(2)Case 2:在步骤(1)的基础上,研究SgT,Case1和Ta,Case1随反应管进口反应混合物温度的变化规律,得到Case 2最优反应器的最小总熵产生率(SgT,Case2),管外热源温度(Ta,Case2)的最优构型和反应管进口反应混合物温度的最优值。
根据优化步骤,表2列出了Case 1,Case 2最优反应器的状态变量和协态变量的边界值。
表2-Case 1,Case 2最优反应器的状态变量和协调变量的边界值
表2中上标“0”均表示参数位于反应管进口的状态,上表“1”表示参数位于反应管出口的状态。
表3列出了参考反应器、Case 1最优反应器、Case 2最优反应器的熵产生率、甲醇产率、CO2转化率和相应的优化变量取值。
表3-参考反应器、Case 1、Case 2最优反应器性能对比
由表3可知,与参考反应器相比,在甲醇产率固定的条件下,Case 1、Case 2最优反应器的总熵产生率分别降低了7.96%,10.89%。其中,粘性流动过程和传质过程的总熵产率变化不大。由于甲醇产率固定,CO2的转化率也基本不变。
图4给出了参考反应器和Case 1最优反应器的反应管管外热源温度、反应混合物温度和渗透管内混合物温度的沿程分布曲线。由图4可见,参考反应器和Case 1最优反应器的反应混合物温度分布曲线(TRM,ref和TRM,Case1)的变化趋势基本一致,均是先增大后减小,TRM,ref在反应器的末端略有升高。这是由于,MSCH 反应是放热反应,在反应器进口附近,反应物料充足,反应速率快,因此反应混合物温度先升高。随着氢气摩尔分数的降低,吸热的水汽变换反应的反应速率提高,因此反应混合物的温度逐渐降低。由于出口附近参考反应器的管外热源温度(Ta,ref)大于TRM,ref,因此TRM,ref在出口处略微升高。Case 1最优反应器的管外热源温度(Ta,Case1)先升高,再降低最后升高,出口处的Ta,Case1与TRM,Case1一致。这是由于,当Ta,Case1与TRM,Case1的变化趋势保持一致时,管外热源和反应混合物之间的温差减小,传热的驱动力减小,传热引起的熵产生率也减小。由管外热源温度(Ta(z))公式和表2中的取值可知,出口处的Ta,Case1等于TRM,Case1,此时能够保证出口处传热过程引起的局部熵产生率为0。参考反应器和Case 1最优反应器的渗透管内混合物温度(TMM,ref和TMM,Case1)均呈线性增大,且变化幅度不大,这是由于渗透管内混合物的流动模式为层流流动,导致渗透管内混合物和反应混合物之间的传热系数较小,进而导致渗透管内混合物的吸热率较小。
图5给出了参考反应器和Case 1最优反应器的传热过程引起的局部熵产生率(σHT,ref和σHT,Case1)的沿程分布曲线。由图5可见,σHT,ref沿反应器无量纲轴向坐标先增大后减小,最后再增大;σHT,Case1沿先减小后增大,最后再减小,在等于0.3到0.8之间,σHT,Case1的变化不明显。σHT,ref和σHT,Case1的沿程变化趋势分别与图4中的Ta,ref与TRM,ref和Ta,Case1与TRM,Case1的温差变化趋势一致。σHT,ref和σHT,Case1与σHT=0围成的面积分别代表参考反应器和Case 1最优反应器的传热过程引起的总熵产生率(SgHT,ref和SgHT,Case1)。由表3可知,与SgHT,ref相比,SgHT,Case1降低24.29%,与参考反应器的总熵产生率(SgT,ref)相比,Case1最优反应器的总熵产生率(SgT,Case1)降低了7.96%。
图6给出了SgT,Case1、SgHT,Case1和Case 1最优反应器的化学反应过程引起的熵产生率(Sgr,Case1)与的关系。如图6所示,随着的增大,SgT,Case1先减小后增大。当等于518 K时,得到了Case 2最优反应器的最小总熵产生率 (SgT,Case2=0.741W·K-1),传热和化学反应过程的熵产生率(SgHT,Case2=0.203W·K-1和Sgr,Case2=0.461 W·K-1)。与SgT,Case1和Sgr,Case1相比,SgT,Case2和Sgr,Case2分别降低了3.14%和22.78%;与SgHT,Case1相比,SgHT,Case2增加了63.71%。因此,在传热过程和化学反应过程引起的熵产生率的综合影响下,SgT,Case1随的增加先减小再增大。由图6可见,随着的增加,Sgr,Case1先减小后增大,当时,Sgr,Case1最小,为0.291。由图6可见,随着的增加,SgHT,Case1逐渐增大。
图7给出了Case 2最优反应器的管外热源温度、反应混合物温度和渗透管内侧混合物温度(Ta,Case2、TRM,Case2和TMM,Case2)的沿程分布曲线。由图4和图7可见, Case 2最优反应器的反应混合物温度(TRM,Case2)、管外热源温度(Tw,Case2)和渗透管内混合物温度(TMM,Case2)与TRM,Case1、Tw,Case1和TMM,Case1的变化趋势基本一致。Tw,Case2和TRM,Case2在端点处重合,这是由于Case 2最优反应器释放了的约束,由表1 可知,等于0,因此等于
图8给出了Case 2最优反应器和分别取565和580K时的Case 1最优反应器的化学反应过程引起的局部熵产生率(σr,Case2,和)的沿程分布曲线。由图可见,σr,Case2和σr,Case1均是先急剧减小,再保持稳定,最后再减小到0。在反应器进口附近(图8中的‘I’区),σr,Case1均大于σr,Case2,当等于0.01-0.2时(图8中的‘II’区)σr,Case1均小于σr,Case2,当大于 0.2后(图8中‘III’区),σr,Case1均大于σr,Case2。随着的增加,在‘II’区内化学反应产生的熵产生率降低,在‘I’区和‘III’区内化学反应产生的熵产生率增大,当大于565K以后,‘I’区和‘III’区内化学反应产生的熵产生率增大的速度明显超过‘II’区内化学反应产生的熵产生率降低的速度,因此Sgr,Case1随着的增加,先减小后增加。
图9给出了Case 1最优反应器和Case 2最优反应器的传热过程与2个子传热过程引起的局部熵产生率(σHT,Case1、σHT1,Case1、σHT2,Case1、σHT,Case2、σHT1,Case2和σHT2,Case2)的沿程分布曲线。如图所示,Case 1和Case 2最优反应器中渗透管内的混合物与反应混合物之间的传热过程引起的熵产生率(σHT2,Case1和σHT2,Case2)较小,这是由于渗透管内混合物和反应混合物之间的换热系数(Um)较小,导致渗透管内混合物与反应混合物之间的换热过程的热流率较小,进而导致σHT2,Case1和σHT2,Case2较小。σHT2,Case2略大于σHT2,Case1,这是由于与TMM,Case2之间的差值大于 TRM,Case1与TMM,Case1之间的差值。Case 2最优反应器中反应混合物与管外热源之间的传热过程引起的熵产生率(σHT1,Case2)与TRM,Case2和Ta,Case2之差的趋势一致。与Case 1最优反应器中反应混合物与管外热源之间的传热过程引起的熵产生率(σHT1,Case1) 相比,σHT1,Case2更大。这是由于TRM,Case2和Tw,Case2之差更大。因此,随着的增大,TRM,Case1与Tw,Case1和TMM,Case1之间的温差增大,进而导致SgHT,Case1增大。
以上所述仅为本发明的较佳实施方式而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,其特征在于,所述CO2加氢合成甲醇膜反应器,包括:
反应管;
渗透管,其同轴设置在反应管内,所述渗透管的材料为羟基方钠石渗透膜,所述反应管与所述渗透管之间填充有催化剂;
所述催化剂为CuO/ZnO/Al2O3催化剂;
所述CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法包括以下步骤:
向反应管中通入合成气,合成气在催化剂的作用下反应生成甲醇,反应中生成的水蒸气进入渗透管,渗透管内的水蒸气由吹扫气带出渗透管;计算反应管管外热源与反应管内混合物之间的传热过程引起的局部熵产生率;
计算反应管内混合物和渗透管内混合物之间的传热过程引起的局部熵产生率;
计算反应管内反应过程引起的局部熵产生率;
计算反应管内混合物粘性流动过程引起的局部熵产生率;
计算反应管内混合物和渗透管内混合物之间传质过程引起的局部熵产生率;
将上述各步骤中的局部熵产生率求和得到局部总熵产生率,局部总熵产生率在反应器轴向积分可得到总熵产生率;
以固定的甲醇收率为约束,以总熵产生率最小为优化目标,得到最优的反应管管外热源温度。
2.如权利要求1所述的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,其特征在于,反应管管外热源与反应管内混合物之间传热过程引起的局部熵产生率的计算方法采用以下公式:
σHT,r=πdriJq,r(1/TRM-1/Ta)
其中,σHT,r为反应管管外热源与反应管内混合物之间传热过程引起的局部熵产生率;dri表示反应管内径;Ta表示反应管管外热源温度;TRM表示反应管内混合物温度;Jq,r表示通过反应管内壁的热流密度,其计算公式为:Jq,r=Ur(Ta-TRM),Ur表示为反应管管外热源和反应管内混合物之间的换热系数。
3.如权利要求2所述的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,其特征在于,反应管内混合物和渗透管内混合物之间的传热过程引起的局部熵产生率的计算方法采用以下公式:
σHT,m=πdmoJq,m(1/TMM-1/TRM)
其中,σHT,m表示反应管内混合物和渗透管内混合物之间的传热过程引起的局部熵产生率;dmo表示渗透管外径;TMM表示渗透管内混合物温度;Jq,m表示通过渗透管的热流密度,其计算公式为:Jq,m=Um(TRM-TMM),Um表示反应管内混合物与渗透管混合物之间的换热系数。
7.如权利要求1所述的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,其特征在于,所述吹扫气为氮气。
8.如权利要求1所述的CO2加氢合成甲醇膜反应器的总熵生产率最小优化方法,其特征在于,以固定的甲醇收率为约束,以总熵产生率最小为优化目标,得到最优的反应管管外热源温度之后还包括:以得到的最优的反应管管外热源温度为基础,进一步以总熵产生率最小为优化目标,得到最优的合成气进入反应管的温度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011553614.4A CN112657435B (zh) | 2020-12-24 | 2020-12-24 | Co2加氢合成甲醇膜反应器及总熵生产率最小优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011553614.4A CN112657435B (zh) | 2020-12-24 | 2020-12-24 | Co2加氢合成甲醇膜反应器及总熵生产率最小优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112657435A CN112657435A (zh) | 2021-04-16 |
CN112657435B true CN112657435B (zh) | 2022-05-31 |
Family
ID=75408474
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011553614.4A Active CN112657435B (zh) | 2020-12-24 | 2020-12-24 | Co2加氢合成甲醇膜反应器及总熵生产率最小优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112657435B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115943194B (zh) * | 2021-06-08 | 2024-06-14 | 日本碍子株式会社 | 液体燃料合成系统 |
JPWO2023120494A1 (zh) * | 2021-12-22 | 2023-06-29 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1147241A (zh) * | 1994-04-08 | 1997-04-09 | 卡萨尔甲醇公司 | 用于气体原料转化的膜反应器 |
US20030066240A1 (en) * | 2001-10-09 | 2003-04-10 | Conoco Inc. | Steam reformer for methane with internal hydrogen separation and combustion |
US20130090394A1 (en) * | 2011-10-07 | 2013-04-11 | Chevron U.S.A Inc. | Process for providing hydrogen to a synthesis gas conversion reaction |
CN103939941A (zh) * | 2014-04-25 | 2014-07-23 | 国家电网公司 | 一种融入了不可逆热力学的锅炉燃烧优化方法 |
CN104298864A (zh) * | 2014-09-30 | 2015-01-21 | 杭州电子科技大学 | 一种合成气一步法制二甲醚动态建模方法 |
CN111196596A (zh) * | 2020-01-14 | 2020-05-26 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 一种流速和浓度均匀分布的微通道甲醇制氢反应器 |
CN111867718A (zh) * | 2018-01-26 | 2020-10-30 | 荷兰应用自然科学研究组织Tno | 用于二氧化碳加氢反应的反应器和工艺 |
-
2020
- 2020-12-24 CN CN202011553614.4A patent/CN112657435B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1147241A (zh) * | 1994-04-08 | 1997-04-09 | 卡萨尔甲醇公司 | 用于气体原料转化的膜反应器 |
US20030066240A1 (en) * | 2001-10-09 | 2003-04-10 | Conoco Inc. | Steam reformer for methane with internal hydrogen separation and combustion |
US20130090394A1 (en) * | 2011-10-07 | 2013-04-11 | Chevron U.S.A Inc. | Process for providing hydrogen to a synthesis gas conversion reaction |
CN103939941A (zh) * | 2014-04-25 | 2014-07-23 | 国家电网公司 | 一种融入了不可逆热力学的锅炉燃烧优化方法 |
CN104298864A (zh) * | 2014-09-30 | 2015-01-21 | 杭州电子科技大学 | 一种合成气一步法制二甲醚动态建模方法 |
CN111867718A (zh) * | 2018-01-26 | 2020-10-30 | 荷兰应用自然科学研究组织Tno | 用于二氧化碳加氢反应的反应器和工艺 |
CN111196596A (zh) * | 2020-01-14 | 2020-05-26 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 一种流速和浓度均匀分布的微通道甲醇制氢反应器 |
Non-Patent Citations (5)
Title |
---|
Dulong-Petit传热规律下甲烷蒸汽重整反应器;陈庆晓等;《节能》;20180331(第3期);第31-40页 * |
Entropy generation minimization for CO2 hydrogenation to light olefins;Lingen Chen etal;《Energy》;20180315;第147卷;第187-196页 * |
Entropy generation rate minimization for hydrocarbon synthesis reactor from carbon dioxide and hydrogen;Lei Zhang etal;《International Journal of Heat and Mass Transfer》;20190731;第137卷;第1112-1123页 * |
Entropy Generation Rate Minimization for Methanol Synthesis via a CO2 Hydrogenation Reactor;Penglei Li etal;《entropy》;20190228;第21卷(第2期);第1-20页 * |
基于熔融盐加热的甲烷蒸汽重整制氢反应器的熵产生率和氢气产率分析;李鹏蕾等;《发电技术》;20200831;第41卷(第4期);第415-428页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112657435A (zh) | 2021-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jin et al. | Experimental and simulation study on a catalyst packed tubular dense membrane reactor for partial oxidation of methane to syngas | |
Li et al. | Multi-objective optimal configurations of a membrane reactor for steam methane reforming | |
Basile et al. | The partial oxidation of methane to syngas in a palladium membrane reactor: simulation and experimental studies | |
Samimi et al. | Green methanol production process from indirect CO2 conversion: RWGS reactor versus RWGS membrane reactor | |
Samimi et al. | Green methanol synthesis process from carbon dioxide via reverse water gas shift reaction in a membrane reactor | |
CN112657435B (zh) | Co2加氢合成甲醇膜反应器及总熵生产率最小优化方法 | |
Alrashed et al. | Comparative analysis of conventional steam methane reforming and PdAu membrane reactor for the hydrogen production | |
De Falco et al. | Hydrogen production by solar steam methane reforming with molten salts as energy carriers: Experimental and modelling analysis | |
Piemonte et al. | Counter-current membrane reactor for WGS process: Membrane design | |
Saithong et al. | Thermodynamic analysis of the novel chemical looping process for two-grade hydrogen production with CO2 capture | |
Castillo et al. | Effect of temperature and pressure on hydrogen production from steam reforming of biogas with Pd–Ag membrane reactor | |
Jeong et al. | Effectiveness factor correlations from simulations of washcoat nickel catalyst layers for small-scale steam methane reforming applications | |
Kong et al. | Performance analysis of hydrogen iodide decomposition membrane reactor under different sweep modes | |
Li | First-principles-based microkinetic rate equation theory for oxygen carrier reduction in chemical looping | |
Tang et al. | An integrated modeling method for membrane reactors and optimization study of operating conditions | |
Talaghat et al. | Mathematical modeling of hydrogen production using methanol steam reforming in the coupled membrane reactor when the output materials of the reformer section is used as feed for the combustion section | |
Klein et al. | Dry methane reforming without a metal catalyst in a directly irradiated solar particle reactor | |
Xie et al. | Reaction/separation coupled equilibrium modeling of steam methane reforming in fluidized bed membrane reactors | |
Rahimpour et al. | Co‐current and countercurrent configurations for a membrane dual type methanol reactor | |
Iliuta et al. | Intensified bio-oil steam reforming for high-purity hydrogen production: Numerical simulation and sorption kinetics | |
Palone et al. | Synthesis of methanol from a chemical looping syngas for the decarbonization of the power sector | |
Chiu et al. | Experimental study on the reaction conditions of a methanol steam reforming process | |
Rathod et al. | Analytical and experimental investigations for hydrogen rich syngas production by biogas reforming processes | |
Chen et al. | Numerical modelling for design of microchannel reactors: Application to hydrogen production from methanol by steam reforming | |
Holt et al. | Numerical modelling of a steam methane reformer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |