多机并联火箭羽流流场及其底部热环境分析

For the problem that the plume flow field structure of a multi engine parallel rocket is complicated and the bottom thermal environment is extremely harsh, which may cause the failure of the engine structural components, the plume flow field and thermal environment at different altitudes are studied through numerical simulation. The result is compared with the measured results in flight which shows that when the rocket is flying at a low altitude, the plume of the engines do not interfere with each other. As the flight altitude increases, the plumes gradually expand and begin to interfere with each other, and finally there is an obvious backflow at the bottom of the rocket. The maximum heat flux at the moment of take off is basically the same as the measured value in flight. Before the backflow occurs, the heat flux mainly consists of radiant heat, the convective heat flow increases as the flight altitude grows, but it is also much smaller than the peak heat flow at takeoff. The result has certain guiding significance for the optimal design of engine structure thermal protection.     

The verification of GNSS tropospheric tomography model in a mountainous area

The GNSS (Global Navigation Satellite System) has not been developed as a meteorological data source provider, but with a careful and sophisticated processing strategy it might be used as one. The term GNSS tomography refers to the usage of the ray traced GNSS signal as scanning rays in the tomographic model input. The model is divided into a number of voxels. The system is inverted and value of refractivity is obtained. Typically, as in the most of the inverse processing, there is a problem of the undetermined system and as a consequence the cofactor matrix is close to singular. To avoid singularity additional conditions or constrains should be added to the system. Here, additional parameters are derived with the help of the air flow analysis in the Sudety mountains (south-west region of Poland), and special Slant Wet Delay (SWD) trimming procedure. The flow’s synthetic parameters like the Bruint-Väisälä frequency and the Froude number are determined. This way the type of the flow is recognized and the analysis of the impact of orographic barrier has been quantified. The SWDs from the GNSS observations were tested against, SWD from raytracing through the COAMPS model field. The modified GNSS tomography model was tested for the real GNSS observations delivered from the GNSS network Karkonosze located in the Sudety mountains and compared with the COAMPS model. The solution shows a considerable improvement in comparison with plain tomographic model results.     

三角翼DBD等离子体流动控制研究进展

现代战机采用较多的三角翼，在大迎角绕流时存在前缘涡破裂等气动问题。作为新型主动流动控制技术，等离子体激励频带宽、响应快、结构简单、便于闭环控制，在解决三角翼气动问题上具有潜力。回顾了介质阻挡放电（DBD）等离子体气动激励的基本原理，及其用于三角翼前缘涡控制的研究进展。从来流条件、几何构型、激励参数等方面分析了DBD等离子体激励对流动控制效果的影响规律；结合不同激励频率下流场演化特性，分析了流动控制机理。最后，从理论研究和工程应用的角度，对三角翼前缘涡控制的发展进行总结展望。     

Numerical simulations of water spray on flame deflector during the four-engine rocket launching

This paper aims to study on the cooling effect with two types of water spray nozzle on the flame deflector during the four-engine launch vehicle take-off. To accurately simulate the two-phase flow of the rocket gas with multispecies and the water spray, the three-dimensional compressible Navier-Stokes equations, discrete ordinates methods and realizable k-Ɛ turbulence model are used to establish the rocket supersonic plumes impact model. The Eulerian dispersed phase (EDP) model was used to simulate the water spray into the exhaust gas. The accuracy and effectiveness of the gas-liquid flow model are verified by a good agreement between simulation results and experimental data. On this basis, a series of numerical simulation studies under different water injection position are performed. The results show that the high temperature regions, along the axis of engines on the deflector plate, have no significant temperature decreasing effect by water spray from the nozzles mounted on the apex of the deflector, and the high temperature converts a large quantity of water into vapor near the plume boundary, which would decrease the flow conductivity. With the cooling spray nozzle fixed directly to the deflector plate, the temperature decrease effect is obvious and the effect of thermal shock on deflector plate induced by exhaust plume is reduced, so that it can prevent the flame deflector from thermal ablation. The study results provide indepth information and engineering guidance for designing the water spray systems and increasing the safety of the launch process.     

基于微型涡喷发动机热喷流的无源流体推力矢量喷管的控制规律

流体推力矢量喷管型面固定、活动部件少、结构重量轻，能够为高机动飞行器提供有效的飞行控制手段，但无源流体推力矢量喷管热喷流的偏转控制规律尚未完全掌握。为了推进无源流体推力矢量技术的实用化，本文设计研制了适用于微型涡喷发动机的耐高温喷管模型，对该喷管在微型涡喷发动机热喷流状态下的控制规律进行研究。利用非接触光学显示和测量手段——红外热成像拍摄和粒子图像测速（PIV）技术对主射流流动特性进行研究，获得流动矢量角随二次流控制阀门闭合度变化的控制规律；利用六分量盒式天平测力实验研究无源流体推力矢量喷管的力学特性，获得推力矢量角随二次流控制阀门闭合度变化的控制规律。研究结果表明：该构型喷管在微型涡喷发动机热喷流下主射流连续可控偏转，最大流动矢量角为-12.3°/12.3°，最大推力矢量角为-12.9°/12.8°，控制规律接近线性，不存在主射流偏转突跳问题。     

A comprehensive study of the 2016 Mw 6.0 Italy earthquake based on high-rate (10 Hz) GPS data

