Numerical simulation of a flow around airfoil using inviscid gas equations

A method for determining a field of inviscid gas equation application in finding the aerodynamic characteristics of bodies is proposed; as an example, use is made of airfoils for which extensive experimental data are available. The comparative analysis of physical and numerical simulation results showed a good verification of the method presented.     

Flow characteristics around airfoils near transonic buffet onset conditions

In transonic flow, buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation. The phenomenon is common in transonic flow, and it has serious impact on the structural strength and fatigue life of aircraft. In this paper, three typical airfoils: the supercritical OAT15A, the high-speed symmetrical NACA64A010, and the thin, transonic/supersonic NACA64A204 are selected as the research objects. The flow fields of these airfoils under pre-buffet and buffet onset conditions are simulated by Unsteady Reynolds Averaged Navier-Stokes (URANS) method, and the mode analysis of numerical results is carried out by Dynamic Mode Decomposition (DMD). Qualitative and quantitative analysis of the shock wave motion, shock wave intensity, shock foot bubble and trailing edge separation, and pressure coefficient fluctuation were performed to attain deep insight of transonic buffet flow features of different airfoils near buffet onset conditions. The results of DMD analysis show that the energy proportion of the steady mode of these airfoils decreases dramatically when approaching the buffet onset angle of attack, while the growth rate of the primary mode increases inversely. It was found that at the onset of buffet, there exist different degrees of merging behavior between shock foot bubble and trailing edge separation during one buffet cycle, and the instability of shock wave and separation induced shear layer are closely related to the merging behavior.     

翼型的气动最优化设计方法和反设计方法＊

本文将数值优化方法同计算流体动力学（ＣＦＤ）相结合,形成两种跨声速翼型的气动设计方法,即：最优化方法和反设计方法。采用求解Ｅｕｌｅｒ方程的有限体积法计算流场,通过结合优化算法,从正反两个方向对跨声速翼型进行气动优化设计。实例证明它们是翼型气动设计的有效方法,有较高的工程应用价值。     

不同热边界对超音速机翼气动热数值分析的影响

飞行器在超音速飞行时受到的气动加热效应给结构强度及热防护设计带来极大影响,且真实状态下的气动热环境需要考虑外流场与结构的耦合及内壁面边界条件的影响。采用S-A湍流模型求解Navier-Stokes方程,通过流场与固体壁面交界处的信息传递,实现外流场与结构场的耦合数值分析。针对三种不同翼型的超音速绕流气动加热进行耦合数值研究,对比翼型内壁面在不同热边界条件下的气动热效应,结果表明：不同翼型具有与气动力相似的的气动热效应;内壁面考虑对流换热的边界条件最接近真实;考虑机翼燃油箱满油时,三维机翼前缘驻点处热流密度最高可达4200w/m2。     

透平跨音速叶栅的优化设计

将遗传算法用于透平跨音速叶栅的优化设计,提出了相应的数学模型和优化方法。优化算例表明,本文所提出的模型及算法能有效的优化叶型方案,避免了一般优化算法的局部最优或维数灾难问题,从而在新的领域内将优化计算在叶轮机械设计中的应用进行了有益的探索。      

旋成面叶栅气动反命题的新解法(Ⅱ)

本文利用映象平面ξ-η的特点, 将旋成面叶栅气动正命题的微分-积分解法与叶型修正的“喷气模型”结合, 建立了一种适用于旋成面叶栅气动反命题的新解法。叶片表面给定的气动参数分布作为反命题的边界条件, 在求解过程中直接得到满足。通过对一例实际透平叶型的计算验证了本文解法的可靠性和实用性。      

重型燃气轮机压气机多截面叶型优化设计

出于对叶型沿叶高方向的不同以及不同叶高处气动边界的不同两方面考虑,开发了一套叶片多截面叶型优化的设计流程.以单级高负荷压气机的进口导叶和大折转角静叶为优化对象,以多截面关联的性能参数为优化目标,采用考虑附面层转捩的MISES计算程序对叶型气动性能进行分析,并通过遗传算法对叶片进行优化.结果显示:在有效工作范围内,各截面叶型经优化后,前部负荷提高,中部及后部负荷降低,性能得到一定提升.三维数值模拟结果显示,在无大尺度二次流动的条件下,多截面叶型优化方法对提升叶片性能切实可行.      

