疏水/超疏水表面防/除冰原理及其研究进展

液体防冰、电热防冰、气热防冰、机械除冰等传统的防/除冰技术具有显著的防/除冰效果,但存在消耗能源大、防冰时间有限以及除冰不彻底等问题,而基于疏水/超疏水表面的防/除冰技术则具有独特的优势。本文首先阐述了利用疏水表面和超疏水表面防/除冰的原理及其影响因素;然后,总结了常用的疏水表面和超疏水表面在防/除冰中存在的问题,并指出利用液滴的反弹特性可更加有效地抑制结冰;最后,提出未来有必要对撞击液滴的动力学过程、传热过程和结冰过程的耦合特性进行系统研究,进而为利用液滴的反弹特性抑制结冰提供一定的理论指导。     

喷雾液滴撞击液膜影响参数及流动机理分析

液滴撞击壁面是喷雾冷却中最常见的现象。使用计算流体力学软件中的Volume of Fluid （VOF）模型对液滴撞击壁面过程进行了数值模拟研究，研究了不同参数的液滴撞击壁面液膜过程的水花形态特点，针对水花高度和直径等参数进行了分析，通过运动间断理论和无量纲参数Re和We的变化分析了不同参数对液膜的流动特征和水花形态的影响及水花产生机理。结果表明：大液滴撞击液膜时能够改善液膜的流动状态；过大的液滴速度会产生大量飞溅及表面干涸；液滴撞击薄液膜时，液膜的流动性较好，当0.6<h*<1.2时，液膜的流动性不随液膜厚度改变。     

超疏水表面的过冷水滴捕获规律研究

飞机在飞行过程中经常面临结冰危险,具有迫切的防冰需求,高效节能的超疏水电热复合防冰蒙皮具有广阔的应用前景,但目前设计流程中尚未考虑超疏水表面过冷水滴收集特性。为充分发挥超疏水表面对飞机防冰系统的作用,本文设计了一套液滴直径与流量可精准控制的喷雾装置,确定了过冷水滴捕获率试验的方案,并结合传热及表面特性相关理论总结了不同试验条件下表面捕获率变化规律。试验结果表明,在一定环境条件中,超疏水表面与聚酰亚胺（PI）表面的过冷水滴捕获率相对值保持稳定,约为25%。结合热传递及受力分析,超疏水表面能够有效减少水滴在表面的停留总量及停留时间,从而降低过冷水滴捕获率。本文为新一代超疏水电热复合蒙皮的防冰功率精确设计与能耗优化提供了依据,对保证恶劣环境下无人机正常工作具有重要意义。     

柱坐标系下准三维液滴撞击壁面的数值分析

液滴撞击壁面现象广泛存在于自然界与工程领域。为研究其运动特性,本文基于柱坐标系与VOF模型,对准三维液滴撞击固体壁面过程进行模拟。该方法实现了三维到二维的合理简化,模拟与实验结果较为吻合,铺展系数约为2.1。在此基础上,本文对单液滴撞击固体壁面过程进行了多工况模拟,模拟表明液滴在速度和直径较大情况下铺展系数较大。相反,在表面张力系数和接触角较大情况下铺展系数较小。最后,对双液滴相继撞击壁面的过程进行了模拟,发现液滴在铺展过程中出现了二次铺展现象,并探究了速度和液滴中心间距对铺展特性的影响。结果表明,双液滴在速度和液滴间距较大的情况下铺展系数较大。     

织构化表面对异质金属润湿性及界面反应的影响

不锈钢与铝合金广泛应用于火箭贮箱中隧道管与补偿器的焊接，熔钎焊已成为其主要连接方式，但接头存在低应力开裂，容易引起燃料泄漏等问题。为提高接头可靠性，首先采用超快激光加工技术在不锈钢表面烧蚀制备了微槽和微坑两种典型的织构化表面。然后，通过铝合金在不锈钢表面的原位润湿铺展试验研究了两种微织构对润湿铺展方向、接触角和界面反应产物的影响。试验结果表明铝合金具有沿微槽加工方向润湿铺展的趋势，而微坑的铺展则接近圆形；在微槽表面表现出更小的接触角，同时基于Wenzel方程对在两种表面润湿铺展的接触角进行数值计算，计算结果略大于实际结果。最后，分析了界面冶金反应对润湿铺展的影响。结果表明界面反应对铝合金的润湿铺展具有一定的促进作用，但过于粗大的界面反应产物则会阻碍润湿铺展。该激光织构化的方法不仅可用于不锈钢与铝合金体系，还可用于航空航天领域中其他金属与金属、金属与非金属之间的连接。     

