基于粒子轨迹的结冰风洞收缩段优化设计数值模拟

为了研究收缩段型面和液滴轨迹特性对结冰风洞流场品质与结冰试验效果的重要影响,获得收缩段的优化设计参数,采用数值模拟方法,使用二维简化模型对3m×2m结冰风洞收缩段的流场特性进行了研究,对几种典型收缩曲线的流场进行了计算,对比了收缩曲线对沿轴向压力分布、出口截面速度分布、不同粒径云雾粒子的极限释放距离及试验段入口云雾粒子分布的影响,结果表明:三次曲线与x_m为0.4的五次曲线不能满足流场合格指标,其他曲线流场特性差别不大,维氏曲线具有最小极限释放距离,在收缩段出口,三次曲线和五次曲线粒子包络面积最大,综合比较,x_m为0.5的五次曲线能有效兼顾流场均匀性、极限释放距离和收缩段出口粒子包络面积比.      

Effects of wing flexibility on aerodynamic performance of an aircraft model

The low-speed wind tunnel experiment is carried out on a simplified aircraft model to explore the influence of wing flexibility on the aircraft aerodynamic performance. The investigation involves the measurements of force, membrane deformation and velocity field at Reynolds number of 5.4 × 10⁴–1.1 × 10⁵. In the lift curves, two peaks are observed. The first peak, corresponding to the stall, is sensitive to the wing flexibility much more than the second peak, which nearly keeps constant. For the optimal case, in comparison with the rigid wing model, the delayed stall of nearly 5° is achieved, and the relative lift increment is about 90%. It is revealed that the lift enhanced region corresponds to the larger deformation and stronger vibration, which leads to stronger flow mixing near the flexible wing surface. Thereby, the leading-edge separation is suppressed, and the aerodynamic performance is improved significantly. Furthermore, the effects of sweep angle and Reynolds number on the aerodynamic characteristics of flexible wing are also presented.     

栅格翼电弧风洞试验的气动参数模拟方法研究

为验证栅格翼防热方案的有效性,针对栅格翼单片级试验件厚度过低无法直接测量试验件冷壁热流密度的难题,提出了一种直接模拟来流气动参数的栅格翼热防护试验设计方法,能够准确考核防热方案,对类似结构的热防护试验具有指导意义。     

飞翼布局飞行器等离子体激励滚转操控试验

飞翼布局飞行器采用多个气动舵面共同作用来控制飞行,常规气动舵面的结构复杂,在大迎角时由于流动分离,舵面操纵效率显著降低。等离子体激励器具有结构简单、重量轻和响应快等优势,常被用在流动控制上。本文利用激励器抑制单侧翼面流动分离产生不对称的气动力,对飞翼布局飞行器滚转通道的控制进行了试验研究,得出了激励器在飞行器上的最优布置位置和最佳控制参数,并和常规副翼舵面滚转操控效果进行了对比。结果表明:布置于内翼、中翼前缘的等离子体激励器能够获得最佳的滚转控制效果;激励器调制频率对飞行器滚转控制效果的影响较大,而激励电压对滚转控制效果的影响较小;与常规副翼相比,等离子体激励器在大迎角时对滚转通道的操控效果优于副翼。     

侧风对舰载直升机悬停性能的影响

针对舰载直升机小速度、大侧滑的飞行状态,采用了旋翼非均匀入流模型,导出了直升机在风场中的运动方程。以某型机为例,对不同风速条件下直升机悬停时的平衡特性进行了计算分析。计算结果表明:左侧风对该机悬停性能影响最大,而尾桨操纵裕度要求左侧风速不能超过30m/s,与同类直升机的飞行数据吻合。分析模型为制订舰载直升机的风险图提供了理论依据。     

Charging of dust grains by anisotropic solar wind multi-component plasma

In this paper we study the charging process of small grain particles by anisotropic multi-component solar wind plasmas (electrons, protons and heavy ions), versus two-component (electron/proton) plasmas. We are focusing attention on the important characteristics of the charging process, namely the charging time, floating potential and current content as functions of plasma parameters such as He⁺⁺/H⁺ (α/p) number density and T_α/T_p temperature ratios of alpha particles to protons, as well as plasma streaming velocity v₀. Measured statistical properties of solar wind plasma parameters at 1 AU show considerable variations in α/p-temperature ratios from 1 to 10, in α/p-number density ratio from 0.01 to 0.35, as well as in values of streaming velocity v₀ from 200 km/s to 1000 km/s and more. Periods of these variations could last for several days each, leading to significant variability in the charging process, according to newly derived general analytical expressions. Numerical calculations performed for protons/alphas plasmas showed large disparity in the charging characteristics. For example, in anisotropic plasma, grain charging time varies up to 90% depending on α/p-particles temperature and number density ratios, whereas changes in floating potential are up to 40%. In contrast, in isotropic plasma, charging characteristic for grains do not change very much for the same plasma parameters variations, with charging time varying about 12% and floating potential only varying about 4%. It is also shown that in highly anisotropic plasma, with all ballistic electrons and ions, dust grains could not hold their charges, and characteristic discharged time is calculated. We note that the analysis is equally applicable to any sized body immersed in solar wind plasma.     

