翼型近波浪水面气动特性研究

采用求解Navier-Stokes方程的数值方法研究了翼型NACA4412近距离经过波浪水面时的气动特性。对数值方法的准确性进行了验证。计算了翼型经过波浪水面和固壁波浪地面2种边界条件下的气动力系数，并进行了对比。研究结果表明:翼型在经过波浪水面时气动力系数会发生周期性的变化，与固壁波浪的情况相比，气动力的变化曲线存在显著差异，波动幅度更大。通过对流场结构的分析，发现了翼型和波浪水面之间的作用机理。波浪水面的质点存在竖直方向上的运动，在小地面间隙时，水面质点向上运动会挤压翼型和水面之间的空气，从而造成翼型气动力大幅波动。同时解释了来流速度越大，气动力系数波动幅度减小的原因。      

不同厚度翼型动态失速涡运动数值研究

在低马赫数下, 对三种不同相对厚度的NACA系列基本厚度翼型在俯仰振荡运动中的动态失速现象进行了数值研究.数值模拟时采用双时间法和LU-SGS隐式解法相结合的模式求解了非惯性坐标系下的非定常纳维-司托克斯(NS)方程组, 空间离散采用Roe格式,并结合刚性动网格生成技术.采用全湍流计算,通过引入BL(Baldwin-Lomax)模型计入湍流影响.计算出的气动力迟滞曲线与实验结果变化趋势符合较好,表明了数值方法的有效性,同时通过比较和分析不同厚度翼型在轻失速和深失速下的流场结构,发现绕翼型的失速旋涡产生和发展规律是明显不同的.      

基于试验气动力的纵向机动飞行载荷分析

提出了一种基于非线性风洞试验数据(简称试验气动力)同步进行纵向机动飞行过程和飞行载荷的计算分析方法.通过曲面样条插值将试验气动力的升力系数和俯仰力矩系数引入纵向机动飞行过程的计算,同时基于试验气动力的压力分布进行全机飞行载荷的静气动弹性修正.弹性载荷的计算采用线性气动力影响系数矩阵处理.采用偶极子格网法提供气动力(简称理论气动力).将基于这两种气动力的计算结果进行了对比.结果表明:基于试验气动力的机动飞行过程更为合理;基于理论气动力的载荷随机动飞行过程的变化趋势有可能与基于试验气动力的载荷变化趋势不一致.      

环量控制翼型非定常气动力建模

针对目前环量控制技术中射流参数与迎角对翼型气动特性的影响高度耦合，对应非定常气动力模型精度较差的研究现状，基于环量控制翼型强迫俯仰振动数值模拟数据，借助Kriging模型实现环量控制翼型的定常气动力插值，借助微分方程模型完成了适用于环量控制翼型的线性微分方程建模，采用两步线性回归参数辨识方法辨识线性微分方程模型中特征时间常数等参数，对高动量系数大振幅流动状态下的非线性影响进行修正。研究结果表明:基于Kriging模型实现的环量控制翼型定常气动力插值精度较传统气动导数模型高，建立的环量控制翼型非定常气动力模型能够精确预测不同流动状态下的气动力和力矩系数变化情况。      

采用欧拉两相流法对翼型表面霜冰的数值模拟

基于欧拉两相流理论提出一种翼型表面霜冰的数值模拟方法.采用同位网格上的SIMPLE (Semi-Implicit Method for Pressure Linked Equations) 算法求解空气控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击,通过求解过冷水滴的控制方程得到翼型表面的水滴收集特性;利用冰层时间推进法模拟积冰过程,使用积冰法向生长假设生成积冰外形.对NACA 0012翼型在0°和4°迎角下的积冰情况进行了研究,冰形预测结果与实验结果一致.通过对结冰翼面的压力分布进行分析,表明积冰对翼型气动性能会产生不利影响.     

飞行器表面三维流场与固壁温度场的耦合分析

