VGs for a normal SBLI with a downstream diffuser

New types of vortex generators for boundary layer control were investigated experimentally in a flow field which contains a Mach 1.4 normal Shockwave followed by a subsonic diffuser. A parametric study of device height and distance upstream of the normal shock was undertaken with two novel devices: ramped-vanes and split-ramps. Flowfield diagnostics included high-speed Schlieren, oil flow visualization, and pitot-static pressure measurements. A number of flowfield parameters including flow separation, pressure recovery, centerline incompressible boundary layer shape factor, and shock stability were analyzed and compared to the baseline. All configurations tested yielded an elimination of centerline flow separation with the presence of the vortex generators. However, the devices also tended to increase the three-dimensionality of the flow with increased side-wall interaction. When located 25δo upstream of the normal shock, the largest ramped-vane device (whose height was about 0.75 the incoming uncontrolled boundary layer thickness, δo) yielded the smallest centerline incompressible shape factor and the least streamwise oscillations of the normal shock. However, additional studies are needed to better understand the corner interaction effects, which are substantial. © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Behaviour of unsteady transonic shock / boundary layer interactions with three-dimensional effects

An experimental investigation into the response of transonic SBLIs to periodic down-stream pressure perturbations in a parallel walled duct has been conducted. Tests have been carried out with a shock strength of M ∞ = 1.5 for pressure perturbation frequencies in the range 16-90 Hz. Analysis of the steady interaction at M∞ = 1.5 has also been made. The principle measurement techniques were high speed schlieren photography and laser Doppler anemometry. The structure of the steady SBLI was found to be highly three-dimensional, with large corner flows and sidewall SBLIs. These aspects are thought to influence the upstream transmission of pressure information through the interaction by affecting the post-shock flow field, including the extent of regions of secondary supersonic flow. At low frequency, the dynamics of shock motion can be predicted using an inviscid analytical model. At increased frequencies, viscous effects become significant and the shock exhibits unexpected dynamic behaviour, due to a phase lag between the upstream transmission of pressure information in the core flow and in the viscous boundary layers. Flow control in the form of micro-vane vortex generators was found to have a small impact on shock dynamics, due to the effect it had on the post-shock flow field outside the viscous boundary layer region. The relationship between inviscid and viscous effects is developed and potential destabilising mechanisms for SBLIs in practical applications are suggested. Copyright © 2009 by Paul Bruce and Holger Babinsky.

Applications of the dynamic mode decomposition

The decomposition of experimental data into dynamic modes using a data-based algorithm is applied to Schlieren snapshots of a helium jet and to time-resolved PIV-measurements of an unforced and harmonically forced jet. The algorithm relies on the reconstruction of a low-dimensional inter-snapshot map from the available flow field data. The spectral decomposition of this map results in an eigenvalue and eigenvector representation (referred to as dynamic modes) of the underlying fluid behavior contained in the processed flow fields. This dynamic mode decomposition allows the breakdown of a fluid process into dynamically revelant and coherent structures and thus aids in the characterization and quantification of physical mechanisms in fluid flow. © 2010 Springer-Verlag.

Measurement-Integrated simulations and Kalman filter applied to a co-flowing jet

This paper deals with the experimental evaluation of a flow analysis system based on the integration between an under-resolved Navier-Stokes simulation and experimental measurements with the mechanism of feedback (referred to as Measurement-Integrated simulation), applied to the case of a planar turbulent co-flowing jet. The experiments are performed with inner-to-outer-jet velocity ratio around 2 and the Reynolds number based on the inner-jet heights about 10000. The measurement system is a high-speed PIV, which provides time-resolved data of the flow-field, on a field of view which extends to 20 jet heights downstream the jet outlet. The experimental data can thus be used both for providing the feedback data for the simulations and for validation of the MI-simulations over a wide region. The effect of reduced data-rate and spatial extent of the feedback (i.e. measurements are not available at each simulation time-step or discretization point) was investigated. At first simulations were run with full information in order to obtain an upper limit of the MI-simulations performance. The results show the potential of this methodology of reproducing first and second order statistics of the turbulent flow with good accuracy. Then, to deal with the reduced data different feedback strategies were tested. It was found that for small data-rate reduction the results are basically equivalent to the case of full-information feedback but as the feedback data-rate is reduced further the error increases and tend to be localized in regions of high turbulent activity. Moreover, it is found that the spatial distribution of the error looks qualitatively different for different feedback strategies. Feedback gain distributions calculated by optimal control theory are presented and proposed as a mean to make it possible to perform MI-simulations based on localized measurements only. So far, we have not been able to low error between measurements and simulations by using these gain distributions.

