水平管内R1234yf的流动沸腾换热特性

化
工进展
Chemical Industry and Engineering Progress
2022年第41卷第7期
水平管内R1234yf 的流动沸腾换热特性
冯龙龙1，钟珂1，张羽森1，商庆春2，贾洪伟1
（1东华大学环境科学与工程学院，上海201620；2山东电力建设第一工程有限公司，山东济南250200）摘要：通过实验研究了环境友好型制冷剂R1234yf 在内径为0.5mm 的水平圆形微通道内的流动沸腾换热特性，测量了不同工况下R1234yf 的沸腾换热系数（HTC ），并与传统制冷剂R134a 进行了对比，分析了质量流速、热流密度和干度对换热系数变化规律的影响。

实验条件为：饱和温度(17±1)℃，质量流速1000~2500kg/(m 2·s)，热流密度25~143kW/m 2。

实验结果表明：R1234yf 的换热系数随着热流密度的增大而显著增大，而质量流速和干度的影响较小，核态沸腾为其主导换热机制。

对比R1234yf 和R134a 在相同工况下的换热特性，发现两种工质的平均换热系数差别较小，并均随着热流密度增大而逐渐增加，但是R1234yf 发生干涸（Dryout ）时的热流密度小于R134a 。

将实验数据与已有文献中的核沸腾主导的经验关联式的预测结果进行了对比，得到了较好的吻合。

关键词：制冷剂R1234yf ；微通道；流动沸腾；关联式中图分类号：TK124
文献标志码：A
文章编号：1000-6613（2022）07-3502-08
Flow boiling heat transfer characteristics of R1234yf in horizontal
microchannel
FENG Longlong 1，ZHONG Ke 1，ZHANG Yusen 1，SHANG Qingchun 2，JIA Hongwei 1
(1College of Environmental Science and Engineering,Donghua University,Shanghai 201620,China;2ESPCO1Electric
Power Construction Co.,Ltd.,Jinan 250200,Shandong,China)
Abstract:The flow boiling heat transfer characteristics of environmentally friendly refrigerant R1234yf in a 0.5mm horizontal circular microchannel were studied experimentally.The heat transfer coefficients (HTCs)of R1234yf were measured and compared with that of R134a,and the effects of mass flux,heat flux and vapor quality on HTC were analyzed.The saturation temperature was (17±1)℃,and the mass fluxes vary from 1000kg/(m 2·s)to 2500kg/(m 2·s)with heat fluxes ranging from 25kW/m 2to 143kW/m 2.The experimental results showed that the HTC of R1234yf in 0.5mm microchannel increases with the increase of heat flux,while the mass flux and vapor quality showed a weak influence on it.The trend indicated that nucleate boiling was the dominant mechanism for flow boiling heat transfer.In addition,the heat transfer performance of R1234yf and R134a were compared under the same working conditions.The HTCs of R1234yf and R134a were almost identical and both increased with the increase of heat flux,but the heat flux for the occurrence of dryout of R1234yf was smaller than that of R134a.Finally,the experimental data for the two refrigerants were compared with two nucleate boiling-dominated correlations from literature and good agreements were obtained.
研究开发
DOI ：10.16085/j.issn.1000-6613.2021-1823
收稿日期：2021-08-25；修改稿日期：2021-11-01。

基金项目：国家自然科学基金（52006030）；上海市青年科技英才“扬帆计划”（18YF1400700）；中国博士后基金（2018M641891）。

第一作者：冯龙龙（1992—），男，博士研究生，研究方向为微通道换热、低GWP 制冷剂换热。

E-mail ：。

通信作者：贾洪伟，博士，讲师，研究方向为多相流数值计算、微尺度热质交换。

E-mail ：。

引用本文：冯龙龙,钟珂,张羽森,等.水平管内R1234yf 的流动沸腾换热特性[J].化工进展,2022,41(7):3502-3509.
Citation ：FENG Longlong,ZHONG Ke,ZHANG Yusen,et al.Flow boiling heat transfer characteristics of R1234yf in horizontal microchannel[J].Chemical Industry and Engineering Progress,2022,41(7):3502-3509.
·
·3502
2022年7月冯龙龙等：水平管内R1234yf的流动沸腾换热特性Keywords:refrigerant R1234yf;microchannel;flow boiling;correlations
《基加利修正案》[1]已于2021年9月15日起对中国生效，具有高全球变暖潜能值（GWP）的
HFCs类制冷剂将会被加速淘汰。

