大攻角流动的非定常特性分析

本文首先通过动态测量手段研究了大攻角流场中的能量成分,分析了各种成分引发的非定常性对大攻角流动的影响。结果表明大攻角流动本身是一种“涡行为”,非定常性对大攻角流动非对称性的影响较弱,影响大攻角流动的本质因素是涡的动力不稳定性对于前体不确定小扰动的响应。     

合成射流大攻角非对称涡控制的试验研究

利用合成射流对细长旋成体大攻角非对称涡控制进行了研究,基于合成射流激励器设计了一频率高达1kHz的非定常小扰动控制机构,并将其成功应用于大攻角非对称涡主动流动控制。应用天平测力和七孔探针流场测试技术,研究了合成射流非定常小扰动电压和频率对非对称涡的控制特性和规律。结果表明,采用合成射流能够完全消除背涡的非对称性,扰动频率是影响非对称涡控制的一个关键参数,高频扰动下模型背风区非对称涡结构趋于无控制流态。且文中结果发现,当攻角α=57.5°、非定常小扰动频率fs=150Hz时,即可将非对称涡完全控制成为对称涡。     

大攻角非对称流动的非定常弱扰动控制

 研制了一种新的大攻角细长旋成体非对称涡的主动控制技术, 即在细长旋成体头部施加非定常弱扰动来控制头部非对称背涡。应用七孔探针测量的空间截面流场揭示了非定常控制下非对称涡变成对称涡的流态特征。测力试验研究结果表明该方法不仅能完全消除背涡的非对称性及其产生的侧向力, 并且有效控制攻角范围从30b直到80b。     

大攻角地效翼非定常流动实验研究

大攻角地效翼非定常流动研究对大攻角地效翼非定常空气动力特性和流动控制具有重要的意义。对攻角为18°的NACA0012地效翼在模型风洞中进行实验,利用三维热线风速仪对地效翼后缘附近尾迹区内测点的速度脉动进行采样,分析不同飞行高度下大攻角地效翼的非定常流动特性。结果表明：地面效应下,大攻角地效翼后缘涡系形成位置后移;随着地效翼高度的降低,非定常流动的脉动频域宽度减小,速度脉动特征频率减小。     

有翼飞行器大攻角分离涡模拟

基于势流方程,采用涡格法和三维离散涡相结合,模拟具有升力面分离涡的大攻角弹冀、弹体、尾翼组合体绕流。该方法适用于高亚音速、小攻角或大攻角有翼飞行器的压力分布、法向力和压力中心等气动特性的数值计算,并具有显示三维涡流场的能力。     

采用DMD方法研究叶栅不同攻角的拟序结构

为分析平面叶栅分离流非定常拟序流动特征,对三个不同攻角下的叶栅进行了单通道的大涡模拟仿真,并采用动力模态分解（DMD）三个工况的流场结构进行了分析。DMD方法对包含复杂时空信息的叶栅分离流流场进行了解耦,剥离出了反映流场主要动力信息的模态,获得了其频率和与之对应的空间结构。并且通过DMD方法,将原本需要研究大量不同时刻的流场,转移到仅需要对少量模态的研究即可,实现了保留主要动力特征的低维近似。通过DMD分析表明:气流经过叶片前缘产生流动分离,形成不稳定的剪切涡结构,它和尾迹区脱落涡相互耦合,并形成新的拟序结构。随着攻角的增大,前缘剪切涡及其与尾迹涡的耦合也同时增强,流场变得更加复杂。      

用非结构动网格技术计算确定再入飞行器配平攻角

提出了用数值计算直接确定飞行器配平攻角的自由振动法.该方法采用非结构动网格技术,使飞行器在气动力作用下绕通过质心的轴线自由振动,最终衰减到配平攻角.构造了一种MUSCL类型有限体积格式,求解考虑了动网格效应的三维非定常Euler方程,采用弹簧近似技术实现非结构网格的变形运动.作为例证,计算确定了载人飞船返回舱的配平攻角,计算结果与实验结果接近.该方法外形适应能力强,效率较高,具有较好的应用前景.     

混合流场控制的大涡模拟

采用大涡模拟方法研究了三维小扰动下,混合流场大尺度拟序结构的产生和演化过程,捕捉了展向涡的卷起、配对、合并,以及二次流向涡的出现等大尺度的三维拟序结构,分析了拟序结构与入口扰动方式之间的内在联系,再现了涡卷自身撕裂而引发转捩的现象。     

