时空守恒CE/SE数值方法在凹槽燃烧流动中的应用研究

应用时空守恒元和解元CE/SE方法,发展了求解化学反应流场的三维并行计算程序。为验证程序,对超声速冲压喷射燃烧室模型进行了计算,计算结果与实测数据以及相关文献的计算结果吻合较好,验证了本文的处理方法和程序的正确性。采用H2和O2的8组分34基元化学反应模型,对不同长深比的凹槽燃烧流动进行了数值模拟,分析了凹槽增强混合对稳定火焰的作用。计算结果表明,CE/SE方法具有模拟凹槽超声速燃烧流动的能力,数值模拟的结果有助于认识凹槽的火焰稳定机理。     

用非结构网格和欧拉方程计算多段翼型绕流

用Ｄｅｌａｕｎａｙ法划分三角形构造非结构网格时，如果布点不当，会产生边界穿透，使网格生成失败，本文分析了几种布点情形，提出不容易产生边界穿透的布点方式。这一原理已用于生成绕多段翼型的非结构网格。     

非结构网格及混合网格复杂无粘流场并行计算方法研究

在分布式存储网络异构编程PVM环境下,本文将原串行计算程序改造为并行计算程序,提高了计算效率.利用非结构网格,对球体及航天飞机外形绕流流场进行了计算,取得了较好的性能结果.最后将该方法应用于混合网格的计算,求解了多种复杂外形的无粘绕流,如翼身组合体和DDF1C导弹,得到了满意的结果.     

非结构网格上N—S方程解法研究及在多段翼型绕流中的应用

在非结构网格上求解N－S方程时，靠近物面的网格伸展比相当大，用有限体积法进行显式求解时不易收敛并且雷诺数有一定的限制。本文采用的人工耗散项中含有网络的方向尺度，因而在不同区域选取不同：在粘性区内该项比物理粘性项小得多，在粘性区外与Jameson的人工耗散项相同。生成非结构网格采用的是层推进技术与传统的阵面推进法相结合的方法，这样能保证在粘中性区内网格有一定的方向规则，而在粘性区外网格完全非结构化。     

非结构网格上Euler方程的区域分裂算法及并行计算

改进了“波阵面”区域分裂算法,并应用于流场区域的划分;对于子区域边界的不光滑现象,为尽量减少通讯消耗,提出了一种边界并行优化策略。利用PVM并行环境,探讨了非结构网格上求解Euler方程的分区并行算法。根据改进的区域分裂算法及优化策略,运用Jameson有限体积法,对二维翼型流场进行了分区并行求解,多区计算的结果与单区计算的结果作了比较,表明了本文研究方法的有效性。     

二维Euler方程的非结构网格再生成自适应方法

给出了提高二维Ｅｕｌｅｒ方程定常解质量的非结构网格自适应方法和模拟结果。流场算法为结点型有限体积法和Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式，网格用阵面推进法生成，生成初始网格时使用均匀的直角坐标背景网格。     

基于非结构网格的多段翼型流场数值模拟

将非结构网格生成技术应用于多段翼型流场数值模拟，针对多连通域可以快捷、迅速地生成非结构网格。通过控制背景网格节点的网格生成参数，进行网格的局部加密，提高多段翼型流场分析的精度，应用有限体积法求解欧拉方程，应用各种加速收敛措施提高流场求解速度，并将数值计算结果与风洞试验结果进行了对比分析，证明该方法具有型号设计应用价值。     

有、无障碍物爆震管起爆过程的CE/SE模拟

用时空守恒元和解元(CE/SE)法模拟有无障碍物情况下爆震管二维爆震波的形成和传播过程,分析爆震波的不同机理,研究在较低能量下,点火区温度、爆震管管径以及障碍物对爆燃-爆震转捩过程(DDT)的影响.控制方程采取无量纲形式,用拓展的隐式梯形法实现刚性源项积分.研究表明,点火能量与点火压力相关,在较低的点火能量下,可以通过提高点火温度、减小爆震管径及加入障碍物加强爆震过程,减小爆燃-爆震转捩过程DDT时间和距离.为研究爆震波点火和DDT过程提供了理论依据.      

三维Euler方程的自适应非结构网络计算

介绍一种三维自适应非结构网格的Euler方程求解技术。采用有限体积中心差分的格点格式对Euler方程进行空间离散,多步龙格-库塔时间推进,并采用当地时间步长、焓阻尼修正等加速收敛。采用Delaunay三角剖分技术生成三维流场的四面体网格。网格自适应技术中的网格局部加密也采用Delaunay三角剖分技术。最后用ONERAM6机翼的跨音速绕流计算说明本文方法的有效性。     

多管脉冲爆震发动机流场数值研究

