竖直圆管内超临界压力低温甲烷传热的数值研究

doi:10.3969/j.issn.1672-6952.2021.05.012

摘要/Abstract

摘要：

为了降低超临界压力甲烷管内对流换热的热阻，提高传热过程的稳定性，对竖直管内超临界甲烷的流动与传热过程进行了数值研究，讨论了热流密度、流动方向对流动换热特性的影响，研究了流场、温度场和湍动能变化对传热不稳定性的影响。结果表明，竖直管内超临界甲烷的传热存在不稳定传热现象，在高热流密度下，管内壁温度和平均温度具有不稳定性，在传热恶化区间震荡；在低热流密度下，局部对流传热系数具有不稳定性，在传热恶化区间震荡；在相同的热流密度下，向上流动时的传热不稳定性大于向下流动时的传热不稳定性，是由于向上流动的热影响大于向下流动的热影响，产生的类气膜是传热不稳定的主要因素。

关键词: 超临界压力, 甲烷, 传热, 不稳定性, 数值研究

Abstract:

In order to reduce the heat resistance and improve stability of convective heat transfer in liquefied natural gas(LNG) vaporizer, the flow and heat transfer processes of supercritical pressure LNG were analyzed by numerical method in a vertical tube. The influence of heat fluxes and flow directions on heat transfer characteristics was discussed, and the effects of the flow field , the temperature field , and the turbulent kinetic energy changes on the heat transfer instability were investigated.The results show that the heat transfer of supercritical methane in vertical tube is unstable,at high heat flux density, the inner wall temperature and the average temperature of the pipe are unstable and and shock in the heat transfer deterioration range; At the low heat flow density , the local convection heat transfer coefficient is unstable , and it shocks in the heat transfer deterioration range; Under the same heat flux density, the heat transfer instability of upward flow is greater than that of downward flow, because the thermal influence of upward flow is greater than that of downward flow. The gas?like film is the main reason of heat transfer instability.

Key words: Supercritical pressure, Methane, Heat transfer, Instability, Numerical studies

中图分类号:

TE088

李琳, 姜文全, 李婷婷, 杨帆, 宿诗雨, 石杰峰. 竖直圆管内超临界压力低温甲烷传热的数值研究[J]. 辽宁石油化工大学学报, 2021, 41(5): 66-71.

Lin Li, Wenquan Jiang, Tingting Li, Fan Yang, Shiyu Su, Jiefeng Shi. Numerical Study on Heat Transfer of Supercritical Pressure Cryogenic Methane in Vertical Pipes[J]. Journal of Liaoning Petrochemical University, 2021, 41(5): 66-71.

图/表 8

图1 竖直圆管内超临界压力低温甲烷传热三维物理模型

图2 基于REFPROP软件的甲烷在5 MPa下的物性参数

表1 各参数的含义

方程	φ	$ρ φ$	$Γ φ$	$S φ$
连续性方程	1	ρ	0	0
动量方程（x方向）	u	ρ	μ	$ρ f x - ∂ p ∂ x$
动量方程（y方向）	v	ρ	μ	$ρ f y - ∂ p ∂ y$
动量方程（z方向）	w	ρ	μ	$ρ f z - ∂ p ∂ z$
能量方程	T	$ρ c p$	λ	$S T$

表1 各参数的含义

方程	φ	$ρ φ$	$Γ φ$	$S φ$
连续性方程	1	ρ	0	0
动量方程（x方向）	u	ρ	μ	$ρ f x - ∂ p ∂ x$
动量方程（y方向）	v	ρ	μ	$ρ f y - ∂ p ∂ y$
动量方程（z方向）	w	ρ	μ	$ρ f z - ∂ p ∂ z$
能量方程	T	$ρ c p$	λ	$S T$

表2 网格独立性验证

序号	网格数	壁面温度/℃	误差/%
1	398 360	91.47	4.89
2	597 540	88.35	1.32
3	796 720	87.97	0.88
4	956 060	87.96	0.87
5	1 195 080	87.20	0

图3 沿程壁面温度的模拟数据及实验数据的对比

图4 不同热流密度下流体向上/向下流动时的温度及对流换热系数

图5 qw=125 kW/m2时竖直向上/向下流动时的温度场和流场云图（a）温度场（b）流场（z方向）

图6 qw=125 kW/m2的工况下向上/向下流动时沿程不同位置的径向流场和湍动能分布

参考文献 35

