Review of computational aeroacoustics for application in electronics cooler noise

A review of computational aeroacoustics (CCA) was made for application in electronics cooler noise. Computational aeroacoustics encompasses all numerical methods where the purposes is to predict the noise emissions from a simulated flow. Numerical simulation of the flow inside and around heat sinks and fans can lead to a prediction of the emitted noise while they are still in the design phase. Direct CCA is theoretically the best way to predict flow-based acoustic phenomena numerically. It is typically used only for low-frequency sound prediction. The boundary element method offers low computational cost and does not use a computational grid, but instead use vortex-surface calculations to determine tonal noise. Axial fans are commonly used to increase the airflow and thus the heat transfer over the heat sinks within the computer cases. Very detailed source simulations in the fan and heat sink region coupled with the use of analogy methods could result in excellent simulation results with a reasonable computational effort.

High-order DES and zonal LES of a labyrinth seal

Hybrid numerical large eddy simulation (NLES) and detached eddy simulation (DES) methods are assessed on a labyrinth seal geometry. A high sixth order discretization scheme is used and is validated using a test case of a two dimensional vortex. The hybrid approach adopts a new blending function and along with DES is initially validated using a simple cavity flow. The NLES method is also validated outside of RANS zones. It is found that there is very little resolved turbulence in the cavity for the DES simulation. For the labyrinth seal calculations the DES approach is problematic giving virtually no resolved turbulence content. It is seen that over the tooth tips the extent of the LES region is small and is likely to be a strong contributor to excessive flow damping in these regions. On the other hand the zonal Hamilton-Jacobi approach did not suffer from this trait. In both cases the meshes used are considered to be hybrid RANS-LES adequate. Fortunately (or perhaps unfortunately) the DES profiles are in agreement with the time mean experimental measurements. It is concluded that for an inexperienced CFD practitioner this could have wider implications particularly if transient results such as unsteady loading are desired. Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc.

Response of two acoustically excited turbulent premixed flames to an imposed phase lag

This paper describes an experimental investigation into the interactions that occur between two acoustically forced lean turbulent premixed flames for an induced phase lag. Phase-averaged FSD from cinematographic OH-PLIF measurements and global heat release measurements were obtained for a range phase lags (ψ S) and amplitudes (A) as a function of flame separation distance, S. The effect of bringing two flames closer together causes jet merging, which alters the vortex flame interactions that drive the thermo-acoustic response. To simulate circumferential modes a phase lag was introduced, which affected the flame dynamics in the region of flame-flame interaction, with these changes dependent on S. For moderate separation distances, the flame structure becomes increasingly asymmetric inducing a very small transverse oscillation. However, for moderate phase lags (φ s ≤ 20) the magnitude of these changes and their subsequent influence on the thermo-acoustic response was found to be slight in comparison with changes in S. Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc.

Winglets for Improved Aerothermal Performance of High Pressure Turbines

This paper investigates the design of winglet tips for unshrouded high pressure turbine rotors, considering aerodynamic and thermal performance simultaneously. A novel parameterization method has been developed to alter the tip geometry of a rotor blade. A design survey of un-cooled, flat-tipped winglets is performed using RANS calculations for a single rotor at engine representative operating conditions. Compared to a plain tip, large efficiency gains can be realized by employing an overhang around the full perimeter of the blade, but the overall heat load rises significantly. By employing an overhang on only the early suction surface, significant efficiency improvements can be obtained without increasing the overall heat transfer to the blade. The flow physics are explored in detail to explain the results. For a plain tip, the leakage and passage vortices interact to create a three-dimensional impingement onto the blade suction surface, causing high heat transfer. The addition of an overhang on the early suction surface displaces the tip leakage vortex away from the blade, weakening the impingement effect and reducing the heat transfer on the blade. The winglets reduce the aerodynamic losses by unloading the tip section, reducing the leakage flow rate, turning the leakage flow in a more streamwise direction and reducing the interaction between the leakage fluid and endwall flows. Generally these effects are most effective close to the leading edge of the tip, where the leakage flow is subsonic.

