火星进入器小攻角飞行的静不稳定性

针对火星探路者号进入段飞行过程,求解三维流体力学Navier-Stokes方程,采用高温真实气体模型和等效比热比完全气体模型,对气动力特性进行预测和分析,对发现的小攻角静不稳定现象开展机理分析.比较海盗号数值计算结果、文献数据和飞行试验数据,结果符合很好,验证了物理化学模型及数值方法;探路者号的升阻特性和配平特性预测数据表明,真实气体模型和LAURA数据符合很好,等效比热比模型可以较准确地计算升力和阻力系数;2°攻角飞行时,进入器沿轨道出现静不稳定,完全气体模型无法模拟,分析认为迎风区声速线在肩部的移动和背风区亚声速区的分布变化以及亚声速区泡的出现,导致激波层压力下降和肩部膨胀区的压力下降向上游的传递过程发生变化,造成表面压力分布改变,最终诱发静不稳定.     

二维叶片襟翼增升的试验研究

本文在雷诺数为0.24×106，0°－40°攻角范围内测定了带有90°，45°，135°三种不同偏转角的位于压力面后缘的襟翼的NACA632—215叶片的表面压力分布．根据压力分布求得的升力系数和阻力系数（压差）表明，襟翼能在不同程度上起到增升的效果，且当襟翼高度相同时，偏转角为90°时增升效果最佳；当偏转角为90°时，襟翼高度越大，升力系数增加越多．     

标模外形化学非平衡流场数值模拟

利用三维并行计算代码求解Navier-Stokes方程,数值模拟标模(ELECTRE)化学非平衡绕流,研究真实气体效应对标模气动热特性的影响,反应模型为Dunn和Kang的7组元7反应化学动力学模型.利用典型弹道点的飞行试验数据验证化学非平衡流计算程序的可靠性.在此基础上,研究不同壁面催化条件下攻角和高度变化对热流的影响.计算表明:真实气体效应主要发生在物面附近很薄的激波层内,并使激波脱体距离减小;完全催化壁驻点热流值高于非催化壁热流值;随着攻角增大,热流分布差异明显,而且攻角越大时,物面电子数密度越小;飞行高度越高,O₂和N₂离解程度越低,驻点热流越低.     

高升阻比模型天平校准安装角误差控制要求

文章针对高升阻比类模型气动力测量的要求,对天平校准安装角误差影响进行了初步分析,天平的单元校准中,主要有两个安装角对天平校准结果产生影响,这两个安装角分别为攻角和滚转角,研究结果表明:天平的不同分量和不同校准系数对安装角误差的要求是不尽相同的,力分量以及滚转力矩校准,安装角误差可以允许在度水平上;俯仰力矩和偏航力矩校准攻角误差只能在分量级上;法向力对轴向力干扰系数校准的攻角误差,以及法向力对侧向力干扰系数校准的滚转角误差均只能允许在十分之几分或百分之几分的水平上.单纯从校准来看,法向力对轴向力的干扰系数要求攻角误差控制,以及法向力对侧向力干扰系数要求的滚转角误差控制,均非常困难;结合风洞实验测量中高升阻比模型的载荷匹配情况,这个误差控制可以适当放宽,攻角误差可以放宽到分量级上,滚转角误差仍为十分之几分量级上,天平校准时采用常规的角度测量仪器进行测量和调整均非常困难,只有配置高精度的角度测量仪器,才能满足高升阻比模型气动力测量对天平校准角度误差控制的要求.      

基于自由涡的风轮失谐气动特性研究

采用一种改进的自由涡方法研究了风切变条件下,叶片安装角存在偏差的失谐风轮的气动特性。计算方法由模拟叶片气动力的非线性升力线法和模拟尾迹涡运动的时间精确自由尾迹法构成。以某2.5 MW风力机为例,计算了功率曲线和叶片气动力分布,并与Bladed软件和CFD计算结果进行了比较。研究了切变风条件下,不同叶片安装角偏差量,风轮推力、扭矩、偏航力矩和俯仰力矩的变化。结果表明,所采用的计算方法是准确有效的。风轮叶片安装角偏差,对风轮推力和扭矩影响较小,对风轮偏航和俯仰力矩影响较大。不同叶片安装角偏差量相反,会显著增加风轮偏航和俯仰力矩波动。     

基于自由涡的风轮失谐气动特性研究EI北大核心CSCD

