高超声速乘波飞行器气动特性研究

用计算流体力学和风洞试验的方法对以锥导乘波体为基础生成的高超声速乘波飞行器的气动性能进行了研究。结果表明，以马赫数6，攻角4度为设计状态的乘波体，在马赫数5～7，攻角4～6度的范围内，都具有良好的气动特性，升阻比接近4。最后，提出了一个简单的以参考温度方法为基础的粘性阻力分析方法。该方法配合使用风洞试验和计算流体的结果，可以用来验证计算流体中难以计算准确的粘性阻力，也可以用来分析在风洞试验难以直接得到的粘性阻力。对于工程上的粘性阻力分析是一个有用的办法。     

基于灰色关联聚类的空中交通复杂度评估方法

为了准确评价空中交通态势的复杂程度,进行有效的空中交通管理,提出了一种空中交通复杂度评估方法。在已有的理论研究成果及实际运行经验的基础上,归纳了定义明确、可定量分析的基础空中交通指标集,使用灰色关联聚类方法对指标集进行聚类分析,利用信息熵理论对聚类指标求解权重,加权求和获得空中交通复杂度。方法弥补了目前使用单一维度指标描述空中交通运行态势的弊端,能更好地体现空中交通复杂度。     

Numerical analysis on noise of rotor with unconventional blade tips based on CFD/Kirchhoff method

A solver is developed aiming at efficiently predicting rotor noise in hover and forward flight. In this solver, the nonlinear near-field solutions are calculated by a hybrid approach which includes the Navier-Stokes and Euler equations based on a moving-embedded grid system and adaptive grid methodology. A combination of the third-order upwind scheme and flux-difference splitting scheme, instead of the second-order center-difference scheme which may cause larger wake dissipation, has been employed in the present computational fluid dynamics (CFD) method. The sound pressure data in the near field can be calculated directly by solving the Navier-Stokes equations, and the sound propagation can be predicted by the Kirchhoff method. A harmonic expansion approach is presented for rotor far-field noise prediction, which gives an analytical expression for the integral function in the Kirchhoff formula. As a result, the interpolation process is simplified and the efficiency and accuracy of the interpolation are improved. Then, the high-speed impulsive (HIS) noise of a helicopter rotor at different tip Mach numbers and on different observers is calculated and analyzed in hover and forward flight, which shows a highly directional characteristic of the rotor HIS noise with a maximum value in the rotor plane, and the HSI noise weakens rapidly with the increasing of the directivity angle. In order to investigate the effects of the rotor blade-tip shape on its aeroacoustic characteristics, four kinds of blade tips are designed and their noise characteristics have been simulated. At last, a new unconventional CLOR-II blade tip has been designed, and the noise characteristics of the presented CLOR-II model rotor have been simulated and measured compared to the reference rotors with a rectangular or swept-back platform blade tip. The results demonstrate that the unconventional CLOR-II blade tip can significantly reduce the HSI noise of a rotor.     

基于MEDLL的分级搜索抗多径算法

以多径估计延迟锁定环路(MEDLL)为代表的参量式基带抗多径算法以其能够同时估计直达信号和多径信号各参数而受到广泛关注,然而该类算法估计精度的提高是以计算复杂度的增加为代价的。为了在保证估计精度的同时降低资源消耗并加快搜索速度,提出了一种MEDLL最大似然(ML)估计搜索方法,将粗搜索后得到的参数估计值作为下一级搜索的先验信息,缩小搜索范围后减小搜索步进,完成逐级精细的参数估计。推导了原始MEDLL及树状分级结构的MEDLL(T-MEDLL)算法计算复杂度表达式;其中在存在一路多径信号的情况下,对于0.01 chip的多径延迟分辨率,复杂度可以降低为传统盲搜索的30%~50%。进一步基于自研的数字中频(IF)GPS信号模拟器和软件GPS接收机平台对不同多径情况下盲搜索法与分级搜索算法的参数估计精度及码伪距多径误差等性能进行了对比,验证了更低计算复杂度的T-MEDLL可以达到与传统MEDLL算法同样的跟踪精度。      

涡轮叶片多重精度MDO方法

考虑到多学科优化设计(multidisciplinary design optimization,MDO)的效率和精度,充分发挥当前涡轮叶片各种精度分析方法的优势,提出了涡轮叶片的多重精度多学科优化设计策略.利用协同优化策略在学科解耦和协调方面的能力,引入包括流固紧密耦合分析、流固热松散耦合分析和近似方程在内的涡轮叶片设计中的多种精度模型,结合两点式标度函数和优化过程阶段性收敛后更新,改进可变复杂度建模方法处理多重精度模型的能力,研究结果表明:仅调用9次高精度流固耦合分析就能确保涡轮带冠叶片多学科优化设计的效率和精度.      

某涡轮级三维流场非定常数值研究

