直升机的气动弹性问题

直升机的气动弹性问题与固定翼飞机不同,不仅要考虑单片桨叶,更要将旋翼视为一个整体,考虑其动态入流、尾迹影响以及旋翼与机身之间的相互耦合等。就单片桨叶而言,在结构动力学上,需要考虑离心力场、几何非线性以及桨叶的非线性挥舞-摆振-扭转耦合;在气动力上,需要考虑动态入流以及桨尖处可能的失速效应,本质上属于非线性气动弹性力学范畴。由于旋翼气动力通常是以周期形式通过旋翼轴传给机身,并引起机身振动,而机身运动又通过改变桨叶根部形态反过来影响旋翼的气动弹性特性,这种旋翼/机身耦合问题,也是近年来直升机气动弹性问题研究中的重要方向和热点之一。此外,随着旋翼流场数值分析方法的日趋成熟,采用动态重叠网格或滑移网格方法来实现桨叶运动,并通过动网格技术来实现桨叶的弹性变形,从而实现弹性旋翼流场的数值模拟,目前正呈现出勃勃生机,成为直升机气动弹性研究的又一重要方向和热点。随着各种新构型直升机的相继出现,如倾转旋翼机、前行桨叶概念旋翼(ABC)直升机和复合式直升机等,也带来了新的气动弹性问题。不断发现问题、解决问题,推动本学科持续发展,永远是气动弹性工作者终身奋斗的目标。     

二元双喉道射流推力矢量喷管流动参数影响的数值研究

采用数值模拟方法研究了不同流动参数对二元双喉道射流推力矢量喷管(Dual-throat fluidic Thrust-vectoring Nozzle,DTN)内流特性和推力矢量控制效果的影响。结果表明,DTN在非推力矢量情况下,NPR在3~4范围时,推力系数较大,达到0.968,而流量系数较小,仅为0.93;NPR再增大,推力系数迅速下降。在推力矢量情况下,落压比一定时,随着次流流量比的增加,推力矢量角增加,而流量系数、推力系数、推力矢量效率减小;次流流量比一定时,随着落压比的增加,推力矢量角减小,系统推力系数先增加后减小,流量系数略微增加。     

柔性扑翼气动结构耦合特性数值研究

微型扑翼体积小、重量轻,其柔性变形对气动特性有显著的影响。通过求解雷诺平均N-S方程(ReynoldsAveraged Navier-Stokes,RANS)和结构动力学方程,对微型柔性扑翼飞行器的气动结构耦合特性进行了数值模拟研究。针对微型扑翼的大幅运动,发展了适用于扑翼的气动结构耦合数值计算方法,研究了微型扑翼的气动结构耦合特性。通过求解雷诺平均Navier-Stokes(RANS)方程得到微型扑翼的非定常气动特性;利用哈密顿原理(Hamilton Principle)推导了扑翼的结构动力学方程,采用结构有限元方法对该动力学方程进行离散并求解,得到扑翼的动态结构特性;采用松耦合方法进行迭代。计算结果与风洞实验结果相比吻合良好,验证了所发展方法的有效性。在此基础上研究了惯性力和关键运动参数对柔性扑翼气动及结构特性的影响规律,有助于比较详细、全面地了解微型扑翼的气动机理,为柔性扑翼的设计提供了参考依据。     

面向未知区域深度测量的序列图像稠密点特征生成算法

对未知着降区平坦度测量是无人机在复杂地形下安全着陆的关键问题。首先,根据小孔成像原理推导出基于单目序列图像的未知区域深度计算方程;其次,针对稀疏匹配存在深度信息重构误差大而稠密匹配在平滑区域误匹配率高的问题,提出一种基于Delaunay三角剖分的稠密点特征生成算法;然后,分别对序列图像中的2帧图像提取亚像素级Harris角点和尺度不变特征变换(SIFT)特征点,并分别进行特征点匹配;再以2种特征点间的欧氏距离作为约束条件将2种特征点进行融合,生成准稠密特征点;最后,将准稠密特征点进行Delaunay三角剖分,并根据每个剖分三角形上3个顶点像素偏差的方差值制定稠密特征点的生成策略,并结合所提出的深度计算方程计算整个未知区域各点的深度信息。通过Vega Prime(VP)搭建仿真演示验证系统,实验结果表明在机载相机距地面400m处计算高度分别为90m和55m的物体深度信息时,其深度测量相对误差不超过0.89%,具有较高的精度。     

(高)超声速流动试验技术及研究进展

近年来,与高速飞行器相关的(高)超声速流动受到了极大的关注。这类流动所具有的非定常性、强梯度和可压缩性对试验方法和风洞设计技术提出了挑战。超声速纳米示踪平面激光散射(NPLS)技术是由作者所在团队研发的非接触光学测试技术。它能够以较高的空间分辨率来揭示超声速三维流场的一个瞬态剖面的时间解析的流动结构。介绍了NPLS技术以及基于NPLS开发的密度场测量、雷诺应力测量和气动光学波前测量等方法,并回顾了这些技术在超声速边界层、超声速混合层、超声速压缩拐角、激波/边界层相互作用和光学头罩绕流等流动中的应用,清晰地再现了边界层、混合层、激波等典型流场结构及其时空演化特性。另外,为了模拟和研究高空大气条件下边界层自然转捩和超声速混合层的转捩特性,介绍了高超声速静风洞、超-超混合层风洞的设计技术以及层流化喷管的设计方法。     

