给定前缘线平面形状的密切锥乘波体设计方法

乘波体因其高超声速阶段的高升阻比性能成为目前研究的热点,但其本身的诸多性能缺陷限制了其在工程中的实际应用. 密切锥乘波体设计 是目前应用较广的乘波体外形设计方法,具有较高的灵活性和生成效率. 本文以弥补乘波体性能缺陷,提高乘波体设计灵活性为目的, 拓展了密切锥乘波体设计方法,推导设计方法中激波出口型线、流线追踪起始线与平面形状轮廓线之间的几何关系,并使用一个微分方程 组给出了具体的数学表达,奠定了定平面形状乘波体设计的理论基础. 通过介绍此微分方程组的数值求解过程,并分析应用此关系的注意 事项,本文提出了给定前缘线平面形状的密切锥乘波体设计方法. 根据此设计几何关系,以渐变前缘、弯曲前缘和双后掠等为例生成定平 面形状乘波体外形,结合计算流体力学方法分析这几类外形的流场,通过流场分布与设计曲线的比较,说明通过此方法设计得到的乘波体 外形保持了高超声速状态的乘波特性,并可以方便的控制平面形状,为提高乘波体的设计灵活性、改善性能缺陷提供了新的途径.      

高超声速飞行器气动防热新概念研究

传统乘波构型的高超声速飞行器尖锐的前缘存在严重的气动加热问题,而简单的前缘钝化气动防热方法由于造成很大的升阻比损失,难以发挥实质性作用. 引入``人工钝前缘(ABLE)'概念,拟以一种新的思路解决这一矛盾. 通过定义ABLE构型的外形参数,并采用CFD数值计算方法研究了各参数对气动力和气动热特性的影响规律,在流场分析的基础上进行了外形优化,最终得到令人满意的新型高超声速飞行器头部外形,总结了运用ABLE概念进行气动防热的相关设计原则和规律.      

上下反翼对双后掠乘波体低速特性的影响

相比于传统乘波体外形, 双后掠乘波体在保持高超声速良好性能的条件下能够提升乘波体低速气动性能, 但其仍存在低速稳定性不好等缺陷. 本文从密切锥乘波体理论提出给定前缘型线的乘波体设计方法, 通过给定三维前缘型线分别生成具有相同平面投影形状的上反和下反机翼双后掠乘波体. 使用CFD技术评估不同上下反程度外翼乘波体的低速性能, 分析升阻特性以及流场涡结构特点. 选取稳定性判据, 研究上下反翼对纵向和横侧向稳定性的影响. 结果表明, 机翼上下反对乘波体低速升阻特性影响较小; 不同外形均为纵向静不稳定的, 且俯仰力矩变化趋势比较类似, 机翼下反可使气动焦点位置后移, 提升纵向稳定性; 机翼上反有助于提升乘波体的横向静稳定性, 而下反则会下降; 机翼上反可以提升侧向稳定性, 且上反程度越大提升效果越明显; 同时机翼上反使乘波体的偏航动态稳定性有明显提升, 下反则会降低, 影响程度与机翼上下反程度呈正相关. 通过结果分析, 说明通过机翼上下反改善乘波体低速稳定性是可行的, 为乘波体在宽速域高超声速飞行器中的应用拓展了途径.      

蜻蜓滑翔时柔性褶皱前翅气动特性分析

