基于BVD原理的高保真空间重构方法

简要综述了一类基于单元边界变差最小化(Boundary Variation Diminishing,BVD)原理,设计双曲守恒律高保真数值格式的空间重构方法。BVD原理要求尽量减少通过重构得到的网格边界两侧物理量之间的差,从而能够有效地控制黎曼求解器中的数值黏性。BVD方法针对数值解的空间分布特征,选择多个函数作为空间重构的候补函数,并根据BVD判定准则从候补函数中选取最合适的函数进行空间重构。BVD判据不需要根据求解对象进行经验参数(阈值)的调整。选用适当的候补函数和BVD准则,可以完全避免现有算法中为抑制数值振荡而必须采用的非线性限制。BVD格式能在抑制数值振荡的同时,有效地控制数值耗散,可以对光滑解与间断解都获得高保真的计算结果。本文概述了BVD方法的基本思想、设计相关格式的基本思路,以及一些具有很强实用价值的BVD格式。并通过单相和两相可压缩流动的一些典型算例验证BVD格式的特点和优势。     

三元二次样条函数及其计算（Ⅰ）——高维数据拟合

半节点二次样条是Ｃ１类连续的，其二阶导函数是阶梯函数，在半节点处产生跳跃。鉴于此，本文利用最小二乘法，获得了一种半节点二次插值样条边界条件确定方法，该方法可以保证二次插值样条在半节点处二阶导数的变化最小，这相当于保证了曲率的变化最小。为了适应四维数据插值的需要，给出了三元二次插值样条的定义及基样条表示，提出了一种边界条件，证明了其存在唯一性。最后把二次插值样条边界条件确定方法推广到三元二次插值样条上去，获得了三元二次插值样条边界条件的确定方法。该方法的意义是，可以直接从插值条件获得边界条件，从而克服了多元样条边界条件难以确定的困难     

基于MILP的多无人机对敌防空火力压制

建立了基于混合整数线性规划（Mixed Integer LinearProgram,MILP）的多无人机编队对敌防空火力压制协同任务分配模型,以0-1决策变量表征无人机一任务指派关系,引入连续时间决策变量来表示任务的执行时间,并通过对决策变量之间的线性等式和不等式的数学描述,建立无人机之间和无人机执行任务之间合理的协同约束关系。采用商用软件CPLEX对模型求解,仿真验证了模型的合理性。     

基于微分进化算法的深空小推力轨道优化设计

针对多目标多任务的深空探测轨道优化设计问题,基于以太阳为主引力体的二体模型,采用电推进小推力系统,提出了将目标行星的探测方式、探测顺序、发射窗口同时进行全局优化的解决方案.依此建立小推力轨道优化设计算法模型,利用微分进化算法进行全局一体化优化设计,得出多目标多任务的一体化深空探测轨道.通过与相同条件下其他设计方法的结果比对分析得出：基于小推力推进方式并采用微分进化算法进行优化的轨道优化设计方案,在多目标多任务深空轨道一体化设计问题上具有可行性及应用价值.     

面向海洋机动目标跟踪监视的卫星重构组网优化与仿真

摘要： 针对海洋机动目标跟踪监视问题，提出通过卫星相位重构调整进行机动目标搜索跟踪的监视策略，构建地理坐标系下机动目标航位预测模型及基于霍曼转移的卫星调相组网方式，并采用NSGA II算法进行优化求解，对比分析卫星重构组网前后的机动目标发现概率.仿真结果表明：在卫星能量消耗允许范围内进行重构组网，能够有效提升对机动目标的发现概率，为海洋机动目标的搜索跟踪问题提供了一定的方法支持.     

边界层内WENO格式特征重构精度的理论分析

