Mach reflection in steady supersonic flow considering wedge boundary-layer correction

Mach reflection in steady supersonic flow is an important phenomenon having received extensive studies, among which simplified theoretical models to predict the size of Mach stem and other flow structure are of particular interest. Past efforts for such models were based on inviscid assumption while in real cases the flow is viscous. Here in this paper we consider the influence of wedge boundary layer on the Mach stem height. This is done by including a simplified boundary layer model into a recently published inviscid model. In this viscous model, the wedge angle and the trailing edge height, which control the Mach stem height, are replaced by their equivalent ones accounting for the displacement effect of the wedge boundary layer, with the boundary layer assumed to be laminar or fully turbulent. This viscous model is shown to compare well with numerical results by computational fluid dynamics and gives a Mach stem height as function of the Reynolds number and Mach number. It is shown that due to the viscous effect, the Mach stem height is increased, through increasing the effective wedge angle.     

Interaction mechanisms of shock waves with the boundary layer and wakes in a highly-loaded NGV using hybrid RANS/LES

Accurate predictions of Shock Waves and Boundary Layer Interaction (SWBLI) and strong Shock Waves and Wake Vortices Interaction (SWWVI) in a highly-loaded turbine propose challenges to the currently widely used Reynolds-Averaged Navier-Stokes (RANS) model. In this work, the SWBLI and the SWWVI in a highly-loaded Nozzle Guide Vane (NGV) are studied using a hybrid RANS/LES strategy. The Turbulence Kinetic Energy (TKE) budget and the Proper Orthogonal Decomposition (POD) method are used to analyze flow mechanisms. Results show that this hybrid RANS/LES method can obtain detailed flow structures for flow mechanisms analysis. Strong shock waves induce boundary layer separation, while the presence of a separation bubble can in turn lead to a Mach reflection phenomenon. The shock waves cause trailing-edge vortices to break clearly, and the wakes, in turn, can change the shocks intensity and direction. Furthermore, the Entropy Generation Rate (EGR) is used to analyze the irreversible loss. It turns out that the SWWVI can reduce the flow field loss. There are several weak shock waves in the NGV flow field, which can increase the irreversible loss. This work offers flow mechanisms analysis and presents the EGR distribution in SWBLI and SWWVI areas in a transonic turbine blade.     

Simulator of pyroshock environment and effect rules of its adjustable parameters

During the launching of spacecraft, the on-board devices will undergo a series of pyroshock environments. In order to verify the reliability of these devices under these pyroshock environments, all of them are needed to take the shock test before launching. This paper has carried out an in-depth research on the simulation method of the pyroshock based on the true explosive excitation. In this study, a simulator containing multiple adjustment parameters is presented and the safety is considered by the design of the protective cover. And the working process of this setup is simulated with the explicit dynamic codes LS-DYNA. What’s more, the effects of the adjustment parameters on the three factors of shock Response Spectrum (SRS) of the resonant board are explored carefully. The rules achieved in this paper are verified by a typical example. The results indicate that the improved simulator can avoid the danger of explosive and make full use of the advantage of actual explosive excitation. And the test condition can be quickly realized at the simulator according to the effect rules of the three adjustable parameters.     

Kink-like mode of a double gradient instability in a compressible plasma current sheet

A linear MHD instability of the electric current sheet, characterized by a small normal magnetic field component, varying along the sheet, is investigated. The tangential magnetic field component is modeled by a hyperbolic function, describing Harris-like variations of the field across the sheet. For this problem, which is formulated in a 3D domain, the conventional compressible ideal MHD equations are applied. By assuming Fourier harmonics along the electric current, the linearized 3D equations are reduced to 2D ones. A finite difference numerical scheme is applied to examine the time evolution of small initial perturbations of the plasma parameters. This work is an extended numerical study of the so called “double gradient instability”, – a possible candidate for the explanation of flapping oscillations in the magnetotail current sheet, which has been analyzed previously in the framework of a simplified analytical approach for an incompressible plasma. The dispersion curve is obtained for the kink-like mode of the instability. It is shown that this curve demonstrates a quantitative agreement with the previous analytical result. The development of the instability is investigated also for various enhanced values of the normal magnetic field component. It is found that the characteristic values of the growth rate of the instability shows a linear dependence on the square root of the parameter, which scales uniformly the normal component of the magnetic field in the current sheet.     

用于爆震引燃的激波聚焦无反应流场数值模拟

