首页/文章/ 详情

管路中甲烷空气扩散混合模拟

11月前浏览1970

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

1 前言

今天,我们研究一个实际应用问题。两段管分别通入空气和甲烷,并在下游的母管汇合排入大气环境。假设某个时刻,空气和甲烷入口突然关闭,显然接下来空气和甲烷将通过扩散进行混合,我们关注甲烷入口什么时候会出现显著量的空气(亦即甲烷浓度显著降低,并且按两种情况考虑,第一种情况:出口未关闭;第二种情况:出口也随即很快关闭。对于第一种情况,相当于还有一个空气源,源源不断进入管道系统,足够久时间管道系统内应该接近于完全由空气占据;第二种情况,相当于形成一个封闭的空间,管道系统内残余的空气和甲烷进行混合至均匀。

2 建模与网格

创建如下的二维平面模型,划分四边形结构化网格。
 
 

3 边界条件与求解设置

湍流模型我们选择k-e两方程模型,按如下设置。
 
 
开启组分输运模型,组分为甲烷和空气,物性参数采用默认设置,需要指出的是,本案例的一个关键参数是扩散系数,这会直接影响到计算结果,甲烷在空气中的扩散系数应该有比较权威的实测值。
 
 
甲烷入口和空气入口均设置为速度入口,速度分别为1m/s和5m/s,组分分别为100%甲烷和100%空气。
 
 
 
 
出口为压力出口,回流组分为100%空气,表压为0Pa,表征外部大气环境。
 
 
开启重力选项,这个设置对于模拟实际情况也是相当重要的,因为把密度差引起的自然对流也考虑进去。需要指出的是,实际上温度场也很重要,倘若计算域内的不同位置具有显著的温度差异,则还需要考虑温度的影响,因为温度影响密度从而影响对流。说到底,这是一个自然对流、扩散的综合问题。当然,考虑的因素越多,计算就越复杂,对于工程实际,这些因素要具体问题具体分析。
 
 
我们采用动网格的event来模拟入口和出口的关闭,即在某个时刻将原先的边界条件改为壁面边界。本案例工况1:在0.5s时将甲烷入口和空气入口改为壁面;工况2:在0.5s将所有入口和出口改为壁面。
 
 
先采用稳态求解获得初始速度场和浓度场,在此基础上开展瞬态计算,时间步长取0.02s。
创建一个监视器,监测甲烷入口位置的甲烷平均浓度。

4 计算结果

首先,获得的初始稳态浓度场如下。
 
 
工况1,一分钟内甲烷入口的甲烷平均浓度曲线如下,可以看出约在10s后空气便扩散到了甲烷入口,后续空气含量越来越多,我们没有监测空气入口的甲烷含量,实际上该位置的甲烷也会越来会多。60s时,甲烷入口的平均浓度为98.69%,空气入口的甲烷含量为7.4%
 
 
 
 
工况2,一分钟内甲烷入口的甲烷平均浓度曲线如下,同样可以看出约在10s后空气便扩散到了甲烷入口,后续空气含量越来越多。60s时,甲烷入口的平均浓度为98.35%,空气入口的甲烷含量为7.8%
 
 
 
 
本案例讨论的问题在SOFC系统测试中具有很强的指导意义,在系统出现急停的时候,如果未在燃料管路继续通入还原性气体,氧气会往阳极扩散,可能会对电堆阳极带来氧化的风险。当然,这里有两个问题值得探讨,第一是在高温下的急停,氧气扩散的同时也会与燃料发生氧化反应,氧分子也在被消耗掉;第二是达到电堆阳极的氧气浓度未必短时间内很高,可能即使被氧化,也可以在后续的启动中进行还原,当然前提是氧化不是很严重,如果很严重则可能氧化镍会将电池片挤。另外,当扩散距离越远时,所需的时间越长,理论上在氧分子到达阳极这段时间内还有时间对系统进行应急处理,而这个时间是可以通过仿真计算进行估计的。总的来说,在系统急停时最好能第一时间通入还原性气体,这样阳极就不会存在氧化的风险了。


来源:仿真与工程
动网格组分输运湍流理论管道
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-09-09
最近编辑:11月前
余花生
签名征集中
获赞 259粉丝 606文章 392课程 0
点赞
收藏
作者推荐

FLUENT直接接触式换热器模拟之二

正文共:1690字13图预计阅读时间:5分钟1前言直接接触式换热器是一种高效的热能传递设备,其核心特点在于冷热流体直接混合接触进行热交换,完全省去了传统间壁式换热器中的固体传热壁面。这种设计消除了壁面热阻和污垢热阻,使得其传热效率极高,单位体积或单位面积的传热能力远大于间壁式换热器。常见的应用形式包括气液系统(如喷淋塔、填料塔用于气体冷却/加热、蒸汽冷凝、增湿/除湿)和液液系统(如搅拌混合罐)。它的工作原理主要依赖于两相流体的充分混合:一种流体被分散成液滴、气泡或薄膜(分散相),与另一种连续相流体在塔器或容器内逆流、并流或错流接触。巨大的、不断更新的相界面是高效传热的关键。其显著优势包括结构相对简单、成本较低、不易结垢、能处理含颗粒或易结垢流体、以及可利用相变潜热(如蒸汽直接接触冷凝)实现极高传热速率。然而,其应用也有严格限制:必须允许两种流体相互混合且不发生有害的化学反应或污染。流体间的互溶性、后续分离的难易程度以及允许的少量夹带是选型的关键考量。此外,流体力学行为(如气液系统中的液泛、压降)需要精确控制。尽管存在限制,在允许流体直接接触的工业场景(如电站冷却塔、化工吸收/解吸伴随传热、海水淡化、地热利用、冶金、食品加工等)中,直接接触式换热器凭借其卓越的效率和解决特殊工况的能力,是不可替代的重要选择。之前我们做一个简单的案例模拟了直接接触式换热器,彼时我们没有考虑相变,今天我们继续做一个案例,本案例考虑了高温空气与冷却水混合后的加湿过程。2建模与网格创建如下的三维混合换热器,具有一定湿度的高温空气(水蒸气质量分数0.00842)从换热器底部进,从顶部出,冷水从中间进,底部出。划分多面体网格,节点数43.6万,最小正交质量0.43。3边界条件与求解设置本案例需要同时用到多相流模型和组分输运模型,其中组分输运模型模拟的是空气的湿度变化,我们采用VOF两相流模型,液态水为主相,空气混合物为次相。同时,需要注意的是组分输运模型的混合物要和两相流的相对应上,空气混合物的组分顺序也要把水蒸气放在前面,因为本案例的UDF需要调用混合物的组分。基本的设置:开启能量项、开启重力项,不再赘述。热空气入口设置为速度入口,速度2m/s,温度100℃,组分为100%空气(其中水蒸气质量分数0.00842)。冷水入口设置为速度入口,速度0.5m/s,温度25℃,组分为100%液态水。热空气出口设置为压力出口,注意回流组分设置为100%空气,这个需要格外注意,如果误设置为液态水,则在迭代计算过程中若出现回流,可能造成计算发散。冷水出口设置为压力出口,回流组分设置为100%液态水。因为湿空气的湿度会随着换热的进行而发生变化,因此其露点温度是变化的,本案例我们采用UDF方法修改湿空气的露点温度,亦即两相流中液态水相变为空气混合物的水蒸气的饱和温度。该UDF代码与我们之前的案例(FLUENT湿空气冷凝模拟)是一致的,读者朋友可以从之前的案例获取,其他的一些设置参考也可以查看该案例。需要注意的是,经过多次调试该代码在本案例中无法用于干空气的加湿,也就是空气混合物的水蒸气质量分数不能太小甚至为零,否则无法收敛。本案例我们将相变系数减小为0.001/s,不然很容易计算发散,当然这个值应该从实验获取,通常设置成较大值,则瞬态计算时间步长需要调小,否则容易发散,反之亦然。为了获得较好的收敛解,本案例的亚松弛因子需做如下设置。采用patch方法设置换热器内初始液位,这个设置很重要,我们没有必要把注水这个过程也模拟,对于给定流量,最终会有一个终了液位,此时入口流量和出口流量相等。瞬态计算,时间步长0.0001s。4计算结果由于计算非常缓慢,本案例并没有计算到稳态。我们主要看一下出口的水蒸气摩尔分数变化情况,可以看出空气的湿度增加了。我们看一下某个时刻的露点温度分布,湿空气的露点温度约在11℃附近,液态水的露点温度亦即该压力下的饱和温度。再看一下计算域截面的水蒸气摩尔分数云图,水面上方附近的摩尔分数略高。最后,我们看一下相变速率云图,显然相变发生在空气和液态水交界面处,因此这些区域的湿度也会如上图所述更高一点。来源:仿真与工程

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