低温室效应工质R1234yf[2]具有与R134a相近的热物性，如表1所示，被认为是一种极具潜力可替代R134a的环保制冷剂[3-5]。

国内外学者也针对R1234yf和R134a的换热性能对比及替代应用展开了实验研究。

例如，Yang等[6]和Longo等[7]通过实验对比了R1234yf与R134a在4mm管内流动沸腾的换热性能。

Yang等[6]的结果表明R1234yf的换热系数比R134a低约17%；而Longo等[7]发现x<0.3时两种制冷剂的换热系数非常接近，而x>0.3时R1234yf的换热系数略低于R134a。

Lu等[8]的实验结果则发现两种制冷剂在3.9mm铜管内的换热系数差别很小。

在应用方面，王锐等[9]的实验研究发现R1234yf可替代R134a和R410A，作为喷雾冷却辅助激光手术治疗的替代制冷剂。

孟照峰等[10]的实验结果表明混合制冷剂R1234yf/R134a（质量分数89%/11%）可以在汽车空调中直接替代R134a。

杨梦等[11]的实验表明R1234yf与R134a的混合制冷剂R513A（质量分数56%/44%）可直接用于替代家用冰箱中的R134a，并且其性能要优于R134a。

近年来，微细/微小通道因其紧凑的结构、较少的工质充注量、较低的材料成本和优异的换热性能等优点[12-14]，在R1234yf的研究中也逐渐受到关注。

De Del Col等[15]对R1234yf在0.96mm管内的流动沸腾换热进行了实验，结果表明，换热系数与热流密度强烈相关，而质量流速几乎没有影响。

Saitoh等[16]研究了R1234yf在2mm管内的流动换热特性，发现低干度时热流密度对两相换热的影响较大，而高干度时质量流速的影响更显著。

Anwar 等[17]通过实验研究了R1234yf在1.6mm微通道内竖直向上流动的换热特征，结果表明换热系数由热流密度控制，而质量流速和干度几乎不影响换热系数变化趋势。

Ribatski等[18]研究了低GWP制冷剂[R1234yf、R1234ze(E)、R600a]和R134a在1.1mm 管内流动的沸腾换热，并基于实验测得的3409个数据点提出了一个新的换热系数预测关联式。

总结前人的研究结果，发现对R1234yf在管径更小（<1mm）的微细通道内流动沸腾换热特性研究以及该尺度条件下与传统制冷剂的换热性能对比均较为缺乏，因此本文对制冷剂R1234yf在0.5mm 微通道内的流动沸腾换热特性进行了实验研究。

实验数据可以扩充R1234yf在微细通道内流动沸腾换热的数据库，主要结论也可为R1234yf的高效应用和推广提供研究基础和实验依据。

1实验装置
1.1实验系统
图1为实验系统示意图，主要由储液罐、齿轮泵、质量流量计、预热段、冷凝器以及实验段组成。

储液罐用于储存和回收系统中的制冷剂。

质量流量（m）通过控制齿轮泵的转速和旁通阀门的开度进行调节，并由质量流量计实时测量。

预热段为一段外壁面缠绕加热丝的铜管，通过调节单相调压器可控制预热段的加热量，并保证制冷剂在实验段入口处保持一定过冷度（<3K）。

在质量流量计和预热段的前端以及实验段的后端分别装有一个冷凝器（与低温恒温槽连接的板式换热器）用来冷却系统中的制冷剂。

测量数据通过安捷伦34972A记录并存储在电脑中。

为了减少系统漏热，所有设备表面覆盖有一层厚度约为9mm、热导率约为0.034 W/(m·K)的保温棉。

1.2实验段介绍
图2为实验段示意图。

实验段为内径0.5mm、外径1.5mm、长度300mm的不锈钢水平光滑圆管，通过光学3D表面轮廓仪（SuperView W1）测得管道内表面平均粗糙度Ra约为1.065μm。

外管壁上等间距布置11只T形热电偶用于测量外壁面温度，热电偶用聚亚酰胺胶带固定。

