Kelvin-Helmholtz-like instabilities in turbulent flows over complex surfaces

Riblets are small surface protrusions aligned with the flow direction, which confer an anisotropic roughness to the surface [6]. We have recently reported that the transitional-roughness effect in riblets, which limits their performance, is due to a Kelvin–Helmholtz-like instability of the overlying mean flow [7]. According to our DNSs, the instability sets on as the Reynolds number based on the roughness size of the riblets increases, and coherent, elongated spanwise vortices begin to develop immediately above the riblet tips, causing the degradation of the drag-reduction effect. This is a very novel concept, since prior studies had proposed that the degradation was due to the interaction of riblets with the flow as independent units, either to the lodging of quasi-streamwise vortices in the surface grooves [2] or to the shedding of secondary streamwise vorticity at the riblet peaks [9]. We have proposed an approximate inviscid analysis for the instability, in which the presence of riblets is modelled through an average boundary condition for an overlying, spanwise-independent mean flow. This simplification lacks the accuracy of an exact analysis [4], but in turn applies to riblet surfaces in general. Our analysis succeeds in predicting the riblet size for the onset of the instability, while qualitatively reproducing the wavelengths and shapes of the spanwise structures observed in the DNSs. The analysis also connects the observations with the Kelvin–Helmholtz instability of mixing layers. The fundamental riblet length scale for the onset of the instability is a ‘penetration length,’ which reflects how easily the perturbation flow moves through the riblet grooves. This result is in excellent agreement with the available experimental evidence, and has enabled the identification of the key geometric parameters to delay the breakdown. Although the appearance of elongated spanwise vortices was unexpected in the case of riblets, similar phenomena had already been observed over other rough [3], porous [1] and permeable [11] surfaces, as well as over plant [5,14] and urban [12] canopies, both in the transitional and in the fully-rough regimes. However, the theoretical analyses that support the connection of these observations with the Kelvin–Helmholtz instability are somewhat scarce [7, 11, 13]. It has been recently proposed that Kelvin–Helmholtz-like instabilities are a dominant feature common to “obstructed” shear flows [8]. It is interesting that the instability does not require an inflection point to develop, as is often claimed in the literature. The Kelvin-Helmholtz rollers are rather triggered by the apparent wall-normal-transpiration ability of the flow at the plane immediately above the obstructing elements [7,11]. Although both conditions are generally complementary, if wall-normal transpiration is not present the spanwise vortices may not develop, even if an inflection point exists within the roughness [10]. REFERENCES [1] Breugem, W. P., Boersma, B. J. & Uittenbogaard, R. E. 2006 J. Fluid Mech. 562, 35–72. [2] Choi, H., Moin, P. & Kim, J. 1993 J. Fluid Mech. 255, 503–539. [3] Coceal, O., Dobre, A., Thomas, T. G. & Belcher, S. E. 2007 J. Fluid Mech. 589, 375–409. [4] Ehrenstein, U. 2009 Phys. Fluids 8, 3194–3196. [5] Finnigan, J. 2000 Ann. Rev. Fluid Mech. 32, 519–571. [6] Garcia-Mayoral, R. & Jimenez, J. 2011 Phil. Trans. R. Soc. A 369, 1412–1427. [7] Garcia-Mayoral, R. & Jimenez, J. 2011 J. Fluid Mech. doi: 10.1017/jfm.2011.114. [8] Ghisalberti, M. 2009 J. Fluid Mech. 641, 51–61. [9] Goldstein, D. B. & Tuan, T. C. 1998 J. Fluid Mech. 363, 115–151. [10] Hahn, S., Je, J. & Choi, H. 2002 J. Fluid Mech. 450, 259–285. [11] Jimenez, J., Uhlman, M., Pinelli, A. & G., K. 2001 J. Fluid Mech. 442, 89–117. [12] Letzel, M. O., Krane, M. & Raasch, S. 2008 Atmos. Environ. 42, 8770–8784. [13] Py, C., de Langre, E. & Moulia, B. 2006 J. Fluid Mech. 568, 425–449. [14] Raupach, M. R., Finnigan, J. & Brunet, Y. 1996 Boundary-Layer Meteorol. 78, 351–382.

