Modelling ignition probability with pairwise mixing-reaction model for flame particle tracking

Reduced order models for ignition analysis can offer insights into ignition processes and facilitate the combustor optimization. In this study, a Pairwise Mixing-Reaction (PMR) model is formulated to model the interaction between the flame particle and the surrounding cell mixture during Lagrangian flame particle tracking. Specifically, the model accounts for the two-way coupling of mass and energy between the flame particle and the surrounding shell layer by modelling the corresponding turbulent mixing, chemical reaction and evaporation process if present. The state of a flame particle, e.g., burnt, hot gas or extinguished, is determined based on particle temperature. This model can properly describe the ignition process with a spark kernel being initiated in a nonflammable region, which is of practical importance in certain turbine engines and has not been rigorously accounted for by the existing models based on the estimation of local Karlovitz number. The model is integrated into an ignition probability analysis platform and is demonstrated for a methane/air bluff-body flame with the flow and fuel/air mixing characteristics being extracted from a non-reacting simulation. The results show that for the spark location being at the extreme fuel-lean outer shear layer of the recirculation zone, PMR can yield ignition events with a significant number of active flame particles. The mechanisms for the survival of the initial flame particles and the entrainment of the survived flame particles into the recirculation zone are analyzed. The results also show that the ignition probability map from PMR agrees well with the experimental observation: a high ignition probability in the shear layer of the recirculation zone near the mean stoichiometric surface, and low ignition probabilities inside the recirculation zone and the top stagnation region of the recirculation zone. The parametric study shows that the predicted shape of the ignition progress factor and ignition probability is in general insensitive to the model parameters and the model is adequate for quantifying the regions with high ignition probabilities.     

The adjusted optical properties for Galileo/BeiDou-2/QZS-1 satellites and initial results on BeiDou-3e and QZS-2 satellites

Solar Radiation Pressure (SRP) is the dominant non-gravitational perturbation for GNSS (Global Navigation Satellite System) satellites. In the absence of precise surface models, the Empirical CODE Orbit Models (ECOM, ECOM2) are widely used in GNSS satellite orbit determination. Based on previous studies, the use of an a priori box-wing model enhances the ECOM model, especially if the spacecraft is a stretched body satellite. However, so far not all the GNSS system providers have published their metadata. To ensure a precise use of the a priori box-wing model, we estimate the optical parameters of all the Galileo, BeiDou-2, and QZS-1 (Quasi Zenith Satellite System) satellites based on the physical processes from SRP to acceleration. Validation using orbit prediction proves that the adjusted parameters of Galileo and QZS-1 satellites exhibit almost the same performance as the corresponding published and “best guess” values. Whereas, the estimated parameters of BeiDou-2 satellites demonstrate an improvement of more than 60% over the initial “guess” values. The resulting optical parameters of all the satellites are introduced into an a priori box-wing model, which is jointly used with ECOM and ECOM2 model in the orbit determination. Results show that the pure ECOM2 model exhibits better performance than the pure ECOM model for Galileo, BeiDou-2 GEO and QZS-1 orbits. Combined with the a priori box-wing model the ECOM model (ECOM+BW) results in the best Galileo, BeiDou-2 GEO and QZS-1 orbits. The standard deviation (STD) of satellite laser ranging residuals reduce by about 20% and 5% with respect to the pure ECOM2 model for Galileo and BeiDou-2 GEO orbits, while the reductions are about 40% and 60% for QZS-1 orbits in yaw-steering and orbit-normal mode respectively. BeiDou-2 IGSO and MEO satellite orbits do not benefit much from the a priori box-wing model. In summary, we suggest setting up a unified SRP model of ECOM+BW for Galileo, QZS-1, and BeiDou-2 orbits based on the adjusted metadata. In addition, we estimate the optical parameters of BeiDou-3e and QZS-2 satellites using a limited number of tracking stations. Results regarding the unified SRP model indicate the same advantages, the STD of satellite laser ranging residuals reduces by about 30% and 20% for QZS-2 and BeiDou-3e orbits respectively over orbit products without a priori model. The estimation procedure is effective and easy to apply to the new emerging satellites in the future.     

Comparison of observed ionospheric vertical TEC over the sea in Indian region with IRI-2016 model

