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

针对火星探路者号进入段飞行过程,求解三维流体力学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₂离解程度越低,驻点热流越低.     

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

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

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

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

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

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

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

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

JF12长实验时间激波风洞10°尖锥气动力实验研究

介绍了在中国科学院力学研究所JF12长实验时间激波风洞上开展的10°尖锥标模的天平测力实验研究结果.JF12激波风洞的实验时间为100~130 ms,名义Mach数为7.0,喷管出口直径为2.5 m,总焓为2.5 MJ/kg,复现了35 km高空的飞行条件.采用六分量应变天平,攻角分别为-5,0,5,10和14°,模型长度为1.5 m,质量为50 kg.实验结果表明,在100~130 ms的实验时间里,应变天平的输出信号含有3~4个完整周期,可以通过对天平的输出信号进行平均直接获得气动力/矩测量结果,而不再需要进行加速度补偿,且气动力系数重复测量的不确定度小于2%.JF12激波风洞气动力系数的测量结果与传统高超声速风洞的结果符合得较好,表明在2.5 MJ/kg的总焓下,真实气体效应对该模型气动力特性的影响不明显.      

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

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

Nonlinear Modeling Study of Aerodynamic Characteristics of an X38-like Vehicle at Strong Viscous Interaction Regions

Strong viscous interaction and multiple flow regimes exist when vehicles fly at high altitude and high Mach number conditions. The Navier–Stokes(NS) solver is no longer applicable in the above situation. Instead, the direct simulation Monte Carlo (DSMC) method or Boltzmann model equation solvers are usually needed. However, they are computationally more expensive than the NS solver. Therefore, it is of great engineering value to establish the aerodynamic prediction model of vehicles at high altitude and high Mach number conditions. In this paper, the hypersonic aerodynamic characteristics of an X38-like vehicle in typical conditions from 70 km to 110 km are simulated using the unified gas kinetic scheme (UGKS), which is applicable for all flow regimes. The contributions of pressure and viscous stress on the force coefficients are analyzed. The viscous interaction parameters, Mach number, and angle of attack are used as independent variables, and the difference between the force coefficients calculated by UGKS and the Euler solver is used as a dependent variable to establish a nonlinear viscous interaction model between them in the range of 70–110 km. The evaluation of the model is completed using the correlation coefficient and the relative orthogonal distance. The conventional viscous interaction effect and rarefied effect are both taken into account in the model. The model can be used to quickly obtain the hypersonic aerodynamic characteristics of X38-like vehicle in a wide range, which is meaningful for engineering design.     

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.     

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.