PSP measurement of flow fields in a supersonic combustor with a backward-facing step

The mixing fields within a SCRAM-jet combustion chamber are visualized using pressuresensitive paint (PSP) as an oxygen sensor. The experiments are performed in a small supersonic wind tunnel at the National Aerospace Laboratory — Kakuda Research Center (NAL-KRC). The main stream Mach number is 2.4, and the dynamic pressure ratios between the injected gas and the main flow are 0.3, 0.7, 1.1 and 1.5. Three fuel injection nozzles are used; oxygen is injected from the central nozzle and air from the two nozzles at either side. The spread of the injected gas is measured to observe the effects of placing the nozzles in different positions. The results show that the jet has its own independent flow structure, and that little mixing of gases occurs between the flow structures created by each nozzle. When the injection dynamic pressure ratio is increased, the oxygen fraction rises in the recirculation zone and falls in the separation zone downstream of the injection.     

Numerical simulations of shock-wave interaction with a boundary layer in the plane supersonic flows with jet injection

A supersonic air flow in a plane channel with a transverse turbulent jet of hydrogen injected through a slot on the bottom wall is simulated. The algorithm for solving the Favre-averaged Navier-Stokes equations for the flow of a perfect multispecies gas on the basis of the WENO scheme is proposed. The main attention is paid to the interaction of the shock-wave structure with the boundary layers on the upper and lower duct walls under the conditions of an internal turbulent flow. Namely, a detailed study of the structure of the flow is done, and separation and mixing depending on the jet slot width are investigated. It is found that in addition to well-known shock-wave structures produced by the interaction of the free stream with the transverse jet and the bow shock interaction with the boundary layers near the walls, an additional system of shock waves and the flow separation arise on the bottom wall downstream at some distance from the jet. The comparison with the experimental data is performed.     

反应性RM不稳定混合特性研究: 反应活性与介质不均匀性的影响

预混火焰界面的RM（Richtmyer-Meshkov）不稳定现象在自然界和工程实践中十分常见,但目前关于反应性RM不稳定的研究主要集中于均匀介质的情况,而实际中的预混气体往往是非均匀的,因此开展非均匀介质中火焰界面演化和混合特性的研究十分必要。采用带单步化学反应的Navier-Stokes方程和高精度数值格式,研究了预混火焰界面在入射激波及反射激波作用下的RM不稳定过程,考察了化学反应活性以及介质非均匀性对RM不稳定过程中火焰界面混合特性的变化规律的影响。结果表明,在入射激波作用后的阶段,在均匀介质中的火焰界面形态呈现典型的"钉-帽-泡"结构,化学反应活性越强,界面的"泡"结构和"钉-帽"结构增长越快;而在非均匀介质中,火焰界面形态则呈现"钉-钉"结构,界面在流向速度差的诱导下被更大程度地拉伸。在第一次反射激波作用后的阶段,混合区的增长速率不依赖于反应活性和均匀性,仅与流动特性有关。时间尺度的研究表明,大尺度流动是反应性RM不稳定的主导因素,其次是化学反应,最后是小尺度混合,化学反应的强化会抑制大尺度流动,非均匀性会强化大尺度流动。      

Flow and mixing fields of turbulent bluff-body jets and flames

The mean structure of turbulent bluff-body jets and flames is presented. Measurements of the flow and mixing fields are compared with predictions made using standard turbulence models. It is found that two vortices exist in the recirculation zone; an outer vortex close to the air coflow and an inner vortex between the outer vortex and the jet. The inner vortex is found to shift downstream with increasing jet momentum flux relative to the coflow momentum flux and gradually loses its circulation pattern. The momentum flux ratio of the jet to the coflow in isothermal flows is found to be the only scaling parameter for the flow field structure. Three mixing layers are identified in the recirculation zone. Numerical simulations using the standard k-? and Reynolds stress turbulence models underpredict the length of the recirculation zone. A simple modification to the C₁ constant in the dissipation transport equation fixes this deficiency and gives better predictions of the flow and mixing fields. The mixed-is-burnt combustion model is found to be adequate for simulating the temperature and mixing field in the recirculation zone of the bluff-body flames.     

