高超声速激波/边界层干扰流动数值模拟研究

在基于Roe格式的全Navier-Stokes方程计算流体力学(CFD)代码中,发展了一种局部熵修正方法,克服了传统熵修正方法在边界层流动模拟中耗散过大的缺点,可用于更加准确的模拟激波/边界层干扰的复杂高超声速流动。对典型高超声速双锥边界层分离与激波干扰的复杂流动进行了模拟,研究了网格密度和熵修正方法耗散性对计算结果的影响。研究表明:高超声速双锥边界层分离与激波干扰流动的数值模拟结果对网格具有很强的敏感性,过稀的网格将产生严重的分离流动预测偏差;低耗散性的局部熵修正方法能更加准确地模拟复杂的高超声速激波与边界层分离流动干扰现象。     

高性能大型金属构件激光增材制造:若干材料基础问题

简要介绍了高性能大型金属构件激光增材制造的技术特点、国内外研究进展及技术发展面临的挑战,分析了大型金属构件激光增材制造的"高性能材料制备"与"复杂结构直接制造"有机融合、"控形/控性"一体化的独特特征。指出高性能大型关键金属构件激光增材制造技术的发展和工程应用,将在很大程度上取决于人们对激光增材制造过程中对激光/金属交互作用行为及能量吸收利用机制、内部冶金缺陷形成机制及力学行为、移动熔池约束快速凝固行为及构件晶粒形态演化规律、非稳态循环固态相变行为及显微组织形成规律、内应力演化规律及构件变形开裂预防方法等材料基础问题的深入研究。     

跨声速压气机转子叶尖流场旋转不稳定现象的数值研究

对某跨声速压气机转子在不同工作流量下的叶尖非定常流场进行了数值研究.结果显示:大流量状态下,该转子叶尖流场几乎不发生振荡.此时,叶尖流场可以按定常流场进行分析;小流量状态下,叶尖泄漏涡大幅振荡,相邻叶片通道内的叶尖泄漏流之间也存在周期性相互干涉.其结果是在稳定状态时出现由于叶尖泄漏涡的振荡及其周向传播造成的"旋转不稳定"现象."旋转不稳定"流场结构主模态旋涡个数大约为40%的叶片通道个数;其周向尺度占据2~3个栅距.     

A study on turbulence transportation and modification of Spalart–Allmaras model for shock-wave/turbulent boundary layer interaction flow

It is of great significance to improve the accuracy of turbulence models in shock-wave/ boundary layer interaction flow. The relationship between the pressure gradient, as well as the shear layer, and the development of turbulent kinetic energy in impinging shock-wave/turbulent boundary layer interaction flow at Mach 2.25 is analyzed based on the data of direct numerical simulation(DNS). It is found that the turbulent kinetic energy is amplified by strong shear in the separation zone and the adverse pressure gradient near the separation point. The pressure gradient was non-dimensionalised with local density, velocity, and viscosity. Spalart–Allmaras(S–A) model is modified by introducing the non-dimensional pressure gradient into the production term of the eddy viscosity transportation equation. Simulation results show that the production and dissipation of eddy viscosity are strongly enhanced by the modification of S–A model. Compared with DNS and experimental data, the wall pressure and the wall skin friction coefficient as well as the velocity profile of the modified S–A model are obviously improved. Thus it can be concluded that the modification of S–A model with the pressure gradient can improve the predictive accuracy for simulating the shock-wave/turbulent boundary layer interaction.     

基于变步长CFD/RBD方法的旋转弹轨迹仿真

采用变步长方法实现非定常Navier-Stokes方程和刚体六自由度运动方程的耦合求解,并对美国ARL旋转弹试验模型的飞行轨迹进行了仿真。旋转弹自由飞行过程中受到空气阻力作用,转速下降,气动力变化趋缓。因此轨迹仿真时先选取较小的时间步长,随着转速的下降逐渐增大步长,这样可以减少计算时间。刚体六自由度运动方程的数值求解采用变步长Adams预估校正法。将定步长和变步长两种方法应用于旋转弹实验模型轨迹仿真,并将计算得到的旋转弹空间位置、姿态角与试验结果进行了对比,结果显示:方法能较精确地仿真旋转弹飞行轨迹。在相同初始步长条件下,变步长方法较定步长方法计算时间减半;在相同计算时间条件下,变步长方法较定步长方法精度更高。     

