跨流域高超声速绕流Boltzmann模型方程并行算法

通过对Boltzmann方程碰撞积分进行模型化处理,提出了统一描述各流域复杂高超声速流动输运现象的气体分子速度分布函数控制方程,使用离散速度坐标法对分布函数方程所依赖的速度空间离散降维,构造出直接求解分子速度分布函数的气体动理论耦合迭代数值格式,研制了复杂飞行器高超声速绕流气动热力学计算模型。基于对气体动理论数值计算方法内在并行性、变量依赖关系、数据通信与并行可扩展性的分析研究,使用区域分解并行化方法提出了新型的气体动理论数值算法并行方案;研究了数据的并行分布与并行执行特征,开展了大规模的并行化程序设计,构造了可稳定运行于成千上万CPU的高性能并行算法,用以模拟各流域复杂飞行器的高超声速绕流问题。以稀薄流到连续流环境下不同Knudsen数、不同马赫数的可重复使用类球锥卫星体及翼身组合复杂飞行器等气动力、热绕流问题为研究对象展开大规模并行计算,并进行算法验证,所得计算结果与理论分析、直接模拟蒙特卡罗方法(DSMC)的模拟值及有关实验数据吻合较好,揭示了飞行器跨流域高超声速下的复杂流动机理与变化规律,提供了一条能够可靠模拟高超声速飞行器跨流域气动力及热问题的统一的算法应用研究途径。     

基于Boltzmann模型方程不同流区复杂三维绕流HPF并行计算

随着现代航空航天快速发展,是否能建立一套可有效模拟稀薄流到连续流各流域绕流问题统一算法,已成为工程应用部门和学术研究领域关心的问题。通过开展基于Boltzmann简化速度分布函数方程数值求解,发展起从稀薄流到连续流统一算法。借助区域分解并行化方法研究建立气体运动论统一算法并行方案,通过对统一算法进行HPF并行化程序设计及算法考验,拟定不同流域三维球体绕流及类"神舟"返回舱外形体绕流算例,进行HPF并行计算,并将计算结果与有关实验数据、DSMC模拟值进行定量比较、分析。研究表明,所发展的统一算法具有很好的并行独立性,基本达到了线性加速的并行效果,算法负载平衡和并行可扩展性较好,可望建立起新型的能可靠模拟不同流域三维绕流问题的HPF并行算法研究方向。     

激波结构内流动问题的气体运动论描述

从介观Boltzmann速度分布函数理论出发,发展计及分子粘性碰撞截面与扩散碰撞截面,可描述各流域一维气体流动问题的Boltzmann简化速度分布函数方程及其气体运动论数值计算方法。通过对不同Knudsen数下非定常激波管流动及不同马赫数定常正激波结构问题数值模拟,研究分析不同流区的激波突跃变化过程以及近连续流、稀薄过渡流特有的分子输运现象,揭示不同马赫数、不同分子模型的激波内流动与传热变化规律,证实基于Bo-ltzmann模型方程的气体运动论数值计算方法用于激波结构内流动研究的准确可靠性。     

跨流域三维复杂绕流问题的气体运动论并行计算

通过研究求解描述跨流域三维绕流问题的Boltzmann模型方程气体运动论耦合迭代数值格式,分析气体运动论数值计算方法的内在并行性;从变量依赖关系、数据通信与并行可扩展性三方面开展基于离散速度空间区域分解计算的研究,发展求解稀薄流到连续流跨流域三维复杂绕流问题的并行算法。通过对不同Knudsen数、不同马赫数、不同攻角三维球体及返回舱绕流的并行计算,计算结果与实验数据和理论分析吻合较好。研究表明,该并行算法负载平衡和并行可扩展性较好,对不同并行计算机系统具有很好适应性,显示直接求解分子速度分布函数的气体运动论计算方法有良好的并行计算特性。     

逆向喷流对双锥导弹外形减阻特性的影响