基于耦合Lantry-Menter转捩模型的SST(Shear stress transport)湍流模型,对某涡轮级在3个雷诺数工况下的全三维流场分别进行了定常、非定常数值模拟,目的在于对比分析两种计算方法所得结果之间的差异。计算结果表明,进口雷诺数为4×104时和5×105时,定常与非定常计算的涡轮级效率的差异较大,最大相差0.6个百分点,而当雷诺数为8×104时,两者差异很小;非定常计算能够较好的模拟叶尖泄漏涡的发展过程、低雷诺数下上游尾迹与下游转子叶片吸力面分离边界层干涉诱导卷起涡的形成、输运及其所导致的叶片吸力面表面压力大幅波动等非定常流动现象。最后,给出了转子出口截面上周向与径向的非定常压力脉动的分布,为基于Lighthill声类比方法的计算气动声学研究奠定了基础。     

A transient bulk flow model with circular whirl motion for rotordynamic coefficients of annular seals

Bulk flow model with perturbation simplification has been used to calculate rotordynamic coefficients in annular seals which have significant influences on the dynamic behavior of rotors in turbomachinery. In this work, a transient bulk flow model with arbitrary rotor motion is developed, and the boundary conditions and friction factor in the model are calibrated with steady Computational Fluid Dynamics (CFD) analysis. The numerical solution scheme is developed based on the finite element method to obtain the transient reaction force in the seal clearance. With a periodic circular rotor orbit, the transient forces at multiple whirling frequencies are used to evaluate the rotordynamic coefficients. The leakage flowrate of CFD analysis has good agreement with experimental results and the calibrated parameters in bulk flow model are dependent on operating conditions. Although CFD calibration improves the accuracy of the perturbed bulk flow model, the direct damping is overestimated and the cross-coupled damping is underestimated. Compared with the perturbed model, the predictions of the transient bulk flow model are more agreeable with the experiment.     

Interface flux reconstruction method based on optimized weight essentially non-oscillatory scheme

Aimed at the computational aeroacoustics multi-scale problem of complex configurations discretized with multi-size mesh, the flux reconstruction method based on modified Weight Essentially Non-Oscillatory (WENO) scheme is proposed at the interfaces of multi-block grids. With the idea of Dispersion-Relation-Preserving (DRP) scheme, different weight coefficients are obtained by optimization, so that it is in WENO schemes with various characteristics of dispersion and dissipation. On the basis, hybrid flux vector splitting method is utilized to intelligently judge the amplitude of the gap between grid interfaces. After the simulation and analysis of 1D convection equation with different initial conditions, modified WENO scheme is proved to be able to independently distinguish the gap amplitude and generate corresponding dissipation according to the grid resolution. Using the idea of flux reconstruction at grid interfaces, modified WENO scheme with increasing dissipation is applied at grid points, while DRP scheme with low dispersion and dissipation is applied at the inner part of grids. Moreover, Gauss impulse spread and periodic point sound source flow among three cylinders with multi-scale grids are carried out. The results show that the flux reconstruction method at grid interfaces is capable of dealing with Computational AeroAcoustics (CAA) multi-scale problems.     

Numerical simulation of broken pin effects on the flow field and cooling performance of a double-wall cooling configuration

In this paper, the flow characteristics of the double wall structure are presented and the effect of the broken pin size on the cooling performance and flow field of the double wall configuration is investigated. A periodic plate model with seven units is adopted, and there are an impingement hole and a film hole in each unit. Under five blowing ratios, six different sizes of the broken pin are compared, and the double wall configuration without broken pins is taken as the baseline. The results show that if the broken pins height is too small, the cooling effectiveness usually cannot be improved. With the presence of broken pins with a height of more than 0.4, the effectiveness is improved due to the enhancement of reattachment and recirculation of coolant. With the increase of the broken pin height, the cooling effectiveness increases. However, the increase of the diameter does not always improve the cooling performance, since the limiting effect of the wall jet. In this study, Case 6 with the largest broken pin always has the best cooling performance, but also the largest flow resistance. In Case 6 temperature is reduced by almost 15 K compared to the baseline, and more areas have relatively higher cooling effectiveness.     

A hybrid multidimensional Riemann solver to couple self-similar method with MULTV method for complex flows

Since proposed, the self-similarity variables based genuinely multidimensional Riemann solver is attracting more attentions due to its high resolution in multidimensional complex flows. However, it needs numerous logical operations in supersonic cases, which limit the method’s applicability in engineering problems greatly. In order to overcome this defect, a hybrid multidimensional Riemann solver, called HMTHS (Hybrid of MulTv and multidimensional HLL scheme based on Self-similar structures), is proposed. It simulates the strongly interacting zone by adopting the MHLLES (Multidimensional Harten-Lax-van Leer-Eifeldt scheme based on Self-similar structures) scheme at subsonic speeds, which is with a high resolution by considering the second moment in the similarity variables. Also, it adopts the MULTV (Multidimensional Toro and Vasquez) scheme, which is with a high resolution in capturing discontinuities, to simulate the flux at supersonic speeds. Systematic numerical experiments, including both one-dimensional cases and two-dimensional cases, are conducted. One-dimensional moving contact discontinuity case and sod shock tube case suggest that HMTHS can accurately capture one-dimensional expansion waves, shock waves, and linear contact discontinuities. Two-dimensional cases, such as the double Mach reflection case, the supersonic shock / boundary layer interaction case, the hypersonic flow over the cylinder case, and the hypersonic viscous flow over the double-ellipsoid case, indicate that the HMTHS scheme is with a high resolution in simulating multidimensional complex flows. Therefore, it is promising to be widely applied in both scholar and engineering areas.