蜻蜓是自然界优秀的飞行家,滑翔是其常见且有效的飞行模式.蜻蜓优异的飞行能力来源于其翅膀的巧妙结构,褶皱是蜻蜓翅膀上最为显著的结构之一,不仅提高了翅膀的刚度,还改变了其气动特性,而飞行过程中柔性翅膀会产生变形是蜻蜓翅膀的另一特性.为揭示蜻蜓在滑翔时,柔性褶皱前翅的变形,探究褶皱和柔性的共同作用对其气动特性的影响,基于逆向工程,依据前人的测量数据和研究成果,通过三维建模软件建立了蜻蜓三维褶皱前翅的计算流体力学(computational fluiddynamics,CFD)模型和计算结构力学(computational structuralmechanics,CSD)模型,并通过模态分析验证了此模型有足够的精度.基于CFD方法和CFD/CSD双向流固耦合计算方法分别对蜻蜓滑翔飞行时刚性和柔性褶皱前翅的气动特性进行了数值模拟,结果表明,柔性褶皱前翅受气动载荷后,翅脉和翅膜产生形变,柔性前翅上下表面压力差相较于刚性前翅减小了,从而其升力和阻力也减小了,而在大攻角时,变形后的前缘脉诱导出比刚性前翅更强的前缘涡.因此在攻角小于10$^\circ$时刚性前翅的气动特性优于柔性前翅,继续增大攻角,柔性前翅的气动特性则优于刚性前翅.前翅受载后气动响应时间短,翅尖的变形最大,仅仅产生了垂直于翅膀所在平面方向上的变形,而没有发生扭转,翼根处受到应力最大,褶皱上凸部分承受蜻蜓滑翔时前翅的主要载荷.      

乘波体压缩面变化对其气动性能影响分析

乘波体是一种利用激波包裹特性获得高升阻比的高速飞行器构型.已有研究中,乘波体气动性能的改善主要依赖于给定源流场条件下的前缘型线优化.本文采用数值优化和计算流体力学模拟为主要手段分析了乘波体压缩面变化对其气动性能的影响,以期有效拓展乘波体的设计空间.主要内容如下:首先给出了一种基于表面局部变形的乘波体设计方法.其次结合运用增量修正参数化方法、计算流体力学分析和微分演化算法构造了乘波体压缩面外形气动优化设计流程,以一种椭圆锥形流场生成的乘波体作为基准构型开展了无黏优化.之后从优化结果中选择升阻比递增的6个典型构型进行前缘钝化处理后,基于N-S方程对其气动性能进行了评估.最后综合依据无黏/黏性计算结果分析了乘波体压缩面变化对其气动性能的影响.结果表明该部分形状的改变对乘波体气动性能影响十分明显,在升力面积不变的条件下,乘波体压缩面形状变化可导致其升阻比出现成倍变化,即使在升力不减条件下,升阻比较基准构型也可获得超过14%的提升.此外,还可导致乘波体相对压心系数出现明显偏移.     

基于实验设计方法的高超声速飞机前缘型线优化分析

为探索前缘线变化对吸气式高超声速飞机气动性能的影响,基于一种旁侧进气布局翼身融合体构型,在飞行马赫数6,攻角4°和高度26 km的巡航飞行条件下,结合运用增量修正参数化设计方法、均匀实验设计方法和计算流体力学模拟,分析了飞行器前缘型线与其升阻力系数及纵向压心等性能参数间的关系.计算结果表明,前缘线形状对飞行器升阻力系数明显高于其对纵向压心影响,设计空间范围内升力系数变化约21.3%,阻力系数变化约31.8%,升阻比变化范围约10.63%,但相对压心变化范围仅为3.87%.在此基础上,通过对典型构型物面压力分布进行分析,发现前缘线形状适当弯曲可利用飞行器下表面侧壁压缩产生的高压气流,利用二者的耦合效应使飞行器获得额外的升力增量.      

Effect of leading edge cut on the aerodynamics of ram‐air parachutes

The effect of the configuration of leading edge cut on the aerodynamic performance of ram‐air parachutes is studied via two‐dimensional flow simulations. The incompressible Reynolds‐averaged Navier–Stokes equations, in primitive variables, are solved using a stabilized finite‐element formulation. The Baldwin–Lomax model is employed for turbulence closure. Flow past an LS(1) 0417 airfoil is investigated for various configurations of the leading edge cut and results are compared with those from a Clark‐Y airfoil section. It is found that the configuration of the leading edge cut affects the lift‐to‐drag ratio (L/D) of the parachute very significantly. The L/D value has strong implications on the flight performance of the parachute. One particular configuration results in a L/D value that is in excess of 25 at 7.5° angle of attack. Results are presented for other angles of attack for this configuration. Copyright © 2004 John Wiley & Sons, Ltd.     

