超声速化学反应流场迎风格式数值模拟及并行计算

基于有限体积迎风格式对超声速燃烧流场进行了的数值模拟.由于超声速燃烧流场绕流的复杂性,要求对多组分Euler/N-S方程求解的数值模拟方法应具有较高的计算精度及效率.本文引用辅助点方法建立了具有空间二阶精度的van Leer迎风矢通量分裂格式,并应用于超声速燃烧流场绕流的数值模拟.化学反应为氢气/空气十反应模型,采用考虑了化学反应特征时间的当地时间步长显式Runge-Kutta时间推进格式.对钝头体模型爆轰现象、后向台阶氢气喷射及二维内外流超声速燃烧流场模型进行了区域分裂技术的并行计算.计算结果与参考文献作了对比,得到了满意的结果.     

非结构网格的高超声速粘性MHD流场数值模拟

应用数值方法在非结构网格上对磁场干扰下的二维高超声速钝头体粘性绕流进行了数值模拟。控制方程为N-S方程耦合Maxwell方程的粘性MHD方程组,空间离散采用有限体积方法,对流项用AUSM格式计算,粘性项用中心格式求解,时间推进用显式5步Runge-Kutta格式,引用双曲型散度清除技术加强▽.B=0的条件。计算模型为二维钝头体,在高超声速来流条件下,分别对有、无均匀磁场干扰下的流动进行了数值模拟。计算结果表明,在均匀磁场干扰下,激波脱体距离显著增加,物体表面压力急剧下降。对比表明文中发展的计算方法可以准确地进行二维粘性MHD流场的数值模拟。     

分布式进化算法及其在翼型气动反设计中的应用

构造了一种新型基于基因算法与博弈论的并行分级多目标优化方法，并应用于多段翼型气动反设计。此方法基于二进制编码的基因算法和博弈论，优化变量被分配给不同的博弈者，因而总体优化问题转变为分裂空间中的局部优化问题。文中给出了一个多段翼型形状，位置可压位流的反设计问题的求解算例，引入了基于非结构网格的分级结构。与传统基因算法数值算例的对比表明了本文构造的并行分级算法具有较高的计算效率，可广泛应用于多目标优化问题。     

带舵旋成体侧向喷流流场特性分析

针对带舵旋成体外形,来流马赫数为2.0,迎角为0°,10°和20°,喷口位于×字舵中间,首先开展了喷主流干扰流场的数值模拟与风洞实验校验,得到相符的喷口附近子午线压力分布,验证了数值模拟方法的可行性;然后根据数值模拟得到的喷主流干扰流场结构和流动参数分布,分析了喷主流干扰流场结构随迎角的变化特点;最后从喷口附近旋成体和尾舵表面压力分布的变化分析了喷主流干扰对气动力/力矩放大因子的影响.研究结果表明:随迎角增大,喷流前激波向喷口靠近,分离区减小,喷主流压比增大,喷流穿透高度增大;气动力/力矩放大因子均大于1.0,随迎角增大,放大因子减小.喷口位于×字舵中间时,喷主流干扰对喷流直接力有增益作用;×字舵的存在减弱了喷流的包裹作用,对流场起到了有利干扰.      

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

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

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法。外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格 库塔方法。采用四阶龙格 库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率。使用非结构动态网格技术对积冰后的翼型网格进行快速调整。在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性。     

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

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

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

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

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

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