Transient response of structures with uncertain properties to nonlinear shock loading

A method is presented to predict the transient response of a structure at the driving point following an impact or a shock loading. The displacement and the contact force are calculated solving the discrete convolution between the impulse response and the contact force itself, expressed in terms of a nonlinear Hertzian contact stiffness. Application of random point process theory allows the calculation of the impulse response function from knowledge of the modal density and the geometric characteristics of the structure only. The theory is applied to a wide range of structures and results are experimentally verified for the case of a rigid object hitting a beam, a plate, a thin and a thick cylinder and for the impact between two cylinders. The modal density of the flexural modes for a thick slender cylinder is derived analytically. Good agreement is found between experimental, simulated and published results, showing the reliability of the method for a wide range of situations including impacts and pyroshock applications. © 2013 Elsevier Ltd. All rights reserved.

Predicting the response of structures to transient shock loading.

This work concerns the prediction of the response of an uncertain structure to a load of short duration. Assuming an ensemble of structures with small random variations about a nominal form, a mean impulse response can be found using only the modal density of the structure. The mean impulse response turns out to be the same as the response of an infinite structure: the response is calculated by taking into account the direct field only, without reflections. Considering the short duration of an impulsive loading, the approach is reasonable before the effect of the reverberant field becomes important. The convolution between the mean impulse response and the shock loading is solved in discrete time to calculate the response at the driving point and at remote points. Experimental and numerical examples are presented to validate the theory presented for simple structures such as beams, plates, and cylinders.

Behaviour of unsteady transonic shock / boundary layer interactions with three-dimensional effects

An experimental investigation into the response of transonic SBLIs to periodic down-stream pressure perturbations in a parallel walled duct has been conducted. Tests have been carried out with a shock strength of M ∞ = 1.5 for pressure perturbation frequencies in the range 16-90 Hz. Analysis of the steady interaction at M∞ = 1.5 has also been made. The principle measurement techniques were high speed schlieren photography and laser Doppler anemometry. The structure of the steady SBLI was found to be highly three-dimensional, with large corner flows and sidewall SBLIs. These aspects are thought to influence the upstream transmission of pressure information through the interaction by affecting the post-shock flow field, including the extent of regions of secondary supersonic flow. At low frequency, the dynamics of shock motion can be predicted using an inviscid analytical model. At increased frequencies, viscous effects become significant and the shock exhibits unexpected dynamic behaviour, due to a phase lag between the upstream transmission of pressure information in the core flow and in the viscous boundary layers. Flow control in the form of micro-vane vortex generators was found to have a small impact on shock dynamics, due to the effect it had on the post-shock flow field outside the viscous boundary layer region. The relationship between inviscid and viscous effects is developed and potential destabilising mechanisms for SBLIs in practical applications are suggested. Copyright © 2009 by Paul Bruce and Holger Babinsky.

Dynamic compressive response of composite square honeycombs

Carbon fibre-epoxy composite square honeycombs, and the parent composite material, were tested in quasi-static compression at a strain rate of 10 -3 s -1 and in dynamic compression at strain rates of 10 3-10 4 s -1 using an instrumented Kolsky bar arrangement. Taken together, these tests provide an assessment of the potential of this composite topology for use as a lightweight sandwich core. The honeycombs had two relative densities, 0.12 and 0.24, and two material orientations, ±45° and 0/90° with respect to the prismatic, loading direction of the honeycomb. Honeycomb manufacture was by slotting, assembling and bonding together carbon fibre/epoxy woven plies of composite sheets of 2 × 2 twill weave construction. The peak value of wall stress in the honeycombs was about one third that of the parent material, for all strain rates. An elastic finite element analysis was used to trace the source of this knock-down in strength: a stress concentration exists at the root of the slots and leads to premature failure by microbuckling. Shock-wave effects were evident at impact velocities exceeding 50 ms -1 for the honeycomb of relative density 0.12. This was traced to stubbing of the buckled cell walls against the face of the Kolsky bar. © 2011 Elsevier Ltd. All rights reserved.

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.

Turbine blade tip heat transfer in low speed and high speed flows

In this paper, high and low speed tip flows are investigated for a high-pressure turbine blade. Previous experimental data are used to validate a CFD code, which is then used to study the tip heat transfer in high and low speed cascades. The results show that at engine representative Mach numbers the tip flow is predominantly transonic. Thus, compared to the low speed tip flow, the heat transfer is affected by reductions in both the heat transfer coefficient and the recovery temperature. The high Mach numbers in the tip region (M>1.5) lead to large local variations in recovery temperature. Significant changes in the heat transfer coefficient are also observed. These are due to changes in the structure of the tip flow at high speed. At high speeds, the pressure side corner separation bubble reattachment occurs through supersonic acceleration which halves the length of the bubble when the tip gap exit Mach number is increased from 0.1 to 1.0. In addition, shock/boundary-layer interactions within the tip gap lead to large changes in the tip boundary-layer thickness. These effects give rise to significant differences in the heat-transfer coefficient within the tip region compared to the low-speed tip flow. Compared to the low speed tip flow, the high speed tip flow is much less dominated by turbulent dissipation and is thus less sensitive to the choice of turbulence model. These results clearly demonstrate that blade tip heat transfer is a strong function of Mach number, an important implication when considering the use of low speed experimental testing and associated CFD validation in engine blade tip design. Copyright © 2009 by ASME.

