翼型前缘变形对动态失速效应影响的数值计算

翼型或机翼的动态失速效应所引起的低头力矩和正气动阻尼限制了飞行器气动性能的提高,甚至可能诱导发生不稳定运动。应用于小尺寸机翼的前缘动态变形（DDLE）技术,通过实时改变前缘形状,能够改善翼型前缘区域的速度梯度,进而抑制动态失速效应。采用转捩剪切应力输运（SST）黏性模型结合分区混合动态网格技术,研究了这种前缘变形对机翼俯仰运动所引起的非定常流动的影响,得到通过小幅度前缘变形抑制和延迟动态失速的方法,从而提高翼型的气动性能。翼型NAC A0012的数值模拟结果与动态失速风洞试验结果比较表明：所使用的数值计算方法能够较为准确地模拟翼型在动态失速过程中升力系数与俯仰力矩系数的变化情况,可用于研究前缘变形对翼型俯仰运动所引起的非定常流动的影响。前缘动态变形翼型俯仰运动过程的非定常流场的数值模拟表明：在大迎角下不同幅度的前缘下垂运动能够抑制流动分离的发生,从而抑制动态失速,但在大迎角下小幅度高频率的前缘下垂变形能更高效地抑制动态失速;前缘变形幅度以及变形沿中弧线的分布对升力系数和俯仰力矩系数的影响并不明显。     

可折叠翼变形飞行器气动特性研究

本文设计了一种可折叠翼概念无人机的变形飞行器模型,提出了一种变形飞行器气动系数的定义,通过数值计算,对比了飞行器变形前后以及不同内翼折角状态下的气动特性,分析了飞行器气动系数与内翼折角、迎角以及飞行马赫数的变化关系。     

涡桨飞机滑流影响的非定常数值模拟验证

通过求解绝对坐标系下的非定常雷诺平均Navier-Stokes方程,计算了螺旋桨非定常滑流对某涡桨飞机起飞构型低速状态的气动特性影响。为描述螺旋桨桨叶的相对运动,采用了动态重叠网格技术,并在并行环境下采用全隐式双时间步方法和多重网格技术来保证时间精度和提高计算效率,计算结果与试验值相当吻合,研究表明:滑流区内襟翼上的动压得到明显增强,同时高速滑流也起到一个吹除在大偏度襟翼上堆积的附面层作用,延迟了襟翼上气流的分离,这两个因素明显提高了襟翼效率,另外经过机翼的导流作用后,滑流对襟翼的洗流不像机翼那么明显。本文还将有无滑流的流场进行对比,通过当地动压增量来定义滑流的加速效应边界,以及通过当地气流角增量来定义滑流的洗流效应边界。该方法能较好地捕捉和解释滑流对飞机部件干扰而使得飞机方向安定性呈现的非线性现象,初步揭示了螺旋桨滑流复杂尾迹流动的特点。     

基于动态混合网格的双鱼串列巡游数值模拟

观察研究发现,在鱼类的群游中,后鱼可以利用前鱼尾迹中有利的涡结构来增大自己的推力,减少功耗.本文利用作者以往发展的动态混合网格生成技术和相应的不可压缩流非定常计算方法,在单鱼巡游数值模拟的基础上,数值模拟了双鱼串列巡游的非定常流动,分析了摆动频率、摆动相位差以及流向间距对推力、功耗等流体动力学特性的影响.     

TSTO级间分离气动特性数值仿真分析

两级入轨（TSTO）重复使用飞行器并联级间分离是一二级飞行器设计必须考虑的关键因素。采用国家数值风洞工程通用CFD软件NNW-FlowStar,基于各向异性非结构混合网格及重叠网格技术,研究了TSTO在侧滑角下的并联级间分离过程。软件通过自主改进的Roe熵修正方法、节点高斯梯度算法实现高超声速流场的精细模拟,运用“物面相交”准则实现多套网格间的并行挖洞。采用双时间步求解非定常分离过程,内迭代计算采用LU-SGS方法。文中首先采用WPFS投放标模验证了软件在常规分离投放方面的数值模拟精度;然后模拟获得TSTO组合体定常气动特性,分析了不同攻角下一二级的气动特性以及可能的安全分离状态;最后数值模拟了TSTO在负攻角、正侧滑角下的并联级间分离过程,与对应的CTS试验结果比较,两级质心位移的计算值和实验值随时间变化具有很好的一致性,偏航角及俯仰角在变化规律和量值上比较吻合,表明软件具备升力体并联级间分离类问题的数值模拟能力。     

