自适应间断装配法模拟弹道靶中超高速弹丸发射

弹道靶自由飞试验是高超声速领域中的一种重要的地面试验手段。试验中厘米级的弹丸/弹托在10米级长度的管道中加速至极高的速度，出膛后通常采用气动力使弹托偏离预定弹道而被拦截器阻挡，仅使弹丸进入试验段开展测试。由于发射过程中弹丸的实际飞行姿态会受到管道内长距离加速、高压驱动下出膛、气动分离干扰等诸多因素的影响，因此基于计算流体力学发展相关的气动预测技术有助于指导试验方案设计，从而保证弹丸运动轨迹的准确性。针对发射问题中涉及的诸如空间尺度变化剧烈、接触-分离、超高速动态分离等数值仿真难点，采用基于非结构动网格技术和格心型有限体积方法发展的自适应间断装配求解器（ADFs），对非定常流场中的运动激波进行装配，通过二维算例对弹丸发射过程数值仿真进行了详细的原理性介绍。一方面，拓展了激波装配方法在工程问题中的应用；另一方面，针对弹道靶中超高速弹丸发射这类问题，建立了一套高效的数值模拟方法，实现了对弹丸从静止到加速、出膛、分离的全动态过程的精细流场刻画。     

一种基于聚类分析的二维激波模式识别算法

在激波捕捉求解器计算的可压缩无黏流场基础上，提出了一种探测并识别二维激波干扰模式的新算法，从3个层面详细介绍了该算法的实施流程。首先，采用基于当地流场参数设计的传统激波探测方法，辨识出激波附近的一系列网格单元；其次，通过经典的K-means聚类算法将这些激波单元划分成许多簇，并根据簇的相邻信息定义每个簇的类别；最后，设定相关准则对某些紧邻的簇进行合并，进而确定各个激波干扰点的位置，记录各条激波分支所对应的簇，采用Bézier曲线拟合算法分别对其聚类中心进行拟合以获取更加光滑的激波线。数值试验表明，该算法不受网格类型的限制，不仅可以保证最终拟合的激波线具有较高的位置精度，还可以清晰地识别出流场中多激波干扰的模式，同时对分析非定常流场中激波的运动与演化过程也提供了一种有效的可视化手段。     

一种改进的非定常激波装配算法

基于非结构动网格技术和格心型有限体积方法,提出一种改进的非定常激波装配算法,进一步拓展了其在包含有运动激波的非定常流场的应用范围。首先,针对激波在直/曲壁面传播这类问题,分别建立了壁面间断节点的运动模型;其次,为保证激波在大范围运动时装配阵面不产生失真,基于Bézier曲线拟合方法实现了间断节点分布的自动重构;接着,通过嵌入局部网格自动重构模块,提高了算法的计算效率和自动化程度;最后,对于激波相交点的运动,设计了一种根据位移推算速度的方法进行装配。数值算例表明,所提算法能够有效地处理激波传播问题,相比激波捕捉方法可以提取更多的流场信息,同时可以获得流场间断更加直观清晰的图谱。     

翼反角对高压捕获翼构型高超气动特性的影响

为研究翼反角变化对高压捕获翼构型高超声速气动特性的影响，基于一种双翼面、单支撑、翼身组合布局的高压捕获翼概念构型，以飞行马赫数6，飞行高度30 km为计算状态，捕获翼和机体三角翼上/下反角为设计变量，结合均匀试验设计方法、数值模拟方法和Kriging建模方法，探寻了升阻特性、纵向和横航向稳定性随翼反角的变化规律。结果表明，升力、阻力及升阻比随翼反角的变化规律基本一致，且对上反角变化更加敏感；小攻角时，翼面上反会明显降低升阻比，而下反会使升阻比先略微增大后缓慢减小；大攻角时，翼反角对升阻比的影响较小；纵向稳定性主要受三角翼反角的影响，三角翼上反时，纵向稳定性降低，下反时，纵向稳定性基本不变；翼面上/下反都会提高航向稳定性，但下反的效果更明显；翼面上反会提高横向稳定性，下反则降低，但大攻角飞行时，三角翼上反角过大可能会导致横向稳定性降低。     

机体尾缘形状对高压捕获翼构型亚声速特性影响

基于圆锥-圆台组合平板捕获翼构型，通过改变尾缘展向扩张角，获得一系列不同外形，在典型亚声速（Ma=0.5）来流条件下开展数值计算，并重点分析了机体尾部截面形状和攻角变化对流动特性和气动特性的影响。结果表明：在0°攻角状态下，机体尾截面展向变宽，机体与捕获翼之间的流场区域对来流的扩张减弱，机体圆台上表面的逆压梯度减小，可有效抑制机-翼之间流场内的流动分离现象，同时整机升力系数增大，阻力系数先减小后增大。随攻角增大，机体圆台上表面压力增大，分离区范围逐渐缩小直至消失，机体尾截面展向变宽可加速分离区消失的进程。当攻角进一步增大时，机体背风面出现横向绕流，但机体尾截面展向变宽可以延缓横向绕流的发展。计算结果还表明，随攻角增大整机升力及阻力主要由捕获翼贡献，机体贡献的气动力随攻角变化不敏感，机体尾截面展向变宽对整机焦点位置影响较小。机体下表面几何形状变化对机体与捕获翼之间的区域内流动特性和捕获翼部件的气动力特性无明显影响。