The structure of turbulent stratified and premixed methane/air flames I: Non-swirling flows

A series of flames in a turbulent methane/air stratified swirl burner is presented. The degree of stratification and swirl are systematically varied to generate a matrix of experimental conditions, allowing their separate and combined effects to be investigated. Non-swirling flows are considered in the present paper, and the effects of swirl are considered in a companion paper (Part II). A mean equivalence ratio of φ=0.75 is used, with φ for the highest level of stratification spanning 0.375-1.125. The burner features a central bluff-body to aid flame stabilization, and the influence of the induced recirculation zone is also considered. The current work focuses on non-swirling flows where two-component particle image velocimetry (PIV) measurements are sufficient to characterize the main features of the flow field. Scalar data obtained from Rayleigh/Raman/CO laser induced fluorescence (CO-LIF) line measurements at 103μm resolution allow the behavior of key combustion species-CH 4, CO 2, CO, H 2, H 2O and O 2-to be probed within the instantaneous flame front. Simultaneous cross-planar OH-PLIF is used to determine the orientation of the instantaneous flame normal in the scalar measurement window, allowing gradients in temperature and progress variable to be angle corrected to their three dimensional values. The relationship between curvature and flame thickness is investigated using the OH-PLIF images, as well as the effect of stratification on curvature.The main findings are that the behavior of the key combustion species in temperature space is well captured on the mean by laminar flame calculations regardless of the level of stratification. H 2 and CO are significant exceptions, both appearing at elevated levels in the stratified flames. Values for surface density function and by extension thermal scalar dissipation rate are found to be substantially lower than laminar values, as the thickening of the flame due to turbulence dominates the effect of increased strain. These findings hold for both premixed and stratified flames. The current series of flames is proposed as an interesting if challenging set of test cases for existing and emerging turbulent flame models, and data are available on request. © 2012 The Combustion Institute.

The structure of turbulent stratified and premixed methane/air flames II: Swirling flows

Experimental results are presented from a series of turbulent methane/air stratified flames stabilized on a swirl burner. Nine operating conditions are considered, systematically varying the level of stratification and swirl while maintaining a lean global mean equivalence ratio of φ̄=0.75. Scalar data are obtained from Rayleigh/Raman/CO laser induced fluorescence (CO-LIF) line measurements at 103μm resolution, allowing the behavior of the major combustion species-CH 4, CO 2, CO, H 2, H 2O and O 2-to be probed within the instantaneous flame front. The corresponding three-dimensional surface density function and thermal scalar dissipation rate are investigated, along with geometric characteristics of the flame such as curvature and flame thickness. Hydrogen and carbon monoxide levels within the flame brush are raised by stratification, indicating models with laminar premixed flame chemistry may not be suitable for stratified flames. However, flame surface density, scalar dissipation and curvature all appear insensitive to the degree of stratification in the flames surveyed. © 2012 The Combustion Institute.

Efficient truncation scheme for modeling deepwater mooring lines

This paper is aimed at enabling the confident use of existing model test facilities for ultra deepwater application without having to compromise on the widely accepted range of scales currently used by the floating production industry. Passive line truncation has traditionally been the preferred method of creating an equivalent numerical model at reduced depth; however, these techniques tend to suffer in capturing accurately line dynamic response and so reproducing peak tensions. In an attempt to improve credibility of model test data the proposed truncation procedure sets up the truncated model, based on line dynamic response rather than quasi-static system stiffness. The upper sections of each line are modeled in detail, capturing the wave action zone and all coupling effects with the vessel. These terminate to an approximate analytical model that aims to simulate the remainder of the line. Stages 1 & 2 are used to derive a water depth truncation ratio. Here vibration decay of transverse elastic waves is assessed and it is found that below a certain length criterion, the transverse vibrational characteristics for each line are inertia driven, hence with respect to these motions the truncated model can assume a linear damper whose coefficient depends on the local line properties and vibration frequency. Stage 3 endeavors to match the individual line stiffness between the full depth and truncated models. In deepwater it is likely that taut polyester moorings will be used which are predominantly straight and have high axial stiffness that provides the principal restoring force to static and low frequency vessel motions. Consequently, it means that the natural frequencies of axial vibrations are above the typical wave frequency range allowing for a quasi-static solution. In cases of exceptionally large wave frequency vessel motions, localized curvature at the chain seabed segment and tangential skin drag on the polyester rope can increase dynamic peak tensions considerably. The focus of this paper is to develop an efficient scheme based on analytic formulation, for replicating these forces at the truncation. The paper will close with an example case study of a single mooring under extreme conditions that replicates exactly the static and dynamic characteristics of the full depth line. Copyright © 2012 by the International Society of Offshore and Polar Engineers (ISOPE).