分析了飞行器表面温度场形成的物理机制,建立了飞行器蒙皮温度场耦合计算的理论模型,较为完整的考虑了各种热源对飞行器蒙皮温度的作用.采用有限体积法和重整化群RNG(Renormalization Group)k-ε湍流模型对流场求解N-S方程,蒙皮固壁中采用热传导模型,采用离散坐标法DOM(Discrete Ordinates Method)计算蒙皮固体壁面对外界大气的辐射.通过数值计算研究了飞行高度、马赫数和内热源对蒙皮温度分布的影响,以及高度及马赫数对飞行器表面气动加热的作用.结果表明,该方法适合考虑飞行器外部流场、固壁等多种热源作用的温度场数值模拟.     

跨音速翼型的气动最优化设计方法

将计算流体动力学(CFD)与最优化技术相结合,通过数值求解欧拉方程,对翼型绕流流场作出了数值模拟.再利用几何、流动和最优化控制方程,反复迭代求得在一定约束条件下气动性能最优的翼型.本文以NACA0012为原始翼型,选取两种设计工况,都取得了满意的结果.     

二维翼型绕流的CE/SE方法

采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础.      

基于降阶模型的翼型结冰冰形预测方法

翼型结冰冰形的数值模拟预测通常比较复杂耗时,为了更加快速准确地预测冰形以减少计算资源消耗,建立了基于本征正交分解（POD）和Kriging模型的冰形快速预测方法。利用CFD数值模拟结果来构建样本空间,以飞行迎角为例详述了降阶模型的冰形预测的实现手段,并结合试验设计方法,完成了多参数的结冰冰形快速预测,同时研究了先进的Blind-Kriging模型的相关方法以及对于预测结果的改进。结果表明,降阶模型预测翼型结冰冰形与CFD数值模拟结果吻合较好,表明降阶模型可以快速、精确地应用于翼型结冰冰形预测。      

某型无人机全机气动力数值计算

使用三维造型软件和计算流体力学软件对某型无人机进行了气动力数值模拟,得到小迎角范围内的升、阻力系数和力矩系数以及无人机表面的压力、密度和来流速度分布情况,同时,对该无人机进行了气动力的工程经验估算,并将估算结果与数值计算结果进行对比分析.研究表明,数值计算能够较好的模拟小尺度低速流场流动,对无人机的研制能起到重要的作用.计算结果为下一步无人机的气动外形优化设计提供了理论依据和重要参考.     

基于NACA0030的波纹状翼型气动特性探索

相对于光滑翼型，波纹状翼型的气动特性呈现出一些独特现象。为了深入探索这种布局的气动特点，在前期风洞试验的基础上，以NACA0030翼型为基础，设计了一组具有不同外形特征的波纹状翼型，开展了非定常数值模拟工作，详细研究了低雷诺数（Re=12×104）流动情况下波纹状外形对流场涡流结构和总体气动特性的影响规律。计算结果表明:相对于光滑翼型，波纹状翼型流动的分离流现象更明显，升力和升力线斜率有明显下降，但推迟了失速现象。波纹状翼型表面越光顺，气动特性越接近于光滑翼型。虽然波纹状翼型的压差阻力大于光滑翼型，但是波纹状外形产生的回流可以减小摩擦阻力。      

积冰几何特性对翼型性能影响的神经网络预测

积冰几何形状对翼型气动系数的影响是复杂的.采用BP(Back Propagation)神经网络的LM(Levenberg-Marguardt)学习算法,建立明冰的典型几何特性(冰角前缘半径、冰角高度和冰角位置)对翼型气动系数影响的神经网络,得到该3种几何参数对气动系数影响的规律;建立了典型冰形参数对最大升力系数影响的神经网络,该网络能很好的预测冰形参数对应的最大升力系数值;此外,建立了冰型位置对舵面铰链力矩系数影响的神经网络.仿真结果表明,BP神经网络仿真结果与实验值具有高度一致性,并能预测非实验值条件下的气动系数;翼型表面积冰位置变化对气动系数影响最大;铰链力矩系数在失速迎角达到之前就发生突变,可以更安全地用来预测失速的发生.      