The 2016 M_w 6.0 Italy earthquake is successfully recorded by the near-field 10?Hz GPS and 200?Hz Strong Motion (SM) stations, providing valuable data for this study. A comprehensive study of this earthquake is carried out based on GPS data, which contains coseismic deformations analysis, noise analysis, seismic wave picking, and magnitude determination. The noise of most GPS-derived displacement waveforms can be described as a combination of white noise, flicker noise, and random walk noise after the earthquake occurrence, and the spectral indices vary significantly for most stations, implying that the seismic signals have affected the noise characteristic of GPS-derived displacement waveforms. S-transform is employed to assess the GPS capability to detect the seismic arrival time. The SM station AMT and the GPS station AMAT are in good agreement in seismic wave picking, and the difference is only 1.2?s in the north component, suggesting that the outcome of seismic wave picking using GPS data is reliable. Then, a classic empirical formula is employed to determine the moment magnitude. A robust moment magnitude (Mw 5.90) can be estimated by the nine GPS stations with about 23.9?s. If four GPS stations near the epicenter is chosen to determine the magnitude, it only take 13.0?s to retrieve a reliable preliminary (Mw 5.82) magnitude, which is 5.4?s ahead of nine stations. In addition, Cross Wavelet Transform (XWT) is adopted to measuring the correlation and phase relationship between GPS and SM records. The result of XWT analysis indicates 10?Hz GPS is capable of capturing reliable and accurate coseismic dynamic deformations, as evidenced by the XWT-based semblance being close to 1 between GPS and SM records. The above results confirm the capability of 10?Hz GPS to capture coseismic dynamic deformations, detect seismic arrival time, and determine earthquake magnitude. Moreover, rapid magnitude determination based on 10?Hz GPS data can be regarded as an important supplement to Earthquake Early Warning (EEW).     

Abnormal changes of dynamic derivatives at low reduced frequencies

A number of studies have found that abnormal changes of dynamic derivatives occurred at very low reduced frequencies, but its inducement mechanism is not very clear. This paper has researched the abnormal changes and analyzed the influence on some parameters by solving the unsteady flow around forced oscillation airfoils based on Navier-Stokes equations. Results indicate that when the reduced frequency approaches to zero, the dynamic derivatives obtained by the numerical method will diverge. We have also proven it in theory that this phenomenon is not physical but completely caused by numerical singularity. Furthermore, the abnormal phenomenon can be effectively mitigated by using the time spectral method to solve the aerodynamic forces and the integral method to obtain the dynamic derivatives. When the reduced frequency is in the range of 0.001–0.01, the dynamic derivative maintains nearly unchanged for the whole speed region. This study can provide a reference for the reasonable choice of the reduced frequency in calculations and experiments of dynamic derivatives.     

Geometric calibration of Colour and Stereo Surface Imaging System of ESA’s Trace Gas Orbiter

There are many geometric calibration methods for “standard” cameras. These methods, however, cannot be used for the calibration of telescopes with large focal lengths and complex off-axis optics. Moreover, specialized calibration methods for the telescopes are scarce in literature. We describe the calibration method that we developed for the Colour and Stereo Surface Imaging System (CaSSIS) telescope, on board of the ExoMars Trace Gas Orbiter (TGO). Although our method is described in the context of CaSSIS, with camera-specific experiments, it is general and can be applied to other telescopes. We further encourage re-use of the proposed method by making our calibration code and data available on-line.     

Numerical investigation on flow nonuniformity-induced hysteresis in scramjet isolator

Numerical simulation and theoretical analysis were conducted to study the hysteresis inside scramjet isolator during the reciprocating process of back pressure variation. It is revealed that only a regular reflection is theoretically possible for two leading shocks when the inflow Mach number is greater than 2.0, and no hysteresis can occur in the transition between shock reflection types. Nevertheless, wall suction, gas injection, and background waves cause non-uniformity of the incoming flow and would make hysteresis possible. Besides the classical hysteresis in the transition between shock reflection, new kinds of hysteresis were found in both the deflection angle of separated boundary layer and the location of the shock train. Moreover, the occurrence of hysteresis in the deflection angle of the separated boundary layer is accompanied with the shock reflection hysteresis. In the case with background waves or gas injection, hysteresis in the starting position of leading shock was observed too. As back pressure decreases, the leading shock does not follow the same path as that as the back pressure increases, and it is anchored at the location where the background shock or the injection interacts with the leading shock. It is inferred that, if two strong adverse pressure gradient regions move towards and interact with each other, hysteresis will take place when they start to separate.     

Flow mechanism between purge flow and mainstream in different turbine rim seal configurations

The rim seal is used to prevent mainstream ingestion to the gap between the vane of a turbine and its blade. In this article, the dolphin lip with a hook configuration and a large seal cavity with hook structures are designed based on the high-pressure turbine datum single shark lip rim seal configuration. The sealing effect and parameters of the flow field are measured by an experiment method and a numerical simulation is used to explain the mechanism. For three configurations, the effect of the leakage slot vortex on the efficiency of the seal and the influence of leakage vortex, generated by the interaction between purge flow and mainstream flow, are discussed in depth. The result shows that the reverse vortex formed by the dolphin lip rim seal with hook structure will increase the sealing efficiency. The seal configuration with a large cavity improves sealing efficiency to a greater extent than the datum structure. At different purge flow rates and with unequal seal structures, the purge flow produces three types of leakage vortices in the passage. Besides, the seal configuration with dolphin lip produces a Kelvin-Helmholtz instability at the interface of the purge and the mainstream flows at a low purge flow rate to induce new leakage vortex branches in the passage of the blade.