后掠机翼横流驻波及其谐波研究

使用Mack方法求解三维不可压Orr-Sommerfeld方程（忽略曲率影响）,得到模型在风洞最佳实验条件下空间最大放大率展向波长为4mm。以此为基础,选择大于此波长的具有代表性的展向扰动波长6mm与7mm,计算了这两个驻波及其谐波空间的放大率。研究在整个层流区主导流动的扰动波长及其发展趋势。使用丝网印刷的方式,在后掠机翼前缘添加粗糙带,人工引入6mm与7mm驻波,通过升华法验证了计算与理论分析。     

基于线化稳定性理论的后掠翼边界层内横流稳定性研究

采用Mack方法求解忽略曲率和压缩性的三维线化稳定性方程（Orr-Sommerfeld方程）,针对无限展长后掠翼,计算分析边界层内不同波长横流驻波沿弦向的放大率,确定最不稳定波的波长。计算结果与参考文献中的实验结论和计算结果符合较好,表明该方法能够准确的预测出最不稳定波的波长,可以用来指导横流驻波主导下的后掠翼流动稳定性问题研究。     

三维机翼可压缩层湍流边界层计算

 从三维时间相依可压缩边界层动量积分方程和平均流动能积分方程出发,将二维可压缩层、湍流边界层积分方程算法推广到求解有限翼展后掠机翼的边界层流动。利用四步Runge-Kutta时间步进方案数值求解积分方程组,并利用当地时间步长加速迭代收敛。文中分析了数值方法的稳定性与收敛性,并考查了横向耗散项对计算结果的影响。算例表明,能获得令人满意的三维机翼定常可压缩层、湍流边界层的计算结果。     

后掠翼可压缩三维边界层稳定性计算

采用可压缩线化稳定性理论对后掠机翼的层流边界层稳定性进行了计算研究。对三维边界层中三维扰动波的可压缩稳定性方程组特征问题采用相关两步格式求解。与常用的初值法相比,具有不需预估特征值和快速收敛的特点。为得到渐近边界条件的准确表达式,导出了外边界上方程的解析解,这对精确的特征值计算和分析至关重要。为避免出现伪不稳定模态,选择了合适的稳定性方程组形式。最后,通过层流控制机翼和自然层流机翼的算例,研究了后掠翼上典型Tollmien-Schlichting(T-S)型扰动和横流(C-F)扰动稳定性问题,其结果与已有计算结果符合甚好。     

一种无人机变后掠翼机构的设计与仿真研究

提出了一种变后掠翼机构的设计方法,能根据机翼平面形状参数,合理布置变掠机构、确定机构部分主要参数,并能得到机翼的运动规律。首先提出了变后掠翼机构方案,将运动机构抽象为数学模型;然后对各参数进行详尽的设计和计算,并对机翼的运动进行了数值分析,在此基础上建立了一种变后掠翼机构设计方法。运用方法对具体的算例,建立机构三维模型,进行运动模拟仿真,CATIA运动仿真结果与MATLAB数值计算结果一致,证明了方法的有效性和合理性。     

后掠翼人工诱导横流驻波的升华法研究

改进了添加粗糙带方法,使用丝网印刷方式在后掠机翼前缘添加特定间距的粗糙带。成功地在边界层内引入了特定波长的驻波,提高了人工诱导驻波波长精度。在西北工业大学低湍流度风洞中使用升华法研究了单一波长驻波的发展趋势,其扰动发展趋势与计算结果吻合得很好。     

大展弦比机翼的静气动弹性分析

大型飞机普遍采用的大展弦比机翼,其气动弹性问题显得尤为突出。通过求解基于三维结构网格的Euler/边界层耦合方程组,采用结构模态分析法,对HIRENASD机翼的静气动弹性进行了数值模拟。计算得到了机翼静弹变形前和变形后的压力分布以及不同攻角下翼尖变形值。计算结果和实验值对比分析,验证了Euler/边界层耦合方程组求解器对静气动弹性数值模拟的准确性,表明所发展的数值模拟方法可以用于大展弦比机翼静气动特性的分析研究。     

