Micro-ramp control for oblique shock wave / boundary layer interactions

Supersonic engine intakes operating supercritically feature shock wave / boundary layer interactions (SBLIs), which are conventionally controlled using boundary layer bleed. The momentum loss of bleed flow causes high drag, compromising intake performance. Micro-ramp sub-boundary layer vortex generators (SBVGs) have been proposed as an alternative form of flow control for oblique SBLIs in order to reduce the bleed requirement. Experiments have been conducted at Mach 2.5 to characterise the flow details on such devices and investigate their ability to control the interaction between an oblique shock wave and the naturally grown turbulent boundary layer on the tunnel floor. Micro-ramps of four sizes with heights ranging from 25% to 75% of the uncontrolled boundary layer thickness were tested. The flow over all sizes of microramp was found to be similar, featuring streamwise counter-rotating vortices which entrain high momentum fluid, locally reducing the boundary layer displacement thickness. When installed ahead of the shock interaction it was found that the positioning of the micro-ramps is of limited importance. Micro-ramps did not eliminate flow separation. However, the previously two-dimensional separation was broken up into periodic three-dimensional separation zones. The interaction length was reduced and the pressure gradient across the interaction was increased.

Predicting shock response in uncertain structures using the Hybrid method

Using the Hybrid method (FE + SEA) it is possible to estimate the frequency response of an uncertain structure. The current work develops the Hybrid method to allow for time domain analysis of the shock response of a structure. Problems to be overcome when taking Hybrid method results into the time domain are a) the Hybrid method frequency response has no phase information, and b) the Hybrid method frequency response is smoothed in frequency and shows no modal peaks. In this paper the first problem has been overcome, using minimum phase reconstruction. Explanation of minimum phase reconstruction and its limitations are described, and application to shock problems described. © 2009 IOP Publishing Ltd.

Contrasting high-resolution characteristic shock-capturing methods in aeroacoustic test problems

Two shock-capturing methods are considered. One is based on a standard conservative Roe scheme with van Leer's MUSCL variable extrapolation method applied to characteristic variables and a Runge-Kutta time stepping scheme. The other is based on the novel CABARET space-time scheme, which uses two sets of staggered variables, one for the conservation step and the other for characteristic splitting into local Riemann invariants. The methods are compared in a range of 2-D inviscid compressible flow test cases. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

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.

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.

An experimental and numerical study of an oscillating transonic shock wave in a duct

A combined experimental and numerical study of a transonic shock wave in a parallel walled duct subject to downstream pressure perturbations has been conducted. Experiments and simulations have been carried out with a shock strength of M∞ = 1.4 for pressure perturbation frequencies in the range 16-90 Hz. The dynamics of unsteady shock motion and the interaction structure between the unsteady transonic shock wave and the turbulent tunnel floor boundary layer have been investigated. It is found that the (experimentally measured) dynamics of shock motion are generally well predicted by the computational scheme, especially at relatively low (≈ 40 Hz) frequencies. However, at higher frequencies (≈ 90 Hz), some subtle differences between the shock dynamics measured in experiments and those predicted by Computational Fluid Dynamics (CFD) exist. There is evidence from experiments that variations in shock / boundary layer interaction (SBLI) structure caused by shock motion are responsible for a change in the nature of shock dynamics between low and high frequency. In contrast, numerical results at low and high frequencies do not differ significantly and this suggests that the numerical method is not fully capturing the physics of the unsteady flow. Possible reasons for this are considered and a number of areas where CFD is unable to replicate experimental observations are identified. Significantly, CFD predicts changes in SBLI structure due to shock motion that are much too large and this may explain why none of the subtle effects on shock dynamics seen in experiments occur in CFD. Further work developing numerical methods that demonstrate a more realistic sensitivity of SBLI structure to unsteady shock motion is required. Copyright © 2010 by P.J.K. Bruce.

Experimental investigation into parameters governing corner interactions for transonic shock wave/ boundary layer interactions

Experiments were conducted investigating the interaction between a normal shock wave and a corner boundary layer in a constant area rectangular duct. Active corner suction and passive blowing were applied to manipulate the natural corner flows developing in the working section of the Cambridge University supersonic wind tunnel. In addition robust vane micro-vortex generators were applied to the corners of the working section. Experiments were conducted at Mach numbers of M∞=1.4 and 1.5. Flow visualisation was carried out through schlieren and surface oil flow, while static pressures were recorded via floor tappings. The results indicate that an interplay occurs between the corner flow and the centre line flow. It is believed that corner flow separation acts to induce a shock bifurcation, which in turn leads to a smearing of the adverse pressure gradient elsewhere. In addition the blockage effect from the corners was seen to result in a reacceleration of the subsonic post-shock flow. As a result manipulation of the corner regions allows a separated or attached centre line flow to be observed at the same Mach number. Copyright © 2010 by Babinsky, Burton, Bruce.