扩压叶栅叶尖非定常流动的涡动力学机理

为了揭示压气机叶尖区旋涡结构与流动非定常性之间的关联，采用URANS对一亚声速平面扩压叶栅在不同攻角下的流场进行了求解，并借助Q判据提取了叶尖瞬态涡系结构。结果表明：泄漏涡的破碎现象能够通过诱导新的涡结构间接作用于相邻通道的叶尖流动，是导致叶尖流场失稳的关键因素。在-0.3°和+0.7°攻角下，叶尖泄漏涡发生了螺旋破碎，并伴随有非定常诱导涡的出现，诱导涡对相邻叶片载荷的影响使得叶尖泄漏涡发生周期性摆动；在+1.7°攻角下，泄漏涡破碎会导致反流涡的形成，反流涡的输运会给叶片载荷和来流攻角带来非定常扰动，反过来又会作用于泄漏涡的破碎和反流涡的生成，最终表现为一种自维持的非定常流动现象。     

压气机二维叶栅涡脱落的数值模拟

对某压气机二维叶栅的非定常分离流场进行了数值模拟,通过对多种工况的计算,进行了频谱分析,对叶栅非定常流动的流场结构和流动机理做了初步的探讨。在来流均匀,定常边界条件下,叶栅内流动仍然表现出强烈的非定常性;在大攻角下有类卡门涡的脱落;旋涡脱落的频率,随着攻角和马赫数的变化而变化;分离点的位置随着攻角和马赫数的变化而变化。     

高超声速钝楔边界层转捩大涡模拟

 以五阶迎风和八阶对称格式混合差分格式求解三维可压缩滤波Navier-Stokes方程，对来流Mach数为6.0、半锥角5°的高超声速空间发展钝楔边界层转捩至完全湍流进行了大涡模拟。时间推进采用紧致存储三阶Runge-Kutta方法，亚格子尺度模型为Driest因子修正的Smagorinsky涡黏性模型。通过定常流场入口边界附近吹／吸引入不稳定扰动斜波的方法数值模拟得到了层流失稳转捩直至完全湍流的空间发展全过程。对扰动的线性、非线性增长以及湍流斑的形成和发展进行了分析，给出了转捩及完全湍流下的速度相关量统计并与实验、DNS结果进行了对比分析，计算结果与理论及实验吻合。     

附面层吸入式进气道内流动损失特性

用数值仿真手段，对进气道内流总压损失采用逐段对比分析的方法，研究了附面层吸入式（BLI）进气道相对于常规S弯进气道的内流总压损失特性。结果表明：进气道出口马赫数不变，较低来流马赫数和较高来流马赫数工况下（方案对应马赫数0.3以下和马赫数0.6以上），BLI进气道产生的流动损失比常规S弯进气道的大，差量达2.4%；中等来流马赫数工况下（本文方案对应马赫数0.3与马赫数0.6之间），BLI进气道产生的流动损失比常规S弯进气道的略小，差量在0.3%以内；流动损失特性间的差异是由于BLI进气道进口前壁面与进口低能附面流改变了进口段流动特性、及在S弯管道内发展的综合作用结果。     

超声速空气介质自由扩散混合流涡的演化

在不同谐波扰动促发方式下,观察二维超声速空气/空气自由扩散混合流中涡的形成及演化.考虑含空气组分扩散效应的可压缩Navier-Stokes方程,对流项采用3阶迎风紧致格式离散,输运项采用6阶对称紧致格式离散,非定常时间推进采用3阶紧致存储显式Runge-Kutta方法.在无相差情形下,获得了大尺度基频涡结构的饱和、一次对并、二次对并、三涡对并等现象.在Ma_c=0.3低对流Mach数下,基频涡较饱满,但流向尺度较小.受空气真实气体特性影响,在Ma_c=0.8高对流Mach数下未发现涡/小激波结构.      

高马赫数下激波湍流边界层干扰数值研究

应用GAO—YONG可压缩湍流方程组数值模拟了入射斜激波／平板湍流边界层相互干扰现象,计算了来流马赫数为5．0,激波入射角度分别为15．876°、23．287°两种不同激波干扰强度下的流场。计算程序中的对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge—Kutta显式时间推进法求解空间离散后的控制方程。计算较好地模拟了高马赫数下的激波／湍流边界层干扰的流场结构,位移边界层厚度,动量损失厚度等,也比较准确地预测了平板壁面压力、摩阻系数等气动力参数的分布。     

