二维舵面绕流的N-S方程数值计算研究

本文基于结构化网格,应用Baldwin-Lomax湍流模型进行绕舵面流动的二维N-S方程数值模拟,将部分计算结果与风洞试验结果进行了对比分析研究,分析了计算方法的可靠性和计算结果的合理性.为开展三维舵面绕流的N-S方程数值模拟和铰链力矩数值计算研究奠定了基础.     

基于非结构嵌套网格方法的旋翼地面效应数值模拟

 建立了一个基于非结构嵌套网格的流场求解器,用来精确模拟复杂的旋翼近地流场,为更好地分析地面效应(IGE)对旋翼气动特性的影响提供一套计算方法。在该求解器中,控制方程采用惯性坐标系下的非定常N-S方程,空间方向上采用二阶迎风格式,并用背景网格的一个面模拟地面作用,以方便地面边界条件的处理。应用所建立的模型,首先针对有实验结果可供对比的旋翼无地面效应（OGE）流场进行了数值模拟,以验证计算方法。然后着重计算了IGE下流场中的桨尖涡空间位置和旋翼拉力增益,并对旋翼在小速度前飞状态下的近地流场及地面涡形成过程进行了模拟和分析。在此基础上,得出了一些有意义的结论。     

翼型扑翼运动流场的数值模拟

基于非定常N-S方程,利用动网格方法和预处理方法对翼型在低速度下扑翼运动进行数值模拟。采用Jameson中心格式的有限体积法、双时间推进法求解非定常N-S方程;紊流模型采用基于SA紊流模型的DES模型。本文计算模拟翼型在低雷诺数下扑翼运动过程,计算所得推力系数结果与实验参考结果吻合较好,对扑翼飞行器设计有一定参考价值。     

多段翼型地面效应数值模拟与分析

通过数值模拟方法研究多段翼型的地面效应,采用有限体积法求解质量加权平均平均Navier-Stokes(N-S)方程,湍流模型采用Spalart-Allmaras(S-A)模型,采用滑移壁面模拟地面的相对运动.计算结果分析表明:随着飞行高度的降低,多段翼型的升力、阻力和低头力矩均减小;升力系数的减小主要是由于地面效应引起...     

二维带动力吸气式高超声速飞行器绕流的PNS-NS混合求解

为了提高吸气式高超声速飞行器绕流的求解效率,采用空间推进求解抛物化PNS(parabolizedNavier-Stokes)方程和时间迭代求解Navier-Stokes(N-S)方程的混合计算流体动力学(CFD)方法来求解高超声速飞行器整机绕流.在超声速占主导的流动区域采用空间推进求解抛物化N-S方程的方法,在亚声速和分离区采用时间迭代求解N-S方程的方法.对于求解二维带化学反应的吸气式高超声速飞行器绕流,混合CFD方法和完全时间迭代方法相比,可得到同等准确的数值模拟结果,并且求解效率提高了数倍.      

多段翼型地效风洞实验的地板边界层控制研究

采用有限体积法求解雷诺平均N-S方程和SSTk-ω两方程湍流模型,数值模拟二维多段翼型在不同风洞地板边界条件下流动,并对固定地板边界层吹吸气控制技术进行了数值研究。计算结果分析了地板边界条件对多段翼型气动力、表面压力分布及流动结构的影响,证明了吹吸气控制技术可以有效减小地板边界层厚度,并对均匀吸气与切向吹气两种方法的特点进行了比较。     

基于旋转爆震三维流场结构分析的计算模型对比研究

为了研究不同计算模型对三维旋转爆震数值模拟的影响,分别基于Euler方程、N-S方程、RANS方法和IDDES方法并结合滑移、无滑移壁面边界条件,耦合氢气/空气的有限化学反应速率模型（7组分8基元反应）,采用高分辨率的五阶有限差分格式WENO-PPM5离散对流项,对圆环筒型燃烧室内的三维旋转爆震波进行数值模拟,得到了不同模型下旋转爆震波的流场结构、传播特性和推力性能。对比了不同计算模型对流场结构的影响。当采用滑移壁面边界条件时,Euler方程和N-S方程的计算结果较为一致。当使用无滑移壁面边界条件时,边界层的存在会导致可燃混气与燃烧产物之间接触面上的爆燃燃烧区域沿壁面向上游渗透,增大爆燃区域的范围;且爆震波锋面非受限侧变窄、严重变形,不同计算方法计算的变形程度有所不同。     

