电场对微细通道内R141b制冷剂流动沸腾压降的影响

罗小平，彭子哲，刘　倩，郭　峰，章金鑫

电场对微细通道内R141b制冷剂流动沸腾压降的影响

罗小平，彭子哲，刘倩，郭峰，章金鑫

（华南理工大学机械与汽车工程学院，广州 510640）

在农业工程领域，微细通道散热技术在农产品培育系统、农业机械、农产品干燥系统中有着广泛的应用。通过施加电场可强化微细通道换热系统的传热效率，为探究电场对微细通道内制冷剂流动沸腾阻力的影响，该文采用了2种电极布置方式（针状和线状），以制冷剂R141b为试验工质，在系统压力为140 kPa，工质入口温度32.5 ℃、质量流率277.35～531.75 kg/(m2·s)、热流密度7.50～21.49 kW/m2、电压0～850 V工况下，在截面尺寸为2 mm×2 mm的矩形微细通道内进行流动沸腾试验，探究直流电场对微细通道内R141b流动沸腾压降特性影响。研究结果表明：在本文试验工况下，电场会增大微细通道内的摩擦压降，针状与线状电极电场作用下的微细通道内摩擦压降分量在总压降中所占比例均比无电极作用下的更大；电场作用下单位长度两相摩擦压降随电压、热流密度的增大而增大，针状电极与线状电极电场作用下平均单位长度两相摩擦压降分别比无电极作用下增加0.7%～15.4%和1.3%～18.7%；电压为0～250 V时，针状电极对压降的影响效果大于线状电极，电压大于400 V后，线状电极对压降的影响效果更为显著。通过COMSOL软件对6 mm长微细通道内2种电场的分布进行了模拟，模拟结果表明相同电压作用下，针状电极产生的电场强度最大值超过线性电极，但线状电极的电场有效作用范围超过针状电极。该文研究结果可为通过施加电场提高微细通道换热器的性能实现微细通道高效节能提供新思路。

制冷剂；电场；电极；微细通道；流动沸腾；压降

0 引 言

微通道散热在集成电路、农业机械、航空等领域有着越来越广泛的应用[1]，在农业工程中，LED光源系统对农作物的生长具有重要意义[2]，然而LED光源的散热却存在一定的问题，将微细通道散热技术应用于LED光源值得深入研究，此外微细通道在太阳能电池的散热、农产品干燥系统的设计等方面也有重要的前景[3-4]。

近年来，国内外学者研究发现电场可以明显地改善散热器传热性能[5-9]，黄烜等[10]研究了电场分布对R123沸腾传热的影响，表明电场的强化传热是电场强度和电场均匀性共同作用的结果。Zhang等[11]研究了竖直微细通道内针状电极和线状电极电场对流动沸腾传热的强化作用，试验结果表明施加电压≤400 V时，针状电极电场强化效果好；电压≥550 V后，线状电极强化效果好。利用电场来强化传热具有简单、快速、低功耗等优点[12]，流动沸腾压降是设计微细通道两相冷却系统时要考虑的关键参数[13]，所以在研究电场的流动沸腾强化传热时，压降是非常重要的性能指标。Bryan等[14]研究了常规直径光滑管和强化管中电场对R134a与R404a的流动沸腾压降的影响，研究表明垂直于流动方向的电场作用力导致了流动阻力增加，从而增加了压降。Cotton等[15]研究了直流和交流电压作用下电场对环形管中HFC-134a流动沸腾压降的影响，表明电场易于增加气态与液态的动量相互作用，使得压力损失沿管道长度的方向，电压越强，压降越大。McGranaghan等[16]研究了60 Hz交流电压对细通道内HFE7000流动沸腾特性的影响，结果表明电压为0到3 kV时，压降随电压的增加而缓慢增加，电压为4～7 kV时，压降随电压增大而迅速增加，电压为7～10 kV时，压降随电压的增加而减小。由于线状与针状电极电场在不同电压范围下对流动沸腾传热强化的规律不同，且目前电场对流动沸腾压降影响的研究较为匮乏，本文以R141b为试验工质，在分别布置针状和线状2种电极的矩形微细通道内进行流动沸腾试验，探究针状与线状2种电极电场对微细通道内R141b流动沸腾压降特性的影响。

