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.

Multi-objective shape optimisation for horizontal-axis wind turbine blades

We are developing a wind turbine blade optimisation package CoBOLDT (COmputa- tional Blade Optimisation and Load De ation Tool) for the optimisation of large horizontal- axis wind turbines. The core consists of the Multi-Objective Tabu Search (MOTS), which controls a spline parameterisation module, a fast geometry generation and a stationary Blade Element Momentum (BEM) code to optimise an initial wind turbine blade design. The objective functions we investigate are the Annual Energy Production (AEP) and the fl apwise blade root bending moment (MY0) for a stationary wind speed of 50 m/s. For this task we use nine parameters which define the blade chord, the blade twist (4 parameters each) and the blade radius. Throughout the optimisation a number of binary constraints are defined to limit the noise emission, to allow for transportation on land and to control the aerodynamic conditions during all phases of turbine operation. The test case shows that MOTS is capable to find enhanced designs very fast and eficiently and will provide a rich and well explored Pareto front for the designer to chose from. The optimised blade de- sign could improve the AEP of the initial blade by 5% with the same flapwise root bending moment or reduce MY0 by 7.5% with the original energy yield. Due to the fast runtime of order 10 seconds per design, a huge number of optimisation iterations is possible without the need for a large computing cluster. This also allows for increased design flexibility through the introduction of more parameters per blade function or parameterisation of the airfoils in future. © 2012 by Nordex Energy GmbH.

Multi-objective shape optimisation for horizontal-axis wind turbine blades

We are developing a wind turbine blade optimisation package CoBOLDT (COmputa- tional Blade Optimisation and Load Deation Tool) for the optimisation of large horizontal- axis wind turbines. The core consists of the Multi-Objective Tabu Search (MOTS), which controls a spline parameterisation module, a fast geometry generation and a stationary Blade Element Momentum (BEM) code to optimise an initial wind turbine blade design. The objective functions we investigate are the Annual Energy Production (AEP) and the apwise blade root bending moment (MY0) for a stationary wind speed of 50 m/s. For this task we use nine parameters which define the blade chord, the blade twist (4 parameters each) and the blade radius. Throughout the optimisation a number of binary constraints are defined to limit the noise emission, to allow for transportation on land and to control the aerodynamic conditions during all phases of turbine operation. The test case shows that MOTS is capable to find enhanced designs very fast and efficiently and will provide a rich and well explored Pareto front for the designer to chose from. The optimised blade de- sign could improve the AEP of the initial blade by 5% with the same apwise root bending moment or reduce MY0 by 7.5% with the original energy yield. Due to the fast runtime of order 10 seconds per design, a huge number of optimisation iterations is possible without the need for a large computing cluster. This also allows for increased design flexibility through the introduction of more parameters per blade function or parameterisation of the airfoils in future. © 2012 AIAA.

Dynamic compression of foam supported plates impacted by high velocity soil

The response of back-supported buffer plates comprising a solid face sheet and foam core backing impacted by a column of high velocity particles (sand slug) is investigated via a lumped parameter model and coupled discrete/continuum simulations. The buffer plate is either resting on (unattached) or attached to a rigid stationary foundation. The lumped parameter model is used to construct maps of the regimes of behaviour with axes of the ratio of the height of the sand slug to core thickness and the normalised core strength. Four regimes of behaviour are identified based on whether the core compression ends prior to the densification of the sand slug or vice versa. Coupled discrete/continuum simulations are also reported and compared with the lumped parameter model. While the model predicted regimes of behaviour are in excellent agreement with numerical simulations, the lumped parameter model is unable to predict the momentum transmitted to the supports as it neglects the role of elasticity in both the buffer plate and the sand slug. The numerical calculations show that the momentum transfer is minimised for intermediate values of the core strength when the so-called "soft-catch" mechanism is in play. In this regime the bounce-back of the sand slug is minimised which reduces the momentum transfer. However, in this regime, the impulse reduction is small (less than 10% of that transferred to a rigid structure). For high values of the core strength, the response of the buffer plate resembles a rigid plate with nearly no impulse mitigation while at low values of core strength, a slap event occurs when the face sheet impinges against the foundation due to full densification of the foam core. This slap event results in a significant enhancement of the momentum transfer to the foundation. The results demonstrate that appropriately designed buffer plates have potential as impulse mitigators in landmine loading situations. © 2013 Elsevier Ltd. All rights reserved.

Multi-objective optimisation of horizontal axis wind turbine structure and energy production using aerofoil and blade properties as design variables

