Unsteady Flow Variability Driven by Rotor-stator Interaction at Rotor Exit

Numerical investigation of the unsteady flow variability driven by rotorstator interaction in a transonic axial compressor is performed. Two models with close and far axial gap between rotor and stator rows are studied in the simulation. Particular attention is attached to the analysis of mechanisms involved in driving rotor wake oscillation, rotor wake skewing and flow angle fluctuation at rotor exit. The results show that smaller axial gap is favorable to enhance the interaction in the region between two adjacent rows, and the fluctuation of the static pressure difference between two sides of rotor wake is improved by potential field from down stator, which is the driving force for rotor wake oscillation. The interaction between rotor and stator is weakened by increasing axial distance, rotor wake shifts to suction side of rotor blade with 5%-10% of rotor pitch, the absolute value of flow angle at rotor exit is less than that in the case of close interspace for every time step, and the fluctuation amplitude is also decreased.     

多块协调变形的网格变形技术及其应用

基于多块网格协调变形的理论,提出了适用于叶轮机内流多物边界的更具一般性的网格变形技术,先获得块边界处的变形量,然后插值得到块内网格各点的变形量。通过二维振荡叶栅和带叶尖间隙的三维跨音风扇叶片的非定常绕流算例,证明网格变形技术在叶轮机气弹稳定性数值模拟中的应用是有效的,能够在保证网格高精度的条件下,提高网格变形速度。     

基于滑移网格的临近空间螺旋桨流场数值仿真

基于滑移网格模型,考虑RNG k-ε湍流模型,通过求解三维非定常N-S方程,仿真研究了螺旋桨非定常旋转流场,比较研究了20km临近空间环境下不同前进速度、不同转速下螺旋桨气动性能的异同.结果显示,拉力随转速的增大而增大,效率随转速的增大先增大后减小;桨叶周围的速度及速度场的分布范围,涡流强度的大小及涡流强度的分布范围都随转速的增大而增大;桨叶旋转运动产生的诱导速度对螺旋桨前面流场的影响范围较小,对螺旋桨后面流场的影响范围较大.     

Computations of unsteady transonic flow over airfoil at low Reynolds number

Transonic flow over a thin airfoil at low Reynolds number was studied numerically by directly solving two-dimensional full Navier-Stokes equations through 5th order weighted essentially non-oscillatory(WENO) scheme without using any turbulence model.A series of distinguished unsteady phenomena for a thin 2-D transonic airfoil flow were presented.Due to continuous adverse pressure gradient in the subsonic flow downstream of the sonic line, the unsteady separated boundary layer with main vortex and secondary vortex was developed at the rear of the airfoil.At the trailing edge,the vortex-shedding was characterized by periodical connection of the main vortex and secondary vortex on the other side of the airfoil.The unsteady separation and vortex-shedding occurred with the same period.On the airfoil surface,the average pulse pressure related to the unsteady supersonic region was obviously smaller than that related to the vortex-shedding at the trailing edge.With the attack angle increasing from 0° to 2°, the frequency of vortex-shedding decreases about 4.2%.At last, the turbulence intensity and many second-order statistics in the wake region were investigated.      

高超声速全动翼面全时域耦合分析方法及应用

在流场-结构温度场同步计算方法的基础上,建立了多物理场全时域耦合分析方法,将方法应用于沿轨道运动的高超声速全动翼面热气动弹性稳定性分析。采用基于有限体积模型的CFD同步计算方法求解高超声速流场和结构温度场,建立映射关系实现结构有限元模型气动载荷加载和温度场赋值。采用移动坐标系和动网格相结合的方式描述变速度飞行和翼面偏转过程。通过坐标系变换将翼面偏转过程和振动过程的网格变形量叠加,考虑翼面振动和偏转的耦合非定常效应。针对沿轨道运动的高超声速飞行器,建立了同步计算方法与全时域耦合分析方法相结合的热气动弹性稳定性分析流程。研究发现,与同步计算方法相比,全时域耦合分析方法能够模拟结构振动对流场和结构温度场的影响,计算得到的监测点热流密度波动幅值占热流峰值的10%左右,而温度变化并不明显,相比于刚体模型,监测点温度只下降了0.3%左右。全时域耦合热气动弹性分析方法得到的颤振临界点在4-5号状态点之间,颤振形式为铰链扭转模态与一阶弯曲模态的耦合颤振,与"冻结"模态的热颤振方法结果一致。     

基于分层思路的动态非线性气动力建模方法