高超声速钝前缘乘波构型优化设计研究

乘波体的高升阻比优势使其在高超声速飞行器设计中极具应用前景. 在实际工程应用中, 为了满足防热要求, 乘波体前缘必须进行钝化处理, 前缘钝化对乘波体气动性能会产生显著影响. 因此, 原始尖前缘最优乘波体并不能保证钝化后仍为最优. 针对这一问题, 首先研究了前缘钝化对不同构型升阻特性的影响程度和作用机理. 结果表明: 前缘钝化会造成乘波体升力小幅度降低, 阻力大幅增加, 升阻比显著降低; 其中钝前缘本身的波阻在阻力增量中起主导作用, 而钝前缘本身的摩阻增加量与物面的摩阻降低量十分接近. 基于上述结果, 提出了一种高效评估钝前缘乘波体气动力的方法, 并结合遗传算法, 开展了直接考虑前缘钝化影响的乘波构型优化设计研究, 获得了钝前缘最优构型. 通过CFD数值模拟对最优构型的气动力特性进行评估, 结果表明: 在不同飞行高度、不同升力和不同钝化半径约束下, 相比尖前缘最优构型, 钝前缘最优构型宽度变窄, 相同纵向位置处的后掠角增大, 且升阻比显著提升. 在M_∞ = 15, H = 50 km, CL = 0.3约束条件下, 钝化半径R = 10 mm的钝前缘最优构型设计点升阻比相比尖前缘最优构型提升量可达9.32%.      

基于神经网络技术的乘波体优化设计

乘波体是高超声速飞行器的主要组成部分,也是飞行器产生升力的主要部分. 针对基 于计算流体动力学(CFD)分析的乘波体优化设计问题,引入人工神经元网络响应面方法. 选 取一定数量的乘波体外形,进行气动性能分析后,利用乘波体的外形控制参数和气动参数做 为训练样本对乘波体进行训练. 利用这些训练样本对人工神经网络进行训练. 在优化计算中 以充分训练的神经网络替代CFD分析,发展了一种基于神经网络技术的乘波体优化设计方法. 利用该方法在马赫数6、雷诺数7\times 10^6条件下,分别对乘波体进行了最大升阻比的单目标和综 合考虑升阻比、容积及表面积的多目标优化. 计算结果表明,采用神经网络响应面技术可在 保证计算稳定性的条件下有效提高计算效率.     

乘波体前体/进气道一体化设计研究

为了研究乘波体几何外形参数和飞行参数对前体/进气道一体化设计的影响,采用理论分析和数值模拟相结合的方法,以马赫数Ma=6和攻角α=0为设计状态、进气道总压恢复系数和前体阻力系数为目标函数,对乘波体前体/进气道进行了优化设计,并在此基础上研究了攻角、马赫数、前缘半径、前体宽度对气动参数的影响。结果表明:该乘波体前体/进气道构型具有良好的攻角特性,总压恢复系数比基准构型提高17.79%,阻力系数比基准构型降低78.5%,符合高超声速飞行器高升力、低阻力的要求,且非常适合小攻角高超声速巡航飞行;为了得到较高升阻比的前体,在前缘半径R≤2mm的范围内进行流场反设计时,可以将设计马赫数的取值比预期低一些。     

基于当地流活塞理论的气动弹性计算方法研究

发展了一种高效、高精度的超音速、高超音速非定常气动力计算 方法------基于定常CFD技术的当地流活塞理论. 运用当地流活塞理论计算非定常 气动力，耦合结构运动方程，实现超音速、高超音速气动弹性的时域模拟. 运用这 种方法计算了一系列非定常气动力算例和颤振算例，并和原始活塞理论、非定 常Euler方程结果作了比较. 由于局部地使用活塞理论假设，这种方法大大地克服 了原始活塞理论对飞行马赫数、翼型厚度和飞行迎角的 限制. 与非定常Euler方程方法相比，当地流活塞理论的效率很高.     