1 材料与方法

1.1 试验系统

试验系统如图1a所示，该试验系统主要由注液与循环流动模块、工质环境控制模块（加热、冷却、压力控制）、试验段模块、试验数据采集模块4部分组成，其中试验段模块如图1b所示，试验段模块由固定盖板、聚四氟乙烯盖板、玻璃板、铝制微细通道板、加热板及铝制基座组成。试验开始前进行氮气充气检漏，而后抽真空进行试验，铝制基座侧面布有测温孔和测压孔。本试验平台的铝制矩形微通道试验段由微细电火花切割而成，规格为240 mm×42 mm×19 mm，由6条矩形微通道组成，其中单个微细通道宽ch为2 mm，微细通道高ch为2 mm，2个相邻微细通道之间的壁厚w为5 mm，整块微细通道板长为240 mm，板宽为42 mm，板高为19 mm，铝制微细通道板示意图如图1c所示。

图1 试验平台

在微通道流动沸腾传热过程中，通道换热表面特性参数对微细通道的流动沸腾阻力特性有一定影响，课题组前期对此进行了大量研究，如文献[17]研究了微细通道表面润湿性对通道内流动沸腾摩擦阻力的影响，研究发现超疏水表面微细通道内的总压降、两相摩擦压降最大，普通光滑表面次之，超亲水表面微细通道内的总压降、两相摩擦压降最小，这是因为润湿性差的表面产生的汽泡受到的表面张力较大，汽泡脱离直径大、脱离频率低，更容易发展成受限汽泡，挤压壁面，导致压降增大。文献[18]研究了微细通道表面粗糙度对通道内流动沸腾摩擦阻力的影响，结果表明壁面粗糙度大的微细通道内两相摩擦压降大，原因是粗糙度大的表面会生成更多的汽泡使通道内流动扰动更加剧烈及通道内形成的液膜表面不平滑，使得压降增大。此外，通道表面特性参数对传热特性也有一定影响，如通道表面粗糙度越大，传热强化效果越好[19]，超疏水表面ONB的过热度最低、超亲水次之、普通光滑表面最高[20]等。目前对于将电场作用于微细通道换热系统的研究较少，电场影响微细通道内压降的研究也十分匮乏，影响机理也不完全清楚，电场影响压降的机理需要进一步的研究。本文在微细通道上分别布置针状、线状2种电极，对电极施加直流电压产生电场，研究针状、线状电极产生的电场对微细通道内流动沸腾压降的影响，并采用COMSOL软件对电场分布进行模拟，从而对电场影响压降的机理进行解释，为设计新型微细通道换热器提供思路。

1.2 电极布置方案

为了探究电场作用下微细通道内R141b的流动沸腾压降特性，本文采用2种不同的电极安装方式，在微细通道的局部相关位置设置针状及线状电极，通过施加不同强度的直流电源，在微细通道内形成不同强度的非均匀电场。同一通道相邻针状电极之间距离为25 mm，针状电极是由不锈钢硬丝加工而成，并将其竖直向下固定于聚四氯乙烯盖板上的电极孔内，针状电极长度为21 mm，电极顶端位于单个矩形微细通道几何中心线上方，针状电极放电顶端距离通道底面1.7 mm，单个矩形通道上针状电极的布置方式如图2a所示，单个矩形通道上方布有7根针状电极，共42根。线状电极由不锈钢软丝加工而成，线状电极的一端从聚四氯乙烯盖板上第一排某个电极孔穿入，另一端由盖板上第六排对应的电极孔穿出，两端拉紧固定使线状电极紧贴在通道上方的盖板内侧，并与矩形通道中心线相平行，单个线状电极长167 mm，线状电极中心轴线距离相应矩形通道底面1.85 mm，每个矩形通道上方布有一根线状电极，共6根。单个线状电极的布置方式如图2b所示。为了尽量减少电极本身对微细通道内所产生的汽泡运动的阻力影响，2种电极的直径均为0.3 mm，且聚四氟乙烯盖板上每个电极孔做了密封处理。