Experimental investigation of dynamic response of acoustically forced turbulent premixed CH4/CO2/air flames

Increasing demand for energy and continuing increase in environmental as well as financial cost of use of fossil fuels drive the need for utilization of fuels from sustainable sources for power generation. Development of fuel-flexible combustion systems is vital in enabling the use of sustainable fuels. It is also important that these sustainable combustion systems meet the strict governmental emission legislations. Biogas is considered as one of the viable sustainable fuels that can be used to power modern gas turbines: However, the change in chemical, thermal and transport properties as well as change in Wobbe index due to the variation of the fuel constituents can have a significant effect on the performance of the combustor. It is known that the fuel properties have strong influence on the dynamic flame response; however there is a lack of detailed information regarding the effect of fuel compositions on the sensitivity of the flames subjected to flow perturbations. In this study, we describe an experimental effort investigating the response of premixed biogas-air turbulent flames with varying proportions of CH4 and CO2 to velocity perturbations. The flame was stabilized using a centrally placed conical bluff body. Acoustic perturbations were imposed to the flow using loud speakers. The flame dynamics and the local heat release rate of these acoustically excited biogas flames were studied using simultaneous measurements of OH and H2CO planar laser induced fluorescence. OH* chemiluminescence along with acoustic pressure measurements were also recorded to estimate the total flame heat release modulation and the velocity fluctuations. The measurements were carried out by keeping the theoretical laminar flame speed constant while varying the bulk velocity and the fuel composition. The results indicate that the flame sensitivity to perturbations increased with increased dilution of CH4 by CO2 at low amplitude forcing, while at high amplitude forcing conditions the magnitude of the flame response was independent of dilution.

Shock dynamics of random structures.

Predicting the response of a structure following an impact is of interest in situations where parts of a complex assembly may come into contact. Standard approaches are based on the knowledge of the impulse response function, requiring the knowledge of the modes and the natural frequencies of the structure. In real engineering structures the statistics of higher natural frequencies follows those of the Gaussian Orthogonal Ensemble, this allows the application of random point process theory to get a mean impulse response function by the knowledge of the modal density of the structure. An ensemble averaged time history for both the response and the impact force can be predicted. Once the impact characteristics are known in the time domain, a simple Fourier Transform allows the frequency range of the impact excitation to be calculated. Experimental and numerical results for beams, plates, and cylinders are presented to confirm the validity of the method.

Shock transmission in a coupled beam system

This paper investigates the circumstances under which high peak acceleration can occur in the internal parts of a system when subjected to impulsive driving on the outside. Motivating examples include the design of packaging for transportation of fragile items. The system is modelled in an idealised form using two beams coupled with point connections. A Rayleigh-Ritz model of such coupled beams was validated against measurements on a particular beam system, then the model was used to explore the acceleration response to impulsive driving in the time, frequency and spatial domains. This study is restricted to linear vibration response and additional mechanisms for high internal acceleration due to nonlinear effects such as internal impacts are not considered. Using Monte Carlo simulation in which the indirectly driven beam was perturbed by randomly placed point masses a wide range of system behaviour was explored. This facilitates identification of vulnerable configurations that can lead to high internal acceleration. The results from the study indicate the possibility of curve veering influencing the peak acceleration amplification. The possibility of veering within an ensemble was found to be dependent on the relative coupling strength of the modes. Understanding of the mechanism may help to avoid vulnerable cases, either by design or by preparatory vibration testing. © 2013 Elsevier Ltd.

Low-order modelling of the response of ducted flames in annular geometries

The adoption of lean premixed prevaporised combustion systems can reduce NOx emissions from gas turbines, but unfortunately also increases their susceptibility to thermoacoustic instabilities. Initially, acoustic waves can produce heat release fluctuations by a variety of mechanisms, often by perturbing the equivalence ratio. If correctly phased, heat release fluctuations can subsequently generate more acoustic waves, which at high amplitude can result in significant structural damage to the combustor. The prediction of this phenomenon is of great industrial interest. In previous work, we have coupled a physics based, kinematic model of the flame with a network model to provide the planar acoustic response necessary to close the feedback loop and predict the onset and amplitude of thermoacoustic instabilities in a lab-scale, axisymmetric single burner combustor. The advantage of a time domain approach is that the modal interaction, the influence of harmonics, and flame saturation can be investigated. This paper extends this approach to more realistic, annular geometries, where both planar and circumferential modes must be considered. In lean premixed prevaporised combustors, fluctuations in equivalence ratio have been shown to be a dominant cause of unsteady combustion. These can occur, for example, due to velocity perturbations in the premix ducts, which can lead to equivalence ratio fluctuations at the fuel injectors, which are subsequently convected downstream to the flame surfaces. Here, they can perturb the heat release by locally altering the flame speed, enthalpy of combustion, and, indirectly, the flame surface area. In many gas turbine designs, particularly aeroengines, the geometries are composed of a ring of premix ducts linking a plenum and an annular combustor. The most unstable modes are often circumferential modes. The network model is used to characterise the flow response of the geometry to heat fluctuations at an appropriate location, such as the fuel injectors. The heat release at each flame holder is determined in the time domain using the kinematic flame model derived, as a function of the flow perturbations in the premix duct. This approach is demonstrated for an annular ring of burners on a in a simple geometry. The approach is then extended to an industrial type gas turbine combustor, and used to predict the limit cycle amplitudes. Copyright © 2012 by ASME.