激波作用不同椭圆氦气柱过程中流动混合研究

在激波与气柱相互作用问题中,压力与密度间断不平行产生的斜压涡量会引起流动的不稳定性,从而促进物质间的混合.本文基于双通量模型,结合五阶加权基本无振荡(WENO)格式,求解多组分二维Navier-Stokes方程,分析激波作用面积相同结构不同的椭圆气柱所致的流动和混合.数值结果清晰地显示了激波诱导Richtmyer-Meshkov不稳定性引起的气柱界面变形和波系演化.同时定量地从界面运动、界面结构参数变化(长度和高度)、气柱体积压缩率、环量及混合率等角度分析激波诱导的流动混合机制,研究椭圆几何构型对氦气混合过程的影响.结果表明,界面及相关参数的演化与气柱初始形状密切相关.当激波沿椭圆长轴作用于气柱时,气柱前端出现空气射流结构,且射流不断增长并渗透到下游界面,致使气柱分离成两个独立涡团,离心率越大,射流发展越快;同时激波作用气柱后在界面处产生不规则反射现象.圆形气柱界面演化与这种作用情形类似.当激波沿椭圆短轴作用于气柱时,界面上游出现类平面结构,随后平面上下缘处产生涡旋,主导流动发展,激波在界面作用产生规则反射,离心率越大,这些现象越明显.界面高度、长度、体积压缩率也因此有所差异.对界面演化、环量和混合率的综合分析表明,激波沿长轴作用于气柱且离心率较大时,流动发展较快,不稳定性导致的流动越复杂,越有利于氦气与环境介质的混合.     

Three-dimensional reactive shock bifurcations

We model interactions of a premixed flame with incident and reflected shocks in a rectangular shock tube using three-dimensional (3D) reactive Navier–Stokes numerical simulations. Shock-flame interactions occur in the presence of boundary layers that cause the reflected shock to bifurcate and form a reactive shock bifurcation (RSB), which contains a flame in the recirculation zone behind the oblique shock. The recirculation zone acts as a flame holder thus attaching the flame to the shock in the vicinity of the wall, and providing a mechanism for a detonationless supersonic flame spread. The accelerated burning induced by an RSB, and Mach stems that may result from RSB–RSB interactions, promote hot-spot formation, and eventually accelerate deflagration-to-detonation transition. Schlieren-type images generated from the simulation results show that the 3D structure of an RSB may not always be easily recognized in experiments if the RSB is attached to the surface of the observation window. The main 3D effect observed in the simulations is caused by the presence of the second no-slip wall in a 3D rectangular channel. Two RSBs that form at adjacent walls interact with each other and produce an oblique Mach stem between two oblique shocks. The oblique Mach stems then interacts with a central Mach stem that forms near symmetry plane, and this interaction creates a hot-spot that leads to a detonation initiation.     

2D mixture fraction measurements in a high pressure and high temperature combustion system using NO tracer-LIF

Mixture fraction measurements in a jet-in-cross flow configuration at high pressures (15 bar) and temperatures (above 1000 K) were performed using planar laser induced fluorescence of nitric oxide (NO-PLIF) as trace species. The goal was the evaluation of this tracer LIF technique for the characterization of the mixing of fuel and hot exhaust gas in the mixing channel. The fuel (natural gas (NG) or H₂/N₂/NG mixture) along with the tracer were injected into the crossflow of the exhaust gas and PLIF measurements were performed in different planes. In order to relate the measured NO-LIF signal to fuel concentration and mixture fraction, effects of pressure, temperature and species concentration were taken into account. Numerical calculations and spectroscopic simulations that mimic the experimental conditions were performed to identify excitation schemes that give optimum correlations between the NO-LIF signal and the mixture fraction. The measured NO-PLIF images were transformed into mixture fraction plots using the computed correlations. The paper reports on the experimental challenges encountered during the measurements and the steps taken to overcome those difficulties. Examples of mixture fraction distributions are presented and discussed. The paper concludes with a detailed analysis on the accuracy of the measured mixture fraction values.     

