Experimental investigation on anti-icing mechanism and characteristics of superhydrophobic electrothermal coupled surface
-
摘要: 作为一种新型防冰技术,超疏水电热复合表面防冰具有良好的防冰效果和较低的能量消耗。基于超疏水表面水滴撞击及润湿特性,依据结冰表面热平衡理论,发展了超疏水电热复合表面防冰热流密度预测模型。在结冰风洞中开展了圆柱模型超疏水电热复合表面防冰实验研究,结果表明,防冰热流密度理论计算值与实验值之间的差别小于6%,验证了该预测模型的准确性。实验结果与能耗分析表明:与传统电加热方法相比,超疏水电热复合表面防冰能够有效降低防冰能耗;在风速10 m/s、液态水含量1 g/m3、水滴平均体积直径65 μm、温度−15 ℃条件下,超疏水表面能够有效防止回流冰形成;对于干、湿表面防冰,超疏水电热复合表面防冰比传统电加热方法能够分别降低约43%和33%的防冰能耗。Abstract: As a novel anti-icing technology, superhydrophobic electrothermal coupled surface anti-icing possesses an excel-lent anti-icing efficiency with low energy consumption. Based on the water droplet impact behaviors and the wetting characteristics of the superhydrophobic surface, a prediction model of the heat flow density of superhydrophobic electrothermal coupled surface anti-icing is developed according to the thermal balance theory of the icing surface. The experimental analysis of the superhydrophobic electrothermal coupled surface anti-icing is carried out in a low-speed icing wind tunnel. The results show that the difference between the theoretical anti-icing heat flux and the experimental results is less than 6%, which verifies the prediction model. The analysis of the experimental results and energy consumption shows that the superhydrophobic electrothermal coupled surface anti-icing effectively reduces the energy consumption compared with the electrothermal method. With the freestream velocity of 10 m/s, liquid water content of 1 g/m3, mean volume diameter of 65 μm, and temperature of −15 ℃, the superhydrophobic coating can effectively prevent the formation of backwater due to its wetting property. For dry and wet surface anti-icing, the superhydrophobic electrothermal coupled surface anti-icing method reduces the energy consumption by about 43% and 33% respectively compared with the electrothermal method.
-
Keywords:
- electrothermal /
- superhydrophobic /
- ice wind tunnel experiments /
- anti-icing /
- low power
-
0 引言
流体壁面剪应力的精确测量是流体力学领域的难题[1]。随着MEMS技术的发展,硅基[2-3]和柔性聚合物基[4]的壁面剪应力微传感器研究发展受到重视,而传感器标定是其测量的基础和前提。根据不同应用需求初步形成了不同类型的动、静态壁面剪应力标定方法和装置[5-9]。本文介绍了平板流动、管道流动和旋转流动等静态标定方法,以及平面波振动、阶梯圆盘旋转运动等动态标定方法,介绍了工作原理和典型结构装置,可为壁面剪应力传感器标定方法的选定提供参考。
1 壁面剪应力静态标定方法
壁面剪应力传感器的静态标定需要施加稳定已知量值的剪应力,读取对应输出。采用流体力学经典模型,通过测量其他物理量从而换算得到对应壁面剪应力值是常用方法。由于之前长期缺乏有效的流体壁面剪应力传感器,相关标定研究不系统,尚未形成统一的标准[10-11]。
目前,常见标定方法主要有平板模型标定法、管道流动标定法[12]和旋转运动标定法[13]。
1.1 平板模型标定法
早在20世纪上叶就有学者和机构对平板模型边界层进行了大量理论和实验研究,并得到平板表面层流和湍流边界层壁面剪应力的理论计算模型。
在空气中,均匀来流绕过一个厚度很薄的平板时,如果平板迎角为0,可以应用普朗特边界层微分方程解决粘性流动问题,平板边界层如图 1所示[14]。
理论上,由于平板绕流压强顺流不变,即
,不会发生边界层分离,也没有压差阻力,只有作用于壁面的摩擦阻力。故在这种情况下壁面切应力可由下式计算:(1) 式中:μ表示流体的动力粘度,ν表示流体的运动粘度,U∞表示平板前端均匀来流流速,x表示距平板前缘的距离。
随着来流速度增大,雷诺数也将增大。当平板表面雷诺数大于临界雷诺数,绕流场可能由层流状态转变为湍流[15],如图 2所示。
通过动量积分方程可以得到湍流边界层中壁面剪应力的计算公式:
(2) 式中:xt为湍流边界层前沿到测量点的距离,ρ为空气密度,δ2为动量厚度,可通过下式计算得到:
(3) 通过确定测量点到转捩位置的距离xt,就可以得到来流速度U∞与剪应力τ的函数关系。
由于平板表面发生转捩的位置不易测得,常常采用强制转捩的方法获得全湍流状态的平板边界层。Tani对转捩带粗糙体高度进行了研究,并得出粗糙体的极限高度计算公式[16]:
(4) 式中:k′ crit为粗糙体的极限高度。
平板模型标定法能产生数百帕量级的壁面剪应力,标定精度受风洞 (或水道) 来流速度测量误差、测试点距前缘距离测量误差及环境温度测量误差等因素的影响。然而,平板边界层的流场状态不易控制,可结合实验方法 (如升华法、油流法等) 初步判断平板表面流动转捩与分离情况。
1.2 管道流动标定法
管道流动因其运动受到管壁约束,管道内的流动参数 (如流量、温度和壁面静压等) 较易测量,也是常见的剪应力标定方法。
1.2.1 管道流动动量损失标定法
当方形管道的宽高比较大时[17],入口处来流稳定的通道内流体运动可近似为二维流动。管道内流体的动量平衡与壁面剪应力作用相关[18]。通过测量管壁沿程压力损失,可以得到壁面剪应力的大小,标定装置如图 3所示。
定义pinlet/poutlet为管道内入口处压力和出口处压力的比值。对于压力比pinlet/poutlet≤1.1,可忽略管道内流体的压缩性,完全发展的不可压缩层流管道内的压力分布是线性的,剪应力表达式可以简单表示为:
(5) 式中:h为管道高度,
表示管道内沿流场方向的压力梯度。通过测量Δx距离沿程的压力变化Δp,就可以得到剪应力的值。马赫数大于0.3的流动问题要考虑流体压缩性。对于压力比pinlet/poutlet≥1.1,管道内压力分布不再是线性变化[19]。通过非线性压力Π=p(x)/poutlet与非线性距离X之间的关系,代入公式 (5) 可以得到大压力比下的管道剪应力计算公式:
(6) 式中:Lc为管道的长度,poutlet为管道出口处的压力,X为沿管道方向x位置与管道长度Lc的比值。常数a和b能够从Π2=aX+b的曲线计算得到[20]。
1997年,美国麻省理工学院研发了一种基于动量损失原理的空气介质剪应力标定装置[21],量程0~1Pa,如图 5所示。在微型气槽内壁面开凿压力检测孔,由于流体在槽腔运动产生摩擦阻力引起流体动量的损失,导致沿气流运动方向有着不同的壁面压力分布,换算可得压差中点剪应力值。此外,还考虑了空气压缩性,给出了不同压差范围的剪应力计算方法并用此装置对自主研发的浮动元件式剪应力传感器进行了标定。
2002年,美国凯斯西储大学Patel等在研究壁面剪应力传感器时采用了此原理搭建起一套小型的扁平管道标定装置,并对浮动式剪应力传感器进行了标定,标定最大剪应力达到80Pa[18]。
2015年,中国空气动力研究与发展中心李建强等研究建立了管道流动动量损失原理的空气介质剪应力标定装置,增设了精密温控系统,可实现从环境温度到55℃的温度控制[22]。
2014年,中国船舶科学研究中心建立了水介质的动量损失管流标定装置,并利用LDV、PIV等对MEMS热膜剪应力传感器进行了测试标定[7]。
1.2.2 微型层流管道标定法
流体在通过极低雷诺数的微型扁薄矩形水道 (具有较大宽高比w/h),可获得可控的槽内纯剪切流[23],如图 6所示。
当管道内的流场形态为均匀层流时,其速度剖面曲线为二次曲线,管道体积流量Q与流动中心处最大流速相关,根据二维定常流动流体壁面剪应力的定义式[24],可以导出管道内壁面剪应力与体积流量之间的关系:
(7) 因矩形通道横截面宽度w无法实现远大于截面高度h,在剪应力计算时引入矩形通道的形状修正因子
,则修正后剪应力计算公式为[23]:(8) 式中:Dh为水力学半径,
;Φ(n) 为矩形通道的修正因子,该修正因子的引入可降低扁薄矩形水道侧壁对中心流速的影响效应[17],
(9) 2005年,美国加州理工学院的Xu Y.等人根据该流动原理搭建起适用于MEMS剪应力传感器标定的微型管道测试平台,管道高度仅0.3mm,并实现对传感器0~5Pa范围内的剪应力标定[23]。
2009年,中国西北工业大学依此原理建立了微型扁薄标定水槽。2014年,李雁冰等根据该装置原理,对恒流驱动的热膜剪应力传感器进行了水下标定实验[25],并给出标定模型。王雷涛、严宇超等增加了溢流系统和温度控制系统,并对系统误差进行了分析。
1.3 旋转运动标定法
平板模型和管道流动标定装置均需要有气源或者水压动力源作为介质流动的动力,尤其搭建专门风洞较为困难,因此有研究者采用圆盘或者圆筒的旋转运动带动周边流场流动从而形成一定的壁面剪切应力。
2004年,美国国家航空航天局 (NASA) Zuckerwar等采用旋转圆轮建立了标定装置[26],如图 8所示。1个电机速度可控的圆形旋转轮和1个与其同轴心的固定圆弧槽,传感器平齐安装在槽上且靠近轮的边缘。当旋转圆轮与传感器之间的间隙h极小时,圆轮与圆弧槽间的流体形成简单库艾特流动,产生的剪应力为:
(10) 式中:r为轮子半径,ω为电机带动的旋转轮的转速,通过调节转速的大小可以获得需要的剪应力。
Zuckerwar通过对圆轮的精确安装,控制圆轮和同轴圆弧槽之间的缝隙最小达50μm,实现最大发生剪应力达到234.49Pa。
1992年,美国麻省理工学院的Shajii J[27]等采用锥板式[28-30]的剪应力标定装置对浮动式剪应力传感器标定进行过相关研究,如图 9(a)所示。2006年,韩国浦项科技大学机械工程系的Kim[31]等采用2个平行放置的圆盘建立了壁面剪应力标定装置,如图 9(b)所示,并对热膜剪应力传感器实现0~0.8Pa范围内的标定。另外,美国加州理工学院的Brown G. L.[32]等采用圆筒式粘度计[33-34]的构造,如图 9(c)所示,搭建了剪应力标定装置。2012年,日本名古屋大学的Terashima O.[35]等对圆筒剪应力发生装置进行了改进,将外圆筒改为平板,如图 9(d)所示,更便于传感器的安装。
旋转运动剪应力标定装置的搭建较为简单,但要保证产生剪应力的精度,对旋转装置偏心度要求较高,且若要产生较大范围的剪应力,圆盘或圆筒的直径和旋转速度要求很大。通常,此类方法仅用于0~2Pa范围内的剪应力标定。此外,该标定方法流场状态复杂,主要是利用其平均效应。
2 动态剪应力标定方法
动态剪应力标定需要加载一个周期性变化的剪应力。通常动态标定是为了研究传感器在中频或者高频激励下的响应情况,因此需要动态标定装置能产生较高频率的激励信号。
2.1 平波管标定方法
平波管实现动态测量原理为用扬声器产生声波,声波驱动流体产生动态剪应力信号,利用麦克风测量产生信号的频率实现对剪应力传感器的动态测量。管道尾部加入消声装置以消除回声对测量的影响。
2000年,美国佛罗里达大学的Chandrasekaran等建立了以平面声波驱动的动态剪应力标定装置,如图 10所示,该装置可产生最大0.1Pa,频率范围100Hz到20kHz的动态剪应力信号,如图 11所示。利用该装置对带有隔热空腔的硅基热敏传感器进行了动态测试[9]。
2001年,美国麻省理工学院的Sheplak等采用相同原理的平波管动态标定装置,并对浮动式剪应力传感器动态性能进行了标定测试[11]。
平面声波将引起振荡的压力梯度,通过测量压力波动,可以换算得到壁面动态变化的剪应力[11]:
(11) 式中:p′为压力波动的幅值。
2.2 圆盘动态标定法
2006年,韩国浦项科技大学机械工程系的Kim等在研究静态圆盘旋转标定法的同时对转动圆盘进行了改造,采用阶梯式转动圆盘,得到最高激励频率400Hz动态变化的剪应力,其结构如图 12所示。阶梯厚度为δh,2个圆盘的间隙的改变为
。在每一个环形库艾特流周期内,圆盘上的m个阶梯提供了m个剪切脉冲。当上部圆盘的转速为n r/min,旋转产生的阶梯脉冲频率为mn/60Hz[31]。3 结论
介绍了国内外流体壁面剪应力标定方法及其原理和装置。通过对比分析,可以得到以下几点启示:
(1) 采用平板模型标定方法可以依托现成风洞或者水道,可以获得较大发生范围的剪应力 (从0帕到上百帕),实验系统可以适配多种类型的剪应力传感器;
(2) 微型管道流槽标定是一种紧凑便捷的专用标定装置,截面尺寸及雷诺数很小,水槽高度经常在毫米级以下 (百微米量级)。发生的剪应力范围通常在10Pa以内;
(3) 动量损失剪应力标定装置也属于专用定制,使用比较方便,剪应力发生范围极大 (可达上百帕);
(4) 旋转运动标定方法的装置较易搭建,对机械装置及装配等精度要求较高,发生的剪应力范围较小,其流场状态复杂,主要取用平均效应;
(5) 动态标定方法中研究及使用较多的是平波管法。
综上,可根据各类流体壁面剪应力传感器的工作要求、尺寸、安装条件 (贴附面为平面或曲面、安装高度等) 和应用环境的不同,选用不同的标定方法或装置。
目前,不同标定方法流场状态及品质影响,误差和不确定度及精度分析等研究较少,尚需关注。
-
表 1 结冰风洞基本性能参数
Table 1 Basic performance parameters of icing wind tunnel
参数 取值范围 实验段尺寸 700 mm(长)×304 mm(宽)×500 mm(高) 风速 5~18 m/s 温度 −25~30 ℃ 液态水含量 0.3~1.0 g/m3 水滴平均体积直径 10~100 μm 表 2 实验条件
Table 2 Experimental condition
风速/(m·s−1) 温度/℃ 液态水含量
/(g·m−3)水滴平均体积
直径/μm10 −15,25 1 65 表 3 不同防冰方法及干/湿表面条件下所需热流密度
Table 3 Heat flux required under different anti-icing methods and dry/wet surface conditions
防冰方法及干/湿表面条件 热流密度/(W·m−2) 电热−干表面防冰 1923 电热−湿表面防冰 824 超疏水电热−干表面防冰 1099 超疏水电热−湿表面防冰 549 超疏水电热−干表面防冰(理论计算值) 1166.79 -
[1] 卜雪琴,彭珑,林贵平,等. 热气防冰系统内表面弦向传热性能衰减规律[J]. 北京航空航天大学学报,2016,42(1):72-78. DOI: 10.13700/j.bh.1001-5965.2015.0065 BU X Q,PENG L,LIN G P,et al. Chordwise attenuation of heat transfer performance on inner surface of hot-air anti-icing system[J]. Journal of Beijing University of Aeronautics and Astronautics,2016,42(1):72-78. doi: 10.13700/j.bh.1001-5965.2015.0065
[2] 郁嘉,赵柏阳,卜雪琴,等. 某型飞机发动机短舱热气防冰系统性能数值模拟[J]. 空气动力学学报,2016,34(3):302-307. DOI: 10.7638/kqdlxxb-2015.0212 YU J,ZHAO B Y,BU X Q,et al. Numerical simulation of the performance of an engine nacelle hot-air anti-icing system[J]. Acta Aerodynamica Sinica,2016,34(3):302-307. doi: 10.7638/kqdlxxb-2015.0212
[3] 卜雪琴,林贵平,郁嘉. 机翼电加热防冰表面内外传热的耦合计算[J]. 航空动力学报,2010,25(7):1491-1496. DOI: 10.13224/j.cnki.jasp.2010.07.030 BU X Q,LIN G P,YU J. Coupled heat transfer calculation on an airfoil electrothermal anti-icing surface[J]. Journal of Aerospace Power,2010,25(7):1491-1496. doi: 10.13224/j.cnki.jasp.2010.07.030
[4] 胡林权. 民用飞机机翼电加热防/除冰应用现状及技术难点[J]. 航空科学技术,2016,27(7):8-11. HU L Q. Application status and technical difficulties for civil aircraft wing electrothermal anti-/ de-icing[J]. Aeronautical Science & Technology,2016,27(7):8-11.
[5] 马莉娅,熊联友,刘立强,等. 用于碳纤维复合材料的电热除冰技术实验研究[J]. 航空学报,2012,33(1):54-61. MA L Y,XIONG L Y,LIU L Q,et al. Experimental study on electro-thermal deicing technique for carbon fiber composite[J]. Acta Aeronautica et Astronautica Sinica,2012,33(1):54-61.
[6] JIANG H,WANG H T,LIU G,et al. Light-weight, flexible, low-voltage electro-thermal film using graphite nanoplatelets for wearable/smart electronics and deicing devices[J]. Journal of Alloys and Compounds,2017,699:1049-1056. doi: 10.1016/j.jallcom.2016.12.435
[7] POURBAGIAN M,HABASHI W G. Surrogate-based optimization of electrothermal wing anti-icing systems[J]. Journal of Aircraft,2013,50(5):1555-1563. doi: 10.2514/1.C032072
[8] MOHSENI M,AMIRFAZLI A. A novel electro-thermal anti-icing system for fiber-reinforced polymer composite airfoils[J]. Cold Regions Science and Technology,2013,87:47-58. doi: 10.1016/j.coldregions.2012.12.003
[9] CORNELL J S,PILLARD D A,HERNANDEZ M T. Comparative measures of the toxicity of component chemicals in aircraft deicing fluid[J]. Environmental Toxicology and Chemistry,2000,19(6):1465-1472. doi: 10.1002/etc.5620190601
[10] LOUCHEZ P, BERNARDIN S, LAFORTE J L. Physical properties of aircraft de-icing and anti-icing fluids[C]//Proc of the 36th AIAA Aerospace Sciences Meeting and Exhibit. 1998: 575. doi: 10.2514/6.1998-575
[11] YAKHYA S, MORENCY F. Numerical study of momentum and heat transfer in propylene glycol jets used for aircraft ground deicing[C]//Proc of the 23rd AIAA Computational Fluid Dynamics Conference. 2017: 4504. doi: 10.2514/6.2017-4504
[12] 李斌. 飞机除冰/防冰液及除冰技术[J]. 清洗世界,2012,28(1):26-31. DOI: 10.3969/j.issn.1671-8909.2012.01.006 LI B. Brief survey of deicing/anti-icing fluid and techniques for aircraft[J]. Cleaning World,2012,28(1):26-31. doi: 10.3969/j.issn.1671-8909.2012.01.006
[13] 李清英,白天,朱春玲. 飞机机械除冰系统的研究综述[J]. 飞机设计,2015,35(4):73-77. DOI: 10.19555/j.cnki.1673-4599.2015.04.015 LI Q Y,BAI T,ZHU C L. Research summary of mechnical de-icing systems of aircrafts[J]. Aircraft Design,2015,35(4):73-77. doi: 10.19555/j.cnki.1673-4599.2015.04.015
[14] DRURY M,SZEFI J T,PALACIOS J. Full-scale testing of a centrifugally powered pneumatic de-icing system for helicopter rotor blades[J]. Journal of Aircraft,2015,54:220-228. doi: 10.2514/1.C033965
[15] SOMMERWERK H, HORST P, BANSMER S. Studies on electro impulse de-icing of a leading edge structure in an icing wind tunnel[C]//Proc of the 8th AIAA Atmospheric and Space Environments Conference. 2016: 3441. doi: 10.2514/6.2016-3441
[16] YEONG Y H, LOTH E, SOKHEY J, et al. Ice adhesion strength on hydrophobic and superhydrophobic coatings[C]//Proc of the 6th AIAA Atmospheric and Space Environments Conference. 2014: 2063. doi: 10.2514/6.2014-2063
[17] CAO L L,JONES A K,SIKKA V K,et al. Anti-icing superhydrophobic coatings[J]. Langmuir,2009,25(21):12444-12448. doi: 10.1021/la902882b
[18] KULINICH S A,FARHADI S,NOSE K,et al. superhydrophobic surface: Are they really ice-repellent?[J]. Langmuir,2011,27(1):25-29. doi: 10.1021/la104277q
[19] FARHADI S,FARZANEH M,KULINICH S A. Anti-icing performance of superhydrophobic surfaces[J]. Applied Surface Science,2011,257(14):6264-6269. doi: 10.1016/j.apsusc.2011.02.057
[20] DE PAUW D,DOLATABADI A. Effect of superhydrophobic coating on the anti-icing and deicing of an airfoil[J]. Journal of Aircraft,2016,54(2):490-499. doi: 10.2514/1.C033828
[21] 彭兰清,卫建勋,陈诺,等. 基于超疏水表层的石墨烯电热除冰实验研究[J]. 科学技术与工程,2021,21(15):6513-6518. DOI: 10.3969/j.issn.1671-1815.2021.15.057 PENG L Q,WEI J X,CHEN N,et al. Experimental study on graphene electrothermal deicing based on superhydrophobic surface[J]. Science Technology and Engineering,2021,21(15):6513-6518. doi: 10.3969/j.issn.1671-1815.2021.15.057
[22] KIMURA S, YAMAGISHI Y, SAKABE A, et al. A new surface coating for prevention of icing on airfoils[R]. SAE Tech-nical Paper, 2007. doi. org/10.4271/2007-01-3315
[23] ZHAO Z H,CHEN H W,LIU X L,et al. Novel sandwich structural electric heating coating for anti-icing/de-icing on complex surfaces[J]. Surface and Coatings Technology,2020,404:126489. doi: 10.1016/j.surfcoat.2020.126489
[24] 朱宝. 低能耗超疏水电热蒙皮设计及防冰性能研究[D]. 西安: 西北工业大学, 2018. DOI: 10.27406/d.cnki.gxbgu.2018.000407 ZHU B. Low power superhydrophobic electrothermal skin and its anti-icing performance[D]. Xi'an: Northwestern Polytechnical University, 2018.
[25] KOLBAKIR C, HU H Y, LIU Y, et al. A hybrid anti-/ de-icing strategy by combining NS-DBD plasma actuator and superhydrophobic coating for aircraft icing mitigation[C]//Proc of the AIAA Scitech 2019 Forum. 2019: 0050. doi: 10.2514/6.2019-0050
[26] ANTONINI C,INNOCENTI M,HORN T,et al. Understanding the effect of superhydrophobic coatings on energy reduction in anti-icing systems[J]. Cold Regions Science and Technology,2011,67(1-2):58-67. doi: 10.1016/j.coldregions.2011.02.006
[27] HU H Y, AL-MASRI F, HU H. An experimental study on ice accretion and anti-/ de-icing of a pitot tube[C]//Proc of the AIAA SCITECH 2022 Forum. 2022: 1913. doi: 10.2514/6.2022-1913
[28] XUE S J,LIU Y H,WANG Y,et al. Variation in anti-icing power of superhydrophobic electrothermal film under different temperatures and wind speeds[J]. International Journal of Aerospace Engineering,2022,2022:3465428. doi: 10.1155/2022/3465428
[29] 朱春玲, 朱程香. 飞机结冰及其防护[M]. 北京: 科学出版社, 2016. ZHU C L, ZHU C X. Aircraft icing and its protection[M]. Beijing: Science Press, 2016.
-
期刊类型引用(5)
1. 陈纪仲,田海平,丁俊飞,马国祯,陈怡纯. 壁湍流摩阻时空信号无接触测量实验研究. 实验流体力学. 2024(04): 66-74 . 本站查看
2. 高南,刘玄鹤. 实用化壁面切应力测量技术的综述与展望. 空气动力学学报. 2023(03): 1-24 . 百度学术
3. 刘丹. 免标定壁面剪应力测量技术研究. 科技创新与应用. 2023(34): 188-192 . 百度学术
4. 陆连山,李栋. 利用双参数非线性拟合求解湍流边界层壁面摩擦速度. 西北工业大学学报. 2020(02): 253-260 . 百度学术
5. Xuanhe Liu,Zhuoyue Li,Nan Gao. An improved wall shear stress measurement technique using sandwiched hot-film sensors. Theoretical & Applied Mechanics Letters. 2018(02): 137-141 . 必应学术
其他类型引用(4)