Vibration of a high-altitude tethered balloon with application to climate engineering

A balloon tethered at an altitude of 20 km could deliver a particulate cloud leading to global cooling. Tethering a balloon at this altitude poses significant problems with respect to vibration and stability, especially in regions of high wind. No-one has ever proposed, yet alone launched, a balloon at an altitude of 20 km tethered to the ground. Owing to wind, the tether needs to be 23 km in length and is to be fixed to a ship at sea or on land in equatorial regions. Whilst the balloon at 20 km is subject to relatively modest wind conditions, at jet stream altitudes (10km) the tether will experience much higher wind loadings, not only because of the high wind speeds of up to 300 km / hr but also because of the high air density. A tether of circular cross section in these high winds will be subject to horizontal and downward drag forces that would bring the aerostat down. For this reason it is advantageous to consider a self-aligning tether of an aerodynamic cross section whereby it is possible to reduce the drag substantially. One disadvantage of a non-circular tether is the possibility of flutter and galloping instabilities. It is reasonably straightforward to model these phenomena for short lengths of aerofoil, but the situation becomes more complex for a 20 km tensioned tether with large deflection and curvature, variable wind speed, variable air density and variable tension. Analysis using models of infinite length are used to establish the stability at a local scale where the tension, aerodynamic and geometric properties are considered constant. Dispersion curve analysis is useful here. But for dynamics on a long-wavelength scale (several km) then a full non-linear analysis is required. This non-linear model can be used to establish the local values of tension appropriate for the dispersion analysis. This keynote presentation will give some insight into these issues.

The modes of a tethered, spherical balloon

Observations of a tethered meteorological balloon show that a strong vibration coupling exists between axial forcing of the tether and ovalling deformations of the balloon. Such coupling may lead to system instabilities and fatigue failure in a tethered-balloon system. This is particularly relevant in the case of a balloon launched from a moving vessel, as is proposed as part of the SPICE geoengineering project. This paper investigates the vibration characteristics of a tethered, spherical balloon using a simple analytical model: a tensioned, spherical membrane attached to a spring. The analytical solution for the natural frequencies and modeshapes of this system is compared to transfer functions obtained by laser vibrometry. These results are then used to determine the most suitable method of modelling the dynamic response of a tethered balloon.

Modeling of soft clays response for seismic triggering of submarine slides

Submarine landslides pose considerable hazards to coastal communities and offshore structures. The difficulty and cost of obtaining undisturbed samples of offshore soils for determining material properties required for slope stability analyses contribute to the complexity of the problem. There are significant advantages in using a simplified model for the seismic response of submarine slopes, compatible with the limited amount of information that can be realistically gathered, but still able to capture the key elements of clay behavior. This paper illustrates the process of parameter determination and calibration of the SIMPLE DSS model, developed for the study of seismic triggering of submarine slope instabilities. The selection of parameters and predictions of monotonic and cyclic simple shear response are carried out for Boston Blue Clay, a marine clay extensively studied and with a large experimental database available in the literature. The results show that the simplified model is able to reproduce the important trends in the response of the soil, especially in accounting for the effect of the slope.

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.

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.

On Mixing and Shaking : Structure and Dynamics of Turbulent Stratified Flames

The drive for low emission combustion systems encourages applications using premixed flames. Yet in many applications, considerations of flame stability or mixing times lead to systems with neither premixed nor diffusion flames, which are often called technically premixed or stratified flames. In this talk we discuss the current state of understanding of the effect of mixing and extent of stratification on the structure, microstructure and dynamics of selected turbulent stratified flames. Over the past few years, a significant database of scalar and velocity data has been built to analyze the effects of unmixedness on local and global flame structure. Microscale studies of the flame structures show in detail how the effect of local stratification affects (or not!) the flame structure, flame surface density and scalar dissipation rates, and production of selected species. The experiments place exacting demands on current spectroscopic diagnostics, and reveal the progress and limits to our understanding of turbulent flames in general. The dynamics of stratified flames with respect to instabilities is also shown to be very rich, as the particular shape of the flames and the stabilization points are is significantly affected by the fuel distribution, modifying the rate and location of heat release, and thus the coupling with the surrounding acoustics and determining the onset of self-excitations.

Carbonating magnesia for soil stabilization

This paper investigates the potential for carbonating reactive magnesia (MgO) to serve as a more sustainable soil stabilization method by providing rapid and significant strength development of the stabilized soil through absorbing substantial quantities of CO2. Gaseous CO2 was forced through laboratory-prepared reactive MgO-treated soil samples in a triaxial cell set-up, and their resulting mechanical and microstructural properties were investigated using unconfined compressive strength, X-ray diffraction, and scanning electron microscopy. The results showed that adequately carbonated MgO-treated soils could, in a few hours, reach a similar strength range to corresponding 28 day Portland cement (PC)-stabilized soils. Hydrated magnesium carbonates, namely nesquehonite and hydromagnesite-dypingite, were the main products of the carbonated MgO in the soil, and were responsible for the significant strength development.