Parallels between wind and crowd loading of bridges

Parallels between the dynamic response of flexible bridges under the action of wind and under the forces induced by crowds allow each field to inform the other.Wind-induced behaviour has been traditionally classified into categories such as flutter, galloping, vortex-induced vibration and buffeting. However, computational advances such as the vortex particle method have led to a more general picture where effects may occur simultaneously and interact, such that the simple semantic demarcations break down. Similarly, the modelling of individual pedestrians has progressed the understanding of human–structure interaction, particularly for large amplitude lateral oscillations under crowd loading. In this paper, guided by the interaction of flutter and vortexinduced vibration in wind engineering, a framework is presented, which allows various human–structure interaction effects to coexist and interact, thereby providing a possible synthesis of previously disparate experimental and theoretical results.

Study of thermal effects in bulk RE-BCO superconductors submitted to a variable magnetic field

When bulk RE-BCO superconductors are used as permanent magnets in engineering applications, they are likely to experience transient variations of the applied magnetic field. The resulting vortex motion may cause a significant temperature increase. As a consequence the initial trapped flux is reduced. In the present work, we first focus on the cause of a temperature increase. The temperature distribution within a superconducting finite cylinder subjected to an alternating magnetic field is theoretically predicted. Results are compared to experimental data obtained by two temperature sensors attached to a bulk YBCO pellet. Second, we consider curative methods for reducing the effect of heat flux on the temperature increase. Hall-probe mappings on YBCO samples maintained out of the thermal equilibrium are performed for two different morphologies : a plain single domain and a single domain with a regularly spaced hole array. The drilled single-domain displays a trapped induction which is weakly affected by the local heating while displaying a high trapped field. © 2006 IOP Publishing Ltd.

Seafloor-riser interaction model

Fatigue stresses associated with extreme storms, vessel movements, and vortex-induced vibrations are critical to the performance of steel catenary risers. The critical location for fatigue damage often occurs within the touchdown zone, where cyclic interaction of the riser with the seabed occurs. Developing a model for seabed stiffness requires characterization of a number of complex nonlinear processes including trench formation, nonlinear soil stiffness, soil suction, and breakaway of the riser from the seafloor. The analytical framework utilized in this research considers the riser-seafloor interaction problem in terms of a pipe resting on a bed of springs, the stiffness characteristics of which are described by nonlinear load-deflection (P-y) curves. The P-y model allows for first penetration and uplift, as well as repenetration and small range motions within the bounding loop defined by extreme loading. The backbone curve is constructed from knowledge of the soil strength, the rate of strength increase with depth, trench width, and two additional parameters, while three parameters are necessary for the cyclic response. © ASCE 2009.

Absolute and convective instability in gas turbine fuel injectors

Hydrodynamic instabilities in gas turbine fuel injectors help to mix the fuel and air but can sometimes lock into acoustic oscillations and contribute to thermoacoustic instability. This paper describes a linear stability analysis that predicts the frequencies and strengths of hydrodynamic instabilities and identifies the regions of the flow that cause them. It distinguishes between convective instabilities, which grow in time but are convected away by the flow, and absolute instabilities, which grow in time without being convected away. Convectively unstable flows amplify external perturbations, while absolutely unstable flows also oscillate at intrinsic frequencies. As an input, this analysis requires velocity and density fields, either from a steady but unstable solution to the Navier-Stokes equations, or from time-averaged numerical simulations. In the former case, the analysis is a predictive tool. In the latter case, it is a diagnostic tool. This technique is applied to three flows: a swirling wake at Re = 400, a single stream swirling fuel injector at Re - 106, and a lean premixed gas turbine injector with five swirling streams at Re - 106. Its application to the swirling wake demonstrates that this technique can correctly predict the frequency, growth rate and dominant wavemaker region of the flow. It also shows that the zone of absolute instability found from the spatio-temporal analysis is a good approximation to the wavemaker region, which is found by overlapping the direct and adjoint global modes. This approximation is used in the other two flows because it is difficult to calculate their adjoint global modes. Its application to the single stream fuel injector demonstrates that it can identify the regions of the flow that are responsible for generating the hydrodynamic oscillations seen in LES and experimental data. The frequencies predicted by this technique are within a few percent of the measured frequencies. The technique also explains why these oscillations become weaker when a central jet is injected along the centreline. This is because the absolutely unstable region that causes the oscillations becomes convectively unstable. Its application to the lean premixed gas turbine injector reveals that several regions of the flow are hydrodynamically unstable, each with a different frequency and a different strength. For example, it reveals that the central region of confined swirling flow is strongly absolutely unstable and sets up a precessing vortex core, which is likely to aid mixing throughout the injector. It also reveals that the region between the second and third streams is slightly absolutely unstable at a frequency that is likely to coincide with acoustic modes within the combustion chamber. This technique, coupled with knowledge of the acoustic modes in a combustion chamber, is likely to be a useful design tool for the passive control of mixing and combustion instability. Copyright © 2012 by ASME.