采用一种改进的自由涡方法研究了风切变条件下,叶片安装角存在偏差的失谐风轮的气动特性。计算方法由模拟叶片气动力的非线性升力线法和模拟尾迹涡运动的时间精确自由尾迹法构成。以某2.5 MW风力机为例,计算了功率曲线和叶片气动力分布,并与Bladed软件和CFD计算结果进行了比较。研究了切变风条件下,不同叶片安装角偏差量,风轮推力、扭矩、偏航力矩和俯仰力矩的变化。结果表明,所采用的计算方法是准确有效的。风轮叶片安装角偏差,对风轮推力和扭矩影响较小,对风轮偏航和俯仰力矩影响较大。不同叶片安装角偏差量相反,会显著增加风轮偏航和俯仰力矩波动。     

垂直轴风力机叶片表面结冰的风洞试验

为研究垂直轴风力机叶片表面结冰的规律以及结冰对其性能的影响,对采用NACA0018翼型的风力机叶片进行了风洞结冰试验研究。在风洞试验段内安装了喷水装置,室外的寒冷空气被吸入风洞后与过冷水滴一起吹向叶片并碰撞结冰。测试了不同水滴流量和叶片攻角下的叶表结冰分布及叶片的升阻力系数变化。在一定攻角范围内,叶表结冰量随翼型迎风面积增加而增加;结冰后的阻力系数增大,升力系数减小,叶片的气动特性降低。     

双锥/双楔流动中的高温气体效应仿真模拟

在高超声速飞行和再入地球大气过程中, 气体分子的振动、电子态激发, 伴随离解、电离反应, 从而产生高温真实气体效应。不同数值方法对高温真实气体效应的模化会造成气体热物性参数的差异, 从而对流场模拟引入不确定度。以高超声速的双锥/双楔流动为例, 通过计算流体力学方法和直接模拟Monte Carlo (DSMC)方法, 研究高温真实气体模型对复杂干扰流动的预测能力。结果表明, 有别于量热完全气体, 若考虑真实物理过程的热化学非平衡过程带来气体热力学性质、输运特性的变化, 会导致激波角、边界层厚度、分离区尺寸等流动结构的改变。因此, 在研究高超声速模拟中应注意数值模型的正确应用。      

光学炮塔气动载荷的非稳态特性

采用耦合J-B模型的IDDES模型与双时间步LU-SGS方法开展炮塔非稳态气动载荷的数值仿真研究。炮塔流动会发生分离,天顶位置的分离角大于90°;当流动绕过炮塔时,形成马蹄涡、脱落涡街等非稳态流场结构,导致气动载荷也具有非稳态特性;炮塔顶点的脉动静压功率谱在1.6~40.0 kHz进入各向同性均匀湍流的惯性子区,基本满足Kolmogrov的-5/3定律;气动力以阻力为主,横向力的脉动幅值大,气动力矩则以俯仰力矩为主,滚转力矩的脉动幅值大,偏航力矩可以忽略不计;气动力和力矩的功率谱主要集中在1 kHz以下,存在多个尖峰频率,主频约为230 Hz (斯特劳哈尔数为0.15)。在ATP系统设计之初,需要考虑光学炮塔所受气动载荷的非稳态特性,并规避尖峰频率尤其是主频的谐振破坏问题。     

高超声速飞行器复杂外形转捩预测

采用k-ω-γ转捩模式对某新型飞行器外形的典型流动特征和边界层失稳特性进行了分析.研究结果表明,横流是影响飞行器大面积转捩的主要因素.随着高度增加,来流Reynolds数减小,迎风面和背风面的转捩起始位置均向下游移动.随着攻角增加,头部附近背风面的展向压力梯度增大,横流效应增强,转捩起始位置向上游移动;另一方面攻角增加导致头部激波增强,波后迎风面密度显著增大,边界层外缘Reynolds数增大,导致迎风面转捩提前发生.0°攻角下背风面中心线附近由压缩面诱导的流动分离导致转捩提前,产生"凸"字型转捩型线,5°攻角时该流动分离发生于转捩之后,"凸"字型转捩型线消失.      

Variation of the temperature factor of models on a ballistic range

A technique for variation of the temperature factor of free-flight models by varying their initial temperature is described. An experiment on a ballistic range is carried out with a free-flying supersonic blunt cone with a half-angle of 15° at a Mach number of 2.3. The flow at the cone base is studied in the transition range (from the laminar to turbulent flow) of Reynolds numbers. The base flow pattern is determined from the shadowgraphs of the flow about the models. The drag coefficient of the blunt cone at a zero angle of attack is found by processing trajectory data. It is found that the near wake geometries and the drag coefficients of the models tested at the laboratory temperature and a temperature of 120 K differ. Explanations of this effect are given.     

Effect of relative camber on the aerodynamic performance improvement of asymmetrical blunt trailing-edge modification

