可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展。但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题。界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因。为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案,本文从数值框架的选择、非守恒方程相容离散、高精度有界格式构造、界面压缩、界面-单介质分区计算方法等多个维度出发,综述近几年我们在可压缩多介质流动问题高精度数值模拟方法方面的研究进展。通过上述多个维度的工作,我们建立了一套适用于可压缩多介质流动问题的低耗散、基本无数值振荡的高精度欧拉数值方法,并成功应用于可压缩多介质大变形流动和界面不稳定性诱导湍流混合等问题的数值模拟。相关数值方法研究成果已集成至武器物理等领域工程数值模拟软件中,为相关工程任务提供了重要技术支撑。     

AIRS观测资料研究全球平流层重力波特性

利用2012—2014年1月和7月AIRS（Atmospheric Infrared Sounder）第79通道的观测数据，分析了平流层重力波活动强弱的全球分布以及重力波发生频率的全球分布；分析了重力波活动随纬度和经度的变化特征，给出了重力波活动在全球范围内的热点区域及其活动强度；对比了白天与夜间的重力波活动强度及发生频率.研究表明重力波活动强度呈现出随纬度变化的特征，在低纬度地区（0°—30°），冬季半球重力波活动强度低，夏季半球重力波活动高；在中高纬度地区，冬季半球重力波活动强度高而夏季半球重力波活动强度低.在1月，全球重力波活动有4个突出的热点区域，分别为50°N附近欧洲大陆与大西洋交接地带、北美洲与大西洋交接地，20°S附近南美洲与大西洋交接地区、非洲与印度洋交接地区.在7月，重力波活动突出的地方为巴塔哥尼亚至南极半岛地区，以及50°S和75°E附近的印度洋区域.重力波活动强度在夜间大于白天，但是夜间的强重力波活动区域小于白天.      

考虑燃料消耗均衡性的低轨通信星座在轨星座构型重构方法研究

针对低轨（low earth orbit,LEO）通信星座在轨重构问题，分析了低轨通信星座的特点，提出了4个重构指标，即全球平均覆盖率、燃料消耗均衡性、重构总时间和重构总速度增量，并分析了提高轨道高度的重构方式。采用基于分解的多目标进化算法（multiobjective evolutionary algorithm based on decomposition，MOEA/D）对重构进行建模分析。结果表明，模型优化结果可以得出多组最优解，构成帕累托前沿，并得出燃料消耗均衡性最优的解，针对两种失效模式得出的重构总时间均为3.6×105s，速度增量方差分别为261.4m2/s2和293.4m2/s2，这些解均可恢复星座原有的覆盖性能，并使燃料消耗最为均衡。     

事件触发机制下的充液航天器姿态控制

针对液体大幅晃动、通信资源受限的充液航天器姿态控制系统，提出一种自适应滑模控制与事件触发机制相结合的控制策略。首先，针对固-液耦合的充液航天器姿态控制系统，选用滑模变结构控制来削弱液体大幅晃动的非线性影响，并设计自适应更新律在线估计不确定参数来提高系统的鲁棒性。然后，考虑星载计算机资源的限制，设计相对阈值的事件触发机制来决定控制输入信号的更新，从而减少控制器与执行器之间的信号更新对通信网络的占用。最后，仿真结果表明，在液体大幅晃动下，所提控制策略不但可以使航天器姿态控制系统最终收敛到任意小的界内，而且可以减少96%的控制信号传输，减轻航天器的通信负载。      

军航无人机与民航航班侧向碰撞风险评估

近年来无人机(UAV)已经成为新型作战力量。为解决军航UAV与民航航班之间的飞行冲突，确保飞行安全，对经典EVENT模型进行改进，提出适用于UAV的碰撞风险模型。着重研究了导航方式、人为因素、高空风对于UAV飞行的影响，并构建了相应的位置偏差模型。首先，利用蒙特卡罗法求解军航UAV与民航航班之间侧向间隔丢失的频率；其次，利用MATLAB进行算例仿真，验证模型的有效性，并得出了侧向碰撞概率随参数变化的关系；最后，通过计算不同安全间隔下的碰撞风险，对空域使用提出建议。      

SBAS电文时序动态编排算法

星基增强系统（SBAS）通过GEO卫星转发SBAS电文实现对GNSS服务性能的提升，以满足民航用户不同飞行阶段的导航需求，因此，合理有效的电文内容及播发时序设计是系统实现高质量服务的重要保证。为提高电文编排的灵活性，避免固定时序填补空余电文引起的播发资源浪费，提出了一种SBAS电文时序动态编排算法，在满足国际标准要求的前提下，综合利用SBAS电文龄期和最大播发间隔实现待播发电文的自动选择。利用NTMF实测数据对当前各主要SBAS的电文进行了特性分析，对所提方法的单双频SBAS电文编排效果进行了评估。结果表明:所提算法可保证电文时序符合国际标准要求，实现了重要电文的优先播发，将空余时隙进行动态分配实现了各类型电文播发间隔的近等比例缩短。与固定时序相比，单频SBAS完好性电文播发间隔缩短约15.0%，首次定位时间缩短约8%，双频SBAS电文首次定位时间缩短约6.5%；与固定时序的BDSBAS B1C电文相比，完好性服务能力提升约14.7%，首次定位时间缩短约16.7%。所提算法有效提升了SBAS电文播发的播发效率，实现了SBAS播发资源的100%有效利用。      

半球谐振子疲劳寿命分析与仿真

半球谐振子的寿命和损伤是直接影响高精度半球谐振陀螺使用时长和安全性的重要因素。目前国内加工的半球谐振子所用的熔融石英玻璃材料主要依靠进口，采用传统的疲劳寿命实验确定方法成本过于昂贵，因此需要利用疲劳分析软件对半球谐振子的疲劳寿命进行分析。文章通过半球谐振陀螺应力分析，采用ANSYS软件对熔融石英半球谐振子进行应力分析仿真，确定因残余应力所引起的疲劳危险部位，并在疲劳部位进行裂纹扩展分析，得到 、 和 型应力强度因子，以使在半球谐振子结构设计和使用过程中对易疲劳的部位进行有效监控和预防。     

逆金字塔式安全着陆点选择方法

针对现有地外天体软着陆安全着陆点选择方法的不足,设计了逆金字塔式安全着陆点选择方法,结合着陆器的落点精度和自主避障及地形坡度适应能力,建立了基于平整点、安全着陆点、可靠避障点和标称着陆点的地形安全性分析与定量搜索模型,该模型呈四层逆金字塔结构,每一层搜索的结果作为下一层搜索的输入,从弱约束到强约束逐步收敛至100%地形安全着陆区,嫦娥五号探测器飞行结果表明该方法能够快速确定安全着陆区及标称着陆点,并有效指导着陆下降过程的标称航迹设计和轨道控制策略制定。     

基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。     

激光捷联惯导量化热设计指标控制方法

随着机载激光捷联惯导向小型化方向发展,惯导结构的热设计正面临着更加严峻的挑战.针对设计中经常出现因无法在设计前期进行热设计指标控制而导致后期设计返工的问题,提出了一种热设计量化控制方法,即根据估算的整机功率计算箱体外表面积,再对箱体外表面积可达性进行判断,如可达则以该参数作为控制指标进行设计,如不可达则采用可达的最大表面积作为输入反算整机功耗,再根据计算结果对各电路功耗进行重新分配.经有限元仿真计算和试验一致证明该方法有效,可以实现惯导设计前期定量化的热设计控制,在很大程度上避免了后期因热可靠性及热试验问题的设计返工.