Effect of Streamline Curvature on Intensity of Streamwise Vortices in the Mixing Layer of Supersonic Jets

The structure of supersonic nonisobaric jets with Mach numbers M_a = 1 and 2 is considered experimentally to find the effect of streamline curvature on the evolution of streamwise vortices in the mixing layer. The spatial development of steady streamwise vortices in the mixing layer of supersonic jets is considered. A method for generation of steady streamwise vortices by applying microroughness elements of controlled size onto the inner surface of the nozzle is developed. Radial profiles and azimuthal variations of total pressure are obtained; the mixinglayer thickness and the curvature of streamlines in supersonic jets are determined. A significant effect of microroughness elements of prescribed shape located on the nozzle surface on the behavior of total pressure in the mixing layer of supersonic jets, as compared to natural disturbances, is obtained.     

超声速来流横向狭缝喷流标量输运的混合RANS/LES模拟

采用基于时均湍流模型的混合RANS/LES模拟方法对超声速来流狭缝横向喷流标量输运过程进行大涡模拟,控制方程对流项采用五阶精度的WENO格式求解,过滤后的组分方程中的亚格子组分对流通量采用梯度扩散模化.模拟得到了超声速来流狭缝横向喷流大尺度涡和组分的演化过程,研究表明喷流具有明显的非定常周期性特点,并且组分随时间的空间分布特征决定于流动的大尺度结构.模拟得到的压力的统计时均结果和实验结果定性一致,但喷流前后的分离区较实验值过大,这表明混合RANS/LES需要进一步改进.     

Decomposition of Mixing Layer Turbulence into Coherent Structures and Background Fluctuations

The eduction of coherent structures from cross-wire rake data in a fully turbulent incompressible mixing layer confirms the feasibility of a decomposition of a turbulent flow field, first suggested by Farge, as non-periodic non-equilibrium coherent structures interacting with a ‘thermalized’ broad-band turbulence. A simple wavelet coefficient decimation algorithm and orthogonalization yields non-periodic dominant flow structures and a background field that has a Gaussian distribution of velocities at the centerline. The coherent structures are classified in terms of their topology. The non-coherent background field has flat energy spectra and normal distribution of velocity components. Most background field statistics depend only weakly on the type of structure on which they are superposed. It may be possible to adapt existing subgrid scale models to this decomposition. This revised version was published online in July 2006 with corrections to the Cover Date.     

Mixing Measurements in a Supersonic Expansion-Ramp Combustor

This paper reports results on molecular mixing for injection via an expansion-ramp into a supersonic freestream with M ₁ = 1.5. This geometry produces a compressible turbulent shear layer between an upper, high-speed “air” stream and a lower, low-speed “fuel” stream, injected through an expansion-ramp at α = 30° to the high-speed freestream. Mass injection is chosen to force the shear layer to attach to the lower guide wall. This results in part of the flow being directed upstream, forming a recirculation zone. Employing the hypergolic hydrogen-fluorine chemical reaction and pairs of “flip” experiments, molecular mixing is quantified by measuring the resulting temperature rise. Initial experiments established the fast-chemistry limit for this flow in terms of a Damköhler number (Da). For Da ≥ 1.4, molecularly mixed fluid effectively reacts to completion. Parameters varied in these experiments were the measurement station location, the injection velocity of the (lower) “fuel” stream, the stoichiometry for the flip experiments, and the density ratio of the fuel and air streams. As expected, mixing increases with increasing distance from the injection surface. The mixed fluid fraction increases by 12% when changing the fuel-to-air stream density ratio from 1 to 0.2. Comparisons with measurements at subsonic (high-speed) “air” stream velocities show that the trend of decreasing mixing with increasing speed documented in free-shear layer flows is also encountered in these flows. The current geometry produces higher mixing levels than do free shear layers.     

Three-Wave Interactions of Disturbances in a Supersonic Boundary Layer

Within the framework of the weakly nonlinear stability theory, group interaction of disturbances in a supersonic boundary layer is considered. The disturbances are represented by two spatial packets of traveling instability waves (wave trains) with multiple frequencies. The possibility of energy redistribution in such wave systems in the case of three-wave resonant interactions of packet constituents is considered. The model is used to test the dynamics of unstable waves arising due to introduction of controlled high-intensity disturbances into a supersonic boundary layer. It is found that this mechanism is not the main one for the features of streamwise dynamics of such nonlinear waves being observed.     