大攻角机翼的气动弹性计算方法

 本文利用大攻角带分离流机翼的完全非定常非线性气动力计算方法,通过与机翼运动方程的同时求解,在时间域内实现了大攻角机翼非线性气动弹性的数值模拟,根据弹性机翼各种状态下的运动过程,可以得到大攻角机翼的颤振速度等重要参数以及亚临界、超临界等飞行状态的运动规律。算例结果表明,大攻角下机翼的气动弹性问题需引起足够的重视。     

基于特征正交分解翼型流场反设计方法研究

基于特征正交分解( Proper Orthogonal Decomposition,POD)方法中模态的思想,利用特征正交分解基,通过改变马赫数或迎角来得到相应的流场快照,对流场参数进行重构,基于快照范围内可以得到高精度的重构流场结果.以NACA0012翼型和RAE2822翼型为研究对象,采用POD方法对流场参数进行了...     

耦合六自由度运动弹性飞机阵风响应数值模拟

耦合飞机刚体六自由度运动,考虑飞机结构气动弹性变形,基于全湍流Navier—Stokes方程流场数值模拟方法．建立了离散阵风作用下弹性飞行器气动载荷与飞行姿态响应的数值模拟技术。其中模态空间中结构动力学方程以及六自由度运动方程分别采用四阶R—K进行时间推进求解,采用RadialBasicFunction（RBF）技术进行气动／结构数据耦合。对应飞机姿态以及弹性变形采用RBF与TFI组合模式的动网格方法进行网格更新;以飞翼布局弹性飞机为数值研究对象,研究了离散阵风作用下飞行器的气动载荷响应以及飞行姿态的变化,对比分析了刚性飞行器与弹性飞行器对离散阵风响应的区别,为进一步系统地分析真实飞机飞行品质提供数值平台。     

基于太阳能飞机应用的低雷诺数翼型研究

以太阳能飞机为背景,对低雷诺数翼型FX 63-137进行了折线型建模以拟合典型晶硅太阳能电池片对气动外形的影响,开展了气动数值模拟分析。首先引用"拟合优度"的定义描述本文折线型翼型轮廓与基准翼型(Baseline)的吻合度,并以此参数为变量建立5种折线型翼型模型;然后,采用计算流体力学(CFD)方法计算分析了不同雷诺数下各折线型翼型的气动特性,并着重研究了低雷诺数下折线型翼型的绕流机理;最后,基于工程应用实际的需求,提出了晶硅太阳能电池片的铺设方法也即折线型翼型设计思想准则,并进行算例验证。研究结果表明:低雷诺数条件下,折线型翼型升阻性能相比光滑翼型在一定程度上表现出了优势,但随着雷诺数的增加,升阻方面的优势逐渐消失;折线型翼型压力分布受各折线段长度影响,前缘吸力峰值、压力平台范围以及压力恢复区分布特征是决定折线型翼型气动性能的主要因素;通过设计的算例验证了本文提出的折线型翼型设计思想的可行性。     

畸变进气下压气机叶片非定常载荷的数值模拟

为了分析畸变进气在轴流压气机转子叶片表面产生的非定常激振力，采用二维无粘线化非定常流理论，建立了周期性气流条件下叶片表面非定常气动力的计算模型，根据实测的畸变实验数据，对某发动机的第一级叶片表面非定激振力进行了计算分析，结果表明在叶片的不同径向位置，非定常压力的分布规律不完全相似，而作用在叶片表面总的非定常气动力在叶尖部位最强。     

大厚度钝后缘翼型气动性能计算研究

大厚度钝后缘翼型由于其气动和结构的优点,近来被作为大型风力机叶片设计时的内侧翼型成为风力机翼型设计的热点之一EULL由于此类翼型厚度和风洞阻塞度的限制,大雷诺数的风洞实验数据很少。耦合RANS方程和基于线性稳定性分析EULL捩预测方法,进行了大厚度钝后缘翼型气动性能计算研究。使用Gridgen程序生成计算网格,湍流模型为S-A模型。对常用的钝后缘修形方式如直接截断、对称增加厚度等并且有风洞实验结果的几种翼型进行了CFD计算并和实验结果及文献计算结果进行了比较。结果表明升力系数在线性段,计算结果和实验结果吻合很好,但是对失速迎角的捕捉能力差;阻力系数计算结果和实验值吻合较好。不同修形方式的计算结果为钝后缘修形设计提供了参考。     

共轴双旋翼非定常流场干扰特性

