Simulation of spherical turbulent premixed flames

Outwardly-propagating spherical turbulent premixed flames are studied using Reynolds-Averaged Navier-Stokes methodology. The reaction rate closure is based on recently developed strained flamelets, where the flamelets are parametrised using the scalar dissipation rate. It is shown that the leading edge speed for the spherical flame is higher than for the planar case, for a given turbulence and thermo-chemical conditions. In addition, it is shown that including the mean curvature effects in the reaction rate closure do not influence the spherical flame speeds when compared with a model that excludes these effects.

Investigation of self excited and acoustically forced combustion instabilities with laser diagnostics

Combustion oscillations in gas turbines can result in serious damage. One method used to predict such oscillations is to analyze the combustor acoustics using a simple linear model. Such a model requires a flame transfer function to describe the response of the heat release to flow perturbations inside the combustor. This paper reports on the application of Planar Laser Induced Fluorescence (PLIF) of OH radicals to analyze the response of a lean premixed flame to oncoming flow perturbations. Both self-excited oscillations and low amplitude forced oscillations at various frequencies are investigated in an atmospheric pressure model combustor rig. In order to visualize fluctuations of local fuel distribution, acetone-PLIF was also applied in non-reacting and acoustically forced flows at oscillation frequencies of 200 Hz and 510 Hz, respectively. OH-PLIF images were acquired over a range of operating parameters. The results presented in this paper originate from data sets acquired at fixed phase angles during the oscillation cycle. Comparative experiments in self excited and forced acoustic oscillations show that the flame and the combustion intensity develop similarly throughout the pressure cycle in both cases. Although the peak fluorescence intensities differ between self excited and the forced instabilities, there is a clear correspondence in the observed frequency and phase information from the two cases. This result encourages a comparison of the OH-PLIF and the acetone-PLIF results. Quantitative measurements of the equivalence ratio in specific areas of the measurement plane offer insight on the complex phenomena coupling acoustic perturbations, i.e. flow velocity fluctuations, to fluctuations in fuel distribution and combustion intensity, ultimately resulting in self excited combustion oscillations.

Task-dependent coordination of rapid bimanual motor responses.

Optimal feedback control postulates that feedback responses depend on the task relevance of any perturbations. We test this prediction in a bimanual task, conceptually similar to balancing a laden tray, in which each hand could be perturbed up or down. Single-limb mechanical perturbations produced long-latency reflex responses ("rapid motor responses") in the contralateral limb of appropriate direction and magnitude to maintain the tray horizontal. During bimanual perturbations, rapid motor responses modulated appropriately depending on the extent to which perturbations affected tray orientation. Specifically, despite receiving the same mechanical perturbation causing muscle stretch, the strongest responses were produced when the contralateral arm was perturbed in the opposite direction (large tray tilt) rather than in the same direction or not perturbed at all. Rapid responses from shortening extensors depended on a nonlinear summation of the sensory information from the arms, with the response to a bimanual same-direction perturbation (orientation maintained) being less than the sum of the component unimanual perturbations (task relevant). We conclude that task-dependent tuning of reflexes can be modulated online within a single trial based on a complex interaction across the arms.

BarraCUDA – a Fast Sequence Mapping Software using Graphics Processing Units

High-throughput DNA sequencing (HTS) instruments today are capable of generating millions of sequencing reads in a short period of time, and this represents a serious challenge to current bioinformatics pipeline in processing such an enormous amount of data in a fast and economical fashion. Modern graphics cards are powerful processing units that consist of hundreds of scalar processors in parallel in order to handle the rendering of high-definition graphics in real-time. It is this computational capability that we propose to harness in order to accelerate some of the time-consuming steps in analyzing data generated by the HTS instruments. We have developed BarraCUDA, a novel sequence mapping software that utilizes the parallelism of NVIDIA CUDA graphics cards to map sequencing reads to a particular location on a reference genome. While delivering a similar mapping fidelity as other mainstream programs , BarraCUDA is a magnitude faster in mapping throughput compared to its CPU counterparts. The software is also capable of supporting multiple CUDA devices in parallel to further accelerate the mapping throughput. BarraCUDA is designed to take advantage of the parallelism of GPU to accelerate the mapping of millions of sequencing reads generated by HTS instruments. By doing this, we could, at least in part streamline the current bioinformatics pipeline such that the wider scientific community could benefit from the sequencing technology. BarraCUDA is currently available at http://seqbarracuda.sf.net