幂律型流体雾化SPH方法数值分析

幂律型流体的雾化过程存在复杂的界面运动,用传统网格法很难精确追踪运动界面.为研究幂律型流体的雾化特性,运用SPH方法对典型的双股幂律型流体撞击雾化问题进行三维数值分析.根据文献实验,采用SPH方法模拟获得与实验条件相应的数值结果,对比后表明,二者雾化角相当吻合,数值结果还成功捕捉到液膜向液丝的破碎过程及网状的液丝分布状态,验证了方法的有效性.分析了射流速度和撞击角度对雾化角的影响,得到雾化角随着射流速度和撞击角度的增加而增大;对雾化后的速度场进行数值分析后表明,撞击点附近,惯性力的作用使速度场变化剧烈;在流体远离撞击点的过程中,粘性耗散作用使速度场趋于稳定,但速度大小小于初始撞击速度.     

飞秒脉冲激光诱导TC4微结构表面抑冰特性实验

飞机发动机进气道前缘唇口积冰将会严重威胁航空安全,仿生研究表明具有微纳结构的疏水表面可以起到良好的抑冰效果。针对飞机唇口材料TC4,采用飞秒脉冲激光诱导制备TC4微结构表面,利用三维形貌仪和扫描电镜对TC4合金表面三维形貌和微纳结构进行观测,应用接触角测量仪分析表面浸润改性,依托结冰特性实验系统测试微结构表面抑冰抑霜性能,并分析飞秒脉冲激光加工工艺参数对表面微观结构和抑霜特性的影响机制。研究结果表明：随着激光扫描速度的增大,TC4合金表面形成的拱形沟壑深度增加,沟壑上方出现干涉条纹以及圆形凸起且微纳凸起的尺寸随扫描速度的增大而增大,接触角先减小后增大再减小;加工后表面液滴冻结时间比未加工表面延迟30 s;扫描速度2 000 mm/s时的液滴冻结时间最长,霜层质量最小,高度最低。飞秒激光加工TC4合金表面形成的微纳结构以及表面吸附的有机物能够改变表面接触角;粗糙度和表面形貌能够影响表面结冰时间和结霜量。     

气泡和液滴在上升或下降时大变形的数值模拟

本文在考虑界面张力的条件下用涡片方法导出了二维无粘气泡和液滴在另一种静止无粘流体中上升或下降时运动和变形的支配方程，提出了求解支配方程的数值方法，并在数值模拟的例子中研究了空气泡和煤油滴在水中上升时及水滴在煤油中下降时的变形规律，研究了数值模拟中界面张力对界面形状的影响，气泡上升的数值模拟结果与相应的试验结果基本一致。数值模拟结果表明本方法可正确模拟饣泡和液滴在上升或下降过程中的大变形，直到界面发     

低速下水滴撞击固体表面运动模拟

以水滴撞击飞机部件造成结冰问题为背景,运用level set界面追踪方法,把不可压黏性方程组中的表面张力用level set函数表示,并利用2阶变密度投影算法求解二维不可压Navier-Stokes方程组,数值模拟了低速下的水气界面追踪问题.为了能够处理密度比为1000的水气界面问题,在界面处引入了一定宽度的光滑过渡带宽,进而消除了由于界面两边大密度比引起的数值震荡,由此较为成功地模拟了大水滴在不同表面张力时撞击平板后的形态变化及其运动轨迹.数值实验结果表明本算法可以有效地处理界面处密度比和黏性很大的不可压多相流问题.考虑界面处的表面张力,本算法可模拟不同表面张力下的界面运动变化.      

飞秒激光制备铜微纳结构表面的润湿及抗结冰特性研究

