基于双滤波器的高精度目标跟踪估计方法

针对连续推力的合作航天器,采用双重无迹卡尔曼滤波(DUKF)算法估计其状态和加速度.通过状态滤波器和参数滤波器的配合,提升滤波精度,完成运动状态和参数的估计,从而实现合作目标的运动轨迹跟踪.与合作航天器相比,非合作航天器存在大小未知、发生时刻未知的机动,无法获得加速度,且信息获取和运动状态的估计难度大.针对非合作航天器...     

星载GPS相位非差低轨卫星定轨及其精度分析

星载GPS相位非差低轨卫星事后精密定轨无需考虑复杂的动力学模型和地面资料,只需低轨卫星上的GPS观测资料和IGS的GPS精密星历产品,而且对于不同高度的卫星定轨都适用,计算简单、方便,能快速、高精度地确定轨道,同时还能确定部分动力学参数。本文在研究相位非差定轨方法基础上,对低轨卫星的误差影响及其处理措施进行探讨,给出了GPS相位非差定轨流程,编写了相应的定轨软件（SHKINE）,并利用CHAMP卫星资料对定轨的可靠性和精度进行分析,表明：利用自行编写的SHKINE定轨软件对CHAMP卫星定轨,3个方向坐标精度为10cm-20cm,点位精度为30cm-40cm,能满足一般定轨要求,是一种简单方便、行之有效的定轨方法。     

温度对磁屏蔽霍尔推力器磁场构型的影响研究

霍尔推力器磁路设计主要通过常温静态磁场仿真得到，并实测推力器非工作状态常温磁场进行复核。大功率霍尔推力器将面临更为严峻的热问题，推力器工作时磁路系统受高温影响，因此在常温下仿真得到的磁场位形会因温度升高而产生偏移，不能反映推力器真实工作时的磁场情况。为研究霍尔推力器工作时热量对磁路系统的影响，通过热磁耦合仿真对10kW磁屏蔽霍尔推力器的热态磁场分布进行研究，并对热态、常温仿真结果进行了对比，发现在阳极附近的径向磁感应强度Br的差异比放电室出口更大。常温设计的磁屏蔽构型在热态时偏离磁屏蔽，磁场和壁面最大不符合度达到13%，通过陶瓷出口型面修正后重新获得磁屏蔽效果，使最大不符合度降低到4.8%以下。合理热设计有助于降低热载荷，热仿真得到磁路系统最高温度低于500℃，低于0.78倍的居里温度Tc磁性急剧转变点，不会出现磁性能急剧下降，但热量对磁屏蔽霍尔推力器磁场构型的影响是应该考虑的。     

基于虚拟慢时间的双基地ISAR成像算法

针对双基地角时变引起的逆合孔径雷达（ISAR）图像畸变和散焦问题，提出了一种基于虚拟慢时间的成像算法。首先，分析了双基地角时变对ISAR成像的影响机理。然后，基于图像对比度最大准则估计等效旋转中心位置，完成初次相位补偿。最后，通过虚拟慢时间构建基于非均匀虚拟采样的补偿系数矩阵，并通过方位向非均匀傅里叶变换得到目标的ISAR像。算法基于图像对比度最大准则解决等效旋转中心位置估计问题，通过虚拟慢时间消除转动相位项的高次项影响，利用非均匀傅里叶变换解决随机虚拟采样的谱估计问题。理论分析和仿真结果验证了算法的有效性和鲁棒性。     

适用于非定常流模拟的分布式并行GMRES方法

为提高计算流体力学方法的收敛性和对高性能并行计算机的适应性,发展了适用于非定常流模拟的GMRES并行全隐式方法,并开展了相应的收敛和并行特性研究。采用变子空间数GMRES方法,减小重启过程计算时间;通过分区并行和Hybrid LU-SGS预处理算子实现方法的分布式并行化;采用鲁棒的Negative-SA湍流模型获得更大CFL数,采取计算和存储雅可比矩阵、网格重排序方法提高计算效率。利用这套方法完成了平面流、NACA0012翼型扰流、翼身组合体扰流、F-16战斗机非定常气动弹性和旋翼前飞流场的数值模拟。结果表明其计算效率较LU-SGS方法提高20%～200%;适用于当代高性能计算机分布式并行结构,并行效率非常高,在240个计算核心上出现了加速比的超线性。     

配置点谱方法求解非均匀介质内辐射传热问题

发动机内的燃气等高温介质随着组分和浓度的变化，会引起折射率在空间上的非均匀分布，从而导致辐射能束沿着曲线传播，其相应的辐射传热过程也更为复杂。为了避免射线追踪方法的复杂计算、提高计算效率，本文提出了配置点谱方法求解二维非均匀介质内辐射传热问题。在求解过程中，角向采用离散坐标法处理，空间采用配置点谱方法处理。通过将三种非均匀介质内辐射传热问题的配置点谱方法结果与文献结果进行对比分析，发现配置点谱方法可以在较少的节点数下，获得准确、有效地计算结果。并且，采用配置点谱方法求解三种算例的计算时间均消耗较少，均在20分钟以内。这将为进一步开展发动机复杂结构内高温燃气辐射快速仿真提供基础。     

跨声速涡轮叶栅两点优化目标函数研究

针对跨声速涡轮叶栅单点优化方法难以获得整体工况性能提升、多点优化方法难以确定合理目标函数形式的问题，提出了两点优化的方法。为了节约优化时间成本，优化过程采用EIF （Equivalent inviscid flow） 模型进行数值模拟，通过添加惩罚函数保证叶栅满足设计流量和负荷要求，并采用叶栅效率线性平均的目标函数形式进行评价。选择两组跨声速涡轮叶栅进行优化设计，并利用CFD方法分析叶型变化对流场马赫数、激波和损失产生的影响。结果显示，所提出的优化设计方法在保证设计工况性能的同时，能够提升叶栅整体工况性能。通过流场分析，揭示了激波结构变化对不同工况损失影响的定性规律。综合全文研究后，给出了一种适用于跨声速涡轮叶栅两点优化设计的目标函数形式。     

非结构网格体心梯度求解方法的精度分析

对于空间二阶精度的非结构网格求解器,无黏通量和黏性通量的计算都需要流场变量在单元体心处的梯度值。体心梯度求解的精度很大程度上决定了求解器整体的精度。本文推导了几种常用的梯度求解方法(包括格林高斯法,节点格林高斯法和最小二乘法)求解体心梯度的理论精度,并从数值方面对上述理论精度进行了验证。在数值测试过程中,通过以当前单元体心为基准进行坐标局部缩放的做法,保证了在非结构网格上做精度测试时,网格拓扑能够严格保持不变。理论推导和数值测试一致发现:使用保线性权插值到节点的节点格林高斯法、使用共面单元的最小二乘法和使用共点单元的最小二乘法,不管网格拓扑关系如何,都能保证梯度求解为一阶精度。而格林-高斯法和使用距离倒数权插值到节点的节点格林-高斯法在一般的网格中只有零阶精度,只有在特定的网格和权函数关系下,才能有一阶精度。最后用NACA 0012翼型亚声速绕流模拟进一步验证了梯度求解方法的收敛精度。