Influence of film cooling hole angles and geometries on aerodynamic loss and net heat flux reduction

Turbine design engineers have to ensure that film cooling can provide sufficient protection to turbine blades from the hot mainstream gas, while keeping the losses low. Film cooling hole design parameters include inclination angle (α), compound angle (β ), hole inlet geometry and hole exit geometry. The influence of these parameters on aerodynamic loss and net heat flux reduction is investigated, with loss being the primary focus. Low-speed flat plate experiments have been conducted at momentum flux ratios of IR = 0.16, 0.64 and 1.44. The film cooling aerodynamic mixing loss, generated by the mixing of mainstream and coolant, can be quantified using a three-dimensional analytical model that has been previously reported by the authors. The model suggests that for the same flow conditions, the aerodynamic mixing loss is the same for holes with different α and β but with the same angle between the mainstream and coolant flow directions (angle κ). This relationship is assessed through experiments by testing two sets of cylindrical holes with different α and β : one set with κ = 35°, another set with κ = 60°. The data confirm the stated relationship between α, β, κ and the aerodynamic mixing loss. The results show that the designer should minimise κ to obtain the lowest loss, but maximise β to achieve the best heat transfer performance. A suggestion on improving the loss model is also given. Five different hole geometries (α =35.0°, β =0°) were also tested: cylindrical hole, trenched hole, fan-shaped hole, D-Fan and SD-Fan. The D-Fan and the SD-Fan have similar hole exits to the fan-shaped hole but their hole inlets are laterally expanded. The external mixing loss and the loss generated inside the hole are compared. It was found that the D-Fan and the SD-Fan have the lowest loss. This is attributed to their laterally expanded hole inlets, which lead to significant reduction in the loss generated inside the holes. As a result, the loss of these geometries is ≈ 50 % of the loss of the fan-shaped hole at IR = 0.64 and 1.44. Copyright © 2011 by ASME.

The influence of turbine stator well coolant flow rate and passage configuration on cooling effectiveness

Market competitiveness for aero engine power plant dictates that improvements in engine performance and reliability are guaranteed a priori by manufacturers. The requirement to accurately predict the life of engine components makes exacting demands of the internal air system, which must provide effective cooling over the engine duty cycle with the minimum consumption of compressor section air. Tests have been conducted at the University of Sussex using a turbine test facility which comprises a two stage turbine with an individual stage pressure ratio of 1.7:1. Main annulus air is supplied by an adapted Rolls-Royce Dart compressor at up to 440 K and 4.8 kg s-1. Cooling flow rates ranging from 0.71 to 1.46 Cw, ent, a disc entrainment parameter, have been used to allow ingress or egress dominated stator well flow conditions. The mechanical design of the test section allows internal cooling geometry to be rapidly re-configured, allowing the effect of jet momentum and coolant trajectory to be investigated. An important facet to this investigation is the use of CFD to model and analyse the flow structures associated with the cavity conditions tested, as well as to inform the design of cooling path geometry. This paper reports on the effectiveness of stator well coolant flow rate and delivery configurations using experimental data and also CFD analysis to better quantify the effect of stator well flow distribution on component temperatures. Copyright © 2011 by Rolls-Royce plc.

Gamma-ray spectroscopy at TRIUMF-ISAC

The 8π spectrometer at TRIUMF-ISAC consists of 20 Compton-suppressed germanium detectors and various auxiliary devices. The Ge array, once used for studies of nuclei at high angular momentum, has been transformed into the world's most powerful device dedicated to radioactive-decay studies. Many improvements in the spectrometer have been made, including a high-throughput data acquisition system, installation of a moving tape collector, incorporation of an array of 20 plastic scintillators for β-particle tagging, 5 Si(Li) detectors for conversion electrons, and 10 BaF2 detectors for fast-lifetime measurements. Experiments can be performed where data from all detectors are collected simultaneously, resulting in a very detailed view of the nucleus through radioactive decay. A number of experimental programmes have been launched that take advantage of the versatility of the spectrometer, and the intense beams available at TRIUMF-ISAC. © 2006 American Institute of Physics.

A review on slip models for gas microflows

