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

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

PMSM电流预测控制参数误差量化分析方法

交流永磁同步电机（PMSM）电流环控制性能是制约交流伺服系统性能的关键。电流预测控制拥有更快的动态响应、更低的电流谐波和优良的转矩响应，但该算法依赖精确的电机模型，参数失配会引发稳态电流误差，无法输出额定转矩，进而导致电机转矩输出效率降低。根据永磁同步电机电流预测模型，详细分析了d、q轴电流静差产生的原因，以及电流预测控制对电机参数误差的敏感性，提出了一种参数误差量化分析方法。该方法引入了电机参数偏差因子，量化描述了电机电感、磁链参数误差与电流静差之间的制约关系。通过仿真分析，验证了所提方法的合理性，为高性能永磁伺服电流预测数字控制技术打下了良好基础。     

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

针对连续推力的合作航天器，采用双重无迹卡尔曼滤波(DUKF)算法估计其状态和加速度。通过状态滤波器和参数滤波器的配合，提升滤波精度，完成运动状态和参数的估计，从而实现合作目标的运动轨迹跟踪。与合作航天器相比，非合作航天器存在大小未知、发生时刻未知的机动，无法获得加速度，且信息获取和运动状态的估计难度大。针对非合作航天器，基于简化的相对运动方程，结合天基平台获得目标的观测信息，采用两个扩展卡尔曼滤波(EKF)及基于半正焦弦的机动检测策略实现多未知脉冲机动的运动状态的估计。仿真结果表明：相比于无迹卡尔曼滤波（UKF），DUKF在对合作航天器的状态和加速度估计方面具有更快的收敛速度和更高的滤波精度；对于存在未知机动的非合作航天器，通过对比验证机动检测策略与滤波器切换策略相结合的方法的有效性，该方法能够检测到多次机动并且减少误判。     

速率模式飞轮姿态控制系统飞轮组合平稳切换方法

对速率模式飞轮姿态控制系统的能控性、稳定性和飞轮组合切换过程进行了研究，揭示了在飞轮组合切换过程中引起卫星姿态抖动的原因，提出了通过对残余角动量进行卸载并将卸载函数同时前馈到系统输入的方法，以保持卫星姿态的稳定度。从理论上证明了可引入残余角动量进行卸载，将卸载函数同时前馈到系统输入，并且不改变原系统的能控性和稳定性。给出了某卫星飞轮切换应用实例，说明了飞轮组合平稳切换方法的有效性。     

脉冲星方位误差估计的两步卡尔曼滤波算法

为了克服钟差和卫星位置误差对脉冲星方位误差估计的影响，设计了两步卡尔曼滤波（TSKF）算法。首先，介绍了脉冲星方位误差估计的传统模型，并通过分析和仿真验证了钟差、卫星位置误差以及2种误差同时存在时会使脉冲星方位误差估计结果产生较大偏差。其次，在传统的估计模型中加入了钟差和卫星位置误差，并将钟差和钟差变化率增广为新的状态量，从而推导出包含2种误差的新模型，并证明了该模型的完全可观测性；根据该模型并按照两步卡尔曼滤波原理，得到了TSKF算法的步骤。最后，通过仿真表明:在钟差和卫星位置误差同时影响下，传统脉冲星方位误差估计算法偏差较大且发散；TSKF算法则能够有效隔离2种误差的影响，使赤经和赤纬误差估计达到0.2 mas之内的精度。      

MIMO窄带加宽带非高斯随机振动试验

提出了一种用于多输入多输出窄带加宽带非高斯随机振动试验方法。首先，详细分析了随机相位与非高斯随机信号偏度和峭度的关系，提出了一种可用于快速生成具有指定偏度和峭度的迭代相位调节法。其次，根据控制通道之间的耦合关系，利用各通道之间随机相位的不相关特性，实现了多通道非高斯随机信号生成的解耦。然后，将生成的非高斯随机信号作为参考输出响应，采用时域逆系统方法生成驱动信号，并利用控制算法分别对响应信号的功率谱、偏度和峭度进行控制。最后，通过三轴向振动台试验验证了本文提出的多输入多输出窄带加宽带非高斯随机振动试验方法的可行性。结果表明，控制点上加速度的功率谱被控制在了±3 dB的容差限内，满足预设的容差要求。同时，时域响应信号的峭度和偏度也满足给定的参考要求。     