针对二维带有收敛扩张喷管的双管脉冲爆震发动机(PDE),采用并行的时空守恒元和解元(CE/SE)方法对单管点火、双管同时点火和不同喉道面积情况下的爆震流场进行了数值模拟。化学反应模型采用H2和O2的8组分34基元反应。计算发现,单管点火情况下传入旁通爆震管的流场可以引爆其中的可燃气体;多管PDE内流场之间的相互作用使管内的压强处于比较高的状态;增大喉道面积削弱了多管PDE流场之间的相互影响,同时减小了封闭端面的压力,从而导致发动机性能下降。     

三维Euler方程组隐式LU分解的高性能并行计算

许多非定常无粘流体力学问题的数值模拟都需要利用Euler方程组来进行计算,而由于在隐格式下,所选取的时间步长可以比在显格式下时大得多,所以隐格式越来越受到重视,其中隐式LU分解是最常用的方法之一。对三维Euler方程组,采用隐式LU分解进行计算时,网格点所在的各个对角阵面之间存在数据依赖关系,本文分析了采用区域分解且边界上用显格式代替隐格式进行计算的高效性,在长方体建筑物内的爆炸模拟表明,在有112个CPU的某MPP巨型机上,并行计算效率超过60%。本文还分析了计算结果与串行计算时的差异,以及利用区域重叠减小这种差异的方法,同时考虑了对处理器进行合理的逻辑组织,将计算网格映射到处理器网格,以最大限度减少通信开销的方法。文中最后以一个爆炸毁伤的例子实际说明了所述方法的可行性与高效性。     

非定常不可压粘性/无粘性耦合方程的一种分步分解方法

给出了数值求解初始变量不可压Navier-Stokes／Euler耦合方程的一种分步块LU分解方法。与传统的时间分裂法不同,该法无需压力中介边条件,从而避免了传统时间分裂法要求的复杂的压力中介边条件逼近。分步块LU分解方法可看做经典的Uzawa算法的改进,后者曾被成功应用于不可压Navier-Stokes／Euler耦合方程的求解。但本文显示分步块LU分解法比经典的Uzawa方法更经济。分析显示该法具有良好的稳定性和高精度,数值结果支持这一理论分析。     

CE/SE法模拟爆震波点火及传播过程

运用时空守恒元和解元(CE/SE)法模拟一维爆震波的形成及传播过程,分析产生爆震波的不同机理和DDT(爆燃-爆震转捩)过程的影响因素.采取激波点火和热点点火两种点火方式,分别对应SDT(激波-爆震转捩)和DDT过程,热点点火初始场采用线性温度分布和阶梯性温度分布,整个管道均为常压.用隐式梯形方法实现刚性源项积分.结果表明,SDT和DDT对应于不同的爆震波机理,并且DDT受点火区不同参数的影响,为研究爆震波点火和DDT过程提供了理论依据.     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法

建立了基于欧拉方程组的二维振荡机翼非定常气动设计反命题方法的数学模型.通过一系列变换,将物理时空中的求解域转换成映象坐标系中的规范区域;并导出了映象坐标下的欧拉方程,结合非定常反命题的边界条件,便可用有限差分方法求解.本方法引入映象坐标系解决反命题边界形状的不确定性,并能利用欧拉方程现有的各种差分格式于反命题求解.     

采用改进的CE/SE方法模拟方管中氢氧爆轰波的稳定传播结构

采用改进的时空守恒元/解元（CE/SE）格式和新型二步化学反应模型,对小扰动情况下方管通道中的氢氧气相爆轰传播过程进行了数值模拟,对不同截面尺寸和两种扰动模式下的爆轰波三波线结构进行了详细计算和分析。数值模拟花费了相对较少的计算代价,清晰地得到了方形截面通道中3种相对稳定的爆轰波传播结构,即直角模式、对角模式和螺旋模式。随着截面尺寸的减小,爆轰波胞格图样与波阵面结构随之产生适应性的变形,乃至出现胞格数目的减少和部分三波线的消失;最终,直角模式和对角模式都会转化为顺时针或逆时针的螺旋模式;临界尺寸下,得到了转换过程中新的过渡结构和压力脉冲振幅的变化,并以此初步讨论了直角和对角传播模式的稳定性。     

基于CE/SE的飞机救生系统动力装置内弹道性能计算

建立某飞机救生系统动力装置套筒式内弹道物理模型,结合CE/SE方法研究粘性二维两相流的内弹道数学模型,计算其性能参数,并与试验数据进行对比。结果表明：套筒内气体压力随时间变化的试验曲线与计算曲线基本吻合,最大压力pm、出口压力pL、工作时间tL和出口速度vL等内弹道性能参数的计算值与试验值最大偏差小于5%。说明运用CE/SE的方法能更加精确的计算内弹道性能,计算具有的准确性和实用性。