自由活塞压缩管ALE方法数值模拟

当前国际上实现高焓气体流动的实验手段之一是自由活塞驱动类脉冲设备,包括自由活塞激波风洞和自由活塞膨胀管.采用自由活塞压缩管作为激波风洞和膨胀管的驱动段时,其驱动能力在很大程度上决定了该类设备的性能.本文采用计算流体力学中任意拉格朗日——欧拉方法(arbitrary Lagrangian Eulerian)数值模拟了压缩管内部的自由活塞运动和气体流动特征.采用移动网格技术来适应活塞运动边界,耦合求解网格运动和气体流动过程,并通过双时间步长方法进行流体运动的时间积分.为了满足几何守恒律(geometric conservation law),对移动网格的法向矢量和表面面积计算进行了修正.不同时刻的活塞位置试验测量结果及欧拉方法预测结果,以及基于简单波理论获得的运动活塞底部气体压力、活塞速度与活塞位置都与当前的ALE方法十分一致.该工作为下一步数值模拟自由活塞激波风洞和自由活塞膨胀管中包括压缩管、激波管和喷管等不同部位的耦合流动提供了基础.      

Lift-Drag Ratio of a Hypersonic Waverider of Delta Planform

The effect of the aperture angle of a V-wing representing a waverider undersurface on the hypersonic waverider lift-drag ratio is investigated for two isoperimetric conditions. The calculation model for the aerodynamic loads on angular configurations makes allowance for viscous interaction on a certain range of flight and geometry parameters on which there is no separation of the boundary layer flow within the shock layer. It is shown that for given values of the lift coefficient and specific volume, the waverider with a V-wing-shaped undersurface can have a lift-drag ratio which considerably exceeds that of a waverider of delta planform with a plane undersurface.     

Diagnostic modelling of an expansion tube operating condition

Computational simulations of an expansion tube were conducted to estimate flow parameters and verify experimental uncertainties. Two types of simulations of the complete facility were undertaken: a one-dimensional simulation, and a hybrid simulation where a one-dimensional simulation of the shock tube section was coupled with a two-dimensional simulation of the acceleration tube. Good agreement between the one-dimensional simulations and experiments were obtained in the shock tube portion of the facility. In the acceleration section, initial two-dimensional simulations did not match the experimentally measured pitot pressure and showed a discrepancy in the shock speed. Further studies examined how the accelerator gas composition affected shock speed, static pressure and pitot pressure levels in expansion tube operation. Subsequent two-dimensional simulations, using an 8% level of air contamination in helium, showed reasonable agreement with experimental data. This prediction of air contamination was later confirmed with experimental measurements of the air partial pressure before operation.      

The numerical and experimental simulation of hypervelocity flow around the HYFLEX vehicle forebody

Numerical and experimental techniques are used to model the flow and pressure distribution around the forebody of the HYFLEX hypersonic flight vehicle. We compare numerical simulation results with modified Newtonian theory and flight data to determine the accuracy of the computational fluid dynamics (CFD) technique used. The numerical simulations closely match the trends in flight data, and show that real gas effects have a small but significant influence on the nose pressure distribution. We also present pressure results from a scale-model tested in a shock tunnel, and compare them with simulation results. For the shock tunnel experiment, the model was placed such that part of the upper surface was in a region of the test flow where nonuniformities were significant, and it was shown that the numerical simulation could adequately capture these experimental flow features. The binary scaling parameter (describing the similarity in species dissociation between flight and model) was used to design the scale-model tests in the shock tunnel, and its effectiveness is discussed. We find that matching the flight Mach number in the shock tunnel experiment is not critical for reproducing flight pressure data, so long as flight velocity is matched, and binary scaling is maintained. Received 11 June 1998 / Accepted 1 September 1998     

基于局部偏转吻切方法的多级压缩乘波体设计

