增强高超声速化学反应流数值计算的稳定性方法

针对高超声速化学反应流数值模拟中包含源项对角化雅可比矩阵的数值格式稳定性很差的问题,在对采用上下对称的高斯-赛德尔(LU-SGS)格式的离散方程进行线性稳定性分析的基础上,提出一种依据流动状态对流场的块网格进行局部加密的方法,以增强流场数值计算稳定性。在高超声速流场中,对包含激波后亚声速高温区域的块网格给予加密,而其他区域中的块网格可以适当进行网格粗化或保持不变。算例结果表明,相比于初始网格,局部加密后的改进网格并不会产生过多的网格量,而且能够采用较大的柯朗-弗里德里奇-列维（CFL）数以达到数值求解过程中稳定计算和快速收敛的目的。     

应用动网格技术计算着陆过程中缓冲发动机喷流流场

利用Fluent软件对着陆过程中缓冲发动机的喷流流场进行数值计算,发动机与地面相对位置的不断变化使计算域形状不断改变.为了保证计算的正确性,必须采用动网格技术.介绍了实现移动网格的基本方法,并结合使用网格弹性平滑移动和网格重构方法实现了着陆过程中缓冲发动机喷流流场的非定常数值计算;得到了随发动机着陆产生的喷流流场的压强...     

某炮射导弹发射过程的阻力系数及升力系数的数值模拟

用数值模拟的方法求解炮射导弹发射过程的主要空气动力系数，以某炮射导弹的气动外形为基础，以优秀的计算流体力学软件FLUENT、为工具。在数值模拟中全部计算网格采用结构化网格，紊流模型采用适于空气动力计算的一方程模型中SPALART-ALLMARAS模型，差分方程的求解采用一阶迎风格式和二阶迎风格式相配合的方式，所有这些保证在数值模拟中的物理模型、数学模型、以及方程的求解等方面提高数值模拟结果的精度，最后给出了主要的空气动力系数，并对其进行了简单的分析。     

应用TTM网格研究冲压增程弹丸进气道内外流场

基于TTM贴体网格生成法,研究了某冲压增程弹丸超音速进气道内外流场数值模拟中单区域网格生成问题,流场求解采用了二阶隐式TVD格式。NACA0012翼型跨音速流场算例结果表明,所采用的算法和编制的程序可靠,能够适用于具有亚、跨、超音速三种流态共存的复杂流场数值计算。对超音速进气道数值模拟结果显示,得到的流场复杂波系结构是合理的,来流马赫数和攻角变化显著影响着进气道性能。     

三维横向流喷与超声速主流干扰流场的数值模拟

本文对三维隅角机翼声速喷流与超声速主流的干扰流扬进行了数值模拟。三维欧拉方程的求解采用非结构网络有限体积伽辽金法（Finite Volume Galerkin Method）。引入了总体结点积分域的概念,简化了从单元矩阵到总体矩阵的汇总过程。通量的分裂采用Osher格式,通过外差使其由一阶精度上升为二阶精度。发展了一种基于线化流量的逆风非结构网格隐式有限元格式以提高求解精度及效率。最后给出了三维隅角机翼流场的算例。     

高超声速飞行器外形热流密度分布计算的高精度方法研究

为提高高超声速条件下飞行器表面热环境的计算精度,发展了高精度算法WCNS-E-5,并开展了热流密度分布的数值计算研究.该算法采用的是五阶精度加权显式非线性格式,结合四阶精度的二阶偏导数差分近似、四阶精度的边界格式,同时还对网格偏导数及温度梯度也离散为四阶精度.通过求解三维Navier-Stokes方程,数值研究了高超声速飞行器外形如钝锥和双椭球体的物面热流密度分布.模拟结果表明,这种保证全流场高精度的WCNS-E-5得到的流场图像清晰、真实.与通常的二阶算法MUSCL对比,WCNS-E-5具有更高的流线分辨率,获得的热流密度更接近于实验测得分布.     

用空间推进算法模拟高超声速进气道流场

在有限体积法框架下，采用空间推进算法SSPNS（Single-Sweep Parabolized Navier-Stokes Algorithm）求解抛物化NS方程（即PNS方程），在流向采用LU—SGS隐式积分，而横向无粘和粘性通量则分别采用AUSM系列格式和中心格式计算。用该方法对1个二维高超声速进气道和2个三维高超声速进气道流场进行了数值模拟，得到的流场波系结构、壁面压力及传热系数分布与文献中相关数值解和实验数据基本一致，表明SSPNS法能够准确地模拟超燃冲发动机进气道内的高超声速流动。对比研究表明，SSPNS法与求解FNS（Full Navier＼lStokes Equations）方程的传统时间迭代法相比，二者计算精度相当，而SSPNS计算速度快1～2个量级，存储量至少低1个量级。本文的研究为CFD在超燃冲压发动机部件及一体化优化设计中的集成，以及大型高超声速工程流动的高效计算，打下了良好的基础。     

