Experimental Investigation of Nozzle Exit Reflector Effect on Supersonic Jet

The present study describes an experimental work to investigate the effect of a nozzle exit reflector on a supersonic jet that is discharged from a convergent–divergent nozzle with a design Mach number of 2.0. An annular reflector is installed at the nozzle exit and its diameter is varied. A high-quality spark schlieren optical system is used to visualize detailed jet structures with and without the reflector. Impact pressure measurement using a pitot probe is also carried out to quantify the reflector’s effect on the supersonic jet which is in the range from an over-expanded to a moderately under-expanded state. The results obtained show that for over-expanded jets, the reflector substantially increases the jet spreading rate and reduces the supersonic length of the jet, compared with moderately under-expanded jets. The reflector’s effect appears more significant in imperfectly expanded jets that have strong shock cell structures, but is negligible in correctly expanded jet.     

一种三阶有限体积法及其在欠膨胀射流激波结构数值模拟中的应用

以喷管出口欠膨胀射流为研究对象,在Lagrange坐标系下建立欠膨胀射流二维积分形式的流动方程。通过在单元交接面处进行三阶ENO(essentially nonoscillatory)格式插值,构造得到一种适用于求解该方程的三阶ENO有限体积法。采用该格式对一维Sod激波管算例和喷管出口欠膨胀射流进行数值计算。计算结果表明,该方法具有高精度、基本无振荡的特点,能很好地捕捉包含激波、滑移线以及三波交点等复杂流场波系结构。计算得到的波系结构中马赫盘的位置与实验结果吻合很好,相对误差小于1.1%。     

Near field shock structure of dual co-axial jets

Experimental results on the shock structure of dual co-axial jets are presented. The effects of the geometric parameters of the inner nozzle, jet static pressure ratio (ratio of the exit plane static pressures of the inner and outer nozzles) and the ratio of outer to inner nozzle throat area on the shock structure were studied. A superimposed outer and inner jet structure was observed in the schlieren photographs. The inner flow is compressed by the outer flow resulting in the formation of a Mach disc and an exit shock. A parameter incorporating the effect of Mach number of the inner nozzle and jet static pressure ratio was found to correlate the observations regarding the Mach disc location.     

Relationship between vortex and sound in under‐expanded supersonic impinging jets

The study of an under‐expanded supersonic jet impinging on a flat plate by using large‐eddy simulation is reported. A third‐order upwind compact difference and a fourth‐order symmetric compact scheme are employed to discretize the nondimensional axisymmetric compressible Favre‐filtered Navier–Stokes equations in space, whereas the third‐order Runge–Kutta method with the total variation diminishing property is adopted to deal with the temporal discretization. The numerical simulation successfully captures the shock wave and vortex structures with different scales in the flow field. Waves with high and low frequencies traveling forward and reflecting back, and sound sources in different locations can be observed. By comparison with the frequency of the impinging tone from the experiment, it can be deduced that the change of pressure and swirling strength in the shear layer, pressure change on the impinging plate, and vortex merging in the jet shear layer are interdependent with the impinging tone. The effects of nozzle lip thickness on the impinging jet flow field have been investigated. The results show that the values of pressure fluctuation and vortex swirling strength in the shear layer near the nozzle have an extremum with the variation of the nozzle lip thickness. The results provide a theoretical foundation for the design of supersonic nozzles. Copyright © 2011 John Wiley & Sons, Ltd.     

Flow field characterization of coaxial conical and serrated (chevron) nozzles

PIV measurements were performed to provide insight into the effect of serrated (chevron) nozzles on the flow field of a coaxial circular jet. The serrations were tested on the primary nozzle. Mean flow results showed that the chevron effectively redistributes momentum from the high velocity center stream outward to the lower velocity secondary stream by creating lateral jets. This leads to a more rapid decay of the peak jet velocity and a consequent reduction in the length of the jet potential core. Local increases of up to 65% in the outer stream velocity were measured. The interaction of the secondary jets with the lower velocity outer stream produces increases in turbulent kinetic energy (TKE) near the center nozzle lip. These flow field effects correlate with the jet’s acoustic emissions: Reduction of low-frequency noise due to large scale mixing and potential core shortening, and increased high-frequency noise due to increased near-field turbulence.     

Active Control of Swirling Coaxial Jet Mixing with Manipulation of Large-Scale Vortical Structures