The design of wind turbine blades is a true multi-objective engineering task. The aerodynamic effectiveness of the turbine needs to be balanced with the system loads introduced by the rotor. Moreover the problem is not dependent on a single geometric property, but besides other parameters on a combination of aerofoil family and various blade functions. The aim of this paper is therefore to present a tool which can help designers to get a deeper insight into the complexity of the design space and to find a blade design which is likely to have a low cost of energy. For the research we use a Computational Blade Optimisation and Load Deflation Tool (CoBOLDT) to investigate the three extreme point designs obtained from a multi-objective optimisation of turbine thrust, annual energy production as well as mass for a horizontal axis wind turbine blade. The optimisation algorithm utilised is based on Multi-Objective Tabu Search which constitutes the core of CoBOLDT. The methodology is capable to parametrise the spanning aerofoils with two-dimensional Free Form Deformation and blade functions with two tangentially connected cubic splines. After geometry generation we use a panel code to create aerofoil polars and a stationary Blade Element Momentum code to evaluate turbine performance. Finally, the obtained loads are fed into a structural layout module to estimate the mass and stiffness of the current blade by means of a fully stressed design. For the presented test case we chose post optimisation analysis with parallel coordinates to reveal geometrical features of the extreme point designs and to select a compromise design from the Pareto set. The research revealed that a blade with a feasible laminate layout can be obtained, that can increase the energy capture and lower steady state systems loads. The reduced aerofoil camber and an increased L/. D-ratio could be identified as the main drivers. This statement could not be made with other tools of the research community before. © 2013 Elsevier Ltd.

Nonlinear thermoacoustics of ducted premixed flames: The influence of perturbation convection speed

When a premixed flame is placed within a duct, acoustic waves induce velocity perturbations at the flame's base. These travel down the flame, distorting its surface and modulating its heat release. This can induce self-sustained thermoacoustic oscillations. Although the phase speed of these perturbations is often assumed to equal the mean flow speed, experiments conducted in other studies and Direct Numerical Simulation (DNS) conducted in this study show that it varies with the acoustic frequency. In this paper, we examine how these variations affect the nonlinear thermoacoustic behaviour. We model the heat release with a nonlinear kinematic G-equation, in which the velocity perturbation is modelled on DNS results. The acoustics are governed by linearised momentum and energy equations. We calculate the flame describing function (FDF) using harmonic forcing at several frequencies and amplitudes. Then we calculate thermoacoustic limit cycles and explain their existence and stability by examining the amplitude-dependence of the gain and phase of the FDF. We find that, when the phase speed equals the mean flow speed, the system has only one stable state. When the phase speed does not equal the mean flow speed, however, the system supports multiple limit cycles because the phase of the FDF changes significantly with oscillation amplitude. This shows that the phase speed of velocity perturbations has a strong influence on the nonlinear thermoacoustic behaviour of ducted premixed flames. © 2013 The Combustion Institute.

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 (a), compound angle (b), 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 a and b but with the same angle between the mainstream and coolant flow directions (angle k). This relationship is assessed through experiments by testing two sets of cylindrical holes with different a and b: one set with k=35 deg, and another set with k=60 deg. The data confirm the stated relationship between α, β, k and the aerodynamic mixing loss. The results show that the designer should minimize k to obtain the lowest loss, but maximize b to achieve the best heat transfer performance. A suggestion on improving the loss model is also given. Five different hole geometries (α=35.0 deg, β=0 deg) 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. © 2013 by ASME.

Optical trapping and manipulation of nanostructures

Optical trapping and manipulation of micrometre-sized particles was first reported in 1970. Since then, it has been successfully implemented in two size ranges: the subnanometre scale, where light-matter mechanical coupling enables cooling of atoms, ions and molecules, and the micrometre scale, where the momentum transfer resulting from light scattering allows manipulation of microscopic objects such as cells. But it has been difficult to apply these techniques to the intermediate-nanoscale-range that includes structures such as quantum dots, nanowires, nanotubes, graphene and two-dimensional crystals, all of crucial importance for nanomaterials-based applications. Recently, however, several new approaches have been developed and demonstrated for trapping plasmonic nanoparticles, semiconductor nanowires and carbon nanostructures. Here we review the state-of-the-art in optical trapping at the nanoscale, with an emphasis on some of the most promising advances, such as controlled manipulation and assembly of individual and multiple nanostructures, force measurement with femtonewton resolution, and biosensors.