基于动态混合网格的多体分离数值模拟方法

在以往发展的动态混合网格技术的基础上,耦合六自由度运动方程计算和非定常流场计算,建立了一种多体分离问题的非定常数值模拟方法。首先采用刚体运动的六自由度运动方程,建立刚体运动轨迹计算模块,然后以该模块为纽带将以往建立的非定常计算方法与改进的基于Delaunay背景网格插值方法和局部网格重构方法相结合的动态混合网格生成方法耦合起来,形成了精度较高、效率较快的多体分离数值模拟方法。利用该方法对典型多体分离问题进行了数值模拟,计算结果与风洞实验结果吻合较好。     

基于非结构网格的气动弹性数值方法研究

基于非结构网格和网格变形技术,研究了耦合求解流体和弹性结构体互相作用问题的数值方法.气动弹性模型采用了耦合求解可压缩雷诺平均的Navier-Stokes方程和线性结构动力学方程的方法.为了能处理在复杂几何域内的问题,采用了基于混合单元的非结构网格有限体积方法.开发的方法在二维翼型简谐振动进行了验证并应用于AGARD 445.6翼型颤振边界的数值模拟.      

两种典型非结构网格变形方法特性对比研究

非定常流动仿真中常常遇到运动边界的情形,基于物理模型的线性弹簧法和基于数学方法的径向基函数(RBFs)插值法是实现计算网格随运动边界变化的两种主要的非结构网格变形方法。以NACA0012翼型的俯仰运动为例,对上述两种方法的变形特性进行较详细地定量分析研究,比较线性弹簧法和RBFs方法的CPU计算时间和最大变形能力,分析两种方法的网格质量随俯仰角度和运动步数的关系,并给出典型状态下网格单元质量云图和流场计算结果。结果表明:基于RBFs插值的网格变形方法计算效率较高,变形能力较强且能保持较高的网格质量,是一种高效的网格变形方法。     

考虑喷流效应的机载导弹发射及气动干扰数值模拟

研究比较了机载导弹发射时是否考虑导弹尾喷流情况下导弹运动特性及机-弹气动干扰影响.通过动态结构嵌套网格技术实现了导弹运动的模拟,结合刚体六自由度运动方程求解三维非定常Euler方程得到流场信息.导弹推力由发动机燃烧室喷流计算得到;当不考虑喷流效应时,推力通过直接在导弹尾部给定非定常作用力实现.应用此方法模拟了基于类"全...     

基于线性弹性体的网格变形方法研究

网格变形方法在计算流体力学的各个方面都发挥着极其重要的作用,是各种应用问题解决方案的基本组成部分.研究和实现了基于线性弹性体的网格变形方法,对基于弹性体法的网格变形方法进行了综述,从弹性方程理论、高效求解方法等方面进行了介绍,给出了方法的最优参数选择,并从基本变形能力、多段翼型设计、非定常动态特性模拟等方面,展示基于线性弹性体网格变形方法的研究和应用进展.结果表明:实现的线弹性体网格变形方法特别善于处理几何大变形,尤其是网格扭转比较大的应用,并且变形后网格品质高、鲁棒性强.     

柔性变后掠飞行器非定常气动特性数值研究

为研究柔性变后掠飞行器变形过程中的非定常气动力,对柔性变后掠飞行器进行了非定常数值仿真。首先分析了柔性变后掠飞行器在特定后掠角下的定常气动特性,接着选用三种变后掠周期进行了非定常计算,分析了不同变后掠速度对飞行器气动特性的影响,以及定常与非定常气动特性的差别,并研究了这种差异产生的原因。结果表明:柔性变后掠飞行器通过后掠角的改变可以使实时气动性能达到最优;不同变后掠速度引起的气动力差异不大;定常气动力与非定常气动力最大差异不超过7%,其差异主要是由于机翼上气动力的差异引起;非定常计算的升力、阻力系数大于定常结果,俯仰力矩系数与定常计算值差异不大。非定常气动力的产生机理是由于机翼的附加速度所引起的,与流场迟滞无关。总体上看,攻角小于14°时,小后掠可以取得较大的升力、阻力系数;大于14°攻角,大后掠的升力、阻力系数较大;所有后掠角均在4°攻角处取得最大升阻比且小后掠角的升阻比较大;当升力系数小于1.28时,小后掠角产生较小的阻力系数,超过这一数值,大后掠角的阻力系数较小。     