由于R141b制冷剂在不同压力下饱和温度不同，本文系统设计压力为140 kPa，工质对应饱和温度为41.6 ℃，为使微细通道内饱和沸腾段以泡状流为主，选择工质入口温度为32.5 ℃，在质量流率范围为277.35～531.75 kg/(m2·s)，热流密度范围为7.50～21.49 kW/m2，电压范围为0～850 V的试验工况下进行2种电极布置的电场作用下微细通道内R141b流动沸腾试验。

注：L1为沿制冷剂流动方向第一个电极孔距微细通道入口的距离，m；L2为单个微细通道上相邻电极孔之间的距离，m。

2 试验数据处理

2.1 微细通道内流动沸腾压降数学模型

试验工质R141b自下而上流经竖直放置的试验段，其入口温度32.5 ℃比相应压力下的R141b饱和温度低，随着微细通道下加热板的不断加热，试验工质R141b从一开始的单相液态逐渐变为两相气液混合态，因此存在单相段长度sp和两相段长度tp，其计算为[21]

式中sp为微细通道内单相段长度，m；tp为微细通道内两相段长度，m；p,l为工质液相定压比热容，kJ/(kg·K)；为试验工质的质量流量，kg/s；sat为制冷剂基于局部压力的饱和温度，K；ave为微细通道试验段的平均有效热流密度，kW/m2；ch为微细通道数量，ch=6；ch为单个矩形微细通道的宽度，m；w为单个矩形通道侧壁宽度的2倍，m；in为试验工质的入口温度，K；为试验段总长，m；e,n为传递给试验段第对测温点的局部有效热流密度，kW/m2，1≤≤4；为铝的导热系数；up,n为第对测温点的上测点温度，K；d,n为第对测温点的下测点温度，K；为上、下测量壁面温度点之间的距离，m。

试验段R141b总压降Δtot的值为工质入口压力in减去出口压力out，由单相摩擦压降Δsp,f、单相重力压降Δsp,g、两相加速度压降Δtp,a、两相摩擦压降Δtp,f、两相重力压降Δtp,g、进口突缩压降Δc以及出口突扩压降Δe组成，其中单相摩擦压降与单相重力压降组成单相流压降Δsp，两相摩擦压降压降、两相加速度压降与两相重力压降组成两相流压降Δtp，对于两相压降的计算，采用分相模型，总压降为4个部分，分别为摩擦压降Δf、重力压降Δg、加速度压降Δa和进出口压降Δce：

重力压降由式（5）计算[22]

式中L、G分别为液相、气相密度，kg/m3；e,out为微细通道出口热力平衡干度[23]。

加速度压降由式（6）计算[22]

通道进出口压降为

通道总摩擦压降为

由于在微细通道流动沸腾试验中的两相摩擦压降是非常重要的参数，而改变热流密度的大小、电场强度等条件会对微细通道内两相段的长度造成影响，因此本文以单位两相流摩擦压降大小Δftp作为重要参考。

其中

式中h为水力直径，m；sp为单相段摩擦系数；为宽高比，chch；LO为全液相雷诺数。

2.2 误差分析

试验采用精度为0.5%的LWGY涡轮流量计，精度为0.5%的HC3160-HVG4压力传感器，精度为0.2%的MIK-ST500温度变送器和WRNK-191热电偶，根据误差传递原理[25]，本文主要物理量误差如表1所示。

表1 主要物理量的误差

3 结果与分析

3.1 电场作用下的压降

3.1.1 无电极作用下的压降

图3为在试验运行压力为140 kPa、工质入口温度为32.5 ℃、质量流率为277.35～531.75 kg/(m2·s)、热流密度为7.50～21.49 kW/m2微细通道内R141b的流动沸腾压降。由图3a可知，在质量流率为277.35～531.75 kg/(m2·s)的范围内，随着质量流率的增加，重力压降、加速度压降、进出口压降的变化趋势皆不明显，而摩擦压降随微细通道内R141b质量流率的增加而增加；由图3b可知，在热流密度为7.50～21.49 kW/m2的范围内，随着热流密度增加，摩擦压降、加速度压降、进出口压降呈现增大的趋势，重力压降随热量流密度的增大而减小。由图3可知，摩擦压降在总压降中占据了最大的比例，为66.2%～73.3%，重力压降、加速度压降、进出口压降在总压降中占比分别为14.9%～28.8%、3.6%～9.4%、1.5%～2.8%。

