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.

Micromachining of glassy carbon toolsets for micro embossing applications

Laser micro machining is fast gaining popularity as a method of fabricating micro scale structures. Lasers have been utilised for micro structuring of metals, ceramics and glass composites and with advances in material science, new materials are being developed for micro/nano products used in medical, optical, and chemical industries. Due to its favourable strength to weight ratio and extreme resistance to chemical attack, glassy carbon is a new material that offers many unique properties for micro devices. The laser machining of SIGRADUR® G grade glassy carbon was characterised using a 1065 nm wavelength Ytterbium doped pulsed fiber laser. The laser system has a selection of 25 preset waveforms with optimised peak powers for different pulsing frequencies. The optics provide spot diameter of 40 μm at the focus. The effect of fluence, transverse overlap and pulsing frequency (as waveform) on glassy carbon was investigated. Depth of removal and surface roughness were measured as machining quality indicators. The damage threshold fluence was determined to be 0.29 J/cm2 using a pulsing frequency of 250 kHz and a pulse width of 18 ns (waveform 3). Ablation rates of 17 < V < 300 μm3/pulse were observed within a fluence range of 0.98 < F < 2.98 J/cm2. For the same fluence variation, 0.6 μm to 6.8 μm deep trenches were machined. Trench widths varied from 29 μm at lower fluence to 47 μm at the higher fluence. Square pockets, 1 mm wide, were machined to understand the surface machining or milling. The depth of removal using both waveform 3 and 5 showed positive correlation with fluence, with waveform 5 causing more removal than waveform 3 for the same fluence. Machined depths varied from less than 1 μm to nearly 40 μm. For transverse overlap variation using waveform 3, the best surface finish with Rz = 1.1 μm was obtained for fluence 0.792 J/cm2 for transverse overlap of 1 μm, 6 μm, and 9 μm at machined depths of 22.9 μm, 6.6 μm, and 4.6 μm respectively. For fluence of 1.426 J/cm2, the best surface finish with Rz = 1.2 μm was obtained for transverse overlap of 6 μm, and 9 μm at machined depths of 12.46 μm, and 8.6 μm respectively. The experimental data was compiled as machining charts and utilised for fabricating a micro-embossing glassy carbon master toolsets as a capability demonstration.

Performance of magnesia cements in porous blocks in acid and magnesium environments

The performance of porous blocks containing three different reactive magnesia-based cements - namely magnesia alone, magnesium oxide: Portland cement (PC) in 1:1 ratio, cured in ambient conditions, and magnesia alone, cured at elevated carbon dioxide conditions, in hydrochloric acid and magnesium sulfate solution - was investigated. Different aggressive chemical solution conditions were used, to which the samples were exposed for up to 12 months and then tested for strength and microstructure. The performance was also compared with that of standard PC-based blocks. The results showed the significant resistance to chemical attack offered by magnesia, both alone and with PC blend in the porous blocks when cured under ambient carbon dioxide conditions, and confirmed the much poorer performance of blocks made from PC alone. The blocks of solely magnesia cured in elevated carbon dioxide conditions, at 20% concentration, showed slightly lower resistance to acid attack than PC; however, the resistance to sulfate attack was much higher. © 2012 Thomas Telford Ltd.

Tsunami wave runup on the clustered islands using boussinesq-type model

This paper describes the implementation of the Boussinesq-type model and extends its application to the tsunami wave runup on the clustered islands (multiple adjacent conical islands), in turn, an extensively validated two-dimensional Boussinesq-type model is employed to examine the interaction between a propagating solitary wave and multiple idealised conical islands, with particular emphasis on a combination effect of two adjustable parameters for spacing interval/diameter ratio between the adjacent conical islands, S/D, and the rotating angle of the structural configuration,θ on maximum soliton runup heights. An extensive parameter study concerning the combination effect of alteringθ and S/D on the maximum soliton runup with the multi-conical islands is subsequently carried out and the distributions of the maximum runup heights on each conical island are obtained and compared for the twin-island cases. The worst case study is performed for each case in respect of the enhancement in the maximum wave runup heights by the multi-conical islands. It is found that the nonlinear wave diffraction, reflection and refraction play a significant role in varying the maximum soliton runup heights on multiconical islands. The comparatively large maximum soliton runups are generally predicted for the merged and bottom mounted clusteredislands. Furthermore, the joints of the clustered-merged islands are demonstrated to suffer the most of the tsunami wave attack. The conical islands that position in the shadow regions behind the surrounding islands are found to withstand relatively less extreme wave impact. Although, these numerical investigations are considerable simplifications of the multi conical islands, they give a critical insight into certain important hydrodynamic characteristics of the interaction between an extreme wave event and a group of clustered conical islands, and thus providing a useful engineering guidance for extreme wave mitigation and coastal development. Copyright © 2012 by the International Society of Offshore and Polar Engineers (ISOPE).

