On the effect of heat release in turbulence spectra of non-premixed reacting shear layers

2009 ◽  
Vol 626 ◽  
pp. 67-109 ◽  
Author(s):  
R. KNAUS ◽  
C. PANTANO

Velocity, mixture fraction and temperature spectra obtained from five direct numerical simulations of non-reacting and reacting shear layers, using the infinitely fast chemistry approximation, are analysed. Two different global chemical reactions corresponding to methane and hydrogen combustion with air, respectively, are considered. The effect of heat release, i.e. density variation, on the inertial and dissipation turbulence subrange of the spectra is investigated. Analysis of the database supports the experimentally available measurements of spectra in turbulent reacting flows showing that heat release effects can be scaled out by utilizing Favre-averaged (density-weighted) large-scale turbulence quantities. This is supported by the simulation results for velocity and mixture fraction in our moderate-Reynolds-number flows but it appears to be less supported in the dissipation subrange of the temperature spectra. The departure from universal scaling using Favre-averaged quantities in the temperature spectrum, which is evident in the dissipation subrange, appears to be caused by the strong nonlinearity of the state relationship relating the mixture fraction to the temperature, as has been suggested previously. These effects are less pronounced at intermediate wavenumbers. Analysis suggests that the nonlinear state relationship and the spectra of mixture fraction moments can be used to reconstruct the temperature spectrum across the flow. Moreover, the governing equation for the temperature variance is analysed to identify a possible surrogate for the overall rate of dissipation of temperature fluctuations and their corresponding dissipation length scale. This scaling analysis is then used to separate planes across the shear layer where the temperature dissipation length scale is alike that of the mixture fraction from regions where smaller length scales are present, and are evidenced in the dissipation subrange using Kolmogorov scaling. In our simulations, these regions correspond to the centre of the shear layer and the mean flame location. The new estimate for the temperature dissipation length scale is able to collapse the compensated spectra profiles at all planes across the shear layer for all simulations.

1978 ◽  
Vol 88 (3) ◽  
pp. 541-562 ◽  
Author(s):  
R. J. Hill

Several models are developed for the high-wavenumber portion of the spectral transfer function of scalar quantities advected by high-Reynolds-number, locally isotropic turbulent flow. These models are applicable for arbitrary Prandtl or Schmidt number, v/D, and the resultant scalar spectra are compared with several experiments having different v/D. The ‘bump’ in the temperature spectrum of air observed over land is shown to be due to a tendency toward a viscous-convective range and the presence of this bump is consistent with experiments for large v/D. The wavenumbers defining the transition between the inertial-convective range and viscous-convective range for asymptotically large v/D (denoted k* and k1* for the three- and one-dimensional spectra) are determined by comparison of the models with experiments. A measurement of the transitional wavenumber k1* [denoted (k1*)s] is found to depend on v/D and on any filter cut-off. On the basis of the k* values it is shown that measurements of β1 from temperature spectra in moderate Reynolds number turbulence in air (v/D = 0·72) maybe over-estimates and that the inertial-diffusive range of temperature fluctuations in mercury (v/D ≃ 0·02) is of very limited extent.


2018 ◽  
Vol 849 ◽  
pp. 1096-1119 ◽  
Author(s):  
D. C. Lander ◽  
D. M. Moore ◽  
C. W. Letchford ◽  
M. Amitay

Scaling characteristics, essential to the mechanisms of transition in square-prism shear layers, were explored experimentally. In particular, the evolution of the dominant instability modes as a function of Reynolds number were reported in the range $1.5\times 10^{4}\lesssim Re_{D}\lesssim 7.5\times 10^{4}$. It was found that the ratio between the shear layer frequency and the shedding frequency obeys a power-law scaling relation. Adherence to the power-law relationship, which was derived from hot-wire measurements, has been supported by two additional and independent scaling considerations, namely, by particle image velocimetry measurements to observe the evolution of length and velocity scales in the shear layer during transition, and by comparison to direct numerical simulations to illuminate the properties of the front-face boundary layer. The nonlinear dependence of the shear layer instability frequency is sustained by the influence of $Re_{D}$ on the thickness of the laminar front-face boundary layer. In corroboration with the original scaling argument for the circular cylinder, the length scale of the shear layer was the only source of nonlinearity in the frequency ratio scaling, within the range of Reynolds numbers reported. The frequency ratio scaling may therefore be understood by the influence of $Re_{D}$ on the appropriate length scale of the shear layer. This length scale was observed to be the momentum thickness evaluated at a transition point, defined where the Kelvin–Helmholtz instability saturates.


Author(s):  
Naman Jain ◽  
Hieu Pham ◽  
Xinyi Huang ◽  
Sutanu Sarkar ◽  
Xiang Yang ◽  
...  

