Shock / boundary-layer interaction control using three-dimensional bumps for transonic wings

Three-dimensional bumps have been developed and investigated, aiming at the two major objectives of shock-wave / boundary-layer interaction control, i.e. drag reduction and suppression of separation, simultaneously. An experimental investigation has been conducted for a default rounded bump in channel now at University of Cambridge and a computational study has been performed for a spanwise series of rounded bumps mounted on a transonic aerofoil at University of Stuttgart. Observed in both cases are wave drag reduction owing to A-shock structures produced by three-dimensional surface bumps and mild control effects on the boundary layer. The effects of rough surface and tall extension have been investigated as well as several geometric variations and multiple bump configurations. A double configuration of narrow rounded bumps has been found to best perform amongst the tested, considerably reducing wave drag through a well-established A-shock structure with little viscous penalty and thus achieving substantial overall drag reduction. Counter-rotating streamwise vortex pairs have been produced by some configurations as a result of local flow separation, but they have been observed to be confined in relatively narrow wake regions, expected to be beneficial in suppressing large-scale separation under off-design condition despite increase of viscous drag. On the whole a large potential of three-dimensional control with discrete rounded bumps has been demonstrated both experimentally and numerically, and experimental investigation of bumps fitted on a transonic aerofoil or wing is suggested toward practical application.

Computational fluid-structure interaction methods for simulation of inflatable aerodynamic decelerators

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. It is theorized that volume-mass characteristics of these decelerators are superior to those of common supersonic/subsonic parachutes and after deployment they may suffer no instabilities at high Mach numbers. A high fidelity computational fluid-structure interaction model is employed to investigate the behavior of tension cone inflatable aeroshells at supersonic speeds up to Mach 2.0. The computational framework targets the large displacements regime encountered during the inflation of the decelerator using fast level set techniques to incorporate boundary conditions of the moving structure. The preliminary results indicate large but steady aeroshell displacement with rich dynamics, including buckling of the inflatable torus that maintains the decelerator open under normal operational conditions, owing to interactions with the turbulent wake. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

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.

RANS-URANS in axial compressor, a design methodology

Turbomachinery flows are inherently unsteady. Until now during the design process, unsteadiness has been neglected, with resort merely to steady numerical simulations. Despite the assumption involved, the results obtained with steady simulations have been used with success. One of the questions arising in recent years is can unsteady simulations be used to improve the design of turbomachines? In this work the numerical simulation of a multi-stage axial compressor is carried out. Comparison of Reynolds averaged Navier-Stokes (RANS) and unsteady Reynolds averaged Navier-Stokes (URANS) calculation shows that the unsteadiness affects pressure losses and the prediction of stall limit. The unsteady inflow due to the wake passing mainly modifies the losses and whirl angle near the endwalls. The computational cost of the fully unsteady compared with a steady simulation is about four times in terms of mesh dimension and two orders of magnitude as number of iterations. A mixed RANS-URANS solution has been proposed to give the designer the possibility to simulate an unsteady stage embedded in a steady-state simulation. This method has been applied to the simulation of a four-stage axial compressor rig. The mixed RANS-URANS approach has been developed using sliding and mixing planes as interface conditions. The rotor-stator interaction has been captured physically while reducing the computational time and mesh size.

Active control of sound radiated by a submarine hull in bending vibration using inertial actuators

Submarines are efficient sources of low frequency radiated noise due to the vibrations induced by the rotation of the propeller in a non uniform wake. In this work the possibility of using inertial actuators to reduce the far field sound pressure is investigated. The submerged vessel is modelled as a cylindrical shell with two conical end caps. Complicating effects such as ring stiffeners, bulkheads and the fluid loading are taken into account. A harmonic radial force is transmitted from the propeller to the hull through the stern end cone and it is tonal at the blade passing frequency (rotational speed of the shaft multiplied by the number of blades). The actuators are attached at the inside of the prow end cone to form a circumferential array. Both Active Vibration Control (AVC) and Active Structural Acoustic Control (ASAC) are analysed and it is shown that the inertial actuators can significantly reduce the far field sound pressure.

Global linear and nonlinear stability of viscous confined plane wakes with co-flow

