Regularities between kinematic and aerodynamic characteristics of flexible membrane wing

The kinematic characteristics of flexible membrane wing have vital influences on its aerodynamic characteristics. To deeply explore the regularities between them, time-resolved aerodynamic forces and deformations at different aeroelastic parameters and angles of attack(α) were measured synchronously by wind tunnel experiments. The membrane motion can be mainly divided into two states at α > 0° with various lift-enhancement regularities: Deformed-Steady State(DSS)at pre-stall, and Dynamic Bala...     

Aeroelastic simulation of the first 1.5-stage aeroengine fan at rotating stall

A massive parallel aeroelastic simulation platform has been built to investigate the first 1.5-stage fan of an aeroengine at rotating stall. The Computational Fluid Dynamics (CFD) solver and Computational Structural Dynamics (CSD) solver are coupled directly by non-matching mesh interfaces. The unsteady rotor/stator interaction is solved by the Sliding Mesh Interface method. The original rotor blades are shrouded by the midspan shrouds. An unshrouded fan is also created to investigate the effects of the midspan shrouds. Both the shrouded fan and unshrouded fan have stable aeroelasticity at the designed state. At rotating stall, the stalled region rotates at 30% of the rotor speed on the absolute reference frame. The energy spectrum of the rotating stalled flow is measured quantitatively. It shows that the first two order excitations are much stronger than the higher order excitations. In the flow of rotating stall, the fifth backward travelling wave mode of shrouded fan is resonated by the fifth excitation of the rotational stalled flow because the rotational speed of the stalled region coincides with the modal rotational speed, while for the unshrouded fan, the first bending mode is resonated by the second excitation of the rotational stalled flow, forming two waves in the circumference of the annulus blades. At rotating stall, the vibration of the shrouded blades is still under control but the vibration of the unshrouded blades is diverged and out of control. A novel tool, i.e., resonance map, is proposed to predict the resonance. It provides a perspective to explain the effects of midspan shrouds at a theoretical level, and it would also be helpful in the structural design of blades.     

Consequences of the application of the streamer fluid model to the study of the sprite inception mechanism

We present the results of a streamer-fluid model used to investigate the electrodynamical coupling between the troposphere and upper atmosphere due to the penetration of lightning electric fields into the mesosphere and the lower ionosphere, generating sprites. The model solves the continuity equation for electrons and ions coupled to Poisson equation. The dominant physical response of the atmosphere is the formation of a screening-ionization wave. The wave shields the atmosphere above it from the action of the lightning field and, together with the conductivity reduction below it due to attachment, the wave amplifies the total field below it, allowing for the penetration of intense electric fields in the mesosphere as it propagates downwards into regions of higher density that compress the wave. This is the key physical mechanism for sprite inception. We evaluated the effects of the thundercloud charge geometry, lightning current waveshape, atmospheric conductivity, via different electron density profiles, and the effect of ionization, attachment and electron mobility coefficients in the electrical breakdown process, related to halo production, and sprite streamer initiation. The results showed that electrons with higher mobility are more efficient in shielding the lightning electric field before breakdown, causing delay, and they contribute to the formation of the streamer seed after breakdown, anticipating the sprite streamer inception. Similarly, a higher effective ionization rate, produced by modifications in the attachment and ionization coefficients, anticipates sprite inception. The simulations with 6 different electron density profiles, and therefore conductivities, spanning 4 orders of magnitude, showed that the altitude of breakdown and sprite initiation, as well as their time delays from the lightning discharge are directly related to atmospheric conductivity: higher conductivities produce halo and sprite inception at lower altitudes with longer delays and may hinder sprite formation. We document that variations of 30 times in the lightning current leads to sprite initiation altitudes in the range 66.0–73.5 km, with delays between 1.550 and 34.500 ms, while variations of 4 orders of magnitude in the conductivity profile lead to initiation altitudes 61.0–70.6 km, with delays in the range 3.825–9.825 ms. Consequently, we suggest that lightning characteristics dominate over atmospheric parameters in determining sprites’ initiation altitude and delay. The simulation of a −CG, with a constant current of 30 kA, did not produce a sprite seed, confirming an asymmetry in the response of the atmosphere to positive and negative lightning. This is due to the free electron drift direction that is away from the screening ionization wave, preventing the formation of the streamer seed for the great majority of −CGs. The same does not apply to halos, which depend on the occurrence of breakdown and can be produced by discharges of both polarities.     

Effect of pitch angle on initial stage of a transport airplane ditching

