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.     

计算机帧时间对半实物仿真的影响

半实物仿真试验中,要求计算机帧时间与计算积分步长相一致,但在实际试验中不容易精确地做到。本文分析了这种不一致给半实物仿真试验带来的影响,给出了解决的办法,文章中还给出了一实例来说明。     

均布压力作用各向异性矩形薄板大挠度问题

基于冯·卡门模型，选择中心挠度为摄动参数，利用摄动方法将正交异性矩形薄板大挠问题的非线性偏微分方程组逐级线性化，导出了各级摄动对应的几个线性偏微分方程组，然后再借助变分法求解，得到了均布载荷作用下正交异性矩形薄板的载荷-挠度曲线。     

轮控小卫星姿态大角度机动递阶饱和控制器设计

针对采用反作用飞轮的小卫星姿态大角度机动控制，在反作用轮输出力矩受限、速率饱和的约束条件下，采用递阶饱和方法，即限制卫星每次姿态机动的最大偏差，对姿态偏差进行逐次消除。在毋需获知最优机动轨迹规划的情况下，可用于卫星任意时刻的姿态捕获和机动控制。数学仿真结果表明，本文设计的控制算．去能够实现快速姿态机动任务，具有良好的鲁棒性。     

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

室火轰燃的大涡模拟

采用一套修正的适用于低马赫数浮力流的N-S方程组来描述火灾场,运用Smagorinsky亚格子模型的大涡模拟方法对室内池火的轰燃发展过程进行了数值模拟,把模拟与实验结果进行了对比和分析。结果表明,该方法能较好地预测轰燃。同时模拟结果显示火焰区向房间后部偏转稍大,近壁处温度预测的误差大一些,随着亚格子模型的不断完善和计算能力的提高,用该技术来预测轰燃现象会更加准确。     

地面模拟静止微重力环境中导线先期着火特性研究

针对在微重力环境中运行的载人航天飞行器上的电缆和导线在工作时由于电流过载导致温度升高而引起着火的情况,提出了"功能模拟"实验原理,并且利用地面实验设备对微重力环境下导线的着火前期特性进行了功能模拟实验研究.通过实验,得到了在微重力情况下由于浮升力的减小使自然对流减弱导致电流过载时导线的热平衡温度高于地面正常重力情况,从而证明了这正是引起航天飞行器着火的潜在点火源.      

超声速凹槽大涡数值模拟研究

在超声速燃烧室内常采用凹槽增强燃料与空气的混合,实现火焰稳定,本文采用大涡模拟方法对开式凹槽流场进行了数值模拟研究,计算采用Smagorinsky 亚网格模型模拟小尺度涡的作用.数值模拟结果给出了凹槽流动的自激振荡的发展过程,以及凹槽流动中拟序结构的演变过程,如涡的卷起、增长,涡与涡之间的合并、破碎等,计算得到的凹槽波动的峰值频率和理论结果一致.     

超声速混合层液滴两相流动的大涡模拟

为研究超声速气流中液滴与气流的混合及液滴蒸发对混合的影响,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构,气相流场采用亚格子(SGS)模型和切应力输运(k-ωSST)湍流模型,液相模拟采用轨道模型和单液滴蒸发模型。在混合层前缘入口处均匀持续地投放液滴,并在液滴入口处下方添加非周期小扰动,并观察液滴蒸发过程对该小扰动造成的影响。分析了入口小扰动在流场中不同的发展情况,发现液滴的蒸发过程使混合层厚度增加并加速混合层的发展,对气相流场扰动较强,可能导致流动失稳,对混合过程有很大的促进作用。     

基于DES方法的三角翼激波-涡干扰流场数值模拟

采用基于Spalart-Allmaras湍流模型的脱体涡模拟（DES）方法,数值求解Navier-Stokes方程,模拟绕尖前缘三角翼的跨音速流动,并对三角翼上翼面的复杂激波-旋涡干扰流场进行了分析。与NASA兰利研究中心的NTF风洞实验结果对比分析表明,DES方法能很好地模拟跨音速三角翼上的旋涡流动。随着攻角由中度攻角增加到大攻角,支架附近的激波越来越强,对主分离涡的干扰作用越来越大,直至出现激波干扰导致的涡破裂。激波的形状、位置及涡破裂位置均与实验结果吻合良好。