Assessment of the close approach frequency and collision probability for satellites in different configurations of large constellations

The Earth orbital environment is drastically changing due to an intensification of the space activities. In particular, several projects of large constellations, proposed for the next years for communications purpose like global internet access, Internet of Things, or for Earth observations, will lead to the deployment of several thousands of new satellites at an unprecedented rate. It is a crucial challenge for space traffic management, which will deal with a great number of satellite conjunctions, potentially causing a collision with damaging consequences for the constellation itself and the space environment sustainability.In this paper, we investigate the close approach frequency and the cumulative collision probability for each referenced constellation. For this purpose, we compute the orbital evolution of satellites in different constellations during the lifecycle, from the deployment to the decommissioning, and we apply the CUBE algorithm and the Foster method to assess the collision probability with the background space debris population assuming a constant uncertainty in position. We show the variation of risk defined by the close approach frequency and the cumulative collision probability as a function of the proposed configuration. In particular, satellites of the Iridium and Kuiper constellation, but also satellite of the Telesat constellation on polar orbits are the most exposed at a collision. Moreover, the decommissioning phase contribute for a major part to the final cumulative collision probability.     

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.     

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

基于PTE的多雷达航迹融合算法

针对杂波背景下的多雷达航迹融合时局部估计误差互协方差矩阵未知的问题，提出基于目标存在概率（PTE）的航迹融合算法，提升了正确航迹率和跟踪精度。首先，通过综合概率数据关联得到单接收站的目标航迹估计集合和对应的目标存在概率。然后，在局部估计误差互协方差矩阵未知的条件下，基于PTE信息提出不带记忆的综合广义凸组合航迹融合算法。进而将前一帧的融合状态进行反馈，提出带记忆的综合广义凸组合航迹融合算法。仿真验证了所提算法的有效性。     

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.     

Fault mode probability factor based fault-tolerant control for dissimilar redundant actuation system

This paper presents a Fault Mode Probability Factor (FMPF) based Fault-Tolerant Control (FTC) strategy for multiple faults of Dissimilar Redundant Actuation System (DRAS) composed of Hydraulic Actuator (HA) and Electro-Hydrostatic Actuator (EHA). The long-term service and severe working conditions can result in multiple gradual faults which can ultimately degrade the system performance, resulting in the system model drift into the fault state characterized with parameter uncertainty. The paper proposes to address this problem by using the historical statistics of the multiple gradual faults and the proposed FMPF to amend the system model with parameter uncertainty. To balance the system model precision and computation time, a Moving Window (MW) method is used to determine the applied historical statistics. The FMPF based FTC strategy is developed for the amended system model where the system estimation and Linear Quadratic Regulator (LQR) are updated at the end of system sampling period. The simulations of DRAS system subjected to multiple faults have been performed and the results indicate the effectiveness of the proposed approach.     

27-day variation in solar-terrestrial parameters: Global characteristics and an origin based approach of the signals

The Earth and the near interplanetary medium are affected by the Sun in different ways. Those processes generated in the Sun that induce perturbations into the Magnetosphere-Ionosphere system are called geoeffective processes and show a wide range of temporal variations, like the 11-year solar cycle (long term variations), the variation of ～27?days (recurrent variations), solar storms enduring for some days, particle acceleration events lasting for some hours, etc.In this article, the periodicity of ～27?days associated with the solar synodic rotation period is investigated. The work is mainly focused on studying the resulting 27-day periodic signal in the magnetic activity, by the analysis of the horizontal component of the magnetic field registered on a set of 103 magnetic observatories distributed around the world. For this a new method to isolate the periodicity of interest has been developed consisting of two main steps: the first one consists of removing the linear trend corresponding to every calendar year from the data series, and the second one of removing from the resulting series a smoothed version of it obtained by applying a 30-day moving average. The result at the end of this process is a data series in which all the signal with periods larger than 30?days are canceled.The most important characteristics observed in the resulting signals are two main amplitude modulations: the first and most prominent related to the 11-year solar cycle and the second one with a semiannual pattern. In addition, the amplitude of the signal shows a dependence on the geomagnetic latitude of the observatory with a significant discontinuity at approx. ±60°.The processing scheme was also applied to other parameters that are widely used to characterize the energy transfer from the Sun to the Earth: F10.7 and Mg II indices and the ionospheric vertical total electron content (vTEC) were considered for radiative interactions; and the solar wind velocity for the non-radiative interactions between the solar wind and the magnetosphere. The 27-day signal obtained in the magnetic activity was compared with the signals found in the other parameters resulting in a series of cross-correlations curves with maximum correlation between 3 and 5?days of delays for the radiative and between 0 and 1?days of delay for the non-radiative parameters. This result supports the idea that the physical process responsible for the 27-day signal in the magnetic activity is related to the solar wind and not to the solar electromagnetic radiation.     

基于Kriging模型梯度解析解的改进一次二阶矩方法

改进一次二阶矩（AFOSM）法是一种基于功能函数梯度的结构可靠性分析方法，鉴于其对隐式函数的梯度较难求解，提出了一种基于Kriging模型梯度解析解的AFOSM方法，利用Kriging代理模型的解析表达式推导得到功能函数对输入变量的梯度解析解，为AFOSM中设计点的确定提供高精度的梯度信息。通过Kriging与AFOSM的结合，很好地解决了基于有限元模型的隐式情况下梯度计算量相当大、可靠性分析难的问题。数值与工程算例验证了所提Kriging梯度解析解的较高精确性，同时验证了所提基于Kriging解析解的AFOSM结构可靠性分析方法的正确性与较高精度。     

动态时频谱分析、探测和跟踪的随机有限集法

动态出现和消失多分量信号的时频分析问题一直是非平稳信号处理的难点之一。为此，提出了一种分析、探测和跟踪多分量信号的随机有限集法。该算法利用时频变换，如短时傅里叶变换或自适应谱估计法，以及多项式预测模型，将多分量信号的时频分析问题归纳成可利用随机有限集进行多目标追踪的问题。分析表明：借助于提出的初始权重赋值算法，以及谱分量幅度和频率的联合似然函数，就可利用高斯混合概率假设密度滤波器来实现对动态时频谱的分析、探测和跟踪。在仿真实验中，所提算法有效提升了动态时频谱的跟踪精度，其对微弱时频谱分量的探测能力，以及对载频差异的分析能力均优于文献报道的算法。     

限带白噪声随机过程的雨流幅值概率密度函数模型

对17种不同带宽系数的限带白噪声随机过程进行了雨流循环计数统计,提出了一个限带白噪声随机过程的雨流幅值概率密度函数(Probability density function,PDF)模型。该模型是Rayleigh分布和Weibull分布的线性组合,其中待定系数均为随机过程谱参数的函数。采用该模型对这17种随机过程的雨流幅值概率密度函数进行了公式拟合,探究了模型中待定参数与随机过程谱参数之间的关系,确定了模型表达式。对照随机过程的雨流计数统计结果,将该模型与Dirlik模型的预测精度进行了比较,结果表明该模型的预测效果优于Dirlik模型。