Abstract Buoyant shear layers encountered in many engineering and environmental applications have been studied by researchers for decades. Often, these flows have high Reynolds and Richardson numbers, which leads to significant/intractable space-time resolution requirements for DNS or LES. On the other hand, many of the important physical mechanisms, such as stress anisotropy, wake stabilization, and regime transition, inherently render eddy viscosity-based RANS modeling inappropriate. Accordingly, we pursue second-moment closure (SMC), i.e., full Reynolds stress/flux/variance modeling, for moderate Reynolds number non-stratified, and stratified shear layers for which DNS is possible. A range of sub-model complexity is pursued for the diffusion of stresses, density fluxes and variance, pressure strain and scrambling, and dissipation. These sub-models are evaluated in terms of how well they are represented by DNS in comparison to the exact Reynolds averaged terms, and how well they impact the accuracy of full RANS closure. For the non-stratified case, SMC model predicts the shear layer growth rate and Reynolds shear stress profiles accurately. Stress anisotropy and budgets are captured only qualitatively. Comparing DNS of exact and modeled terms, inconsistencies in model performance and assumptions are observed, including inaccurate prediction of individual statistics, non-negligible pressure diffusion, and dissipation anisotropy. For the stratified case, shear layer and gradient Richardson number growth rates, and stress, flux and variance decay rates, are captured with less accuracy than corresponding flow parameters in the non-stratified case. These studies lead to several recommendations for model improvement.


2021 ◽  
Author(s):  
Naman Jain ◽  
Xinyi Huang ◽  
Xiang Yang ◽  
Robert Kunz ◽  
Hieu T. Pham ◽  
...  

Abstract Buoyant shear layers are encountered in many engineering and environmental applications, and have been studied by researchers in the context of experiments and modeling for decades. Often, these flows have high Reynolds and Richardson numbers, and this leads to significant/intractable space-time resolution requirements for DNS or LES modeling. On the other hand, many of the important physical mechanisms in these systems, such as stress anisotropy, wake stabilization, and regime transition, inherently render eddy viscosity-based RANS modeling inappropriate. Accordingly, we pursue second-moment closure (SMC), i.e., full Reynolds stress/flux/variance modeling, for moderate Reynolds number non-stratified, and stratified shear layers for which DNS is possible. A range of sub-model complexity is pursued for the diffusion of stresses, density fluxes and variance, pressure strain and scrambling, and dissipation. These sub-models are evaluated in terms of how well they are represented by DNS in comparison to the exact Reynolds averaged terms, and how well they impact the accuracy of the full RANS closure. For the non-stratified case, the SMC model predicts the shear layer growth rate and Reynolds shear stress profiles accurately. Stress anisotropy and budgets are captured only qualitatively. Comparing DNS of exact and modeled terms, inconsistencies in model performance and assumptions are observed, including inaccurate prediction of individual statistics, non-negligible pressure diffusion, and dissipation anisotropy. For the stratified case, shear layer and gradient Richardson number growth rates, and stress, flux and variance decay rates, are captured with less accuracy than corresponding flow parameters in the non-stratified case. These studies lead to several recommendations for model improvement.


2006 ◽  
Vol 128 (4) ◽  
pp. 874-879 ◽  
Author(s):  
Roberto C. Aguirre ◽  
Jennifer C. Nathman ◽  
Haris C. Catrakis

Flow geometry effects are examined on the turbulent mixing efficiency quantified as the mixture fraction. Two different flow geometries are compared at similar Reynolds numbers, Schmidt numbers, and growth rates, with fully developed turbulence conditions. The two geometries are the round jet and the single-stream planar shear layer. At the flow conditions examined, the jet exhibits an ensemble-averaged mixing efficiency which is approximately double the value for the shear layer. This substantial difference is explained fluid mechanically in terms of the distinct large-scale entrainment and mixing-initiation environments and is therefore directly due to flow geometry effects.


2018 ◽  
Vol 30 (9) ◽  
pp. 4141-4154
Author(s):  
Abbas Ebrahimi ◽  
Majid Hajipour ◽  
Kamran Ghamkhar

PurposeThe purpose of this paper is to control flow separation over a NACA 4415 airfoil by applying unsteady forces to the separated shear layers using dielectric barrier discharge (DBD) plasma actuators. This novel flow control method is studied under conditions which the airfoil angle of attack is 18°, and Reynolds number based on chord length is 5.5 × 105.Design/methodology/approachLarge eddy simulation of the turbulent flow is used to capture vortical structures through the airfoil wake. Power spectral density analysis of the baseline flow indicates dominant natural frequencies associated with “shear layer mode” and “wake mode.” The wake mode frequency is used simultaneously to excite separated shear layers at both the upper surface and the trailing edge of the airfoil (dual-position excitation), and it is also used singly to excite the upper surface shear layer (single-position excitation).FindingsBased on the results, actuations manipulate the shear layers instabilities and change the wake patterns considerably. It is revealed that in the single-position excitation case, the vortices shed from the upper surface shear layer are more coherent than the dual-position excitation case. The maximum value of lift coefficient and lift-to-drag ratio is achieved, respectively, by single-position excitation as well as dual-position excitation.Originality/valueThe paper contributes to the understanding and progress of DBD plasma actuators for flow control applications. Further, this research could be a beneficial solution for the promising design of advanced low speed flying vehicles.