自由流单位Re数对边界层转捩的影响

利用激波管实验原理和实验技术,在截面为94 ×94mm~2的激波管实验段中,测量了自由流单位Re数单一变化时壁面边界层转捩Re数的变化规律。为此研制了一种可直接显示边界层从层流向湍流过渡全过程的热电模拟网络。根据大量测量结果,并以自由流单位Re数为自变量,整理出预测平板边界层转捩Re数的经验公式。     

气动声学直接数值模拟的空间差分格式分析

流场内声波的直接数值模拟对数值格式的耗散、色散特性提出了严格要求,首先介绍了两种适于声学模拟空间差分格式——紧致差分格式与色散关系保持(DRP)格式的基本思想,然后分析了格式的色散及其与四阶龙格-库塔法结合时的耗散特性。根据Fourier误差分析和数值模拟发现:无论从格式可有效求解的波数范围还是从与龙格-库塔时间推进结合后的耗散特性来看六阶紧致格式的可有效求解范围要大于DRP格式;不过,数值计算表明紧致格式的计算量远远超过DRP格式,对计算机资源提出了较高的要求。      

来流情况下入口状态对吸气式激光推力器冲量耦合系数的影响

激光推力器在高空环境下飞行不可避免地会受到来流的影响,而其入口状态直接影响激光聚焦点火之前的定常流场状态,从而影响推力器的推进性能。在来流马赫数为5的条件下,对某种吸气式推力器模型的入口在打开、关闭以及半开3种状态下的内外流场演化过程进行了数值模拟。结果表明：入口打开时,轴向阻力低于同工况下入口关闭和半开的情况;入口关闭时,聚焦区的激波传播速度明显增加,正推力峰值明显增大,抵消了入口关闭带来的额外负推力;入口半开时,正推力较入口关闭没有得到有效提高,轴向阻力较入口打开却大幅增加,导致其冲量耦合系数最低,推进性能最差。通过计算和分析预测,入口构形的改进设计可大幅度提高来流情况下推力器的推进性能,对入口构形优化还需进行大量工作。     

Transonic wing stall of a blended flying wing common research model based on DDES method

Numerical simulation of wing stall of a blended flying wing configuration at transonic speed was conducted using both delayed detached eddy simulation(DDES) and unsteady Reynolds-averaged Navier-Stokes(URANS) equations methods based on the shear stress transport(SST) turbulence model for a free-stream Mach number 0.9 and a Reynolds number 9.6 × 10~6. A joint time step/grid density study is performed based on power spectrum density(PSD) analysis of the frequency content of forces or moments, and medium mesh and the normalized time scale0.010 were suggested for this simulation. The simulation results show that the DDES methods perform more precisely than the URANS method and the aerodynamic coefficient results from DDES method compare very well with the experiment data. The angle of attack of nonlinear vortex lift and abrupt wing stall of DDES results compare well with the experimental data. The flow structure of the DDES computation shows that the wing stall is caused mainly by the leeward vortex breakdown which occurred at x/x_(cr)= 0.6 at angle of attack of 14°. The DDES methods show advantage in the simulation problem with separation flow. The computed result shows that a shock/vortex interaction is responsible for the wing stall caused by the vortex breakdown. The balance of the vortex strength and axial flow, and the shock strength, is examined to provide an explanation of the sensitivity of the breakdown location. Wing body thickness has a great influence on shock and shock/vortex interactions, which can make a significant difference to the vortex breakdown behavior and stall characteristic of the blended flying wing configuration.     

大攻角风洞自由飞动导数实验研究

实验是在0.2×0.2m~2超音速风洞中进行的,其来流Mach数为2.0。模型为带翼导弹模型。实验时用气体发射装置将模型向实验段上游发射;并以另一气流横吹模型使其呈大攻角状态飞行。最大攻角可达100°。采用高速鼓轮相机记录模型飞行姿态随时间的变化。用改进的数值积分法进行数据处理,求得俯仰阻尼导数随攻角的变化。     

Improvement of forebody/inlet integration for hypersonic vehicle

Results of a numerical (CFD) study of the influence of the forebody shape on local flow parameters at a bottom-mounted inlet entrance are presented. The free-stream Mach number is assumed to be 3.5–7.0. Some recommendations on forebody shape optimization are provided. Main characteristics of the air inlet are evaluated.