Instability of the steady-state regime of generation in laser diodes with a broad gain spectrum

The laser-diode parameters at which the steady-state regime of generation becomes unstable are analyzed within the framework of the mode-locking model. The crucial role of the transverse inhomogeneity of the field, pumping intensity, and spectrum width in developing the instabilities of the steady-state regime of generation is demonstrated. The calculated values of the instability threshold are shown to be consistent with the experimental results. © 2008 Springer Science+Business Media, Inc.

CFD investigation of turbulent premixed flame response to transverse forcing

Screech is a high frequency oscillation that is usually characterized by instabilities caused by large-scale coherent flow structures in the wake of bluff-body flameholders and shear layers. Such oscillations can lead to changes in flame surface area which can cause the flame to burn unsteadily, but also couple with the acoustic modes and inherent fluid-mechanical instabilities that are present in the system. In this study, the flame response to hydrodynamic oscillations is analyzed in a controlled manner using high-fidelity Computational Fluid Dynamics (CFD) with an unsteady Reynolds-averaged Navier-Stokes approach. The response of a premixed flame with and without transverse velocity forcing is analyzed. When unforced, the flame is shown to exhibit a self-excitation that is attributed to the anti-symmetric shedding of vortices in the wake of the flameholder. The flame is also forced using two different kinds of low-amplitude out-of-phase inlet velocity forcing signals. The first forcing method is harmonic forcing with a single characteristic frequency, while the second forcing method involves a broadband forcing signal with frequencies in the range of 500 - 1000 Hz. For the harmonic forcing method, the flame is perturbed only lightly about its mean position and exhibits a limit cycle oscillation that is characteristic of the forcing frequency. For the broadband forcing method, larger changes in the flame surface area and detachment of the flame sheet can be seen. Transition to a complicated trajectory in the phase space is observed. When analyzed systematically with system identification methods, the CFD results, expressed in the form of the Flame Transfer Function (FTF) are capable of elucidating the flame response to the imposed perturbation. The FTF also serves to identify, both spatially and temporally, regions where the flame responds linearly and nonlinearly. Locking-in between the flame's natural self-excited frequency and the subharmonic frequencies of the broadband forcing signal is found to alter the dynamical behaviour of the flame. Copyright © 2013 by ASME.

Properties of two model soils stabilized with different blends and contents of GGBS, MgO, lime, and PC

This paper addresses the use of ground granulated blast furnace slag (GGBS) and reactive magnesia (MgO) blends for soil stabilization, comparing them with GGBS-lime blends and Portland cement (PC) for enhanced technical performance. A range of tests were conducted to investigate the properties of stabilized soils, including unconfined compressive strength (UCS), permeability, and microstructural analyses by using X-ray diffraction (XRD) and scanning electron microscopy (SEM). The influence of GGBS:MgO ratio, binder content, soil type, and curing period were addressed. The UCS results revealed that GGBS-MgO was more efficient than GGBS-lime as a binder for soil stabilization, with an optimum MgO content in the range of 5-20% of the blends content, varying with binder content and curing age. The 28-day UCS values of the optimum GGBS-MgO mixes were up to almost four times higher than that of corresponding PC mixes. The microstructural analyses showed the hydrotalcite was produced during the GGBS hydration activated by MgO, although the main hydration products of the GGBS-MgO stabilized soils were similar to those of PC. © 2014 American Society of Civil Engineers.

Global linear stability of the boundary-layer flow over a rotating sphere

We consider the linear global stability of the boundary-layer flow over a rotating sphere. Our results suggest that a self-excited linear global mode can exist when the sphere rotates sufficiently fast, with properties fixed by the flow at latitudes between approximately 55°-65° from the pole (depending on the rotation rate). A neutral curve for global linear instabilities is presented with critical Reynolds number consistent with existing experimentally measured values for the appearance of turbulence. The existence of an unstable linear global mode is in contrast to the literature on the rotating disk, where it is expected that nonlinearity is required to prompt the transition to turbulence. Despite both being susceptible to local absolute instabilities, we conclude that the transition mechanism for the rotating-sphere flow may be different to that for the rotating disk. © 2014 Elsevier Masson SAS. All rights reserved.

Properties and microstructure of GGBS-magnesia pastes

