高超音速化学非平衡钝劈绕流数值计算

 自1970年Davis提出粘性激波层方法以来,用数值方法求解高超音速轴对称钝体绕流问题国内外已做了大量工作,但至今仍未见到关于平面问题的计算结果。对于平面问题,虽然方程在形式上比轴对称简单,但由于二维效应,激波层较厚,用文献[2]的方法向下游区推进有困难。另外,在驻点线上采用极限关系式虽能克服方程的奇性,但驻点解对流向步长Δξ有依赖。     

湍流燃烧模型对氢燃料超燃室流场模拟的影响

采用化学平衡的假定概率密度函数(PDF)模型和火焰面模型计算了德国宇航研究中心的超燃室反应流,计算结果与有限速率反应模型的和实验的结果进行了对比.使用有限体积法离散Favre平均的N-S方程,湍流模型采用k-ε模型.研究表明:(1)有限速率反应模型在喷氢孔近场,化学平衡的假定PDF模型在喷氢孔远场不能准确捕捉流场的细致结构,而火焰面模型对全流场预测较好,后两种模型的计算时间较有限速率反应模型节省约38%;(2)超燃室内湍流和燃烧相互作用不可忽略,从预测精度和计算效率来看,火焰面模型有较好的工程应用前景.      

回流燃烧室流动特性试验

为了揭示有/无燃烧状态下燃烧室热态和冷态流场的特征和流动特性,针对某型回流燃烧室单头部试验件,使用粒子图像测速仪（PIV）,测量燃烧室燃烧状态下不同截面处的热态流场,以及没有燃烧状态下不同截面处的冷态流场,探讨不同总压损失系数对回流燃烧室热态/冷态流场特征及流动特性的影响。研究表明:随着总压损失系数的增大,冷态条件下各截面流场结构基本保持不变,如射流孔穿透深度、射流角度、回流区位置及大小、流线等基本保持一致,但是各位置点速度大小逐渐增大。热态条件下各截面流场随着总压损失系数增大,流场结构也基本保持不变;相同总压损失系数时,热态流场与冷态流场存在差异,燃油喷射与气流的相对运动将会对燃烧室头部的流场结构造成影响,速度较冷态流动时略微增大。      

回流环形燃烧室反应流场数值模拟

采用自编计算程序,在任意曲线坐标系下对带扩压器、涡流杯、火焰筒、内外环冷却通道与弯管的回流环形燃烧室三维两相化学反应整体流场进行计算.数值预测火焰筒头部结构参数对燃烧流场的影响.所得计算结果与实验数据符合较好,结果表明:采用的数学模型和计算方法合理,该软件可供回流燃烧室优化设计与研制使用.      

模型燃烧室中两相化学反应流动的大涡模拟

对带V型火焰稳定器的模型燃烧室中气液两相化学反应流动进行大涡模拟(LES),研究中采用了欧拉一拉格朗日方法。利用大涡模拟来研究气相流场,k方程亚网格尺度模型来模拟亚网格涡粘性。化学反应速率由亚网格EBU燃烧模型来确定,两步亚网格紊流燃烧模型被用来确定CO浓度,NO亚网格污染物生成模型被用来评估NO生成速率,辐射热量传递采用了热通量辐射模型。并在交错网格体系下气相采用SIMPLE算法和混合差分格式求解,液相在拉格朗日坐标系下采用PSIC算法对其进行求解。通过大涡模拟,能得到气相和液相之间的质量、动量和能量交换。通过计算值和实验值的比较可知,出口温度和污染物成分浓度实验结果与计算结果吻合较好。因此采用大涡模拟方法来模拟两相化学反应流动和污染物生成过程是可行的。     

双组元姿控发动机喷管化学反应流场数值模拟

本文对混合比为０．９、１．０、１．１三种状态下工作的双组元自然推进剂（肼／四氧化二氮）姿控发动机喷管内化学反应流动进行了数值模拟。数值模拟时采用了弱耦合点隐式方法的数值方法及肼／四氧化二氮的十二组分、十三个基元反应的有限速率化学反应模型。得到了三种混合比下反应流及混合比为１．０时冻结流发动机的推力和比推力、喷管中的流动参数及各组分的质量分数。分析表明，数值模拟的结果与理论分析一致，结果可靠。本文工作为姿控发动机的喷管设计提供了理论依据。     

气相爆轰波传播特性的数值模拟及实验对照

本文应用基元反应模型和频散可控耗散格式(DCD)对氢氧爆轰波进行了二维数值模拟.氢氧混合物的化学反应模型考虑了8种组分20个反应方程式.在处理化学反应引起的刚性问题时采用了时间算子分裂的方法.本文首先对爆轰波数值结果和实验结果进行了对照验证,然后对爆轰波在楔面反射由马赫反射向规则反射转变的过程进行了数值分析,得到了反射转变临界角,并和实验结果及理论分析进行了比较,结果是令人满意的;本文还对爆轰波的多波结构进行了初步的数值分析.     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

旋流器后火焰流场的试验研究

 本文给出了用LDA测量旋流器后火焰流场的试验研究结果。试验表明,中心开孔的平底碟形火焰稳定器后方形成双回流区,湍流脉动的均方根速度在回流区边界区域上发生突变,且湍流脉动速度的概率密度分布图上存在双峰;用热电偶测出的温度分布表明内外两支火焰在双回流区的内外边界附近开始,并向下游伸展。试验结果为应用旋流器燃烧室性能改进提供了依据,也为发展带回流的湍流反应流数值计算方法提供了验证依据。     

超燃发动机燃烧室三维化学反应流场的数值模拟

采用ＬＵ隐式方法强耦合求解有限体积法离散的完全N-S方程及组份方程,数值模拟超燃发动机燃烧室内的化学反应流场。计算中采用了AUSM通量分裂格式,两方程及八方程的化学反应模型。计算得到了流场的结构与组份分布,并对两种化学反应模型的结果进行了比较与分析。