The vertical ionospheric TEC values obtained from GAGAN grid based ionospheric delay correction values over the sea in the Indian equatorial region have been compared with the corresponding values derived from the International Reference Ionosphere model, IRI-2016. The objective of this work is to study the deviation of the vertical TEC derived from the IRI model from ground truths over the sea for different conditions. This will serve the basic intention of assessing the candidature of the IRI model as an alternative ionospheric correction model in navigation receivers in terms of accuracy. We have chosen different solar activity periods, seasons, geomagnetic conditions, locations etc. for our comparison and analysis. The TEC values by the IRI-2016 were compared with the actual measured values for the given conditions and errors were obtained. The measured vertical TEC values at the ionospheric grid points were derived from the GAGAN broadcast ionospheric delay data and used as reference. The IRI model with standard internal functions was used in estimating the TEC at the same ionospheric grid points. The errors in the model derived values are statistically analysed. Broadly, the results show that, for the Indian sector over the sea, the IRI model performs better on quiet days in off equatorial regions, particularly in the northern region. The overall performance degrades for other conditions with the model generally underestimating the true TEC values and most severely in the equatorial region. The performance is worst in this region for the disturbed days of the equinoctial period. The comparison study is also done with the TEC data measured directly by dual frequency GPS receivers. The results were found to be in general agreement with those obtained by comparing the model with GAGAN broadcast data as reference. This study will be useful in considering the IRI-2016 model for real time estimates of TEC as an alternative to the current parametric model in a satellite navigation receiver in absence of other options.     

A modern robust approach to remotely estimate chlorophyll in coastal and inland zones

The chlorophyll concentration of a water body is an important proxy for representing the phytoplankton biomass. Its estimation from multi or hyper-spectral remote sensing data in natural waters is generally achieved by using (i) the waveband ratioing in two or more bands in the blue-green or (ii) by using a combination of the radiance peak position and magnitude in the red-near-infrared (NIR) spectrum. The blue-green ratio algorithms have been extensively used with satellite ocean color data to investigate chlorophyll distributions in open ocean and clear waters and the application of red-NIR algorithms is often restricted to turbid productive water bodies. These issues present the greatest obstacles to our ability to formulate a modern robust method suitable for quantitative assessments of the chlorophyll concentration in a diverse range of water types. The present study is focused to investigate the normalized water-leaving radiance spectra in the visible and NIR region and propose a robust algorithm (Generalized ABI, GABI algorithm) for chlorophyll concentration retrieval based on Algal Bloom index (ABI) which separates phytoplankton signals from other constituents in the water column. The GABI algorithm is validated using independent in-situ data from various regional to global waters and its performance is further evaluated by comparison with the blue-green waveband ratios and red-NIR algorithms. The results revealed that GABI yields significantly more accurate chlorophyll concentrations (with uncertainties less than 13.5%) and remains more stable in different waters types when compared with the blue-green waveband ratios and red-NIR algorithms. The performance of GABI is further demonstrated using HICO images from nearshore turbid productive waters and MERIS and MODIS-Aqua images from coastal and offshore waters of the Arabian Sea, Bay of Bengal and East China Sea.     

Simulating planetary wave propagation to the upper atmosphere during stratospheric warming events at different mountain wave scenarios

Parameterization schemes of atmospheric normal modes (NMs) and orographic gravity waves (OGWs) have been implemented into the mechanistic Middle and Upper Atmosphere Model (MUAM) simulating atmospheric general circulation. Based on the 12-members ensemble of runs with the MUAM, a composite of the stratospheric warming (SW) has been constructed using the UK Met Office data as the lower boundary conditions. The simulation results show that OGW amplitudes increase at altitudes above 30 km in the Northern Hemisphere after the SW event. At altitudes of about 50 km, OGWs have largest amplitudes over North American and European mountain systems before and during the composite SW, and over Himalayas after the SW. Simulations demonstrate substantial (up to 50–70%) variations of amplitudes of stationary planetary waves (PWs) during and after the SW in the mesosphere-lower thermosphere of the Northern Hemisphere. Westward travelling NMs have amplitude maxima not only in the Northern, but also in the Southern Hemisphere, where these modes have waveguides in the middle and upper atmosphere. Simulated variations of PW and NM amplitudes correspond to changes in the mean zonal wind, EP-fluxes and wave refractive index at different phases of the composite SW events. Inclusion of the parameterization of OGW effects leads to decreases in amplitudes (up to 15%) of almost all SPWs before and after the SW event and their increase (up to 40–60%) after the SW in the stratosphere and mesosphere at middle and high northern latitudes. It is suggested that observed changes in NM amplitudes in the Southern Hemisphere during SW could be caused by divergence of increased southward EP-flux. This EP-flux increases due to OGW drag before SW and extends into the Southern Hemisphere.     

