跨音速喷管化学反应流数值模拟

用时间相关的半隐格式有限差分数值方法求解了化学非平衡反应跨音速喷管流场，在喷管收敛段，流动接近化学平衡状态，控制方程的刚性问题严重，数值积分困难。通过对时间差分项隐式离散、对空间差分项显式离散，流场边界采用参考平面上的特征线法计算等，成功地解决了由于化学反应有限速率带来的数值解不稳定问题。该格式简单、需要计算机存贮空间少。本文完成了一维和轴对称非平衡化学反应喷管流动计算，并与化学平衡流和冻结流的计算结果做了比较。     

斜切反喷管性能计算

本文提出了斜切反喷管流场计算方法,内部点用MacCormack二步显格式数值解控制方程,入口边界和固壁边界的参数用物理边界条件和可应用的沿双特征线的相容性方程计算.固壁边界计算中采用了激波装配技术和开关格式,有效地捕获了流场内的激波系列,计算中得到的激波情况和物理分析一致.本文在流场计算的基础上,进行了反喷管的性能计算.本文提出的方法可应用于反喷管的性能预估.     

火箭发动机塞式喷管流场的数值研究

用时间相关法和张涵信的无振荡、无自由参数ＮＮＤ显式格式,对火箭发动机塞式喷管二维轴对称定常流场进行了数值模拟。计算得到了有粘和无粘中两种情况下,喷管流场的马赫数压强、温度、流线以及速度矢量分布图,计算网络采用代数方法生成,并通过求解Ｌａｐａｃｅ方程对网格进行修正。计算结果表明：ＮＮＤ格式用于计算塞式喷管流场具有较高的效率和良好的分辨率;无粘和有粘情况得到的流场结构有较大的差别;基础喷管的出口马赫数     

粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

外流对塞式喷管性能的影响

为了了解外流对塞式喷管性能的影响。本文从N-S方程出发，采用NND格式对塞式喷管有外流的流场进行了数值模拟，重点研究了外流马赫数和外流迎角对塞式喷管的影响，研究表明，在一定飞行高度下，外流会使塞式喷管的性能降低，并且外流为亚声速时对塞式喷管的影响较小，在跨声速和超声速时影响较大；当外流有迎角时，会改变塞式喷管推力矢量，正迎角会使塞式喷管性能得到补偿，负迎角则会降低塞式喷管的性能，随着背压减小，外流对塞式喷管的影响也会减小。     

固体火箭发动机喷管喉衬动边界条件下的温度场计算

本文探讨了固体火箭发动机喷管喉村在烧蚀退移、燃气中固体微粒的沉积以及沉积物吹除所引起的动边界条件下的温度场数值计算模型。模型的建立以有限体积能量平衡的基本原理为依据,采用有限差分的数值计算方法及稳定性较好的交替方向隐式格式求解。为确保数值计算结果在边界退移及网格脱落时的稳定性,采用了插值与迭代相结合的方法,以全隐式格式求解表层网格节点和受热壁面的温度。模型还考虑了多层背村材料结构、变物性材料特性等诸种因素,使之具有一定的普遍性,并可应用于工程设计和研究中。数次发动机实验证明,本文提出的温度场计算模型能够较真实地反映喉衬温度响应,具有较为满意的计算精度。     

硅基内衬喷管扩张段烧蚀和温度场的耦合计算方法

本文叙述固体火箭喷管扩张段硅基内衬烧蚀和温度场的耦合计算方法.烧蚀采用非定常的液体层烧蚀模型;温度场计算采用圆柱坐标系下的一维径向瞬态导热模型,考虑材料在高温下形成的炭化层和热解层,而由烧蚀产生的边界移动问题用坐标变换的方法处理.利用烧蚀和导热计算方程组中二者的结合部,通过充分耦合,获得了更准确的预示计算结果.     

气动塞式喷管底部二次流特性的数值模拟

发展了一套可以对底部有二次流的气动塞式喷管的性能进行预示的数值方法，采用二阶精度的NND格式求解二维层流N-S方程，对底部二次流的入口边界条件进行了特殊处理。以一种实验用的直排式气动塞式喷管为对象，针对两种不同的背压条件，对二次流流量分别为主流的0%，1%，3%和5%的几种工况研究了塞式喷管底部二次流流量变化对喷管性能影响。计算结果表明，底部二次流的加入使得气动塞式喷管性能有比较明显的提高。数值方法可以用于气动塞式喷管的设计和性能预报。     

固体火箭发动机喷管结构瞬态温度场的理论预估

本文采用目前广泛应用、效果显著的有限元法,对固体火箭发动机复合喷管结构的轴对称瞬态温度场进行了理论预估.在计算中考虑了对流换热系数沿喷管轴向的变化及材料的方向性.所编程序通用性较好.     

塞式喷管流场数值模拟

采用ＮＮＤ２Ｍ差分格式,从考虑к－ε湍流模型的理想气体三维薄层近似Ｎ－Ｓ方程出发,对塞式喷管内外流场进行了数值模拟。计算结果与国外文献十分一致。通过流场分析,对塞式喷管流动机理有了清楚的认识。计算表明ＮＮＤ２Ｍ差分格式及相应计算软件适合于模拟复杂喷流流场,可用于塞式喷管发动机的研究。     

SLAU格式在大涡模拟方法中的应用