1	何依，邹斌，张丽，等.基于LNG冷能的双循环⁃卡琳娜冷电联供系统［J］.辽宁石油化工大学学报，2020，40（1）：43⁃51.
2	方旭，谢禹钧，常佩琛，等.基于断桥结构LNG罐式集装箱内接管传热分析［J］.石油化工高等学校学报，2019，32（2）：98⁃102.
3	韩昌亮，任婧杰，王焱庆，等. SCV耦合传热特性实验研究与数值模拟［J］.化工学报，2017，68（3）：854⁃863.
4	王博杰，匡以武，齐超，等.中间介质气化器中超临界LNG换热过程分析［J］.化工学报，2015，66（S2）：220⁃225.
5	李仲珍，郭少龙，陶文铨.超临界LNG管内流动与换热特性的研究［J］.工程热物理学报，2013，34（12）：2314⁃2317.
6	潘杰，吕涛，李冉，等.超临界压力下SuperORV传热管的传热特性［J］.工程热物理学报，2016，37（9）：1926⁃1934.
7	王亚洲，华益新，孟华.超临界压力下低温甲烷的湍流传热数值研究［J］.推进技术，2010，31（5）：606⁃611.
8	李玮哲，林文胜.超临界甲烷在印刷电路板换热器中加热过程模拟［J］.低温工程，2017（5）：60⁃64.
9	贾丹丹，赵忠超，张永，等.超临界LNG在印刷板式汽化器微细流道内的流动与换热性能数值研究［J］.船舶工程，2017，39（5）：35⁃40.
10	Gu H F，Li H Z，Wang H J，et al. Experimental investigation on convective heat transfer from horizontal miniature tube to methane at supercritical pressures［J］.Applied Thermal Engineering，2013，58（1⁃2）：490⁃498.
11	郭占魁，忻晔晨，李凌.超临界甲烷在竖直圆管内加热的数值模拟［J］.上海理工大学学报，2017，39（1）：30⁃34.
12	Huang D，Wu Z，Sunden N B，et al. A brief review on convection heat transfer of fluids at supercritical pressures in tubes and the recent process［J］.Applied Energy，2016，162：494⁃505.
13	白万金，徐肖肖，吴杨杨.低质量流速下超临界CO₂在管内冷却换热特性［J］.化工学报，2016，67（4）：1244⁃1250.
14	Zhang Q，Li H X，Kong X F，et al.Special heat transfer characteristics of supercritical CO₂ flowing in a vertically⁃upward tube with low mass flux［J］.International Journal of Heat and Mass Transfer，2018，122（1）：469⁃482.
15	Shiralkar B S，Griffith P.Deterioration in heattransfer to fluids at supercritical pressure and high heat fluxes［J］.Journal of Heat Transfer，1969，91（1）：27⁃36.
16	曲默丰，梁梓宇，赵云杰，等.低质量流速光管内超超临界水传热特性的实验与数值模拟［J］.西安交通大学学报，2018，52（7）：52⁃59.
17	Xu R N，Luo F，Jiang P X.Buoyancy effects on turbulent heat transfer of supercritical CO₂ in a vertical mini⁃tube based on continuous wall temperature measurements［J］.International Journal of Heat and Mass Transfer，2017，110（3）：576⁃586.
18	Yang Z，Cheng X，Zheng X H，et al.Numerical investigation on heat transfer of the supercritical fluid upward in vertical tube with constant wall temperature［J］.International Journal of Heat and Mass Transfer，2019，128（3）：875⁃884.
19	Kim D E，KimM H.Experimental study of the effects of flow acceleration and buoyancy on heat transfer in a supercritical fluid in a circular tube［J］.Nuclear Engineering and Design，2010，240（10）：3336⁃3349.
20	Jiang P X，Wang Z C，Xu R N.A modified buoyancy effect correction method on turbulent convection heat transfer of supercritical CO₂ pressure fluid bases on RANS models［J］.International Journal of Heat and Mass Transfer，2018，127（8）：257⁃267.
21	Zhao C R，Liu Q F，Zhang Z，et al.Investigation of buoyancy⁃enhanced heat transfer of supercritical in upward and downward tube flows［J］.Journal of Supercritical Fluids，2018，138：154⁃166.
22	Tian R，ZhangY，Ma Y Z.Experimental study of buoyancy effect and its criteria for heat transfer of supercritical R134a in horizontal tubes［J］.International Journal of Heat and Mass Transfer，2018，127（8）：555⁃567.
23	刘生晖，黄彦平，刘光旭，等.浮升力因子和流动加速因子改进及其在超临界流体混合对流传热中的应用［J］.中国科学：技术科学，2017，47（2）：176⁃189.
24	Zhu B G， Xu J L， Wu X M， et al. Supercritical “boiling” number， a new parameter to distinguish two regimes of carbon dioxide heat transfer in tubes［J］.International Journal of Thermal Sciences，2019，136（20）：254⁃266.
25	Song J H，Kim H Y，Bae Y Y.Heat transfer characteristics of supercritical fluid flow in a vertical tube［J］.Journal of Supercritical Fluids，2008，44（10）：164⁃171.
26	浦航，李素芬，东明，等.超临界压力RP⁃3起始加热段传热恶化数值研究［J］.工程热物理学报，2017，38（10）：2242⁃2248.
27	贾洲侠，徐国强，闻洁，等.超临界压力RP⁃3在竖直细圆管内混合对流研究［J］.北京航空航天大学学报，2016，42（1）：152⁃157.
28	王淑香，张伟，牛志愿，等.超临界压力下CO₂在螺旋管内的混合对流换热［J］.化工学报，2013，64（11）：3917⁃3926.
29	Yildiz S，Groeneveld D C.Diameter effect on supercritical heat transfer［J］.International Communications in Heat and Mass Transfer，2014，54（2）：27⁃32.
30	范辰浩，王海军，宋子琛，等.小管径超临界流体传热恶化特性研究［J］.工程热物理学报，2018，39（9）：2032⁃2039.
31	李志刚，淮秀兰，陶毓伽，等.变物性对微通道内流动与传热的影响［J］.工程热物理学报，2007，28（4）：640⁃643.
32	代宝民，李敏霞，吕佳桐，等.超临界CO₂/R41小通道内的换热特性［J］.化工学报，2015，66（3）：924⁃931.
33	王彦红，李素芬，赵星海.超临界压力航空煤油不稳定流动实验［J］.航空动力学报，2018，33（12）：2838⁃2844.
34	宇波.数值传热学实训［M］.北京：科学出版社，2018.
35	李志辉，姜培学.高雷诺数条件下超临界压力CO₂在竖直圆管内换热特性的实验研究［J］.核动力工程，2008，29（6）：41⁃45.