The unidirectional emptying box

A theoretical description of the turbulent mixing within and the draining of a dense fluid layer from a box connected to a uniform density, quiescent environment through openings in the top and the base of the box is presented in this paper. This is an extension of the draining model developed by Linden et al. (Annu. Rev. Fluid Mech. vol. 31, 1990, pp. 201-238) and includes terms that describe localized mixing within the emptying box at the density interface. Mixing is induced by a turbulent flow of replacement fluid into the box and as a consequence we predict, and observe in complementary experiments, the development of a three-layer stratification. Based on the data collated from previous researchers, three distinct formulations for entrainment fluxes across density interfaces are used to account for this localized mixing. The model was then solved numerically for the three mixing formulations. Analytical solutions were developed for one formulation directly and for a second on assuming that localized mixing is relatively weak though still significant in redistributing buoyancy on the timescale of the draining process. Comparisons between our theoretical predictions and the experimental data, which we have collected on the developing layer depths and their densities show good agreement. The differences in predictions between the three mixing formulations suggest that the normalized flux turbulently entrained across a density interface tends to a constant value for large values of a Froude number FrT, based on conditions of the inflow through the top of the box, and scales as the cube of FrT for small values of FrT. The upper limit on the rate of entrainment into the mixed layer results in a minimum time (tD) to remove the original dense layer. Using our analytical solutions, we bound this time and show that 0.2tE ≈tD tE, i.e. the original dense layer may be depleted up to five times more rapidly than when there is no internal mixing and the box empties in a time tE. © 2010 Cambridge University Press.

The design of multi-element airfoils through multi-objective optimization techniques

This paper presents the development and the application of a multi-objective optimization framework for the design of two-dimensional multi-element high-lift airfoils. An innovative and efficient optimization algorithm, namely Multi-Objective Tabu Search (MOTS), has been selected as core of the framework. The flow-field around the multi-element configuration is simulated using the commercial computational fluid dynamics (cfd) suite Ansys cfx. Elements shape and deployment settings have been considered as design variables in the optimization of the Garteur A310 airfoil, as presented here. A validation and verification process of the cfd simulation for the Garteur airfoil is performed using available wind tunnel data. Two design examples are presented in this study: a single-point optimization aiming at concurrently increasing the lift and drag performance of the test case at a fixed angle of attack and a multi-point optimization. The latter aims at introducing operational robustness and off-design performance into the design process. Finally, the performance of the MOTS algorithm is assessed by comparison with the leading NSGA-II (Non-dominated Sorting Genetic Algorithm) optimization strategy. An equivalent framework developed by the authors within the industrial sponsor environment is used for the comparison. To eliminate cfd solver dependencies three optimum solutions from the Pareto optimal set have been cross-validated. As a result of this study MOTS has been demonstrated to be an efficient and effective algorithm for aerodynamic optimizations. Copyright © 2012 Tech Science Press.

LES of heat transfer in electronics

Numerical methods based on the Reynolds Averaged Navier-Stokes (RANS) and Large Eddy Simulation (LES) equations are applied to the thermal prediction of flows representative of those found in and around electronics systems and components. Low Reynolds number flows through a heated ribbed channel, around a heated cube and within a complex electronics system case are investigated using linear and nonlinear LES models, hybrid RANS-LES and RANS-Numerical-LES (RANS-NLES) methods. Flow and heat transfer predictions using these techniques are in good agreement with each other and experimental data for a range of grid resolutions. Using second order central differences, the RANS-NLES method performs well for all simulations. © 2011 Elsevier Inc.

Lift and the leading edge vortex