The global stability of confined uniform density wakes is studied numerically, using two-dimensional linear global modes and nonlinear direct numerical simulations. The wake inflow velocity is varied between different amounts of co-flow (base bleed). In accordance with previous studies, we find that the frequencies of both the most unstable linear and the saturated nonlinear global mode increase with confinement. For wake Reynolds number Re = 100 we find the confinement to be stabilising, decreasing the growth rate of the linear and the saturation amplitude of the nonlinear modes. The dampening effect is connected to the streamwise development of the base flow, and decreases for more parallel flows at higher Re. The linear analysis reveals that the critical wake velocities are almost identical for unconfined and confined wakes at Re ≈ 400. Further, the results are compared with literature data for an inviscid parallel wake. The confined wake is found to be more stable than its inviscid counterpart, whereas the unconfined wake is more unstable than the inviscid wake. The main reason for both is the base flow development. A detailed comparison of the linear and nonlinear results reveals that the most unstable linear global mode gives in all cases an excellent prediction of the initial nonlinear behaviour and therefore the stability boundary. However, the nonlinear saturated state is different, mainly for higher Re. For Re = 100, the saturated frequency differs less than 5% from the linear frequency, and trends regarding confinement observed in the linear analysis are confirmed.

The local and global stability of confined planar wakes at intermediate Reynolds number

At high Reynolds numbers, wake flows become more globally unstable when they are confined within a duct or between two flat plates. At Reynolds numbers around 100, however, global analyses suggest that such flows become more stable when confined, while local analyses suggest that they become more unstable. The aim of this paper is to resolve this apparent contradiction by examining a set of obstacle-free wakes. In this theoretical and numerical study, we combine global and local stability analyses of planar wake flows at $\mathit{Re}= 100$ to determine the effect of confinement. We find that confinement acts in three ways: it modifies the length of the recirculation zone if one exists, it brings the boundary layers closer to the shear layers, and it can make the flow more locally absolutely unstable. Depending on the flow parameters, these effects work with or against each other to destabilize or stabilize the flow. In wake flows at $\mathit{Re}= 100$ with free-slip boundaries, flows are most globally unstable when the outer flows are 50 % wider than the half-width of the inner flow because the first and third effects work together. In wake flows at $\mathit{Re}= 100$ with no-slip boundaries, confinement has little overall effect when the flows are weakly confined because the first two effects work against the third. Confinement has a strong stabilizing effect, however, when the flows are strongly confined because all three effects work together. By combining local and global analyses, we have been able to isolate these three effects and resolve the apparent contradictions in previous work.

Computational fluid-structure interaction of DGB parachutes in compressible fluid flow

This paper describes large-scale simulations of compressible flows over a supersonic disk-gap-band parachute system. An adaptive mesh refinement method is used to resolve the coupled fluid-structure model. The fluid model employs large-eddy simulation to describe the turbulent wakes appearing upstream and downstream of the parachute canopy and the structural model employed a thin-shell finite element solver that allows large canopy deformations by using subdivision finite elements. The fluid-structure interaction is described by a variant of the Ghost-Fluid method. The simulation was carried out at Mach number 1.96 where strong nonlinear coupling between the system of bow shocks, turbulent wake and canopy is observed. It was found that the canopy oscillations were characterized by a breathing type motion due to the strong interaction of the turbulent wake and bow shock upstream of the flexible canopy. Copyright © 2010 by ASME.

Cinematographic OH-PLIF measurements of two interacting turbulent premixed flames with and without acoustic forcing

