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.

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.

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.

Effects of the biodiesel fuel physical properties on the swirl stabilised spray combustion characteristics

An increasin g interest in biofuel applications in modern engines requires a better understanding of biodiesel combustion behaviour. Many numerical studies have been carried out on unsteady combustion of biodiesel in situations similar to diesel engines, but very few studies have been done on the steady combustion of biodiesel in situations similar to a gas turbine combustor environment. The study of biodiesel spray combustion in gas turbine applications is of special interest due to the possible use of biodiesel in the power generation and aviation industries. In modelling spray combustion, an accurate representation of the physical properties of the fuel is a first important step, since spray formation is largely influenced by fuel properties such as viscosity, density, surface tension and vapour pressure. In the present work, a calculated biodiesel properties database based on the measured composition of Fatty Acid Methyl Esters (FAME) has been implemented in a multi-dimensional Computational Fluid Dynamics (CFD) spray simulation code. Simulations of non-reacting and reacting atmospheric-pressure sprays of both diesel and biodiesel have been carried out using a spray burner configuration for which experimental data is available. A pre-defined droplet size probability density function (pdf) has been implemented together with droplet dynamics based on phase Doppler anemometry (PDA) measurements in the near-nozzle region. The gas phase boundary condition for the reacting spray cases is similar to that of the experiment which employs a plain air-blast atomiser and a straight-vane axial swirler for flame stabilisation. A reaction mechanism for heptane has been used to represent the chemistry for both diesel and biodiesel. Simulated flame heights, spray characteristics and gas phase velocities have been found to compare well with the experimental results. In the reacting spray cases, biodiesel shows a smaller mean droplet size compared to that of diesel at a constant fuel mass flow rate. A lack of sensitivity towards different fuel properties has been observed based on the non-reacting spray simulations, which indicates a need for improved models of secondary breakup. By comparing the results of the non-reacting and reacting spray simulations, an improvement in the complexity of the physical modelling is achieved which is necessary in the understanding of the complex physical processes involved in spray combustion simulation. Copyright © 2012 SAE International.

A numerical study into the effects of Bio-synthetic Paraffinic Kerosine blends with Jet-A fuel for civil aircraft engine

Growing concerns regarding fluctuating fuel costs and pollution targets for gas emissions, have led the aviation industry to seek alternative technologies to reduce its dependency on crude oil, and its net emissions. Recently blends of bio-fuel with kerosine, have become an alternative solution as they offer "greener" aircraft and reduce demand on crude oil. Interestingly, this technique is able to be implemented in current aircraft as it does not require any modification to the engine. Therefore, the present study investigates the effect of blends of bio-synthetic paraffinic kerosine with Jet-A in a civil aircraft engine, focusing on its performance and exhaust emissions. Two bio-fuels are considered: Jatropha Bio-synthetic Paraffinic Kerosine (JSPK) and Camelina Bio-synthetic Paraffinic Kerosine (CSPK); there are evaluated as pure fuels, and as 10% and 50% blend with Jet-A. Results obtained show improvement in thrust, fuel flow and SFC as composition of bio-fuel in the blend increases. At design point condition, results on engine emissions show reduction in NO x, and CO, but increases of CO is observed at fixed fuel condition, as the composition of bio-fuel in the mixture increases. Copyright © 2012 by ASME.

Spray flame structure of rapeseed biodiesel and Jet-A1 fuel

The spray combustion characteristics of rapeseed methyl esters (RME) were compared to Jet-A1 fuel using a gas turbine type combustor. The swirling spray flames for both fuels were established at a constant power output of 6 kW. The main swirling air flow was preheated to 350 C prior to coaxially enveloping the airblast-atomized liquid fuel spray at atmospheric pressure. Investigation of the fundamental spray combustion was performed via measurements of the fuel droplet sizes and velocities, gas phase flow fields and flame reaction zones. The spray flame droplets and flow fields in the combustors were characterised using phase Doppler anemometry (PDA) and particle imaging velocimetry (PIV) respectively. Flame chemiluminescence imaging was employed to identify the flame reaction zones. The highest droplet concentration zone extends along a 30 angle from the symmetry axis, inside the flame zone. Only small droplets(<17 μ) (<17 μm)are found around the centreline region, while larger droplets are found at the edge of the spray outside the flame reaction zone. RME exhibits spray characteristics similar to Jet-A1 but with droplet concentration and volume fluxes four times higher, consistent with the expected longer droplet evaporation timescale. The flow field characteristics for both RME and Jet-A1 spray flames are very similar despite the significantly different visible characteristics of the flame reaction zones. © 2013 Elsevier Ltd. All rights reserved.