Large-scale vortical structures and associated mixing in methane/air swirling coaxial jets are actively controlled by manipulating the outer shear layer of the outer swirling coaxial jet with miniature flap actuators. In order to investigate the control mechanisms, stereoscopic particle image verocimetry (stereo-PIV) and plannar laser-induced fluorescence (PLIF) techniques are employed. It is found that intense vortex rings are produced in the outer shear layer in phase with the periodic flap motion regardless of the swirl number examined. The vortical structures in the inner shear layer, however, are strongly dependent on the swirl rate. This is because the central methane jet is accelerated by the negative axial pressure gradient, of which strength is determined by the swirl. As a result, the inner vortex formation is significantly suppressed at a higher swirl rate. On the other hand, at a relatively low swirl rate, the inner vortices are shed continuously and the methane jet is pinched off. This particular mode promotes the mixing of methane and air, so that the flammable mixture can be formed at an earlier stage of the jet flow development. In addition, the evolution of secondary streamwise vortices is prompted by the combination of the periodic vortex ring shedding and the swirl. They also contribute to the mixing enhancement in the downstream region.     

LES and RANS modelling of under-expanded jets with application to gaseous fuel direct injection for advanced propulsion systems

A density-based solver with the classical fourth-order accurate Runge-Kutta temporal discretization scheme was developed and applied to study under-expanded jets issued through millimetre-size nozzles for applications in high-pressure direct-injection (DI) gaseous-fuelled propulsion systems. Both large eddy simulation (LES) and Reynolds-averaged Navier-Stokes (RANS) turbulence modelling techniques were used to evaluate the performance of the new code. The computational results were compared both quantitatively and qualitatively against available data from the literature. After initial evaluation of the code, the computational framework was used in conjunction with RANS modelling (k-ω SST) to investigate the effect of nozzle exit geometry on the characteristics of gaseous jets issued from millimetre-size nozzles. Cylindrical nozzles with various length to diameter ratios, namely 5, 10 and 20, in addition to a diverging conical nozzle, were studied. This study is believed to be the first to provide a direct comparison between RANS and LES within the context of nozzle exit profiling for advanced high-pressure injection systems with the formation of under-expanded jets. It was found that reducing the length of the straight section of the nozzle by 50% resulted in a slightly higher level of under-expansion (∼2.6% higher pressure at the nozzle exit) and ∼1% higher mass flow rate. It was also found that a nozzle with 50% shorter length resulted in ∼6% longer jet penetration length. At a constant nozzle pressure ratio (NPR), a lower nozzle length to diameter ratio resulted in a noticeably higher jet penetration. It was found that with a diverging conical nozzle, a fairly higher penetration length could be achieved if an under-expanded jet formed downstream of the nozzle exit compared to a jet issued from a straight nozzle with the same NPR. This was attributed to the radial restriction of the flow and consequently formation of a relatively smaller reflected shock angle. With the conical nozzle used in this study and a 30 bar injection pressure, an under-expanded hydrogen jet exhibited ∼60% higher penetration length compared to an under-expanded nitrogen jet at 100 μs after start of injection. Moreover, the former jet exhibited ∼22% higher penetration compared to a nitrogen jet issued through the conical profile with 150 bar injection pressure.     

Numerical study of hysteresis in annular swirling jets with a stepped‐conical nozzle

This study investigates the experimentally observed hysteresis in the mean flow field of an annular swirling jet with a stepped‐conical nozzle. The flow is simulated using the Reynolds‐averaged Navier–Stokes (RANS) approach for incompressible flow with a k–ε and a Reynolds stress transport (RSTM) turbulence model. Four different flow structures are observed depending on the swirl number: ‘closed jet flow’, ‘open jet flow low swirl’, ‘open jet flow high swirl’ and ‘coanda jet flow’. These flow patterns change with varying swirl number and hysteresis at low and intermediate swirl numbers is revealed when increasing and subsequently decreasing the swirl. The influence of the inlet velocity profile on the transitional swirl numbers is investigated. When comparing computational fluid dynamics with experiments, the results show that both turbulence models predict the four different flow structures and the associated hysteresis and multiple solutions at low and intermediate swirl numbers. Therefore, a good agreement exists between experiments and numerics. Copyright © 2006 John Wiley & Sons, Ltd.     

Experiments on dynamic behavior of near-field region in variable property jet with swirling flow

