首页/文章/ 详情

FLUENT高压氢气泄漏自燃模拟

2天前浏览9

正文共: 1297字 13图     预计阅读时间: 4分钟

1 前言

高压氢气泄漏到空气中自动点燃(也称为自燃或自发点火)是一个复杂但已被实验和理论证实的现象。它主要发生在高压(通常几十兆帕以上)氢气突然泄漏到大气中的情况下。发生自动点燃的核心原因是泄漏瞬间局部区域达到了氢气的自燃温度(约550°C)当高压氢气从容器或管道中突然泄漏(例如阀门失效、管道破裂)时,氢气会以极高的速度喷射出来。这股高速喷射的氢气流会猛烈地冲击和压缩前方的静止空气。根据热力学原理(理想气体定律和绝热压缩过程),气体被快速压缩时,其温度会急剧升高。计算表明,在足够高的初始压力下(例如35MPa或更高),被压缩的空气团温度可以瞬间超过550°C,达到甚至超过1000°C。当这股高温压缩空气与随后喷出的新鲜氢气接触时,只要混合物的浓度在可燃范围内(氢气在空气中的可燃极限是4%-75%体积比),就会立即发生点燃。这个现象是氢安全研究的一个重要方向,今天我们以FLUENT模拟一下这个过程。

2 建模与网格

创建如下的二维模型[1]划分四边形结构化网格,注意下游管道和高压储氢室之间的爆破片要创建实体边,便于后续模拟爆破过程,也就是说整个计算域被爆破片分为两个子域。
 
 

3 边界条件与求解设置

为了观察管道中湍流的精细结构,本案例我们采用二维LES大涡模型模拟湍流,氢气的燃烧采用EDC模型,这两种模型被验证了用于高压氢气泄漏自燃工况的可行性。注意,对于二维模型,LES模型需要用TUI启动,语句为(rpsetvar 'les-2d? #t)。氢气的燃烧我们采用氢气-空气单步模型,燃烧反应的指前因子、活化能、温度指数等采用软件默认值。这种高压氢气泄漏流动特别还包含燃烧、爆炸现象等通常必须考虑气体可压缩性,因此混合物密度采用理想气体模型,而不能用默认的不可压缩理想气体模型。
 
 
 
 
 
 
 
 
用于采用了理想气体模型模拟可压缩流动,通常参考压力设置为0Pa。那么,我们在输入压力边界条件时,压力值应为绝对压力。
 
 
下游管道端口为压力出口边界,表征外界环境,亦即回流组分为空气(21%氧气+79%氮气)
 
 
 
 
我们采用动网格的event来模拟爆破过程,即在某个时刻将边界条件由壁面修改为内部边界,注意这个模拟还是偏激进的,因为表征的爆破耗时为0s,但实际上即使爆破再快,也是需要时间的,可能是毫秒或者微秒级别。不过,这样的设置通常对于模拟研究是可以接受的。
 
 
以下列条件初始化,注意必须还要用patch功能修改两个子域的组分、压力初始值,毕竟初始时刻下游管道内全部是空气,高压储氢室压力为15MPa(本案例假定值)。
 
 
瞬态求解,时间步长1e-8s

4 计算结果

我们看一下某个时刻下游管道内的温度、水蒸气分布,可以看出氢气发生了自燃。
 
 
 
 
我们看一下下游管道距离爆破片40mm、80mm和120mm处的压力曲线,可以看到压力波的传播特性,我们将两点距离除以压力传播的时间间隔即可得到压力波的传播速度,本案例压力波的速度约为1900m/s。
 
 
参考文献
[1] 不同泄放条件下高压氢气泄漏激波传播行为及诱导自燃特性
[2] 高压氢气泄漏自燃现象的模拟

来源:仿真与工程
Fluent燃烧动网格湍流理论爆炸管道
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-07-17
最近编辑:2天前
余花生
签名征集中
获赞 234粉丝 515文章 337课程 0
点赞
收藏
作者推荐

FLUENT河流弯道二次流模拟

正文共: 662字 8图 预计阅读时间: 2分钟1 前言河流在弯曲的河道中流动时会产生二次流,使得弯曲河岸外侧(凹岸)被水冲刷,内侧(凸岸)有沉积物堆积[1]。今天我们用一个三维案例简单模拟一下这个现象。2 建模与网格根据参考文献[1]创建如下的三维河道模型,划分六面体结构化网格,单元数19.2万,最小正交质量0.8。 3 边界条件与求解设置水沙运动是两相流运动,本案例的核心设置是采用欧拉双流体模型将水和沙子均看作是流体,认为二者之间相互穿透,用体积法分数来描述每一相的存在[2]。主相为水,次相为沙子,我们偷懒直接修改空气的物性参数,比如沙子密度2650kg/m3,直径0.1mm。其他的参数我们先不修改,采用默认值,将来有精力再做系统性研究。 考虑水深方向的重力,务必不能忘了。 河道入口设置为速度入口,两相的速度均为0.5m/s,其中沙子的体积分数0.1%。 出口设置为压力出口,表压0Pa。两侧和底面均为无滑移壁面。顶部为对称边界。这是一种简化的理想模型处理,这样就省略了顶部的空气的影响。本案例的湍流模型设置如下,可能不是最合适的,可以考虑RSM模型。 欧拉模型的稳态求解可能较难达到良好的收敛结果,比如本案例,如果采用稳态计算,残差有点难看。可以这样处理,先稳态求解获取初始值再进行瞬态求解,时间步长0.5s,最后可以获得很好的收敛解,残差图如下。 4 计算结果我们看一下弯头几个截面的沙子体积分数分布,可以明显看出河道凸岸这一侧的沙子沉积现象。这一点也可以从沙子体积分数等值面看出。 参考文献[1] 基于CFD的河流弯道二次流的三维数值模拟[2] 风沙运动的欧拉双流体模型参数研究来源:仿真与工程

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