Instability analysis and drag coefficient prediction on a swept RAE2822 wing with constant lift coefficient

Swept wing is widely used in civil aircraft,whose airfoil is chosen,designed and optimized to increase the cruise speed and decrease the drag coefficient.The parameters of swept wing,such as sweep angle and angle of attack,are determined according to the cruise lift coefficient requirement,and the drag coefficient is expected to be predicted accurately,which involves the instability characteristics and transition position of the flow.The pressure coefficient of the RAE2822 wing with given constant lift coefficient is obtained by solving the three-dimensional Navier-Stokes equation numerically,and then the mean flow is calculated by solving the boundary layer(BL) equation with spectral method.The cross-flow instability characteristic of boundary layer of swept wing in the windward and leeward is analyzed by linear stability theory(LST),and the transition position is predicted by eNmethod.The drag coefficient is numerically predicted by introducing a laminar/turbulent indicator.A simple approach to calculate the lift coefficient of swept wing is proposed.It is found that there is a quantitative relationship between the angle of attack and sweep angle when the lift coefficient keeps constant;when the angle of attack is small,the flow on the leeward of the wing is stable.when the angle of attack is larger than 3°,the flow becomes unstable quickly;with the increase of sweep angle or angle of attack the disturbance on the windward becomes more unstable,leading to the moving forward of the transition position to the leading edge of the wing;the drag coefficient has two significant jumping growth due to the successive occurrence of transition in the windward and the leeward;the optimal range of sweep angle for civil aircraft is suggested.     

脉冲爆震发动机最大净推力和最佳工作频率

为研究气动阀式脉冲爆震发动机净推力和工作频率之间的变化规律,建立了PDE净推力的近似计算模型。通过计算和试验的方法,研究了PDE能够达到的最大净推力及相应的最佳工作频率,以及爆震管长度L和PDE的阻力经验系数K1,K2等因素对PDE的最大净推力、最佳工作频率的影响。结果表明:L,K1,K2一定,PDE存在一个最佳工作频率fopt,在该频率下,PDE净推力最大。随着L增加,fopt下降;随着K1,K2增加,fopt和最大净推力都下降。研究结果对于PDE优化设计,提高PDE的工作性能具有重要的参考价值。     

测定飞行推力和极曲线新方法

根据飞机在纵向运动时沿风轴方向的运动方程,重新定义飞机的推力和阻力,提出了一种新的测定飞行推力和阻力的方法,并将此法用于某型飞机的风洞与飞行试验极曲线相关性研究,取得了好的结果.     

Study of riblet drag reduction for an infinite span wing with different sweep angles

This paper studies the riblet drag reduction effect for an infinite swept wing under a low Reynolds number using a large-eddy simulation. The results show that the drag reduction ratio is not linear under different sweep angles. The maximum drag reduction ratio in this study is 9.5% for a wing with a 45° sweep angle. The local surface streamline angle and turbulence quantities are calculated to analyze the drag reduction mechanism. The results demonstrate that the riblets considerably suppress the Reynolds stresses above the wing upper surface, while the turbulence kinetic energy in the near wake is increased. A possible relaminarization phenomenon is observed at the middle part of the wing. Quasi-two-dimensional flow structures are observed near the wall, and a peak frequency is considered as the dominant frequency of the region.     

Drag reduction research in supersonic flow with opposing jet

A CFD study on drag reduction in supersonic flow with opposing jet has been conducted. Flowfield parameters, reattachment point position and surface pressure distributions are obtained and validated with experiments. From the analysis on the physical mechanism of drag reduction, it shows the phenomenon that, when the opposing jet blows, the high pressure region is located between the bow shock wave and the Mach disk, which makes the nose region much lower pressure. As the pressure ratio increases, the high pressure region is gradually pushed away from the surface. Larger the total pressure ratio is, the lower of the drag coefficient is. To study the effect of the intensity of opposing jet more reasonably, a new parameter R_PA has been introduced by combining the flux and the total pressure ratio. The study shows that the same shock wave position and drag coefficient can be obtained with the same R_PA with different fluxes and the total pressures, which means the new parameter could stand for the intensity of opposing jet and could be used to analyze the influence of opposing jet on flow field and aerodynamic force.     

从动态飞行获得极曲线及相关性研究

本文研究飞机极曲线试飞及其与风洞的相关性问题.主要通过动态飞行获得飞机极曲线(C_y≤0.3)的飞行数据.飞行条件:高度为3～12 km;M数为0.8,0.9,1.2,1.8.结果表明:动态飞行技术是可行的,飞行极曲线与风洞极曲线的相关性较好.     