WENO格式因其高精度以及良好的激波捕捉能力而被广泛地用于可压缩湍流的数值模拟过程中。为了在有强间断的流场中仍然能保证数值稳定，通常将流动物理量投影到特征空间后再进行WENO重构。然而这种特征重构方法在边界层内的精度问题仍需要进一步分析与研究。通过可压缩边界层自相似解构造平板边界层内法向的流动参数，并在该方向上进行WENO重构，得到求解半点数值通量的WENO非线性权重。通过分析实际得到的权重与理想权重的偏差大小发现在边界层流动中特征重构的精度明显低于分量重构。为了更深入地分析特征重构精度降低的原因，通过对曲线坐标系下的流动方程进行理论分析，推导得到特征重构时的特征变量形式。根据特征变量表达式以及自相似解中权重误差分布发现引起特征重构精度降低的主要原因是特征投影过程在半点处出现了额外的极值点，同时这种因特征投影过程而产生极值的特性将不会随着左右特征矩阵的不同选取而发生改变。进一步地，基于理论分析的普适性，对于任何光滑流场，只要使用特征重构的投影过程，都会在半点处产生极值从而使得WENO系列格式产生精度降低的现象。     

可压缩平板边界层转捩的大涡模拟

入口利用层流速度型叠加二维和三维Tollmien-Schliting(T-S)波,应用动态亚格子模型、预处理技术以及高精度算法,采用大涡模拟(Large eddy simulation,LES)数值计算方法,对可压缩平板边界层空间模式的转捩机理及其拟序结构进行了研究。计算结果展示了T-S波经历线性、非线性发展阶段破碎形成发卡涡,最终演化成具有4个涡环的涡链结构的整个过程,并揭示了发卡涡涡腿并非完整对称及低速条带在涡环产生的位置不连续的现象,探讨了主流向涡与壁面相互作用诱导生成次流向涡的机理。     

基于引射器函数法的程序优化设计及验证

针对考虑多因素综合影响的引射器优化问题,基于引射器函数法,深入研究了静压协调函数的数学特性并分析了其数学曲线上的奇异点。在此基础上,通过程序设计和基准推进归纳法得到了多分支工作特性曲线,分析了不同分支下解的特性。另外结合工作特性曲线重点研究了混合室背压、主次流总压比和混合不均匀度对引射性能的影响,提出了静压特性曲面、临界曲线的设计概念。通过试验数据对比分析验证了引射器函数法的可靠性。结论表明:靠近临界曲线工况时引射器设计性能较好;混合不均匀度对引射性能影响重大,不均匀度为1.5时对比理想状态,引射系数最大误差达到32.73％;考虑壁面摩擦设计时需对摩擦因数公式模型进行修正。研究结果为引射器优化设计提供了重要指导。      

一种基于故障统计数据的发动机风险预测

给出了发动机故障危险等级划分 、危险系数确定、风险因子与飞行风险的计算方法,建立了基于蒙特卡罗仿真的发动机多故 障模式风险评估流程,预测发动机故障发生情况,评估发动机使用阶段的故障风险。在已有 航空发动机历史故障数据的基础上,对多种故障模式进行了风险预测,制定了两种可降低风 险的方案,分析了维修周期对风险因子的影响,给出了合理的降低风险方案。     

基于行为的非合作目标多航天器编队轨迹规划

目前航天器编队轨迹规划的研究并未综合考虑目标非合作性、构型复杂性及控制协同性等问题,针对非合作目标多航天器编队系统,提出了一种基于行为的相对运动轨迹规划方法。首先,基于聚集行为模型设计了期望运动场,将期望速度视为一系列具有不同行为特性的速度矢量和,通过构造平衡态构型公式计算行为调节参数。其次,采用高斯型环境函数描述了非合作目标特性,并通过C-W方程对期望运动场进行了改进,使得该方法在多航天器编队系统与非合作目标形成期望构型的过程中,既能避免编队成员之间发生碰撞,又能保证与非合作目标之间的安全距离。仿真结果表明,该算法可以使多航天器编队系统在3000s内运动至期望构型。该方法在复杂的、多个控制目标的编队飞行轨迹规划中优势明显,具有自主性、协同性、鲁棒性强等特点,同时可应用于不同类型的编队构型相对运动规划。