飞机机翼表面霜冰的三维数值模拟

基于欧拉两相流理论对三维情况下飞机机翼表面的霜冰进行了数值模拟.根据水滴拟流体模型建立三维水滴控制方程;提出一套水滴控制方程的数值求解方法;由三维水滴流场的求解结果计算机翼表面的水滴收集特性,提出一种三维积冰外形的生成方法,完成了对飞机机翼表面霜冰的三维数值模拟.对ONERA M6机翼在不同迎角下霜冰的积冰情况进行了数值预测,并分析了结冰条件对积冰的影响.      

基于双eN方法的翼身组合体流动转捩自动判断

发展翼身组合体复杂外形流动转捩自动判断方法,对高亚声速民机自然层流（NLF）机翼设计具有重要意义。使用多块结构化网格和三维雷诺平均Navier-Stokes（RANS）方程求解器,耦合边界层方程求解和基于线性稳定性理论（LST）的完全双e^N方法,发展了一套可同时计及Tollmien-Schlichting波和横流不稳定性扰动诱导转捩的翼身组合体流动转捩自动判断方法。对DLR-F4翼身组合体绕流进行了转捩自动判断,将得到的转捩位置与试验结果进行比较,验证了所发展方法的正确性。使用上述方法对配置自然层流机翼的中短程民机翼身组合体外形进行了数值模拟,并将结果与单独机翼的转捩位置进行了对比,结果表明机身三维位移效应增强了自然层流后掠机翼边界层的横流不稳定性强度,导致翼根转捩位置提前至前缘区。     

机翼前缘后掠角对飞机RCS影响的数值模拟

为了在飞机总体设计时改善其隐身性能,对机翼前缘后掠角参数化可调的飞机三维数字样机的RCS特性进行了研究。使用CATIA软件,建立机翼前缘后掠角参数化可调的飞机三维数字样机;基于物理光学法和等效电磁流法,采用RCSAnsys软件,使用X波段雷达对飞机进行探测,雷达入射波的俯仰角在-15＆#176;、0＆#176;和15＆#176;条件下,数值模拟机翼前缘后掠角在-30＆#176;~＋60＆#176;之间变化时飞机的RCS特性,并对数值模拟结果进行数理统计分析。在机翼前缘后掠角变化的条件下,飞机RCS特性数值模拟结果表明：飞机头向RCS峰值之一的方位角与机翼前缘后掠角的角度相等;飞机头向RCS算术平均值特性为直机翼大、前掠翼和后掠翼小、大后掠翼更小;飞机侧向和尾向的RCS算术平均值变化相对不大。     

转捩模式与转捩准则预测高超声速边界层流动

对原始的k-ω-γ转捩模式和"层流+转捩准则"模型进行了改进,在2种方法中分别增加了横流模态时间尺度和横流转捩准则用于预测横流失稳诱导转捩。通过对网格预处理可并行计算获得边界层外缘信息以及边界层内横流速度。采用不同雷诺数条件下的0°攻角尖锥以及HIFiRE-5外形对2种方法预测高超声速边界层转捩的性能进行了比较分析。研究结果表明,2种方法均能正确反映高超声速边界层转捩起始位置和转捩区长度随雷诺数的变化趋势,但不能捕捉转捩区热流峰值;"层流+转捩准则"模型计算得到的传热系数在全湍流区较k-ω-γ转捩模式偏高。对于同时存在流向不稳定和横流不稳定的HIFiRE-5外形,改进的k-ω-γ转捩模式和改进的"层流+转捩准则"模型相比于原始的模型均能更加准确地预测中心线两侧横流失稳诱导形成的转捩;对于中心线附近因速度剖面拐点引起的边界层转捩,"层流+转捩准则"模型由于与边界层厚度相关,预测得到的转捩位置较试验结果靠前,k-ω-γ转捩模式与试验结果吻合很好。     