三维烧蚀内部热响应数值计算研究

对热解型碳化复合材料三维烧蚀内部热响应数值计算关键技术进行了研究。采用碳化层—热解面—原始材料层模型,将热解气体与碳化层之间的对流换热处理为源项,通过有限元法建立移动边界条件下温度场求解方程组,采用Gauss Seidel迭代法计算热解气体质量流量和温度场。同时,研究和分析了三维烧蚀移动边界处理方法以及动网格生成方法。由于每个时间步都需要网格重划,烧蚀热防护数值计算对存储效率和计算效率要求较高,本文研究了内部热响应计算中影响存储效率和计算效率的主要因素,并提出了相应的压缩存储方案和求解方案。计算结果表明,移动边界处理方法准确合理;存储方案的存储效率较高;保持刚度矩阵和形函数矩阵正定对称性可以加快温度场计算的收敛速度。     

残量分布格式初步研究与应用

介绍了一种新的流场计算数值格式——残量分布格式。它是一种多维迎风格式,适于三角形网格和四面体网格。给出了该格式的标量方程和Euler方程。算例显示:该格式能处理稳态与非稳态问题,计算精度高,有较强的捕捉激波能力,且结构紧凑,易实现编码。     

发动机喷管内流场对流换热系数影响因素的数值分析

基于Roe-FDS格式,采用标准κ-e湍流模型对某液氧煤油发动机喷管内部燃气流对喷管内壁的对流传热换热系数进行了数值仿真与研究,分析了网格雷诺数以及恒温壁壁温对喷管内壁对流换热系数的影响.当近壁面网格第一层法向高度为10-3 mm数量级时,仿真结果与工程估算公式的计算结果相吻合.同时,在一定范围内,对流换热系数随来流速...     

分析插值误差和斜率的轨迹优化网格细化方法

提出一种基于插值误差和斜率分析的轨迹优化自适应网格细化方法，包括节点插入算法和节点删除算法。节点插入算法分析各个离散节点的控制变量的插值误差。若插值误差较大，则在该节点周围增加节点细化网格；否则，不进行细化。节点删除算法分析各个离散节点处的控制变量斜率。若某个节点的左斜率和右斜率都为零，那么删除该节点；否则，保留该节点。采用三个典型的轨迹优化算例验证了所提出的方法的有效性和特色，并且与其它几种网格细化方法进行了对比。仿真结果表明，本文方法生成的网格规模较小，需要的网格迭代次数较少，能够快速、高精度求解非光滑轨迹优化问题。     

一种求解轨迹优化问题的改进多分辨率技术

针对轨迹优化多分辨率技术可能发生细化遗漏或者细化失败的缺陷，设计了一种改进多分辨率技术。该方法通过在每次网格细化迭代过程中增加两层高分辨率节点并且引入节点检测算法，确保在轨迹的非光滑区域能够高效、持续地细化网格。针对原始多分辨率技术的初始网格的节点数目存在二进限制的弊端，定义了一种新的广义二分网格，使得多分辨率技术的初始网格节点数目可以为任意奇数。采用多个轨迹优化算例验证了本文方法的有效性并且与其它方法进行了对比。此外，仿真结果还表明，网格充分细化对目标函数影响不大，但是能够更加准确地满足路径约束和终端状态约束。     

一种基于图割的稠密三维场景重建算法

立体视觉技术是自主探测机器人在未知环境中获取信息的重要手段,通过对可视场景的稠密三维重建实现导航、定位及路径规划等一系列工作。本文在候选点匹配的基础上结合图割理论,首先在世界坐标系建立代表深度信息的网格节点,接着依据区域匹配算法对候选点进行初步筛选,去除大部分相关值较低的节点,建立简化的网格图,最后通过寻找图中最小割来实现能量函数的全局最小,完成稠密的三维场景重建。实验证明,相关阈值γ设为0.6时,简化网格图的重建精度和计算效率达到相对平衡。图割算法解决了候选点测量时潜在的匹配歧义问题,且对低纹理区域有较好的匹配效果。     

