水平集方法模拟沸腾的comsol实现（1）——计算模型

澹之

已于 2024-06-20 11:18:33 修改

阅读量347

点赞数 5

分类专栏：相变文章标签：算法

于 2024-06-20 11:02:15 首次发布

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接：https://blog.csdn.net/weixin_67664336/article/details/139365583

版权

相变专栏收录该内容

1 篇文章 0 订阅

订阅专栏

冷凝是一个复杂的气液相变问题，涉及两相流模拟，热量传递，质量传递。

相变模型

质量传递考虑温度饱和模型[1]：

$m=C\rho_l\ \ (T-T_s)/T_s$

m表示质量传递速率， $\rho_l$ 表示液相（水）的密度， $T_s$ 为饱和温度，C是一个比例常数，通过调整C维持界面温度为饱和温度。

控制方程[2]

连续性方程

质量蒸发（冷凝）速率可以公式化为界面法向通量，其用于界面上的连续性方程源项。

$m=\rho_f\ (U-u_f)\bullet\mathbf{n}$

$\rho_f$ 为流体密度，U为界面速度， $u_f$ 为流体速度。

液相气相的速度表示为

$u_l=U-m/\rho_l$

$u_v=U-m/\rho_v$

在界面处取液相和气相速度的差值，就可以得到由相变和相应的密度变化所产生的散度。

$u_l-u_v=(1/\rho_v\ -1/\rho_l\ )m$

连续性方程的源项等于上式乘以一个只允许源项在界面上非零的界面delta函数。

$\nabla\bullet u=\delta(u_l-u_v)$

$\nabla\bullet u=\delta(1/\rho_v\ -1/\rho_l\ )m$

动量方程直接考虑不可压缩流体的NS方程

传热方程

$C_p\ \ \partial T/\partial t+\rho C_p\ (u\bullet\nabla)T-\nabla\bullet k\nabla T=S_h$

其中Sh=-δmhlg 是热源，hlg 是潜热

水平集方程

因为Level set 函数默认情况下只考虑界面的平流，如下图，它不考虑生成液相时的界面移动，所以也需要对其进行修改。

水平集方法（comsol）模拟气泡上升

水平集方程[3]

$\partial\phi/\partial t+U\bullet\nabla\phi=0$

考虑相变后

$U=u_f+m/\rho_l$

所以在level set函数中添加的源项为 $m/\rho_f\ \ \mathrm{\nabla\phi}$ ，修改之后的水平集方程便可以考虑蒸发冷凝带来的界面变化。

参考文献：

[1]Tsujimoto, K. , Kambayashi, Y. , Shakouchi, T. & Ando, T. (2009). Numerical simulation of gas-liquid two-phase flow with phase change using cahn-hilliard equation. Turbulence, Heat and Mass Transfer, 6, 1-12.

[2]Forster C J, Smith M K. The Transient Modeling of Single-Bubble Nucleate Boiling in a Sub-Cooled Liquid Using an ALE Moving Mesh[C]//Excerpt from the Proceeding of the 2011 COMSOL Conference in Boston. 2011.

[3] Olsson E, Kreiss G. A conservative level set method for two phase flow[J]. Journal of computational physics, 2005, 210(1): 225-246.

关注

5
点赞
踩
4

收藏

觉得还不错? 一键收藏
0
评论
水平集方法模拟沸腾的comsol实现（1）——计算模型

冷凝是一个复杂的气液相变问题，涉及两相流模拟，热量传递，质量传递。
复制链接

扫一扫

专栏目录

评论

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。