This paper describes an experimental investigation into the interactions that occur between two lean turbulent premixed flames stabilised on conical bluff-bodies when they are moved closer together. Cinematographic OH-PLIF measurements were acquired to investigate adjacent flame front interactions as a function of flame separation distance (S). Flame surface density (FSD) and curvature were determined to characterise the unforced flames. Acoustic forcing was then applied to explore the amplitude dependent thermo-acoustic response. Phase-averaged FSD and global heat release measurements in the form of OH * chemiluminescence were obtained for a range of forcing frequencies (f) and amplitudes (A) as a function of S. As the flames were brought closer together the adjacent annular jets were found to merge into a single jet structure. This caused adjacent flame fronts to merge above the wake region between the two flames at a location determined by the jet efflux (flame angle) and S. This region of flame-flame interaction we refer to as 'interacting region'. In the unforced flames, a trend of increasingly negative curvature for decreasing S produced a small net increase in flame surface area via cusp formation. When subjected to acoustic forcing, S-dependent regimes were found in the global heat release response as a function A. The overall trend showed that the occurrence of jet/flame merging reduces the value of A at which non-linear response occurs. In support of previous findings for flames stabilised along shear layers, the phase-averaged FSD showed that the flame dynamics that drive the thermo-acoustic response result from the roll-up of vortices which generate large-scale vortex-flame interactions. Compared with axisymmetric flames, the occurrence of jet merging alters the vortex-flame interactions resulting in an asymmetric contribution to the heat release between the wall and interacting regions. The majority of the heat release was found to occur in the interacting region through the rapid production and destruction of flame surface area. The occurrence of jet merging and large-scale interactions between adjacent flames result in different physical mechanisms that drive the thermo-acoustic response compared with single axisymmetric flames. © 2011.

Experimental study of homogenisation in grid turbulence

Turbulence statistics have been measured immediately downstream of a regular grid made of round rods with rod spacing M. 2D-2C PIV was used to analyse a measurement area of 14M x 4M in the down and cross-stream directions respectively. The relevant Reynolds number span the range Re M = U ∞M/ν = 5 500 - 16 500. The Reynolds shear stresses recorded on two parallel measurement planes differently located relative to the grid exhibit significant discrepancies over the first 5M, but have completely homogenised in the cross-stream direction by x/M = 7. The downstream evolution of the two-point velocity correlation functions shows a progressive loss of coherence and a clear trend towards the expected isotropic behavior. The same conclusions apply to measurements taken in the wake of another regular grid made of square rods. Changes in the vortex shedding pattern from the grid were observed at the lowest Reynolds number, with two of the four rod wakes captured shedding in phase with each other but in anti-phase with a third one. The impact of this early flow coherence on the turbulence statistics did not persist due to the homogenisation process.

Thermal imaging study of scalar transport in shallow wakes

The thermal imaging technique relies on the usage of infrared signal to detect the temperature field. Using temperature as a flow tracer, thermography is used to investigate the scalar transport in the shallow-water wake generated by an emergent circular cylinder. Thermal imaging is demonstrated to be a good quantitative flow visualization technique for studying turbulent mixing phenomena in shallow waters. A key advantage of the thermal imaging method over other scalar measurement techniques, such as the Laser Induced Fluorescence (LIF) and Planar Concentration Analysis (PCA) methods, is that it involves a very simple experimental setup. The dispersion characteristics captured with this technique are found to be similar to past studies with traditional measurement techniques. © 2012 Publishing House for Journal of Hydrodynamics.

2D PIV measurements in the near field of grid turbulence using stitched fields from multiple cameras

We present measurements of grid turbulence using 2D particle image velocimetry taken immediately downstream from the grid at a Reynolds number of Re M = 16500 where M is the rod spacing. A long field of view of 14M x 4M in the down- and cross-stream directions was achieved by stitching multiple cameras together. Two uniform biplanar grids were selected to have the same M and pressure drop but different rod diameter D and crosssection. A large data set (10 4 vector fields) was obtained to ensure good convergence of second-order statistics. Estimations of the dissipation rate ε of turbulent kinetic energy (TKE) were found to be sensitive to the number of meansquared velocity gradient terms included and not whether the turbulence was assumed to adhere to isotropy or axisymmetry. The resolution dependency of different turbulence statistics was assessed with a procedure that does not rely on the dissipation scale η. The streamwise evolution of the TKE components and ε was found to collapse across grids when the rod diameter was included in the normalisation. We argue that this should be the case between all regular grids when the other relevant dimensionless quantities are matched and the flow has become homogeneous across the stream. Two-point space correlation functions at x/M = 1 show evidence of complex wake interactions which exhibit a strong Reynolds number dependence. However, these changes in initial conditions disappear indicating rapid cross-stream homogenisation. On the other hand, isotropy was, as expected, not found to be established by x/M = 12 for any case studied. © Springer-Verlag 2012.

