高超声速化学非平衡绕流分布式并行算法

讨论了非结构网格高超声速化学非平衡绕流流场的分布式并行数值模拟方法。控制方程采用带化学反应源项的N-S方程,数值离散格式采用Jameson有限体积法。化学反应模型为五组元五反应模型,没有考虑电离效应,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法。通过调用METIS库函数将整体计算区域进行分裂,使用PVM并行机制在PCs-Cluster机群上对二维钝头体高超声速化学非平衡绕流进行了分布式并行计算,得到了结果,并对算法的加速特性进行了对比分析。     

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

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

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

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

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

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

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

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

高超飞行器内流道激波振荡问题的数值研究及试验验证

采用数值模拟和风洞试验的方法,对吸气式高超声速飞行器盲腔状态的流动特性进行了研究.采用“双时间步”方法进行了内外流一体化的非定常数值模拟,利用彩色纹影系统对高超声速飞行器前体流场进行显示,并采用动态压力传感器测得了飞行器内流道的壁面压力随时间的变化.结果表明:盲腔构型在高超声速飞行中会出现周期性的激波振荡现象.数值模拟所得流场变化特征、内流道壁面压力振荡周期和壁面压力变化趋势与试验结果吻合良好.     

基于混合网格上的高温非平衡流的数值模拟

新一代再入飞行器及空间传输系统需要了解其飞行时的气动问题及热力学问题。数值模拟高空、低密度、高熵及非平衡流动是一具有挑战性的问题。本文针对带座舱飞船高超声速再入大气过程中存在的严重气动加热现象,利用混合网格及Osher逆风格式,数值求解了三维化学非平衡Navier-Stokes方程,其中化学模型是5组分1 7个化学反应的空气化学模型,对带座舱飞船再入高度为40 km和马赫数为2 0 ,1 0的化学非平衡流场进行了数值模拟,给出了迎角为0°和2 0°情况下的各个组分的密度分布、压力等参数,并与量完全气体的计算结果进行了比较。     

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

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

带减阻杆高超声速飞行器外形气动特性研究

采用高超声速风洞测力试验方法测量钝头飞行器头部减阻杆的高超声速气动特性,研究减阻杆的气动减阻原理,分析了多组不同构型减阻杆的减阻效果.结果表明,减阻杆显著减少了钝头飞行器高超声速的阻力,最大的减阻率达到60%之多;减阻效果与减阻杆构型和迎角状态密切相关;减阻杆会诱发稳定性、“热斑”以及非定常脉动等不利问题.     

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

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

热化学非平衡流场辐射特性数值计算

设计了耦合辐射的三维热化学非平衡流场计算方法,以适用于任意形状的网格。采用Jameson有限体积法求解耦合辐射源项的三维N-S方程。化学模型包含11个组元,20个化学反应。辐射源项通过直接求解辐射输运方程RTE获得。在空间和方向上分别离散后,利用有限体积法求解不同方向上的辐射输运方程,精确计算了波长为100-1500nm所有主要组元的吸收系数,逐条计算了目前已观测到的空气组元的原子谱线,而对分子的带状谱采用了带模型计算,并最终得到了驻点处辐射密度随波长变化的分布情况。     

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

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

盒式机翼布局气动特性研究

采用数值模拟的方法研究了盒式机翼布局的气动特性。在数值模拟过程中,控制方程为雷诺平均N-S方程,采用基于非结构网格的有限体积方法进行求解。空间离散采用了Osher逆风格式,时间离散则采用五步Runge-Kutta格式。紊流模型采用了S-A方程模型。对盒式机翼布局外形在不同状态下的流场进行了数值模拟,给出了每个状态下的气动参数,并进行了分析,结果显示盒式机翼布局具有较好的气动特性。     

含化学反应流动的TVD格式及其在燃气射流中的应用

以Harten标准TVD格式为基础，结合固体火箭燃气射流的特点，以数学方法系统地推导了适合于高温、高压和高速流体流动的数值格式，给出了压力偏导数的合理计算公式；利用化学动力学知识，对火箭燃气射流流场中存在的多组分、含有限速率的化学反应问题进行了论述，阐明了化学反应质量源项的求解方法。以12组分9反应方程模型为例，利用编制的计算程序，对某火箭燃气自由射流流场进行了模拟。通过对结果的分析，肯定了数值格式的正确性。     

高超声速圆球模型飞行流场的数值模拟和实验验证

对高超声速圆球模型飞行流场进行数值模拟 ,分别采用空气完全气体模型、平衡气体模型以及热化学非平衡 1 1组元气体模型求解非定常轴对称N -S方程组。使用有限差分时间相关法捕捉激波 ,得到了定常流场的解。差分方程隐式部分采用了LU -SGS方法以避免矩阵运算 ,对化学反应和振动能量源项采用预处理矩阵以解决刚性问题。由计算结果处理得到的阴影图和干涉条纹图与再入物理弹道靶实验照片进行了对比分析 ,验证了实验中圆球飞行流场大部分区域接近于平衡状态。     

非结构网格二维理想MHD绕流逆风格式解法

在非结构网格上发展了针对二维理想磁流体方程组的逆风格式求解方法.控制方程中的对流项采用AUSM格式处理,时间推进采用显式5步龙格-库塔方法.为了消除计算中产生的磁场散度的影响,引入了双曲型散度清除方法.通过对磁流体激波管问题的求解验证了该方法对激波的捕捉能力,对有均匀磁场干扰下的喷管流动情况进行了数值模拟,并与文献中结果进行了对比.计算结果显示了磁场对磁流体流动的干扰效应,该结果与参考文献中的数值模拟结果相吻合.     

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。