The closed thorium-transuranic fuel cycle in reduced-moderation PWRs and BWRs

Multiple recycle of long-lived actinides has the potential to greatly reduce the required storage time for spent nuclear fuel or high level nuclear waste. This is generally thought to require fast reactors as most transuranic (TRU) isotopes have low fission probabilities in thermal reactors. Reduced-moderation LWRs are a potential alternative to fast reactors with reduced time to deployment as they are based on commercially mature LWR technology. Thorium (Th) fuel is neutronically advantageous for TRU multiple recycle in LWRs due to a large improvement in the void coefficient. If Th fuel is used in reduced-moderation LWRs, it appears neutronically feasible to achieve full actinide recycle while burning an external supply of TRU, with related potential improvements in waste management and fuel utilization. In this paper, the fuel cycle of TRU-bearing Th fuel is analysed for reduced-moderation PWRs and BWRs (RMPWRs and RBWRs). RMPWRs have the advantage of relatively rapid implementation and intrinsically low conversion ratios, which is desirable to maximize the TRU burning rate. However, it is challenging to simultaneously satisfy operational and fuel cycle constraints. An RBWR may potentially take longer to implement than an RMPWR due to more extensive changes from current BWR technology. However, the harder neutron spectrum can lead to favourable fuel cycle performance. A two-stage TRU burning cycle, where the first stage is Th-Pu MOX in a conventional PWR feeding a second stage continuous burn in RMPWR or RBWR, is technically reasonable, although it is more suitable for the RBWR implementation. In this case, the fuel cycle performance is relatively insensitive to the discharge burn-up of the first stage. © 2013 Elsevier Ltd. All rights reserved.

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-hit ballistic damage characterisation of 304 stainless steel plates with finite elements

Multi-impact of projectiles on thin 304 stainless steel plates is investigated to assess the degradation of ballistic performance, and to characterise the inherent mechanisms. Assessment of ballistic degradation is by means of a double-impact of rigid spheres at the same site on a circular clamped plate. The limiting velocity of the second impact, will be altered by the velocity of the antecedent impact. Finite element analyses were used to elucidate experimental results and understand the underlying mechanisms that give rise to the performance degradation. The effect of strength and ductility on the single and multi-impact performance was also considered. The model captured the experimental results with excellent agreement. Moreover, the material parameters used within the model were exclusively obtained from published works with no fitting or calibration required. An attempt is made to quantify the elevation of the ballistic limit of thin plates by the dynamic mechanism of travelling hinges. Key conclusions: The multi-hit performance scales linearly with the single-hit performance; and strength is a significantly greater effector of increased ballistic limit than ductility, even at the expense of toughness. © 2014 Elsevier Ltd.

Interaction of a premixed flame with a liquid fuel film on a wall

© 2004 The Combustion Institute. Published by Elsevier Inc. All rights reserved. In piston engines and in gas turbines, the injection of liquid fuel often leads to the formation of a liquid film on the combustor wall. If a flame reaches this zone, undesired phenomena such as coking may occur and diminish the lifetime of the engine. Moreover, the effect of such an interaction on maximum wall heat fluxes, flame quenching, and pollutant formation is largely unknown. This paper presents a numerical study of the interaction of a premixed flame with a cold wall covered with a film of liquid fuel. Simulations show that the presence of the film leads to a very rich zone at the wall in which the flame cannot propagate. As a result, the flame wall distance remains larger with liquid fuel than it is for a dry wall, and maximum heat fluxes are smaller. The nature of the interaction of flame wall interaction with a liquid fuel is also different from the classical flame/dry wall interaction: it is controlled mainly by chemical mechanisms and not by the thermal quenching effect observed for flames interacting with dry walls: the existence of a very rich zone created above the liquid film is the main mechanism controlling quenching.

GPU-accelerated optimization of fuel treatments for mitigating wildfire hazard

Fuel treatment is considered a suitable way to mitigate the hazard related to potential wildfires on a landscape. However, designing an optimal spatial layout of treatment units represents a difficult optimization problem. In fact, budget constraints, the probabilistic nature of fire spread and interactions among the different area units composing the whole treatment, give rise to challenging search spaces on typical landscapes. In this paper we formulate such optimization problem with the objective of minimizing the extension of land characterized by high fire hazard. Then, we propose a computational approach that leads to a spatially-optimized treatment layout exploiting Tabu Search and General-Purpose computing on Graphics Processing Units (GPGPU). Using an application example, we also show that the proposed methodology can provide high-quality design solutions in low computing time. © 2013 The Authors. Published by Elsevier B.V.