高压储氢罐充放气过程的热效应模拟与性能预测

发布时间：2020-07-15 09:30

【摘要】：本文针对高压储氢罐充放气过程中的热效应展开机理研究,主要研究充放气过程中罐内氢气温度、氢气压力和罐壁温度等随着时间的变化规律,以及不同充气参数对热效应的影响。针对这一问题,本文运用质量守恒方程和能量守恒方程建立了高压储氢罐充放气过程热效应的理论模型,该模型将整个高压储氢罐系统分为两个区域,一个是气态区;另一个是罐壁区。针对罐壁是否进行简化处理,可分别建立相应的单区单温模型和双区双温模型,对模型进行推导可得到相应的单区解析解和双区解析解;另外利用Matlab/Simulink仿真平台可以搭建相应的集总参数模型,分别得到相应的单区数值解和双区数值解,同时Matlab/Simulink还可将罐壁中的内衬和外套分开考虑,从而得到相应的三区(气态区、内衬区和外套区)数值模型。数值解模拟的结果可以和解析解计算的结果进行对比,从而验证理论模型的有效性。基于单区模型的解析解,本文对预冷氢气温度、充气质量和充气状态(SOC)进行了相应的参数研究。为了避免氢气温度超过85℃,从而保证储氢罐的安全性能,需要对氢气进行预冷处理,预冷氢气温度的选择成为一个关键问题。为了解决这一问题,本文在单区模型的基础上,成功推导了预冷氢气温度的解析解,给出了预冷氢气温度和各充气参数(初始温度、初始压力和充气时间)之间的函数关系,并计算出了当最高氢气温度恰好达到358K时所对应的临界预冷条件,研究发现,当实际初始温度高于临界温度、实际初始压力小于临界压力,实际充气时间比临界时间短时,均需要对氢气进行预冷处理。本文还研究了充气质量随充气参数之间的变化规律,充气质量的解析解可以由氢气温度的解析解推导得到,充气质量可以分别表示为关于充气速率、初始氢气压力和进气温度的函数关系,利用充气质量的解析解对相应的数据进行拟合,可以确定充气质量与各充气参数之间的定量关系。基于储氢罐中氢气密度均匀分布的假设以及充气质量的解析解,本文还提出了一个计算SOC的数学模型,模型中相应的参数值可以用最小二乘法确定;在研究充气速率和进气温度同时对SOC的影响时,可以在换热系数的表达式中引入努塞尔数和雷诺数来对原有的模型进行修正,通过研究发现,新的模型具有更好的拟合效果以及表达数据的能力。
【学位授予单位】：武汉理工大学
【学位级别】：硕士
【学位授予年份】：2018
【分类号】：U469.7
【图文】：

第 2 章 高压储氢罐充/放气过程理论模型及解析解本章主要阐述理论模型的建立过程以及单区和双区解析解的推导过程。2.1 理论模型如图 2-1 所示，为高压储氢罐充放气过程的传热传质示意图。氢气从进气口处进入储氢罐，同时伴随着能量的变化，从而引起气体温度和压力的改变，整个储氢罐系统的热效应可以用相应的质量守恒方程和能量守恒方程进行表示。

建立了合适的数学物理模型，然而在对应微分方程求解的过的结果似乎无法精确得到，上文虽然通过解耦微分方程的办法，得和罐壁温度的近似解析解，但也对实际结果带来了一定的误差。为描述相应的物理规律，精确得到氢气温度和罐壁温度等随充气参数关系，这里还可利用一些数值方法对物理模型进行求解，通过b/Simulink 仿真平台，得到了相应的数值解。数值解的出现，为实际明了新的方向，同时，它还可跟近似解析解进行对比，用来验证近效性和适用性。单区模型数值解 Matlab/Simulink 仿真模拟平台中描述单区模型的质量守恒方程，如。其中充气速率和放气速率（如图 3-2）可用平台中内置的常数项过一定的数学运算模块复合之后，再利用积分模块可得氢气质量随系，积分运算中需要输入初始质量 m0的值。

图 3-2 放气速率，在 Matlab/Simulink 平台中描述能量守恒方程，如图 3-3 为待求项，显然是未知的，而在能量守恒方程中，又需要内置的 from 和 goto 模块表示。

【相似文献】

相关期刊论文 前10条

1 蔡睿贤;;热力机械基本方程的严格解析解[J];机械工程学报;2009年03期

2 汪懋骅;方程u_(tt)-cu_(xx)+au+bu~3=0的解析解答[J];应用数学和力学;1988年07期

3 冯登泰;磨耗形轮轨弹塑性静力接触的解析解[J];铁道学报;1988年02期

4 王文信,田石,冯庆典;Ⅰ型裂纹HRR模型的解析解[J];应用力学学报;1989年04期