Leading edge vortices are considered to be important in generating the high lift coefficients observed in insect flight and may therefore be relevant to micro-air vehicles. A potential flow model of an impulsively started flat plate, featuring a leading edge vortex (LEV) and a trailing edge vortex (TEV) is fitted to experimental data in order to provide insight into the mechanisms that influence the convection of the LEV and to study how the LEV contributes to lift. The potential flow model fits the experimental data best with no bound circulation, which is in accordance with Kelvin's circulation theorem. The lift-to-drag ratio is well approximated by the function 'cot α' for α > 15°, which supports the tentative conclusion that shortly after an impulsive start, at post-stall angles of attack, lift is caused non-circulatory forces and by the action of the LEV as opposed to bound circulation. Copyright © 2012 by C. W. Pitt Ford.

Stiffness and strength of tridimensional periodic lattices

This paper presents a method for the linear analysis of the stiffness and strength of open and closed cell lattices with arbitrary topology. The method hinges on a multiscale approach that separates the analysis of the lattice in two scales. At the macroscopic level, the lattice is considered as a uniform material; at the microscopic scale, on the other hand, the cell microstructure is modelled in detail by means of an in-house finite element solver. The method allows determine the macroscopic stiffness, the internal forces in the edges and walls of the lattice, as well as the global periodic buckling loads, along with their buckling modes. Four cube-based lattices and nine cell topologies derived by Archimedean polyhedra are studied. Several of them are characterized here for the first time with a particular attention on the role that the cell wall plays on the stiffness and strength properties. The method, automated in a computational routine, has been used to develop material property charts that help to gain insight into the performance of the lattices under investigation. © 2012 Elsevier B.V.

Lift and the leading-edge vortex

Flapping wings often feature a leading-edge vortex (LEV) that is thought to enhance the lift generated by the wing. Here the lift on a wing featuring a leading-edge vortex is considered by performing experiments on a translating flat-plate aerofoil that is accelerated from rest in a water towing tank at a fixed angle of attack of 15°. The unsteady flow is investigated with dye flow visualization, particle image velocimetry (PIV) and force measurements. Leading-and trailing-edge vortex circulation and position are calculated directly from the velocity vectors obtained using PIV. In order to determine the most appropriate value of bound circulation, a two-dimensional potential flow model is employed and flow fields are calculated for a range of values of bound circulation. In this way, the value of bound circulation is selected to give the best fit between the experimental velocity field and the potential flow field. Early in the trajectory, the value of bound circulation calculated using this potential flow method is in accordance with Kelvin's circulation theorem, but differs from the values predicted by Wagner's growth of bound circulation and the Kutta condition. Later the Kutta condition is established but the bound circulation remains small; most of the circulation is contained instead in the LEVs. The growth of wake circulation can be approximated by Wagner's circulation curve. Superimposing the non-circulatory lift, approximated from the potential flow model, and Wagner's lift curve gives a first-order approximation of the measured lift. Lift is generated by inertial effects and the slow buildup of circulation, which is contained in shed vortices rather than bound circulation. © 2013 Cambridge University Press.

On the effects of turbulence modelling in design optimisation for high-lift devices

Aerodynamic shape optimisation is being increasingly utilised as a design tool in the aerospace industry. In order to provide accurate results, design optimisation methods rely on the accuracy of the underlying CFD methods applied to obtain aerodynamic forces for a given configuration. Previous studies of the authors have highlighted that the variation of the order of accuracy of the CFD solver with a fixed turbulence model affects the resulting optimised airfoil shape for a single element airfoil. The accuracy of the underlying CFD model is even more relevant in the context of high-lift configurations where an accurate prediction of flow is challenging due to the complex flow physics involving transition and flow separation phenomena. This paper explores the effect of the fidelity of CFD results for a range of turbulence models within the context of the computational design of aircraft configurations. The NLR7301 multi-element airfoil (main wing and flap) is selected as the baseline configuration, because of the wealth of experimental an computational results available for this configuration. An initial validation study is conducted in order to establish optimal mesh parameters. A bi-objective shape optimisation problem is then formulated, by trying to reveal the trade-off between lift and drag coefficients at high angles of attack. Optimisation of the airfoil shape is performed with Spalart-Allmaras, k - ω SST and k - o realisable models. The results indicate that there is consistent and complementary impact to the optimum level achieved from all the three different turbulence models considered in the presented case study. Without identifying particular superiority of any of the turbu- lence models, we can say though that each of them expressed favourable influence towards different optimality routes. These observations lead to the exploration of new avenues for future research. © 2012 AIAA.