The interaction of riblets with wall-bounded turbulence

The purpose of this thesis is to give answer to the question: why do riblets stop working for a certain size? Riblets are small surface grooves aligned in the mean direction of an overlying turbulent flow, designed specifically to reduce the friction between the flow and the surface. They were inspired by biological surfaces, like the oriented denticles in the skin of fastswimming sharks, and were the focus of a significant amount of research in the late eighties and nineties. Although it was found that the drag reduction depends on the riblet size scaled in wall units, the physical mechanisms implicated have not been completely understood up to now. It has been explained how riblets of vanishing size interact with the turbulent flow, producing a change in the drag proportional to their size, but that is not the regime of practical interest. The optimum performance is achieved for larger sizes, once that linear behavior has broken down, but before riblets begin adopting the character of regular roughness and increasing drag. This regime, which is the most relevant from a technological perspective, was precisely the less understood, so we have focused on it. Our efforts have followed three basic directions. First, we have re-assessed the available experimental data, seeking to identify common characteristics in the optimum regime across the different existing riblet geometries. This study has led to the proposal of a new length scale, the square root of the groove crosssection, to substitute the traditional peak-to-peak spacing. Scaling the riblet dimension with this length, the size of breakdown of the linear behavior becomes roughly universal. This suggests that the onset of the breakdown is related to a certain, fixed value of the cross-section of the groove. Second, we have conducted a set of direct numerical simulations of the turbulent flow over riblets, for sizes spanning the full drag reduction range. We have thus been able to reproduce the gradual transition between the different regimes. The spectral analysis of the flows has proven particularly fruitful, since it has made possible to identify spanwise rollers immediately above the riblets, which begin to appear when the riblet size is close to the optimum. This is a quite surprising feature of the flow, not because of the uniqueness of the phenomenon, which had been reported before for other types of complex and porous surfaces, but because most previous studies had focused on the detail of the flow above each riblet as a unit. Our novel approach has provided the adequate tools to capture coherent structures with an extended spanwise support, which interact with the riblets not individually, but collectively. We have also proven that those spanwise structures are responsible for the increase in drag past the viscous breakdown. Finally, we have analyzed the stability of the flow with a simplified model that connects the appearance of rollers to a Kelvin–Helmholtz-like instability, as is the case also for the flow over plant canopies and porous surfaces. In spite of the model emulating the presence of riblets only in an averaged, general fashion, it succeeds to capture the essential attributes of the breakdown, and provides a theoretical justification for the scaling with the groove cross-section.

Enhanced integrated modelling approach to reconfiguring manufacturing enterprises