可压缩混合层中细光束畸变效应的实验研究

利用双曝光CCD相机成像技术,对小直径光束穿越可压缩混合层流场存在的大尺度结构特性及其引起的气动光学效应进行了实验研究。实验研究结果表明:在不同对流马赫数(Mc=0.17,0.45)下,可压缩混合层中的大尺度结构形状不同,因而具有不同的折射率分布场,从而导致光束投影发生了不同程度的变形和偏移。利用光束质心的概念定量分析了不同对流马赫数下光束偏移和抖动的统计特性。统计结果表明在高对流马赫数下光束投影的偏移量减小,光束沿流场展向的平均偏移角度减小,但摆动幅度加大。     

Study on the Fuel Air Mixing Induced by a Shock Wave Propagating into a H2-Air Interface

2nd-order upwind TVD scheme was used to solve the laminar, fully Navier-Stokes equations. The numerical simulations were done on the propagation of a shock wave with Ma _S = 2 and 4 into a hydrogen and air mixture in a duct and a duct with a rearward step. The results indicate that a swirling vortex may be generated in the lopsided interface behind the moving shock. Meanwhile, the complex shock system is also formed in this shear flow region. A large swirling vortex is produced and the fuel mixing can be enhanced by a shock wave at low Mach number. But in a duct with a rearward step, the shock almost disappears in hydrogen for Ma _S = 2. The shock in hydrogen will become strong if Ma _S is large. Similar to the condition of a shock moving in a duct full of hydrogen and air, a large vortex can be formed in the shear flow region. The large swirling vortex even gets through the reflected shock and impacts on the lower wall. Then, the distribution of hydrogen behind the rearward step is divided into two regions. The transition from regular reflection to Mach reflection was observed as well in case Ma _S = 4.     

Research on Coherent Structures in a Mixing Layer of the Fene-P Polymer Solution

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅＭｉｘｉｎｇｌａｙｅｒｉｓｃｏｍｍｏｎｌｙｅｎｃｏｕｎｔｅｒｅｄｉｎｖａｒｉｏｕｓｃｈｅｍｉｃａｌｐｒｏｃｅｓｓｅｓ .Ｉｔｉｓｔｈｅｒｅｆｏｒｅｓｉｇｎｉｆｉｃａｎｔｔｏｓｔｕｄｙｔｈｅｂｅｈａｖｉｏｒｓｏｆｃｏｈｅｒｅｎｔｓｔｒｕｃｔｕｒｅｓｉｎａｍｉｘｉｎｇｌａｙｅｒｉｎｏｒｄｅｒｔｏｐｒｅｄｉｃｔａｎｄｃｏｎｔｒｏｌｔｈｅｅｖｏｌｕｔｉｏｎｏｆｔｈｅｆｌｏｗ ,ｔｈｅｎｉｍｐｒｏｖｅｔｈｅｍｉｘｉｎｇａｎｄｃｈｅｍｉｃａｌｒｅａｃｔｉｏｎｅｆｆｉｃｉｅｎｃｙ .Ｓ…     

混合层流动拟序结构的大涡模拟

采用大涡模拟方法对空间发展的二维平面混合层进行了数值模拟 ,动量方程采用分步投影法求解 ,亚格子项采用标准Smagorinsky亚格子模式模拟 ,压力泊松方程采用修正的循环消去法快速求解 ,同时求解了标志物输运方程以实现数值流场显示。模拟结果给出了混合层流动的瞬态发展过程以及流动中拟序结构的发展演变过程 ,成功地模拟了混合层发展中的各种瞬态细节过程 ,如涡的卷起、增长 ,涡与涡之间的配对、合并过程 ,以及大涡破碎为小涡的级联过程 ,为各种以混合层流动为原型流动的射流、尾流等工业流动的控制和优化提供了理论基础。     

Vortex-In-Cell and Probability Density Function Approach for a Passive Scalar Field in a Mixing Layer

