基于欧拉方程的跨声速翼型设计

给出了一种基于欧拉方程的跨声速翼型设计方法。方法以Ｔａｋａｎａｓｈｉ提出的“正反迭代余量修正”设计原理为基础，在气动力分析模块中，以欧拉方程为基本控制方程，并采用Ｗａｌｚ边界层方法对其进行粘性修正；反设计模块采用经过改进的二维翼型设计方法。方法的程序经过几个设计实例主题证明是有效而实用的。     

关于三维Navier－Stokes方程的粘性项计算

借助于张量分析和张量计算，在贴体曲线坐标系下本文讨论了不同求解变量导致了粘性项个数上的重大差异和不同大小的计算量，并提出了便于粘性计算的最佳形式。文中借助于有限体积离散技术，通过引进两个对称辅助矩阵［Ａ］和［Ｂ］，使粘性项的计算量大大减少，这对完成三维粘性流的数值计算具有重要的指导意义。借助于上述方法，本文完成了某型真实进气道两种工况的三维粘性Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程计算（即Ｍ∞＝３．０，α＝０°和设计工况Ｍ∞＝２．６５，α＝０°），获得了满意的全场结果；对于Ｍ∞＝２．６５的设计工况，同实验数据作了比较，符合良好。由于本文的方法明显的减少了粘性项的计算量且节省了大量内存，以致于使三维流场的Ｎ－Ｓ求解能在普通微机上进行。     

应用流固耦合数值方法研究机翼的射流减振技术

采用流固耦合的数值计算方法研究了NACA 0 0 15翼型在大迎角 ( 15°～ 6 0°)范围的颤振 ,以及在翼型背部部分引入射流的减振技术。在流体区域用高精度、高分辨率算法求解Farve平均的Navier Stokes方程 ,在固体区域用四阶Runge Kutta法求解振动方程 ,并且每一个时间步后都在两个区域之间传递边界条件。计算结果表明在翼型背部引入适当射流会降低翼型的振动 ,并提高升力。但如果引入射流的速度过高 ,会在叶背处形成新的分离流 ,升力反而会下降     

翼型动失速和非定常载荷的综合翼型数据方法

本文给出受迫谐振翼型的动失速工程估算方法。本方法基于凤洞试验，综合分析翼型动、稳态特性之间的差别与动、稳态条件之间的关系，建立一套估算动、稳态特性之间差别的经验公式，修正稳态特性，得到相应的动态特性。用本方法计算了三种翼型不同动态条件（包括深失速和后掠）的动失速特性，并与测量和文献结果进行了比较，结果符合得相当好。     

求解大型对称特征值问题的迭代Chebyshev-Lanczos方法

本文研究计算大型对称矩阵极端(几个最大或最小)特征值及相应特征向量的问题,讨论了Chebyshev迭代法对Lanczos方法的应用,提出了Chebyshev-Lanczos方法。计算实践表明迭代Chebyshev-Lanczos方法比迭代Lanczos方法优越。     

关于Davidson—Lanczos方法的收敛率

本文对文[1]提出的求解大型对称矩阵A的极端(几个最大或最小)特征值及相应特征向量的Davidson—Lanczos方法,用Rayleigh—Ritz逼近理论,研究了该方法的收敛率。证明了由该方法产生的规范正交向量{v_i}_i~m=1是Krylov子空间K_m≡Span(v_1,Av_1,…,A~(m-1)v_1)的一组基。设A的k个最大特征值为又,λ_1>λ_2>…>λ_k,相应的近似特征值为λ_i~(m)(i=1,…,k),得到 这里γ_i(γ_i>1),W_i和W_i~(m)是常数。     

Transonic buffet control research with two types of shock control bump based on RAE2822 airfoil

Current research shows that the traditional shock control bump (SCB) can weaken the intensity of shock and better the transonic buffet performance.The author finds that when SCB is placed downstream of the shock,it can decrease the adverse pressure gradient.This may prevent the shock foot separation bubble to merge with the trailing edge separation and finally improve the buffet performance.Based on RAE2822 airfoil,two types of SCB are designed according to the two different mechanisms.By using Reynolds-averaged Navier-Stokes (RANS) and unsteady Reynolds-averaged Navier-Stokes (URANS) methods to analyze the properties of RAE2822 airfoil with and without SCB,the results show that the downstream SCB can better the buffet performance under a wide range of freestream Mach number and the steady aerodynamics characteristic is similar to that of RAE2822 airfoil.The traditional SCB can only weaken the intensity of the shock under the design condition.Under the off-design conditions,the SCB does not do much to or even worsen the buffet performance.Indeed,the use of backward bump can flatten the leeward side of the airfoil,and this is similar to the mechanism that supercritical airfoil can weaken the recompression of shock wave.     

Aerodynamics of smart flap under ground effect

Aerodynamic characteristics of two-dimensional smart flap under the ground effect have been assessed by a numerical simulation. In this process, a pressure-based implicit procedure to solve Navier–Stokes equations on a nonorthogonal mesh with collocated finite volume formulation is used. The boundedness criteria for this procedure are determined from the Normalized Variable Diagram (NVD) scheme. The procedure incorporates the k–ε eddy–viscosity turbulence model. Cantilever beam with uniformly varying load with roller support at the free end is considered for the configuration of the smart flap. The method is first validated against experimental data. Then, the algorithm is applied for turbulent aerodynamic flows around airfoil with smart and conventional flaps for different attack angle, flap angle and ground clearance where the results of two flaps are compared. The comparisons show that the quality of the solution is considerable.     

常气压辉光放电等离子体控制翼型失速的数值模拟研究

基于Shyy提出的大气压下均匀辉光放电等离子体与空气干扰的物理模型,通过求解电位势方程得到电场分布,得到了作用于流体上的电场力.通过数值求解考虑等离子体作用的流体运动控制方程,以NACA0015翼型低速绕流为对象,研究了常压辉光放电等离子体位置和个数控制翼型绕流分离的影响.位于分离点上游的等离子体能够有效地抑止翼型分离,而在分离区的等离子体对流动影响很弱,同实验结果一致,并给出了等离子体对翼型气动力影响的规律.     

结冰翼型表面明冰的压力分布和应力计算

在不同来流速度和攻角下,比较分析了冰型外部压力分布的变化规律,计算分析了外部压力分布对表面冰型的力学特性影响规律,研究了冰型和蒙皮界面间最大剪切应力的变化规律,并结合剪切强度随界面温度变化的实验曲线,分析了最大剪切应力对冰层脱落的贡献.计算结果表明:最大剪切应力是剪切强度的5.48%,因此气动力对冰层的脱落贡献不大;当开启热除冰装置后,由于蒙皮表面温度的不断上升,冰型和蒙皮界面间剪切强度的下降,当表面温度为-2℃时,最大剪切应力超过剪切强度的10%,此位置可能率先产生破裂,从而加速冰层的脱落.