The unsteady development of a turbulent wake through a downstream low-pressure turbine blade passage

This paper presents two-dimensional LDA measurements of the convection of a wake through a low-pressure (LP) turbine cascade. Previous studies have shown the wake convection to be kinematic but have not provided details of the turbulent field. The spatial resolution of these measurements has facilitated the calculation of the production of turbulent kinetic energy and this has revealed a mechanism for turbulence production as the wake converts through the bladerow. The measured ensemble-averaged velocity field confirmed the previously reported kinematics of wake convection while the measurements of the turbulence quantities showed the wake fluid to be characterised by elevated levels of turbulent kinetic energy (TKE) and to have an anisotropic structure. Based on the measured mean and turbulence quantities, the production of turbulent kinetic energy was calculated. This highlighted a TKE production mechanism that resulted in increased levels of turbulence over the rear suction surface where boundary layer transition occurs. The turbulence production mechanism within the bladerow was also observed to produce more nearly isotropic turbulence. Production occurs when the principal stresses within the wake are aligned with the mean strains. This coincides with the maximum distortion of the wake within the blade passage and provides a mechanism for the production of turbulence outside of the boundary layer.

Reynolds number influence on turbulent blocking effects of a rigid plane boundary (2nd report, direct numerical simulation using townsend's eddy forcing)

The Reynolds number influence on turbulent blocking effects by a rigid plane boundary is studied using direct numerical simulation (DNS). A new forcing method using 'simple model eddies' (Townsend 1976) for DNS of stationary homogeneous isotropic turbulence is proposed. A force field is obtained in real space by sprinkling many space-filling 'simple model eddies' whose centers are randomly but uniformly distributed in space and whose axes of rotation are random. The method is applied to a shear-free turbulent boundary layer over a rigid plane boundary and the blocking effects are investigated. The results show that stationary homogeneous isotropic turbulence is generated in real space using the present method. By using different model eddies with different sizes and rotation speeds, we could change the turbulence properties such as the integral and micro scales, the turbulent Reynolds number and the isotropy of turbulence. Turbulence intensities near the wall showed good agreements with the previous measurement and the linear analysis based on a rapid distortion theory (RDT). The splat effect (i.e., turbulence intensities of the components parallel to the boundary are amplified) occurs near the boundary and the viscous effect prohibits the splat effect at the quasi steady state at low Reynolds number.

Secondary flow near an undulatory surface induced by wall blocking effect

Prandtl's secondary mean motions of the second kind near an undulating surface were explained in terms of turbulent blocking effect and kinematic boundary conditions at the surface, and its order of magnitude was estimated. Isotropic turbulence is distorted by the undulating surface of wavelength λ and amplitude h with a low slope, so that h « λ. The prime mechanism for generating the mean flow is that the far-field Isotropic turbulence is distorted by the non-local blocking effect of the surface to become anisotropic axisymmetric turbulence near the surface with principal axis that is not aligned with the local curvature of the undulation. Then the local analysis can be applied and the mechanism is similar to the mean flow generation mechanism for homogeneous axisymmetric turbulence over a planer surface, i.e. gradients of the Reynolds stress caused by the turbulent blocking effect generate the mean motions. The results from this simple analysis are consistent with previous exact analysis in which the effects of curvature are strictly taken into account. The results also qualitatively agree with flow visualization over an undulating surface in a mixing-box.

Contrasting code performances for computational aeroacoustics of jets

The numerical propagation of subcritical Tollmein-Schlichting (T-S), inviscid vortical and cut-on acoustic waves is explored. For the former case, the performances of the very different NEAT, NTS, HYDRA, FLUXp and OSMIS3D codes is studied. A modest/coarse hexahedral computational grid that starkly shows differences between the different codes and schemes used in them is employed. For the same order of discretization the five codes show similar results. The unstructured codes are found to propagate vortical and acoustic waves well on triangular cell meshes but not the T-S wave. The above code contrasting exercise is then carried out using implicit LES or Smagorinsky LES for and Ma = 0.9 plane jet on modest 0.5 million cell grids moving to circa 5 million cell grids. For this case, even on the coarse grid, for all codes results were generally encouraging. In general, the spread in computational results is less than the spread of the measurements. Interestingly, the finer grid turbulence intensity levels are slightly more under-predicted than those of the coarse grid. This difference is attributed to the numerical dispersion error having a favourable coarse grid influence. For a non-isothermal jet, HYDRA and NTS also give encouraging results. Peak turbulence values along the jet centreline are in better agreement with measurements than for the isothermal jets. Copyright © 2006 by University of Wales.

Random lasing in a bistable smectic A liquid crystal

In this work, we examine the phenomenon of random lasing from the smectic A liquid crystal phase. We summarise our results to date on random lasing from the smectic A phase including the ability to control the output from the sample using applied electric fields. In addition, diffuse random lasing is demonstrated from the electrohydrodynamic instabilities of a smectic A liquid crystal phase that has been doped with a low concentration of ionic impurities. Using a siloxane-based liquid crystal doped with ionic impurities and a laser dye, nonresonant random laser emission is observed from the highly scattering texture of the smectic A phase which is stable in zero-field. With the application of a low frequency alternating current electric field, turbulence is induced due to motion of the ions. This is accompanied by a decrease in the emission linewidth and an increase in the intensity of the laser emission. The benefit in this case is that a field is not required to maintain the texture as the scattering and homeotropic states are both stable in zero field. This offers a lower power consumption alternative to the electric-field induced static scattering sample.

Investigation of extinction in unsteady flames in turbulent combustion by 2D-LIF of OH radials and flamelet analysis

An experimental and theoretical investigation of premixed turbulent combustion in an engine simulator is presented. The distribution of hydroxyl radicals formed in the combustion of propane/air mixtures was visualized by 2D-LIF and used to monitor the progress of the combustion process. For stoichiometric mixtures, images showed a continuous wrinkled flame front, while in lean (λ=1.5) mixtures, local flame extinction was observed as discontinuities in the reaction zone. A bright active reaction zone was still observed in flame inlets and closed concave structures. The effects of self-absorption and of collisional quenching on the fluorescence signal are considered and appear to have only a minor net influence on the shape and width of the flame front. The images are evaluated and interpreted in terms of the Lewis number effect and the laminar flamelet model. Analysis was performed by determining the contour lines of the images (specifically, the ratios of average maximum to equilibrium OH concentration) and comparing with corresponding ratios from unstrained flame simulations. The results show that although the degree of turbulence is not high enough for straining effects to be important, flamelet curvature does play a significant role in the combustion of lean mixtures; this is manifested by a mean effective flame velocity that is less than the laminar burning velocity. © 1991 Combustion Institute.

Compressor wake/leading-edge interactions at off design incidences

In this paper, the effects of wake/leading-edge interactions were studied at off-design conditions. Measurements were performed on the stator-blade suction surface at midspan. The leading-edge flow-field was investigated using hotwire micro-traverses, hotfilm surface shear-stress sensors and pressure micro-tappings. The trailing-edge flow-field was investigated using hotwire boundary-layer traverses. Unsteady CFD calculations were also performed to aid the interpretation of the results. At low flow coefficients, the time-averaged momentum thickness of the leading-edge boundary layer was found to rise as the flow coefficient was reduced. The time-resolved momentum-thickness rose due to the interaction of the incoming rotor wake. As the flow coefficient was reduced, the incoming wakes increased in pitch-wise extent, velocity deficit and turbulence intensity. This increased both the time-resolved rise in the momentum thickness and the turbulent spot production within the wake affected boundary-layer. Close to stall, a drop in the leading-edge momentum thickness was observed in-between wake events. This was associated with the formation of a leading-edge separation bubble in-between wake events. The wake interaction with the bubble gave rise to a shedding phenomenon, which produced large length scale disturbances in the surface shear stress. Copyright © 2008 by ASME.

The prediction of hot streak migration in a high-pressure turbine

Accurate predictions of combustor hot streak migration enable the turbine designer to identify high-temperature regions that can limit component life. It is therefore important that these predictions are achieved within the short time scales of a design process. This article compares temperature measurements of a circular hot streak through a turning duct and a research turbine with predictions using a three-dimensional Reynolds-averaged Navier-Stokes solver. It was found that the mixing length turbulence model did not predict the hot streak dissipation accurately. However, implementation of a very simple model of the free stream turbulence (FST) significantly improved the exit temperature predictions on both the duct and research turbine. One advantage of the simple FST model described over more complex alternatives is that no additional equations are solved. This makes the method attractive for design purposes, as it is not associated with any increase in computational time.

Velocity distributions for low pressure turbines

A parametric set of velocity distributions has been investigated using a flat plate experiment. Three different diffusion factors and peak velocity locations were tested. These were designed to mimic the suction surfaces of Low Pressure (LP) turbine blades. Unsteady wakes, inherent in real turbomachinery flows, were generated using a moving bar mechanism. A turbulence grid generated a freestream turbulence level that is believed to be typical of LP turbines. Measurements were taken across a Reynolds number range of 50,000-220,000 at three reduced frequencies (0.314, 0.628, 0.942). Boundary layer traverses were performed at the nominal trailing edge using a Laser Doppler Anemometry system and hot-films were used to examine the boundary layer behaviour along the surface. For every velocity distribution tested, the boundary layer separated in the diffusing flow downstream of the peak velocity. The loss production is dominated by the mixing in the reattachment process, mixing in the turbulent boundary layer downstream of reattachment and the effects of the unsteady interaction between the wakes and the boundary layer. A sensitive balance governs the optimal location of peak velocity on the surface. Moving the velocity peak forwards on the blade was found to be increasingly beneficial when bubblegenerated losses are high, i.e. at low Reynolds number, at low reduced frequency and at high levels of diffusion. Copyright © 2008 by ASME.

Detached-eddy simulation of transonic flow past a fan-blade section

Detached-eddy simulation of transonic flow past a thin section of a fan blade has been carried out. The inflow Mach number is 1.03, and a bow shock forms upstream of the blade. The shock (corresponding to an adjacent blade) impinges on the suction-side boundary layer which causes separation and rapid transition to turbulence. The boundary layer later re-attaches near the trailing edge. The pressure-side boundary layer transitions near the leading edge and remains attached. Mean surface pressure shows basic agreement with a steady RANS calculation; strong shock motion in the DES is the major cause of discrepancy. Surface pressure spectra are investigated, and low-frequency two-dimensional disturbances associated with the shock motion are dominant. Removing the two-dimensional component from the spectra, the pressure-side three-dimensional spectra reproduce the spectral shape given by a correlation for flat-plate boundary layer wall-pressure spectra developed by Goody. 1 The suction-side disturbances produce similar high- and intermediate-frequency scalings despite substantially different boundary layer development. Near-wake results show that disturbance kinetic energy peaks at the suction-side inflection point of the mean profile, and that the energy is concentrated at low frequencies relative to the near-trailing edge surface pressure. Copyright © 2009 by the authors.

Bypass transition to sustained thermoacoustic oscillations in a linearly stable rijke tube

In this paper we examine triggering in a simple linearly-stable thermoacoustic system using techniques from flow instability and optimal control. Firstly, for a noiseless system, we find the initial states that have highest energy growth over given times and from given energies. Secondly, by varying the initial energy, we find the lowest energy that just triggers to a stable periodic solution. We show that the corresponding initial state grows first towards an unstable periodic solution and, from there, to the stable periodic solution. This exploits linear transient growth, which arises due to nonnormality in the governing equations and is directly analogous to bypass transition to turbulence. Thirdly, we introduce noise that has similar spectral characteristics to this initial state. We show that, when triggering from low noise levels, the system grows to high amplitude self-sustained oscillations by first growing towards the unstable periodic solution of the noiseless system. This helps to explain the experimental observation that linearly-stable systems can trigger to self-sustained oscillations even with low background noise. © 2010 by University of Cambridge. Published by the American Institute of Aeronautics and Astronautics, Inc.

Assessment of CFD models for shock boundary-layer interaction

A workshop on the computational fluid dynamics (CFD) prediction of shock boundary-layer interactions (SBLIs) was held at the 48th AIAA Aerospace Sciences Meeting. As part of the workshop numerous CFD analysts submitted solutions to four experimentally measured SBLIs. This paper describes the assessment of the CFD predictions. The assessment includes an uncertainty analysis of the experimental data, the definition of an error metric and the application of that metric to the CFD solutions. The CFD solutions provided very similar levels of error and in general it was difficult to discern clear trends in the data. For the Reynolds Averaged Navier-Stokes methods the choice of turbulence model appeared to be the largest factor in solution accuracy. Large-eddy simulation methods produced error levels similar to RANS methods but provided superior predictions of normal stresses.