A Reynolds-averaged simulation based on the vortex-in-cell (VIC) and the transport equation for the probability density function (PDF) of a scalar has been developed to predict the passive scalar field in a two-dimensional spatially growing mixing layer. The VIC computes the instantaneous velocity and vorticity fields. Then the mean-flow properties, i.e. the mean velocity, the root-mean-square (rms) longitudinal and lateral velocity fluctuations, the Reynolds shear stress, and the rms vorticity fluctuations are computed and used as input to the PDF equation. The PDF transport equation is solved using the Monte Carlo technique. The convection term uses the mean velocities from the VIC. The turbulent diffusion term is modeled using the gradient transport model, in which the eddy diffusivity, computed via the Boussinesq's postulate, uses the Reynolds shear stress and gradients of mean velocities from the VIC. The molecular mixing term is closed by the modified Curl model.

The computational results were compared with two-dimensional experimental results for passive scalar. The predicted turbulent flow characteristics, i.e. mean velocity and rms longitudinal fluctuations in the self-preserving region, show good agreement with the experimental measurements. The profiles of the mean scalar and the rms scalar fluctuations are also in reasonable agreement with the experimental measurements. Comparison between the mean scalar and the mean velocity profiles shows that the scalar mixing region extends further into the free stream than does the momentum mixing region, indicating enhanced transport of scalar over momentum. The rms scalar profiles exhibit an asymmetry relative to the concentration centerline, and indicate that the fluid on the high-speed side mixes at a faster rate than the fluid on the low-speed side. The asymmetry is due to the asymmetry in the mixing frequency cross-stream profiles. Also, the PDFs have peaks biased toward the high-speed side.     

An Experimental Study of the Near-Field Mixing Characteristics of a Swirling Jet

The present experimental investigation is devoted to the mixing characteristics of a passive scalar in the near-field region of a moderately swirling jet issuing from a fully developed axially rotating pipe flow. Instantaneous streamwise and azimuthal velocity components as well as the temperature were simultaneously accessed by means of a combined X-wire and cold-wire probe. The results indicate a modification of the turbulence structures to that effect that the swirling jet spreads, mixes and evolves faster compared to its non-swirling counterpart. The high correlation between streamwise velocity and temperature fluctuations as well as the streamwise passive scalar flux are even more enhanced due to the addition of swirl, which in turn shortens the distance and hence time needed to mix the jet with the ambient air.     

气固两相混合层二维涡配对的数值研究

采用双向耦合模型中的速度耦合模型,数值模拟了气固两相混合层中涡的配对、合并过程,文中采用拟谱方法对流场进行了直接数值模拟,用Lagrange模型跟踪固粒,结果发现,在双向耦合过程中Stokes数仍然是重要的参数,但不是唯一影响流场的参数,流场的发展还与固粒的体积浓度、固粒的相对密度以及固粒大小等因素有关。固粒抑制流场的发展,阻碍涡的配对、合并,加快了涡量的扩散;小St数的固粒仍然跟随流体运动,大St数的固粒趋向于低涡量区的趋势减弱。     

Determining the Optimal Contour of a Three-Dimensional Supersonic Jet Engine Outlet System

The problem of determining the optimal contour of a two-dimensional jet engine outlet system was solved in [1] using the direct method of the calculus of variations. In this study the method is developed for designing the optimal contour of a three-dimensional outlet system providing maximum thrust in a given direction.     

Evolution of Artificial Perturbations in the Boundary Layer on a Plate and in Its Wake at Supersonic Free-Stream Velocity

The evolution of artificial perturbations in the boundary layer on the flat section of a plate, on the backward-facing wedge behind the rarefaction wave fan, and in the wake is studied experimentally at the Mach number M=2.     

Mixing Characteristics of a Flapping Jet from a Self-Exciting Nozzle