The dynamic behavior of the near-field region in a coaxial variable property jet has been experimentally investigated under a swirling flow produced by rotating cylindrical inner and outer tubes, focusing on how the swirl of the outer jet affects the formation of a stagnation point in the swirling inner jet. The inner and outer jets rotate in the same direction. Air, CO₂, or He is issued from the inner tube as a variable property jet, and air is issued from the outer tube in this work. In the case of a CO₂ jet (a high-density, low-viscosity gas jet), a stagnation point flow is more easily formed than in the case of an air jet, and the stagnation point location is significantly lower than in that of the air jet. When the swirl of the outer jet is introduced, a stagnation point flow is more easily formed than in the case of a nonswirling outer jet, and the stagnation point location is much lower than in the case of a nonswirling outer jet. In the case of a He jet (a low-density and high-viscosity gas jet), the inner jet does not have a stagnation point flow, and its overall behavior remains nearly unchanged even under high swirl numbers of the inner and outer jets. These results clearly show that the density and viscosity differences between the inner and outer jets have a significant impact on the dynamic behavior of the near-field region in the coaxial swirling jet. The significant lowering of the stagnation point location can be physically explained by considering the theoretical equation obtained in this work.     

Nonlinear instability of an annular liquid sheet exposed to gas flow

Nonlinear instability and breakup of an annular liquid sheet has been modeled in this paper. The liquid sheet is considered to move axially and is exposed to co-flowing inner and outer gas streams. Also, the effect of outer gas swirl on sheet breakup has been studied. In the developed model a perturbation expansion method has been used with the initial magnitude of the disturbance as the perturbation parameter. This is a comprehensive model in that other geometries of planar sheet and a coaxial jet can be obtained as limiting cases of very large inner radius and inner radius equal to zero, respectively. In this temporal analysis, the effect of liquid Weber number, initial disturbance amplitude, inner gas-to-liquid velocity ratio, outer gas-to-liquid velocity ratio and outer gas swirl strength on the breakup time is investigated. The model is validated by comparison with earlier analytical studies for the limiting case of a planar sheet as well as with experimental data of sheet breakup length available in literature. It is shown that the linear theory cannot predict breakup of an annular sheet and the developed nonlinear model is necessary to accurately determine the breakup length. In the limiting case of a coaxial jet, results show that gas swirl destabilizes the jet, makes helical modes dominant compared to the axisymmetric mode and decreases jet breakup length. These results contradict earlier linear analyses and agree with experimental observations. For an annular sheet, it is found that gas flow hastens the sheet breakup process and shorter breakup lengths are obtained by increasing the inner and the outer gas velocity. Axially moving inner gas stream is more effective in disintegrating the annular sheet compared to axially moving outer gas stream. When both gas streams are moving axially, the liquid sheet breakup is quicker compared to that with any one gas stream. In the absence of outer gas swirl, the axisymmetric mode is the dominant instability mode. However, when outer gas flow has a swirl component higher helical modes become dominant. With increasing outer gas swirl strength, the maximum disturbance growth rate increases and the most unstable circumferential wave number increases resulting in a highly asymmetric sheet breakup with shorter breakup lengths and thinner ligaments.     

激波诱导环形SF6气柱演化的机理

基于可压缩多组分Navier-Stokes控制方程,结合5阶加权本质无振荡格式以及网格自适应加密技术和level-set方法,数值模拟了平面激波(Ma=1.23)与环形SF₆气柱(内外半径分别为8和17.5 mm)界面的相互作用过程。相比于之前的实验结果,数值模拟结果揭示了入射激波在界面内4次透射过程中的复杂波系结构,观察到透射激波在内部界面传播时形成自由前导折射结构并向自由前导冯诺依曼折射结构转换的波系演变过程;另外,界面内的复杂激波结构诱导内部下游界面上的涡量发生了3次反向;在界面演化后期,内部界面形成的“射流”结构与下游界面相互作用,诱导界面形成一对主涡、一对次级涡以及一个反向“射流”结构。定量分析了环形界面长度、宽度、位移、环量以及混合率的变化情况,结果表明,内部气柱的存在减弱了前期小涡结构合并形成大涡结构过程中对界面高度与长度的影响,同时提高了重质气体与环境气体的混合率。     

Prediction of momentum and scalar transport in turbulent swirling flows with an objective Reynolds-stress transport closure