一种孔径和频率二维稀疏的步进频SAR成像方法

 步进频率信号(SFWs)在不增加雷达系统瞬时带宽的情况下能够获得高的距离向分辨率的同时,也存在着抗干扰能力较差及其等效重复频率较低的问题,并且在方位向积累时间内由于雷达载机工作状态的变化,会导致方位向的数据录取不完整。针对上述问题,提出一种孔径和频率二维稀疏的步进频合成孔径雷达(SAR)成像方法。首先,分析了稀疏步进频率信号(SSFWs)的SAR成像模型,然后基于压缩感知理论完成距离向成像处理。其次,针对稀疏孔径的回波数据,通过构造成像算子和压缩感知重建模型的方法实现其距离徙动校正和方位压缩处理,进而获得二维成像结果。相比于传统的步进频率信号SAR成像,利用所提方法能够在少量的频率资源和雷达回波数据情况下实现准确的SAR成像。最后,通过对仿真和实测的步进频率雷达数据进行成像处理,验证了所提方法的有效性和可行性。     

运动油滴/固体壁面斜碰撞的状态辨识及特征分析

 航空发动机轴承腔油气两相流动与换热研究是轴承腔润滑与换热设计的重要基础工作,而运动油滴与轴承腔壁面的碰撞状态辨识及特征分析是轴承腔油气两相流动与换热研究的组成部分。基于VOF方法建立了运动油滴与固体壁面斜碰撞数值仿真模型,实现了运动油滴与固体壁面斜碰撞状态的数值模拟,分析了碰撞油滴直径、入射角和碰撞速度对沉积油膜铺展长度、油膜堆积厚度和溅射油滴数目等碰撞状态特征参数的影响规律,提出了运动油滴与固体壁面斜碰撞条件下的状态判断准则,并得到了文献的物理试验工作的支持。研究结果表明:随着油滴入射角减小,沉积油膜铺展长度和油膜堆积厚度增大,溅射油滴数目减少;随着碰撞速度和油滴直径的增大,沉积油膜铺展长度和溅射油滴数目均增大,但油膜堆积厚度对应前者呈减小趋势、对应后者仍呈增大趋势。     

重型燃气轮机燃烧过程的反应动力学数值模拟

为了阐明重型燃气轮机燃烧过程的反应动力学特性,采用Gri_3.0,NUI_Galway与USC_2.0动力学模型对甲烷燃料在定容燃烧反应器中的燃烧特性进行了数值计算,并与实验结果进行了对比分析,确定了甲烷燃料的反应动力学模型;耦合该反应动力学模型与CFD计算软件,对燃用甲烷燃料的重型燃气轮机燃烧室单个火焰筒的燃烧过程进行了反应动力学分析,并与相应实验结果进行了对比分析。结果表明:与NUI_Galway及USC_2.0动力学模型相比,采用Gri_3.0动力学模型计算得到的甲烷燃料多工况下的燃烧特性与实验值吻合较好;在该重型燃气轮机燃烧室单个火焰筒中,在火焰筒头部与主燃区同时存在两个明显的呈对称状态的回流区;在该火焰筒中的高温区域,O,OH与H等活性组分以及CO2的摩尔分数达到最大,而在低温区域CO的摩尔分数达到最大;同时,与实验值相比,采用Gri_3.0动力学模型计算得到的单个火焰筒的出口平均温度略高约4K,热点温度高约197K。     

燃气涡轮无气膜冷却动叶参数化设计系统应用

为进行燃气涡轮冷却结构的设计,在考虑总的温度分布系数(OTDF)及不考虑OTDF两种情况下,采用一套设计方法,完成了燃气涡轮第一级动叶冷却结构的整体设计。通过设计表明:采用管网计算,并通过三维导热计算进行热分析,最后通过气热耦合计算能够快速地设计出较佳的冷却结构;不考虑OTDF设计时,第一腔流量为16.9g/s,第二腔流量为40.8g/s,最大无量纲温度0.700,采用双进口蛇形通道,换热效果较佳,且结构设计较为简单;考虑OTDF时,为达到设计要求,第一腔流量为18.2g/s,第二腔流量为25.4g/s,第三腔流量为5.3g/s,最大无量纲温度0.752。通过多方案设计,得出在无气膜情况下,采用三个冷气进口多回转通道能够达到较好的冷却效果。     

乙醇-空气稀液雾值班火焰的大涡模拟研究

为加深对气液两相湍流燃烧现象的认知,检测火焰面模型在液雾燃烧中的适用性,本文在欧拉-拉格朗日架构下使用火焰面/反应进度变量模型(FPV)数值模拟了乙醇-空气稀液雾值班火焰。欧拉坐标系下的气相湍流场使用大涡模拟方法模拟,离散液相则使用拉格朗日颗粒轨道模型进行描述,考虑了相间质量、动量、能量交换。燃烧模型中采用碳元素定义混合物分数,在混合物分数方程源项中体现液相对燃烧模型的影响。模拟得到的气相温度分布和液相统计值均和实验数据较好吻合,验证了该燃烧模型对稀液雾扩散类型火焰的适用性。分析瞬时图发现,该稀液雾火焰的最高燃烧温度往往位于当量混合物分数附近,在出口下游20倍直径处火焰完全点着,此处上游FPV模型能给出局部点火熄火现象。蒸发作用在剪切层和点火区域较强,而液滴和火焰作用较弱,单一液滴很少被火焰包围。