Angle-of-attack effects on counter-rotating propellers at take-off

Due to their potential for significant fuel consumption savings, Counter-Rotating Open Rotors (CRORs) are currently being considered as an alternative to high-bypass turbofans. When CRORs are mounted on an aircraft, several 'installation effects' arise which are not present when the engine is operated in isolation. This paper investigates how flow features arising from one such effect - The angle-of-attack of the engine centre-line relative to the oncoming flow - can influence the design of CROR engines. Three-dimensional full-annulus unsteady CFD simulations are used to predict the time-varying flow field experienced by each rotor and emphasis is put on the interaction of the frontrotor wake and tip vortex with the rear-rotor. A parametric study is presented that quantifies the rotorrotor interaction as a function of the angle-of-attack. It is shown that angle-of-attack operation significantly changes the flow field and the unsteady lift on both rotors. In particular, a frequency analysis shows that the unsteady lift exhibits sidebands around the rotor-rotor interaction frequencies. Further, a non-linear increase in the total rear-rotor tip unsteadiness is observed for moderate and high angles-of-attack. The results presented in this paper demonstrate that common techniques used to mitigate CROR noise, such as modifying the rotor-rotor axial spacing and rear-rotor crop, can not be applied correctly unless angle-of-attack effects are taken into account. Copyright © 2012 by ASME.

Vortical structures on three-dimensional shock control bumps

The vortical wake structure produced by a three-dimensional shock control bump (SCB) is thought to be useful for controlling transonic buffet on airfoils. However, at present the vorticity produced is relatively weak and the production mechanism is not well understood. Using a combined experimental and computational approach, a preliminary investigation on the wake vorticity for different bump geometries has been carried out. The structure of the wake for on and off-design conditions are considered, and the effects on the downstream boundary layer demonstrated. Three main vortical structures are observed: a primary vortex pair, weak inter-bump vortices and shear flow in the lambda-shock region. The effect of pressure gradients on vortex strength is examined and it is found that spanwise pressure gradients on the front section of the bump are the most significant parameter influencing vortex strength. © 2013 by S.P. Colliss et al.

Experimental studies of an accelerating, pitching, flat plate at low Reynolds number

The understanding of low Reynolds number aerodynamics is becoming increasingly prevalent with the recent surge in interest in advanced Micro-Air Vehicle (MAV) technology. Research in this area has been primarily stimulated by a military need for smaller, more versatile, autonomous, surveillance aircraft. The mechanism for providing the high lift coefficient required forMAV applications is thought to be largely influenced by the formation of a Leading Edge Vortex (LEV). This paper analyses experimentally, the influence of the LEV effect for a flat plate wing (AR = 4) under fast and slow pitch-up motions at Re =10,000 using a combination of dye flow visualisation and PIV measurements. It is found that a fast pitch over 1c shows a flow topology dominant LEV, while for a slow pitch case over 6c, the flow is largely separated. The development of the suction surface flow and the LEV was strongly correlated with the kinematics of the leading edge, suggesting that the effective local angle of incidence at the Leading Edge (LE) is of considerable significance in unsteady pitching motions. © 2013 by P.R.R.J Stevens.

