非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。     

激波诱导燃烧模拟中轴对称边界数值异常研究

对钝体激波诱导燃烧模拟中引起对称轴附近出现数值异常现象的几个影响因素进行了细致讨论.在飞行速度低于爆轰速度情形的模拟中,在一定条件下靠近对称轴的位置会出现明显的数值异常——燃烧阵面扭曲.分别从边界处理方法、计算格式、化学反应机制和计算网格等4个方面对该异常产生的原因进行了详细地讨论.研究表明4个因素都会导致数值异常产生,深入的分析表明计算格式是其中最根本的因素,而这种燃烧阵面扭曲的数值异常应该与粉刺问题(Carbuncle)类似.针对如何克服轴对称边界数值异常提出了一些合理的措施,并给出了验证算例的对比与分析.      

高超声速飞行器内流道流态及其对全机气动力影响的实验研究

对吸气式高超声速飞行器而言,内流道的流态对推进系统乃至飞行器的工作特性有重要影响.本文用风洞实验的方法研究了升力体(类X-43)、二元进气道和轴对称进气道三种构型在不同来流、内部几何和燃烧室反压条件下的内流道工作情况,着重分析了起动/不起动状态对内流道流动特征及全机气动力的影响,提出并发展了通过在燃烧室内喷射高压气流的内流道反压模拟技术,并运用这一技术成功实现了内流道由起动向不起动、由不起动向起动/再起动的切换.     

预处理AUSM+格式在全速域反应流模拟中的应用

对预处理修正的AUSM+(advection upwind splitting method+)格式进行了详细讨论, 并将其应用于全速域范围内的化学反应流场数值模拟.基于多组分Navier-Stokes方程, 结合控制方程时间导数项的预处理, 采用预处理修正的AUSM+格式发展了一种模拟全速域范围内化学反应流场的有效算法.验证算例包括超声速、亚声速以及极低速情形的化学反应流场, 计算结果表明基于预处理AUSM+格式发展的算法能够有效地用于全速域范围内化学反应流动问题的模拟.      

先进飞行器气动隐身一体化设计方法研究现状与趋势

介绍了先进飞行器气动隐身一体化设计的研究现状与发展趋势。首先介绍了设计中的常用的数值模拟技术。针对气动力计算,重点分析对比了位势方程、Euler方程和N—S方程等气动力计算方法;针对隐身性能计算,重点分析对比了高频近似方法、时域法、频域法等雷达散射截面计算方法。在此基础上。讨论了先进飞行器气动隐身一体化设计的策略和步骤,最后对该设计方法的发展趋势进行了分析。     

高超声速轴对称流道冷流特征及气动力特性研究

对一种轴对称形式的高超声速飞行器全流道开展了风洞实验和数值模拟研究, 分析了不同来流总压、飞行攻角全流道的流场结构和气动力特性.研究结果表明:(1)一定范围内雷诺数的变化对全流道的流动结构和模型的气动力特性无显著影响, 因此所获得的风洞实验结果有望通过某种形式推广到飞行状态下使用;(2)飞行攻角对全流道的流动结构和升力系数有着显著影响, 但阻力系数的影响并不明显;(3)研究范围内来流马赫数的变化对全流道的流动结构有着一定影响, 但研究范围内, 阻力系数随马赫数的变化幅度较小;(4)由于轴对称流道的浸润面积较大, 研究范围内该类飞行器的摩擦阻力在全机阻力中占据了较大的比重, 设计状态下达全机气动力的62%;(5)与实验结果的对照表明, 所采用的数值模拟方法具有较高的精度.      

多级压缩锥导乘波体设计与分析

提出了多级压缩锥导乘波体的设计方法,该方法应用吻切锥理论和零攻角圆锥绕流基准流场通过流线追踪生成具有多个压缩面的乘波体。对以吸气式冲压发动机为动力的高超声速飞行器,应用多级压缩乘波前体可充分发挥前体的预压缩作用,为进气道的正常工作提供所需的均匀流场。以二级压缩乘波体为例阐述了该设计方法,设计方法通过对二级压缩基准流场进行重构,使其符合Taylor-Maccoll流动模型以获得新的二级压缩基准流场。同时编写设计程序生成了一级、二级和三级压缩乘波体,通过数值模拟结果校验设计方法的正确性,并对其压缩性、升阻比、总压恢复系数等性能进行了对比分析。     

机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

基于动态非结构重叠网格法的直升机前飞非定常流场数值模拟研究

 发展了非结构重叠网格技术，提出适合于直升机旋翼运动的动态非结构重叠网格算法，该算法使得在采用非结构重叠网格模拟直升机流场时实现程序自动处理网格，不需任何人工干预。流动控制方程采用非定常可压缩Euler方程，空间离散采用Jameson有限体积法，时间推进采用双时间步长法，为加速收敛，在双时间步长内循环迭代中采用预处理技术。为验证本发展的预处理解算器以及动态非结构重叠网格的算法的可信性和有效性，开展一系列的数值实验，模拟了不同速度范围的NACA0012翼型绕流,Caradonna旋翼悬停和前飞流场以及Robin直升机机身-旋翼组合的前飞非定常流场，结果表明采用本文算法得到的结果是可信的，网格处理是自动高效的，可以应用于直升机全机流场模拟，为直升机设计提供可信指导数据。