颤振边界预测的系统稳定性分析方法

对颤振边界预测方法中的颤振裕度法与自回归滑动平均模型（ARMA）稳定性分析方法进行研究,采用数值仿真算例研究这两种方法受参数识别误差的影响,采用风洞颤振试验算例对比分析其预测颤振边界的有效性与准确性,结果表明:①颤振裕度法的判据受阻尼比识别误差的影响较小,相对于阻尼比识别误差更容易受频率识别误差的影响,ARMA稳定性分析方法的判据相对于颤振裕度法更易受阻尼比识别误差的影响,颤振裕度法比ARMA稳定性分析方法鲁棒性更高;②对于机翼弯扭耦合颤振,两方法的判据变化趋势相近,相差采样频率4次方的数量级,两者随着风速的增加具有平缓的下降趋势,有助于较早地预测颤振边界;③对于面内弯曲为主的颤振,ARMA稳定性分析方法比颤振裕度法适用性更好。      

受限混合层的线性稳定性

针对超声速边界层/混合层组合流动,利用可压缩线性稳定性理论研究了流动的线性失稳特性。基本流场选取了具有不同速度特征的两股来流,采用双曲正切的混合层剖面叠加可压缩边界层自相似性解剖面构造。重点考察了混合层中心与壁面距离、对流马赫数等参数对组合流动稳定性特征的影响,其中壁面采用绝热壁面。混合层中心与壁面的距离为5~15倍的边界层厚度,混合层的对流马赫数为0.6~1.2。结果表明:该组合流动中存在独特的多重不稳定模态,并相互影响;且其不稳定模态随着壁面距离及对流马赫数的变化呈现出不同的主导行为。      

下单翼布局飞机不同整流减阻特性研究

下单翼常规布局飞机的机翼机身结合处存在严重的相互干扰,该处的附面层堆积是造成气流分离的主要原因,分离会致使飞机的阻力增加。通过整流将机翼机身结合处的附面层外推,在远前方来流的作用下吹除,通过减小附面层厚度减小阻力;对比分析三种不同整流方案的特点,并进行数值分析。结果表明:本文的三种整流减阻方案能够达到整流减阻的效果,改善了气动性能,可为下单翼常规布局飞机整流减阻提供参考。     

高超声速二维钝楔边界层参数分析

为了指导二元高超声速进气道型面设计,基于两类高超声速二维钝楔外形,采用高精度CFD方法计算了四种典型工况下二维钝楔边界层参数,分析了前缘半径、壁面温度、压力梯度对边界层参数的影响,得出了边界层位移厚度、边界层总压恢复及壁面压力梯度沿流向的变化规律。研究表明,前缘半径及壁面温度对边界层参数影响较大,在两者一定的情况下,逆压梯度使边界层位移厚度增速放缓,顺压梯度会使边界层位移厚度快速增长。     

无网格法耦合RNGk-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

 基于移动最小二乘无网格方法,耦合RNG (Re-Normalisation Group) k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM (Advection Upstream Splitting Method)⁺-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP (Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNG k-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。     

利用涡流发生器抑制S弯进气道旋流畸变的数值模拟研究

S弯进气道由于其良好的隐身性能为现代作战飞机所青睐,但是由于进气道是弯曲的,较大的表面曲率容易导致边界层分离,降低了进气道的总压恢复,增加了气流的湍流脉动,导致旋流畸变的发生。对边界层流动实施控制可以推迟、扼制流动分离,改变进气道的气动性能。本研究在某机型进气道的CAD模型的基础上,利用FLUENT软件对采用不同类型的边界层扰流片(涡流发生器)进行数值模拟,控制边界层分离,使大S弯进气道出口的总压恢复和旋流畸变指标得以改善。     

超声速燃烧数值模拟中的湍流与化学反应相互作用模型

高精度数值模拟有助于理解超声速湍流燃烧中湍流与化学反应的相互作用,可为发动机燃烧室等工程应用设计提供可靠的预测模型。除直接数值模拟外,目前在湍流燃烧应用中使用的大涡模拟和雷诺平均Navier-Stokes模拟均需要借助模型模化发生在湍流小尺度上的流动与化学反应过程对湍流大尺度运动的影响。现有的湍流与化学反应相互作用模型大致可分为:火焰面类模型和概率密度函数类模型,2类模型在不同的应用中各自具有优势和局限性。此外,现有模型大都基于低马赫数燃烧,而超声速燃烧中通常会伴随快速混合、局部熄火和再着火以及激波等复杂过程,这为发展其中的湍流与化学反应相互作用模型提出了更多的挑战。     

Mechanism of unconventional aerodynamic characteristics of an elliptic airfoil

Abstract The aerodynamic characteristics of elliptic airfoil are quite different from the case of conventional airfoil for Reynolds number varying from about 10~4to 10~6.In order to reveal the fundamental mechanism,the unsteady flow around a stationary two-dimensional elliptic airfoil with 16%relative thickness has been simulated using unsteady Reynolds-averaged Navier–Stokes equations and the γ-Reθt transition turbulence model at different angles of attack for flow Reynolds number of 5×10~5.The aerodynamic coefficients and the pressure distribution obtained by computation are in good agreement with experimental data,which indicates that the numerical method works well.Through this study,the mechanism of the unconventional aerodynamic characteristics of airfoil is analyzed and discussed based on the computational predictions coupled with the wind tunnel results.It is considered that the boundary layer transition at the leading edge and the unsteady flow separation vortices at the trailing edge are the causes of the case.Furthermore,a valuable insight into the physics of how the flow behavior affects the elliptic airfoil’s aerodynamics is provided.     

超声速H_2/Air湍流扩散燃烧RANS数值模拟

为了研究修正的火焰面反应进度变量燃烧模型在超声速湍流扩散燃烧问题中的适用性,对德国宇航中心(DLR)超声速燃烧室开展RANS数值模拟。基于Open Foam软件平台中密度求解器分别对三维冷态场和燃烧场进行模拟分析。将网格自适应加密技术用于流场的计算;燃烧场计算中,通过分析不同压力下层流火焰面数据库,引入了反应进度变量源项的压力修正系数,压力修正系数α等于2.2。计算结果表明,冷态场中压力分布、波系分布、速度分布以及燃烧场中波系分布、速度分布、温度分布结果均与实验值符合较好。压力修正方法能够较好地解决超声速湍流扩散燃烧问题。湍流Schmidt数敏感性分析表明,湍流Schmidt数Sct对湍流火焰结构有较大影响,文中Sct等于0.7时能得到与实验值较为一致的分布。     

航空发动机主轴承热分析边界条件处理方法

为了提高航空发动机主推力球轴承热分析的计算精度,对轴承的摩擦发热和对流换热边界条件进行了分类及研究。应用ANSYS有限元分析软件,采用将摩擦热按体积生热率处理和将摩擦热按热流密度处理的2种不同方式,对边界条件进行了加载,分别对试验器状态的发动机主轴承进行了热分析计算,并与试验测量结果进行了对比。计算结果表明:采用表面效应单元加载热流密度的方式得到的轴承温度分布更理想,内部热点温度更集中,热点温度比按体积生热率加载的高。2种边界条件处理方法均已应用到航空发动机润滑系统热分析中,提高了航空发动机润滑系统热分析的准确性。