使用界面追踪对膜态沸腾的数值研究

计算物理 ›› 2008, Vol. 25 ›› Issue (5): 576-584.

使用界面追踪对膜态沸腾的数值研究

袁明豪, 杨燕华, 李天舒, 胡志华

上海交通大学核科学与工程学院, 上海 200240

收稿日期:2007-04-05 修回日期:2007-09-14 出版日期:2008-09-25 发布日期:2008-09-25
作者简介:袁明豪(1980-),male,wujiang,Jiangsu,PhD,interested in the numerical simulation of multiphase flow.
基金资助:
National Natural Science Foundation of China(No.50376036)

Numerical Study on Film Boiling with Interface Tracking

YUAN Minghao, YANG Yanhua, LI Tianshu, HU Zhihua

School of Nuclear Science and Engineering, Shanghai Jiaotong University, Shanghai 200240, China

Received:2007-04-05 Revised:2007-09-14 Online:2008-09-25 Published:2008-09-25
Supported by:
National Natural Science Foundation of China(No.50376036)

摘要/Abstract

摘要： 提出一种计算带相变的自由界面的数值算法.基于分段线性界面重构(PLIC)的VOF(volume-of-fluid)方法用于追踪自由界面,并对汽液交界面上的相变导致的不连续速度场给出处理方法.此方法容易实施且被证明是有效的.流场的求解使用SIMPLE方法,表面张力使用连续表面力模型(CSF)进行计算.在三维直角坐标系下,模拟了水平壁面上的膜态沸腾,在二维适体坐标系下,模拟了竖直圆头柱体表面的自然对流膜态沸腾.计算结果与理论关系式符合较好.

关键词: 膜态沸腾, volume-of-fluid, 相变

Abstract: A numerical method for simulation of interfacial flows with phase change is presented. A volume-of-fluid (VOF) method based on piecewise linear interface construction (PLIC) is used to track interface. Spocial treatment is taken to deal with discontinuous velocity field due to phase change at liquid-vapor interface. The treatment is easy and effective. SIMPLE method is used for flow field and a CSF model for surface tension. The method is used to simulate three-dimensional film boiling on a horizontal plate with Cartesian coordinates and axisymmetric two-dimensional natural convection film boiling on a hemisphericallyended rod with body fitted coordinates. They are in good agreement with analytical correlations.

Key words: film boiling, volume-of-fluid, phase change

中图分类号:

TL33

袁明豪, 杨燕华, 李天舒, 胡志华. 使用界面追踪对膜态沸腾的数值研究[J]. 计算物理, 2008, 25(5): 576-584.

YUAN Minghao, YANG Yanhua, LI Tianshu, HU Zhihua. Numerical Study on Film Boiling with Interface Tracking[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2008, 25(5): 576-584.

[1]	付韬, 邬龙, 李晨光. 基于生成函数方法的现实网络座键渗流建模[J]. 计算物理, 2022, 39(2): 212-222.
[2]	王子墨, 李凌. 超短激光打孔中快速相变的格子玻尔兹曼模拟[J]. 计算物理, 2020, 37(3): 299-306.
[3]	李琼, 刘海风, 张弓木, 张其黎. 模拟退火算法在化学自由能模型中的应用[J]. 计算物理, 2019, 36(3): 259-264.
[4]	刘昌宇, 李仔武, 李栋. 含相变材料玻璃类围护结构传热系数分析[J]. 计算物理, 2016, 33(4): 427-433.
[5]	吴海娜, 公卫江, 易光宇, 魏国柱. 纵向横向晶体场中自旋为1的量子伊辛二聚化链的基态[J]. 计算物理, 2015, 32(3): 374-378.
[6]	刘超, 石艺娜, 秦承森, 梁仙红. α铁冲击相变的离散元方法[J]. 计算物理, 2014, 31(1): 51-58.
[7]	胡安杰, 李隆键, 曾建邦. 作用力求取方法对多相伪势模型的影响[J]. 计算物理, 2012, 29(6): 828-836.
[8]	赵元松, 梁卫华, 陈文振. 围绕变壁温水平圆柱热源接触熔化的理论分析[J]. 计算物理, 2012, 29(2): 257-262.
[9]	黄勇, 宣益民, 李强. 水平圆管内磁性潜热型功能流体对流换热特性的数值模拟[J]. 计算物理, 2012, 29(1): 87-94.
[10]	赵元松, 梁卫华, 陈文振. 圆管内自由固体相变材料定热流接触熔化[J]. 计算物理, 2011, 28(4): 529-534.
[11]	周志坚, 胡燕飞. 金红石→萤石TiO₂相变的第一性原理计算[J]. 计算物理, 2011, 28(1): 111-118.
[12]	郝军华, 吴志强, 王铮, 金庆华, 李宝会, 丁大同. 加压下ZnO结构相变和电子结构的第一性原理计算[J]. 计算物理, 2010, 27(5): 759-764.
[13]	黄纯青, 邓绍军. 三维Ising模型的蒙特卡罗模拟[J]. 计算物理, 2009, 26(6): 937-941.
[14]	解立强. 高亏格膜泡形状的数值计算[J]. 计算物理, 2009, 26(5): 763-766.
[15]	杨培培, 刘红(小), 刘红. 双轴向列相液晶的Monte Carlo模拟[J]. 计算物理, 2009, 26(4): 597-602.