很多非线性气动力模型难以精确预测系统的小扰动线性特征。针对这一局限,提出了一种非线性分层模型,用于辨识跨声速非线性非定常气动力。分层建模需要同时提供微幅振荡和大幅振荡两套训练样本,首先通过线性模型(如带外输入的自回归(ARX)模型)对微幅振荡样本进行建模,而后采用非线性模型(如径向基函数神经网络(RBFNN))对大幅振荡的样本与线性模型的差量进行建模,进而把线性模型和非线性模型的输出相叠加,得到分层非线性动力学模型。算例表明建立的分层气动力模型与单一自回归径向基函数(AR-RBF)神经网络模型相比不仅具有更高的数值精度,可以精确预测大幅运动中的非线性特征,而且在小扰动状态下自动退化为线性模型,能够精确刻画结构微幅振荡下的线性动力学特性。     

高温燃气进入舵机舱过程仿真与流动机理分析

采用计算流体力学(CFD)方法研究了火箭发动机工作拖尾段高温发动机燃气进入舵机舱的物理现象。结合导弹实际飞行弹道参数变化特点和超声速流场扰动不向前传递的空气动力学理论,提出了简化而不失真的非定常流场仿真方案,显著缩短了仿真周期;复现了某型导弹实际飞行时舵机舱先被"抽气"再进高温燃气的动态过程,并分析了高温发动机燃气进入舵机舱的流动机理,即在发动机工作段,导弹底端面压强低于舵机舱内压强,舵机舱被"抽气",在拖尾段随着燃烧室总压降低,喷口附近的马赫盘向导弹底端面移动,使导弹底端面压强增大且高于舵机舱内压强,高温燃气进入舵机舱烧毁电路致使导弹折断;明确了某型导弹折断故障产生的诱因,提出了改进措施和检测方法,并得到了大量飞行靶试的验证,解决了舵机舱热防护结构可靠性问题。     

突然启动流动问题:从不可压到高超声速流动

突然启动流动问题在气动弹性、扑翼飞行和机动飞行中有重要应用价值。鉴于此,对突然启动流动问题涉及的流动与升力演化机制和理论分析方法进行统一介绍并指出目前的理论空白。比较性地介绍了升力随时间变化的原因、预测升力演化的理论方法和主要变化规律。从中发现,不可压缩和可压缩突然启动问题存在各自独立的研究方法与结论,因此进行统一介绍与分析对建立二者之间的关联有指导意义。同时指出,大迎角可压缩突然启动问题尚无理论分析方法和研究结果。作为一项补充研究,采用数值计算发现一个之前尚未报道的现象,即升力系数首先从较低值快速增加至一个峰值,接着快速下降,趋于定常值。该文综述介绍的方法可用于分析突然启动问题升力演化规律。     

对转开式转子非定常气动干扰特性分析

采用动态面搭接技术,求解非定常雷诺平均Navier-Stokes方程实现对转开式转子的非定常气动数值模拟.对10片桨叶单个转子和10×10对转开式转子构型分别进行模拟,对比分析前后两个转子间的气动干扰和滑流流动干扰.结果表明:与单个转子相比,对转开式转子前转子的拉力系数和功率系数减小,后转子的拉力系数和功率系数则都增大,并在一个旋转周期内都呈现20次周期性波动.前后转子拉力系数的频谱分析显示振荡发生在桨叶通过频率的偶数倍,且在2倍时后转子拉力的振幅最大.前转子桨尖涡与后转子的周期性干扰,引起后转子桨叶力分布的变化和非定常性.与单个转子相比,前转子后滑流轴向速度偏大,周向速度偏小.后转子对滑流有二次加速及旋涡恢复作用.      

地面效应中垂直起降状态旋翼的气动特性计算

基于时间步进自由尾迹和地面面元模型,构建了一种新的非定常空气动力学模型,用于计算有地效垂直飞行状态下的旋翼气动性能。为提高瞬态尾迹结构的求解稳定性和计算速度,本文使用了精度阶数更高的显式CB3D格式。计算证实,该格式能有效控制地面干扰所产生的数值误差,得到的尾迹结构也更符合物理实际。同时引入一种适用于非定常状态的"等体积修正"方法,以对计算中容易出现的涡线落入地面下方的非物理现象进行修正。在有地效状态下,将新构建模型的计算结果同试验数据进行了对比,验证了模型的准确性。进而计算了平面与斜面上有地效垂直起降状态下旋翼气动特性的变化特点,分析了其影响因素。