Investigation on the flameholding mechanisms in supersonic flows: backward-facing step and cavity flameholder

Abstract  

As effective devices to extend the fuel residence time in supersonic flow and prolong the duration time for hypersonic vehicles cruising in the near-space with power, the backward-facing step and the cavity are widely employed in hypersonic airbreathing propulsive systems as flameholders. The two-dimensional coupled implicit RANS equations, the standard k-ε turbulence model, and the finite-rate/eddy-dissipation reaction model have been used to generate the flow field structures in the scramjet combustors with the backward-facing step and the cavity flameholders. The flameholding mechanism in the combustor has been investigated by comparing the flow field in the corner region of the backward-facing step with that around the cavity flameholder. The obtained results show that the numerical simulation results are in good agreement with the experimental data, and the different grid scales make only a slight difference to the numerical results. The vortices formed in the corner region of the backward-facing step, in the cavity and upstream of the fuel injector make a large difference to the enhancement of the mixing between the fuel and the free airstream, and they can prolong the residence time of the mixture and improve the combustion efficiency in the supersonic flow. The size of the recirculation zone in the scramjet combustor partially depends on the distance between the injection and the leading edge of the cavity. Further, the shock waves in the scramjet combustor with the cavity flameholder are much stronger than those that occur in the scramjet combustor with the backward-facing step, and this causes a large increase in the static pressure along the walls of the combustor.     

激波与层流/湍流边界层相互作用实验研究

在超声速风洞中,分别对层流和湍流来流条件下的边界层和斜激波(激波强度足以引起流动分离)相互干扰进行了实验研究,利用纳米粒子示踪平面激光散射(NPLS)技术获得了两种条件下流场的精细结构图像;利用粒子图像测速(PIV)技术获得了两种条件下流场的速度场和涡量场;综合运用NPLS结果和PIV结果对比分析了两种流动的瞬时流动结构和时间相关性,实验结果表明:层流边界层内的分离区呈现出狭长的条状,而湍流边界层内分离区呈现出较规则的椭圆;在入射激波上游距入射点较远的位置,层流边界层外围拟序结构会诱导出一系列压缩波系,进而汇聚成空间位置不稳定的诱导激波,而湍流边界层则是在入射激波上游较近的地方直接形成较强且稳定的诱导激波;在入射激波下游,层流边界层内的膨胀区域较小且急促,膨胀后产生的再附激波很弱,而湍流边界层内的膨胀区域较大,膨胀后产生的激波较强。     

Direct Numerical Simulation of Three-Dimensional Richtmyer--Meshkov Instability

Direct numerical simulation （DNS） is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interracial instability is known as the Richtmyer-Meshkov （R-M） instability. The compressible Navier-Stoke equations are discretized with group velocity control （GVC） modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interracial instability on produced vorticity and mixing is discussed.     

激波冲击火焰的流动拓扑研究

为深入研究激波冲击火焰现象的内在机制,采用二维带化学反应的Navier-Stokes方程对现象进行数值研究,通过对速度梯度张量特征方程的分析证明Okubo-Weiss函数适用于可压缩流动,并重点分析火焰区的流动拓扑特性.结果表明,波后火焰区内Okubo-Weiss函数积分量基本守恒,但在火焰区内部和表面具有截然不同的流动状态,且火焰发展基本不受流场可压缩性的影响;波后火焰区的流动拓扑分类主要以焦点和鞍点为主,意味着流场中变形占主导.     

二维拉格朗日坐标系下气粒混合双向耦合对激波流场影响的计算

