车用液氢气瓶快充过程数值研究.pdf
《车用液氢气瓶快充过程数值研究.pdf》由会员分享,可在线阅读,更多相关《车用液氢气瓶快充过程数值研究.pdf(8页珍藏版)》请在咨信网上搜索。
1、车用液氢气瓶快充过程数值研究金树峰,王崇宇,姚淑婷,谭风光(兰州理工大学石油化工学院,兰州730050)摘要:车用液氢气瓶加注过程中,因加注速度快且气瓶容积小,极易造成过充现象,导致安全事故发生。针对设有气包防过充装置的车载液氢气瓶建立三维对称模型,采用 CFD 方法对气瓶快速加注过程开展数值仿真。分析了开孔孔径、加注速率、初始充满率和防过充装置容积对防过充装置性能的影响规律。结果表明:开孔孔径是影响防过充装置性能的主要因素,其限制了液体的流通能力;采用 6mm 孔径和初始充装量为 50%的状况下,充装结束后气瓶都将过充,最大充装率可达 98.4%;采用大容积的防过充装置可以为气瓶预留更多的气
2、相安全空间,但同时也会导致充装量不足等问题,当防过充装置容积为 100L 时,最终充装量仅有 77.2%。关键词:车载液氢气瓶;液氢加注;防过充装置;充装率中图分类号:TB657.9文献标志码:A文章编号:10067086(2024)02018008DOI:10.12446/j.issn.1006-7086.2024.02.011Numerical Study on Fast Filling Process of Vehicle Liquid Hydrogen CylinderJIN Shufeng,WANG Chongyu,YAO Shuting,TAN Fengguang(School o
3、f Petrochemical Engineering,Lanzhou University of Technology,Lanzhou730050,China)Abstract:Thevehicleliquidhydrogencylinderiseasytobeoverfilledduetothefastfillingspeedandthesmallsizeofhydrogencylinderduringthefillingprocessandsafetyaccidentsmightbecaused.Athree-dimensionalsymmetricmodelisbuiltforvehi
4、cleliquidhydrogencylinderthatisequippedwithanullagetank,andanumericalstudyusingCFDsimulationmethodisperformedonthefastfillingprocessofvehicleliquidhydrogencylinder.Theeffectsofholesizes,fillingrate,ini-tialfillingrateandvolumeofoverfillpreventiondeviceontheperformanceoftheoverfillpreventiondevicesar
5、estudied.Theresultshavebeenobtained.Theholesizeshaveasignificanteffectontheperformanceoftheoverfillpreventiondevice,whichlimitstheflowcapacityofliquid.Whentheholesizeis6mmandtheinitialfillingrateis50%,thecylinderswillbeoverfilledattheendoffilling,andthemaximumfillingratecanreach98.4%.Theuseoflargevo
6、lumeoverfillpreventionde-vicescanreservemorephaseforcylinders,butitwillalsoleadtoissuessuchasinsufficientfilling.Whenthevolumeofover-fillpreventiondeviceis100L,thefinalfillingrateisonly77.2%.Key words:vehicleliquidhydrogencylinder;liquidhydrogenfilling;overfillpreventiondevices;fillingrate0引言低温容器加注过
7、程与一般的常温液体容器加注有许多不同之处。低温液体受热容易膨胀,而液体具有不可压缩性,所以气瓶必须留出一定的气相空间,使介质受热后有膨胀空间1-2。如果加注过量,气相空间不足,随着外部热量进入气瓶,瓶内的压力将急剧升高。一旦气瓶内的压力大于气瓶的爆破压力,则会造成安全事故3。因此,须确保气瓶加注完成后留有一定的气相空间,通常会在瓶内设置防过充装置,如溢流阀和气包装置等,确保气瓶的安全充装。收稿日期:2023-11-08基金项目:甘肃省青年基金(21JR7RA269)作者简介:金树峰,副教授,主要从事低温储运技术研究。E-mail:引文信息:金树峰,王崇宇,姚淑婷,等.车用液氢气瓶快充过程数值研
8、究J.真空与低温,2024,30(2):180187.JINSF,WANGCY,YAOST,etal.NumericalstudyonfastfillingprocessofvehicleliquidhydrogencylinderJ.VacuumandCryogenics,2024,30(2):180187.真空与低温第 30 卷第 2 期180VacuumandCryogenics2024年3月针对低温液体的加注过程已有大量研究。赵金鑫等4建立了液氦罐箱的一维非稳态热力学模型,研究了加注流量对大型液氦罐箱加注过程的影响,并将其加注过程分为预冷、积液和静置三个阶段。Kang 等5研究了车用液
9、氢气瓶的初期加注过程,发现在气瓶进出口管道处出现了涡流。马原等6建立了液氮立式容器模型,并对其加注过程中产生的降温及热力特性进行研究,所建立模型能够较好地预测容器加注过程产生的瞬态降温及热应力。Ma 等7利用数值模拟的方法研究了微重力条件下液氢储罐的无排气充注过程,构建了二维轴对称模型,并通过实验数据进行验证,对比分析了变重力条件下的无排气填充过程,并分别讨论了入口结构、入口液体温度、初始壁温和入口流量四个因素的影响。Wang 等8还使用相同的实验台对储罐的垂直和水平放置进行了无排气填充对比实验,结果表明在相同工况下卧式储罐在无排气填充中的表现优于立式储罐。Wang 等9针对 LNG 无排气加
10、注过程中储罐的热力学过程进行了建模和分析,模拟了低温储罐内部的复杂传热和传质现象。Kim 等10对四氟甲烷的无排气加注进行了实验研究,实验表明加注液体温度和初始壁面温度是影响无排气加注性能的主要因素。Ferrn 等11针对 0.5m3立式低温储罐建立了二维 CFD 模型,探究了不同初始充装率和绝热层厚度对低温储罐内 LNG 的压力、温度以及蒸发率的影响。Choi 等12建立了 4.9m3低温储罐的 CFD 模型,存储介质为 LNG 和液氮,研究了不同充注率情况下储罐内部低温流体出现的热分层现象。Jo 等13建立了非平衡状态 LNG储罐数学模型,并进行动态模拟计算,探究了不同初始充装率对储罐内压
11、力、温度等热力学参数的影响。现有文献关于车载液氢加注系统,以及车载液氢气瓶内气包防过充装置的作用过程和机理的相关研究和实验报道较少。本文针对设有气包防过充装置的车载液氢气瓶加注过程进行数值模拟分析,以期了解液氢加注过程气瓶内部的热力学变化过程,以及气包防过充装置的作用过程,为气包防过充装置的安全设计提供理论依据和支撑。1模型建立1.1物理模型以车载液氢气瓶为研究对象,总容积 500L,具体结构如图 1 所示14。为了便于分析与计算,本文对车用液氢气瓶进行了合理简化,只针对气瓶内胆及其内部流体域进行建模。入口底部开有小孔防过充装置图 1车用液氢气瓶结构示意图Fig.1Schematicofveh
12、icleliquidhydrogencylinder气包防过充装置由防过充筒体和防过充封头组成,与后封头焊为一体,形成一个独立的腔室,其容积为气瓶容积的 10%,且在防过充筒体底部开有小孔。其作用原理是首次充装内筒腔室满液时,由于气包防过充筒底部开孔孔径很小,在短时间内流入气包的液体量很少,从而保留一定的气相安全空间。加注完成后随着时间延长,内筒腔室内液体会通过小孔逐渐流入气包内,达到内外液位平衡,从而确保气瓶的气相安全空间15。以气包底部开孔处为断面,液体流入气包的速率可由式(1)计算:Qm=AC02p p0+ghL(1)式中:为泄漏液体密度;A 为开孔孔径面积;C0为液体泄漏系数;p 为气
13、瓶内的压力;p0为防过充装置内的压力;g 为重力加速度;hL为内筒体与气包的液位差。1.2CFD 数学模型鉴于气瓶结构的对称性,建立三维对称模型。采用 Fluentmeshing 软件对车用液氢气瓶进行网格划分,对气包开孔区域作网格加密处理,如图 2 所示,并进行网格无关性验证,最终选择总数约为 45 万的网格进行计算。图 2CFD 模型与网格划分Fig.2CFDmodelandmeshgeneration对于车载液氢气瓶加注过程中流体区内的气液两相流动,采用 VOF 多相流模型获取气液相界金树峰等:车用液氢气瓶快充过程数值研究181面情况。为了描述气瓶内的相变现象,采用了 Lee模型,相变强
14、度控制因子 C 取为 0.1。在该模型中,将每个流体单元温度 T 与饱和温度 Tsat进行比较,如果 TTsat,液体蒸发,而如果 TTsat,蒸汽冷凝。相变过程中的传质可通过以下公式计算,能量源项Q 的绝对值为相变质量与气化潜热的乘积。Sg=Sl=Cll|T Tsat|Tsat(2)Sl=Sg=Cgg|T Tsat|Tsat(3)Q=Slhfg(4)式中:T 和 Tsat分别为流体温度和饱和温度;和 hfg为相变的体积分数和潜热;为流体密度;S 和 Q 分别为质量转移速率和能量源项;C 为相变强度控制因子;下标l 为液相,g 为气相。在液氢气瓶加注过程中,气瓶内存在明显的温度梯度,因此需考虑
15、流体物性参数随温度变化。对于气体,密度采用理想气体模型,其他物性设为定值。液体视为不可压缩,其密度随温度的变化,通过查询 NIST 数据库中液氢的物性数据在 Fluent 软件中进行插值分段线性拟合。气瓶加注过程入口边界条件选择质量流量入口,入口液氢温度为 20K,初始时刻气液相处于饱和状态。因加注过程时间较短,且气瓶绝热性能较好,气瓶外壁设置为绝热边界条件,内壁设置为耦合边界类型。采用 PISO 算法计算,动量、能量方程采用二阶迎风格式,时间步长设置为 0.001s。2结果与讨论2.1加注过程瓶内流场分布及压力变化加注流量取 0.05kg/s,进液温度 20K,初始充装率为 5%(液位处于防
16、过充装置下方),瓶内初始压力 0.1MPa,瓶内初始液氢温度为该压力下的饱和温度 20.3K,对气瓶加注过程进行数值模拟。以防过充装置开孔孔径 1mm 和 6mm 为例,图 3 给出了气瓶加注过程的气液相分布云图,当体积分数为 1 时,表示此区域全为氢蒸气,当体积分数为 0时表示此区域全为液氢。加注约 20s,液面开始漫过防过充装置底部,逐渐有液体通过小孔流入防过充装置。当孔径为 1mm 时,因开孔孔径较小,流入防过充装置内的液体较少。但随着加注进行,由于气瓶内腔室液体不断蒸发及低温液体液位不断升高对气相空间的挤压,气瓶内腔室压力不断升高。其与防过充装置内部将存在很大压差,液体进入防过充装置。
17、但因孔径较小,液体流通能力较弱,加注完成时只有少部分液体进入防过充装置。当开孔孔径为 6mm 时,由于孔径较大,开孔处液体流通能力较强,加注完成时液体几乎充满防过充装置。1.00.90.80.70.60.40.30.20.1015 s150 s(a)孔径 1 mm(b)孔径 6 mm280 s15 s150 s280 s气相体积1.00.90.80.70.60.40.30.20.10气相体积分数分数图 3加注过程中气瓶气液相分布云图Fig.3Graphiccontoursofphasedistributionincylinderduringfilling在加注结束后防过充装置顶部气体受挤压温度
18、和压力升高,防过充装置内部和气瓶内部的低温液氢具有一定过冷度,分别通过热对流与热传导的方式对防过充装置顶部气体进行冷凝降温降压。随着时间持续,受冷凝作用影响,防过充装置顶部气体质量变少,压力和温度也逐渐降低,气相空间182真空与低温第30卷第 2 期体积减小。随着时间进一步持续,氢气与液氢的温差将变得很小,液氢的冷凝作用减弱,且伴随着外部环境热量的渗入及瓶内压力降低饱和温度降低,液氢也在蒸发,气瓶内部流体会处于不断变化的状态,防过充装置顶部将始终存在一小部分气相空间。气瓶加注过程压力场分布云图如图 4 所示,开孔孔径为 1mm 时,在加注过程中气瓶内筒体部分的压力受液体的蒸发及气相空间的压缩影
19、响不断上升,气相区域各处压力值差异较小。而防过充装置底部因开孔孔径较小,气液流通能力较弱,其内部压力上升较慢,在加注结束时气瓶内部与防过充装置内部的压差达到了 1.06MPa。当防过充装置的开孔孔径为 6mm 时,由于开孔孔径较大,液体流通能力较强,在加注过程中大量液体不断地流入防过充装置内部,气瓶内筒体部分与防过充装置内部虽然存在压差但其值很小。在加注结束时,由于液柱静压的原因,液相区域压力值向自上而下递增,低温液体底部出现最大压力值,气相区域压力与液相区域压力最大差值仅约为 360Pa。1.002105压力/Pa压力/Pa压力/Pa压力/Pa压力/Pa压力/Pa1.011051.02610
20、51.0381051.051051.151051.251051.351051.451052.01054.01056.01058.01051.01061.21061.0331051.0351051.0371051.0391051.411051.4221051.4341051.4461051.2011061.201 11061.201 21061.201 310615 s150 s(a)孔径 1 mm(b)孔径 6 mm280 s15 s150 s280 s图 4加注过程瓶内压力分布云图Fig.4Graphiccontoursofpressuredistributionincylinderduri
21、ngfilling2.2孔径及初始充装率对最终充装率的影响加注流量取 0.05kg/s,进液温度 20K,初始充装率为 5%(液位位于防过充装置下方)和 50%,瓶内初始压力 0.1MPa,瓶内初始液氢温度为该压力下饱和温度 20.3K,开孔孔径分别取 1mm、2mm、3mm、6mm,对气瓶加注过程进行数值模拟。图 5和图 6 给出了加注过程气瓶与防过充装置内部压差随时间变化曲线和加注过程防过充装置内液体容积变化曲线。由式(1)可知液体流入防过充装置的速率是由气瓶内筒体和防过充装置内部的压差及液位差决定的。由图 6 中防过充装置内液体容积变化曲线可以看出,随着加注进行,曲线上的点切线斜率逐渐增
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 液氢 气瓶快充 过程 数值 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。