The heliospheric modulation of cosmic rays: Effects of a latitude dependent solar wind termination shock

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock and that it has begun its merger with the local interstellar medium. The compression ratio of this shock affects galactic cosmic rays when they enter the heliosphere. Hydrodynamic (HD) models show that the compression ratio can have a significant latitude dependence; with the largest value in the nose direction of the heliosphere, becoming significantly less towards the polar regions. The modulation effects of such large latitude dependence are studied, using a well-established numerical drift and shock modulation model. We focus on computing the modulated spectra for galactic protons with emphasis on the radial and polar gradients in the equatorial plane and at a polar angle of θ = 55°, corresponding to the heliolatitude of Voyager 1. Two sets of solutions are computed and compared each time; with and without a latitude dependence for the compression ratio. All computations are done for the two magnetic field polarity cycles assuming solar minimum conditions. Including the termination shock in the model allows the study of the re-acceleration of galactic protons in the outer heliosphere. We find that for the A < 0 polarity cycle the intensity between ∼200 MeV and ∼1 GeV in the vicinity of the shock in the heliospheric equatorial plane may exceed the local interstellar value specified at the heliopause. Unfortunately, at θ = 55°, the effect is reduced. This seems not possible during an A > 0 cycle because significant modulation is then predicted between the heliopause and the termination shock, depending on how strong global gradient and curvature drifts are in the heliosheath. The overall effect of the shock on galactic protons in the equatorial plane is to reduce the total modulation as a function of radial distance with respect to the interstellar spectrum. Making the compression ratio latitude dependent enhances these effects at energies E < 200 MeV in the equatorial plane. At larger heliolatitudes these effects are even more significant. The differences in the modulation between the two drift cycles are compelling when the compression ratio is made latitude dependent but at Earth this effect is insignificant. A general result is that the computed radial gradient changes for galactic protons at and close to the TS and that these changes are polarity dependent. In line with previous work, large polarity dependent effects are predicted for the inner heliosphere and also close to the shock’s position in the equatorial plane. In contrast, at θ = 55°, the largest polarity effect occurs in the middle heliosphere (50 AU), enhanced by the latitude dependence of the compression ratio. At this latitude, the amount of proton modulation between the heliopause and the termination shock is much reduced. If galactic cosmic rays were to experience some diffusive shock acceleration over the 100–1000 MeV range at the shock, the radial gradient should change its sign in the vicinity of the shock, how large, depends on the compression ratio and the amount of drifts taking place in the outer heliosphere. The effective polar gradient shows a strong polarity dependence at Earth but this dissipates at θ = 55°, especially with increasing radial distance. This tendency is enhanced by making the compression ratio latitude dependent.     

Kinetic theory of sheath formation in solar wind plasma

We present a general self-consistent kinetic theory for plasma sheath formation in solar wind plasma. The theory could be applied to anisotropic, as well as to isotropic collisionless plasma without resorting to any simplifications, limitations, or assumptions, such as the necessary existence of a ‘pre-sheath’ region of ions acceleration to ensure the Bohm criterion. The kinetic framework is first applied to sheath formation around an arbitrary oriented planar absorbing surface, charged by solar wind anisotropic plasma, under the condition of negligible photoelectric effect. We then make use of our kinetic approach for the plane geometry in isotropic collisionless plasma, as a particular case of a planar electrode orientation parallel to plasma streaming velocity, also analyzing the sheath structure around spherical and cylindrical absorbing electrodes submerged in isotropic collisionless plasma. Obtained results demonstrate principal differences in spatial charge distributions in sheath regions between spherical or cylindrical electrodes of large size and an unbound planar surface submerged in isotropic plasma. In the case of a planar electrode, we directly compare results obtained in our kinetic and hydrodynamic theories and conventional hydrodynamic theory of plasma sheath formation. The outcome from the present study have direct implications to the analysis of plasma sheath structure and associated distribution in space of charged dust grains, which is relevant to the moon exploration near the optical terminator region or in shadowed craters in the moon.     

低速风洞绳牵引并联支撑系统的模型气动试验研究(英文)

根据绳牵引并联机构的优点,建造了一种新型的用于低速风洞试验的飞机模型绳牵引并联支撑系统,研究该系统中模型的空气动力参数的测量与计算方法。文中对机构进行了静力学分析,建立了描述实验模型气动载荷的数学模型,提出了通过测量绳系拉力求得模型空气动力参数的解算方法;设计并构建了绳系拉力测量及数据采集系统;将该系统置于开口式回流低速风洞中进行了吹风试验,采集了模型在不同姿态和不同风速下的各牵引绳的拉力数据,并对数据进行了处理分析,通过解算得到了不同吹风条件下模型的气动载荷参数曲线。研究结果表明,绳牵引并联机构用于低速风洞试验的支撑系统是可行的。     

内埋式模型操纵面自动变角度系统研制

针对风洞试验模型操纵面角度的改变,开发了内埋式自动变角机构及控制系统,对控制角度进行了地面校准,并采用三角翼模型进行了风洞常规测力实验、重复性实验和操纵面效率试验。结果表明,采用自动变角机构改变操纵面角度的试验精准度均较高,完全满足国军标的要求,说明本操纵面自动变角系统的研制是成功的,在风洞试验中可以代替人工更换角度片方式,大大提高实验效率。