入流周期性温度激励下的二维超声速混合层混合增强

对入流施加周期性温度激励下的二维超声速混合层的混合增强进行大涡模拟研究     

基于不同PDF的超声速扩散燃烧火焰面模型对比

为了研究混合分数概率密度函数对湍流扩散燃烧的影响,采用稳态火焰面模型描述超声速扩散燃烧过程,分别采用β-PDF和δ-PDF方法对层流火焰面数据库积分得到两种湍流火焰面数据库,并对比分析了湍流火焰面数据库。结合混合RANS/LES程序,利用DLR氢燃料超燃燃烧室算例进一步对比验证采用不同PDF方法生成的湍流火焰面数据库。研究结果表明,混合分数β分布数据库中间组分随混合分数方差变化大,采用混合分数β分布的算例计算结果脉动速度较大,与实验结果符合得更好。     

基于超声速气流中凹腔主动喷注的强迫点火方案研究

为了进一步研究在入口来流Ma=2.1和T0=846K,p0=0.7MPa的条件下凹腔主动喷注对点火的影响,通过分析高速摄影拍摄的瞬态图像和合成后的统计图像对比了采用凹腔前壁面喷注、后壁面喷注等组合方式的三种乙烯点火方案,并对采用凹腔前壁和后壁面喷注相结合的组合喷注方式开展点火前冷态流场数值仿真分析。研究表明,在此工况条件下,凹腔前壁面喷注容易使燃料在凹腔角回流区内和火花塞处聚集,而凹腔后壁面喷注更加利于点火产生的初始火核在凹腔内传播,但在相同喷前压力下凹腔后壁面的喷注压降要更低;为了使燃料主动喷注后在凹腔内分布更加均匀,凹腔前壁面喷注位置应该设置在靠近剪切层处,而凹腔后壁面喷注位置应该设置在靠近凹腔底壁处;采用凹腔前壁喷注和后壁面喷注相结合的喷注方式相比于单采用前壁喷注和单采用后壁面喷注的喷注方式,凹腔的燃料分布更加均匀,点火后凹腔内火焰更加均匀稳定。当凹腔主动喷注当量比达到0.03、全局当量比达到0.14时,火焰就能穿过剪切层引燃凹腔下游横向射流,在串联凹腔中稳定地燃烧。     

基于平行层推移的含表观裂纹缺陷固体发动机装药燃面计算

将固体火箭发动机装药燃面视为运动界面,利用Level Set界面追踪方法追踪燃面变化过程,计算燃面肉厚曲线;空间、时间离散分别采用五阶高精度Weighted-ENO、三阶TVD Runge-Kutta格式;基于平行层推移规律,对某型号星孔装药进行燃面计算,对比Level Set方法与实体造型法计算结果,二者吻合较好;利...     

基于主动喷注方式的后缘突扩凹腔激光诱导等离子体点火实验研究

针对马赫数2.92,总压2.6MPa和总温1530K的超声速来流条件,基于主动喷注乙烯的燃料喷注方案,在后缘突扩凹腔燃烧室中开展了激光诱导等离子体点火实验研究。通过采集50kHz的CH*基自发辐射图像,详细观测了火焰传播过程并进而研究了主动式燃料喷注方案对点火过程的影响。研究表明,在激光诱导等离子体点火以后,CH*基强度沿时间的变化会呈现出不同的阶段性特征,大体可以分为四个不同的发展变化阶段(激光激发阶段、初始火焰阶段、过渡阶段和全局火焰阶段)。对于采用凹腔主动式燃料喷注方案,都会经历一个类似的初始火焰形成、减弱、增长和迅速发展成全局火焰的过程,区别在于采用凹腔后壁面燃料喷注方案的初始火焰要更加微弱,而且要经历一个难以观测到CH*基信号的初始火焰沉寂阶段。在全局当量比0.03～0.07的条件下,凹腔前壁面燃料喷注方案要比凹腔后壁面燃料喷注方案更加利于初始火焰的发展,但在形成全局火焰以后会引起较大的波动不利于火焰稳定。     

超声速反应流火焰面/进度变量模型

基于Pierce低马赫数条件下的火焰面/进度变量模型,考虑可压缩性和激波的影响,发展了用于超声速湍流反应流的火焰面/进度变量模型.该模型通过求解混合分数和进度变量两个标量方程来描述化学反应,且采用火焰面方程解的完整的三个分支:稳态燃烧分支、具有局部熄火状态的非稳态分支和完全熄火分支.利用该模型对德国宇航中心(DLR)的超燃冲压发动机进行模拟,与稳态火焰面模型得到的附着火焰不同,计算结果捕捉到了熄火和火焰抬举现象;且数值阴影和平均温度分布与实验结果吻合较好.      

一种改进的类DES湍流模拟方法

构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法.      

三维超声速混合层内标量混合的大涡模拟

对三维对流马赫数0.62的超声速混合层流动的标量混合进行大涡模拟, 控制方程对流项采用五阶精度的WENO格式求解, 小尺度涡的作用采用一方程LDKM亚格子模型处理, 过滤后的组分方程中的亚格子组分对流通量采用梯度扩散模化.模拟得到了混合层流场大尺度拟序结构以及标量场的演化过程, 研究表明标量混合过程受混合层内涡系演化所控制, 标量场具有明显的三维特征.模拟得到的速度、组分及其脉动的统计时均结果和实验结果相符较好.      

超声速燃烧室氢气强迫点火过程实验

模拟飞行Ma=4.0的超燃冲压发动机的燃烧室进口条件,在直连式试验台上对超燃冲压发动机氢气火花塞点火过程进行了研究。通过高速摄影与高速纹影捕捉了点火过程初始火焰发展与流程波系演化过程。结果表明,初始火焰传播受高速对流、湍流火焰传播和自点火多种机制控制,凹腔结构不仅提供低速环境,还在初始火焰传播过程中具有能量反馈和质量反馈(活化分子)作用,促进剪切层内火焰的发展。     

基于超声速边界层转捩的WENO格式对比分析

基于高精度weighted essentially non-oscillation(WENO)格式对马赫数为4.5,雷诺数为10000的超声速来流条件下平板边界层转捩过程进行了大涡模拟.无黏通量分别用5阶、7阶、9阶WENO格式进行离散,黏性通量离散采用4阶中心差分格式,时间推进采用具有总变差递减(TVD)性质的3阶精度Runge-Kutta方法.通过在入口边界叠加一对等幅最不稳定第一模态谐波扰动,分别采用3种WENO格式计算得到了平板层流边界层失稳转捩的演化过程.结果表明:5阶WENO格式的数值耗散明显高于7阶和9阶WENO格式,在捕捉湍流涡和流场脉动特性上存在明显不足.9阶WENO格式的耗散小,能够捕捉到流场更小尺度的涡和高频脉动.研究具有高频脉动特性的问题或者欲捕捉精细涡结构时,建议采用7阶以上的高精度格式.