Author(s):  
Kiran Manoharan ◽  
Travis Smith ◽  
Benjamin Emerson ◽  
Christopher M. Douglas ◽  
Tim Lieuwen ◽  
...  

This study is motivated by the necessity to develop a low order prediction approach for unsteady heat release response characteristics in lean premixed gas turbine combustors. This in turn requires an accurate description of the coherent hydrodynamic oscillations induced in the combustor flow by acoustic forcing. Time resolved velocity and flame position fields are obtained using sPIV and OH-PLIF measurements on a single nozzle, swirl-stabilized, premixed, methane-air flame in a model “unwrapped” annular combustor rig. A natural acoustic oscillation in the rig at 115 Hz results in a coherent flow oscillation that is concentrated primarily within the shear layer between the annular jet flow and the central recirculation zone. A linear stability analysis performed about time averaged base flow fields shows that the flow does not have any self-excited hydrodynamic modes. We then compare predictions from a forced response analysis at a forcing frequency of 115 Hz, based on the linearized Navier-Stokes equations for this coherent response. Good qualitative agreement between linear forced response analysis predictions and experimental response results, is seen for the spatial variation of velocity oscillation amplitude fields, away from the burner centerline. Further, good quantitative agreement between predictions and the experimental response is seen for the phase speed of velocity oscillations along the shear layer between the annular jet and the central recirculation zone. This phase velocity is an important flow field characteristic that has a significant impact on the heat release response that results from these coherent velocity oscillations. Present methods for forced response analysis assume uniform forcing amplitude along the radial direction at the forcing location, as well as, open flows along the streamwise direction. Both these assumptions are not strictly true for the present burner which has a center body on its axis. This maybe the reason for somewhat poor qualitative and quantitative agreement between experiments and predictions at the centerline.


2005 ◽  
Vol 53 (620) ◽  
pp. 408-413
Author(s):  
Mikiya Araki ◽  
Jun Osaka ◽  
Osamu Imamura ◽  
Mitsuhiro Tsue ◽  
Michikata Kono

2016 ◽  
Vol 811 ◽  
pp. 37-50 ◽  
Author(s):  
Giuseppe A. Rosi ◽  
David E. Rival

A constantly accelerating circular plate was investigated towards understanding the effect of non-stationarity on shear-layer entrainment and topology. Dye visualizations and time-resolved particle image velocimetry measurements were collected for normalized accelerations spanning three orders of magnitude. Increasing acceleration acts to organize shear-layer topology. Specifically, the Kelvin–Helmholtz instabilities within the shear layer better adhered to a circular path and exhibited consistent and repeatable spacing. Normalized starting-vortex circulation was observed to collapse with increasing acceleration, which one might not expect due to increased levels of mixing at higher instantaneous Reynolds numbers. The entrainment rate was shown to increase nonlinearly with increasing acceleration. This was attributed to closer spacing between instabilities, which better facilitates the roll-up of fluid between the shear layer and vortex core. The shear-layer organization observed at higher accelerations was associated with smaller spacings between instabilities. Specifically, analogous point-vortex simulations demonstrated that decreasing the spacing between instabilities acts to localize and dampen perturbations within an accelerating shear layer.


Author(s):  
C. P. Premchand ◽  
Manikandan Raghunathan ◽  
Midhun Raghunath ◽  
K. V. Reeja ◽  
R. I. Sujith ◽  
...  

Abstract The tonal sound production during thermoacoustic instability is detrimental to the components of gas turbine and rocket engines. Identifying the root cause and controlling this oscillatory instability would enable manufacturers to save in costs of power outages and maintenance. An optimal method is to identify the structures in the flow-field that are critical to tonal sound production and perform control measures to disrupt those “critical structures”. Passive control experiments were performed by injecting a secondary micro-jet of air onto the identified regions with critical structures in the flow-field of a bluff-body stabilized, dump, turbulent combustor. Simultaneous measurements such as unsteady pressure, velocity, local and global heat release rate fluctuations are acquired in the regime of thermoacoustic instability before and after control action. The tonal sound production in this combustor is accompanied by a periodic flapping of the shear layer present in the region between the dump plane (backward-facing step) and the leading edge of the bluff-body. We obtain the trajectory of Lagrangian saddle points that dictate the flow and flame dynamics in the shear layer during thermoacoustic instability accurately by computing Lagrangian Coherent Structures. Upon injecting a secondary micro-jet with a mass flow rate of only 4% of the primary flow, nearly 90% suppression in the amplitude of pressure fluctuations are observed. The suppression thus results in sound pressure levels comparable to those obtained during stable operation of the combustor. Using Morlet wavelet transform, we see that the coherence in the dominant frequency of pressure and heat release rate oscillations during thermoacoustic instability is affected by secondary injection. The disruption of saddle point trajectories breaks the positive feedback loop between pressure and heat release rate fluctuations resulting in the observed break of coherence. Wavelet transform of global heat release rate shows a redistribution of energy content from the dominant instability frequency (acoustic time scale) to other time scales.


Sign in / Sign up

Export Citation Format

Share Document