The accurate prediction of turbulent swirling flows requires the use of a differential Reynolds-stress transport model to close the time-averaged Navier–Stokes equations. The performance of such model is largely determined by the way in which the fluctuating pressure–strain correlations are approximated. A number of alternative approximations are available, all of which depend explicitly on the mean vorticity tensor. Such dependence renders a constitutive relation inconsistent with the principle of Material Frame Indifference (MFI). In this paper, an objective model (i.e. one which is consistent with MFI) for the pressure–strain correlations is presented. This model, which was developed using Tensor Representation Theory, has fewer terms than the conventional alternatives and is therefore easier to implement in computational codes. Moreover, the model was calibrated to correctly reproduce the relative stress levels in both free and wall-bounded flows without the need to employ wall-damping corrections. The performance of this model is assessed using experimental data from both weakly- and strongly-swirled jets. Comparisons are also made with results obtained using three widely-used alternative models for the pressure–strain correlations. It is found that the objective model, although simpler in formulation than the others, yields results that are generally in closer correspondence with the data. The paper also reports on the prediction of mass transfer in a swirling jet. The case considered was that of a co-axial, strongly-swirled flow with an outer annular air stream and an inner helium jet. Swirl was imparted to the outer stream only. The concentration of helium was predicted using a differential scalar-flux transport closure. Close agreement was obtained with the measured concentrations. Analysis of the predicted mass fluxes revealed that the turbulent diffusivity is strongly anisotropic in this flow.     

可压缩流向涡与激波轴对称干扰的数值模拟

用ＮＳ方程数值模拟了可压缩流向涡和激波轴对称相互作用现象。数值模拟包括定常和非定常两种情况,计算结果分别与相应的实验进行了比较,结果表明数值模拟成功地捕捉到了激波和旋涡相互作用过程中发生的激波波面变形,激波振荡,涡核变大以及激波波后出现驻点、回流区等流场特征。提出了判断流向涡与运动激波相互作用中旋涡破碎的准则。     

Multiple pulsed hypersonic liquid diesel fuel jetsdriven by projectile impact

Further studies on high-speed liquid diesel fuel jets injected into ambient air conditions have been carried out. Projectile impact has been used as the driving mechanism. A vertical two-stage light gas gun was used as a launcher to provide the high-speed impact. This paper describes the experimental technique and visualization methods that provided a rapid series of jet images in the one shot. A high-speed video camera (10⁶ fps) and shadowgraph optical system were used to obtain visualization. Very interesting and unique phenomena have been discovered and confirmed in this study. These are that multiple high frequency jet pulses are generated within the duration of a single shot impact. The associated multiple jet shock waves have been clearly captured. This characteristic consistently occurs with the smaller conical angle, straight cone nozzles but not with those with a very wide cone angle or curved nozzle profile. An instantaneous jet tip velocity of 2680 m/s (Mach number of 7.86) was the maximum obtained with the 40 nozzle. However, this jet tip velocity can only be sustained for a few microseconds as attenuation is very rapid.Received: 13 December 2003, Accepted: 11 April 2004, Published online: 11 February 2005[/PUBLISHED]K. Pianthong: Correspondence to:      

An experimental study on the vortical structures and behaviour of jets issuing from inclined coaxial nozzles

An experimental study on inclined coaxial jets using laser-induced fluorescence and particle image velocimetry is presented here. The Reynolds numbers of the inner primary jet and outer secondary jet were Re = 2,500 and between Re = 500 and 2,000 (based on gap size), respectively, which corresponded to secondary-to-primary jet velocity ratios (VR) of VR = 0.5–2.0. The secondary-to-primary jet area ratio was 2.25, and 45° and 60° incline-angles were studied. Flow visualizations show that relatively independent inclined primary and secondary jet vortex roll-ups were formed at VR = 0.5. At VR = 1.0, regular pairings and mergings between primary and secondary jet vortex roll-ups led to large-scale entrainment of secondary jet and ambient fluids into the primary jet column and conferred a “serpentile”-shaped outline upon it. While the “serpentile”-shaped outline continued to exist at VR = 2.0, it was a result of stronger secondary jet inner vortex roll-ups which “pinched” the primary jet column regularly. These flow behaviours are observed to intensify with an increase in the incline-angle used. Velocity measurements demonstrate that inclined coaxial nozzles promoted vectoring of the primary jet momentum towards the longer nozzle lengths when velocity-ratio and/or incline-angle were increased. Lastly, peak velocity and higher turbulence intensity levels due to augmented vortical interactions are also detected along shorter nozzle lengths.     

LES study of global instability in annular jets

