首页/文章/ 详情

关于沸腾问题汽化核心的模拟思考

9月前浏览475

正文共: 747字 2图     预计阅读时间: 2分钟

1 前言

在笔者看来,沸腾的模拟是非常困难的,因为没有一个比较完美的沸腾模型,通常各种模型参数多多少少都带有经验性。比如最为常见的Lee模型,有一个非常关键的相变系数(1/s),这个系数需要通过实验进行调整,通常的使用范围会相差数个量级。还有一个非常重要的问题,那就是汽化核心的模拟,我们知道汽化核心是产生气泡的种子,其尺寸通常在1~100μm级别。如果要模拟真实的汽化核心,就需要将这个几何细节也创建出来,这对于计算资源而言是一个挑战,而且当尺度小到一定程度时,可能NS方程都不适用了,这个因素我们暂且不考虑。Lee模型在处理源项时,源项的产生仅仅与温度有关,当温度高于(汽化)或低于(冷凝)饱和温度时,边产生了源项。而对于沸腾而言,相变是出现在固液界面处,也就是说在汽化核心处只要温度达到了一定条件就产生源项即产生了气泡,其他的传热传质发生在气液界面处。因此,笔者想通过这样一种不成熟的方式来处理,在汽化核心处采用常规的Lee模型(沸腾),而在其他的位置采用修改后的Lee模型(蒸发),如下图所示。这两个模型的UDF都在我们以往的案例中演示了虽然我们知道通常采用初始气泡的方式来模拟汽化核心,这种方式总体来说还是具有作弊性质。注意,对于下图的案例,如果在整个流体域直接采用常规的Lee模型,则在整个流体域特别是靠近固体的温度相高的位置随时出现气泡,当我们想关注气泡的产生和发展,这种方式就无法满足模拟要求了。
 
 

2 计算结果

采用上述的方法,我们做了一个二维的模拟案例,气泡的生成动画如下,看起来效果还不错,当然了由于计算资源有限,网格不够细,气液界面不够光滑漂亮。这里需要着重指出,汽化核心的尺寸为100μm,整个计算域宽度只有3mm,计算时间步长1e-7s,迭代计算收敛效果还是不佳。
 
 


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

蒸汽蓄热器放汽闪蒸过程模拟

正文共:917字6图预计阅读时间:3分钟1前言蒸汽蓄热器是一种以水为介质存储蒸汽热能的压力容器,其工作原理主要基于饱和水的显热储存和压力变化下的汽液相变转化。释放能量时,饱和水因压力骤降成为“过热水”,部分水瞬间闪蒸为低压蒸汽排放。今天,我们模拟一下这个闪蒸过程。2建模与网格创建如下的二维平面模型,包括罐体和一个排放口。划分四边形结构化网格,节点数5459。3边界条件与求解设置我们只描述一下几个关键的设置。闪蒸现象涉及到两相流,我们采用VOF模型,主相为液态水,次相为水蒸气。考虑重力影响。我们采用理想气体模型模拟水蒸气的密度变化,这个很重要,必须能体现压力和密度的对应关系。由于采用了理想气体模型,我们将参考压力设置为0Pa。我们采用软件自动的蒸发凝结模型来模拟闪蒸相变,饱和温度通过UDF定义为压力的函数,代码如下。DEFINE_PROPERTY(saturation,c,t){realp_ab;realp_op=RP_Get_Real("operating-pressure");realt_sat;p_ab=p_op+C_P(c,t);t_sat=1730.63/(8.07131-log10(p_ab/133.322))-233.426+273.15;returnt_sat;}排放口设置为压力出口,表压通过UDF定义为指数下降规律,代码如下,回流相为100%水蒸气,回流温度为罐内压力对应的饱和温度。DEFINE_PROFILE(out_pressure,thread,i){face_tf;realt=RP_Get_Real("flow-time");begin_f_loop(f,thread){F_PROFILE(f,thread,i)=1.27e6*exp(-5*t);}end_f_loop(f,thread)}罐内初始压力1.27MPa,初始温度为对应压力下的饱和温度463.6921K,采用patch设置罐内一半液位。4计算结果模拟时间内,罐内的液位变化动画如下,可以看到由于饱和水的闪蒸产生的蒸汽推动了液位上升,这也是所谓的“虚假水位”现象。当然,由于我们的模型是随意创建的,再计算下去可能液态水都排出了。但是,上述方法可适应实际工程模拟。来源:仿真与工程

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