直升机共轴双旋翼相遇过程的非定常气动干扰在其周期性干扰中最为强烈。为了研究该过程中的干扰特征,建立了一个基于非定常雷诺平均Navier-Stokes方程的气动干扰数值方法,双旋翼的反转运动采用运动嵌套网格进行模拟。为了探究双旋翼相遇时的气动干扰机理,以两个相对运动的双翼型系统来模拟双旋翼特定展向截面的非定常相遇过程。分析了双旋翼和双翼型相遇时的气动特性和流场特征,并对双翼型系统进行了参数影响研究,结果表明:双旋翼上下桨叶相遇时,上下旋翼拉力均会出现先增后减的变化趋势且波动幅值分别为其对应总拉力的30%和22%,双翼型系统的升力波动趋势与之相似;间距增大使上下翼型间的气动干扰减小,且上下翼型升力波动对间距变化的敏感时间不同;翼型相对厚度增大使双翼型升力波动幅值增大,并出现二次波动;有附加来流时,相对来流速度较小的翼型升力波动的幅值与范围更大。      

多段翼型低速流场与气动特性数值计算研究

在结构化网格求解RANS方程的基础上,将嵌套网格技术应用在低速绕流气动特性计算中,并以NLR-7301两段翼型和某运输机5段翼型在起飞构型下对低速绕流流场进行数值模拟,将气动特性及流场特性的数值模拟数据与对接网格计算以及风洞实验数据进行了比较.结果表明,该气动特性计算方法有着较高的精度和可靠性.     

跨声速无粘流绕振动翼型流动的数值模拟

本文从Euler方程出发,利用LU-ADI格式对绕二元振动翼型跨声速流场进行了数值模拟。文章中给出了绕NACA0012和NACA 64A 010翼型不同状态的几组数值结果,在同某些实验和其它文献结果基本吻合的前提下,研究了振动对翼型表面压力、激波位置的非定常效应以及其它振动参数对流场非定常效应的影响。     

自由翼的空气动力特性研究

自由翼具有广泛的应用前景.本文研究了自由翼的翼型的几何特征,表明它们的中弧线都具有反折的特点;用数值计算方法对于同一种平面形状,翼型分别为G0741、MH45、MH49和MH78的四种不同的自由翼的气动特性进行了研究,明确了中弧线反折的程度对气动特性的影响.本文还对自由翼在非均匀气流环境中的气动特性进行了研究.     

基于CST参数化的翼型优化遗传算法研究

用外形建模的CST参数化方法,构建翼型几何可直接利用其后缘角、前缘半径等几何特征,其控制参数更能反映翼型特有的气动敏感性,有助于遗传算法搜索寻优。基于这一参数化方法,结合遗传算法,构造了用于翼型优化的设计方法。算法中,CST控制参数作为设计变量,采用二进制编码,并通过引入精英策略,提高了遗传算法的收敛性能。算法适应度评估涉及的流场求解则采用了基于Jameson有限体积法的Euler方程解算程序。先以NACA0012翼型为例,以其某一已知的表面压力分布为目标,进行了遗传算法的重构运算,给出了重构的翼型几何外形,验证了方法。在此基础上,进行了带约束的跨音速翼型优化设计,给出了升力系数极大化和阻力系数极小化等设计算例,展示出翼型优化设计的效果。     

多目标翼型优化设计基因算法研究

本文把基因算法用于多目标翼型优化设计,针对基因算法的搜索原理和气功外形特点,提出了与之相适应的Ｂｅｚｉｅｒ曲线基因表达,构造了用于优化设计的适应度函数,发展了一种小种群演化方法。针对小种群基因算法中出现的近亲繁殖问题,提出了远亲配对策略、构造了自适应非均匀变异算子,提高了算法的收敛速度。文中最后给出了基于全位势方程求解程序的翼型重构反设计以及单目标和多目标翼型优化算例。     

Modelling and design of wing tip devices at various flight conditions using a databased aerodynamic prediction tool

A design study of wing tip devices at high and low speeds is described. The basis of the design study is an equivalent drag approach containing both aerodynamic drag gain and structural weight penalty. A comprehensive parameter study is carried out using a rapid aerodynamic prediction tool named Lift and Drag Component Analysis (LIDCA). Adding to an available lifting-line method a databased module for airfoil data is employed that uses results of two-dimensional flow simulations by multidimensional interpolation. Detailed validation studies of the method at high lift and high speed have demonstrated good accuracy. RANS computations of the selected wing tip designs confirm the predicted benefits at cruise condition. The results of the most effective wing tip designs are analysed at both flight conditions. Finally, options for improving the performance at take-off are suggested.