The paper presents results of numerical investigations of annular non–swirling jets. The research is performed with the use of large eddy simulations (LES) and a high–order code based on the compact differences and the Fourier approximation combined with a projection method for pressure–velocity coupling. We consider various Reynolds numbers and inner shear layer thicknesses and observe that both parameters affect the dynamics of the flow in such a way that the formation of the turbulent structures follows two different scenarios. For the particular group of the inlet conditions the flow pattern is found to be similar to low–swirling jets. The instantaneous flow–fields reveal the formation of the spiral structures located in the inner and outer mixing layers. The analysis shows that the structures are the result of the instability that leads to the precession of the recirculation region formed in the near–field. The winding sense of the spirals are clockwise or counterclockwise depending on the case and azimuthal distribution of perturbations.     

Calculation of hysteresis and flow-rate oscillations of base pressure in supersonic annular nozzles

The interaction of the turbulent axisymmetric near wake behind the face of the central body of an annular nozzle with the supersonic annular jet discharging from this nozzle is analyzed. The flow in the monoparametric near wake is calculated by the integral method [1] while the flow in the nonviscous jet is calculated by the method of through calculation using a monotonic explicit difference system of the first order of accuracy [2]. The interaction between the nonviscous and turbulent streams is determined by the displacement thickness of the wake. The initial conditions of the wake are determined from the integral conditions of attachment with the mixing flow in the isobaric base region. The interaction flow is described by the particular solution of the equations which passes through the singular saddle point — the throat of the wake. The near wake and base pressure in different modes of discharge from an annular nozzle at the exit cross section of which the ratio of outer and inner radii is y₂/y₁ = 1.3 and the Mach number is M = 2.54 are calculated as an example. The region of hysteresis of the base pressure, connected with the ambiguity of the interaction flow owing to the formation of the throat of the wake within the first or second barrel of the jet, and the parameters of the low-frequency flow-rate oscillations of base pressure in this region are determined. The results of the calculations are in satisfactory agreement with experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 125–130, January–February, 1977.     

Manipulation of Large-Scale Vortical Structures and Mixing in a Coaxial Jet with Miniature Jet Actuators Toward Active Combustion Control

Manipulation of large-scale vortical structures and associated mixing in a methane-air coaxial jet is carried out by using miniature jet actuators installed on the inner surface of the annular nozzle. The periodic radial miniature jet injections are achieved with a rapid-response servo-valve. The spatio-temporal primary jet structures are investigated through phase-locked 2C-PIV (2 Component Particle Image Velocimetry) and stereoscopic-PIV. It is found that intense ring-like vortices are produced perfectly in phase with the periodic miniature jet injections regardless of the valve-driven frequency f_v examined. When the Strouhal number St_v, which is defined with f_v, is larger than unity, the ring-like vortices are densely formed and thus methane/air mixing is prompted with low periodic fluctuation. The diameter of the vortices becomes small as St_v is increased, so that the transport range of the inner methane and outer air fluids can be controlled by changing St_v. In addition, the evolution of counter-rotating vortex pair is also observed in the cross-sectional plane. These streamwise vortices are directly formed as a result of the radial miniature jet injection, which leads to entrainment of the ambient fluid near the primary jet shear layer, and they also contribute to the mixing enhancement. Moreover, it is demonstrated that coaxial jet flame characteristics such as carbon monoxide (CO) emission and flame holding can be drastically improved under different equivalence ratios by the present jet control scheme.     

环形通道内湍流旋流流动的数值模拟

本文应用一种考虑湍流－旋流相互作用及湍流脉动各向异性的新的代数Ｒｅｙｎｏｌｄｓ应力模型,对环形通道内的湍流旋流流动进行了数值模拟,研究了改主为旋流流数,进口轴向速度及半径比等参数对环形通道内湍流流动的影响,以及对强化环形通道内传热的作用。     

Coherent structures in unsteady swirling jet flow

An LDA technique and phase-averaging analysis were used to study unsteady precessing flow in a model vortex burner. Detailed measurements were made for Re=15,000 and S=1.01. On the basis of the analysis of phase-averaged data and vortex detection by the λ₂-technique of Joeng and Hussain (1995), three precessing spiral vortex structures were identified: primary vortex (PV), inner secondary vortex (ISV), and outer secondary vortex (OSV). The PV is the primary and most powerful structure as it includes primary vorticity generated by the swirler; the ISV and OSV are considered here as secondary vortical structures. The jet breakdown zone is the conjunction of a pair of co-rotating co-winding spiral vortices, PV and ISV. The interesting new feature described is that the secondary vortices form a three-dimensional vortex dipole with a helical geometry. The effect of coupling of secondary vortices was suggested as a mechanism of enhanced stability reflected in their increased axial extent.