On the effects of turbulence modelling in design optimisation for high-lift devices

Aerodynamic shape optimisation is being increasingly utilised as a design tool in the aerospace industry. In order to provide accurate results, design optimisation methods rely on the accuracy of the underlying CFD methods applied to obtain aerodynamic forces for a given configuration. Previous studies of the authors have highlighted that the variation of the order of accuracy of the CFD solver with a fixed turbulence model affects the resulting optimised airfoil shape for a single element airfoil. The accuracy of the underlying CFD model is even more relevant in the context of high-lift configurations where an accurate prediction of flow is challenging due to the complex flow physics involving transition and flow separation phenomena. This paper explores the effect of the fidelity of CFD results for a range of turbulence models within the context of the computational design of aircraft configurations. The NLR7301 multi-element airfoil (main wing and flap) is selected as the baseline configuration, because of the wealth of experimental an computational results available for this configuration. An initial validation study is conducted in order to establish optimal mesh parameters. A bi-objective shape optimisation problem is then formulated, by trying to reveal the trade-off between lift and drag coefficients at high angles of attack. Optimisation of the airfoil shape is performed with Spalart-Allmaras, k - ω SST and k - ε realisable models. The results indicate that there is consistent and complementary impact to the optimum level achieved from all the three different turbulence models considered in the presented case study. Without identifying particular superiority of any of the turbu- lence models, we can say though that each of them expressed favourable influence towards different optimality routes. These observations lead to the exploration of new avenues for future research. © 2012 by the authors.

Detection reliability analysis of passive DAS RFID system

Recently, it has been shown that improved wireless communication coverage can be achieved by employing distributed antenna system (DAS). The DAS RFID system is based on a novel technique whereby two or more spatially separated transmit and receive antennas are used to enable greatly enhanced tag detection performance over longer distances using antenna diversity combined with frequency and phase hopping. In this paper, we present a detection reliability evaluation of the DAS RFID in a typical lab environment. We conduct an extensive experimental analysis of passive RFID tag detection with different locations and orientations. The tag received signal strengths corresponding to various tag locations on one of the six different sides of a cube, and for different reader transmit power are collected and analyzed in this study.

On high-frequency noise scattering by aerofoils in flow

Abstract A theoretical model is developed for the sound scattered when a sound wave is incident on a cambered aerofoil at non-zero angle of attack. The model is based on the linearization of the Euler equations about a steady subsonic flow, and is an adaptation of previous work which considered incident vortical disturbances. Only high-frequency sound waves are considered. The aerofoil thickness, camber and angle of attack are restricted such that the steady flow past the aerofoil is a small perturbation to a uniform flow. The singular perturbation analysis identifies asymptotic regions around the aerofoil; local 'inner' regions, which scale on the incident wavelength, at the leading and trailing edges of the aerofoil; Fresnel regions emanating from the leading and trailing edges of the aerofoil due to the coalescence of singularities and points of stationary phase; a wake transition region downstream of the aerofoil leading and trailing edge; and an outer region far from the aerofoil and wake. An acoustic boundary layer on the aerofoil surface and within the transition region accounts for the effects of curvature. The final result is a uniformly-valid solution for the far-field sound; the effects of angle of attack, camber and thickness are investigated. © 2013 Cambridge University Press.

Macro - A tool to support distributed MDO

This paper addresses a new way for handling distributed design know as the Macro concept. It is based round the assumption that future design teams will become more distributed in nature as industry exploits the Internet and other integrated communication and data exchange systems. The paper notes that this concept is part of an attack on the problems associated with the total process of Distribute Multi-Disciplinary design and Optimisation. The concepts rely on the creation of distributed self-building and self-organising teams made up from members who are globally distributed. The paper describes both the approach adopted and its implementation in a prototype software system operating over the Internet. In essence the work presented is describing a novel method for implementing a distributed design process which is far from complete but which is producing challenging ideas. © 2000 by Cranfield University.