使用高阶逆风通量差分裂格式的悬停旋翼流场数值模拟

为减少Jameson二阶中心差分有限体积法导致的旋翼尾迹的数值耗散,将三阶逆风格式(MUSCL)与通量差分裂方法相结合,建立了一个三维雷诺平均N-S方程数值模拟悬停状态旋翼流场的方法。为了充分考虑旋翼尾迹对流场的影响,采用周期性边界条件和由动量理论导出的远场边界条件。为进一步减少尾迹数值耗散和便于添加上述边界条件,采用了嵌套网格方法。然后,进行了算例计算,给出了桨叶表面的压强分布,与可得到的试验数据及二阶中心差分方法的计算结果进行了对比,并针对几种不同桨尖形状的旋翼悬停流场进行计算,数值结果显示:后掠桨尖可减弱超临界流动。此外,还计算和分析了旋翼下方不同位置上的涡量分布。计算结果表明,本文方法能够有效地减少旋翼尾迹的数值耗散。      

扑翼微飞行器飞行时的数值模拟与分析

采用数值方法求解三维非定常不可压缩的N-S方程,生成了相应的扑翼微飞行器的三维动网格;模拟了扑翼微飞行器在不同频率、不同攻角、不同振幅、不同飞行速度、不同的翅膀初始位置下的非稳态流场;计算了扑翼微飞行器在不同运动模态下的升阻力系数;分析了影响扑翼微飞行器空气动力学的特性的不同因素;为扑翼微飞行器的飞行机制和建造提供了理论依据。     

微型飞行器低雷诺数三维流场数值模拟

针对当前所研究的微型飞行器(MAV)流动特征,采用时间相关法和有限体积法求解N-S方程,选择了Van Albada插值限量因子的高分辨率空间差分格式,建立了模拟MAV三维低雷诺数流场的数值计算方法。通过熵校正技术改进和粘性计算准确性与网格密度关系的分析,确保了流场计算的准确可靠。使用该数值算法,对某常规气动布局的MAV进行流场模拟和分析,并比较了其流场与常规尺寸飞行器的不同。     

地面效应下前/后掠翼气动特性研究

针对飞行器在地效区飞行时复杂的流场特性,通过求解定常可压N-S方程,改变机翼后掠角和地效区飞行高度,研究不同前/后掠角机翼在地效区内的气动特性。结果表明:在地效区内,随着后掠角的增大,机翼的升力系数和阻力系数呈现先增后减的变化规律,后掠角在0°附近时升力系数达到最大值,阻力系数在10°附近达到最大值;俯仰力矩系数随着后掠角增加而减小;展向流动在后掠角为0°时最小,展向流动随着后掠角增大或减小剧烈变化;机翼下洗角随着后掠角增大而减小,随着离地高度的减小而减小。研究结果可为地效飞行器的概念方案设计和优化提供理论依据。     

基于结构动网格的无人机地面效应研究

为了更好地分析地面效应对无人机气动特性的影响,采用一种新的结构动网格生成技术,运用CFD手段对某无人机的地面效应进行了研究。在一套静态网格CFD计算基础上,运用结构动网格生成技术,获得非定常的网格移动,得到一系列的不同攻角、地面效应高度的CFD网格及初始流场。在此基础上进行定常流场求解,计算得到地面效应的定常CFD结果。最后对某双尾撑布局无人机的地面效应进行了CFD计算研究,对相应的全机气动特性的影响进行了分析。计算结果表明：地面效应对平尾的影响大于对机翼的影响并使平尾的失速迎角提前;该方法是可行的,具有一定的工程应用价值。     

多操纵面飞翼构型飞机舵面故障在线诊断方法

舵面故障会影响飞机的操纵性,为实现飞机的控制重构需在线诊断出舵面的故障信息.针对多操纵面飞翼构型飞机,提出了一种在线诊断其舵面故障的方法.飞翼构型飞机舵面故障在线诊断需辨识的操纵导数个数较多,直接使用最小二乘法难以精确辨识得到各操纵导数.基于限定记忆最小二乘法,采用分组辨识法对飞机的操纵导数进行辨识,减少了每次辨识操纵...     