模型的转捩流动计算分析简

 为了在黏性流动数值模拟中实现边界层转捩的自动预测,将γ-Re_θt转捩模型引入到三维非结构混合网格的雷诺平均Navier-Stokes方程求解程序（HUNS3D）。该转捩模型由两个依赖当地变量定义的关于间歇因子和当地化转捩起始动量厚度雷诺数的输运方程组成,其数值求解算法与流场求解程序中湍流模型的求解方法相同。为了考察和验证HUNS3D程序中γ-Re_θt转捩模型对航空工程中的常见附面层自由转捩问题的预测精度,对低速平板流动、Aerospatial-A翼型、NLR 7301超临界翼型和NASA Trap wing 高升力构型等典型外形的自由转捩流动进行了计算,并将计算结果与相关试验结果进行了对比分析。算例结果表明：γ-Re_θt转捩模型对于转捩位置具有很好的敏感性,能比较准确地预测自然转捩和分离转捩,可以有效提高HUNS3D程序对实际流动的模拟能力和预测精度。     

C-γ-Reθ高超声速三维边界层转捩预测模型

针对国家数值风洞（NNW）工程高超声速三维边界层转捩预测需求,开展了高超声速边界层横流转捩判据及模型研究。采用线性稳定性分析e^N方法对高超声速转捩数据进行扩展,结合横流强度与表面粗糙度构造当地化的高超声速横流转捩判据,对基于Chant 2.0计算平台的高超声速修正γ-Re_θ转捩模型进行了横流模式拓展,建立了适用于高超声速三维边界层横流转捩预测的C-γ-Re_θ转捩模型。采用构建的转捩模型对多状态下的高超声速尖锥进行横流转捩预测,取得了与试验结果符合较好的预测效果。     

用主动柔顺壁运动控制边界层转捩

应用主动柔顺壁运动抑制平板边界层内扰动的增长和推迟边界层转捩的实验研究结果表明,边界转捩过程中的扰动即使已发展到非线性阶段,也能被主动柔顺壁的运动所控制。通过控制柔顺壁下方空腔内的输入声激励强度所产生的控制扰动,可以明显地抑制边界层内扰动的非线性增长率和推迟边界层转捩过程。改变控制扰动与初始扰动间的相位差,对于二维初始扰动,可以明显地提前和推迟边界层转捩过程,但对三维初始扰动,发迹相位差对边界层转     

高超声速边界层流动稳定性实验研究

高超声速条件下边界层转捩相关研究是近年来空气动力学领域的研究热点。通过采用基于纳米粒子示踪平面激光散射(Nano-tracer-based planar laser scattering,NPLS)技术以及温敏漆等测试技术,对高超声速边界层的流动稳定性开展了实验研究。通过NPLS技术,对圆锥边界层中第二模态波精细结构进行了测量,并基于时间相关的测量结果,对第二模态波的波长和频率进行了分析。针对裙锥边界层NPLS结果,计算了特定位置上功率谱空间分布结果,测量得到了高次谐波成分。通过温敏漆和NPLS结果,发现主导三角翼前缘边界层转捩的模态为行进横流模态,分析了该模态的特性,并且与Kulite传感器测量得到的频率进行比较。     

翼型几何模型对边界层转捩预测精度的影响

对飞行器表面的边界层转捩现象进行准确地预测,是设计高性能层流翼型的关键。基于边界层数值模拟和现代CFD技术的特点,选取与SSTk-ω湍流模型相耦合的γ-Reθt转捩模型,对不同类型的翼型边界层转捩进行预测。结果表明:在转捩计算过程中出现的摩阻曲线震荡,是由翼型表面计算网格不光顺所引起的;提高翼型几何模型光顺程度能够在不同程度上提高翼型气动特性的计算精度,并能得到更加准确的转捩预测结果。     