基于微结构表面形成的疏水特性可以有效应用于水腐蚀防护、电缆防护、飞机防冰等重要领域。研究中,不借助其他化学处理手段,仅通过飞秒激光直接制备了多种疏水功能特性微纳结构铜表面,并系统研究了温度对微纳结构表面的疏水特性影响及抗结冰性能。结果表明,这些疏水微纳结构表面在低温下与原表面相比都能有效延缓水滴的结冰过程,这一现象的产生与其表面润湿性能及微纳结构形貌有关。     

ν-gap度量及其在飞行控制律评估中的应用

 传统控制律评估方法主要用于单输入单输出（SISO）系统，且对模型参数摄动考虑不够全面,针对这些不足，研究了ν-gap度量方法。在介绍系统广义稳定裕度相关概念的基础上，给出了ν-gap度量的定义、特点和性质以及近似摄动模型的计算，提出ν-gap度量评估控制律的步骤。实例结果表明，该方法不仅克服了上述传统评估方法的缺陷，而且还有根据所求的各摄动影响情况忽略影响小的元素，以减少计算量及可以找到最坏情况下的参数摄动组合等优点。     

Integrated Entry Guidance for Reusable Launch Vehicle

A method for the implementation of integrated three-degree-of-freedom constrained entry guidance for reusable launch vehicle is presented. Given any feasible entry conditions, terminal area energy management interface conditions, and the reference trajectory generated onboard then, the method can generate a longitudinal guidance profile rapidly, featuring linear quadratic regular method and a proportional-integral-derivative tracking law with time-varying gains, which satisfies all the entry corridor constraints and meets the requirements with high precision. Afterwards, by utilizing special features of crossrange parameter, establishing bank-reversal corridor, and determining bank-reversals according to dynamically adjusted method, the algorithm enables the lateral entry guidance system to fly a wide range of missions and provides reliable and good performance in the presence of significant aerodynamic modeling uncertainty. Fast trajectory guidance profiles and simulations with a reusable launch vehicle model for various missions and aerodynamic uncertain-ties are presented to demonstrate the capacity and reliability of this method.     

基于飞机油箱模型形状特征油量测量切片步长选择方法研究

 在分析飞机数字式油量测量过程中目前广泛使用的切片法油量测量原理的基础上,针对现有的定步长切片法无法得到准确、可靠的燃油质量特性数据库的缺陷,结合对飞机油箱模型形状特征的分析,提出了基于飞机油箱模型形状特征的油量测量切片步长选择方法。此方法包括切片步长整体和局部选择两个过程,整体选择以实现相邻两切片平面所夹油箱模型体积近似相等为目的来确定切片步长,以体现油箱模型截面整体变化规律;局部选择以设计切片平面与截面突变平面重合或尽可能接近的方式,突出油箱截面的局部变化特征。实验结果表明：该切片步长选择方法较定步长方法能够建立更为合理、可靠的燃油质量特性数据库,从而提高了油量测量精度。     

The effect of inlet conditions on the flow and heat transfer in multiple rotating cavity with axial throughflow

This paper discusses experimental results from two different build configurations of a heated multiple rotating cavity test rig.Measurements of heat transfer from the discs and tangential velocities are presented.The test rig is a 70% full scale version of a high pressure compressor stack of an axial gas turbine engine.Of particular interest are the internal cylindrical cavities formed by adjacent discs and the interaction of these with a central axial throughflow of cooling air.Tests were carried out for a range of non-dimensional parameters representative of high pressure compressor internal air system flows(Re up to 5×106 and Rez up to 2×105).Two different builds have been tested.The most significant difference between these two build configurations is the size of the annular gap between the(non-rotating) drive shaft and the bores of the discs.The heat transfer data were obtained from thermocouple measurements of surface temperature and a conduction solution method.The velocity measurements were made using a two component,LDA system.The heat transfer results from the discs show differences between the two builds.This is attributed to the wider annular gap allowing more of the throughflow to penetrate into the cavity.There are also significant differences between the radial distributions of tangential velocity in the two builds of the test rig.For the narrow annular gap,there is an increase of non-dimensional tangential velocity V/Ωr with radial location to solid body rotation V/Ωr=1.For the wider annular gap,the non-dimensional velocities show a decrease with radial location to solid body rotation.      