In this paper, the aerodynamic performance of the S series of wind turbine airfoils with different relative cambers and their modifications is numerically studied to facilitate a greater understanding of the effects of relative camber on the aerodynamic performance improvement of asymmetrical blunt trailing-edge modification. The mathematical expression of the blunt trailing-edge modification profile is established using the cubic spline function, and S812, S816 and S830 airfoils are modified to be asymmetrical blunt trailing-edge airfoils with different thicknesses. The prediction capabilities of two turbulence models, the k-ω SST model and the S-A model, are assessed. It is observed that the k-ω SST model predicts the lift and drag coefficients of S812 airfoil more accurately through comparison with experimental data. The best trailing-edge thickness and thickness distribution ratio are obtained by comparing the aerodynamic performance of the modifications with different trailing-edge thicknesses and distribution ratios. It is, furthermore, investigated that the aerodynamic performance of original airfoils and their modifications with the best thickness of 2% c and distribution ratio being 0:4 so as to analyze the increments of lift and drag coefficients and lift–drag ratio. Results indicate that with the increase of relative camber, there are relatively small differences in the lift coefficient increments of airfoils whose relative cambers are less than 1.81%, and the lift coefficient increment of airfoil with the relative camber more than 1.81% obviously decreases for the angle of attack less than 6.3°. The drag coefficient increment of S830 airfoil is higher than that of S816 airfoil, and those of these two airfoils mainly decrease with the angle of attack. The average lift–drag ratio increment of S816 airfoil with the relative camber of 1.81% at different angles of attack ranging from 0.1° to 20.2° is the largest, closely followed by S812 airfoil. The lift–drag ratio increment of S830 airfoil is negative as the angle of attack exceeds 0.1°. Thus, the airfoil with medium camber is more suited to the asymmetrical blunt trailing-edge modification.     

再入飞行器失稳的分叉理论分析与数值仿真验证

飞行器再入大气层时的姿态稳定性事关飞行安全, 是气动设计的关键问题之一.文章采用非线性自治动力系统分叉理论, 耦合求解非定常Navier-Stokes方程和俯仰运动方程, 研究了钝体和细长体两类航天飞行器再入过程单自由度俯仰运动失稳问题.研究表明, 航天飞行器再入时, 如果仅有1个配平攻角, 随Mach数降低, 其配平攻角处的俯仰姿态失稳一般对应于Hopf分叉, 并存在亚临界Hopf分叉和超临界Hopf分叉两种失稳形态; 如果再入时随着Mach数的降低, 其配平攻角由1个演化至多个(一般为3个), 其配平攻角处的俯仰姿态失稳形态将更为复杂, 可能发生鞍结点分叉形态的刚性失稳行为;随Mach数的进一步降低, 其俯仰运动还可能进一步发生Hopf分叉和同宿分叉.      

Numerical simulation of the thermal non-equilibrium flow-field characteristics of a hypersonic Apollo-like vehicle

In order to investigate the relationship between the flow-field parameters outside the vehicle and the altitude, this paper takes the Atmospheric Reentry Demonstrator (ARD) with an angle of attack of -20° as the research object and adopts a two-temperature model coupled with the shear-stress transport k-ω turbulence model to focus on the variation of flow-field parameters including flow-field pressure, Mach number and temperature with the reentry altitude. It is found that the flow-field high-pressure region and low-Mach region both appear in the shock layer near the head of the ARD, while the maximum pressure of the surface appears on the windward side of the ARD's head with a toroidal distribution, and the numerical magnitude is inversely proportional to the radius of the torus. With fluid through the shoulder of the ARD flow expansion plays a dominant role, the airflow velocity increases, the Mach number of the windward side of the rear cone increases and the flow-field pressure and surface pressure rapidly decrease. When the fluid passes through the shock layer, the translational-rotation temperature will increase before the vibration-electron temperature, there is a thermal non-equilibrium effect and the two temperatures will rapidly decrease again when approaching the surface of the ARD due to the existence of temperature gradient. At the same time, both the windward side of the shoulder and the back cover of the ARD suffer from a large thermal load and require thermal protection.     

Investigation of hypersonic flows through a cavity with sweepback angle in near space using the DSMC method

