首页/文章/ 详情

FLUENT气瓶压力热胀冷缩效应模拟

10月前浏览645
正文共:1064字 11图    预计阅读时间:3分钟
 
 

1 前言

 
 
    本案例灵感源自现场实验现象,本周冷空气来袭,气温发生了明显的降低。笔者在为发电系统更换气源(气瓶)时,按照以往的常规操作方法,发现新的气瓶换上后,压力在持续缓慢下降,怀疑有泄漏,但是怎么都没找到漏点。当时正值下午,气瓶间开始摆脱太阳光的照射,因此我们怀疑是气温造成的。由于环境温度下降,使得气瓶内的温度下降,根据理想气体状态方程pv=nRT,温度下降对应压力也降低。经过了十几个小时的观察,最终确认是温度影响造成该现象。下图为观察期间的压力变化情况(2020年12月30日20点至次日早上8点)。本案例模拟的是密闭气瓶热胀冷缩对内部压力的影响。
 
 

2 建模与网格

 
 
    建立如下的二维气瓶模型,气瓶内径400mm,壁厚5mm,考虑对称性,只建一半模型。划分全六面体网格,节点数9119,最小正交质量0.66。
 
 

3求解设置

 
 
    开启能量方程。
    由于本案例的计算结果对流场不敏感,湍流模型选择默认即可。
    气体介质选择甲烷,密度模型选择理想气体(这个设置非常重要)。另外,本案例并不考虑气体本身的温度分布,即采用集总参数思想处理气体温度,因此把导热系数设置成很大值,这样气体的温度就是均匀的(这样设置是为了减小计算量,做了简化)。
    气瓶初始温度10℃,气瓶外壁为等温壁面,按照如下的升温速率升温。
    DEFINE_PROFILE(mytemp, thread, position){face_t f;real t;t=RP_Get_Real("flow-time");begin_f_loop(f,thread){F_PROFILE(f,thread,position)=283.15+100*t;}end_f_loop(f,thread)}
        本案例不考虑重力场,这样内部的气体基本保持静止,这样设置也是为了减小计算量,做了简化。参考压力1MPa,也即气瓶内部初始压力1MPa,初始化时表压为0Pa即可。
        按如下进行初始化
        瞬态求解,时间步长0.001s,计算时间步数100,这样气瓶表面温度将从10℃经过0.1s上升到20℃。
     
     

    4计算结果

     
     
        气瓶内甲烷的平均温度随时间变化如下,从10℃上升到11.504373℃。
        气瓶内甲烷的平均压力(绝对压力)随时间变化如下,从1000000Pa上升到1005313Pa。
        根据理想气体状态方程,初始状态和终了状态的气体压力和温度满足P1/T1=P2/T2,据此可计算得到温度为11.504373℃时的气体压力为
    P2=T2/T1*P1=(284.654373÷283.15)×1000000Pa=1005313.99Pa,与FLUENT计算结果一致。
    来源:仿真与工程
    湍流
    著作权归作者所有,欢迎分享,未经许可,不得转载
    首次发布时间:2023-07-05
    最近编辑:10月前
    余花生
    签名征集中
    获赞 165粉丝 267文章 278课程 0
    点赞
    收藏
    未登录
    还没有评论

    课程
    培训
    服务
    行家

    VIP会员 学习 福利任务 兑换礼品
    下载APP
    联系我们
    帮助与反馈