5 车维毅;;框架结构的动力微分方程和静力解析解[J];辽宁大学学报(自然科学版);1989年04期

6 吴开谡 ,陈永寿 ,舒能川;重元素s-过程解析解[J];中国原子能科学研究院年报;2002年00期

7 李江飞;贺晓;徐康泰;翟克平;李智博;罗璇;邹振春;王振业;;基于局部解析解的非稳态导热计算[J];科技通报;2017年04期

8 张学元;两类时滞微分方程的初等解析解[J];湖南纺织高等专科学校学报;1999年01期

9 马宝东,司建国;广义Jabotinsky方程的解析解[J];滨州师专学报;2001年02期

10 司建国;两类线性泛函微分方程解析解的存在性[J];应用数学学报;1991年02期

相关会议论文 前10条

1 黄拳章;武文明;王哲君;强洪夫;李艾华;;关于含单个流体夹杂无限大问题解析解的讨论[A];中国力学大会-2015论文摘要集[C];2015年

2 薛定宇;白鹭;;多解非线性矩阵方程的准解析解方法[A];系统仿真技术及其应用（第16卷）[C];2015年

3 陶震;倪明玖;;矩形槽道热电磁流体流动解析解[A];中国力学大会——2013论文摘要集[C];2013年

4 王耀;杨在林;;非均匀介质波动问题解析解答的等效转化关系研究[A];中国力学大会-2015论文摘要集[C];2015年

5 邱吉宝;宋海燕;程昊;秦朝红;刘张君;;热过屈曲梁振动的解析解[A];中国力学大会-2015论文摘要集[C];2015年

6 王虎;胡长顺;王秉纲;;弹性地基上无限长板受局部荷载作用的解析解[A];第八届全国结构工程学术会议论文集（第Ⅰ卷）[C];1999年

7 傅向荣;赵东;蒋东英;王宏志;赵阳;田歌;张浩;;楔形板的解析解及其数值解答探讨[A];北京力学会第十六届学术年会论文集[C];2010年

8 靳庆华;张明利;;Isogrid稳定性解析解分析[A];北京力学会第15届学术年会论文摘要集[C];2009年

9 周海兵;;拉氏数值模拟程序验证过程的人为解析解方法研究[A];中国工程物理研究院科技年报/2011年版[C];2011年

10 任青文;徐炜;曹茂森;;基于Lipschitz指数的梁表面裂缝量化[A];中国力学大会——2013论文摘要集[C];2013年

相关博士学位论文 前10条

1 苟晨华;能源利用问题的代数显式解析解[D];中国科学院研究生院（工程热物理研究所）;2007年

2 卢萌盟;复杂条件下复合地基固结解析理论研究[D];浙江大学;2009年

3 许韬;符号计算与若干非线性发展方程的解析研究[D];北京邮电大学;2011年

4 王慧平;联合荷载作用下海底管道的弯曲破坏特性研究[D];大连理工大学;2016年

5 陈晨;火电机组自动发电控制系统非线性特性与控制研究[D];北京交通大学;2016年

6 郭霄;复杂条件下竖井地基固结解析理论研究[D];浙江大学;2016年

7 周明;氯盐环境下混凝土耐久性分析与定量设计[D];广西大学;2016年

8 程瑶;基于结构元理论的煤层气输运模糊数学模型及解析解[D];辽宁工程技术大学;2011年

9 王晖;无线传感器网络生命周期跨层优化研究[D];上海交通大学;2008年

10 王玉林;承压地下水开采井流模型及其渗流理论研究[D];浙江大学;2011年

相关硕士学位论文 前10条

1 王旭;高压储氢罐充放气过程的热效应模拟与性能预测[D];武汉理工大学;2018年

2 李元媛;固体力学、流体力学及管理学方面的代数显式解析解[D];西北大学;2009年

3 程胜华;二维Eshelby夹杂问题解析解的数值验证和拓展研究[D];南昌大学;2016年

4 赵云;一类Volterra积分—微分方程的解析解研究[D];西安理工大学;2017年

5 李远;任意形状隧洞的粘弹性位移解析解[D];华北电力大学(北京);2017年

6 杜红梅;带有微扰动的Schr（？）dinger方程的解析解和数值解的比较[D];东北师范大学;2005年

7 白淑婷;两个新非线性可积方程的解析解[D];内蒙古师范大学;2017年

8 杨静;缓坡方程与修正型缓坡方程模拟一维分片光滑地形上波浪反射的准确解析解[D];广西民族大学;2011年

9 任晔;单裂隙岩体渗流与传热耦合的解析解与参数敏感度分析[D];北京交通大学;2009年

10 刘媛媛;超静定曲线梁面内变形解析解答[D];大连海事大学;2016年

本文编号：2756311