Reactive magnesia (MgO) was used as an alkali activator for ground granulated blast-furnace slag (GGBS) and its activating efficiency was investigated compared with hydrated lime. GGBS-MgO and GGBS-hydrated lime paste samples with different compositions and different water to solid ratios were prepared and cured for different periods. A range of tests was conducted to investigate the properties and microstructure of the pastes, including compressive strength, X-ray diffraction, scanning electron microscopy, Fourier transform infrared spectroscopy, energy dispersive X-ray and thermogravimetric analysis. The results showed that the reactive MgO acts as an effective alkali activator of GGBS, achieving higher 28-day compressive strength than that of the corresponding GGBS-hydrated lime system. The extensive microstructural investigation indicated that the main hydration product of reactive MgO-activated GGBS and hydrated lime-activated GGBS systems was hydrated calcium silicate, but there was much more hydrotalcite present in the former, which contributed to its superior 28-day compressive strength.

Comparison of acoustic velocity perturbation measurements using PIV vs. two-microphone technique

Understanding combustion instabilities requires accurate measurements of the acoustic velocity perturbation into injectors. This is often accomplished via the use of the two microphone technique, as this only requires two pressure transducers. However, measurements of the actual velocities emerging from the injectors are not often taken, leaving questions regarding the assumptions about the acoustic velocity. A comparison of velocity measured at downstream of the injector with that of two-microphone technique can show the accuracy and limitations of two-microphone technique. In this paper, velocity measurements are taken using both particle image velocimetry (PIV) and the two-microphone technique in a high pressure facility designed for aeroengine injector measurements. The flow is excited using an area modulation device installed on the choked end of the combustion chamber, with PIV measurements enabled by optical access downstream of the injector through a quartz tube and windows. Acoustic velocity perturbations at the injector are determined by considering the Fourier transformed pressure fluctuations for two microphones installed at a known distance upstream of the injector. PIV measurements are realized by seeding the air flow with micrometric water particles under 2.5 bar pressure at ambient temperature. Phase locked velocity fields are realized by synchronizing the acquisition of PIV images with the revolution of the acoustic modulator using the pressure signal measured at the face of injector. The mean velocity fluctuation is calculated as the difference between maximum and minimum velocities, normalized by the mean velocity of the unforced case. Those values are compared with the peak-to-peak velocity fluctuation amplitude calculated by the two-microphone technique. Although the ranges of velocity fluctuations for both techniques are similar, the variation of fluctuation with forcing frequencies diverges significantly with frequency. The differences can be attributed to several limitations associated with of both techniques, such as the quality of the signal, the signal/noise ratio, the accuracy of PIV measurements and the assumption of isentropic flow of the particle velocity from the plenum through the injector. We conclude that two-microphone methods can be used as a reference value for the velocity fluctuation in low order applications such as flame transfer functions, but not for drawing conclusions regarding the absolute velocity fluctuations in the injector. Copyright © 2013 by ASME.

Monitoring and preventing numerical oscillations in 3D simulations with coupled Monte Carlo codes

Previous studies have reported that different schemes for coupling Monte Carlo (MC) neutron transport with burnup and thermal hydraulic feedbacks may potentially be numerically unstable. This issue can be resolved by application of implicit methods, such as the stochastic implicit mid-point (SIMP) methods. In order to assure numerical stability, the new methods do require additional computational effort. The instability issue however, is problem-dependent and does not necessarily occur in all cases. Therefore, blind application of the unconditionally stable coupling schemes, and thus incurring extra computational costs, may not always be necessary. In this paper, we attempt to develop an intelligent diagnostic mechanism, which will monitor numerical stability of the calculations and, if necessary, switch from simple and fast coupling scheme to more computationally expensive but unconditionally stable one. To illustrate this diagnostic mechanism, we performed a coupled burnup and TH analysis of a single BWR fuel assembly. The results indicate that the developed algorithm can be easily implemented in any MC based code for monitoring of numerical instabilities. The proposed monitoring method has negligible impact on the calculation time even for realistic 3D multi-region full core calculations. © 2014 Elsevier Ltd. All rights reserved.

Thermal shock resistance of air plasma sprayed thermal barrier coatings

The spallation resistance of an air plasma sprayed (APS) thermal barrier coating (TBC) to cool-down/reheat is evaluated for a pre-existing delamination crack. The delamination emanates from a vertical crack through the coating and resides at the interface between coating and underlying thermally grown oxide layer (TGO). The coating progressively sinters during engine operation, and this leads to a depth-dependent increase in modulus. Following high temperature exposure, the coating is subjected to a cooling/reheating cycle representative of engine shut-down and start-up. The interfacial stress intensity factors are calculated for the delamination crack over this thermal cycle and are compared with the mode-dependent fracture toughness of the interface between sintered APS and TGO. The study reveals the role played by microstructural evolution during sintering in dictating the spallation life of the thermal barrier coating, and also describes a test method for the measurement of delamination toughness of a thin coating. © 2014 Elsevier Ltd.