Development of a tool for rapidly assessing the implementation difficulty and emissions benefits of innovations

Consumer goods manufacturers aiming to reduce the environmental impact associated with their products commonly pursue incremental change strategies, but more radical approaches may be required if we are to address the challenges of sustainable consumption. One strategy to realize step change reductions is to prepare a portfolio of innovations providing different levels of impact reduction in exchange for different levels of organizational resource commitment. In this research a tool is developed to support this strategy, starting with the assumption that through brainstorming or other eco-innovation approaches, a long-list of candidate innovations has been created. The tool assesses the potential greenhouse gas benefit of an innovative option against the difficulty of its implementation. A simple greenhouse gas benefit assessment method based on streamlined LCA was used to analyze impact reduction potential, and a novel measure of implementation difficulty was developed. The predictions of implementation difficulty were compared against expert opinion, and showed similar results indicating the measure can be used sensibly to predict implementation difficulty. The assessment of the environmental gain versus implementation difficulty is visualized in a matrix, showing the trade-offs of several options. The tool is deliberately simple with scalar measures of CO 2 emissions benefits and implementation difficulty so tool users must remain aware of other potential environmental burdens besides greenhouse gases (e.g. water, waste). In addition, although relative life cycle emissions benefits of an option may be low, the absolute impact of an option can be high and there may be other co-benefits, which could justify higher levels of implementation difficulty. Different types of consumer products (e.g. household, personal care, foods) have been evaluated using the tool. Initial trials of the tool within Unilever demonstrate that the tool facilitates rapid evaluation of low-carbon innovations. © 2011 Elsevier Ltd. All rights reserved.

Measurements of triggering and transient growth in a model lean-premixed gas turbine combustor

Instability triggering and transient growth of thermoacoustic oscillations were experimentally investigated in combination with linear/nonlinear flame transfer function (FTF) methodology in a model lean-premixed gas turbine combustor operated with CH 4 and air at atmospheric pressure. A fully premixed flame with 10kW thermal power and an equivalence ratio of 0.60 was chosen for detailed characterization of the nonlinear transient behaviors. Flame transfer functions were experimentally determined by simultaneous measurements of inlet velocity fluctuations and heat release rate oscillations using a constant temperature anemometer and OH */CH * chemiluminescence emissions, respectively. The phase-resolved variation of the local flame structure at a limit cycle was measured by planar laser-induced fluorescence of OH. Simultaneous measurements of inlet velocity, OH */CH * emission, and acoustic pressure were performed to investigate the temporal evolution of the system from a stable to a limit cycle operation. This measurement allows us to describe an unsteady instability triggering event in terms of several distinct stages: (i) initiation of a small perturbation, (ii) exponential amplification, (iii) saturation, (iv) nonlinear evolution of the perturbations towards a new unstable periodic state, (v) quasi-steady low-amplitude periodic oscillation, and (vi) fully-developed high-amplitude limit cycle oscillation. Phase-plane portraits of instantaneous inlet velocity and heat release rate clearly show the presence of two different attractors. Depending on its initial position in phase space at infinitesimally small amplitude, the system evolves towards either a high-amplitude oscillatory state or a low-amplitude oscillatory state. This transient phenomenon was analyzed using frequency- and amplitude-dependent damping mechanisms, and compared to subcritical and supercritical bifurcation theories. The results presented in this paper experimentally demonstrate the hypothesis proposed by Preetham et al. based on analytical and computational solutions of the nonlinear G-equation [J. Propul. Power 24 (2008) 1390-1402]. Good quantitative agreement was obtained between measurements and predictions in terms of the conditions for the onset of triggering and the amplitude of triggered combustion instabilities. © 2011 The Combustion Institute.