Accurate modeling of gas microflow is crucial for the microfluidic devices in MEMS. Gas microflows through these devices are often in the slip and transition flow regimes, characterized by the Knudsen number of the order of 10-2∼100. An increasing number of researchers now dedicate great attention to the developments in the modeling of non-equilibrium boundary conditions in the gas microflows, concentrating on the slip model. In this review, we present various slip models obtained from different theoretical, computational and experimental studies for gas microflows. Correct descriptions of the Knudsen layer effect are of critical importance in modeling and designing of gas microflow systems and in predicting their performances. Theoretical descriptions of the gas-surface interaction and gas-surface molecular interaction models are introduced to describe the boundary conditions. Various methods and techniques for determination of the slip coefficients are reviewed. The review presents the considerable success in the implementation of various slip boundary conditions to extend the Navier-Stokes (N-S) equations into the slip and transition flow regimes. Comparisons of different values and formulations of the first- and second-order slip coefficients and models reveal the discrepancies arising from different definitions in the first-order slip coefficient and various approaches to determine the second-order slip coefficient. In addition, no consensus has been reached on the correct and generalized form of higher-order slip expression. The influences of specific effects, such as effective mean free path of the gas molecules and viscosity, surface roughness, gas composition and tangential momentum accommodation coefficient, on the hybrid slip models for gas microflows are analyzed and discussed. It shows that although the various hybrid slip models are proposed from different viewpoints, they can contribute to N-S equations for capturing the high Knudsen number effects in the slip and transition flow regimes. Future studies are also discussed for improving the understanding of gas microflows and enabling us to exactly predict and actively control gas slip. © Springer-Verlag 2012.

Thermal stratification produced by plumes and jets in enclosed spaces

The airflow and thermal stratification produced by a localised heat source located at floor level in a closed room is of considerable practical interest and is commonly referred to as a 'filling box'. In rooms with low aspect ratios H/R ≲ 1 (room height H to characteristic horizontal dimension R) the thermal plume spreads laterally on reaching the ceiling and a descending horizontal 'front' forms separating a stably stratified, warm upper region from cooler air below. The stratification is well predicted for H/R ≲ 1 by the original filling box model of Baines and Turner (J. Fluid. Mech. 37 (1968) 51). This model represents a somewhat idealised situation of a plume rising from a point source of buoyancy alone-in particular the momentum flux at the source is zero. In practical situations, real sources of heating and cooling in a ventilation system often include initial fluxes of both buoyancy and momentum, e.g. where a heating system vents warm air into a space. This paper describes laboratory experiments to determine the dependence of the 'front' formation and stratification on the source momentum and buoyancy fluxes of a single source, and on the location and relative strengths of two sources from which momentum and buoyancy fluxes were supplied separately. For a single source with a non-zero input of momentum, the rate of descent of the front is more rapid than for the case of zero source momentum flux and increases with increasing momentum input. Increasing the source momentum flux effectively increases the height of the enclosure, and leads to enhanced overturning motions and finally to complete mixing for highly momentum-driven flows. Stratified flows may be maintained by reducing the aspect ratio of the enclosure. At these low aspect ratios different long-time behaviour is observed depending on the nature of the heat input. A constant heat flux always produces a stratified interior at large times. On the other hand, a constant temperature supply ultimately produces a well-mixed space at the supply temperature. For separate sources of momentum and buoyancy, the developing stratification is shown to be strongly dependent on the separation of the sources and their relative strengths. Even at small separation distances the stratification initially exhibits horizontal inhomogeneity with localised regions of warm fluid (from the buoyancy source) and cool fluid. This inhomogeneity is less pronounced as the strength of one source is increased relative to the other. Regardless of the strengths of the sources, a constant buoyancy flux source dominates after sufficiently large times, although the strength of the momentum source determines whether the enclosure is initially well mixed (strong momentum source) or stably stratified (weak momentum source). © 2001 Elsevier Science Ltd. All rights reserved.

Sensitivity analysis of a time-delayed thermo-acoustic system via an adjoint-based approach