折叠翼变体飞行器非定常气动特性实验研究

折叠翼变体飞行器是一种可以在飞行中改变自身气动外形的新型飞行器。研制出了一种折叠翼变体飞行器的风洞实验模型,在风洞实验中测得了模型不同变体位置下的气动力以及进行变体运动时气动力的动态变化过程,并通过PIV实验手段获得模型周围的流场在变体运动过程中的变化情况。结果表明：在机翼变形过程中,折叠翼模型有明显的非定常气动现象产生,而且折叠变形的速度越大,非定常现象越明显。出现非定常现象的主要原因是变体运动对机翼前缘涡的影响。     

折叠机翼变体飞机纵向操纵性与稳定性研究

针对一种飞翼布局折叠机翼变体飞机方案,建立了相应的研究模型,使用涡格估算方法计算出机翼折叠角度对全机纵向静稳定性的影响.利用飞机纵向小扰动运动方程,分别得到了机翼展开和折叠状态的长、短周期模态,并对其动稳定性进行了讨论.结合工程估算方法和风洞试验方法,计算出内段机翼折叠过程全机力矩系数、升降舵操纵导数以及升降舵配平偏角...     

滑动蒙皮变后掠无人机非定常气动特性研究

为了探索飞机变形过程中的非定常气动特性及其机理,基于滑动蒙皮变后掠无人机这一可变形飞行器大尺度全局变形研究平台,通过风洞实验对其进行了非定常气动特性研究。结果表明:无人机变后掠过程中非定常气动特性曲线在相应的准定常曲线周围形成滞回环;无人机变后掠速率越快,滞回效应越显著;无人机以约7.5°/s的速率变后掠,会使非定常气动特性数值偏离相应的准定常数值5%以上;变后掠无人机非定常气动特性产生的主要原因可能在于流场结构迟滞和机翼附加速度。     

旋转机翼飞机旋翼模式前飞状态干扰气动特性

和传统直升机相比,旋转机翼（CRW）飞机在旋翼模式前飞时各部件之间存在更为严重的气动干扰。为了获得旋转机翼/机身/鸭翼/平尾之间的非定常气动干扰规律,基于运动嵌套网格技术,通过求解三维非定常雷诺平均Navier-Stokes（URANS）方程,建立了旋翼前飞流场数值模拟方法。首先对传统直升机旋翼/机身干扰模型进行了计算,验证了方法的可靠性,然后对某旋转机翼飞机全机在旋翼模式前飞状态下的非定常流场进行了数值模拟,并对各个气动部件上的非定常气动力和力矩的变化进行了分析。结果表明：飞机在旋翼模式前飞时,机身部件对旋转机翼的干扰较弱,在经过机身上方时拉力峰值仅略有增加;旋转机翼对鸭翼和垂尾干扰较弱,对机身和平尾干扰较强,随着前飞速度增大,旋转机翼对平尾的干扰会产生较大的升力损失和抬头力矩,需要引起重视。计算结果为该类飞行器的总体综合设计提供了参考。     

多段翼型局部主动变形流动控制的非定常数值模拟

对30P30N三段翼型失速攻角附近的分离流动进行了数值模拟研究。为了抑制大攻角时背风区的流动分离,在主翼段上表面引入了行波壁变形模型和抛物型局部主动振动模型,利用作者以往发展的动态混合网格技术和相应的非定常计算方法,对变形过程中的非定常分离流动进行了数值模拟,分析了各种变形参数对流动分离的影响。计算结果表明,在适当的条件下,局部主动变形能够抑制翼型背风区的分离,由此可以起到增升减阻的作用,改善翼型的气动性能。     

非定常预处理方法在倾转旋翼飞行器悬停状态气动干扰模拟中的应用