Growth and shape control of disks by bending and extension

Differential growth of thin elastic bodies furnishes a surprisingly simple explanation of the complex and intriguing shapes of many biological systems, such as plant leaves and organs. Similarly, inelastic strains induced by thermal effects or active materials in layered plates are extensively used to control the curvature of thin engineering structures. Such behaviour inspires us to distinguish and to compare two possible modes of differential growth not normally compared to each other, in order to reveal the full range of out-of-plane shapes of an initially flat disk. The first growth mode, frequently employed by engineers, is characterised by direct bending strains through the thickness, and the second mode, mainly apparent in biological systems, is driven by extensional strains of the middle surface. When each mode is considered separately, it is shown that buckling is common to both modes, leading to bistable shapes: growth from bending strains results in a double-curvature limit at buckling, followed by almost developable deformation in which the Gaussian curvature at buckling is conserved; during extensional growth, out-of-plane distortions occur only when the buckling condition is reached, and the Gaussian curvature continues to increase. When both growth modes are present, it is shown that, generally, larger displacements are obtained under in-plane growth when the disk is relatively thick and growth strains are small, and vice versa. It is also shown that shapes can be mono-, bi-, tri- or neutrally stable, depending on the growth strain levels and the material properties: furthermore, it is shown that certain combinations of growth modes result in a free, or natural, response in which the doubly curved shape of disk exactly matches the imposed strains. Such diverse behaviour, in general, may help to realise more effective actuation schemes for engineering structures. © 2012 Elsevier Ltd. All rights reserved.

Stiffness of sands through a laboratory test database

Deformations of sandy soils around geotechnical structures generally involve strains in the range small (0·01%) to medium (0·5%). In this strain range the soil exhibits non-linear stress-strain behaviour, which should be incorporated in any deformation analysis. In order to capture the possible variability in the non-linear behaviour of various sands, a database was constructed including the secant shear modulus degradation curves of 454 tests from the literature. By obtaining a unique S-shaped curve of shear modulus degradation, a modified hyperbolic relationship was fitted. The three curve-fitting parameters are: an elastic threshold strain γe, up to which the elastic shear modulus is effectively constant at G0; a reference strain γr, defined as the shear strain at which the secant modulus has reduced to 0·5G0; and a curvature parameter a, which controls the rate of modulus reduction. The two characteristic strains γe and γr were found to vary with sand type (i.e. uniformity coefficient), soil state (i.e. void ratio, relative density) and mean effective stress. The new empirical expression for shear modulus reduction G/G0 is shown to make predictions that are accurate within a factor of 1·13 for one standard deviation of random error, as determined from 3860 data points. The initial elastic shear modulus, G0, should always be measured if possible, but a new empirical relation is shown to provide estimates within a factor of 1·6 for one standard deviation of random error, as determined from 379 tests. The new expressions for non-linear deformation are easy to apply in practice, and should be useful in the analysis of geotechnical structures under static loading.

Multiply conditioned analyses of stratification in highly swirling methane/air flames

The effects of stratification on a series of highly swirling turbulent flames under globally lean conditions (φg=0.75) are investigated using a new high-spatial resolution multi-scalar dataset. This dataset features two key properties: high spatial resolution which approaches the 60 micron optical limit of the measurement system, and a wavelet oversampling methodology which significantly reduces the influence of noise. Furthermore, the very large number of realizations (30,000) acquired in the stratified cases permits statistically significant results to be obtained even after aggressive conditioning is applied. Data are doubly conditioned on equivalence ratio and the degree of stratification across the flame in each instantaneous realization. The influence of stoichiometry is limited by conditioning on the equivalence ratio at the location of peak CO mass fraction, which is shown to be a good surrogate for the location of peak heat release rate, while the stratification is quantified using a linear gradient in equivalence ratio across the instantaneous flame front. This advanced conditioning enables robust comparisons with the baseline lean premixed flame. Species mass fractions of both carbon monoxide and hydrogen are increased in temperature space under stratified conditions. Stratification is also shown to significantly increase thermal gradients, yet the derived three-dimensional flame surface density is shown to be relatively insensitive to stratification. Whilst the presence of instantaneous stratification broadens the curvature distribution relative to the premixed case, the degree of broadening is not significantly influenced by the range of global stratification ratios examined in this study. © 2012 The Combustion Institute.

