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

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

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

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

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

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

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

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

空间目标碰撞概率的等效矩形域快速算法

针对线性相对运动假设下的碰撞概率计算问题，提出了等效矩形域方法。利用该方法将概率积分计算进行近似处理，推导出概率密度积分的解析表达式。针对空间目标误差椭球形状固定的情况，给出了最大碰撞概率的计算方法。对目标在交会时刻相距较远和接近碰撞的情况进行了碰撞概率以及最大碰撞概率计算。通过将计算结果与空间压缩无穷级数法的对比，验证了等效矩形域方法近似计算碰撞概率积分的可靠性，同时对于不同的碰撞概率计算情况，等效矩形域方法显示出更高精确度和更好的估计偏差稳定性。     

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

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

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

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

多支点柔性转子系统临界转速稳健设计方法

基于响应面法和容差模型提出了定量考虑参数变差影响的转子系统临界转速稳健设计方法。该方法利用响应面模型获得多参数、多目标之间的函数关系,并以多阶临界转速分布为约束条件、以临界转速对支承刚度变差敏感度最低为设计目标。对小涵道比涡扇发动机低压转子系统临界转速稳健设计结果表明:在考虑各支点支承刚度变差情况下,采用多参数、多目标稳健设计方法,可保证多支点柔性转子系统的临界转速特性在满足设计要求的同时,转子系统临界转速对支承刚度变差的敏感度最低,验证了该方法的有效性。      

强电磁脉冲下柴油发动机系统薄弱环节识别

为识别强电磁脉冲环境下柴油发动机系统的薄弱环节，提出了一种加权故障树和分层贝叶斯网络相结合的柴油发动机系统薄弱环节识别方法。该方法综合考虑同层单元失效的相关性，加权故障树的局部应用解决了部分条件转移概率表不易获取问题。运用贝叶斯网络双向推理功能，首先，通过柴油发动机辐照试验和电磁仿真软件获得的各部件敏感度阈值及电磁应力数据，计算出强电磁脉冲下部件级到系统级的先验失效概率；然后，依据贝叶斯概率公式计算在发动机失效条件下各部件故障的后验概率，并排序以识别其薄弱环节，为电磁防护方案的设计提供参考和建议。以宽带高功率微波(WBHPM)辐照为例，说明了柴油发动机系统分层贝叶斯网络故障模型参数获取与概率计算过程。结果表明:执行器和凸轮轴位置传感器、曲轴位置传感器既为柴油发动机系统的重要部件，也为较薄弱环节，是需要重点防护的对象。