-
工业过程中经常遇到需要快速冷却高温物体或者在有限空间内持续冷却运动高温板带的场合,如高负荷微电子元件在狭小空间内需快速排出电流产生的焦耳热以维持正常工作温度[1];带材连续生产线上需要高且均匀的热流密度以对运动带材实施气体淬火,保证马氏体组织转变均匀[2-3];航空发动机涡轮叶片内部通道通过实施强制冷却以使其在高温环境下正常运转[4-6]。从传热的角度讲,上述冷却过程都需要比较大的换热系数,气体冲击射流冷却可以适应这种需求。冲击冷却通常是气体从圆孔或狭缝型喷嘴喷射到固体表面,由于其流动行程短,边界层薄,因而换热系数比常规的管内流动换热系数高出几倍直至一个数量级。对于给定的换热系数,采用冲击换热原理设计的装置比常规装置在尺寸上通常小1~2个数量级[7-9]。
在圆孔气体射流换热方面,前人多是采用实验方法研究冲击换热系数,获得传热准数方程。Gardon等[10-11]揭示了冲击射流局部换热系数非单调性变化是由于湍流度变化的缘故,并以喷嘴间距Xn为定性尺寸定义努谢尔数(Nu),引入到达速度定义雷诺数(Re),获得的准数方程在目前的工业设计中被广泛采用。Martin[12]以Krotzsch[13]的传质实验数据为基础,通过传热传质类比,总结出适用性更广的单股圆孔射流和圆孔阵列射流准数方程。Martin公式以喷嘴直径D为定性尺寸,综合考虑了各种几何因素(如喷射高度、开孔率)、喷气速度与介质种类等的影响,并讨论了喷射系统的最佳结构参数问题。Kuntikana[14]和Viskanta[15]通过实验分析了气体冲击射流和火焰冲击射流的传热过程,研究结果均针对特殊的几何结构。Zuckerman[16]列举了多位研究者的单孔和多排圆孔射流的传热准数方程,Meola[17]探讨了对多排圆孔射流建立新的通用传热关联式的方式。
通常冲击射流的传热准数影响因素可归因于流体特性和结构参数,流体特性主要指雷诺数和普朗特数(Pr),结构参数通常指量纲为一喷射高度(H/D)和量纲为一喷嘴间距(Xn/D)或者开孔率(Af)。但是对比已有文献研究可知,这两个量纲为一参数并不能涵盖阵列圆孔射流的所有结构参数,如采用平板开孔、短圆管和长圆管的结果就有显著差异[18],此外射流废气的排出方式,如强制排气、排气不畅、排气良好对换热系数也会影响显著[19]。因此Martin的回归公式与Glaser[20]的实验值存在明显偏差,当Re>20 000时,Glaser的努谢尔数高出回归值约20%~30%,当Re<20 000时,Ott[21]的实验结果也高出Martin公式计算值的25%左右。Meola等[22]也曾指出Gardon的实验值与Martin公式的计算值差异较大,Gardon的实验值较大。可见,目前还难以获得通用的圆孔阵列射流换热准数方程。
近年来,随着电子工业和制造业等的发展,气体冲击射流冷却需要更高的换热系数和更好的均匀性,如带钢热处理领域生产高强度钢板所需的平均换热系数可达1 000 W/(m2·K)。这意味着一方面需要进一步提高雷诺数,另一方面需要进一步优化冲击射流阵列结构参数,例如采用密集型的气体射流装置。王磊等[23]研究了涡轮叶片用密集冲击射流冷却的换热情况,考察了小喷嘴间距条件下的换热性能。Meola[22]测量了方形、圆孔等喷嘴的密集冲击射流换热系数,结果表明换热系数比Martin公式计算值高出20%左右。
本文针对刀片式密集圆孔气体冲击换热进行实验研究,以得到平均对流换热系数,用于指导冷却装置的开发、设计与运行。
-
图1所示为刀片式圆孔射流装置,喷孔布置在狭长型的刀片状肋管上,类似于狭缝喷嘴,肋管之间形成规则的狭长形逸气通道,排气空间较大;在肋管头部布置间距为Yn、直径为D的圆孔,相邻肋管的间距为Xn。这种结构使得刀片式圆孔气体射流装置兼具狭缝喷嘴和圆孔喷嘴的优点,一方面肋管的集气功能可以大幅降低供风和集气环节的阻力损失,在风机功率一定的情况下可以提高气体流量;另一方面良好的排气空间,使得喷孔可以密集布置,提高了开孔率。
刀片式圆孔气体射流换热均匀,其测量原理如图2所示。射流速度为ug,温度为Tg的气体从冲击高度为H的刀片式多排圆孔喷嘴喷射出来,垂直向下喷射到温度为Ts的高温薄铜板上进行强化换热。由于铜板厚度极薄,毕渥数(Bi)不超过0.01,忽略其厚度方向的温度差别。高温薄铜板背部黏贴有电压为U、电流为I的薄膜电阻加热器供热。射流气体与薄铜板换热后,废气流向相邻肋管间的逸气空间,最终从两侧排出。当薄铜板表面温度稳定在Ts后,薄铜板与气体的换热满足如式(1)所示的热平衡式:
式中:qe为薄膜电阻加热产生的热流密度,W/m2;A为冲击板面积,m2;hav为平均对流换热系数,W/(m2·K);qloss为包含薄铜板对环境的辐射散热和加热器背面导热散热引起的热流密度损失,W/m2。由此平均对流换热系数为
定义射流雷诺数(ReD)和平均努谢尔数(NuDav)分别为
式中:kg和νg分别为气体的导热系数和运动黏性系数,单位分别为W/(m·K)和m2/s;定性温度采用气体喷口处的温度值Tg。改变射流速度和喷射高度,测量各工况下的平均对流换热系数,可以获得NuDav与ReD的关系。
-
图3所示为圆孔气体射流实验原理,主要分为空气管路系统、刀片式喷嘴喷射系统、冲击板与热流密度测试系统、位置调节系统和信号采集系统。
室温空气经过滤器和消音器后由压头为12 000 Pa,风量为3 600 m3/h的离心式风机鼓入直径400 mm的风管,经主管道上的调节阀和孔板流量计后进入风箱。在风箱前,管道逐渐由圆管缓慢过渡到风箱上矩形接口,以保证气体流量在喷嘴组上均匀分配。
刀片式圆孔喷嘴的具体尺寸如图4所示,结构参数包括喷孔直径D、肋管深度E、喷嘴纵向间距Xn、横向间距Yn(见图5)
、喷射高度H和喷嘴个数n。开孔率(Af)表示喷嘴出口总面积与换热面总面积的比值,图5中分别示出了狭缝喷嘴(5(a))、等间距圆孔喷嘴(5(b))和刀片式圆孔喷嘴(5(c))的开孔率定义。刀片式圆孔喷嘴的具体几何尺寸如表1所示,其开孔率高达0.036 6,是常规管式风箱开孔率的3倍,接近狭缝喷嘴的开孔率。当 D=8 mm时,如果带钢与喷嘴间距过小,例如 H/D <5,则带材振动容易损坏喷嘴,将不具备实际用途;如果带钢与喷嘴间距过大(H/D >15),则会造成风机的动力浪费。由于刀片式圆孔射流的应用场景为带材连续退火,考虑到带材顺利运行和喷嘴效率,选取6<H/D<13,具体为6.25、7.50和12.5。 D/mm Xn/D Yn/D H/D E/mm n Af 8 6.25 3.125 6.25, 7.50, 12.5 150 182 0.036 6 表 1 刀片式圆孔喷嘴的几何尺寸表
Table 1. Geometric data for blade-type circular nozzle
冲击板的作用是模拟气体射流冲击后形成的流场,其材质为不锈钢,几何尺寸为600 mm×600 mm×8 mm。在冲击板的中央镶嵌着边长250 mm的方形薄铜板,其作用是模拟高温带材表面与气体射流的换热,并测定平均热流密度。高温薄铜板与冲击板之间填充厚约10 mm的绝热材料,以减小薄铜板边缘与冲击板之间的接触导热,安装时使薄铜板上表面与不锈钢冲击板面保持在同一平面上。
总管冷却气体的流量采用孔板流量计测量,误差约为3%。铜板温度采用T型热电偶测量,测温范围0~350 ℃,最大温度误差为0.5 ℃;气体温度采用水银温度计测量,测温范围0~100 ℃,最大温度误差为0.2 ℃。电压表为1.0级,测量范围0~200 V,最大误差为2 V。根据误差分析[24],可得雷诺数的相对误差为3%,平均努谢尔数或平均对流换热系数的相对误差为
误差的主要来源为薄膜电阻热流密度的测量误差。
为了确定薄铜板的热流密度损失(qloss),在正式实验前,将加热器通电,使薄铜板升温,但不启动气体射流风机。经过一段时间后,电加热器的供热量与热流计表面的辐射散热,以及未计入的导热散热量达到平衡,并使热流计表面维持温度不变。根据薄铜板的热平衡可得
${q_{\rm{e}}} = {q_{\rm{loss}}}{\rm{ = }}{h_{\rm{loss}}}({T_{\rm{s}}} - {T_{\rm{g}}})$ ,其中${h_{\rm{loss}}}$ 为热损失换热系数。通过数次实验可以确定${h_{\rm{loss}}} \approx 9.2\;{\rm{W/}}({{\rm{m}}^2}\cdot {\rm{K}})$ ,由于${h_{\rm{loss}}}$ 远小于${h_{{\rm{av}}}}$ ,实验过程中热流计表面温度控制在80~110 ℃,可认为${h_{\rm{loss}}}$ 不随热流计表面温度变化而改变。 -
采用的刀片式圆孔喷嘴组横向与纵向间距较小(Xn/D=6.25, Yn/D=3.125),局部对流换热系数可视为均匀分布,采用薄铜板热流计表面的平均温度来计算平均努谢尔数NuDav。但实际上,换热系数在空间上仍然存在一定的不均匀性,如果采用表面温度的最高值和最低值估算努谢尔数的最大值NuDav_MAX和最小值NuDav_MIN,它们与平均努谢尔数NuDav的偏差可以表征换热系数在空间上的不均匀程度。图6示出了H/D =7.50时,以NuDav_MAX−NuDav为正偏差、NuDav−NuDav_MIN为负偏差的NuDav与ReD的对应关系。由图6可以看出,当雷诺数较高时,NuDav的偏差较大,但是最大偏差仍小于5%,平均偏差为2%。由此可见,可以采用热流计的平均温度估算NuDav。
图7示出了不同H/D条件下ReD对NuDav的影响,主要研究了H/D分别为6.25,7.50和12.5时的3种情况。由图7可以看出,NuDav随雷诺数的增加而增大,即换热系数随着喷气速度的增大而增加。表明提高射流速度会使气体射流与冲击面的动量交换加剧,从而增强换热效果。由图7还可以看出NuDav对ReD的依赖程度。3种H/D条件下,在双对数坐标下,实验点均落在3条近乎平行的线附近,表明3种H/D条件下准数方程中ReD的幂次接近。当H/D<7.50时,喷射高度的影响并不明显,当H/D>7.50时,增加喷射高度会显著降低换热系数。
图8示出了
$N{u_{D{\rm{av}}}}{\rm{/}}Re_{D}^{m}$ (其中m为指数)与量纲为一喷射高度$H{\rm{/}}D$ 的关系,按照$(N{u_{D{\rm{av}}}}{\rm{/}}Re_{D}^{m}) \propto {(H/D)^p}$ (其中p为指数)的形式进行回归,可得到$p = - 0.606$ ,说明NuDav随$H{\rm{/}}D$ 的−0.606次幂递减。 -
在6 000≤Rex≤600 000,8≤Xn/D≤32条件下,Gardon获得的准数方程如下[11]
式中,ua为到达冲击面的速度,m/s,考虑了喷射高度的影响。Martin总结的准数方程为[12]
公式的适用范围为2 000≤ReD≤100 000,2≤H/D≤12,0.004≤Af≤0.04。
图9比较了H/D =7.5时函数的测量值与Gardon公式的计算值,可以看出两者吻合很好。
图 9 H/D=7.5测量值与Gardon公式比较
Figure 9. Comparision of measured data with values reduced by Gardon’s equation when H/D =7.5
图10(a)比较了H/D =7.5时函数的测量值与Martin公式的计算值,可以看出测量值明显高于Martin公式的计算值,前者大约是后者的1.4倍左右。这与Glaser[20]和Ott[21]实验结果与Martin公式的比较结果是一致的。
图 10 H/D =7.5时实验值与Martin公式(a)和修正Martin公式(b)结果比较
Figure 10. Comparision of measured data with values reduced by Martin’s equation (a) and corrected Martin’s equation (b) when H/D=7.5
图11示出了H/D=7.5时,函数
$K\left( {{A_{\rm{f}}},H/D} \right) $ 、$G\left( {{A_{\rm{f}}},H/D} \right) $ 及修正函数KA$\left( {{A_{\rm{f}}},H/D} \right) $ 分别与Af 的关系。由图11(a)可以看出,当H/D=7.5、Af>0.03时函数$K\left( {{A_{\rm{f}}},H/D} \right)$ 和$G\left( {{A_{\rm{f}}},H/D} \right)$ 几乎不变化;由图11(b)可以看出,当Af>0.02时,函数$K\left( {{A_{\rm{f}}},H/D} \right)$ ×$G\left( {{A_{\rm{f}}},H/D} \right)$ 几乎不变化。这主要是Martin公式中函数$K\left( {{A_{\rm{f}}},H/D} \right)$ 和$G\left( {{A_{\rm{f}}},H/D} \right)$ 过度压制了开孔率Af的作用。图 11 H/D =7.5时函数K(Af,H/D)、G(Af,H/D)和修正函数KA(Af,H/D)以及函数K(Af,H/D)× G(Af,H/D)与Af的关系
Figure 11. Relationship between function K (Af,H/D),G (Af,H/D) and modified KA (Af,H/D)、function K (Af,H/D) × G (Af,H/D) with Af when H/D=7.5
如果引入修正函数
来代替函数
$K\left( {{A_{\rm{f}}},H/D} \right)\times G\left( {{A_{\rm{f}}},H/D} \right)$ ,${K_{\rm{A}}}\left( {{A_{\rm{f}}},H/D} \right) $ 与Af的关系如图11(b)所示,可得到Martin公式的修正式如式(5)所示。图10(b)也对比了H/D =7.5时的实验数据与Martin修正公式结果,可以看出二者吻合得很好。 -
(1)以空气为模化介质,考察了刀片式多排圆孔冲击射流换热情况,结果表明平均努谢尔数随雷诺数的增加而增大,与量纲为一喷射高度
$H{\rm{/}}D$ 的0.606次幂成反比。(2)刀片式多排圆孔气体射流采用密集布置喷孔后,提高了开孔率,换热系数显著提高,约为Martin公式计算值的1.4倍。
(3)实验结果与Gardon公式非常吻合,与Martin公式存在显著差别,引入修正的KA(Af, H/D)函数后两者吻合较好。
刀片式多排密集圆孔气体冲击射流冷却的实验研究
Convective Cooling of Blade-Type Multiply Intensive Circular Gas Impinging Jets
-
摘要: 刀片式多排圆孔气体射流可用于受限表面的密集冷却。考察了喷射高度和雷诺数对密集圆孔气体射流冲击光滑平板换热效果的影响。应用热薄膜法测量壁面温度并利用等热流加热获得平均对流换热系数,在射流雷诺数(ReD)的取值范围为5 000≤ReD≤48 500、量纲为一喷射高度(H/D)分别为6.25、7.50、12.5时,平均努谢尔数随ReD增加而
递增,并与H/D的0.606次幂成反比。建立了平均努谢尔数关于ReD、普朗特数(Pr)和H/D的准则关系式,以及密集气体冲击射流冷却的换热准数方程。实验结果与Martin公式的模拟结果比较表明:刀片式圆孔喷嘴密集布置喷孔后,开孔率提高,换热系数大幅提高,约为采用Martin公式得到的换热系数计算值的1.4倍。实验数据与Gardon公式的计算结果非常吻合,与Martin公式的计算结果差别显著,但引入修正函数后与后者能很好地拟合。 -
关键词:
- 刀片式 /
- 多排密集圆孔气体射流 /
- 冲击冷却
Abstract: Blade-type multiply circular gas jets can be utilized for intensive cooling within limited area.The effects of impingement height H and jet Reynolds numbers on the heat transfer when multiply intensive circular air jets normally impinge on a smooth and flat surface are investigated experimentally. The average convective heat transfer coefficients are measured with heated thin foil technique which involves the detection of surface temperature and apples constant heat flux heating. When the jet Reynolds numbers (ReD) vary from 5 000 to 48 500 and normalized impingement height (H/D) is 6.25. 7.50, 12.5, respectively, an increase in jet Reynolds numbers will produce higher average Nusselt numbers, and normalized impingement height with an exponent of −0.606 also have a strong impact on it. The experimental results are compared with Martin’s correlation and about 40% higher is observed. It is indicated that intensive circular nozzle array can significantly enhance heat transfer rate due to its high open porosity. The experimental data are quite well fitted by the corrected Gardon's equation. Conversely, the correlation proposed by Martin is not suitable to deduce present data and a modified geometrical function is developed. -
表 1 刀片式圆孔喷嘴的几何尺寸表
Table 1. Geometric data for blade-type circular nozzle
D/mm Xn/D Yn/D H/D E/mm n Af 8 6.25 3.125 6.25, 7.50, 12.5 150 182 0.036 6 -
[1] 桂进乐. 基于微小通道与冲击射流的电子器件冷却机理研究[D]. 武汉: 华中科技大学, 2015.
[2] 郑海燕. 冷轧连续退火炉冷却技术的发展和应用[J]. 轧钢, 2018, 35(1): 52-56.
[3] 刘志伟, 郭晓宏, 张瑞琦, 等. 超高强度冷轧耐候马氏体钢的组织与性能[J]. 上海金属, 2017, 39(5): 47-50. doi: 10.3969/j.issn.1001-7208.2017.05.011
[4] 柳焕楠. 涡轮叶片冲击冷却仿真与实验研究[D]. 哈尔滨: 哈尔滨工程大学, 2017.
[5] XING Y G, ZHANS L, ZHAO Q, et al. Experimental and numerical investigations of overall cooling effectiveness on a vane end wall with jet impingement and film cooling[J]. Applied Thermal Engineering, 2019, 148(5): 1148-1163. [6] 刘亮亮, 竺晓程, 刘昊, 等. 纵向涡发生器对涡轮叶片前缘冲击冷却的影响[J]. 动力工程学报, 2018, 38(9): 719-724.
[7] 顾维藻, 神家锐, 马重芳, 等. 强化传热[M]. 北京: 科学出版社, 1990: 189-216.
[8] 杨世铭, 陶文铨. 传热学[M]. 第4版. 北京: 高等教育出版社, 2006: 277-280.
[9] INCROPERA F P, DEWITT D P, BERGMAN T L, et al. 葛新石, 叶宏译. 传热和传质基本原理[M]. 北京: 化学工业出版社, 2007: 279-283
[10] GARDON R, AKFIRAT J C. The role of turbulence in determining the heat transfer characteristics of impinging jets[J]. International Journal of Heat and Mass Transfer, 1965, 8(10): 1261-1272. doi: 10.1016/0017-9310(65)90054-2 [11] GARDON R, COBONPUE J. Heat transfer between a flat plate and jets of air impinging on it[C]//International Developments in Heat Transfer.[s.l.]: ASME, 1962: 454-460. [12] MARTIN H. Heat and mass transfer between impinging gas jet and solid surface[C]//Hartnett J P, Irvine J T F. Advances in Heat Transfer, New York: Academic Press Inc, 1977:1-60. [13] KROTZSCH P. Warme-und sroffubergang bei prallstromung aus dusenund blendenfeldern[J]. Chemie Ingenieur Technik, 1968, 40(7): 339-344. doi: 10.1002/cite.330400706 [14] KUNTIKANA P, PRABHU S V. Isothermal air jet and premixed flame jet impingement: Heat transfer characterisation and comparison[J]. International Journal of Thermal Sciences, 2016, 100(2): 401-415. [15] VISKANTA R. Heat transfer to impinging isothermal gas and flame jets[J]. Experiment Thermal and Fluid Science, 1993, 6(2): 111-134. doi: 10.1016/0894-1777(93)90022-B [16] ZUCKERMAN N, LIOR N. Jet impingement heat transfer: Physics, correlations, and numerical modeling[J]. Advances in Heat Transfer, 2006, 39: 565-631. [17] MEOLA C. A new correlation of Nusselt number for impinging plate[J]. Heat Transfer Engineering, 2009, 30(3): 221-228. doi: 10.1080/01457630802304311 [18] ETEMOGLU A B, KEMAL I . Investigation into the effect of nozzle shape on the nozzle discharge coefficient and heat and mass transfer characteristics of impinging air jets[J]. Heat and Mass Transfer, 2010, 46(11): 1395-1410. [19] DONG-HO R, PIL-HYUN Y, HYUNG H C. Local heat/mass transfer and flow characteristics of array impinging jets with effusion holes ejecting spend air[J]. International Journal of Heat and Mass Transfer, 2003, 46(6): 1049-1061. doi: 10.1016/S0017-9310(02)00363-0 [20] GLASER H. Untersuchungen an schlitz- und mehrdusenanord-nungen bei der trocknung feuchter oberflachen durch warm-luftstrahlen[J]. Chemie Ingenieur Technik, 1962, 34(7): 200-207. [21] OTT H. Warmeubergang an einerdutchLuftstrahlengekuhlten Platte[J]. Schweizer Bauzeitung, 1961, 79(46): 834-840. [22] MEOLA C, CARLOMAGNO G M. Intensive cooling of large surfaces with arrays of jets[C]// Quantitative InfraRed Thermography Conference, 2006, Padova, Italy:[s.n.]. http://dx.doi.org/10.21611/qirt.2006.064. [23] 王磊, 张靖周, 杨卫华. 密集型阵列冲击射流换热特性实验[J]. 航空动力学报, 2009, 24(6): 1264-1269.
[24] KLINE S J, MCCLINTOCK F A. Describing uncertainties in single sample experiments[J]. Mechanical Engineering, 1953, 75(1): 3-8. -