Computation of transient and steady turbulent isothermal flows with an impinging rectangular jet

The computation of both transient and steady turbulent incompressible isothermal flows is studied. The flow is very complex, having streamline curvature, large vortex structures and stagnation resulting from an impinging rectangular jet. For transient computations, the standard k-ε model is adopted. For steady flows, the k-ε, high and low Reynolds number k-l and mixing length models are tried. Zonal approaches combining the above turbulence models are also investigated. None of the models are found to give satisfactory agreement with velocity measurements.

Bifurcations and instability in the adhesion of intrinsically curved rods

Motivated by applications such as gecko-inspired adhesives and microdevices featuring slender rod-like bodies, there has been an increase in interest in the deformed shapes of elastic rods adhering to rigid surfaces. A central issue in analyses of the rod-based models for these systems is the stability of the predicted equilibrium configurations. Such analyses can be complicated by the presence of intrinsic curvatures induced by fabrication processes. The results in the present paper are used to show how this curvature can lead to shear-induced bifurcations and instabilities. To characterize potential instabilities, a new set of necessary conditions for stability are employed which cater to the possible combinations of buckling and delaminating instabilities. © 2013 Elsevier Ltd. All rights reserved.

Blade Loading and Its Application in the Mean-Line Design of Low Pressure Turbines

In order to minimize the number of iterations to a turbine design, reasonable choices of the key parameters must be made at the preliminary design stage. The choice of blade loading is of particular concern in the low pressure (LP) turbine of civil aero engines, where the use of high-lift blades is widespread. This paper considers how blade loading should be measured, compares the performance of various loss correlations, and explores the impact of blade lift on performance and lapse rates. To these ends, an analytical design study is presented for a repeating-stage, axial-flow LP turbine. It is demonstrated that the long-established Zweifel lift coefficient (Zweifel, 1945, "The Spacing of Turbomachine Blading, Especially with Large Angular Deflection" Brown Boveri Rev., 32(1), pp. 436-444) is flawed because it does not account for the blade camber. As a result the Zweifel coefficient is only meaningful for a fixed set of flow angles and cannot be used as an absolute measure of blade loading. A lift coefficient based on circulation is instead proposed that accounts for the blade curvature and is independent of the flow angles. Various existing profile and secondary loss correlations are examined for their suitability to preliminary design. A largely qualitative comparison demonstrates that the loss correlations based on Ainley and Mathieson (Ainley and Mathieson, 1957, "A Method of Performance Estimation for Axial-Flow Turbines," ARC Reports and Memoranda No. 2974; Dunham and Came, 1970, "Improvements to the Ainley-Mathieson Method of Turbine Performance Prediction," Trans. ASME: J. Eng. Gas Turbines Power, July, pp. 252-256; Kacker and Okapuu, 1982, "A Mean Line Performance Method for Axial Flow Turbine Efficiency," J. Eng. Power, 104, pp. 111-119). are not realistic, while the profile loss model of Coull and Hodson (Coull and Hodson, 2011, "Predicting the Profile Loss of High-Lift Low Pressure Turbines," J. Turbomach., 134(2), pp. 021002) and the secondary loss model of (Traupel, W, 1977, Thermische Turbomaschinen, Springer-Verlag, Berlin) are arguably the most reasonable. A quantitative comparison with multistage rig data indicates that, together, these methods over-predict lapse rates by around 30%, highlighting the need for improved loss models and a better understanding of the multistage environment. By examining the influence of blade lift across the Smith efficiency chart, the analysis demonstrates that designs with higher flow turning will tend to be less sensitive to increases in blade loading. © 2013 American Society of Mechanical Engineers.

Prediction of FRP debonding using the global-energy-balance approach

A major research program was carried out to analyze the mechanism of FRP debonding from concrete beams using global-energy-balance approach (GEBA). The key findings are that the fracture process zone is small so there is no R-curve to consider, failure is dominated by Mode I behavior, and the theory agrees well with tests. The analyses developed in the study provide an essential tool that will enable fracture mechanics to be used to determine the load at which FRP plates will debond from concrete beams. This obviates the need for finite element (FE) analyses in situations where reliable details of the interface geometry and crack-tip stress fields are not attainable for an accurate analysis. This paper presents an overview of the GEBA analyses that is described in detail elsewhere, and explains the slightly unconventional assumptions made in the analyses, such as the revised moment-curvature model, the location of an effective centroid, the separate consideration of the FRP and the RC beam for the purposes of the analysis, the use of Mode I fracture energies and the absence of an R-curve in the fracture mechanics analysis.

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.