The evolution of a stratified turbulent cloud

Localized regions of turbulence, or turbulent clouds, in a stratified fluid are the subject of this study, which focuses on the edge dynamics occurring between the turbulence and the surrounding quiescent region. Through laboratory experiments and numerical simulations of stratified turbulent clouds, we confirm that the edge dynamics can be subdivided into materially driven intrusions and horizontally travelling internal wave-packets. Three-dimensional visualizations show that the internal gravity wave-packets are in fact large-scale pancake structures that grow out of the turbulent cloud into the adjacent quiescent region. The wave-packets were tracked in time, and it is found that their speed obeys the group speed relation for linear internal gravity waves. The energetics of the propagating waves, which include waveforms that are inclined with respect to the horizontal, are also considered and it is found that, after a period of two eddy turnover times, the internal gravity waves carry up to 16 % of the cloud kinetic energy into the initially quiescent region. Turbulent events in nature are often in the form of decaying turbulent clouds, and it is therefore suggested that internal gravity waves radiated from an initial cloud could play a significant role in the reorganization of energy and momentum in the atmosphere and oceans.©2013 Cambridge University Press.

Water droplet flow paths and droplet deposition in low pressure steam turbines

The complex three-dimensional two-phase flow in a low pressure steam turbine is investigated with comprehensive numerical flow simulations. In addition to the condensation process, which already takes place in the last stages of steam turbines, the numerical flow model is enhanced to consider the drag forces between the droplets and the vapour phase. The present paper shows the differences in the flow path of the phases and investigates the effect of an increasing droplet diameter. For the flow simulations a performance cluster is used because of the high effort for such multi-momentum two-phase flow calculations. In steam turbines the deposition of small water droplets on the stator blades or on parts of the casing is responsible for the formation of large coarse water droplets and these may cause additional dissipation as well as damage due to blade erosion. A method is presented that uses detailed CFD data to predict droplet deposition on turbine stator blades. This simulation method to detect regions of droplet deposition can help to improve the design of water removal devices. © Springer-Verlag Berlin Heidelberg 2013.

A flame surface density approach to large-eddy simulation of premixed turbulent combustion

The flame surface density approach to the modeling of premixed turbulent combustion is well established in the context of Reynolds-averaged simulations. For the future, it is necessary to consider large-eddy simulation (LES), which is likely to offer major advantages in terms of physical accuracy, particularly for unsteady combustion problems. LES relies on spatial filtering for the removal of unresolved phenomena whose characteristic length scales are smaller than the computational grid scale. Thus, there is a need for soundly based physical modeling at the subgrid scales. The aim of this paper is to explore the usefulness of the flame surface density concept as a basis for LES modeling of premixed turbulent combustion. A transport equation for the filtered flame surface density is presented, and models are proposed for unclosed terms. Comparison with Reynolds-averaged modeling is shown to reveal some interesting similarities and differences. These were exploited together with known physics and statistical results from experiment and from direct numerical stimulation in order to gain insight and refine the modeling. The model has been implemented in a combustion LES code together with standard models for scalar and momentum transport. Computational results were obtained for a simple three-dimensional flame propagation test problem, and the relative importance of contributing terms in the modeled equation for flame surface density was assessed. Straining and curvature are shown to have a major influence at both the resolved and subgrid levels.

Structure of cryogenic flames at elevated pressures

This paper presents new experimental results on cryogenic jet flames formed by a coaxial injector at a pressure of 70 bar, which approaches the pressures found in rocket engines. This element, fed with liquid oxygen and gaseous hydrogen, is placed in a square combustion chamber equipped with quartz windows. The flame is examined via spectroscopy, OH* emission, and backlighting, the aim being to provide basic information on the flame structure. It is found that some of the OH* emission is absorbed by the OH radicals present in the flame. A detailed examination of this effect is presented, in which it is shown that, for this turbulent flame, the Abel transform gives the position of the intense reaction region, whether or not absorption is signficant. The flame is attached to the oxygen injector, as at low pressure. At high pressure, flame expansion is reduced compared with low pressure and is also less dependent on the momentum flux ratio between the hydrogen and the oxygen streams. An analysis of the relevant Damköhler numbers suggests that this is because the rate of combustion is mainly controlled by large-scale turbulent mixing at high pressure, and it is dominated by jet break-up, atomization, and vaporization at low pressures. Jet break-up is particularly dependent on the momentum flux ratio. Finally, the mean volumetric heat release rates and flame surface density in the experimental facility are estimated.