Study of the turbulence in the air-side and water-side boundary layers in experimental laboratory wind induced surface waves

This study detailed the structure of turbulence in the air-side and water-side boundary layers in wind-induced surface waves. Inside the air boundary layer, the kurtosis is always greater than 3 (the value for normal distribution) for both horizontal and vertical velocity fluctuations. The skewness for the horizontal velocity is negative, but the skewness for the vertical velocity is always positive. On the water side, the kurtosis is always greater than 3, and the skewness is slightly negative for the horizontal velocity and slightly positive for the vertical velocity. The statistics of the angle between the instantaneous vertical fluctuation and the instantaneous horizontal velocity in the air is similar to those obtained over solid walls. Measurements in water show a large variance, and the peak is biased towards negative angles. In the quadrant analysis, the contribution of quadrants Q2 and Q4 is dominant on both the air side and the water side. The non-dimensional relative contributions and the concentration match fairly well near the interface. Sweeps in the air side (belonging to quadrant Q4) act directly on the interface and exert pressure fluctuations, which, in addition to the tangential stress and form drag, lead to the growth of the waves. The water drops detached from the crest and accelerated by the wind can play a major role in transferring momentum and in enhancing the turbulence level in the water side.On the air side, the Reynolds stress tensor's principal axes are not collinear with the strain rate tensor, and show an angle α σ≈=-20°to-25°. On the water side, the angle is α σ≈=-40°to-45°. The ratio between the maximum and the minimum principal stresses is σ a/σ b=3to4 on the air side, and σ a/σ b=1.5to3 on the water side. In this respect, the air-side flow behaves like a classical boundary layer on a solid wall, while the water-side flow resembles a wake. The frequency of bursting on the water side increases significantly along the flow, which can be attributed to micro-breaking effects - expected to be more frequent at larger fetches. © 2012 Elsevier B.V.

A Physical Interpretation of Stagnation Pressure and Enthalpy Changes in Unsteady Flow

This paper provides a physical interpretation of the mechanism of stagnation enthalpy and stagnation pressure changes in turbomachines due to unsteady flow, the agency for all work transfer between a turbomachine and an inviscid fluid. Examples are first given to illustrate the direct link between the time variation of static pressure seen by a given fluid particle and the rate of change of stagnation enthalpy for that particle. These include absolute stagnation temperature rises in turbine rotor tip leakage flow, wake transport through downstream blade rows, and effects of wake phasing on compressor work input. Fluid dynamic situations are then constructed to explain the effect of unsteadiness, including a physical interpretation of how stagnation pressure variations are created by temporal variations in static pressure; in this it is shown that the unsteady static pressure plays the role of a time-dependent body force potential. It is further shown that when the unsteadiness is due to a spatial nonuniformity translating at constant speed, as in a turbomachine, the unsteady pressure variation can be viewed as a local power input per unit mass from this body force to the fluid particle instantaneously at that point. © 2012 American Society of Mechanical Engineers.

Comparison of bleed and micro-vortex generator effects on supersonic boundary-layers

An experimental study of bleed and vortex generators in supersonic ow has been conducted. Methods were developed to analyze and directly compare the two systems' effects on turbulent boundary layers to better understand their potential to mitigate ow separation. LDA was used to measure two components of velocity in the boundary-layer for three cases|baseline, with bleed, or with a VG|at Mach numbers of 1.3, 1.5 and 1.8. The bleed system was comprised of a series of 2mm diameter normal holes operated at different suction rates, removing up to 10% of the incoming boundary layer. Three VG shapes were tested only at Mach 1.5 and 1.8. Measurements of the evolution of Hi and Cf downstream of each device indicate that Hi is not an appropriate parameter to gauge the effectiveness of vortex generators due to boundary layer wake distortion. The skin friction coeficient Cf may be a more appropriate measure. Similar increases in Cf were generated by VGs and bleed. The recovery to baseline conditions downstream of bleed was sensitive to Mach number, and more investigation of that effect will be required. Copyright © 2012 by University of Cambridge.