Dynamism and uncertainty are real challenges for present day manufacturing enterprises (MEs). Reasons include: an increasing demand for customisation, reduced time to market, shortened product life cycles and globalisation. MEs can reduce competitive pressure by becoming reconfigurable and change-capable. However, modern manufacturing philosophies, including agile and lean, must complement the application of reconfigurable manufacturing paradigms. Choosing and applying the best philosophies and techniques is very difficult as most MEs deploy complex and unique configurations of processes and resource systems, and seek economies of scope and scale in respect of changing and distinctive product flows. It follows that systematic methods of achieving model driven reconfiguration and interoperation of component based manufacturing systems are required to design, engineer and change future MEs. This thesis, titled Enhanced Integrated Modelling Approach to Reconfiguring Manufacturing Enterprises , introduces the development and prototyping a model-driven environment for the design, engineering, optimisation and control of the reconfiguration of MEs with an embedded capability to handle various types of change. The thesis describes a novel systematic approach, namely enhanced integrated modelling approach (EIMA), in which coherent sets of integrated models are created that facilitates the engineering of MEs especially their production planning and control (PPC) systems. The developed environment supports the engineering of common types of strategic, tactical and operational processes found in many MEs. The EIMA is centred on the ISO standardised CIMOSA process modelling approach. Early study led to the development of simulation models during which various CIMOSA shortcomings were observed, especially in its support for aspects of ME dynamism. A need was raised to structure and create semantically enriched models hence forming an enhanced integrated modelling environment. The thesis also presents three industrial case examples: (1) Ford Motor Company; (2) Bradgate Furniture Manufacturing Company; and (3) ACM Bearings Company. In order to understand the system prior to realisation of any PPC strategy, multiple process segments of any target organisation need to be modelled. Coherent multi-perspective case study models are presented that have facilitated process reengineering and associated resource system configuration. Such models have a capability to enable PPC decision making processes in support of the reconfiguration of MEs. During these case studies, capabilities of a number of software tools were exploited such as Arena®, Simul8®, Plant Simulation®, MS Visio®, and MS Excel®. Case study results demonstrated effectiveness of the concepts related to the EIMA. The research has resulted in new contributions to knowledge in terms of new understandings, concepts and methods in following ways: (1) a structured model driven integrated approach to the design, optimisation and control of future reconfiguration of MEs. The EIMA is an enriched and generic process modelling approach with capability to represent both static and dynamic aspects of an ME; and (2) example application cases showing benefits in terms of reduction in lead time, cost and resource load and in terms of improved responsiveness of processes and resource systems with a special focus on PPC; (3) identification and industrial application of a new key performance indicator (KPI) known as P3C the measuring and monitoring of which can aid in enhancing reconfigurability and responsiveness of MEs; and (4) an enriched modelling concept framework (E-MUNE) to capture requirements of static and dynamic aspects of MEs where the conceptual framework has the capability to be extended and modified according to the requirements. The thesis outlines key areas outlining a need for future research into integrated modelling approaches, interoperation and updating mechanisms of partial models in support of the reconfiguration of MEs.

Properties of electrons and holes during femtosecond cooperative emission

Spectral and energetic characteristics of the cooperative recombination of high density electrons and holes in bulk GaAs are experimentally studied at room temperature. It is shown that the properties and parameters of femtosecond superradiant pulses are conditioned by the collective properties of electrons and holes. Electron-hole pairing and the formation of a short-living coherent e-h BCS state distinguish strongly the regime of cooperative emission from all radiative e-h recombination regimes, which have been observed earlier. The dependences of the energy gap (the order parameter), the Fermi energy, and the band gap of the coherent e-h BCS state on the concentration of electron-hole pairs are obtained.

Superfluorescence in semiconductor lasers

An analysis is made of the conditions for the generation of superfluorescence pulses in an inverted medium of electron-hole pairs in a semiconductor. It is shown that strong optical amplification in laser semiconductor amplifiers characterised by αL ≫ 1 leads to suppression of phase re-laxation of the medium during the initial stages of evolution of superfluorescence and to formation of a macroscopic dipole from electron-hole pairs. Cooperative emission of radiation in this system results in generation of a powerful ultrashort pulse of the optical gain, which interacts coherently with the semiconductor medium. It is shown that coherent pulsations of the optical field, observed earlier by the author in Q-switched semiconductor lasers, are the result of superfluorescence and of the coherent interaction between the optical field and the medium.

Superfluorescence in semiconductor lasers

An analysis is made of the conditions for the generation of superfluorescence pulses in an inverted medium of electron-hole pairs in a semiconductor. It is shown that strong optical amplification in laser semiconductor amplifiers characterised by αL ≫ 1 (α is the small-signal gain and L is the amplifier length) leads to suppression of phase relaxation of the medium during the initial stages of evolution of superfluorescence and to formation of a macroscopic dipole from electron - hole pairs. Cooperative emission of radiation in this system results in generation of a powerful ultrashort pulse of the optical gain, which interacts coherently with the semiconductor medium. It is shown that coherent pulsations of the optical field, observed earlier by the author in Q-switched semiconductor lasers, are the result of superfluorescence and of the coherent interaction between the optical field and the medium.

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.