注：图3a中，P=140 Pa，Tin=32.5 ℃，qave=12.46 kW·m-2；图3b中，P=140 Pa，Tin=32.5 ℃，G=277.35 kg·m-2·s-1；ΔPf、ΔPg、ΔPa、ΔPce分别为摩擦压降、重力压降、加速度压降、进出口压降，Pa。

3.1.2 针状和线状电极作用下的压降

图4a和图4b为针状和线状电极在不同电压作用下微细通道内R141b流动沸腾各压降组成的变化情况。

Note: P=140Pa, Tin=32.5℃.

由图4a和图4b可知，在同一热流密度和质量流率下，针状电极电场作用下的摩擦压降会随着电压增加而增加，重力压降、加速度压降与进出口压降无明显变化规律；线状电极电场作用下，摩擦压降在0～250 V电压范围内无明显变化，在250～850 V电压范围内随电压的增加而增加，重力压降、加速度压降与进出口压降在0～850 V电压范围内无明显变化规律。针状电极电场作用下，摩擦压降在总压降中占比67.8%～76.3%；线状电极作用下，总压降中摩擦压降的占比为69.8%～78.6%，与无电极作用下摩擦压降在总压降中占比例对比可知，2种电极作用下的摩擦压降在总压降中所占比例均有增大，其原因是由于电场对汽泡产生作用，使汽泡运动变得更加复杂，汽泡运动对流体产生扰动，引起摩擦压降增大。

3.2 电场对单位长度两相摩擦压降的影响

图5为2种电极作用下的微细通道内R141b流动沸腾单位长度两相摩擦压降随电压变化曲线图。由图5可知，在本文试验的热流密度范围内，2种电极作用下的单位长度两相摩擦压降随着热流密度的增加而增大，这是因为高热流密度下，微细通道内沸腾更加剧烈，流体的质量含气率与汽泡

速度加快[26]，使得气液、液体与微细通道壁面摩擦加强；在同一热流密度下，单位长度两相摩擦压降随电压的增大而增大。针状电极作用下，施加0、250、400、550、700和850 V电压下的微细通道内随热流密度变化的R141b流动沸腾平均单位长度两相摩擦压降比无电极作用下的平均单位长度两相摩擦压降分别提升0.7%、6.6%、9.2%、11.8%、13.5%、15.4%；线状电极作用下施加0、250、400、550、700和850 V电压下的平均单位长度两相摩擦压降比无电极作用下分别提升1.3%、2.0%、10.7%、13.6%、16.8%、18.7%。由此可知，在运行压力140 kPa、工质入口温度为32.5 ℃、质量流率为277.35 kg/(m2·s)工况下，2种电极本身对压降的影响非常微弱。由于在微细通道内施加电压后，电场对汽泡的作用力使汽泡朝着低强度场强的方向运动，此方向指向加热壁面[27]，汽泡在电场力的作用下会被沿轴向拉长[28]，在受到指向加热壁面力的作用下，微细通道内拉长的汽泡会挤压微细通道壁面的液膜，随着电压的升高，电场作用下汽泡受到的指向加热壁面的力增大，汽泡与液膜之间、液体与微细通道壁面之间的摩擦也随之增大，使得本文试验条件下，单位长度两相摩擦压降随电压的增大而增大。

Note: P=140 Pa, Tin=32.5 ℃, G=277.35 kg·m-2·s-1.

3.3 电极对压降的影响

图6为2种电极作用下各压降随电压的变化图。由图6可知，在运行压力140 kPa、工质入口温度为32.5 ℃、质量流率为277.35 kg/(m2·s)时，高热流密度时2种电极作用下的总压降、两相流压降、单位长度两相摩擦压降均明显高于低热流密度下的对应值。同一热流密度下施加电压超过250 V后，2种电极作用下各压降随电压的增大而增大，线状电极作用下的微细通道内R141b流动沸腾总压降、两相流压降、单位长度两相摩擦压降皆比对应针状电极下的高，2种电极的总压降、两相流压降、单位长度两相摩擦压降相对于0 V时皆有明显的提升。

Note: P=140 Pa, Tin=32.5 ℃, G=277.35 kg·m-2·s-1.