Experimental results of the mixing characteristics of a low-frequency flapping jet from a self-exciting nozzle are presented. The simple fluidic device used to generate the flapping motion is also described. The nozzle contains no external trigger and, unlike the flip-flop nozzle of Viets, contains no external feedback path. Both conventional and conditional averaging schemes are employed to characterise the turbulent mixing characteristics of the jet using data obtained from hot-wire anemometry. Flow-visualisation is used to characterise the flapping motion. It is revealed that the dynamic flapping motion enhances the large-scale mixing of the jet while concurrently suppressing the generation of the fine-scale turbulence. The results also indicate that high turbulence intensities, initiated by the flapping motion, are sustained even in the far-field flow region. This revised version was published online in July 2006 with corrections to the Cover Date.     

Non-Equilibrium Mixing of Turbulence Scales Using a Continuous Hybrid RANS/LES Approach: Application to the Shearless Mixing Layer

A new subgrid-scale (SGS) model based on partially integrated transport method (PITM) is applied to the case of a turbulent spectral non-equilibrium flow created by the mixing of two turbulence fields of differing scales: the shearless mixing layer. The method can be viewed as a continuous hybrid RANS/LES approach. In this model the SGS length scale is no longer given by the size of the discretization step, but is dynamically estimated using an additional transport equation for the dissipation rate. The results are compared to those corresponding to the classical model of Smagorinsky and to the experimental data of Veeravalli and Warhaft. A method for creating an anisotropic analytical pseudo-random field for inflow conditions is also proposed. This approach based on subgrid-scale transport modelling combined with anisotropic inlet conditions gives better results for the prediction of the shearless mixing layer.     

Numerical Studies on the Mixing of CH4 and Kerosene Injected into a Supersonic Flow with H2 Pilot Injection

Two-fluid model and divisional computation techniques were used. The multi-species gas fully N-S equations were solved by upwind TVD scheme. Liquid phase equations were solved by NND scheme. The phases-interaction ODE equations were solved by 2nd Runge-Kutta approach. The favorable agreement is obtained between computational results and PLIF experimental results of iodized air injected into a supersonic flow. Then, the numerical studies were carried out on the mixing of CH ₄ and kerosene injected into a supersonic flow with H ₂ pilot injection. The results indicate that the penetration of kerosene approaches maximum when it is injected from the second injector. But the kerosene is less diffused compared with the gas fuels. The free droplet region appears in the flow field. The mixing mechanism of CH ₄ with H ₂ pilot injection is different from that of kerosene. In the staged duct, H ₂ can be entrained into both recirculation zones produced by the step and injectors. But CH ₄ can only be carried into the recirculation between the injectors. Therefore, initiations of H ₂ and CH ₄ can occur in those regions. The staged duct is better in enhancing mixing and initiation with H ₂ pilot flame.     

混合层中粒子扩散的实验研究

本文应用流动显示和LDA测量,对混合层中粒子的扩散问题进行了实验研究。结果表明：粒子的扩散强烈地依赖于St数,具有较小St的粒子的扩散近似地与流体示踪粒子相同,随着St的增大,扩散角度增大,也就是说粒子扩散加快,并且随着粒子惯性的增大,它对湍流脉动的响应将会减弱。     

Large Eddy Simulation (2D) of a Reacting Plane Mixing Layer Using Filtered Density Function Closure

The filtered density function (FDF) is implemented for a two-dimensional, large eddy simulation (LES) of a gas phase, spatially developing, reacting and non-reacting, constant-density, plane mixing layer in a flow regime prior to the mixing transition where the flow is mainly two-dimensional. The unresolved scalar fluctuations are taken into account by considering the probability density function (PDF) of subgrid scale (SGS) scalar quantities following the FDF approach. In the derived FDF transport equation, the effect of chemical reactions appears in a closed form. The Lagrangian Monte Carlo scheme is used to solve the FDF transport equation. The applicability and performance of the FDF for LES of a reacting plane mixing layer are assessed by comparisons with experimental measurements. In non-reacting flow, the calculated mean streamwise velocity profiles and mean mixture fraction profiles relax to self-similarity, which is in satisfactory agreement with the measurements. In reacting flow, the FDF calculation provided a satisfactory accuracy in comparison with measurements of mean reactant and product concentration. The increase in the total amount of product formation in the flip case demonstrates the asymmetric characteristics of the entrainment and mixing characteristics in the mixing layer. This revised version was published online in July 2006 with corrections to the Cover Date.