[1]	葛韶华, 艾贵琳, 杨帆, 陈昊楠, 屈思雨, 姜文全. 超临界CO₂及其混合物在竖直圆管内的传热特性数值研究[J]. 辽宁石油化工大学学报, 2023, 43(3): 75-80.
[2]	翟郑佳，李国龙，朱恒宣，杨历，王进. 磁场对Fe₃O₄⁃水纳米流体传热特性的影响[J]. 辽宁石油化工大学学报, 2020, 40(2): 22-29.
[3]	彭雪松，马贵阳. LAD⁃40对甲烷水合物生成影响研究[J]. 辽宁石油化工大学学报, 2020, 40(2): 53-57.
[4]	李宏达，张连红，蒋林时，肖宏. 超声波合成三（4⁃乙酰氨基苯基）甲烷[J]. 辽宁石油化工大学学报, 2019, 39(5): 26-29.
[5]	卜英南，陈祥彬，张莹莹，高磊. 折流板结构对换热器壳程流动和传热性能的影响[J]. 辽宁石油化工大学学报, 2018, 38(05): 72-76.
[6]	赵波,王阔,隋成玉. 三叶草型催化剂体相温度分布数值模拟[J]. 辽宁石油化工大学学报, 2017, 37(2): 10-16.
[7]	王云,马贵阳,徐子健. 热油管道停输安全出站温度数值研究[J]. 辽宁石油化工大学学报, 2015, 35(1): 29-32.
[8]	官学源,刘德俊,范欣,毕潆,王秋莎,赵志强. 含硫天然气管道泄漏危险区域的数值模拟[J]. 辽宁石油化工大学学报, 2013, 33(1): 39-42.
[9]	王玉福，申龙涉. 集肤效应电伴热管道再启动三维非稳态传热过程数值模拟[J]. 辽宁石油化工大学学报, 2011, 31(4): 42-44.
[10]	许丹, 申龙涉, 杜明俊, 陈忠华. 埋地热油管道停输三维非稳态传热过程的数值模拟[J]. 辽宁石油化工大学学报, 2010, 30(4): 47-50.
[11]	李小玲,马贵阳. 土壤源热泵单U 型埋管换热器短期运行换热分析[J]. 辽宁石油化工大学学报, 2010, 30(1): 51-54.
[12]	仇涤凡，李帆. 柴油机配气机构动力学特性的数值研究[J]. 辽宁石油化工大学学报, 2008, 28(2): 44-47.
[13]	李　力, 高文艺, 等. 苯和1 - 十二烯在甲烷磺酸催化剂中的液相烷基化[J]. 辽宁石油化工大学学报, 2007, 27(1): 17-20.
[14]	杨　军, 陈保东,孙成家. 螺旋与弓形折流板换热器性能对比及螺旋角优化[J]. 辽宁石油化工大学学报, 2005, 25(2): 59-62.
[15]	苏婷婷, 张连惠, 宫　红, 王　锐. 甲烷磺酸锌的热分析及脱水动力学[J]. 辽宁石油化工大学学报, 2005, 25(1): 20-22.