电压250 V时，线状电极电场作用下微细通道内R141b流动沸腾总压降、两相流压降、单位长度两相摩擦压降与线状电极施加0 V电压时相比并无明显的变化，平均变化绝对值分别为0.8%、0.5%、0.2%，而针状电极作用下的总压降、两相流压降、单位长度两相摩擦压降与针状电极施加0 V电压时的相比都有明显的增加，甚至超过对应线状电极作用下的压降数值。

为探究在本文试验工况下，针状与线状电极作用下微细通道内流动沸腾压降随电压不同变化的原因，使用COMSOL软件对微细通道内电场的分布进行模拟。由于针状电极的电场强度在其3 mm距离处开始出现严重的衰减，所以选取6 mm试验段可以有效地表征通道内电场分布，因此本文计算区域大小设为2 mm×2 mm×6 mm。边界条件设置如下：通道底面和2个侧壁面为固体域，并设置为零电势边界；电极表面为高电势边界。本文高电势的值设置为各试验电压，电极材料为304不锈钢，有限元计算区域介质为R141b，与文献[29]的研究类似，本文将R141b的介电常数假设为恒定值，且将任意时刻的电场状态看作为稳态，根据高斯定理提出的静电场的方程[30]为

式中为介电常数；为电势，V；e为电荷密度，C/m3。

图8为针状电极与线状电极布置方式下6 mm长微细通道内电场分布模拟结果。图8中的模拟结果中电场强度大小的数量级与文献[11]中一致，可本文的模拟结果具有一定的准确性。由图8可知，与线状电极相比，同一电压下针状电极在微细通道内形成的电场强度更大，微细通道内电场强度随着与电极之间距离的增大而减小，施加电压越大，微细通道内同一位置电场强度也随之越大，250 V时，针状电极电场强度最大值为2.75×106V/m，线状电极电场强度最大值为1.96×106V/m。250 V之前，针状电极电场作用下微细通道内R141b流动沸腾总压降、两相流压降、单位长度两相摩擦压降随电压增加而增加，线状电极电场作用下的压降无明显变化，电压超过250 V后，线状电极作用下微细通道内R141b流动沸腾总压降、两相流压降、单位长度两相摩擦压降才开始增加，所以250 V时线状电极作用下对汽泡产生的指向加热壁面电场力的大小为临界值[11]，电压低于250 V时，线状电极电场对汽泡的作用力没有对汽泡产生有效的作用，只有当电场对汽泡产生指向加热壁面的作用力w超过临界值，汽泡受到电场的作用明显，汽泡对流体的扰动作用才会导致压降明显增大。250 V电压时，针状电极由于其结构的特殊性，针状电极形成的电场强度高，针状电极作用下的电场对汽泡产生指向加热壁面的作用力w超过临界值，电场对汽泡产生显著的影响，这使得250 V电压下针状电极对压降的影响比线状电极显著。

注：A点为线状电极边缘处任意一点，B点为针状电极边缘处任意一点。

图8 针状与线状电极布置方式下微细通道内电场分布模拟结果

图9为2种电极作用下电场分布的有效范围对比，由图9可知，由于针状电极的设置较为分散，在微细通道内所形成的电场有效作用范围相对于线状电极更小，当电压达到400 V后，线状电极作用下的超过临界值，2种电场作用下的电场力超过临界值后，线状电极电场作用范围大，微细通道中有更多的汽泡受到电场力的作用，针状电极电场作用范围小，受到电场作用的汽泡少，电场是通过对汽泡的作用影响压降，所以线状电极的总压降、两相流压降、单位长度两相摩擦压降均大于相应条件下的针状电极。

图9 400 V时2种电极电场分布的有效范围

4 结 论

本文设置2种电极布置方式（针状和线状），以R141b为试验工质，在2种电场作用下的微细通道内进行流动沸腾传热试验，通过研究电场作用下的微细通道内R141b的流动沸腾压降特性，研究结果可为施加电场提高微细通道换热器的性能，降低电场对微细通道流动沸腾压降的影响，为实现微细通道高效节能提供新思路，并得出以下主要结论：

1）在本文试验工况下，针状电极电场作用下摩擦压降会随着电压增加而增加；线状电极电场作用下的摩擦压降在0～250 V电压范围内无明显变化，在250～850 V电压范围内随电压增加而增加。针状电极电场作用下摩擦压降在总压降中占比约为67.8%～76.3%，线状电极电场作用下的占比约为69.8%～78.6%，2种电极电场作用下的摩擦压降在总压降中的占比相对于无电极作用下均有升高。

