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, 2]。以斜激波与圆柱/圆球弓形激波之间的相互作用为典型模型,根据激波干扰位置和激波强度不同,Edney[3]总结出6种具有普适性的激波干扰类型(Ⅰ~Ⅵ类)。其中,第Ⅳ类激波干扰因其产生的超声速射流冲击壁面会带来严峻的热/力载荷而受到设计和研究人员的广泛关注[4-7]。在某些情况下,激波干扰产生的超声速射流高度弯曲,不会直接冲击壁面,这种激波干扰被称为第Ⅳa类[8]。
高超声速进气道流动存在多种形式的激波干扰现象,尤其以设计工况下前体激波入射进气道唇口[9]形成的激波干扰最为突出。针对前体激波与二元进气道唇口弓形激波的干扰问题,Wieting等[10]进行了系统性研究,获得了丰富的压力和热流数据。近年来,前缘高度后掠的内转式进气道[11-13]逐渐成为研究热点之一。然而,其普遍采用的V形溢流口(唇口)与二元进气道唇口具有明显的区别。最近,Xiao等[14]在不考虑前体激波的情况下,将内转式进气道溢流口简化为V形钝前缘,研究发现:在一定的几何参数下,V形钝前缘自身就会形成类似于Edney第Ⅳ类激波干扰的流动结构。显然,当前体激波入射V形钝前缘溢流口时,将产生更加复杂的波系结构,很可能加剧溢流口的气动热。
另一方面,根据来流条件的不同,进气道会出现亚额定、额定和超额定等3种工作状态[15]。针对二元进气道的研究表明[16-17]:在超额定状态下,前体激波入射在唇口内部,当前体激波与唇口激波发生Edney第Ⅱ类激波干扰时,形成的亚声速区可能导致进气道局部不起动。相对而言,内转式进气道唇口更加复杂,前体激波与其相互作用的研究尚未见系统报道,亟需针对性的机理揭示,以补充认识上的不足。
采用完整的内转式进气道对V形溢流口激波干扰特性进行研究,不仅会耗费大量的计算和实验资源,更会带来多种复杂因素的影响。抓住内转式进气道V形溢流口的主要流动特征,有侧重地进行简化,继而揭示其流动机理和关键影响因素,当属高效简洁的研究方式。因此,本文将内转式进气道V形溢流口简化为V形钝前缘平板。鉴于内转式进气道的前体激波为复杂的三维曲面激波,为便于准确描述和控制,本文暂不考虑入射激波的展向曲率,而是采用简化的斜激波,以突出V形钝前缘溢流口自身带来的主要特征。继而采用激波风洞实验和数值模拟相结合的方法,揭示入射激波与V形钝前缘溢流口相对位置变化引起的激波干扰形态演变,以期为内转式进气道溢流口设计提供有价值的参考。
1 实验和计算方法
1.1 实验设备和模型
V形钝前缘平板由较长的平直段和交叉位置倒圆部分组成,如图 1所示。交叉位置倒圆半径R=6.5mm,前缘钝化半径r=2mm,前缘扩张角β=48°。前缘平直段的长度L=31r,以使其产生的脱体激波充分发展。以V形前缘驻点位置为坐标原点,定义流向为x方向,平板展向为z方向,垂直于平板方向为y方向。
为了更加清楚地观察波系结构,对V形钝前缘平板实验模型进行了2.5倍放大。实验在中国科学技术大学的KDJB330反射型激波风洞[18-20]中进行。实验来流马赫数为6,静压为891Pa,静温为101K。纹影系统采用Z字形光路布置,如图 2所示。通过高速摄影记录流场纹影图像, 拍摄帧率为15000Hz,曝光时间为2μs。入射斜激波由气流偏转角为12°的斜劈产生,斜劈流向长度为200mm,展向宽度为260mm。为了获得丰富的流场信息,考虑到模型对光路的遮挡,分别从平行于平板(z方向,见图 2)和垂直于平板(y方向,见图 3)2个方向对流场进行拍摄。通过前后移动斜劈的方式,改变入射斜激波与V形钝前缘溢流口的相对位置。为了便于分析,斜激波入射位置采用不考虑激波干扰时无粘斜激波在平板上的入射点坐标。
1.2 计算方法
数值计算采用Fluent软件求解Navier-Stokes方程,并使用具有二阶精度的格式离散。对流通量采用AUSM格式,湍流模型采用k-ω SST模型,气体分子粘性采用Sutherland公式,气体定压比热使用CHEMKIN温度多项式进行计算。对于V形钝前缘平板流场的计算,由于来流和几何的对称性,为了节省计算资源,采用1/4计算域进行模拟。采用结构化六面体网格,计算域及局部网格如图 4和5所示,近壁面采用等比加密,壁面y+ < 1,总网格量约3000万。有斜激波入射时的计算域如图 6所示,此时六面体网格总数在5000万左右。为了与风洞实验进行对比,来流条件设置与实验一致。在数值模拟中监视各方程残差,监测出口位置的流量及质量平均马赫数、驻点位置的壁面压力和热流量。待残差收敛并且各监测参数稳定,认为流场收敛。
由于三维流动的复杂性,数值模拟获得的流场结构难以直接与具有光路积分效应的实验纹影进行对比。为此,本文按照纹影法的原理对数值模拟流场进行处理。如为了获得z方向穿过全流场的数值纹影,将流场沿z方向划分为k个等间距的区域,假设每个小区域内的流场参数沿z方向均匀分布,则经过流场的光线偏转角εx由式(1)[21]得到:
(1) 式中,n为折射率,ξ1和ξ2分别为光线进入流场和离开流场位置的z坐标。由此获得的光强变化可以与z方向穿过流场的实验纹影进行直观对比。
2 结果与分析
V形钝前缘自身的三维绕流会导致复杂的激波干扰,因此在研究斜激波入射V形钝前缘平板激波干扰前,有必要对V形钝前缘基本波系结构进行研究。
2.1 V形钝前缘平板基本流场
从y方向拍摄的流场纹影与流场x-z对称面数值纹影的对比如图 7(a)所示,可以看到数值模拟较好地捕捉到了基本的波系结构。图 7(b)为上半部分波系结构的示意图,其中DS代表V形钝前缘斜直段产生的脱体激波,BS代表V形前缘交叉位置形成的弓形激波,脱体激波DS入射弓形激波BS,相交于三波点T,在三波点位置产生透射激波TS和滑移线Σ。脱体激波DS波后的超声速气流在V形前缘交叉倒圆位置被压缩,压缩波汇聚形成弯曲激波CS。透射激波TS与弯曲激波CS发生异侧激波干扰,形成马赫杆MS。激波BS以及马赫杆MS波后气流均为亚声速,而透射激波TS波后气流为超声速,超声速气流在亚声速区域之间不断地压缩和膨胀,交替出现激波和膨胀波,形成超声速射流结构。超声速射流沿壁面运动,在驻点前发生碰撞,并卷起形成反向旋转的涡对结构。涡对结构与弓形激波BS相互作用,使得弓形激波波面向上游凸起,并出现Kink。将脱体激波DS视为入射激波,则V形钝前缘激波干扰结构与Edney第Ⅴ类激波干扰类似。流场中存在射流失稳和剪切层脱涡等流动过程,数值模拟虽然不能捕捉这些非定常过程,但能够刻画基本波系结构,流动的非定常性没有显著影响本文对波系干扰类型的研究。
为了探讨流场波系的三维结构,从z方向也进行了纹影拍摄。使用本文1.2节所描述的三维流场后处理方法得到的数值纹影与实验纹影图像进行了对比,如图 8(a)所示。由于流场的三维性,在z方向纹影图像上可以看到3道激波,从左到右依次为弓形激波BS上凸起的激波面(Bulge)、弓形激波BS和弯曲激波CS。图 8(b)进一步定量比较了数值模拟与实验获得的激波BS及其凸起的激波面的位置。结果表明:激波BS的位置符合较好,而凸起的激波面(Bulge)略有差别。这是由于激波BS相对稳定,而凸起的激波面受涡对结构的影响(见图 7),具有一定的非定常性。整体而言,数值模拟与实验纹影吻合较好,可用于本文后续的激波干扰研究。
图 9比较了相同钝化半径下二维圆柱产生的弓形激波和V形钝前缘平板弓形激波。其中,Δ1和Δ2分别代表V形钝前缘和圆柱弓形激波的脱体距离,V形钝前缘弓形激波由x-y对称面的数值模拟给出,圆柱产生的弓形激波形状由Billig[22]公式给出,而弓形激波面上声速点位置的激波角β*由以下公式[23]给出:
(2) 式中,Ma∞为来流马赫数, γ为气体比热比。从图 9可以看出,由于自身的激波干扰,V形钝前缘交叉位置弓形激波的脱体距离较大,并且驻点前存在大范围的亚声速区。进一步采用马赫数为5.95的等值面给出V形钝前缘平板的三维激波结构(图 10),可以看出:前缘脱体激波DS和弓形激波BS均为三维曲面激波,并且BS波面有较大范围的凸起。这表明V形钝前缘弓形激波的近正激波区域较大,当斜激波与弓形激波干扰时,出现第Ⅲ和第Ⅳ类干扰的范围更大[24]。
2.2 斜激波入射V形钝前缘溢流口激波干扰流场
当入射斜激波距离V形溢流口上游较远时,对应于亚额定工况下内转式进气道前体激波与唇口激波相交于唇口外的情况。图 11(a)为斜激波无粘入射平板位置坐标为x/r=-10.5时的流场纹影图片,图 11(b)为波系示意图。斜激波IS从平板下侧入射弓形激波BS超声速区域,发生同侧激波相交,出现Edney第Ⅵ类激波干扰。激波相互作用导致弓形激波的激波角减小,入射激波干扰后形成透射激波TS和剪切层Σ。
向下游移动斜劈,当斜激波无粘入射位置为x/r=-3时,斜激波入射在V形前缘驻点前的亚声速区,此时对应于内转式进气道额定工况下前体激波入射唇口的情况。仅从z方向拍摄纹影时,由于模型遮挡,难以获得干扰点附近激波干扰形态。因此,需要结合y方向的纹影来判断干扰类型。从图 12(a)流场z方向的纹影可以看出, 入射斜激波IS分为了两道斜激波,分别为IS和IS′。其中,IS′是入射斜激波IS与V形钝前缘的脱体激波DS干扰产生的透射激波。IS′入射V形前缘平板壁面形成反射激波RS,而入射斜激波IS与弓形激波BS的干扰点被模型遮挡。从图 12(b)流场y方向的纹影可以看到,在驻点前存在两道弓形激波BS和BS′。结合图 12(a)和(b)两个方向的纹影,推测激波干扰类型可能为Edney第Ⅳ类激波干扰,其流场结构如图 13所示。此时,由于斜激波IS入射,原有的弓形激波BS转变为通过透射激波TS连接的两段弓形激波BS和BS′,并伴随着超声速射流的产生,该流场结构很复杂,将在2.3节进一步分析。因此,在内转式进气道额定工况下,前体激波入射溢流口驻点前亚声速区域,可能带来复杂的激波干扰流场,对进气道流动带来负面影响。
进一步向下游移动斜劈,使得斜激波无粘入射位置为x/r=0.2,此时对应于内转式进气道超额定状态下前体激波入射溢流口内部的情况。流场纹影如图 14(a)所示,图 14(b)为其流场示意图。入射激波IS与弓形激波BS干扰形成的剪切层Σ并没有撞击壁面,而是形成一种有别于经典的Edney第Ⅲ类激波干扰的结构。其原因在于:V形钝前缘脱体激波DS与弓形激波BS干扰形成的脱体距离和亚声速区域较大,使得发生第Ⅲ类激波干扰的区域增大[23]。
当斜激波无粘入射位置为x/r=3时,斜激波与弓形激波BS的超声速区域干扰,流场纹影如图 15(a)所示,流场如图 15(b)所示。入射斜激波IS与弓形激波BS相交,形成长度较短的马赫杆MS、透射激波TS以及剪切层Σ1和Σ2。这种由Edney第Ⅱ类激波干扰形成的亚声速区,可能会带来进气道局部不起动[16-17],应避免出现。
2.3 斜激波/弓形激波第Ⅳ类激波干扰流场分析
当斜激波入射V形钝前缘时,三维波系干扰的情况如图 16所示(这里采用斜激波在上)。根据激波位置和特点,可能会出现3种干扰情况:第一种是入射斜激波IS与脱体激波DS相交,第二种是入射斜激波IS与弓形激波BS相交,第三种是入射斜激波与激波BS和DS的交线相交。鉴于三维激波干扰的复杂性,为了深入理解斜激波入射V形前缘时的激波干扰流场,本节以斜激波无粘入射壁面位置坐标为x/r=-5时作为典型工况,结合数值模拟进行分析。与入射位置为x/r=-3的工况(图 12)相比,此时激波干扰点位置更接近BS的正激波部分。
图 17(a)采用流向截面和展向对称面(z/r=0)展示了斜激波IS与V形钝前缘波系干扰三维结构的演化情况。其中,流向截面的马赫数云图仅显示了一半流场,而展向对称面通过马赫数等值线显示。在x/r=-5截面上(局部放大见图 17(b)),入射斜激波IS与脱体激波DS的较弱部分发生异侧相交,形成Edney第Ⅰ类干扰,产生的透射激波TS相对于激波IS向DS激波面偏折。在x/r=-2.5截面上(局部放大见图 17(c)),入射斜激波IS与脱体激波DS的较强部分发生异侧相交,形成Edney第Ⅱ类激波干扰,并产生马赫杆MS,此时的透射激波TS1入射钝前缘壁面。图 17沿流向发展的波系结构表明:由于入射斜激波IS与前缘脱体激波DS相交处DS激波强度的变化,激波干扰类型由Edney第Ⅰ类转变为第Ⅱ类。而在V形钝前缘溢流口处,如图 17(a)中的蓝色线框所示,入射激波IS与弓形激波BS发生干扰并产生超声速射流,这一区域的流场结构将在图 18中进行分析说明。
在斜激波入射条件下,多股气流共同汇集在V形钝前缘交叉位置,形成更加复杂的三维激波干扰流场。图 18给出了V形钝前缘交叉位置z/r=0和y/r=-0.4两个正交截面上的流场结构。从z/r=0截面可以看出,入射斜激波与弓形激波BS干扰后,形成通过透射激波TS4连接的两段弓形激波BS和BS′以及超声速射流。该超声速射流并未直接入射壁面,而是发生高度弯曲沿壁面方向溢流,此干扰类型是第Ⅳ a类干扰。而y/r=-0.4截面展示了斜激波IS同时与弓形激波BS和脱体激波DS发生多波系干扰的情况,由于IS与DS干扰已经形成了马赫杆MS(见图 17(c)),MS再与弓形激波BS干扰形成透射激波TS3以及超声速射流区域。可见,入射激波IS、弓形激波BS与脱体激波DS三者之间的干扰互相耦合。
将数值模拟流场沿y方向进行三维纹影处理,如图 19(a)所示,与图 19(b)给出的实验纹影图像符合较好。从三维纹影结果中可以看到,V形钝前缘溢流口弓形激波BS被分为两道激波BS和BS′。实验纹影图像与数值计算结果相互印证了斜激波与V形钝前缘弓形激波发生了第Ⅳa类激波干扰。
3 结论
通过简化构型对内转式进气道V形钝前缘溢流口激波干扰特性进行了研究,得到以下结论:
(1) V形钝前缘自身形成的激波干扰,导致其驻点前弓形激波的脱体距离较大,波后存在大范围的亚声速区,因此在估算内转式进气道溢流口激波脱体距离时,需慎用圆柱脱体激波距离计算公式。
(2) 斜激波入射在V形钝前缘驻点前弓形激波接近正激波的位置时,发生Edney第Ⅳa类激波干扰。此时,入射激波与弓形激波BS的干扰以及弓形激波BS与脱体激波DS的干扰互相耦合,形成多处超声速射流区域。因此,在内转式进气道额定工况下,前体激波入射V形溢流口驻点前亚声速区域可能带来复杂的流场结构,在实际应用中应当予以关注。
(3) 当斜激波入射在V形钝前缘弓形激波亚声速区的声速点附近时,出现一种有别于Edney第Ⅲ类激波干扰的波系结构,此时干扰产生的剪切层未入射壁面。
(4) 入射斜激波与V形钝前缘脱体激波DS、弓形激波BS的超声速区域干扰时的形态与Edney第Ⅰ、Ⅱ、Ⅵ类激波干扰类似,其中第Ⅱ类干扰形成的马赫杆后局部亚声速区域,可能造成内转式进气道局部不起动。
鉴于V形溢流口三维流动的复杂性,在激波干扰引起的气动热以及三维前体激波入射方面,还有待进一步的细致研究。
-
表 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.