Detailed measurements of spike formation in an axial compressor

This paper presents new experimental measurements of spike-type stall inception. The measurements were carried out in the single stage Deverson compressor at the Whittle Laboratory. The primary objective was to characterize the flow field in the tip clearance gap during stall inception using sufficient instrumentation to give high spatial and temporal resolution. Measurements were recorded using arrays of unsteady pressure transducers over the rotor tips and hot-wires in the tip gap. Pre-stall ensemble averaged velocity and pressure maps were obtained as well as pressure contours of the stall event. In order to study the transient inception process in greater detail, vector maps were built up from hundreds of stalling events using a triggering system based on the stalling event itself. The results show an embryonic disturbance starting within the blade passage and leading to the formation of a clear spike. The origins of the spike and its relation to the tip leakage vortex are discussed. It has also been shown that before stall the flow in the blade passage which is most likely to stall is generally more unsteady, from revolution to revolution, than the other passages in the annulus. Copyright © 2012 by ASME.

Mixing and internal dynamics of droplets impacting and coalescing on a solid surface.

The coalescence and mixing of a sessile and an impacting liquid droplet on a solid surface are studied experimentally and numerically in terms of lateral separation and droplet speed. Two droplet generators are used to produce differently colored droplets. Two high-speed imaging systems are used to investigate the impact and coalescence of the droplets in color from a side view with a simultaneous gray-scale view from below. Millimeter-sized droplets were used with dynamical conditions, based on the Reynolds and Weber numbers, relevant to microfluidics and commercial inkjet printing. Experimental measurements of advancing and receding static contact angles are used to calibrate a contact angle hysteresis model within a lattice Boltzmann framework, which is shown to capture the observed dynamics qualitatively and the final droplet configuration quantitatively. Our results show that no detectable mixing occurs during impact and coalescence of similar-sized droplets, but when the sessile droplet is sufficiently larger than the impacting droplet vortex ring generation can be observed. Finally we show how a gradient of wettability on the substrate can potentially enhance mixing.

On the two classes of global primary modal instability in laminar separation bubbles

The self-excited global instability mechanisms existing in flat-plate laminar separation bubbles are studied here, in order to shed light on the causes of unsteadiness and three- dimensionality of unforced, nominally two-dimensional separated flows. The presence of two known linear global mechanisms, namely an oscillator behavior driven by local regions of absolute inflectional instability and a centrifugal instability giving rise to a steady three- dimensionalization of the bubble, is studied in a series of model separation bubbles. Present results indicate that absolute instability, and consequently a global oscillator behavior, does not exist for two-dimensional bubbles with a peak reversed-flow velocity below 12% of the free-stream velocity. However, the three-dimensional instability becomes active for recirculation levels as low as urev ≈ 7%. These findings suggest a route to the three-dimensionality and unsteadiness observed in experiments and simulations substantially different from that usually found in the literature, in which two-dimensional vortex shedding is followed by three-dimensionalization.

High-order detached eddy simulation, zonal LES and URANS of cavity and labyrinth seal flows

Hybrid numerical large eddy simulation (NLES), detached eddy simulation (DES) and URANS methods are assessed on a cavity and a labyrinth seal geometry. A high sixth-order discretization scheme is used and is validated using the test case of a two-dimensional vortex. The hybrid approach adopts a new blending function. For the URANS simulations, the flow within the cavity remains steady, and the results show significant variation between models. Surprisingly, low levels of resolved turbulence are observed in the cavity for the DES simulation, and the cavity shear layer remains two dimensional. The hybrid RANS-NLES approach does not suffer from this trait.For the labyrinth seal, both the URANS and DES approaches give low levels of resolved turbulence. The zonal Hamilton-Jacobi approach on the other had given significantly more resolved content. Both DES and hybrid RANS-NLES give good agreement with the experimentally measured velocity profiles. Again, there is significant variation between the URANS models, and swirl velocities are overpredicted. © 2013 John Wiley & Sons, Ltd.