管道两端焊接有四通
表1R1234yf和R134a在饱和温度为17℃时的部分热物性参数
制冷剂R1234yf
R134a ODP
GWP
4
1430
饱和压力p sat
/MPa
0.542
0.521
液体密度ρl
/kg·m-3
1120.3
1236.2
气体密度ρv
/kg·m-3
30.01
25.305
汽化潜热i lv
/kJ·kg-1
151.56
184.89
定压比热容c
p，l
/kJ·kg-1·K-1
1.3563
1.3939
黏度μl
/μPa·s
169.69
215.24
热导率λl
/mW·m-2·K-1
66.011
84.579
表面张力σ
/mN·m-1
7.182
9.093
注：ODP为消耗臭氧潜能值；GWP为全球变暖潜能值。

··3503
化工进展，2022，41（7）
结构，并在其垂直方向上分别装有一个铠装热电偶和一个压力传感器，分别用于测量制冷剂在实验段进出口处的温度和压力。

实验段两端直接连接一个直流电源（5V ，200A ），根据焦耳效应（Joule effect ）对实验管段进行加热，并通过调节输出电
压改变加热量。

2数据分析
2.1
数据处理及单相实验验证
进行两相沸腾实验之前，需要通过单相流动换热实验测定系统的加热效率（ε）[19-20]。

由于R1234yf 和R134a 在常压下的沸点较低，为了更接近R1234yf 和R134a 的沸腾实验中实验段外壁面的
温度条件，以便更准确地估算有效加热量，单相实验采用了常压沸点较高的R245fa 作为流动工质。

为了避免工质在实验段内发生相变，在单相测试过程中应保持实验段温度始终低于制冷剂在对应压力下的饱和温度。

由热平衡可知，直流电源输出功率（Q ）等于有效加热量（Q eff ）与热损失（Q loss ）之和，如式(1)。

Q =Q eff +Q loss
(1)
式中，Q=U ×I ，为直流电源的输出功率，W ；Q eff 是工质流过实验段吸收的热量，单相换热仅为显热变化，可写为式(2)。

Q eff =mc p ,l (T l,out -T l,in )
(2)
式中，c p ,l 是定压比热容，J/(kg·K)；T l,in 和T l,out
分别是实验段进出口过冷制冷剂流体的温度。

ε定义为Q 与Q eff 的比值，如式(3)。

ε=
Q eff Q
(3)
通过上述步骤可获得加热效率，结果如图3所示。

可以看到本实验台的平均加热效率约为78.91%。

此外，将获得的单相流动换热的Nusselt 数与Gnielinski [21]关联式(4)的预测结果进行比较，如图4所示。

可以看到实验测得的单相换热Nusselt 数与预测关联式吻合较好。

并且从图4(b)可以看出，97.9%的单相流动数据点都落在±15%的误差范围内，这也表明本文实验装置具有较好的可靠性。

Nu Gnielinski =
(f /8)(Re -1000)
1+12.7f /8(Pr 2/3-1)(4)
根据实验获得的加热效率可以估算实验段加热的有效热流密度（q eff ），如式(5)。

q eff =
εQ
πD in L
(5)
同时，沿程各热电偶测点对应的内壁面温度（T wall,in,z ）可通过Fourier 定律估算，如式(6)。

T wall,in,z =T wall,out,z
+εQ 4πλL éëêêζ2()
1-ln ζ2
-1ζ2-1ùû
ú
ú,ζ=D out D in (6)式中，T wall,out,z 为外壁面温度；D out 为实验段管外径；λ为不锈钢热导率。

于是可获得局部换热系数（h z ），如式(7)。

h z =
q eff
T wall,in,z -T f,z
(7)
质量流速（G ）为质量流量与通道截面积的比值，如式(8)。

G =
4m πD in 2
(8)
各热电偶测点处对应的干度（x z ）可通过热平衡计算，如式(9)。

x z =
i l,in +
εQ m ⋅z
L -i sat
i l,v
(9)
式中，i L,in 是进口液体制冷剂在进口温度T in 和
图1实验装置示意图
1—齿轮泵；2—旁通阀；3—冷凝器；4—质量流量计；5
—预热段；6—实验段；7—数据采集仪&电脑；8—储液罐；9—干燥过滤器
图2实验段及热电偶贴附示意图
·
·3504
2022年7月冯龙龙等：水平管内R1234yf 的流动沸腾换热特性
进口压力p in 下的焓值；i sat 是饱和段起始点处的焓值；i l,v 是制冷剂气化潜热。