开裂角对阻力方向舵颤振特性的影响

对于无垂尾飞翼式布局飞机,阻力方向舵是一种十分有效的偏航控制装置。偏航控制时,舵面大角度偏转,引起气流分离、涡等复杂流场运动,非定常气动力复杂。计算气动弹性学科提出流固耦合求解方法,以期用于阻力方向舵气动弹性问题的求解。本文采用基于CFD技术的流固耦合方法求解阻力方向舵二维气动弹性问题,计算结果表明,随着开裂角的增大,阻力方向舵的颤振速度增加。对阻力方向舵气动特性进行了计算分析,结果表明阻力方向舵开裂,舵面背风区形成死水区,舵面非定常气动力影响系数减小,阻力方向舵开裂角越大,其颤振速度越大。     

弹性飞机阵风减缓研究

阵风响应减缓可有效抑制飞机由于离散突风和紊流引起的附加过载。以某民机为研究对象,基于飞机结构动力学方程和非定常空气动力学模型,建立弹性飞机在突风作用和控制面偏转联合激励下的开环气动弹性方程。引入坐标变换,建立飞机刚体模态对应的广义坐标与飞行动力学运动变量之间的关系,从而实现气动弹性方程与飞行动力学方程的融合。并以副翼和升降舵为突风减缓控制面。采用俯仰角速率,翼尖加速度和质心加速度作为反馈信号,针对连续突风和离散突风,分别设计相应的控制律。仿真结果表明,减缓系统对翼尖过载、翼根弯矩具有较好的抑制能力。     

高压涡轮全环非定常流动数值模拟

为了进一步准确地掌握涡轮内部复杂流动特征,采用全环非定常数值模拟方法研究了某型高压涡轮内部非定常流动,重点分析了上游导叶尾迹和激波在动叶流动中的非定常作用机制,探讨了引起动叶通道内及叶片表面压力分布周期性变化的非定常扰动。数值结果表明:与尾迹作用相比,由导叶激波引起的非定常效应更为明显,尾迹引起的非定常扰动在动叶0.5倍轴向弦长位置下游较为突出。在动叶流场中及叶片表面上的静压变化显示,在动叶通道内还存在低频扰动,对局部流场的非定常扰动明显。最后指出,在非定常计算时采用区域缩放法对流场中非定常扰动的预测存在一定误差。     

Contribution of geometry of interaction between interplanetary and terrestrial magnetic fields into global magnetospheric state and geomagnetic activity

The paper presents results of our study of dependence of geomagnetic activity from geoeffective parameters taking into account mutual orientation of the interplanetary magnetic field, electric field of the solar wind and geomagnetic moment. We attract a reconnection model elaborated by us made allowance for changes of geometry of the solar wind–magnetosphere interaction during annual and diurnal motions of the Earth. We take as our data base the interplanetary magnetic field and solar wind velocity measured at 1 a.u. at ecliptic plane for the period of 1963–2005 and K_p, D_st, a_m indices. Taken as a whole a geoeffective parameter suggested by us explains 95% of observed variations of the indices. Changes of the geometric factor determined by mutual orientation of the solar wind electric field and geomagnetic moment explain larger than 75% of observed statistical variations of D_st and a_m indices. Based on our results we suggest a new explanation of semi-annual and UT variation of geomagnetic activity.     

风力机翼型阵风响应计算与分析

基于网格速度方法,对阵风条件下风力机翼型运动的气动特性进行数值计算。中心格式有限体积法、双时间推进法求解预处理后的非定常N-S方程。计算结果表明:在弦向阵风作用下,翼型的运动能有效降低响应波动峰值,即弦向阵风对运动影响较小;而在法向阵风作用下,气动响应波动很大。计算翼型在不同参数阵风作用下的响应过程,阵风的频率变化对响应影响不大。弦向阵风的幅值增大,叶片运动对阵风响应峰值减缓效果增加;法向阵风的幅值增大1倍,响应峰值也增加1倍。