Airworthiness regulations require that the transport airplane should be proved to ensure the survivability of the ditching for the passengers. The planned ditching of a transport airplane on the calm water is numerically simulated. The effect of pitch angle on the impact characteristics is especially investigated by a subscaled model. The Reynolds-averaged Navier-Stokes (RANS) equations of unsteady compressible flow are solved and the realizable j-e equations are employed to model the turbulence. The transformation of the air-water interface is tracked by volume of fluid (VOF) model. The motion of the rigid body is modeled by dynamic mesh method. The initial ditching stage of the transport airplane is analyzed in detail. The numerical results show that as the pitching angle increases, the maximal normal force decreases and the pitching motion becomes much gentler. The aft fuselage would be sucked down by the water and lead to pitching up, whereas the low horizontal tail prevents this trend. Consequently, the transport aircraft with low horizontal tail should ditch on the water at an angle between 10 and 12 as a recommendation.     

表面张力贮箱匹配特性研究

本文研究两个以上表面张力贮箱在并联使用中存在的问题，发现加注和挤出过程均存在不可忽视的偏差；将贮箱液阻特性作为产生偏差的根本原因，进行了较深入的探讨；并给出大幅度降低这些偏差的有效方法。     

Numerical simulations of high enthalpy flows around entry bodies

Ablation flows around entry bodies are at highly nonequilibrium states. This paper presents comprehensive computational fluid dynamics simulations of such hypersonic flow examples.The computational scheme adopted in this study is based on the Navier–Stokes equations, and it is capable of simulating multiple-dimensional, non-equilibrium, chemically reacting gas flows with multiple species. Finite rate chemical reactions, multiple temperature relaxation processes, and ionizations phenomena with electrons are modeled. Simulation results of several hypersonic gas flows over axisymmetric bodies are presented and compared with results in the literature. It confirms that some past treatment of adopting less species for hypersonic flows is acceptable, and the differences from more species and more chemical reactions are not significant.     

叶轮机非定常气动设计的缘线匹配(Ⅰ)理论与实施方法

针对缘线匹配展开理论与实施方法探讨，分析了缘线匹配潜在效益，并给出缘线匹配的两种实施方法：基于非定常流动束缚于拟流回转面假设的基元流动展向积分方法和全优化方法。分析指出：缘线匹配具有通过沿展高调协叶轮机相邻叶排非定常流动相位关系实现气动、气弹、气动噪声和热传导性能全方位、多目标优化潜力，为叶轮机非定常设计提供了新的设计自由度。     

重力对单气泡生长特性的影响

采用VOF (Volume of Fluid)多相流模型, 通过用户自定义函数UDF (User Defined Function)实现相变过程中质量和能量的输运, 对微重力条件下尺寸为10mm × 10mm × 25mm的矩形通道的池沸腾现象进行数值模拟, 得到了微重力及常重力作用下单个气泡生长特性的差异. 模拟结果表明, 微重力条件下气泡周围的流线与温度场的分布有显著差异; 由于表面张力作用, 微重力下的气泡脱离特性与常重力下不同; 在微重力条件下, 气泡直径的变化比较复杂, 并与重力加速度的大小有关; Marangoni流对微重力下的流动影响很大, 使换热系数波动, 而且波动的幅度随重力加速度的减小而增大.      

二元喉道倾斜矢量喷管调节方法

利用数值模拟方法,对二元喉道倾斜矢量喷管的调节方法进行了研究。研究结果表明,当喉道注气流量不变,随着扩张段注气流量逐渐增大,喷管内出现三种典型流动状态。当扩张段注气流量比较小时,喉道处声速线的形状和位置没有明显变化,在扩张段注气口前面形成亚声速区域;当喷管扩张段注气流量比较大时,喷管中出现了强激波系和上下贯通的亚声速区域,这时喷管内的流动损失最大,推力系数最低;进一步增大扩张段注气流量,声速线会从喷管喉道移动到喉道注气口和扩张段注气口之间,出现典型的喉道倾斜现象。单独控制流量系数可通过保持扩张段注气流量不变,改变喉道注气流量来实现;矢量角的控制可通过保持喉道注气流量不变,改变扩张段注气流量来实现。     

Don Pettit 国际空间站微重力流体实验与若干现象的解析

Don Pettit 是美国NASA 的资深宇航员, 2002 年11 月至2003 年5 月, 在国际空间站工作期间, 进行了一系列名为Saturday Morning Science的微重力实验. Don Pettit 实验揭示了微重力下一些令人惊奇的现象, 涉及流体力学、气泡动力学、晶体生长以及运动学的诸多基础问题, 是能够获得的国际空间站为数不多的微重力实验资料, 具有重要研究意义. 实验包括14 组流体实验和7 组力学实验, 而流体实验又大致分为水膜实验和水球实验, 包括大尺寸水膜制备、扩散和对流、Marangoni 对流、薄膜结晶、水球/气泡旋转等. 本文对Don Pettit 实验的流体部分进行了介绍, 并对其中的Marangoni对流和薄膜结晶实验现象进行了解析.