通过对用于脉冲爆震引燃的内楔表面激波聚焦冷态流场的数值模拟，分析了不同角度内楔表面在激波聚焦时流场结构的不同特点。数值计算的结果表明：对于大内楔角的情况，激波进入内楔角之后产生规则反射，温度的大幅度升高发生在激波于内楔顶点反射之后；对于小内楔角的情况，激波进入内楔角之后产生马赫反射，温度的升高在入射激波到达内楔顶点之前就已经发生，而在入射激波于内楔顶点反射之后得到进一步的提升。同时，还对不同角度内楔的反射激波后气流温度进行了比较。     

激光推进中飞行器结构响应对聚焦性能的影响

利用ABAQUS软件对激光推进飞行器进行了单、多脉冲结构冲击响应的数值模拟;建立了变形抛物面焦点区域的数值计算方法,计算出结构变形后的焦点区域的变化程度:发现因结构变形引起的聚焦区域的变化不可忽视,是飞行器设计时必须考虑的因素。     

聚心火焰与诱导激波相互作用及爆燃转爆轰过程

基于带化学反应的二维轴对称Eu ler方程,利用带限制函数的波传播算法,在两端均敞开和一端封闭一端敞开的两种边界情况下,对圆柱型燃烧室中的氢气-空气预混气的聚心燃烧进行了数值模拟,探讨了聚心火焰引发爆轰过程中的波系结构以及火焰的变形。数值结果表明,在两端均敞开的情况下,燃烧诱导的激波在轴心、火焰和壁面之间来回反射,使火焰失稳,加速了火焰的燃烧,但没有转变成爆轰。而在一端封闭的情况下,封闭端反射产生的激波,不断穿越火焰,使火焰严重失稳,加剧了燃烧速度,最终导致爆轰的形成。稳定爆轰阵面的参数与CJ理论计算的结果进行了对比,两者符合的较好。同时,火焰在与激波的作用过程中,形状扭曲变形,对于两端敞开的情况呈对称的扁平头部蘑菇云,而对于一端封闭的情况,则呈封闭端小敞口端大的扁平头部蘑菇云。     

连续式高焓超声速预混加热器设计与分析

为研究超声速燃烧和爆轰相关的机理问题,提出了一种结合燃烧型加热器和阵列喷管的超声速预混加热器设计思想。通过预热燃烧室来提供总温可变的高焓富氧气流,经过特征线型面喷管膨胀降温后,在喷管扩张段的适当位置以一定角度喷入燃料,经过混合段后形成所需的连续高焓总温和当量比可调的预混气流。通过对混合过程的数值模拟和预混气体的着火延迟时间分析了当前的预混高焓加热器的混合和自燃问题。在超声速气流中加入斜劈采用纹影技术进行激波点火实验,并验证了当前的预混加热器设计是成功。     

Propagation of interplanetary shocks into the Earth’s magnetosphere

This paper is devoted to the study of propagation of disturbances caused by interplanetary shocks (IPS) through the Earth’s magnetosphere. Using simultaneous observations of various fast forward shocks by different satellites in the solar wind, magnetosheath and magnetosphere from 1995 till 2002, we traced the interplanetary shocks into the Earth’s magnetosphere, we calculated the velocity of their propagation into the Earth’s magnetosphere and analyzed fronts of the disturbances. From the onset of disturbances at different satellites in the magnetosphere we obtained speed values ranging from 500 to 1300 km/s in the direction along the IP shock normal, that is in a general agreement with results of previous numerical MHD simulations. The paper discusses in detail a sequence of two events on November 9th, 2002. For the two cases we estimated the propagation speed of the IP shock caused disturbance between the dayside and nightside magnetosphere to be 590 km/s and 714–741 km/s, respectively. We partially attributed this increase to higher Alfven speed in the outer magnetosphere due to the compression of the magnetosphere as a consequence of the first event, and partially to the faster and stronger driving interplanetary shock. High-time resolution GOES magnetic field data revealed a complex structure of the compressional wave fronts at the dayside geosynchronous orbit during these events, with initial very steep parts (10 s). We discuss a few possible mechanisms of such steep front formation in the paper.     

惯导系统两类冲击环境条件之间的等效转换

GJB150A采用冲击响应谱的形式来描述复杂冲击环境条件,有助于提升惯导系统等航天电子设备地面试验的真实性。但冲击响应谱试验结果的较大差异性也给产品研制过程中的试验研究带来很大困扰。为了解决试验真实性与试验结果较大差异性之间的矛盾,提出了一种将复杂冲击条件转换为经典冲击条件的等效方法,用经典冲击波形等效冲击响应谱,给出了等效转换的基本准则,基于数值试验结果导出了等效转换公式。经试验验证,按照该方法转换得到的经典波形冲击试验结果接近多次冲击响应谱冲击试验结果的平均水平,表明该方法有效。