采用格心格式的有限体积法求解三维大涡模拟控制方程,成功地将最新改进的SLAU格式推广到大涡数值模拟方法中。SLAU格式的空间重构方式采用五阶WENO格式,时间推进采用LU-SGS隐式时间推进。湍流模式为最新提出的基于随机分析的动态亚格子模型。计算了圆柱绕流、NACA4412翼型绕流、NACA0012翼型绕流三个算例,分别从可压缩较大迎角翼型绕流、可压缩跨音速翼型绕流以及不可压缩低马赫数圆柱绕流三个方面验证了将SLAU格式应用到湍流大涡数值模拟计算中的可行性。结果表明,结合SLAU格式的大涡模拟方法可以达到较好的数值模拟效果,同时间接检验了基于随机分析的动态亚格子模型的数值模拟效果。     

跨音速大迎角Euler方程数值分析

应用有限体积法离散三维Euler方程,选用三步显式格式进行时间推进求解;并通过当地时间步长、残值光顺、焓修正等方法加速收敛。对一维波动方程的von Neumann稳定性分析说明这种三步格式的最大Courant数是2。对三角翼跨音速大迎角流动的数值分析表明文中的三步格式及焓修正方法优于人们经常采用的四步、五步格式以及焓阻尼技术。     

求解可压缩流动的低耗散E-CUSP迎风格式性能分析

采用对比分析的方法,对求解可压缩流动的低耗散E-CUSP(LDE)迎风格式的性能进行测评,目的是进一步更好地将其用于分离涡数值模拟。对比格式采用广泛使用的Roe格式。控制方程为Euler和雷诺平均Navier-Stokes方程耦合Spalart-Allmaras一方程湍流模型,采用全隐式无分裂Gauss-Seidel线松驰时间推进方法。结果表明,LDE格式具有强的捕捉激波能力,计算结果与Roe格式一致,和实验数据吻合良好。相对于Roe格式,LDE格式具有较快的收敛速度。     

圆截面U形管流有限元模拟

弯曲圆管内流动具有很强的工业应用背景,管内二次涡的大小、位置及发生的频率会对流动系统的阻力及热质交换产生重大作用,研究弯曲圆管内的流动具有重要的理论意义.本文应用湍流大涡模拟方法对弯道流动进行了数值模拟.用时间推进法,给出了隐含流线迎风耗散作用的二阶完全展开Taylor-Galerkin有限元离散格式,数值计算的结果与实验结果符合较好.表明大涡模拟方法适用于边界形状复杂,存在各向异性的大尺度涡的流动仿真.     

小能量点火脉冲爆震发动机工作过程数值模拟

为完善脉冲爆震发动机数值模拟模型,采用有限速率化学反应模型,通过二阶迎风格式差分逼近二维欧拉方程,模拟以煤油和氧气为可爆混合物的小能量点火触发爆震的脉冲爆震发动机工作过程。采用CFD软件中spark ignition模型模拟实验中电火花塞点火,实现了50mJ小能量点火触发爆震。数值模拟结果和试验结果吻合,说明所采用的计...     

高超声速边界层感受性问题高精度数值模拟

通过五阶WENO格式和六阶对称紧致格式以及三阶TVD R-K结合的方法,对存在强激波和小扰动相互干扰的高超声速边界层感受性问题进行研究.结果表明:此方法能够模拟边界层内不稳定波的产生和发展以及小扰动和激波相互干扰等现象,因此能广泛应用于含激波的感受性等问题的可压缩湍流直接数值模拟以及具有激波和边界层干扰等复杂流场的计算.     

锯齿形转捩片触发高超声速进气道边界层转捩的大涡模拟

为促进锯齿形转捩片在高超声速进气道中的应用,以地面风洞条件下的二元进气道为研究对象,采用高精度大涡模拟方法对锯齿形转捩片在三级压缩楔面上触发的边界层转捩现象开展了研究。数值方法基于隐式亚格子模型,空间离散采用高精度通量限制型紧致格式,时间推进采用显式Runge-Kutta方法。数值模拟清晰捕捉到了边界层转捩的空间发展演化过程,并获得了统计平均流场以及流场脉动特征。数值模拟结果表明转捩片能够有效触发进气道压缩面边界层转捩;通过与等熵压缩面及单楔面数值模拟结果的对比分析,获得了转捩片触发边界层转捩的内在机理,为后续研究工作奠定了基础。     

航空发动机滑油通风系统性能计算仿真

采用自主研发的航空发动机滑油通风系统性能分析软件,对某型发动机通风流路进行了数值仿真。数值仿真以流动换热的网络算法为基础,经过结构分析,构建了通风流路仿真元件与节点组成的网络仿真模型。在给定实验条件下,数值仿真了稳态情况下通风流路腔室节点的压力,各密封装置泄漏量的分配。最终的仿真结果与实验参考值的误差在l％以内,验证了数值仿真方案的精度。     

用虚拟黏性法构造Navier-Stokes方程黏性项的隐式格式

提出的虚拟黏性法是一种关于Navier-Stokes方程黏性项(黏性应力项和传热项)的隐式新方法。通过引入虚拟时间和虚拟黏性项,将隐式格式的构造大大简化,从而避免了大型复杂隐式差分方程组的常规求解。在虚拟时间推进过程中所需求解的方程组的系数矩阵是一个三对角矩阵,它具有计算简单且计算量小的优点。用模型方程和Navier-Stokes方程进行了数值仿真,研究各种参数对计算的影响,并在精度和效率上与显式方法进行比较,证实了算法的正确性和优势。最后分析了本文格式的适用范围。