上述数据处理过程均由MATLAB 编程实现，
计算过程中所需要的制冷剂的热物性参数通过调用物性软件REFPROP 9.1动态数据库直接获得。

2.2
不确定性分析
实验所用的热电偶都经过标准热电偶校准，精度为±0.2K 。

压力传感器的满量程精度为±0.25%；质量流量计的精度为0.1%；数字电流表和电压表分别用来显示直流电源的输出电流及电压，不确定度均为±0.5%。

数据处理过程中涉及的计算参数的不确定度可通过Moffat [22]的方法计算得到，如
式(10)。

δR =
(10)
式中，δR 表示参数R 的绝对不确定度；X i 为自
变量；δX i 为其绝对不确定度。

各参数的相对不确定度分别为δX i /X i 和δR/R 。

不确定度计算结果如表2所示。

3结果与讨论
3.1
热流密度对换热的影响
图5是R1234yf 在不同热流密度下局部换热系数（h local ）随着沿程干度（x ）的变化。

由图5可知，在同一干度条件下，h local 均随着热流密度的增大而增大。

这是因为随着热流密度的增加，管道内壁面上更多的核化点被激活，气泡数量增多，同时流体扰动也更加剧烈，换热能力得到显著增强[23-24]。

热流密度对两相流动换热具有显著的促进作用，核态沸腾为沸腾换热的主导机制[25]。

x <0.7时，h local 的变化趋势基本保持一致，干度对流动换热的影响较小。

热流密度较低时，随着干度的增加
换热系数略微下降。

在低干度区，干度的增大使更多的核化点被激活，管内流型多为泡状流；在中高干度区，随着小气泡的生长和大气泡的合并，管内流型过渡到拉长气泡流。

热流密度较高时，换热系数整体呈略微上升的趋势。

较高的热流密度使管内成核点在干度较低时就已经被高度激活，换热系数处在较高水平；随着干度的增加，管内流型更早的向拉长气泡流和环状流转变。

x>0.7时，实验段后端的h local
出现明显的下降，表明管内制冷剂出现干涸。

此时管内流型多为弥散流[26]，即制冷剂气体与
管道内壁面直接接触，制冷剂液滴以弥散状态分散在管道中心位置，气态制冷剂的导热热阻成为阻碍换热的主要原因。

由于热量不能被迅速带走而不断积聚使得壁面过热度迅速增加，导致换热系数显著
图3
实验段加热效率
图4
单相Nusselt 数与Gnielinski [21]关联式预测值的对比
表2
实验参数的不确定度
测量值
参数名m /kg·s D ，L/m U /V ，I /A T （T-type ）/℃p /Pa
不确定度±0.1%±2%
±0.5%±1.17%±0.25%
计算参数参数名G /kg·m -2·s -1q /kW·m -2h /kW·m -2·K -1
x
Q /W
不确定度±2.83%±2.12%±6.76%±0.87%±0.71%
·
·3505
化工进展，2022，41（7）
下降。

此外，当热流密度增大时，管内流动的干涸
点位置会提前。

3.2质量流速对换热的影响
图6为R1234yf 在G =1000~2500kg/(m 2·s)范围内局部换热系数随着沿程干度的变化。

当q =50kW/m 2时，数据点的整体平均换热系数（h mean ）约为6.52kW/(m 2·K)；当q 增加到58kW/m 2时，h mean 约为7.45kW/(m 2·K)。

相同热流密度下的局部换热系数基本在同一水平，质量流速对换热的影响较小。

当q =58kW/m 2、G =1000kg/(m 2·s)时，实验段后端已经出现换热恶化，而其他质量流速工况还未出现，因此高质量流速能延缓干涸的发生[图6(b)]。

图7为不同质量流速下R1234yf 的平均换热系数（h ave ）随热流密度的变化。

由图7可知，在干涸发生之前，沿程平均换热系数均随热流密度的增加而增加，同一热流密度下，不同质量流速的平均换热系数差别很小。