High order schemes on three-dimensional general polyhedral meshes - Application to the level set method

In this article, we detail the methodology developed to construct arbitrarily high order schemes - linear and WENO - on 3D mixed-element unstructured meshes made up of general convex polyhedral elements. The approach is tailored specifically for the solution of scalar level set equations for application to incompressible two-phase flow problems. The construction of WENO schemes on 3D unstructured meshes is notoriously difficult, as it involves a much higher level of complexity than 2D approaches. This due to the multiplicity of geometrical considerations introduced by the extra dimension, especially on mixed-element meshes. Therefore, we have specifically developed a number of algorithms to handle mixed-element meshes composed of convex polyhedra with convex polygonal faces. The contribution of this work concerns several areas of interest: the formulation of an improved methodology in 3D, the minimisation of computational runtime in the implementation through the maximum use of pre-processing operations, the generation of novel methods to handle complex 3D mixed-element meshes and finally the application of the method to the transport of a scalar level set. © 2012 Global-Science Press.

Effects of flow geometry and Damkohler number on turbulent premixed flame

DNS of planar turbulent flame and turbulent V-flame has been conducted to investigate turbulence-scalar interaction in relatively practical turbulent combustion. Several turbulence quantities are examined for the understandings of fundamental characteristics of flow field in V-flame. Due to the additional turbulence production by the hot-rod, turbulence does not simply decay in V-flame. Turbulence-scalar interaction, scalar alignments with the principal strain rate in other words, is then clarified. The competition of turbulence and dilatation can be found in the conditional PDF of flame normal alignment. The results suggests that the alignment characteristics in high Da flames are applicable to low Da flames in the region of intense heat release.

In vivo dynamics of the internal fibrous structure in smooth adhesive pads of insects.

Many insects with smooth adhesive pads can rapidly enlarge their contact area by centripetal pulls on the legs, allowing them to cope with sudden mechanical perturbations such as gusts of wind or raindrops. The short time scale of this reaction excludes any neuromuscular control; it is thus more likely to be caused by mechanical properties of the pad's specialized cuticle. This soft cuticle contains numerous branched fibrils oriented almost perpendicularly to the surface. Assuming a fixed volume of the water-filled cuticle, we hypothesized that pulls could decrease the fibril angle, thereby helping the contact area to expand laterally and longitudinally. Three-dimensional fluorescence microscopy on the cuticle of smooth stick insect pads confirmed that pulls significantly reduced the fibril angle. However, the fibril angle variation appeared insufficient to explain the observed increase in contact area. Direct strain measurements in the contact zone demonstrated that pulls not only expand the cuticle laterally, but also add new contact area at the pad's outer edge.

Validation of a lattice Boltzmann model for gas-solid reactions with experiments

A lattice Boltzmann method is used to model gas-solid reactions where the composition of both the gas and solid phase changes with time, while the boundary between phases remains fixed. The flow of the bulk gas phase is treated using a multiple relaxation time MRT D3Q19 model; the dilute reactant is treated as a passive scalar using a single relaxation time BGK D3Q7 model with distinct inter- and intraparticle diffusivities. A first-order reaction is incorporated by modifying the method of Sullivan et al. [13] to include the conversion of a solid reactant. The detailed computational model is able to capture the multiscale physics encountered in reactor systems. Specifically, the model reproduced steady state analytical solutions for the reaction of a porous catalyst sphere (pore scale) and empirical solutions for mass transfer to the surface of a sphere at Re=10 (particle scale). Excellent quantitative agreement between the model and experiments for the transient reduction of a single, porous sphere of Fe 2O 3 to Fe 3O 4 in CO at 1023K and 10 5Pa is demonstrated. Model solutions for the reduction of a packed bed of Fe 2O 3 (reactor scale) at identical conditions approached those of experiments after 25 s, but required prohibitively long processor times. The presented lattice Boltzmann model resolved successfully mass transport at the pore, particle and reactor scales and highlights the relevance of LB methods for modelling convection, diffusion and reaction physics. © 2012 Elsevier Inc.