We apply adjoint-based sensitivity analysis to a time-delayed thermo-acoustic system: a Rijke tube containing a hot wire. We calculate how the growth rate and frequency of small oscillations about a base state are affected either by a generic passive control element in the system (the structural sensitivity analysis) or by a generic change to its base state (the base-state sensitivity analysis). We illustrate the structural sensitivity by calculating the effect of a second hot wire with a small heat-release parameter. In a single calculation, this shows how the second hot wire changes the growth rate and frequency of the small oscillations, as a function of its position in the tube. We then examine the components of the structural sensitivity in order to determine the passive control mechanism that has the strongest influence on the growth rate. We find that a force applied to the acoustic momentum equation in the opposite direction to the instantaneous velocity is the most stabilizing feedback mechanism. We also find that its effect is maximized when it is placed at the downstream end of the tube. This feedback mechanism could be supplied, for example, by an adiabatic mesh. We illustrate the base-state sensitivity by calculating the effects of small variations in the damping factor, the heat-release time-delay coefficient, the heat-release parameter, and the hot-wire location. The successful application of sensitivity analysis to thermo-acoustics opens up new possibilities for the passive control of thermo-acoustic oscillations by providing gradient information that can be combined with constrained optimization algorithms in order to reduce linear growth rates. © Cambridge University Press 2013.

Structural sensitivity of spiral vortex breakdown

Previous numerical simulations have shown that vortex breakdown starts with the formation of a steady axisymmetric bubble and that an unsteady spiralling mode then develops on top of this. We investigate this spiral mode with a linear global stability analysis around the steady bubble and its wake. We obtain the linear direct and adjoint global modes of the linearized Navier-Stokes equations and overlap these to obtain the structural sensitivity of the spiral mode, which identifies the wavemaker region. We also identify regions of absolute instability with a local stability analysis. At moderate swirls, we find that the m=-1 azimuthal mode is the most unstable and that the wavemaker regions of the m=-1 mode lie around the bubble, which is absolutely unstable. The mode is most sensitive to feedback involving the radial and azimuthal components of momentum in the region just upstream of the bubble. To a lesser extent, the mode is also sensitive to feedback involving the axial component of momentum in regions of high shear around the bubble. At an intermediate swirl, in which the bubble and wake have similar absolute growth rates, other researchers have found that the wavemaker of the nonlinear global mode lies in the wake. We agree with their analysis but find that the regions around the bubble are more influential than the wake in determining the growth rate and frequency of the linear global mode. The results from this paper provide the first steps towards passive control strategies for spiral vortex breakdown. © 2013 Cambridge University Press.

An Arbitrary Lagrangian Eulerian Method for Three-Phase Flows with Triple Junction Points

We present a moving mesh method suitable for solving two-dimensional and axisymmetric three-liquid flows with triple junction points. This method employs a body-fitted unstructured mesh where the interfaces between liquids are lines of the mesh system, and the triple junction points (if exist) are mesh nodes. To enhance the accuracy and the efficiency of the method, the mesh is constantly adapted to the evolution of the interfaces by refining and coarsening the mesh locally; dynamic boundary conditions on interfaces, in particular the triple points, are therefore incorporated naturally and accurately in a Finite- Element formulation. In order to allow pressure discontinuity across interfaces, double-values of pressure are necessary for interface nodes and triple-values of pressure on triple junction points. The resulting non-linear system of mass and momentum conservation is then solved by an Uzawa method, with the zero resultant condition on triple points reinforced at each time step. The method is used to investigate the rising of a liquid drop with an attached bubble in a lighter liquid.

Local equilibrium of the Gibbs interface in two-phase systems

We analyze the local equilibrium assumption for interfaces from the perspective of gauge transformations, which are the small displacements of Gibbs' dividing surface. The gauge invariance of thermodynamic properties turns out to be equivalent to conditions for jumps of bulk densities across the interface. This insight strengthens the foundations of the local equilibrium assumption for interfaces and can be used to characterize nonequilibrium interfaces in a compact and consistent way, with a clear focus on gauge-invariant properties. Using the principle of gauge invariance, we show that the validity of Clapeyron equations can be extended to nonequilibrium interfaces, and an additional jump condition for the momentum density is recognized to be of the Clapeyron type. © 2012 Europhysics Letters Association.

The impact of sand slugs against beams and plates: Coupled discrete particle/finite element simulations