Detached-eddy simulation of supersonic flow past a spike-tipped blunt nose

Space vehicle in atmosphere travels mostly at supersonic speed and generates a very strong bow shockwave around its blunt nose. Oblique shock and conical separated flow zone generated by a forward disk-tip spike significantly reduce the drag by reducing the high pressure area on the blunt nose. This study employs improved delayed detached eddy simulation to investigate the characteristic flow structures around a spike-tipped blunt nose at Mach number of 3 and Reynolds number(based on the blunt-body diameter) of 2.72x10~8. The calculated time-averaged quantities agree well with experimental data. Characteristic frequencies in different flow regions are extracted using fast Fourier transform. It is found that two distinct instability modes exist: oscillation mode and pulsation mode. The former is related to the foreshock/turbulence interaction with nondimensional frequency at around 0.004. The latter corresponds to the interaction between turbulence and shock structures around the blunt nose, with a typical coherent structure shedding frequency at 0.092.     

Drag and energy accommodation coefficients during sunspot maximum

Conditions appropriate to gas-surface interactions on satellite surfaces in orbit have not been successfully duplicated in the laboratory. However, measurements by pressure gauges and mass spectrometers in orbit have revealed enough of the basic physical chemistry that realistic theoretical models of the gas-surface interaction can now be used to calculate physical drag coefficients. The dependence of these drag coefficients on conditions in space can be inferred by comparing the physical drag coefficient of a satellite with a drag coefficient fitted to its observed orbital decay. This study takes advantage of recent data on spheres and attitude stabilized satellites to compare physical drag coefficients with the histories of the orbital decay of several satellites during the recent sunspot maximum. The orbital decay was obtained by fitting, in a least squares sense, the semi-major axis decay inferred from the historical two-line elements acquired by the US Space Surveillance Network. All the principal orbital perturbations were included, namely geopotential harmonics up to the 16th degree and order, third body attraction of the Moon and the Sun, direct solar radiation pressure (with eclipses), and aerodynamic drag, using the Jacchia-Bowman 2006 (JB2006) model to describe the atmospheric density. After adjusting for density model bias, a comparison of the fitted drag coefficient with the physical drag coefficient has yielded values for the energy accommodation coefficient as well as for the physical drag coefficient as a function of altitude during solar maximum conditions. The results are consistent with the altitude and solar cycle variation of atomic oxygen, which is known to be adsorbed on satellite surfaces, affecting both the energy accommodation and angular distribution of the reemitted molecules.     

Command generator tracker based direct model reference adaptive tracking guidance for Mars atmospheric entry

In order to accurately deliver an entry vehicle through the Martian atmosphere to the prescribed parachute deployment point, active Mars entry guidance is essential. This paper addresses the issue of Mars atmospheric entry guidance using the command generator tracker (CGT) based direct model reference adaptive control to reduce the adverse effect of the bounded uncertainties on atmospheric density and aerodynamic coefficients. Firstly, the nominal drag acceleration profile meeting a variety of constraints is planned off-line in the longitudinal plane as the reference model to track. Then, the CGT based direct model reference adaptive controller and the feed-forward compensator are designed to robustly track the aforementioned reference drag acceleration profile and to effectively reduce the downrange error. Afterwards, the heading alignment logic is adopted in the lateral plane to reduce the crossrange error. Finally, the validity of the guidance algorithm proposed in this paper is confirmed by Monte Carlo simulation analysis.     

Thermosphere density and wind measurements in the equatorial region using a constellation of drag balance nanospacecraft

A mission for in situ thermosphere density and winds measurement is described, based on nanospacecraft equipped with a drag balance instrument (DBI) and a GPS receiver. The mission is based on nanosatellite clusters deployed in three orbital planes. In this study, clusters of 10 nanospacecraft are considered, leading to a mission based on a total of 30 nanospacecraft. The geometry analyzed is a symmetrical one, including an equatorial orbit and two orbits with the same inclination and opposing ascending nodes. The main idea is that, by combining the accurate information on the satellite inertial position and velocity provided by the GPS receiver and the drag acceleration intensity provided by the DBI, due to the orbits’ geometrical configuration, both atmospheric drag and wind can be resolved in a region close to the orbit nodes. Exploiting the Earth oblateness effect, a complete scan of the equatorial regions can be accomplished in the short mission lifetime typical of very low Earth orbit satellites, even in high solar activity peaks, when the expected nanospacecraft lifetime is about 40 days.