On the use of combined compact scheme for numerical solution of the two-layer oceanic shallow water equations

Many types of oceanic physical phenomena have a wide range in both space and time. In general, simplified models, such as shallow water model, are used to describe these oceanic motions. The shallow water equations are widely applied in various oceanic and atmospheric extents. By using the two-layer shallow water equations, the stratification effects can be considered too. In this research, the sixth-order combined compact method is investigated and numerically implemented as a high-order method to solve the two-layer shallow water equations. The second-order centered, fourth-order compact and sixth-order super compact finite difference methods are also used to spatial differencing of the equations. The first part of the present work is devoted to accuracy assessment of the sixth-order super compact finite difference method (SCFDM) and the sixth-order combined compact finite difference method (CCFDM) for spatial differencing of the linearized two-layer shallow water equations on the Arakawa's A-E and Randall's Z numerical grids. Two general discrete dispersion relations on different numerical grids, for inertia-gravity and Rossby waves, are derived. These general relations can be used for evaluation of the performance of any desired numerical scheme. For both inertia-gravity and Rossby waves, minimum error generally occurs on Z grid using either the sixth-order SCFDM or CCFDM methods. For the Randall's Z grid, the sixth-order CCFDM exhibits a substantial improvement , for the frequency of the barotropic and baroclinic modes of the linear inertia-gravity waves of the two layer shallow water model, over the sixth-order SCFDM. For the Rossby waves, the sixth-order SCFDM shows improvement, for the barotropic and baroclinic modes, over the sixth-order CCFDM method except on Arakawa's C grid. In the second part of the present work, the sixth-order CCFDM method is used to solve the one-layer and two-layer shallow water equations in their nonlinear form. In one-layer model with periodic boundaries, the performance of the methods for mass conservation is compared. The results show high accuracy of the sixth-order CCFDM method to simulate a complex flow field. Furthermore, to evaluate the performance of the method in a non-periodic domain the sixth-order CCFDM is applied to spatial differencing of vorticity-divergence-mass representation of one-layer shallow water equations to solve a wind-driven current problem with no-slip boundary conditions. The results show good agreement with published works. Finally, the performance of different schemes for spatial differencing of two-layer shallow water equations on Z grid with periodic boundaries is investigated. Results illustrate the high accuracy of combined compact method.

Study of spike and modal stall phenomena in a low-speed axial compressor

This paper presents a study of stall inception mechanisms a in low-speed axial compressor. Previous work has identified two common flow breakdown sequences, the first associated with a short lengthscale disturbance known as a `spike', and the second with a longer lengthscale disturbance known as a `modal oscillation'. In this paper the physical differences between these two mechanisms are illustrated with detailed measurements. Experimental results are also presented which relate the occurrence of the two stalling mechanisms to the operating conditions of the compressor. It is shown that the stability criteria for the two disturbances are different: long lengthscale disturbances are related to a two-dimensional instability of the whole compression system, while short lengthscale disturbances indicate a three-dimensional breakdown of the flow-field associated with high rotor incidence angles. Based on the experimental measurements, a simple model is proposed which explains the type of stall inception pattern observed in a particular compressor. Measurements from a single stage low-speed compressor and from a multistage high-speed compressor are presented in support of the model.

Enhancing the stability of subsonic compressors using casing grooves

Casing grooves are known to increase the stable operating range of axial compressors. The mechanism by which this stability enhancement occurs is poorly understood. This paper develops a better understanding of the behavior of casing grooves through analysis of new data. An experimental parametric study is used to demonstrate the effect of varying the axial location of a single casing groove on the stability and efficiency of the compressor. The effect that the groove has on rotor outflow blockage, blade loading, and the near-casing flow field is then investigated using both experimental and computational methods. It is found that the interaction of the groove with the flow field is different when the groove is positioned forward or aft relative to the blade. The interaction of the groove with the flow in the tip region in both of these positions is presented in detail. Finally, the implications of these findings for the design of casing grooves of different depths are discussed. © 2011 American Society of Mechanical Engineers.

Enhancing the stability of subsonic compressors using casing grooves