The impact of a slug of dry sand particles against a metallic sandwich beam or circular sandwich plate is analysed in order to aid the design of sandwich panels for shock mitigation. The sand particles interact via a combined linear-spring-and-dashpot law whereas the face sheets and compressible core of the sandwich beam and plate are treated as rate-sensitive, elastic-plastic solids. The majority of the calculations are performed in two dimensions and entail the transverse impact of end-clamped monolithic and sandwich beams, with plane strain conditions imposed. The sand slug is of rectangular shape and comprises a random loose packing of identical, circular cylindrical particles. These calculations reveal that loading due to the sand is primarily inertial in nature with negligible fluid-structure interaction: the momentum transmitted to the beam is approximately equal to that of the incoming sand slug. For a slug of given incoming momentum, the dynamic deflection of the beam increases with decreasing duration of sand-loading until the impulsive limit is attained. Sandwich beams with thick, strong cores significantly outperform monolithic beams of equal areal mass. This performance enhancement is traced to the "sandwich effect" whereby the sandwich beams have a higher bending strength than that of the monolithic beams. Three-dimensional (3D) calculations are also performed such that the sand slug has the shape of a circular cylindrical column of finite height, and contains spherical sand particles. The 3D slug impacts a circular monolithic plate or sandwich plate and we show that sandwich plates with thick strong cores again outperform monolithic plates of equal areal mass. Finally, we demonstrate that impact by sand particles is equivalent to impact by a crushable foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally, yet give predictions to within 5% of the full discrete particle calculations for the monolithic and sandwich beams and plates. These foam projectile calculations suggest that metallic foam projectiles can be used to simulate the loading by sand particles within a laboratory setting. © 2013 Elsevier Ltd.

Nonlinear phenomena in thermoacoustic systems with premixed flames

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies, amplitudes and stability of limit cycles. Limit cycles in thermoacoustic systems are reached when the energy input from driving processes and energy losses from damping processes balance each other over a cycle of the oscillation. In this paper an integral relation for the rate of change of energy of a thermoacoustic system is derived. This relation is analogous to the well-known Rayleigh criterion in thermoacoustics, but can be used to calculate the amplitudes of limit cycles, as well as their stability. The relation is applied to a thermoacoustic system of a ducted slot-stabilized 2-D premixed flame. The flame is modelled using a nonlinear kinematic model based on the G-equation, while the acoustics of planar waves in the tube are governed by linearised momentum and energy equations. Using open-loop forced simulations, the flame describing function (FDF) is calculated. The gain and phase information from the FDF is used with the integral relation to construct a cyclic integral rate of change of energy (CIRCE) diagram that indicates the amplitude and stability of limit cycles. This diagram is also used to identify the types of bifurcation the system exhibits and to find the minimum amplitude of excitation needed to reach a stable limit cycle from another linearly stable state, for single- mode thermoacoustic systems. Furthermore, this diagram shows precisely how the choice of velocity model and the amplitudedependence of the gain and the phase of the FDF influence the nonlinear dynamics of the system. Time domain simulations of the coupled thermoacoustic system are performed with a Galerkin discretization for acoustic pressure and velocity. Limit cycle calculations using a single mode, as well as twenty modes, are compared against predictions from the CIRCE diagram. For the single mode system, the time domain calculations agree well with the frequency domain predictions. The heat release rate is highly nonlinear but, because there is only a single acoustic mode, this does not affect the limit cycle amplitude. For the twenty-mode system, however, the higher harmonics of the heat release rate and acoustic velocity interact resulting in a larger limit cycle amplitude. Multimode simulations show that in some situations the contribution from higher harmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Copyright © 2012 by ASME.

The dynamic response of edge clamped plates loaded by spherically expanding sand shells

The dynamic deformation of both edge clamped stainless steel sandwich panels with a pyramidal truss core and equal mass monolithic plates loaded by spherically expanding shells of dry and water saturated sand has been investigated, both experimentally and via a particle based simulation methodology. The spherically expanding sand shell is generated by detonating a sphere of explosive surrounded by a shell of either dry or water saturated synthetic sand. The measurements show that the sandwich panel and plate deflections decrease with increasing stand-off between the center of the charge and the front of the test structures. Moreover, for the same charge and sand mass, the deflections of the plates are significantly higher in the water saturated sand case compared to that of dry sand. For a given stand-off, the mid-span deflection of the sandwich panel rear faces was substantially less than that of the corresponding monolithic plate for both the dry and water saturated sand cases. The experiments were simulated via a coupled discrete-particle/ finite element scheme wherein the high velocity impacting sand is modeled by interacting particles while the plate is modeled within a Lagrangian finite element setting. The simulations are in good agreement with the measurements for the dry sand impact of both the monolithic and sandwich structures. However, the simulations underestimate the effect of stand-off in the case of the water saturated sand explosion, i.e. the deflections decrease more sharply with increasing stand-off in the experiments compared to the simulations. The simulations reveal that the momentum transmitted into the sandwich and monolithic plate structures by the sand shell is approximately the same, consistent with a small fluid-structure interaction effect. The smaller deflection of the sandwich panels is therefore primarily due to the higher bending strength of sandwich structures. © 2013 The Authors. Published by Elsevier Ltd. All rights reserved.