R1234yf 的两相流动换热受热流密度的影响较大，而质量流速的影响作用较小。

在较高热流密度时，由于靠近实验段出口处出现换热恶化，从而导致平均换热系数随热流密度变化的斜
率变小，并且质量流速越大，转折点对应的热流密度越大。

这是因为达到干涸所需的热流密度与质量流速成正比[27]，在高质量流速下，壁面热量可以被迅速带走，使壁面温度在较高热流密度时仍保持较低水平。

3.3R1234yf 与R134a 的换热特性对比
图8对比了R1234yf 和R134a 的局部换热系数随干度的变化。

可以看出两种制冷剂具有相似的换热趋势。

但是由于R1234yf 的汽化潜热比R134a 低约18%
，通过相变换热带走相同热量所需的制冷剂
图6
质量流速对局部换热系数的影响
图5
热流密度对局部换热系数的影响
图7沿程平均换热系数随热流密度的变化
·
·3506
2022年7月冯龙龙等：水平管内R1234yf 的流动沸腾换热特性
蒸发量更多，所以在同一热流密度下，相同测点处R1234yf 的干度略大于R134a 。

此外，由于R1234yf 的表面张力比R134a 低约21%，在高热流密度下，R1234yf 的近壁面液膜更易破碎，也更易出现换热恶化。

当q =58.68kW/m 2时，R1234yf 在测点P-10和P-11处出现换热系数急剧下降，当q 进一步增加至66.96kW/m 2时，在测点P-9处就已出现换热恶化，
而R134a 在热流密度达到66.96kW/m 2时仍未出现换热恶化现象。

图9对比了R1234yf 和R134a 的平均换热系数（h ave ）随热流密度的变化。

在干涸发生前，R1234yf 和R134a 的h ave 分布基本重合，相同热流密度下的h ave 差别很小。

而随着热流密度进一步升高，
两种工质的管内流动均会出现换热恶化，平均换热系数的斜率变小，并且转折点的热流密度均随质量流速的增大而增大，但R1234yf 干涸出现时的热流密度显著小于R134a 。

3.4
实验数据与预测模型对比
实验所得的R1234yf 与R134a 的平均换热系数与Cooper [28]
和Tran 等[29]
提出的核沸腾主导的换热系数预测模型进行了对比，并采用实验数据与关联式预测值的平均绝对误差（MAE ）来评估预测模型的准确性，如式(11)。

MAE =1N ∑i =1N ||X pre,i -X exp,i X exp,i
×100%(11)
图10(a)是R1234yf 与R134a 的平均换热系数与
Cooper [28]模型的对比结果。

两种制冷剂的实验数据
与关联式预测值的MAE 分别为10.32%和9.96%。

R1234yf 有98.25%的实验值与预测值的偏差可以维持在±30%以内；R134a 有94.74%的实验值与预测值的偏差可以维持在±30%以内。

R1234yf 与R134a 的平均换热系数与Tran 等[29]的模型对比结果
如图10(b)所示，两种制冷剂的实验数据与关联式
预测值的MAE 分别为5.3%和25.97%。

R134a 有70.18%的实验值与预测值的偏差可以维持在±30%以内，而R1234yf 的全部实验值与预测值的偏差都
在±30%以内。

总体来说，两种模型均对R1234yf
在微通道内沸腾换热具有较好的预测。

4结论
（1）在0.5mm 微通道内，R1234yf 的沸腾换热
系数随热流密度的增大而增大，热流密度对两相换热具有显著的促进作用，而质量流速和干度的影响较小，核态沸腾为其主导换热机制。

（2）微尺度条件下，R1234yf 和R134a 呈现相
似的换热趋势。

相同工况下，两种工质的换热系数
图9R1234yf 与R134a
的平均换热系数对比
图8R1234yf 与R134a 的局部换热系数对比
·
·3507
化工进展，2022，41（7）
几乎一致，但R1234yf 出现干涸时的热流密度小于R134a 。

（3）实验测得的平均换热系数与Cooper 和Tran
等的模型预测结果进行对比，发现两种模型都能对R1234yf 进行较为准确的预测，但对于R134a ，
Cooper 模型的预测结果优于Tran 等模型。