Casing grooves are known to increase the stable operating range of axial compressors. The mechanism by which this stability enhancement occurs is poorly understood. This paper develops a better understanding of the behaviour of grooves through analysis of new data. An experimental parametric study is used to demonstrate the effect of varying the axial location of a single casing groove on the stability and efficiency of the compressor. The effect that the groove has on rotor outflow blockage, blade loading and the near-casing flow field is then studied using both experimental and computational methods. It is found that the interaction of the groove with the flow field is different when the groove is positioned forward or aft relative to the blade. The interaction of the groove with the flow in the tip region in both of these positions is presented in detail. Finally, the implications of these findings for the design of casing grooves of different depths are discussed. Copyright © 2009 Rolls-Royce plc.

A Study of Fan-Distortion Interaction Within the NASA Rotor 67 Transonic Stage

The performance of a transonic fan operating within nonuniform inlet flow remains a key concern for the design and operability of a turbofan engine. This paper applies computational methods to improve the understanding of the interaction between a transonic fan and an inlet total pressure distortion. The test case studied is the NASA rotor 67 stage operating with a total pressure distortion covering a 120-deg sector of the inlet flow field. Full-annulus, unsteady, three-dimensional CFD has been used to simulate the test rig installation and the full fan assembly operating with inlet distortion. Novel post-processing methods have been applied to extract the fan performance and features of the interaction between the fan and the nonuniform inflow. The results of the unsteady computations agree well with the measurement data. The local operating condition of the fan at different positions around the annulus has been tracked and analyzed, and this is shown to be highly dependent on the swirl and mass flow redistribution that the rotor induces ahead of it due to the incoming distortion. The upstream flow effects lead to a variation in work input that determines the distortion pattern seen downstream of the fan stage. In addition, the unsteady computations also reveal more complex flow features downstream of the fan stage, which arise due to the three dimensionality of the flow and unsteadiness. © 2012 American Society of Mechanical Engineers.

Stall warning by blade pressure signature analysis

At low mass flow rates axial compressors suffer from flow instabilities leading to stall and surge. The inception process of these instabilities has been widely researched in the past - primarily with the aim of predicting or averting stall onset. In recent times, attention has shifted to conditions well before stall and has focussed on the level of irregularity in the blade passing signature in the rotor tip region. In general, this irregularity increases in intensity as the flow rate through the compressor is reduced. Attempts have been made to develop stall warning/avoidance procedures based on the level of the flow irregularity, but little effort has been made to characterise the irregularity, or to understand its underlying causes. Work on this project has revealed for the first time that the increase in irregularity in the blade passing signature is highly dependent on both tip-clearance and eccentricity. In a compressor with small, uniform, tip-clearance, the increase in blade passing irregularity which accompanies a reduction in flow rate will be modest. If the tip-clearance is enlarged, however, there will be a sharp rise in irregularity at all circumferential locations. In a compressor with eccentric tip-clearance, the increase in irregularity will only occur in the part of the annulus where the tip-clearance is largest, regardless of the average clearance level. In this paper, some attention is also given to the question of whether this irregularity observed in the pre-stall flow field is due to random turbulence, or to some form of coherent flow structure. Detailed flow measurements reveal that the latter is the case. From these findings, it is clear that a stall warning system based on blade passing signature irregularity will not be viable in an aero-engine where tip-clearance size and eccentricity change during each flight cycle and over the life of the compressor. Copyright © 2011 by ASME.

The structure of turbulent stratified and premixed methane/air flames I: Non-swirling flows