Collective motion of self-propelled particles: Stabilizing symmetric formations on closed curves

We provide feedback control laws to stabilize formations of multiple, unit speed particles on smooth, convex, and closed curves with definite curvature. As in previous work we exploit an analogy with coupled phase oscillators to provide controls which isolate symmetric particle formations that are invariant to rigid translation of all the particles. In this work, we do not require all particles to be able to communicate; rather we assume that inter-particle communication is limited and can be modeled by a fixed, connected, and undirected graph. Because of their unique spectral properties, the Laplacian matrices of circulant graphs play a key role. The methodology is demonstrated using a superellipse, which is a type of curve that includes circles, ellipses, and rounded rectangles. These results can be used in applications involving multiple autonomous vehicles that travel at constant speed around fixed beacons. ©2006 IEEE.

α-Synuclein senses lipid packing defects and induces lateral expansion of lipids leading to membrane remodeling.

There is increasing evidence for the involvement of lipid membranes in both the functional and pathological properties of α-synuclein (α-Syn). Despite many investigations to characterize the binding of α-Syn to membranes, there is still a lack of understanding of the binding mode linking the properties of lipid membranes to α-Syn insertion into these dynamic structures. Using a combination of an optical biosensing technique and in situ atomic force microscopy, we show that the binding strength of α-Syn is related to the specificity of the lipid environment (the lipid chemistry and steric properties within a bilayer structure) and to the ability of the membranes to accommodate and remodel upon the interaction of α-Syn with lipid membranes. We show that this interaction results in the insertion of α-Syn into the region of the headgroups, inducing a lateral expansion of lipid molecules that can progress to further bilayer remodeling, such as membrane thinning and expansion of lipids out of the membrane plane. We provide new insights into the affinity of α-Syn for lipid packing defects found in vesicles of high curvature and in planar membranes with cone-shaped lipids and suggest a comprehensive model of the interaction between α-Syn and lipid bilayers. The ability of α-Syn to sense lipid packing defects and to remodel membrane structure supports its proposed role in vesicle trafficking.

Compact multimode polymer waveguide bends for board-level optical interconnects

Multimode polymer waveguides are promising for use in board-level optical interconnects. In recent years, various on-board optical interconnection architectures have been demonstrated making use of passive routing waveguide components. In particular, 90° bends have played important roles in complex waveguide layouts enabling interconnection between non co-linear points on a board. Due to the dimensions and index step of the waveguides typically used in on-board optical interconnects, low-loss bends are typically limited to a radius of ∼ 10 mm. This paper therefore presents the design and fabrication of compact low-loss waveguide bends with reduced radii of curvature, offering significant reductions in the required areas for on-board optical circuits. The proposed design relies on the exposure of the bend section to the air, achieving tighter light confinement along the bend and reduced bending losses. Simulation studies carried out with ray tracing tools and experimental results from polymer samples fabricated on FR4 are presented. Low bending losses are achieved from the air-exposed bends up to 4 mm of radius of curvature, while an improvement of 14 μm in the 1 dB alignment tolerances at the input of these devices (fibre to waveguide coupling) is also obtained. Finally, the air-exposed bends are employed in an optical bus structure, offering reductions in insertion loss of up to 3.8 dB. © 2013 IEEE.

The shape of helically creased cylinders

Creasing in thin shells admits large deformation by concentrating curvatures while relieving stretching strains over the bulk of the shell: after unloading, the creases remain as narrow ridges and the rest of the shell is flat or simply curved. We present a helically creased unloaded shell that is doubly curved everywhere, which is formed by cylindrically wrapping a flat sheet with embedded foldlines not axially aligned. The finished shell is in a state of uniform self-stress and this is responsible for maintaining the Gaussian curvature outside of the creases in a controllable and persistent manner. We describe the overall shape of the shell using the familiar geometrical concept of a Mohr's circle applied to each of its constituent features-the creases, the regions between the creases, and the overall cylindrical form. These Mohr's circles can be combined in view of geometrical compatibility, which enables the observed shape to be accurately and completely described in terms of the helical pitch angle alone. Copyright © 2013 by ASME.