可压缩涡环和平面激波相互作用的数值研究

本文从非定常的Ｅｕｌｅｒ方程出发,进行了可压缩涡环和平面激波相互作用的数值研究。     

超燃冲压发动机不同燃烧室结构流场数值模拟

对带不同长深比凹腔结构（燃烧室）的二维超燃冲压发动机流场进行数值模拟,采用迎风二阶精度NND格式求解二维欧拉方程,对流场进行冷态数值计算,能够更加清楚地捕捉、分析激波的特点,反映流场的状态和旋流特点;结果证明凹腔可以使流场中形成低速高压回流区,以达到燃料混合燃烧的目的。凹腔长深比对流场特性、总压恢复、燃烧室阻力影响显著。随着长深比的越大,凹腔旋流强度越强,总压损失越大,燃烧室阻力越大。     

超声速弱欠膨胀冲击射流流场结构

为了对超声速弱欠膨胀冲击射流的流场结构细节进行研究,使用大涡模拟方法对其进行了数值模拟。利用三阶迎风和四阶对称紧致格式对无量纲化轴对称可压缩滤波N-S方程进行空间离散,时间上推进采用的是三阶精度的TVD型Rugge-kutta法。亚格子尺度模型采用的是修正Sm agorinsky涡粘性模型。通过与经典的冲击射流实验比较,证明了程序的可靠性。数值模拟得到了剪切层以及壁面射流中的涡结构和主射流中的激波结构,并且在此基础上对涡合并和板前激波和涡干扰现象进行了深入研究。发现涡合并现象主要出现在流场的上游,越往下游出现的几率越小;涡和板前激波的相互作用会引起激波位置和强度以及冲击平板上冲击区的压强的显著变化,同时也会导致涡的变形。     

非平衡流解耦方法及其计算激波诱导燃烧的应用验证

介绍并发展了一种改进的化学非平衡流动解耦方法.将流动控制方程与化学反应生成源项解耦处理,组分对流项的矢通量分裂基于流动方程密度对流项的分裂形式给出.对流项采用五阶WENO格式求解,化学反应源项的刚性采用简化的隐式方法进行处理.将该方法应用于预混H₂/O₂和H₂/Air激波诱导燃烧计算,得到的定常及非定常燃烧过程与实验观测相符,振荡频率与实验测量较吻合,结果表明该方法应用于多组分多步反应流计算是可行的.      

激波管在酚醛树脂高温热解动力学研究中的应用

使用激波管作为加热手段,利用其加热速率快的优点,突破了传统方法在加热速率上的限制,研究了酚醛树脂在1400~1700K温度范围内的热解动力学.主要碳氢产物是甲烷、乙烯、乙炔和苯.通过对反应扩散过程的分析,考察了扩散对热解过程的影响.结果表明,实验中酚醛树脂的反应扩散过程迅速达到稳态,扩散影响可以忽略,首次获得了酚醛树脂在芳环开环热解机制下的热解速率常数.     

超声速风扇叶型设计研究

为了设计低流动损失超声风扇叶型,研制出基于数值最优化超声叶型设计软件.以叶型主要几何参数为叶型设计参数,针对轴向进气第一级风扇转子采用给定轮缘功、损失最小为目标函数.设计的优化叶型后面大部分呈平直因此气流转角小、减速增压完全通过激波实现.沿弦向叶片厚度逐渐增加使叶栅通道呈收缩,总体上使流过叶栅超声气流得到减速增压.优化叶型所构成叶栅流场具有多道激波、流动损失低.设计结果表明设计方法的可行性,并可为超声风扇叶型设计提供新思路.      

振动叶栅中尾迹/激波相干的数值研究

研究了尾迹/激波相干对振动叶栅非定常流动的影响.数值方法利用了先进的谐波方法.该方法基于傅立叶变换,把叶栅非定常流变量分解为一时间平均量和若干谐波分量,继而,时域非定常N-S方程被转化为频域定常傅立叶系数方程.定常傅立叶系数方程的求解采用了伪时间推进方法.为加速数值解收敛,采用了本地时间步长和多重网格技术.通过NASA67转子算例结果表明:上游尾迹对振动叶栅非定常流动具有显著影响.当存在尾迹相干时,叶片通道内的激波在轴向方向上表现出较大的漂移运动,叶片所受的非定常力幅值变大,同时失去了时间周期性.      

一种缩短隔离段简易措施的实验研究

为了提高超燃冲压发动机隔离段耐反压能力以及缩短隔离段长度,设计了一种带等宽度平直斜楔的隔离段,斜楔放置在隔离段进口下壁面厚附面层一侧,在M1.98非对称的隔离段进口来流条件下完成吹风实验。实验结果表明,隔离段加等宽度平直斜楔后可在同样的反压下使激波串相对长度从13.87缩短至11.12,缩短了2.75。同样的隔离段加斜楔后能够承受的最大反压从来流静压的3.55倍上升到3.85倍,提高了8.45%左右。     

激波诱导燃烧模拟中轴对称边界数值异常研究

对钝体激波诱导燃烧模拟中引起对称轴附近出现数值异常现象的几个影响因素进行了细致讨论.在飞行速度低于爆轰速度情形的模拟中,在一定条件下靠近对称轴的位置会出现明显的数值异常——燃烧阵面扭曲.分别从边界处理方法、计算格式、化学反应机制和计算网格等4个方面对该异常产生的原因进行了详细地讨论.研究表明4个因素都会导致数值异常产生,深入的分析表明计算格式是其中最根本的因素,而这种燃烧阵面扭曲的数值异常应该与粉刺问题(Carbuncle)类似.针对如何克服轴对称边界数值异常提出了一些合理的措施,并给出了验证算例的对比与分析.      

非对称交叉激波和湍流边界层相互作用的数值研究

针对7°×11°双尖鳍外形的非对称交叉激波与湍流边界层相互作用,采用Navier-Stokes方程和5种湍流模型进行了计算.主要考察对壁面压强、热传导、绝热壁面温度和壁面摩擦力线分布的计算精度.计算结果表明:两道斜激波相交后的区域的壁面压强和热传导都比较高;计算的压强和壁面摩擦力线与试验吻合很好,绝热壁面温度次之,热传导最差,峰值高达试验的3倍左右.湍流模型对壁面压强和壁面摩擦力线影响很小,对绝热壁面温度和热传导影响很大.在5种湍流模型中,TNT(turbulent/non-turbulent)和SST(shear-stress transport)模型表现较好.