A series of flames in a turbulent methane/air stratified swirl burner is presented. The degree of stratification and swirl are systematically varied to generate a matrix of experimental conditions, allowing their separate and combined effects to be investigated. Non-swirling flows are considered in the present paper, and the effects of swirl are considered in a companion paper (Part II). A mean equivalence ratio of φ=0.75 is used, with φ for the highest level of stratification spanning 0.375-1.125. The burner features a central bluff-body to aid flame stabilization, and the influence of the induced recirculation zone is also considered. The current work focuses on non-swirling flows where two-component particle image velocimetry (PIV) measurements are sufficient to characterize the main features of the flow field. Scalar data obtained from Rayleigh/Raman/CO laser induced fluorescence (CO-LIF) line measurements at 103μm resolution allow the behavior of key combustion species-CH 4, CO 2, CO, H 2, H 2O and O 2-to be probed within the instantaneous flame front. Simultaneous cross-planar OH-PLIF is used to determine the orientation of the instantaneous flame normal in the scalar measurement window, allowing gradients in temperature and progress variable to be angle corrected to their three dimensional values. The relationship between curvature and flame thickness is investigated using the OH-PLIF images, as well as the effect of stratification on curvature.The main findings are that the behavior of the key combustion species in temperature space is well captured on the mean by laminar flame calculations regardless of the level of stratification. H 2 and CO are significant exceptions, both appearing at elevated levels in the stratified flames. Values for surface density function and by extension thermal scalar dissipation rate are found to be substantially lower than laminar values, as the thickening of the flame due to turbulence dominates the effect of increased strain. These findings hold for both premixed and stratified flames. The current series of flames is proposed as an interesting if challenging set of test cases for existing and emerging turbulent flame models, and data are available on request. © 2012 The Combustion Institute.

Stall Warning by Blade Pressure Signature Analysis

At low mass flow rates, axial compressors suffer from flow instabilities leading to stall and surge. The inception process of these instabilities has been widely researched in the past---primarily with the aim of predicting or averting stall onset. In recent times, attention has shifted to conditions well before stall and has focused on the level of irregularity in the blade passing signature in the rotor tip region. In general, the irregularity increases in intensity as the flow rate through the compressor is reduced. Attempts have been made to develop stall warning/avoidance procedures based on the level of flow irregularity, but little effort has been made to characterize the irregularity itself, or to understand its underlying cause. Work on this project has revealed for the first time that the increase in irregularity in the blade passing signature is highly dependent on both tip-clearance size and eccentricity. In a compressor with small, uniform, tip-clearance, the increase in blade passing irregularity that accompanies a reduction in flow rate will be modest. If the tip-clearance is enlarged, however, there will be a sharp rise in irregularity at all circumferential locations. In a compressor with eccentric tip-clearance, the increase in irregularity will only occur in the part of the annulus where the tip-clearance is largest, regardless of the average clearance level. In this paper, some attention is also given to the question of whether the irregularity observed in the prestall flow field is due to random turbulence or to some form of coherent flow structure. Detailed flow measurements reveal that the latter is the case. From these findings, it is clear that a stall warning system based on blade passing signature irregularity would be difficult to implement in an aero-engine where tip-clearance size and eccentricity change during each flight cycle and over the life of the compressor. © 2013 American Society of Mechanical Engineers.

Non-reflective boundary conditions for non-uniform mean velocity inlet and outlet flows

The numerical solution of problems in unbounded physical space requires a truncation of the computational domain to a reasonable size. As a result, the conditions on the artificial boundaries are generally unknown. Assumptions like constant pressure or velocities are only valid in the far field and lead to spurious reflections if applied on the boundaries of the truncated domain. A number of attempts have been made over the past decades to design conditions that prevent such reflections. One approach is based on characteristics. The standard analysis assumes a spatially uniform mean flow field but this is often impractical. In the present paper we show how to extend the formulation to the more general case of a non-uniform mean velocity field. A number of test cases are provided and our results compare favourably with other boundary conditions. In principle the present approach can be extended to include non-uniformities in all variables.

Droplet Generator

A method and device for periodically perturbing the flow field within a microfluidic device to provide regular droplet formation at high speed.

Can fundamental shock-wave/boundary-layer interaction research be relevant to inlet aerodynamics?

The important influence of shock waves on supersonic inlet performance has led to much time and effort being expended in the area of shock wave/boundary layer interaction research (SWBLI) and SWBLI control. In this short review, the impact of SWBLIs on supersonic inlet aerodynamic research is discussed and is contrasted with fundamental SWBLI research. Inlet research focussed on internal flow performance is reviewed, based on the salient results, conclusions, and the limitations of such work. The role of fundamental SWBLI research in relation to supersonic inlet research is considered, and the possible positive impact of improving the link between fundamental SWBLI research and inlet design is considered. A simple flow-field is discussed which is thought to be able to simulate at least some more of the flow physics found in a typical inlet. A brief review of real inlet parameters is then given to help determine appropriate fundamental experimental parameters such as incoming Mach number, incoming boundary-layer thickness and subsonic difiuser angle. Copyright © 2012 by N. Titchener, H. Babinsky, and E. Loth.