离散粗糙元诱发边界层转捩的实验研究

针对直升机转子叶片模型表面离散粗糙元诱发边界层转捩问题开展了实验研究,分析不同雷诺数下粗糙元尺寸参数对转捩位置的影响。实验在中航工业气动院直升机转子叶片模拟装置进行,在模型转子叶片表面布置不同尺寸的离散柱状粗糙元,利用红外热像技术探测边界层转捩,并提出一种基于湍流／层流区域面积比的转捩位置判定准则,目的是实现边界层转捩位置自动识别,进而分析粗糙元尺寸参数对转捩位置的影响。实验转速为300至600 r／m,对应叶尖切向速度为25～40 m／s。实现了对旋转叶片边界层转捩位置的定量测量,通过实验验证,转捩位置判定算法正确可靠,初步得到了不同高度 DRE 诱发转捩位置与雷诺数之间的关系,随着粗糙元高度的增加,转捩位置逐渐靠前。     

高超声速锥柱裙模型边界层转捩的弹道靶实验

为研究高超声速边界层转捩现象、给边界层计算提供可靠的对比数据,在中国空气动力研究与发展中心趟高速弹道靶上开展了锥柱裙模型高超声速边界层转捩的自由飞实验。所采用的锥柱裙模型全长105mm,飞行速度1．94km／s（Ma=5．65）,单位雷诺数4．32×10^7～1．20×10^8m-1。使用激光阴影成像技术,获得了锥柱裙模型边界层转捩和湍流边界层发展的图像,测得的湍流边界层厚度在0．6～2．2mm之间,湍流涡的流向尺寸与边界层厚度的比值介于0．3～0．8之间且沿流向呈下降趋势。实验结果表明：弹道靶实验能够获得给定飞行环境下的高超声速边界层转捩图像,从图像中可以清晰判断转捩位置或区域、测量边界层厚度和分析湍流涡的尺寸。     

前缘锯齿对边界层不稳定噪声的影响

为探索仿生学前缘锯齿结构的降噪规律,试验研究了低雷诺数到中等雷诺数（Re=（2~8）×10⁵）不同攻角状态下9种前缘锯齿结构对叶片层流边界层不稳定噪声的影响。研究表明：前缘锯齿可以减弱甚至完全抑制边界层不稳定噪声,降噪效果对锯齿振幅和锯齿波长均比较敏感,锯齿振幅越大、波长越小,降噪效果越好,降噪量可达30 dB;前缘锯齿结构可以诱导产生流向涡,影响叶片下游边界层流动,破坏声学反馈回路;前缘锯齿对边界层不稳定噪声峰值频率没有影响。     

激波风洞边界层转捩测量技术及应用

高超声速边界层转捩对摩阻、传热等有重要影响。在高超声速飞行器研制中,迫切希望能精确预测和控制边界层转捩。激波风洞作为高超声速气动热环境试验的主要地面模拟设备,是研究高超声速边界层转捩的重要设备。但激波风洞原有测量技术适用于工程型号试验,需要依据高超声速边界层转捩特点进行适应性改造和升级。依据高超声速边界层转捩过程中的热流、压力、密度等物理参数变化,发展了薄膜热流传感器测热技术、温敏热图测量技术、高频脉动压力测量技术、高清晰度纹影显示技术等适用于激波风洞的边界层转捩测量技术。并针对头部钝度0.05 mm的半锥角7°尖锥模型,在中国空气动力研究与发展中心Ø2 m激波风洞（FD-14A）马赫数10、单位雷诺数1.2×10⁷/m的流场条件下开展了边界层转捩试验。采用多种转捩测量技术同时进行测量,获得尖锥模型表面边界层转捩情况、边界层脉动压力频谱特征、边界层内清晰的第2模态波和湍流斑纹影图像,不同测量技术获取的试验结果可相互印证,线性稳定性理论分析结果与试验结果相吻合。