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

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

AUFS格式在多组分反应流场模拟中的应用

对AUFS格式进行了讨论,并将其应用于高超声速多组分化学反应流场的数值模拟.超声速燃烧流场复杂性使得发展基于非结构网格的高精度逆风格式数值求解方法成为必然.本文中控制方程为轴对称多组分Eu-ler方程,空间分裂采用基于非结构网格的有限体积逆风格式AUFS方法,时间推进为考虑了化学反应特征时间效应的显式5步Runge-Kutta方法.化学反应模型采用经典的H2/Air的7组分8反应模型.对钝头体激波诱导亚爆轰及超爆轰流场进行了数值模拟,得到了与实验及文献相符合的计算结果.     

Roe格式在多组元燃烧流场数值模拟中的应用

对Roe通量分裂格式进行了详细讨论,并将其推广应用到基于非结构网格的高超声速多组元燃烧流场的数值模拟中。控制方程采用三维多组元Euler方程,重点针对Roe格式开展研究;考虑到燃烧流场求解过程中出现的刚性问题,对化学反应源项采用点隐式方法处理。对超声速条件下钝头体激波诱导燃烧流场进行了数值模拟,比较分析了3种经典化学动力学模型对流场结构的影响,获得与数值模拟、试验数据相符的结果。结果表明发展的方法可以用于多组元燃烧流场的数值分析,且化学动力学模型的选取直接影响激波诱导燃烧流场中诱导区厚度。     

不同进口马赫数下超声速燃烧流场特性研究

对有无楔板超燃冲压发动机模型内横向氢气喷流超声速燃烧流场进行了数值模拟,分析了进口马赫数对超声速燃烧流场特性及特征参数分布的影响特性。采用有限体积法求解多组元Navier-Stokes(N-S)方程,对不同进口马赫数下的燃烧流场进行了数值模拟,细致对比了流场激波结构、喷流穿透深度、燃烧阵面,燃烧效率及总压恢复系数等参数随进口条件的变化特征。结果表明:无论是否存在楔板结构,喷口后流场压强均随着进口马赫数的增加而减小,并且随进口马赫数的增加,氢气喷流穿透深度减小,楔板对喷流穿透深度基本无影响。较无楔板结构而言,设置楔板结构可以缓解燃烧室内流场对马赫数变化的敏感度,使燃烧更为稳定。在同一进口马赫数条件下,楔板布局有明显的促燃作用及激波点火效果,在一定程度上可增加此类发动机工作的马赫数范围,但以总压恢复系数略微降低为代价。     

HLLE+格式在高超声速热化学非平衡流场中的应用

对HLLE+格式进行了讨论,并将其应用于高超声速热化学非平衡反应流场的数值模拟.本文控制方程为带化学反应源项的多组元轴对称Navier-Stokes方程,为了保证数值模拟的精确度,空间分裂采用基于结构网格的有限体积迎风格式HLLE+方法,时间推进为隐式的LU-SGS方法.化学反应模型采用经典的Dunn-Kang反应模型,并且采用Park的双温模型模拟流场的热力学非平衡效应.对高超声速条件下的RAM-CII模型进行了数值模拟,得到了与实验以及文献相符合的结果,验证了对HLLE+格式扩展的正确性.     

弹射座椅系统外高速流场的数值模拟

对绕弹射座椅系统的高速流场进行了数值模拟.控制方程为雷诺平均的Navier-Stokes方程,空间离散为Osher逆风格式和五步龙格库塔格式,湍流模拟采用了基于M-SST两方程紊流模型的DES方法,计算中采用了区域分裂的并行技术.对来流马赫数分别为0.6,1.2的流场进行了计算,计算结果表明该方法可以较好地模拟弹射座椅系统外的高速流场.     

双级涡流器环形燃烧室整体流场数值模拟

在任意曲线坐标系下对包括扩压器、双级涡流器及火焰筒在内的环形燃烧室三维两相燃烧整体流场进行数值模拟。由于环形燃烧室形状复杂,采用偏微分方程法与整体分区结合法生成环形燃烧室整体网格。所用的数学模型有:k-ε紊流模型、EBU-Arrhenius紊流燃烧模型、六通量热辐射模型以及颗粒群轨道模型等。在非交错网格系下,气相采用SIMPLE算法求解,液相采用PSIC算法求解。数值分析不同燃烧室进口气流参数以及涡流器几何尺寸对燃烧室流场的影响,并将计算结果与实验数据进行了比较。结果表明本文的计算方法合理性可用于环形燃烧室的研制与优化设计。     

三维非结构混合网格高超声速流场并行计算

在非结构混合网格上对三维高超声速化学非平衡粘性绕流进行了基于PC-Cluster的分布式并行数值模拟.本采用区域分裂思想,研究了三维非结构混合网格区域自动分解技术,并以此为基础对高超声速化学非平衡绕流进行了并行数值计算.控制方程为多组分N-S方程,空间离散采用有限体积格心格式,时间推进为显式Runge-Kutta格式.化学非平衡动力学模型为七组元带电离反应模型,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法.在PC-Cluster上对三维双椭球模型的高超声速绕流流场进行了基于区域分解技术的并行数值模拟,所得数值结果与参考献中的结果作了对比验证。     

辐射及化学非平衡流耦合场计算方法

把一种新的三维热辐射计算方法应用到高超声速非平衡流场的数值模拟过程中:在非结构网格上对高超声速化学非平衡流场的数值模拟过程中,耦合了辐射输运方程来求解辐射热流。控制方程为含有化学反应源项和辐射源项的三维N-S方程,数值离散格式采用了Jam eson有限体积法,辐射输运方程采用有限体积法求解。化学反应模型为7组元7反应模型,并采用了化学反应特征时间步长。本文分别对RAM C-Ⅱ和MU SES-C模型进行了计算,计算结果与参考文献的结果吻合。     

边界条件对三维空腔流动振荡的影响

基于全隐式无分裂方法,数值求解三维非定常雷诺平均的N av ier-S tokes方程,紊流模型采用修正的B-L代数模型。运用此方法,本文数值模拟了三维开式空腔的超声速、跨声速、亚声速流动。研究了空腔口处剪切层的非定常运动现象,着重探讨了入口边界条件对三维空腔流自持性振荡的影响。数值计算结果表明,与层流速度型相比,入口条件为紊流速度型和特征边界型得到的结果能够吻合。