乘波体因优异的气动特性,被认为是突破现有"升阻比屏障"的有效途径之一,已成为高超声速飞行器气动设计的研究热点.针对常规单级压缩乘波前体压缩量不足的问题,基于局部偏转吻切方法提出一种多级压缩乘波体设计方法,实现了多道非轴对称激波的逆向乘波设计.通过引入多道非轴对称激波,可充分发挥乘波前体的预压缩效果,并为复杂外形条件下的...     

Supersonic Flow Around a V-Shaped Wing at Incidence and Yaw

A solution of the problem of the flow around a V-wing with supersonic leading edges at low angles of attack and yaw is obtained within the framework of the linear theory. Possible patterns of nonsymmetric flow around the wing are analyzed as functions of the wing geometry and the freestream velocity direction, and the ranges of angles of attack and yaw on which these patterns are realized are established. Some previously undescribed shock wave configurations are found to exist in the wing-induced conical flows.     

Numerical simulation of flapping‐wing insect hovering flight at unsteady flow

A computational fluid dynamics (CFD) analysis was conducted to study the unsteady aerodynamics of a virtual flying bumblebee during hovering flight. The integrated geometry of bumblebee was established to define the shape of a three‐dimensional virtual bumblebee model with beating its wings, accurately mimicking the three‐dimensional movements of wings during hovering flight. The kinematics data of wings documented from the measurement to the bumblebee in normal hovering flight aided by the high‐speed video. The Navier–Stokes equations are solved numerically. The solution provides the flow and pressure fields, from which the aerodynamic forces and vorticity wake structure are obtained. Insights into the unsteady aerodynamic force generation process are gained from the force and flow‐structure information. The CFD analysis has established an overall understanding of the viscous and unsteady flow around the virtual flying bumblebee and of the time course of instantaneous force production, which reveals that hovering flight is dominated by the unsteady aerodynamics of both the instantaneous dynamics and also the past history of the wing. A coherent leading‐edge vortex with axial flow and the attached wingtip vortex and trailing edge vortex were detected. The leading edge vortex, wing tip vortex and trailing edge vortex, which caused by the pressure difference between the upper and the lower surface of wings. The axial flow, which include the spanwise flow and chordwise flow, is derived from the spanwise pressure gradient and chordwise pressure gradient, will stabilize the vortex and gives it a characteristic spiral conical shape. Copyright © 2006 John Wiley & Sons, Ltd.     

基于POD和代理模型的静气动弹性分析方法

静气动弹性问题考虑弹性结构与定常气动力间的相互耦合作用,对飞行器的性能和安全具有显著的影响.在现代飞行器设计阶段,计算流体力学(CFD)/计算结构力学(CSD)直接耦合方法是精确考察静气动弹性影响的重要手段.然而,基于CFD技术的气动力仿真手段在耦合过程中计算量大且耗时长,难以满足设计阶段的需求.因此,为了兼顾计算精度与效率,文章采用本征正交分解(POD)和Kriging代理模型相结合的模型降阶方法,替代CFD求解过程并耦合有限元分析(FEA)方法,建立了高效、准确的静气动弹性分析框架.相较于传统的以模态法为主的静气动弹性分析方法,该方法能够解决更为复杂的静气动弹性问题以及提供静气动弹性变形过程中的气动分布载荷.针对典型三维跨声速HIRENASD机翼模型开展的马赫数、迎角变化的算例验证表明:由建立的静气动弹性分析方法与CFD/CSD直接耦合方法计算得到机翼翼梢处的静变形量间的相对误差在5%以内;同时该方法预测静平衡位置处的气动分布载荷的误差在5%以内,静气动弹性分析的计算效率至少提升了6倍.     

基于雨燕翅膀的仿生三角翼气动特性计算研究

针对低雷诺数微型飞行器的气动布局,设计出类似雨燕翅膀的一组具有不同前缘钝度的中等后掠(Λ=50?)仿生三角翼.为了定量对比研究三角翼后缘收缩产生的气动效应,设计了一组具有同等后掠的普通三角翼.为了深入研究仿生三角翼布局的前缘涡演化特性以及总体气动特性,采用数值模拟方法详细地探索了低雷诺数(Re=1.58×104)流动条...