Microvortex generators applied to flowfield containing a normal shock wave and diffuser

The flow through a terminating shock wave and the subsequent subsonic diffuser typically found in supersonic inlets has been simulated using a small-scale wind tunnel. Experiments have been conducted at an inflow Mach number of 1.4 using a dual-channel working section to produce a steady near-normal shock wave. The setup was designed so that the location of the shock wave could be varied relative to the diffuser. As the near-normal shock wave was moved downstream and into the diffuser, an increasingly distorted, three-dimensional, and separated flow was observed. Compared with the interaction of a normal shock wave in a constant area duct, the addition of the diffuser resulted in more prominent corner interactions. Microvortex generators were added to determine their potential for removing flow separation. Although these devices were found to reduce the extent of separation, they significantly increased three-dimensionality and even led to a large degree of flow asymmetry in some configurations. Copyright © 2011 by Neil Titchener and Holger Babinsky.

CFD-DEM simulation of propagation of sound waves in fluid particles fluidised medium

In this work, speed of sound in 2 phase mixture has been explored using CFD-DEM (Computational Fluid Dynamcis - Discrete Element Modelling). In this method volume averaged Navier Stokes, continuity and energy equations are solved for fluid. Particles are simulated as individual entities; their behaviour is captured by Newton's laws of motion and classical contact mechanics. Particle-fluid interaction is captured using drag laws given in literature.The speed of sound in a medium depends on physical properties. It has been found experimentally that speed of sound drops significantly in 2 phase mixture of fluidised particles because of its increased density relative to gas while maintaining its compressibility. Due to the high rate of heat transfer within 2 phase medium as given in Roy et al. (1990), it has been assumed that the fluidised gas-particle medium is isothermal.The similar phenomenon has been tried to be captured using CFD-DEM numerical simulation. The disturbance is introduced and fundamental frequency in the medium is noted to measure the speed of sound for e.g. organ pipe. It has been found that speed of sound is in agreement with the relationship given in Roy et al. (1990). Their assumption that the system is isothermal also appears to be valid.

Shock transmission and eigenvalue veering within a coupled system.

The operation of dynamical systems in harsh environments requires continuous monitoring. Internal sensors may be used to monitor the conditions in real time. A typical example is the sensor and electronic components used in space structures which, especially during launch, are subject to huge g force. The paper will present an experimental and theoretical study on a simplified model used to analyze the possible cause of high acceleration on the enclosed sensors and equipments due to impulsive loading. The model system consists of two beams coupled using compliant connections. An impulse hammer excites one beam, and vibrations are transmitted to the indirectly driven beam. A theoretical model is developed using a Rayleigh-Ritz approach and validated using experimental results in both the frequency and time domains. Monto Carlo simulation was done with random masses positioned on the indirectly driven beam to determine the worst-case conditions for maximum peak acceleration. Highest acceleration levels were found when mode matching in the two beams led to veering behavior in the coupled modes. The results suggest guidelines for the detailed design of internal components of a structure exposed to shock loading from its environment. [The authors thank Schlumberger Cambridge Research for financial support.].

Baroclinic waves in an air-filled thermally driven rotating annulus.

In this study an experimental investigation of baroclinic waves in air in a differentially heated rotating annulus is presented. Air has a Prandtl number of 0.707, which falls within a previously unexplored region of parameter space for baroclinic instability. The flow regimes encountered include steady waves, periodic amplitude vacillations, modulated amplitude vacillations, and either monochromatic or mixed wave number weak waves, the latter being characterized by having amplitudes less than 5% of the applied temperature contrast. The distribution of these flow regimes in parameter space are presented in a regime diagram. It was found that the progression of transitions between different regimes is, as predicted by recent numerical modeling results, in the opposite sense to that usually found in experiments with high Prandtl number liquids. No hysteresis in the flow type, with respect to variations in the rotation rate, was found in this investigation.

Synchronization of modulated traveling baroclinic waves in a periodically forced, rotating fluid annulus.

Frequency entrainment and nonlinear synchronization are commonly observed between simple oscillatory systems, but their occurrence and behavior in continuum fluid systems are much less well understood. Motivated by possible applications to geophysical fluid systems, such as in atmospheric circulation and climate dynamics, we have carried out an experimental study of the interaction of fully developed baroclinic instability in a differentially heated, rotating fluid annulus with an externally imposed periodic modulation of the thermal boundary conditions. In quasiperiodic and chaotic amplitude-modulated traveling wave regimes, the results demonstrate a strong interaction between the natural periodic modulation of the wave amplitude and the externally imposed forcing. This leads to partial or complete phase synchronization. Synchronization effects were observed even with very weak amplitudes of forcing, and were found with both 1:1 and 1:2 frequency ratios between forcing and natural oscillations.