2）本文试验工况下，电场作用下的单位长度两相摩擦压降随热流密度的增大而增大，同一热流密度下，单位长度两相摩擦压降随电压的增大而增大。在质量流率为277.35 kg/(m2·s)、热流密度为21.49 kW/m2、电压为850 V条件下，针状电极作用下的平均单位长度两相摩擦压降相对于无电极作用下提升15.4%，线状电极电场作用下提升18.7%。

3）本文试验工况下，施加电压为0～250 V时，线性电极电场对汽泡产生指向加热壁面的作用力小于临界值，线状电极电场对压降无明显作用，针状电极由于其结构的特殊性而产生尖端效应，在低电压下会产生高电场力，使针状电极电场在0～250 V电压内对压降的影响效果大于线状电极电场；电压达到400 V后，线状电极电场对汽泡产生指向加热壁面的作用力大于临界值，且线状电极电场的有效作用范围比针状电极电场更大，线状电极电场对压降的作用更为显著。

[1]罗小平，王文，张超勇，等. 换热器铝基微细通道微纳结构表面制备及其传热特性[J]. 农业工程学报，2018，34(20)：216－224.

Luo Xiaoping, Wang Wen, Zhang Chaoyong, et al. Surface preparation and heat transfer characteristics of aluminum-based microchannel micro-nanostructures in heat exchangers[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(20): 216－224. (in Chinese with English abstract)

[2]刘晓英，焦学磊，徐志刚，等. 可调LED光源系统设计及其对菠菜生长的影响[J]. 农业工程学报，2012， 28(1)：208－212.

Liu Xiaoying, Jiao Xuelei, Xu Zhigang, et al. Design on LED flexible light system and its effect on growth of spinach[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(1): 208－212. (in Chinese with English abstract)

[3]闫素英，李洪阳，史志国，等. 太阳能电池冷却用微通道散热器内纳米流体换热特性[J]. 农业工程学报，2016，32(13)：212－217.

Yan Suying, Li Hongyang, Shi Zhiguo, et al. Heat transfer characteristics of nanofluid in microchannel applied on solar cell cooling[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(13): 212－217. (in Chinese with English abstract)

[4]李伟钊，盛伟，张振涛，等. 热管联合多级串联热泵玉米干燥系统性能试验[J]. 农业工程学报，2018，34(4)：278－284.

Li Weizhao, Sheng Wei, Zhang Zhentao, et al. Experiment on performance of corn drying system with combination of heat pipe and multi-stage series heat pump equipment[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(4): 278－284. (in Chinese with English abstract)

[5]Ogata J, Yabe A . Basic study on the enhancement of nucleate boiling heat transfer by applying electric fields[J]. International Journal of Heat and Mass Transfer, 1993, 36(3): 775－782.

[6]Liu X, Hu C, Li H, et al. Effects of an inhomogenous electric field on an evaporating thin film in a microchannel[J]. International Journal of Thermophysics. 2018, 39(3): 43.

[7]Nangle Smith S, Cotton J S. EHD based load controllers for R134a convective boiling heat exchangers[J]. Applied Energy, 2014, 134: 125－132.

[8]Grassi W, Testi D. Electrohydrodynamic convective heat transfer in a square duct[J]. Annals of the New York Academy of Sciences. 2009, 1161(1): 452－462.

[9]Cotton J S, Chang J S, Shoukri M, et al. Electrohydrodynamically enhanced flow boiling in an eccentric horizontal cylindrical channel[J]. College of Electrical Engineering and Computer, 2008, 15(3): 220－223.

[10]黄烜，李瑞阳，郁鸿凌，等. 电场分布对R123沸腾换热的影响[J]. 化工学报，2007，58(8)：1926－1930.

Huang Xuan, Li Ruiyang, Yu Hongling, et al. Effects of electric field distribution on R123 boiling heat transfer enhancement[J]. Journal of Chemical Industry and Engineering (China), 2007, 58(8): 1926－1930. (in Chinese with English abstract)

[11]Zhang J, Luo X, Feng Z, et al. Effects of pin and wire electrodes on flow boiling heat transfer enhancement in a vertical minichannel heat sink[J]. International Journal of Heat and Mass Transfer. 2019, 136: 740－754.

[12]Seyed Yagoobi J, Bryan J E. Enhancement of heat transfer and mass transport in single-phase and two-phase flows with electrohydrodynamics[J]. Advances in Heat Transfer, 1999: 33, 95－186.

[13]Costa Patry E, Olivier J, Michel B, et al. Two-phase flow of refrigerants in 85m-wide multi-microchannels: Part II-Heat transfer with 35 local heaters[J]. International Journal of Heat and Fluid Flow. 2011, 32(2): 464－476.

[14]Bryan J E, Seyed Yagoobi J. Electrohydrodynamically enhanced convective boiling of alternative refrigerants: fundamental understanding and applicability[J]. HVAC&R Research. 2002, 8(4): 337－355.

[15]Cotton J, Robinson A J, Shoukri M, et al. A two-phase flow pattern map for annular channels under a DC applied voltage and the application to electrohydrodynamic convective boiling analysis[J]. International Journal of Heat and Mass Transfer. 2005, 48(25): 5563－5579.

[16]Mcgranaghan G, Robinson A J. EHD augmented convective boiling: Flow regimes and enhanced heat transfer[J]. Heat Transfer Engineering. 2014, 35(5): 517－527.

[17]罗小平，廖政标，周建阳，等. 极端润湿性微细通道内R141b的流动沸腾压降特性[J]. 华南理工大学学报：自然科学版，2018，46(2)：109－117.

Luo Xiaoping, Liao Zengbiao, Zhou Jianyang, et al. Influence of two-phase flow physical properties on emulsification using novel parallelized multi-scale channels[J]. Journal of South China University of Technology: Natural Science Edition, 2018, 46(2): 109－117. (in Chinese with English abstract)

[18]刘波，罗小平，王二利. 粗糙度对微通道内两相流摩擦压降的影响[J]. 中南大学学报：自然科学版，2015(11)：4334－4340.

Liu Bo, Luo Xiaoping, Wang Erli, et al. Effect of surface roughness on two-phase fractional pressure drop through microchannels[J]. Journal of Central South University: Science and Technology, 2015(11): 4334－4340. (in Chinese with English abstract)

[19]周建阳，罗小平，冯振飞，等. DMLS微换热器粗糙度对Al2O3/R141b流动沸腾传热影响[J]. 农业机械学报，2016，47(11)：421－428.

Zhou Jianyang, Luo Xiaoping, Feng Zhenfei, et al. Effect of roughness on flow boiling heat transfer of Al2O3/R141b in DMLS micro heat exchanger[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016,47(11): 421－428. (in Chinese with English abstract)

[20]罗小平，王文，廖政标，等. 基于不同润湿性微细通道过冷沸腾起始点（ONB）的实验研究[J]. 化工进展，2018，37(3)：884－892.

Luo Xiaoping, Wang Wen, Liao Zhengbiao, et al. Experimental study on onset of nucleate boiling (ONB) in different wettability micro-channels[J]. Chemical Industry and Engineering Progress, 2018, 37(3): 884－892. (in Chinese with English abstract)

[21]Lee H, Park I, Mudawar I, et al. Micro-channel evaporator for space applications-1. Experimental pressure drop and heat transfer results for different orientations in earth gravity[J]. International Journal of Heat and Mass Transfer. 2014, 77(4): 1213－1230.

[22]Kim S M, Mudawar I. Consolidated method to predicting pressure drop and heat transfer coefficient for both subcooled and saturated flow boiling in micro-channel heat sinks[J]. International Journal of Heat and Mass Transfer, 2012, 55(13/14): 3720－3731.

[23]阎昌琪. 气液两相流[M]. 哈尔滨：哈尔滨工程大学出版社，2007.

[24]Geiger G E. Sudden Contraction Losses in Single and Two-phase flow[D]. Pittsburgh: University of Pittsburgh, 1964.

[25]Mofat R J. Describing the uncertainties in experimental results[J]. Experimental Thermal & Fluid Science, 1988, 1(1): 3－17.

[26]Choi K, Pamitran A S, Oh J, et al. Pressure drop and heat transfer during two-phase flow vaporization of propane in horizontal smooth minichannels[J]. International Journal of Refrigeration, 2009, 32(5): 837－845.

[27]Cheung K, Ohadi M M, Dessiatoun S, et al. EHD-enhanced boiling coefficients and visualization of R-134a over enhanced tubes[J]. Journal of Heat Transfer, 1997, 119(2): 332－338.

[28]Zu Y Q, Yan Y Y. A numerical investigation of electrohydrodynamic (EHD) effects on bubble deformation under pseudo-nucleate boiling conditions[J]. International Journal of Heat and Fluid Flow, 2009, 30(4): 761－767.

[29]Diao Y H, Liu Y, Zhang J, et al. Effect of electric field on the enhanced heat transfer characteristic of an evaporator with multilayered sintered copper mesh[J]. Journal of Electrostatics, 2015, 73: 26－32.

[30]Kasayapanand Nat. Numerical modeling of the effect of number of electrodes on natural convection in an EHD fluid[J]. Journal of Electrostatics, 2007, 65(7): 465－474.

Effect of electric field on flow boiling pressure drop characteristics of R141b in microchannel

Luo Xiaoping, Peng Zizhe, Liu Qian, Guo Feng, Zhang Jinxin

(,,510640,)

Microchannel heat exchanger is highly efficient in heat exchange due to its microscopic structure and large heat exchange surface area ratio. Compared with conventional heat exchanger, it not only improves heat exchange efficiency but also saves space and reduces the weight of heat dissipation equipment. In agriculture, microchannel heat dissipation technology has been widely used in cultivation, machinery, and product drying. Microchannel increases resistance to fluid flow and could thus result in an increase in energy consumption. It could also give rise to a pressure-drop oscillation affecting operation of the system. Previous work showed that electric field has a beneficial impact on heat transfer efficiency in the microchannel heat exchanger. The purpose of this paper is to investigate how different electrodes (needle type and linear types) affect pressure drop of the boiling R141b coolant when it flows in the microchannel. The pressure in the system was kept at 140 kPa, the working inlet temperature was 32.5 ℃, the mass flow rate ranged from 277.35 to 531.75 kg/(m2·s), the heat flux density ranged from 7.50 to 21.49 kW/m2, the voltage ranged from 0 to 850 V, and the microchannel was rectangular with a cross section of 2 mm×2 mm. The results showed that the electric field under both electrodes increased the frictional pressure drop in the microchannel, compared with that without electric field. The frictional pressure drop of the two-phase along per unit length increased with both voltage and heat flux. The average two-phase frictional pressure drop along per unit length under the needle and the linear electrode was increased 0.7% to 15.4% and 1.3% to 18.7%, respectively, that without electric field. When the voltage changed from 0 to 250 V, the effect of the needle electrode on the pressure drop was greater than that of the linear electrode, while when the voltage was higher than 400 V, the effect of the linear electrode on the pressure drop was more significant. We simulated the electric field distribution in a microchannel 6 mm long using COMSOL for both electrodes. The results showed that under the same voltage, the maximum of the electric field intensity in the needle electrode was higher than that in the linear electrode, but the effective range of the electric field in the linear electrode was higher than that in the needle electrode. The results presented in this paper provide an alternative to improve performance of microchannel heat exchanger and reduce its energy consumption.

refrigerant; electric field; electrode; microchannel; flowing boiling; pressure drop

罗小平，彭子哲，刘倩，郭峰，章金鑫. 电场对微细通道内R141b制冷剂流动沸腾压降的影响[J]. 农业工程学报，2020，36(1)：257－265.doi：10.11975/j.issn.1002-6819.2020.01.031 http://www.tcsae.org

Luo Xiaoping, Peng Zizhe, Liu Qian, Guo Feng, Zhang Jinxin. Effect of electric field on flow boiling pressure drop characteristics of R141b in microchannel[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(1): 257－265. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.01.031 http://www.tcsae.org

2019-07-27

2019-10-29

国家自然科学基金资助项目（21776096）

罗小平，博士，教授，博士生导师，主要研究方向为微通道换热器相变传热、分子动力学模拟等研究。Email：mmxpluo@scut.edu.cn

10.11975/j.issn.1002-6819.2020.01.031

TK124

A

1002-6819(2020)-01-0257-09