卫星动平台下动态门控调度的设计与实现

当运动节点构成的大规模动态平台执行分布式协作任务时,如何保证关键数据通信的实时性和确定性,这是一项新的挑战。文章以时间敏感网络（time sensitive networking,TSN）技术为背景,结合卫星构成的分布式动态通信平台上的链路特性,基于TSN中的门控调度技术,提出了一种适用于卫星链路场景下的动态调度算法。通过改进TAS（time aware shaper）调度算法,减少卫星动平台因通信距离动态变化以及空间摄动等因素造成的传输时延,从而增加动平台上卫星通信系统的确定性。同时根据动态卫星链路周期性运动特征,提出周期性门控调度机制,保证所有周期内相同时间点的门控调度一致。最后通过现场可编程门阵列（FPGA）设计的验证测试结果表明,提出的TAS改进方案可以适用于卫星网络的确定性调度,充分发挥时间敏感网络的优势。     

基于候选点稠密匹配的三维场景重构方法

针对星球探测机器人在未知环境中三维场景重建存在的计算复杂度问题,提出基于网格候选点的三目立体匹配算法。在空间中建立代表深度信息的网格节点,并对深度方向的节点分布进行合理规划,确保候选点稠密匹配的准确性和高效性。候选点匹配解决了传统立体匹配算法中图像校正带来的实时性问题,同时采用三目视觉系统代替双目,通过另一组对应点的相似性测度对潜在歧义的少量候选点进行二次判决。实验证明,由于处理每组图像对不再需要进行极线校正,因此计算代价与传统的匹配算法相比有一定降低,而第三台摄像机有效消除了匹配歧义,计算量相对于双目系统也没有明显增加。

     

自由飞行空间机器人捕捉运动目标的力矩控制算法及其仿真

本文研究了自由飞行空间机器人（ＦＦＳＲ）捕捉运动目标的力矩控制算法及其仿真。首先，建立了ＦＦＳＲ关节驱动力矩的求解算法；其次结合基于广义雅可比矩阵的分解运动速度控制方法提出了一种ＦＦＳＲ捕捉运动目标的力矩控制算法；最后，计算机仿真验证了本文提出算法的正确性。     

高超声速滑翔飞行器再入轨迹快速、高精度优化

针对高超声速滑翔飞行器再入轨迹优化问题，提出一种基于稀疏差分法和网格细化技术的快速、高精度求解方法。该方法应用局部配点法将再入轨迹优化问题转化为非线性规划(NLP)问题，从两方面提高轨迹优化的效率和精度。一方面，引入一种高效的稀疏差分法计算NLP的一阶偏导数，提高NLP的求解效率；另一方面，提出一种基于新型广义二分网格的网格细化算法调整离散节点的数量和分布，使得方法能够采用较少的节点数目取得较高的优化精度，从而减小NLP的规模和计算量。应用该方法求解了高超声速滑翔再入轨迹优化问题，仿真结果表明所述方法能够快速生成一条严格满足各种约束的最优三维再入轨迹。在此基础上，研究了滑翔飞行器的再入落点区范围，进一步检验了该方法的有效性。     

一种有撞击角和视场角约束的运动目标的偏置比例导引算法

针对以一定角度攻击运动目标问题,提出一种采用偏置比例导引的间接撞击角度控制方法。定义满足攻击角度约束的碰撞三角形概念,在这基础上给出了具有偏置比例导引形式的制导律。根据碰撞三角形要求的偏置项积分值设计两阶段和一阶段这两种计算偏置项的方法,其中两阶段方法不需要计算剩余时间,而一阶段方法能够攻击机动目标。考虑视场角和过载等约束条件限制,对偏置项进行三阶段改进,使导弹在满足这些约束条件的同时满足撞击角要求的偏置项积分值。为了适应机动目标的情况,考虑当前时刻导弹和目标运动信息,在线地偏置项进行实时更新。最后对运动目标和机动目标两种情况进行了仿真实验,验证了本文方法的有效性。     

喉栓式推力可调固体火箭发动机动态响应特性数值分析

基于ALE(Arbitrary Lagrangian-Eulerian)描述的N-S方程,利用动网格方法适应边界移动,对喉栓式推力可调固体火箭发动机在推力调节过程中发动机的内流场进行了非稳态数值模拟,分析了喉栓运动速度、发动机自由容积对推力调节性能的影响规律,揭示了喉栓式发动机推力调节过程中发动机的动态响应特性.所得结论可为喉栓式推力可调发动机的设计、试验提供依据.