A novel satellite image encryption algorithm based on hyperchaotic systems and Josephus problem

With the rapid growth of the number of Earth observation satellite (EOS) supporting critical applications, it is required to improve the security techniques to protect the sensitive data and images during the transmission between the satellites and the ground stations. This paper introduces a new satellite image encryption algorithm based on the Linear Feedback Shift Register (LFSR) generator, SHA 512 hash function, hyperchaotic systems, and Josephus problem. LFSR generates a matrix that is used to construct the 512-bits value of the hash function. These bits are used to set the initial values and parameters of the proposed encryption algorithm. Firstly, the six dimensions (6-D) hyperchaotic system is divided into three parts, where every two equations are considered as one part. Secondly, the 1-D hyperchaotic logistic-tent system is considered as the controller to select one part. The selected part is used to generate a matrix that is XORed with the original image. Thirdly, the scrambling operation by Josephus sequences is applied to the output of the previous step by scrambling the rows and the columns according to the selected part to produce the pre-encrypted image. Finally, if the number of iterations is less than the required number which is considered as a parameter of the secret key, the previous operations will be repeated in the pre-encrypted image; otherwise, the pre-encrypted image is considered as the final cipher image. Experimental and analyses results show that the proposed algorithm has good performance in terms of high level of security, large enough key-space, tolerance to Single Event Upsets (SEU) as well as low time complexity.     

航天器姿态预设性能控制方法综述

针对存在不确定惯量矩阵及其它未知不确定性(如执行器故障、航天器结构发生突变等)下的航天器姿态控制问题，综述了预设性能控制及其应用的研究进展。首先阐释了预设性能控制方法的基本内涵及其关键步骤；然后分析总结了现有航天器姿态控制以及预设性能控制研究的基本概况与发展趋势；最后面向未来复杂化、智能化的空间任务对航天器控制系统提出的新需求，提出航天器预设性能控制值得研究的问题和方向。     

基于数字孪生的智能脉动管控

为提升飞机总装脉动生产线智能化水平,开展了基于数字孪生的智能脉动管控研究。基于数字孪生五维模型理论构建了飞机总装脉动生产线数字孪生应用架构,从反应式计划调度、物流精准配送、智能作业指导、产线健康监视、资源迭代优化配置5个方面开展了技术研究,为复杂产品装配的智能制造应用提供参考。     

Maximal admissible mode decision delay in terminal guidance

Mode decision-maker is a critical component in the logic-based Integrated Estimation and Guidance(IEG) system. For the best possible estimation and guidance performance, the mode decision delay of the mode decision-maker should be limited to a range as small as possible. This paper presents a numerical method for computing the maximal admissible mode decision delay that varies with time-to-go. Particular attention has been paid to highly maneuvering target interception in terminal guidance. The results of this research offer useful guidelines for the design of the mode decision-maker in IEG systems.     

Orbit determination of BDS-3 satellite based on regional ground tracking station and inter-satellite link observations

The BeiDou global navigation satellite system (BDS-3) has established the Ka-band inter-satellite link (ISL) to realize a two-way ranging function between satellites, which provides a new observation technology for the orbit determination of BDS-3 satellites. Therefore, this study presents a BDS satellite orbit determination model based on ground tracking station (GTS) observations and ISL ranging observations firstly to analyze the impact of the ISL ranging observations on the orbit determination of BDS-3 satellites. Subsequently, considering the data fusion processing, the variance component estimation (VCE) algorithm is applied to the parameter estimation process of the satellite orbit determination. Finally, using the measured data from China’s regional GTS observations and BDS-3 ISL ranging observations, the effects of ISL ranging observations on the orbit determination accuracy of BDS-3 satellites are analyzed. Moreover, the impact of the VCE algorithm on the fusion data processing is evaluated from the aspects of orbit determination accuracy, Ka-band hardware delay parameter stability, and ISL ranging observation residuals. The results show that for China’s regional GTSs, the addition of BDS-3 ISL ranging observations can significantly improve the orbit determination accuracy of BDS-3 satellites. The observed orbit determination accuracy of satellite radial component is improved from 48 cm to 4.1 cm. In addition, when the initial weight ratio between GTS observations and ISL ranging observations is not appropriate, the various indicators which include orbit determination accuracy, ISL hardware delay, and ISL observation residuals were observed to have improved after the adjustment of the VCE algorithm. These results validate the effectiveness of the VCE algorithm for the fusion data processing of the GTS observations and ISL ranging observations.