基于COSMIC掩星精密定轨数据的等离子体层电子含量研究

等离子体层是日地环境重要的组成部分.本文利用COSMIC掩星精密定轨数据经处理后得到的podTec文件获取等离子体层电子含量（PEC）对等离子体层进行研究.将podTec数据进行处理后获得的PEC（pod-PEC）和IRI-Plas经验模型提供的PEC （IRI-PEC）进行对比，发现pod-PEC与IRI-PEC符合得较好.在低（0°—20°）、中（20°—50°）、高（50°—90°）修正地磁纬度带下，分析了COSMIC在太阳活动极大年（2014年）3，6，9和12月的pod-PEC，得到如下结论:PEC随着纬度升高而逐渐减少，且3，9月PEC在中低纬关于磁赤道的南北对称性较好，6月北半球各纬度带的PEC明显高于南半球同一纬度带的值，而12月情况则完全相反，南半球中纬的PEC甚至会等于北半球低纬的PEC值；PEC在白天高而晚上低，高纬地区的PEC昼夜变化不明显；PEC具有明显的季节性.对于北半球，一年中PEC最大值出现在春季，冬秋季次之，夏季最低，具有明显的年度异常现象.      

SABRE预冷器计算模型及其在整机模型中的应用

预冷器的性能对SABRE等预冷组合循环发动机具有重要影响,为实现发动机方案设计阶段预冷器的快速设计与评估,建立了预冷器准二维快速评估模型。将SABRE预冷器的几何结构简化为一个扇环形区域,沿径向和周向将该区域划分为二维节点。应用守恒方程及传热关联式完成单个节点计算,再求解节点矩阵的平衡方程组,计算内外流体特定节点上参数的二维分布,得到预冷器出口参数。将模型嵌入发动机总体性能程序中,实现了发动机设计及非设计点的预冷器性能计算的功能。与文献数据对比结果表明,预冷器模型传热计算误差小于5%,摩擦阻力误差小于10%。整机计算结果显示,Ma0~5范围内,预冷器空气侧温降范围为143K至932K,温降随飞行马赫数升高单调上升。预冷器传热有效度范围为0.896-0.945,空气侧总压恢复系数范围为0.852-0.904。     

基于光学测角数据的风云四号同步轨道卫星精密定轨

针对风云四号同步卫星的精密定轨和精度评估需求，首先利用地面光学测角数据对FY-4A卫星进行精密定轨，定轨后方位角和高度角的残差rms分别为0.25＂和0.45＂。与基于测距数据的轨道相比，位置精度在有测角数据的弧段内小于50m。进一步联合测角数据和测距数据对FY-4A卫星进行联合定轨，定轨后轨道重叠精度优于15m。利用联合定轨结果评估了基于测距数据的实时轨道产品精度，可以明显发现轨道精度随着测距数据的积累而逐步提高。     

一种调焦机构运动方向与光轴平行性测试方法

调焦机构运动方向与光轴平行性的好坏直接影响调焦过程中光学系统主点位置、像质均匀性、像面照度均匀性等多项指标。为了在安装过程中保证调焦机构运动方向与光轴平行,需要测量调焦机构运动方向与光轴之间的角度关系。文章通过三坐标测量系统和自准直经纬仪,设计了一种对感光元件非接触式的测量方法,来测量调焦机构运行方向与焦平面法线方向之间的角度关系,通过结构设计与加工补偿实现两者的平行装调。试验结果表明,调焦机构运动方向与焦平面法线方向的角度测量模型准确可靠,数据处理算法合理可行。