Spatial correlation of heat release rate and sound emission from turbulent premixed flames

The two-point spatial correlation of the rate of change of fluctuating heat release rate is central to the sound emission from open turbulent flames, and a few attempts have been made to address this correlation in recent studies. In this paper, the two-point correlation and its role in combustion noise are studied by analysing direct numerical simulation (DNS) data of statistically multi-dimensional turbulent premixed flames. The results suggest that this correlation function depends on the separation distance and direction but, not on the positions inside the flame brush. This correlation can be modelled using a combination of Hermite-Gaussian functions of zero and second order, i.e. functions of the form (1-Ax2)e-Bx2 for constants A and B, to include its possible negative values. The integral correlation volume obtained using this model is about 0.2δL3 with the length scale obtained from its cube root being about 0.6δ L, where δ L is the laminar flame thermal thickness. Both of the values are slightly larger than the values reported in an earlier study because of the anisotropy observed for the correlation. This model together with the turbulence-dependent parameter K, the ratio of the root-mean-square (RMS) value of the rate of change of reaction rate to the mean reaction rate, derived from the DNS data is applied to predict the far-field sound emitted from open flames. The calculated noise levels agree well with recently reported measurements and show a sensitivity to K values. © 2012 The Combustion Institute.

Nested sampling for materials: the case of hard spheres

The recently introduced nested sampling algorithm allows the direct and efficient calculation of the partition function of atomistic systems. We demonstrate its applicability to condensed phase systems with periodic boundary conditions by studying the three dimensional hard sphere model. Having obtained the partition function, we show how easy it is to calculate the compressibility and the free energy as functions of the packing fraction and local order, verifying that the transition to crystallinity has a very small barrier, and that the entropic contribution of jammed states to the free energy is negligible for packing fractions above the phase transition. We quantify the previously proposed schematic phase diagram and estimate the extent of the region of jammed states. We find that within our samples, the maximally random jammed configuration is surprisingly disordered.

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.

On representing chemical environments

We review some recently published methods to represent atomic neighbourhood environments, and analyse their relative merits in terms of their faithfulness and suitability for fitting potential energy surfaces. The crucial properties that such representations (sometimes called descriptors) must have are differentiability with respect to moving the atoms, and invariance to the basic symmetries of physics: rotation, reflection, translation, and permutation of atoms of the same species. We demonstrate that certain widely used descriptors that initially look quite different are specific cases of a general approach, in which a finite set of basis functions with increasing angular wave numbers are used to expand the atomic neighbourhood density function. Using the example system of small clusters, we quantitatively show that this expansion needs to be carried to higher and higher wave numbers as the number of neighbours increases in order to obtain a faithful representation, and that variants of the descriptors converge at very different rates. We also propose an altogether new approach, called Smooth Overlap of Atomic Positions (SOAP), that sidesteps these difficulties by directly defining the similarity between any two neighbourhood environments, and show that it is still closely connected to the invariant descriptors. We test the performance of the various representations by fitting models to the potential energy surface of small silicon clusters and the bulk crystal.

Self-similar breakup of near-inviscid liquids.

The final stages of pinchoff and breakup of dripping droplets of near-inviscid Newtonian fluids are studied experimentally for pure water and ethanol. High-speed imaging and image analysis are used to determine the angle and the minimum neck size of the cone-shaped extrema of the ligaments attached to dripping droplets in the final microseconds before pinchoff. The angle is shown to steadily approach the value of 18.0 ± 0.4°, independently of the initial flow conditions or the type of breakup. The filament thins and necks following a τ(2/3) law in terms of the time remaining until pinchoff, regardless of the initial conditions. The observed behavior confirms theoretical predictions.

Density of photon states in dye-doped chiral nematic liquid crystal cells in the presence of losses and gain.

We calculate the density of photon states (DOS) of the normal modes in dye-doped chiral nematic liquid crystal (LC) cells in the presence of various loss mechanisms. Losses and gain are incorporated into the transmission characteristics through the introduction of a small imaginary part in the dielectric constant perpendicular and along the director, for which we assume no frequency dispersion. Theoretical results are presented on the DOS in the region of the photonic band gap for a range of values of the loss coefficient and different values of the optical anisotropy. The obtained values of the DOS at the photonic band gap edges predict a reversal of the dominant modes in the structure. Our results are found to be in good agreement with the experimentally obtained excitation thresholds in chiral nematic LC lasers. The behavior of the DOS is also discussed for amplifying LC cells providing additional insight to the lasing mechanism of these structures.

Flexoelectric and elastic coefficients of odd and even homologous bimesogens.

It is known that bimesogenic liquid crystals exhibit a marked "odd-even" effect in the flexoelastic ratio (the effective flexoelectric coefficient to the average elastic coefficient), with the ratio being higher for the "odd-spaced" bimesogens (those with an odd number of alkyl groups in the spacer chain) than their neighboring even-spaced counterparts. To determine the contribution of each property to the flexoelastic ratio, we present experimental results on the flexoelectric and elastic coefficients of two homologous nonsymmetric bimesogens which possess odd and even alkyl spacers. Our results show that, although there are differences in the flexoelectric coefficients, there are substantially larger differences in the effective elastic coefficient. Specifically, the odd bimesogen is found to have both a low splay elastic coefficient and a very low bend elastic coefficient which, when combined, results in a significantly lower effective elastic coefficient and consequently a higher flexoelastic ratio.