喷射颗粒与气体混合是内爆压缩领域的热点和难点. 针对喷射混合中的气粒双向耦合问题, 开展了理论建模、离散算法以及颗粒反馈对激波流场的影响研究. 建立了拉格朗日计算框架下的数学模型; 给出了耦合源项的离散算法; 开展了平面及汇聚构型条件下, 气粒双向耦合的数值模拟研究; 发现了颗粒反馈导致气体激波提速现象以及气区流场物理量分布形态的改变, 初步获得了量化分析结果. 本文建立的数学模型、计算方法和获得的新的物理认识, 为深入理解喷射混合现象、解决相关工程应用问题提供了重要理论支撑.     

激波与SF_6梯形气柱相互作用的数值模拟

基于大涡模拟,结合五阶加权基本无振荡格式以及沉浸边界法对平面入射激波与两种SF_6梯形重气柱的相互作用过程进行了数值模拟,数值结果清晰地显示了激波诱导Richtmyer-Meshkov不稳定性所导致的两种梯形重气柱的变形过程,详细分析了入射激波在两种梯形重气柱界面发生反射、折射、绕射以及折射激波与透射激波在气柱内部来回反射的过程,并研究了该过程中所产生的复杂波系结构,对两种梯形气柱变形过程中与周围空气的混合过程进行了分析;通过记录气柱界面四个特征尺寸随时间的变化对两种梯形气柱界面的不同演化过程进行了定量分析。     

受限空间内旋流回流区特性

为了研究受限空间内旋流回流区的三维结构特性,采用realizable k-ε模型模拟了旋流数等于0.884时,不同受限空间内的旋流流场.受限率是影响回流区形态的重要因素,受限率大于6时,中心回流区与下游回流区是两个独立的区域,有两对涡结构;受限率在3～6之间时,中心回流区与下游回流区合并到一起,存在两对独立的涡结构;受...     

激波作用下SF6气泡界面演化和射流发展的数值模拟

 数值研究了平面激波冲击氮气环境中SF₆气泡界面的Richtmyer-Meshkov不稳定性,重点关注其中的激波聚焦及射流的产生和发展过程。在入射激波马赫数为1.23的情况下,给出了压力、密度、数值纹影和涡量等物理量的演化图像,定量分析了流场中压力最大值、密度最大值、射流速度、环量和斜压力矩随时间的变化关系。计算结果表明,平面激波冲击SF₆气泡过程有很强的聚能效应,在气泡内部靠近下游极点处发生激波近似理想聚焦和点爆炸现象,直接导致出现二次波系以及向下游运动的细长射流结构。相比入射激波,二次波系产生斜压力矩和涡量的能力要弱得多。     

带喷流超声速光学头罩流场气动光学畸变试验研究

超声速光学头罩在大气层内飞行时, 需要在光学窗口表面顺来流方向进行喷流冷却, 致使窗口上方流场更为复杂. 目标光线穿过窗口上方流场, 受到激波、膨胀波、混合层、湍流边界层等流场结构引起的变密度场影响而产生波前畸变, 导致成像出现偏移、抖动、模糊等气动光学效应. 本文对马赫数3.8来流条件下有无喷流时超声速光学头罩流场引起的气动光学波前畸变进行了试验研究. 基于纳米示踪的平面激光散射技术, 首先对流场图像进行密度校准获得高时空分辨率密度场, 然后采用光线追迹法计算得到波长532 nm平面光波垂直于光学窗口穿过流场后的光程差(optical path difference, OPD)分布, 并对窗口上方近壁区有无喷流状态的流场结构引起的 OPD分布进行了研究. 发现无喷流时, 流场结构相对较为简单, 窗口上方有较长的回流区和层流区, 而有喷流时窗口上方出现复杂的剪切层、混合层及湍流边界层, 流动很快就转捩为湍流结构, 其引起的气动光学畸变要明显高于无喷流状态. 无喷流状态相隔5 μs的流场引起的光程差均方根值分别为0.0348 和0.0356 μm, 有喷流状态的光程差均方根值分别为0.0462 和0.0485 μm. 关键词： 超声速 喷流 气动光学 光程差     