近水面飞行气动特性数值计算方法研究

根据对水面波的数学描述，给出了其对流场影响的精确解析解。针对地(水)面效应飞行器的特点，分析了近水面波飞行时的流场特征，利用线性叠加原理，提出了飞行器与水面波耦合影响的非定常速度势修正项概念。分析了非定常气动特性的数值计算方法。对某实用型地效飞行器进行了计算，并对计算结果与试验结果进行了分析总结。     

三维多段机翼地面效应数值模拟

 通过数值模拟方法研究多段机翼的地面效应,采用有限体积法求解质量加权平均Navier-Stokes方程,湍流模型选用Spalart-Allmaras模型,利用运动壁面边界模拟地面的相对运动。计算结果分析表明:随着飞行高度的降低,多段机翼的升力、阻力和低头力矩均减小;迎角、展弦比越大,地面效应越明显,升力损失越大;升力的减小主要是由于地面效应导致机翼下方静压增大的气流通过缝隙进入机翼上表面流场,使得机翼下翼面压力的增加量小于上翼面吸力的减小量;地面效应使机翼上翼面翼尖容易发生分离;翼尖涡沿着展向方向向外移动,机翼诱导阻力减小。该文研究结果可以为大型飞机的增升装置地面效应设计提供参考依据。     

考虑地面效应的翼型动态特性数值模拟

为了研究考虑地面干扰的飞行器非定常气动特性,基于滑移网格技术,通过求解二维非定常Euler方程,对不同近地高度下的NACA0012翼型的动态振动特性进行了数值模拟。分析了存在地面干扰时的升力系数以及俯仰力矩系数迟滞环变化,并结合强迫振动法,计算了翼型各个近地高度下的纵向组合动导数。结果表明,地面效应不仅对定常流场产生影响,更显著地影响了非定常气动力及力矩,近地高度越小,这种效应越明显,使得升力系数的迟滞环面积变小,幅值减小,而俯仰力矩系数的迟滞环变化不规律,力矩系数的动导数随着近地高度减小而增大(绝对值减小),纵向的动态稳定性产生损失。地面效应干扰在飞行器非定常气动研究中应该引起重视。     

螺桨调速系统压力膜盒破损故障分析

 某型飞机螺桨调速系统中的压力膜盒曾产生破损事故，导致发动机空中停车，影响飞机正常飞行。对该螺桨调速系统进行了理论分析和实验研究。求取了系统各组成元件及管路结构数学模型，进行了系统的液压脉动响应数字仿真计算，并在飞机地面试车条件下对实际系统进行了压力脉动测量及频谱分析。计算所得的压力膜盒前压力脉动响应的峰－峰值，与飞机地面试车实测结果吻合。进而对压力膜盒破损事故的产生原因进行了仿真计算及分析，并提出消除故障的对策，经有关部门采用，取得了较好效果     

地面效应中悬停旋翼的自由尾迹计算

发展了一种具有良好收敛性的自由尾迹模型,用于对地面效应(IGE)状态下的旋翼悬停尾迹进行数值模拟。其中使用直涡元对桨尖涡丝进行离散,使用面元法模拟地面对流场的影响。针对以往类似研究中旋翼高度较低时自由尾迹迭代不易收敛的问题,引入了地下尾迹等比例修正并提出了环式地面面元分布等多项修改办法来改善计算模型的收敛性,同时这些改进也起到了削减计算量的作用。为了对新构建模型进行验证,以多个高度下的旋翼有地效状态为算例,对其尾迹几何结构进行了计算并与实验值进行了对比,结果显示二者符合得较好,证实了该模型具有较高的准确性。     

导弹地效区飞行空气动力和飞行力学特性研究

采用CFD和风洞试验技术,研究了导弹地效区飞行的空气动力和飞行力学相关特性。研究结果表明,导弹地效区飞行,由于受到地面效应的作用,其空气动力和飞行力学特性同常规自由空间飞行相比,发生了较大变化,常规导弹气动布局很难实现地效区飞行。由于地效区飞行运动响应十分复杂,加之飞行高度很低,这使现有飞行控制技术难以适应。这要求在气动布局设计中,必须将气动力、飞行力学和控制响应一体化考虑,气动布局技术是协调解决运动响应和动力学耦合的核心。