Low-order modelling of ducted flames with temporally varying equivalence ratio in realistic geometries

The interaction between unsteady heat release and acoustic pressure oscillations in gas turbines results in self-excited combustion oscillations which can potentially be strong enough to cause significant structural damage to the combustor. Correctly predicting the interaction of these processes, and anticipating the onset of these oscillations can be difficult. In recent years much research effort has focused on the response of premixed flames to velocity and equivalence ratio perturbations. In this paper, we develop a flame model based on the socalled G-Equation, which captures the kinematic evolution of the flame surfaces, under the assumptions of axisymmetry, and ignoring vorticity and compressibility. This builds on previous work by Dowling [1], Schuller et al. [2], Cho & Lieuwen [3], among many others, and extends the model to a realistic geometry, with two intersecting flame surfaces within a non-uniform velocity field. The inputs to the model are the free-stream velocity perturbations, and the associated equivalence ratio perturbations. The model also proposes a time-delay calculation wherein the time delay for the fuel convection varies both spatially and temporally. The flame response from this model was compared with experiments conducted by Balachandran [4, 5], and found to show promising agreement with experimental forced case. To address the primary industrial interest of predicting self-excited limit cycles, the model has then been linked with an acoustic network model to simulate the closed-loop interaction between the combustion and acoustic processes. This has been done both linearly and nonlinearly. The nonlinear analysis is achieved by applying a describing function analysis in the frequency domain to predict the limit cycle, and also through a time domain simulation. In the latter case, the acoustic field is assumed to remain linear, with the nonlinearity in the response of the combustion to flow and equivalence ratio perturbations. A transfer function from unsteady heat release to unsteady pressure is obtained from a linear acoustic network model, and the corresponding Green function is used to provide the input to the flame model as it evolves in the time domain. The predicted unstable frequency and limit cycle are in good agreement with experiment, demonstrating the potential of this approach to predict instabilities, and as a test bench for developing control strategies. Copyright © 2011 by ASME.

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.

Flux pumping, fluctuations and climbing fields

This paper describes the behaviour of bulk superconductors when subjected to a varying magnetic field. A magnetic model is described together with experimental results which explain and describe the behaviour of superconducting bulks when subjected to varying magnetic fields. We demonstrate how the behaviour is dependent on the magnitude and period of the perturbations in the fields. The model which we use has been implemented using the Comsol™pde solver. It is a fully integrated model which uses a variable heat source to regulate the magnetic circuit and thereby to achieve flux pumping. Comsol™is used for post solution visualization and the model is presented alongside experimental results which support and confirm the conclusions from the model. © 2012 IOP Publishing Ltd.

Conditional Moment Closure for Turbulent Premixed Flames

The conditional moment closure (CMC) method has been successfully applied to various non-premixed combustion systems in the past, but its application to premixed flames is not fully tested and validated. The main difficulty is associated with the modeling of conditional scalar dissipation rate of the conditioning scalar, the progress variable. A simple algebraic model for the conditional dissipation rate is validated using DNS results of a V-flame. This model along with the standard k- turbulence modeling is used in computations of stoichiometric pilot stabilized Bunsen flames using the RANS-CMC method. A first-order closure is used for the conditional mean reaction rate. The computed non reacting and reacting scalars are in reasonable agreement with the experimental measurements and are consistent with earlier computations using flamelets and transported PDF methods. Sensitivity to chemical kinetic mechanism is also assessed. The results suggest that the CMC may be applied across the regimes of premixed combustion.