针对倾转旋翼飞行器的多种速度尺度并存的复杂流动问题,发展了能够统一求解同时包含不可压缩流动和可压缩流动的非定常预处理方法和一种考虑了流场物理信息传播方向的预处理远场边界条件。结合适用于倾转旋翼飞行器的复杂外形及大位移问题的非结构混合网格技术和重叠插值技术,建立了倾转旋翼飞行器的数值模拟方法。在此基础上,首先通过悬停算例验证了数值方法的有效性,确定了选择预处理截断参数的原则。然后探究了倾转旋翼飞行器直升机模式悬停状态的尾迹涡结构演化过程,研究了“喷泉效应”和“地面效应”对倾转旋翼气动性能的影响。计算结果表明,旋翼尾迹涡的耗散在不同方位角存在较大差异。“地面效应”抑制了旋翼尾迹的发展,提高了旋翼整个方位角区域的升力系数,而不仅限于机翼上方方位角区域。 “地面效应”作用强于“喷泉效应”作用,使旋翼相对于孤立旋翼拉力系数略微增大了5.26%,但同时整个倾转旋翼飞行器的总的向上的力只有孤立旋翼拉力的82.46%。结论表明预处理方法和预处理远场边界条件能够较好地模拟倾转旋翼飞行器悬停状态的干扰流场,兼顾了数值计算的耗散性和稳定性。气动干扰对倾转旋翼飞行器整体气动性能有重大影响。      

对转螺旋桨流场气动干扰数值模拟

基于结构网格动态面搭接技术,通过求解三维非定常雷诺平均纳维-斯托克斯(RANS)方程,对单独螺旋桨(SRP)滑流流场进行了数值模拟,结果与风洞试验数据吻合良好,在此基础上深入研究了对转螺旋桨（CRP）的非定常滑流效应及前后桨之间的气动干扰现象。计算结果表明:所采用的非定常RANS方法能够良好刻画对转螺旋桨滑流流场的发展变化特征,适用于分析前后桨气动干扰特性。前后桨桨尖涡的复杂相互作用是导致对转螺旋桨气动干扰的直接原因,后桨桨尖涡对于前桨桨尖涡具有一定耗散作用。在气动干扰影响下,对转螺旋桨气动力会产生周期性波动,波动周期与螺旋桨桨叶数相关, 6×6对转螺旋桨在1个旋转周期内产生了12次波动。气流穿过对转螺旋桨会发生两次加速,因此在相同工况下,对转螺旋桨的拉力系数和功率系数相比于单独螺旋桨分别增大约1倍,效率提升约1.5%。      

涵道螺旋桨与孤立螺旋桨气动特性的 数值模拟对比

基于非结构动态嵌套网格方法,对涵道螺旋桨与孤立螺旋桨的气动特性进行了非定常Euler方程数值模拟对比研究和分析.将计算区域分成旋转区与静止区两个子域,生成相互重叠的嵌套网格,有效地解决了螺旋桨与涵道以及静止部分桨毂之间的相对旋转问题.计算结果与实验结果相一致,验证了计算方法的正确性.涵道入口前缘形成的负压区是产生涵道附加拉力的原因所在.对比计算结果表明:与孤立螺旋桨相比,涵道的存在改善了螺旋桨桨尖区域的绕流特性,减小了桨尖损失;相同转速情况下,涵道螺旋桨产生更大的拉力,而所需功率略小,涵道螺旋桨系统具有更高的气动效率.      

Numerical study of aerodynamic characteristics of FSW aircraft with different wing positions under supersonic condition

This paper investigates the influence of forward-swept wing(FSW) positions on the aerodynamic characteristics of aircraft under supersonic condition(Ma = 1.5). The numerical method based on Reynolds-averaged Navier–Stokes(RANS) equations, Spalart–Allmaras(S–A) turbulence model and implicit algorithm is utilized to simulate the flow field of the aircraft. The aerodynamic parameters and flow field structures of the horizontal tail and the whole aircraft are presented. The results demonstrate that the spanwise flow of FSW flows from the wingtip to the wing root, generating an upper wing surface vortex and a trailing edge vortex nearby the wing root.The vortexes generated by FSW have a strong downwash effect on the tail. The lower the vertical position of FSW, the stronger the downwash effect on tail. Therefore, the effective angle of attack of tail becomes smaller. In addition, the lift coefficient, drag coefficient and lift–drag ratio of tail decrease, and the center of pressure of tail moves backward gradually. For the whole aircraft,the lower the vertical position of FSW, the smaller lift, drag and center of pressure coefficients of aircraft. The closer the FSW moves towards tail, the bigger pitching moment and center of pressure coefficients of the whole aircraft, but the lift and drag characteristics of the horizontal tail and the whole aircraft are basically unchanged. The results have potential application for the design of new concept aircraft.