UEGO based measurement of EGR rate and residual gas fraction

A technique is presented for measuring the exhaust gas recirculation (EGR) and residual gas fraction (RGF) using a fast UEGO based O2 measurement of the manifold or in-cylinder gases, and of the exhaust gases. The technique has some advantages over the more common CO2-based method. In the case of an RGF measurement, fuel interference must be eliminated and special fuelling arrangements are is required. It is shown how a UEGO-based measurement, though sensitive to reactive species in the exhaust (such as H 2), as a system reports EGR/ RGF rates faithfully. Preliminary tests showed that EGR and RGF measurements using the O2 approach agreed well with CO2-based measurements. © 2011 SAE International.

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.

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

This paper investigates the use of inertial actuators to reduce the sound radiated by a submarine hull under excitation from the propeller. The axial forces from the propeller are tonal at the blade passing frequency. The hull is modeled as a fluid-loaded cylindrical shell with ring stiffeners and equally spaced bulkheads. The cylinder is closed at each end by circular plates and conical end caps. The forces from the propeller are transmitted to the hull by a rigid foundation connected to the propeller shaft. Inertial actuators are used as the structural control inputs. The actuators are arranged in circumferential arrays and attached to the internal end plates of the hull. Two active control techniques corresponding to active vibration control and discrete structural acoustic sensing are implemented to attenuate the structural and acoustic responses of the submarine. In the latter technique, error information on the radiated sound fields is provided by a discrete structural acoustic sensor. An acoustic transfer function is defined to estimate the far field sound pressure from a single point measurement on the hull. The inertial actuators are shown to provide control forces with a magnitude large enough to reduce the sound due to hull vibration. © 2012 American Society of Mechanical Engineers.

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.

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.

Oscillatory sensitivity patterns for global modes in wakes

Globally unstable wakes with co-flow at intermediate Reynolds numbers are studied, to quantify important spatial regions for the development and control of the global instability. One region of high structural sensitivity is found close to the inlet for all wakes, in agreement with previous findings for cylinder wakes. A second, elongated region of high structural sensitivity is seen downstream of the first one for unconfined wakes at Re = 400. When base flow modifications are considered, a spatially oscillating sensitivity pattern is found inside the downstream high structural sensitivity region. This implies that the same change in the base flow can either destabilize or stabilize the flow, depending on the exact position where it is applied. It is shown that the sensitivity pattern remains unchanged for different choices of streamwise boundary conditions and numerical resolution. Actual base flow modifications are applied in selected configurations, and the linear global modes recomputed. It is confirmed that the linear global eigenvalues move according to the predicted sensitivity pattern for small amplitude base flow modifications, for which the theory applies. We also look at the implications of a small control cylinder on the flow. Only the upstream high sensitivity region proves to be robust in terms of control, but one should be careful not to disturb the flow in the downstream high sensitivity region, in order to achieve control. The findings can have direct implications on the numerical resolution requirements for wakes at higher Reynolds numbers. Furthermore, they provide one more possible explanation to why confined wakes have a more narrow frequency spectrum than unconfined wakes.

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.

Hybrid passive/active inertial actuators to attenuate the structural response of a submerged hull

This paper theoretically investigates the application of tuned vibration absorbers and hybrid passive/active inertial actuators to reduce the vibrational responses of plates and shells. The passive/active actuators are initially applied to a simple plate. A model of a submerged hull consisting of a ring stiffened finite cylinder with bulkheads and external fluid loading is then considered. The fluctuating forces from the propeller result in excitation of the low frequency global hull modes. Inertial actuators and tuned vibration absorbers are located at each end of the hull and in circumferential arrays to reduce the hull structural response at its axial resonances. The control performance of the hybrid passive/active inertial actuator, where the passive component is tuned to a structural resonance, is compared to the attenuation achieved by a fully passive tuned vibration absorber. This work shows the potential of using hybrid passive/active inertial actuators to attenuate the global structural responses of a submerged vessel.

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.