Single sublattice endotaxial phase separation driven by charge frustration in a complex oxide.

Complex transition-metal oxides are important functional materials in areas such as energy and information storage. The cubic ABO3 perovskite is an archetypal example of this class, formed by the occupation of small octahedral B-sites within an AO3 network defined by larger A cations. We show that introduction of chemically mismatched octahedral cations into a cubic perovskite oxide parent phase modifies structure and composition beyond the unit cell length scale on the B sublattice alone. This affords an endotaxial nanocomposite of two cubic perovskite phases with distinct properties. These locally B-site cation-ordered and -disordered phases share a single AO3 network and have enhanced stability against the formation of a competing hexagonal structure over the single-phase parent. Synergic integration of the distinct properties of these phases by the coherent interfaces of the composite produces solid oxide fuel cell cathode performance superior to that expected from the component phases in isolation.

Optimal recall from bounded metaplastic synapses: predicting functional adaptations in hippocampal area CA3.

A venerable history of classical work on autoassociative memory has significantly shaped our understanding of several features of the hippocampus, and most prominently of its CA3 area, in relation to memory storage and retrieval. However, existing theories of hippocampal memory processing ignore a key biological constraint affecting memory storage in neural circuits: the bounded dynamical range of synapses. Recent treatments based on the notion of metaplasticity provide a powerful model for individual bounded synapses; however, their implications for the ability of the hippocampus to retrieve memories well and the dynamics of neurons associated with that retrieval are both unknown. Here, we develop a theoretical framework for memory storage and recall with bounded synapses. We formulate the recall of a previously stored pattern from a noisy recall cue and limited-capacity (and therefore lossy) synapses as a probabilistic inference problem, and derive neural dynamics that implement approximate inference algorithms to solve this problem efficiently. In particular, for binary synapses with metaplastic states, we demonstrate for the first time that memories can be efficiently read out with biologically plausible network dynamics that are completely constrained by the synaptic plasticity rule, and the statistics of the stored patterns and of the recall cue. Our theory organises into a coherent framework a wide range of existing data about the regulation of excitability, feedback inhibition, and network oscillations in area CA3, and makes novel and directly testable predictions that can guide future experiments.

The rocking response of large flexible structures to earthquakes

The rocking response of structures subjected to strong ground motions is a problem of 'several scales'. While small structures are sensitive to acceleration pulses acting successively, large structures are more significantly affected by coherent low frequency components of ground motion. As a result, the rocking response of large structures is more stable and orderly, allowing effective isolation from the ground without imminent danger of overturning. This paper aims to characterize and predict the maximum rocking response of large and flexible structures to earthquakes using an idealized structural model. To achieve this, the maximum rocking demand caused by different earthquake records was evaluated using several ground motion intensity measures. Pulse-type records which typically have high peak ground velocity and lower frequency content caused large rocking amplitudes, whereas non-pulse type records caused random rocking motion confined to small rocking amplitudes. Coherent velocity pulses were therefore identified as the primary cause of significant rocking motion. Using a suite of pulse-type ground motions, it was observed that idealized wavelets fitted to velocity pulses can adequately describe the rocking response of large structures. Further, a parametric analysis demonstrates that pulse shape parameters affect the maximum rocking response significantly. Based on these two findings, a probabilistic analysis method is proposed for estimating the maximum rocking demand to pulse-type earthquakes. The dimensionless demand maps, produced using these methods, have predictive power in the near-field provided that pulse period and amplitude can be estimated a priori. Use of this method within a probabilistic seismic demand analysis framework is briefly discussed. © 2013 Springer Science+Business Media Dordrecht.