水下超声速气体射流气液两相复杂流动研究

本文实验研究了水下超声速气体射流气液两相复杂流动。利用高速摄影仪和电子相机分别实时记录了过膨胀超声速工况水下气体射流的喷射状态和整体形貌,显示了不同工况水下高速气体垂直射流的演化过程和动态不稳定性形貌。研究结果表明:在射流的初始段存在与射流内部复杂波系相关的激波反馈特性,激波反馈特性发生之前存在能量积聚的高频低幅的胀鼓过程,二者均随机发生;在射流的主体段,在气水掺混和卷吸大规模能量交换作用下,射流呈现随机的偏摆效应,并且偏摆受环境流场影响明显。     

在冷气干扰下叶栅绕流的控制方程及求解方法

1控制方程高温涡轮的冷却空气量随着燃气轮机性能的提高而在不断增加，因而在数值模拟气冷涡轮流场时则须考虑喷射冷气对燃气主流的影响。本文采用文献［1］的“质量源”模型，推导出具有统一形式的适用于叶栅三维和周向平均混合流动的控制方程。在绝热和忽略粘性应力作功率条件下，将文献出中的通用形式的控制方程在柱坐标系中展开，然后再经坐标变换，将其转换到任意曲线坐标系《，刀J中，得到了在冷气干扰下的叶栅三维绕流的双曲守恒型控制方程（用矩阵形式表示）：aU／＆＋aE／a＋OF／的十OG／0（一H（1）式中仅三项含有冷气参数，…     

Numerical simulations of the flow field ahead of an accelerating flame in an obstructed channel

The development of the unburned gas flow field ahead of a flame front in an obstructed channel was investigated using large eddy simulation (LES). The standard Smagorinsky–Lilly and dynamic Smagorinsky–Lilly subgrid models were used in these simulations. The geometry is essentially two-dimensional. The fence-type obstacles were placed on the top and bottom surfaces of a square cross-section channel, equally spaced along the channel length at the channel height. The laminar rollup of a vortex downstream of each obstacle, transition to turbulence, and growth of a recirculation zone between consecutive obstacles were observed in the simulations. By restricting the simulations to the early stages of the flame acceleration and by varying the domain width and domain length, the three-dimensionality of the vortex rollup process was investigated. It was found that initially the rollup process was two-dimensional and unaffected by the domain length and width. As the recirculation zone grew to fill the streamwise gap between obstacles, the length and width of the computational domain started to affect the simulation results. Three-dimensional flow structures formed within the shear layer, which was generated near the obstacle tips, and the core flow was affected by large-scale turbulence. The simulation predictions were compared to experimental schlieren images of the convection of helium tracer. The development of recirculation zones resulted in the formation of contraction and expansion regions near the obstacles, which significantly affected the centerline gas velocity. Oscillations in the centerline unburned gas velocity were found to be the dominate cause for the experimentally observed early flame-tip velocity oscillations. At later simulation times, regular oscillations in the unburned streamwise gas velocity were not observed, which is contrary to the experimental evidence. This suggests that fluctuations in the burning rate might be the source of the late flame-tip velocity oscillations. The effect of the obstacle blockage ratio (BR) on the development of the unburned gas flow field was also investigated by varying the obstacle height. Simulation predictions show favorable agreement with the experimental results and indicate that turbulence production increases with increasing obstacle BR.     

Numerical study on spinning detonations

Spinning detonations in both a circular tube and a square tube are presented in order to reveal characteristics of spinning modes by using three-dimensional simulations with a detailed chemical reaction model. The present results show a feature of a single spinning detonation which was discovered in 1926. The shock patterns in both cases are similar except the pressure trail, however, the shock wave angles and the shock wave lengths are shown to be dependent on the cross section configuration of the tube. The pitch angle, the track angle, the Mach stem angle, and the incident shock angle on the tube wall in the numerical results agree well with those in the experimental ones, and they are independent of the compositions of mixture, tube diameters, and initial pressures.