符号说明
G ——质量流速，kg/(m 2·s)
GWP ——全球变暖潜能值（global worming potential ）
HTC ——换热系数（heat transfer coefficient ），kW/(m 2
·K)h ——换热系数，kW/(m 2·K)i ——焓，J/kg
MAE ——平均绝对误差（mean absolute error ），%m ——质量流量，kg/s
ODP ——消耗臭氧潜能值（ozone depletion potential ）p ——压力，Pa q ——热流密度，kW/m 2T ——温度，℃x ——干度ε——加热效率，%下角标
ave ——平均值exp ——实验值l ——液态
local ——局部值pre ——预测值sat ——饱和值v ——气态参考文献
[1]UNEP.[Ｃ]//38th Meeting of the Open-ended Working Group (OEWG 38)of the Parties to the Montreal Protocol on Substances that Deplete the Ozone Layer and 28th Meeting of the Parties to the Montreal Protocol (MOP 28),Kigali,2016.
[2]NIELSEN O J,JAVADI M S,SULBAEK ANDERSEN M P,et al.Atmospheric chemistry of CF 3CF ==CH 2:kinetics and mechanisms of gas-phase reactions with Cl atoms,OH radicals,and O 3[J].Chemical Physics Letters,2007,439(1/2/3):18-22.
[3]MCLINDEN M O,KAZAKOV A F,STEVEN BROWN J,et al.A thermodynamic analysis of refrigerants:possibilities and tradeoffs for Low-GWP refrigerants[J].International Journal of Refrigeration,2014,38:80-92.[4]QI Z G.Performance improvement potentials of R1234yf mobile air conditioning system[J].International Journal of Refrigeration,2015,58:35-40.
[5]GARCIA PABON J,KHOSRAVI A,NUNES R,et al.Experimental investigation of pressure drop during two-phase flow of R1234yf in smooth horizontal tubes with internal diameters of 3.2mm to 8.0mm [J].International Journal of Refrigeration,2019,104:426-436.
[6]YANG C Y,NALBANDIAN H,LIN F C.Flow boiling heat transfer and pressure drop of refrigerants HFO-1234yf and HFC-134a in small
circular tube[J].International Journal of Heat and Mass Transfer,2018,121:726-735.
[7]LONGO G A,MANCIN S,RIGHETTI G,et al.R1234yf and R1234ze (E)as environmentally friendly replacements of R134a:Assessing flow boiling on an experimental basis[J].International Journal of Refrigeration,2019,108:336-346.
[8]LU M C,TONG J R,WANG C C.Investigation of the two-phase convective boiling of HFO-1234yf in a 3.9mm diameter tube[J].
International Journal of Heat and Mass Transfer,2013,65:545-551.[9]王锐,陈斌,王嘉丰,等.R1234yf 瞬态喷雾冷却及过热度影响的实验研究[J].化工学报,2018,69(2):595-601.WANG Rui,CHEN Bin,WANG Jiafeng,et al.Experimental research of R1234yf transient spray cooling and influence of cryogen superheat degree[J].CIESC Journal,2018,69(2):595-601.
[10]孟照峰,张华,秦延斌,等.R1234yf/R134a 混合物在汽车空调中替
代R134a 的实验研究[J].化工学报,2018,69(6):2396-2403.
MENG Zhaofeng,ZHANG Hua,QIN Yanbin,et al.Experimental study on R1234yf/R134a mixture as alternative to R134a in automobile air conditioner[J].CIESC Journal,2018,69(6):2396-2403.
[11]杨梦,张华,秦延斌,等.混合制冷剂R134a/R1234yf(R513A)与R134a
热力学性能对比及实验[J].化工进展,2019,38(3):1182-1189.YANG Meng,ZHANG Hua,QIN Yanbin,et al.Thermodynamic performance comparison and experimental study of mixed refrigerant R134a/R1234yf(R513A)and R134a[J].Chemical Industry and Engineering Progress,2019,38(3):1182-1189.
[12]ANWAR Z,PALM B,KHODABANDEH R.Flow boiling heat transfer
and dryout characteristics of R152a in a vertical mini-channel[J].Experimental Thermal and Fluid Science,2014,53:
207-217.图10换热系数实验结果与模型计算结果的对比
·
·3508
2022年7月冯龙龙等：水平管内R1234yf的流动沸腾换热特性
[13]LI W,WU Z.A general criterion for evaporative heat transfer in micro/ mini-channels[J].International Journal of Heat and Mass Transfer, 2010,53(9/10):1967-1976.
[14]YANG Q C,SHU B F,WANG J,et al.Experimental investigation on flow boiling heat transfer and flow patterns in a single micro-channel with large mass velocity[J].Experimental Thermal and Fluid Science, 2018,91:283-291.
[15]DE DEL COL D,BORTOLIN S,TORRESIN D,et al.Flow boiling of R1234yf in a1mm diameter channel[J].International Journal of Refrigeration,2013,36(2):353-362.
[16]SAITOH S,DANG C B,NAKAMURA Y,et al.Boiling heat transfer of HFO-1234yf flowing in a smooth small-diameter horizontal tube[J]. International Journal of Refrigeration,2011,34(8):1846-1853. [17]ANWAR Z,PALM B,KHODABANDEH R.Flow boiling heat transfer, pressure drop and dryout characteristics of R1234yf:experimental results and predictions[J].Experimental Thermal and Fluid Science, 2015,66:137-149.
[18]SEMPÉRTEGUI-TAPIA D F,RIBATSKI G.Flow boiling heat transfer of R134a and low GWP refrigerants in a horizontal micro-scale channel[J].International Journal of Heat and Mass Transfer,2017, 108:2417-2432.
[19]ONG C L,THOME J R.Flow boiling heat transfer of R134a,R236fa and R245fa in a horizontal1.030mm circular channel[J].Experimental Thermal and Fluid Science,2009,33(4):651-663.
[20]党超,贾力,黄浅.矩形微槽道内R134a流动沸腾换热特性的实验
研究[J].工程热物理学报,2017,38(6):1327-1332. DANG Chao,JIA Li,HUANG Qian.Experimental study on flow boiling heat transfer characteristics of R134a in a rectangular micro-channel[J].Journal of Engineering Thermophysics,2017,38(6):1327-1332.
[21]GNIELINSKI V.New equations for heat and mass transfer in turbulent pipe and channel flow[J].International Journal of Chemical
Engineering,1976,16:359-368.
[22]MOFFAT R J.Describing the uncertainties in experimental results[J]. Experimental Thermal and Fluid Science,1988,1(1):3-17. [23]杞卓玲,贾力,党超.微细通道内非共沸工质流动沸腾换热实验研
究[J].工程热物理学报,2018,39(9):2005-2011.
QI Zhuoling,JIA Li,DANG Chao.Experimental investigation on flow boiling heat transfer characteristics of zeotropic mixture in microchannels[J].Journal of Engineering Thermophysics,2018,39(9): 2005-2011.
[24]刘梦悦,柳建华,张良.R1234yf在小管径管内流动沸腾换热特性研
究[J].热能动力工程,2020,35(9):65-70.
LIU Mengyue,LIU Jianhua,ZHANG Liang.Study on flow boiling heat transfer characteristics of R1234yf in small-diameter tubes[J].Journal of Engineering for Thermal Energy and Power,2020,35(9):65-70.
[25]KAEW-ON J,SAKAMATAPAN K,WONGWISES S.Flow boiling heat transfer of R134a in the multiport minichannel heat exchangers [J].Experimental Thermal and Fluid Science,2011,35(2):364-374.
[26]TORFEH S,KOUHIKAMALI R.Numerical investigation of mist flow regime in a vertical tube[J].International Journal of Thermal Sciences, 2015,95:1-8.
[27]MARTI´N-CALLIZO C,ALI R,PALM B.Dryout incipience and critical heat flux in saturated flow boiling of refrigerants in a vertical uniformly heated microchannel[C]//Proceedings of ASME20086th International Conference on Nanochannels,Microchannels,and Minichannels,June23-25,2008,Darmstadt,Germany.2009: 705-712.
[28]COOPER M G.Saturation nucleate pool boiling—A simple correlation
[C]//First U.K.National Conference on Heat Transfer.Amsterdam: Elsevier,1984:785-793.
[29]TRAN T N,WAMBSGANSS M W,FRANCE D M.Small circular-and rectangular-channel boiling with two refrigerants[J].International Journal of Multiphase Flow,1996,22(3):485-498.
·
·3509。