逆向喷流是一种主动流动控制技术，具有减阻降热特性，可用于高超声速飞行器设计。以典型双锥导弹外形的球头、单锥、双锥（全弹）为研究对象，将喷流发生器和弹体固连，采用CFD方法对逆向喷流的减阻特性进行了数值研究，对比分析了喷流马赫数、喷流压比等参数对不同对象减阻效果的影响。结果表明：逆向喷流流场存在长、短射流穿透两种模态；球头在小压比长射流模态时的减阻效果最佳；单锥和双锥在大压比短射流模态时的减阻效果更好。存在一个最佳压比，使得逆向喷流的减阻效果最佳；喷流压力过大，减阻效果变差，甚至出现阻力系数不降反增情形。逆向喷流减阻效果对控制体选取敏感，若将逆向喷流对头部的减阻特性（超过40%）直接推广至飞行器整机（6%左右），评估结果过于乐观。综合最佳减阻效果、最佳喷流压比、流量与所需储气瓶体积等影响因素，工程应用时逆向喷流应优先选用声速喷流。     

红外热图在高超声速低密度风洞测热试验中的应用概述

围绕大极角情况下模型表面温度的测量、红外热图测热精度的提高、模型物面坐标与红外热图像素位置对应关系等问题,对近年来在高超声速偌氏密度风洞开展的红外热图测热工作进行了总结,给出了一些实用、简单、经济的方法.为验证上述技术,在高超声速低密度风洞开展了用红外热图技术与热电偶同时测量一平板带劈薄壁模型表面的气动加热率分布以及半球圆柱模型红外测热数据与DMSC计算结果的比较.不同手段获得的数据与红外测热数据相互验证的结果表明:这些技术的解决,有利于红外热图技术向工程化实用化迈进.     

离散速度方法及其在跨流域问题中的应用研究进展

往返大气层跨流域飞行器的气动力和气动热计算一直是当前流体力学研究的难点和热点之一。由于该类飞行器面临的流动场景不再是单一的连续流或稀薄流,采用Navier-Stokes方程求解器或DSMC方法均不能获得全流域的准确结果。近年来,以不依赖于连续性假设的Boltzmann模型方程为基础,通过在位置空间和速度空间同时离散求解该方程,实现了跨流域问题的统一求解。本文对该类算法的研究进展进行回顾和分析,着重介绍气体动理论统一算法（Gas kinetic unified algorithm, GKUA）、统一气体动理学格式（Unified gas kinetic scheme, UGKS）和改进离散速度方法（Improved discrete velocity method, IDVM）3种数值途径,分析它们的基本假设和实现方式,关注它们目前取得的进展和应用情况。同时,本文还将IDVM进一步扩展到非定常情形,以便用于非定常跨流域问题求解。最后,本文对该类算法存在的一些问题进行讨论,期望在后续的研究中能予以解决。     

阿波罗指令舱稀薄气体动力学特征的蒙特卡罗数值模拟

本文对阿波罗指令舱稀薄气体动力学特征进行了蒙特卡罗直接模拟（ＤＳＭＣ）。采用中点数值逼近法确定模拟分子与物面碰撞点，利用网格编号将三维计算问题一维化，把流场入口处的分子补充与分子排序有效结合，利用分子横穿流场次数确定流场稳定前的时间步长计算次数。简化了计算过程，大大节约了计算机时和内存。给出了指令舱从自由分子流到过渡区不同Ｋｎｕｄｓｅｎ数、各攻角下的气动特征。模拟结果与有关理论、实验数据相比，证实     

翼伞后缘偏转过程的流固耦合动力学特性

深入研究翼伞后缘偏转过程的气动与结构耦合动力学问题是解决大型翼伞精确空投系统机动转弯和雀降等操纵动作设计分析的重点内容。首先基于ALE算法和罚函数耦合方法对翼伞后缘偏转过程进行流固耦合动力学建模,之后基于结构化的ALE求解方法和瞬态非线性求解器对翼伞三维模型的单个气室后缘偏转进行仿真验证,预测了后缘偏转运动引起的周围流场流动分离现象。分别针对翼伞后缘单侧下偏和双侧下偏过程的流固耦合行为进行仿真分析,获得全时域内翼伞结构场和周围流场特性动态演化结果,以及翼伞气动性能参数时间历程曲线,发现了后缘下偏过程的操纵延迟现象。最后通过风洞试验对仿真结果进行验证,证明了方法的有效性,为大型冲压翼伞的设计和应用提供理论和技术支撑。     

基于模糊关系模型的电火花加工工艺数据库

]