二维翼段颤振的μ控制

 采用超声电机作为作动器来实现含控制面的翼段颤振鲁棒抑制。针对作者设计的二维翼段颤振主动抑制系统，通过理论与实验相结合的方法，建立了考虑沉浮方向阻尼和作动器模型参数不确定性的控制系统模型，设计了μ控制器，并对控制器做了降阶处理。数值仿真和风洞试验表明，μ控制器可有效地抑制颤振的发生，将颤振临界速度提高23.4%。相对于H_∞控制器，μ控制器的控制效果和鲁棒性更好。     

航空发动机管路测量数据分割方法

 为了提高航空发动机管路测量数据的反求建模效率,提出了一种区域增长分割算法。该方法主要是利用管路表面在曲面索引系数映射、主曲率映射、高斯映射上的特性,并基于均值漂移算法和遗传算法提出了区域分割算法中种子区域的选择策略。然后利用区域增长分割方法实现了对管路测量数据的分割。经过仿真和实测数据验证,所提管路分割算法具有较好的分割质量和效率。     

Abnormal Shape Mould Winding

为解决网格化芯模的缠绕问题,本文提出了复合材料面片缠绕机理;接着详细分析了面片缠绕过程中的芯模凹曲面上纤维滑线和架空现象,应用微分几何曲面理论和空间几何理论,提出判据及其解决方案;最后,针对飞机发动机进气道的缠绕成型,编制缠绕控制程序并进行相应的实验,验证了面片缠绕方法的可行性。     

Interfacial Properties Modification of Carbon Fiber/ Polyarylacetylene Composites

This work was dedicated to performing surface oxidation and coating treatments on carbon fibers (CF) and investigating the changes of fiber surface properties after these treatments, including surface composition, relative volume of functional groups, and sur- face topography with X-ray photoelectron spectroscopy (XPS) and atom force microscopy (AFM) technology. The results show that, after oxidation treatments, interfacial properties between CF and non-polar polyarylacetylene (PAA) resin are remarkably modified by removing weak surface layers and increasing fiber surface roughness. Coating treatment by high char phenolic resin solution after oxida- tion makes interface of CF/PAA composites to be upgraded and the interfacial properties further bettered.     

Thermal Stresses and Cracks During the Growth of Large-sized Sapphire with SAPMAC Method

本文利用有限元法对冷心放肩微量提拉法(SAPMAC)生长大尺寸蓝宝石单晶过程中的热应力进行了模拟计算,应用临界缺陷理论解释了裂纹的萌生与扩展机理。研究表明,晶体内的热应力主要与晶体生长速率,环境温度和结晶取向有关;较大热应力多出现在结晶界面,放肩、收尾以及直径急剧变化等位置;最大热应力总是出现在籽晶与新生晶体的界面附近。裂纹将在临界缺陷位置产生,并在应力作用下沿a或m面扩展。计算结果与实验结果基本吻合,通过晶体生长系统和生长工艺的改进,对晶体的开裂问题得到了有效的抑制。     

Optimization Design System for Composite Structures Based on Grid Technology

To solve the topology optimization of complicated multi-objective continuous/discrete design variables in aircmit structure design, a Parallel Pareto Genetic Algorithm （PPGA） is presented based on grid platform in this paper. In the algorithm, the commercial finite element analysis （FEA） software is integrated as the calculating tool for analyzing the objective functions and the filter of Pareto solution set based on weight information is introduced to deal with the relationships among all objectives. Grid technology is utilized in PPGA to realize the distributed computations and the user interface is developed to realize the job submission and job management locally/remotely. Taking the aero-elastic tailoring of a composite wing for optimization as an example, a set of Pareto solutions are obtained for the decision-maker. The numerical results show that the aileron reversal problem can be solved by adding the limited skin weight in this system. The algorithm can be used to solve complicated topology optimization for composite structures in engineering and the computation efficiency can be improved greatly by using the grid platform that aggregates numerous idle resources.