Near space has been paid more and more attentionin recent years due to its military application value.However,flow characteristics of some fundamental configurations(e.g.,the cavity) in near space have rarely been investigated due to rarefied gas effects,which make the numerical simulation methods based on continuous flow hypothesis lose validity.In this work,the direct simulation Monte Carlo(DSMC),one of the most successful particle simulation methods in treating rarefied gas dynamics,is employed to explore flow characteristics of a hypersonic cavity with sweepback angle in near space by considering a variety of cases,such as the cavity at a wide range of altitudes 20-60 km,the cavity at freestream Mach numbers of 6-20,and the cavity with a sweepback angle of 30°-90°.By analyzing the simulation results,flow characteristics are obtained and meanwhile some interesting phenomena are also found.The primary recirculation region,which occupies the most area of the cavity,causes pressure and temperature stratification due to rotational motion of fluid inside it,whereas the pressure and temperature in the secondary recirculation region,which is a small vortex and locates at the lower left corner of the cavity,change slightly due to low-speed movement of fluid inside it.With the increase of altitude,both the primary and secondary recirculation regions contract greatly and it causes them to separate.A notable finding is that rotation direction of the secondary recirculation region would be reversed at a higher altitude.The overall effect of increasing the Mach number is that the velocity,pressure,and temperature within the cavity increase uniformly.The maximum pressure nearby the trailing edge of the cavity decreases rapidly as the sweepback angle increases,whereas the influence of sweepback angle on velocity distribution and maximum temperature within the cavity is slight.     

Investigation of magnetohydrodymamic interaction in a supersonic air flow at <Emphasis Type="Italic">M</Emphasis> = 8

The influence of magnetohydrodynamic interaction localized before a model on the position of a shock wave attached to a wedge is experimentally and numerically investigated. The investigation is carried out in an air flow with a Mach number of 8. It is shown that, for a hydromagnetic interaction parameter on the order of 0.1, the slope angle of the shock wave can be increased by 10°. Experiments are conducted for the case when the flow is ionized by an electron beam or by a pulsed electric discharge. Good agreement between experimental and numerical results is obtained for both ways of ionization if the Joule heating of the gas is insignificant. The conclusion is drawn that the way of providing a nonequilibrium conductivity of the flow has a minor effect on the position of the oblique shock wave near the wedge with the hydromagnetic interaction parameter being the same.     

钝锥再入飞行器黏性干扰效应建模

针对典型的钝锥外形, 采用统一气体动理学格式(UGKS)模拟了高度70~110 km下不同Mach数和攻角的流场, 进行了流场特性的分析, 并基于黏性干扰的理论成果, 将气动力特性与第3黏性干扰参数、攻角和Mach数等参数进行关联, 建立了气动力系数的黏性干扰模型, 给出了模型预测结果的相关性分析和准确性评估。经初步测试, 该模型预测结果与UGKS直接模拟结果具有良好的一致性, 对工程应用快速获取高空气动特性具有重要意义。      

Flow characteristics of hypersonic inlets with different cowl-lip blunting methods

Under hypersonic flight conditions,the sharp cowl-lip leading edges have to be blunted because of the severe aerodynamic heating.This paper proposes four cowl-lip blunting methods and studies the corresponding flow characteristics and performances of the generic hypersonic inlets by numerical simulation under the design conditions of a flight Mach number of 6 and an altitude of 26 km.The results show that the local shock interference patterns in the vicinity of the blunted cowl-lips have a substantial influence on the flow characteristics of the hypersonic inlets even though the blunting radius is very small,which contribute to a pronounced degradation of the inlet performance.The Equal Length blunting Manner(ELM)is the most optimal in that a nearly even reflection of the ramp shock produces an approximately straight and weak cowl reflection shock.The minimal total pressure loss,the lowest cowl drag,maximum mass-capture and the minimal aeroheating are achieved for the hypersonic inlet.For the other blunting manners,the ramp shock cannot reflect evenly and produces more curved cowl reflection shock.The Type V shock interference pattern occurs for the Cross Section Cutting blunting Manner(CSCM)and the strongest cowl reflection shock gives rise to the largest flow loss and drag.The cowl-lip blunted by the other two blunting manners is subjected to the shock interference pattern that transits with an increase in the blunting radius.Accordingly,the peak heat flux does not fall monotonously with the blunting radius increasing.Moreover,the cowl-lip surface suffers from severe aerothermal load when the shear layer or the supersonic jet impinges on the wall.     

基于等离子体气动激励的斜劈诱导激波控制

基于弧光等离子体气动激励,采用不同的放电通道间距、放电通道数、放电直流输入电压、斜劈劈角、有无磁场作用等激励条件,实验研究了在超音速来流条件下(马赫数为2.2)斜激波位置、角度、强度的变化规律。结果表明:施加等离子体气动激励后,激波的起始位置平均前移1～8 mm,激波角平均减小4%～8%,激波强度平均减弱8%～26%。这主要是由于等离子体气动激励产生高温高压的表面等离子体层,使边界层分离点逆气流前移,改变了原有激波系结构,使原有的激波位置前移,激波角减小;同时由于局部的高温导致当地音速增大,使得当地马赫数减小,上述两个原因均可导致激波强度减弱。