The breakdown of the viscous regime in riblets

We investigate the mechanisms involved in the breakdown of the viscous regime in riblets, with a view to determining the point of optimum performance, where drag reduction ceases to be proportional to the riblet size. This occurs empirically for a groove cross-section $A_g^+ \approx 120^+$. To study the interaction of the riblets with the overlaying turbulent flow, we systematically conduct DNSes in a ribbed turbulent channel with increasing riblet size. The conditionally averaged crossflow above and within the grooves reveals a mean recirculation bubble that exists up to the point of viscous breakdown, isolating the groove floor from the overlying crossflow, and preventing the high momentum fluid from entering the grooves. We do not find evidence of outside vortices lodging within the grooves until $A_g^+ \approx 400$, which is well past the drag minimum, and already into the drag increasing regime. Interestingly, as the bubble breaks down, we observe that quasi-two-dimensional spanwise structures form just above the riblets, similar to those observed above porous surfaces and plant canopies, which appear to be involved in the performance degradation.

Sand column impact onto a Kolsky pressure bar

A laboratory-based methodology to launch cylindrical sand slugs at high velocities is developed. The methodology generates well-characterised soil ejecta without the need for detonation of an explosive; this laboratory-based tool thereby allows for the experimental investigation of the soil-structure events. The experimental set-up comprises a launcher with a cylindrical cavity and a piston to push out the sand slug. The apparatus is used to launch both dry and water-saturated sand slugs. High speed photography is used to characterise the evolution of the sand slugs after launch. We find that the diameter of the slugs remains unchanged, and the sand particles possess only an axial component of velocity. However, the sand particles have a uniform spatial gradient of axial velocity and this results in lengthening of the slugs as they travel towards their target. Thus, the density of the sand slugs remains spatially homogenous but decreases with increasing time. The velocity gradient is typically higher in the dry sand slugs than that of the water-saturated slugs. The pressure exerted by the slugs on a rigid-stationary target is measured by impacting the slugs against a direct impact Kolsky bar. After an initial high transient pressure, the pressure reduces to a value of approximately ρv 2 where ρ is the density of the impacting sand slug and v is the particle velocity. This indicates that loading due to the sand is primarily inertial in nature. The momentum transmitted to the Kolsky bar was approximately equal to the incident momentum of the sand slugs, regardless of whether they are dry or water-saturated. © 2013 Elsevier Ltd. All rights reserved.

A laboratory-scale buried charge simulator

An experimental technique has been developed in order to mimic the effect of landmine loading on materials and structures to be studied in a laboratory setting, without the need for explosives. Compressed gas is discharged beneath a sand layer, simulating the dynamic flow generated by a buried explosive. High speed photography reveals that the stages of soil motion observed during a landmine blast are replicated. The effect of soil saturation and the depth of the sand layer on sand motion are evaluated. Two series of experiments have been performed with the buried charge simulator to characterise subsequent impact of the sand. First, the time variation in pressure and impulse during sand impact on a stationary target is evaluated using a Kolsky bar apparatus. It is found that the pressure pulse imparted to the Kolsky bar consists of two phases: an initial transient phase of high pressure (attributed to wave propagation effects in the impacting sand), followed by a lower pressure phase of longer duration (due to lateral flow of the sand against the Kolsky bar). Both phases make a significant contribution to the total imparted impulse. It is found that wet sand exerts higher peak pressures and imparts a larger total impulse than dry sand. The level of imparted impulse is determined as a function of sand depth, and of stand-off distance between the sand and the impacted end of the Kolsky bar. The second study uses a vertical impulse pendulum to measure the momentum imparted by sand impact to a target which is free to move vertically. The effect of target mass upon imparted momentum is investigated. It is concluded that the laboratory-scale sand impact apparatus is a flexible tool for investigating the interactions between structures and dynamic sand flows. © 2013 Elsevier Ltd. All rights reserved.