连续变弯度翼型动态气动特性数值模拟

针对翼型后缘连续变弯度运动中后缘边界精确数值模拟的问题，提出了基于二维多项式的时空曲面拟合方法，实现了对后缘边界时空位置的精确模拟。在此基础上，基于OpenFOAM发展了翼型后缘连续变弯度与大幅度俯仰运动耦合的运动边界数值模拟，并计算了翼型耦合运动的气动力，讨论了后缘线性/非线性变形对翼型大迎角动态气动特性的影响规律。结果表明:后缘运动对翼型俯仰运动的升阻特性影响显著，特别在翼型大幅度俯仰时后缘非线性变形对升阻特性改善效果比线性变形大6%~10%。同时还研究了翼型俯仰与后缘变形运动相位差对气动特性的影响。特别地，当相位差为180°时，后缘运动使动态失速时的最大升力提高50.3%，平均升力提高34.6%；当2种运动相位差为0°时，后缘运动使动态失速时的最大阻力降低39.7%，平均阻力降低30.2%，最大升阻比提高22.3%，平均升阻比提高16.8%；同时，翼型在动态俯仰过程中出现负阻力现象，对产生负阻力的原因进行了分析。这些结果可用于指导连续变弯度后缘控制律的设计。      

BLI效应下整流罩设计对翼型气动特性的影响

边界层吸入(BLI)效应对飞行器气动特性的影响比较显著,而整流罩的设计会进一步影响BLI效应下的翼型气动特性。为了揭示BLI效应下整流罩的主要设计参数对翼型气动特性的影响及其原因,本文采用计算流体力学(CFD)和Morris敏感度分析相结合的方法对该问题进行了详细研究,得到了整流罩主要设计参数对翼型气动特性的敏感度排序和耦合影响程度排序;对敏感度较高和耦合影响较大的参数进行了流动分析。结果表明:在巡航和起飞2种状态下,对气动系数影响相对较大的设计参数是整流罩最大厚度和进气边界弦向位置,整流罩最大厚度对翼型气动特性影响的主要原因是整流罩背风面会发生局部分离,且其还会改变阻力-流量系数曲线的趋势;整流罩最大厚度和进气边界弦向位置对翼型气动特性的耦合影响作用较强。      

机翼防冰过程中冰脊问题的数值分析

对机翼防冰过程中冰脊的形成特点及冰脊对机翼气动特性的影响进行了计算分析.基于经典的Messinger结冰模型开发了多步结冰程序,对不同条件下过冷水滴的撞击结冰进行了热质耦合计算,计算结果和文献试验结果吻合较好,表明该热质耦合算法的正确性.在此基础上对不同环境温度、飞行速度和不同加热功率等条件下冰脊的生长特点和机翼的气动特性进行了计算分析.结果表明:在非霜冰条件下,冰脊主要在热防护极限外并紧挨着热防护极限的位置处形成和发展,而在霜冰条件下,冰脊主要在机翼下表面形成,但是在热防护区域内有显著的结冰出现,该条件下的结冰对机翼的气动特性具有较大的破坏性.     

鸭式旋翼/机翼飞行器转换末段气动特性

采用数值模拟方法研究鸭式旋翼/机翼(CRW,Canard Rotor/Wing)飞行器在转换过程末段,旋翼转速极低时全机气动特性变化规律及其产生原因.给出了旋翼旋转一周时,全机气动力、气动力矩、焦点位置变化规律,对此布局形式,转换过程末段全机升力、阻力变化幅度可达10.7%,3.7%,焦点可移动0.6 m.研究显示:旋翼处于前后不对称流场及旋翼处于不同方位角时对机体的不对称干扰是气动力与气动力矩变化原因,旋翼与平尾升力线斜率变化、旋翼自身焦点位置变化导致了全机焦点移动.