Influence of insoluble surfactant on the deformation and breakup of a bubble or thread in a viscous fluid

The influence of surfactant on the breakup of a prestretched bubble in a quiescent viscous surrounding is studied by a combination of direct numerical simulation and the solution of a long-wave asymptotic model. The direct numerical simulations describe the evolution toward breakup of an inviscid bubble, while the effects of small but non-zero interior viscosity are readily included in the long-wave model for a fluid thread in the Stokes flow limit. The direct numerical simulations use a specific but realizable and representative initial bubble shape to compare the evolution toward breakup of a clean or surfactant-free bubble and a bubble that is coated with insoluble surfactant. A distinguishing feature of the evolution in the presence of surfactant is the interruption of bubble breakup by formation of a slender quasi-steady thread of the interior fluid. This forms because the decrease in surface area causes a decrease in the surface tension and capillary pressure, until at a small but non-zero radius, equilibrium occurs between the capillary pressure and interior fluid pressure. The long-wave asymptotic model, for a thread with periodic boundary conditions, explains the principal mechanism of the slender thread's formation and confirms, for example, the relatively minor role played by the Marangoni stress. The large-time evolution of the slender thread and the precise location of its breakup are, however, influenced by effects such as the Marangoni stress and surface diffusion of surfactant. © 2008 Cambridge University Press.

Continuous spectrum growth and modal instability in swirling duct flow

We present results on the stability of compressible inviscid swirling flows in an annular duct. Such flows are present in aeroengines, for example in the by-pass duct, and there are also similar flows in many aeroacoustic or aeronautical applications. The linearised Euler equations have a ('critical layer') singularity associated with pure convection of the unsteady disturbance by the mean flow, and we focus our attention on this region of the spectrum. By considering the critical layer singularity, we identify the continuous spectrum of the problem and describe how it contributes to the unsteady field. We find a very generic family of instability modes near to the continuous spectrum, whose eigenvalue wavenumbers form an infinite set and accumulate to a point in the complex plane. We study this accumulation process asymptotically, and find conditions on the flow to support such instabilities. It is also found that the continuous spectrum can cause a new type of instability, leading to algebraic growth with an exponent determined by the mean flow, given in the analysis. The exponent of algebraic growth can be arbitrarily large. Numerical demonstrations of the continuous spectrum instability, and also the modal instabilities are presented.

Analytically-based approach to rotor-stator interaction noise in mean swirling flow

A method for modelling and predicting the noise generated by the interaction between the unsteady wake shed from the rotor and a downstream row of stators in a modern ultra-high bypass ducted turbofan engine is described. An analytically-based model is developed to account for three main features of the problem. First, the way in which a typical unsteady wake disturbance from the rotor interacts and is distorted by the mean swirling flow as it propagates downstream. The analysis allows for the inclusion of mean entropy gradients and entropy perturbations. Second, the effects of real stator-blade geometry and proper representation of the genuinely three-dimensional nature of the problem. Third, to model the propagation of the resulting noise back upstream in mean swirling flow. The analytical nature of the problem allows for the inclusion of all wake harmonics and enables the response at all blade passing frequencies to be determined. Example results are presented for an initial wake distribution corresponding to a genuine rotor configuration. Comparisons between numerical data and the asymptotic model for the wake evolution are made. Copyright © 2004 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Curvature and wrinkling of premixed flame Kernels comparison between planar laser induced fluorescence (PLIF) and 3D direct numerical simulations (DNS)

A direct comparison between time resolved PLIF measurements of OH and two dimensional slices from a full three dimensional DNS data set of turbulent premixed flame kernels in lean methane/air mixture was presented. The local flame structure and the degree of flame wrinkling were examined in response to differing turbulence intensities and turbulent Reynolds numbers. Simulations were performed using the SEGA DNS code, which is based on the solution of the compressible Navier Stokes, species, and energy equations for a lean hydrocarbon mixture. For the OH PLIF measurements, a cluster of four Nd:YAG laser was fired sequentially at high repetition rates and used to pump a dye laser. The frequency doubled laser beam was formed into a sheet of 40 mm height using a cylindrical telescope. The combination of PLIF and DNS has been demonstrated as a powerful tool for flame analysis. This research will form the basis for the development of sub-grid-scale (SGS) models for LES of lean-premixed combustion systems such as gas turbines. This is an abstract of a paper presented at the 30th International Symposium on Combustion (Chicago, IL 7/25-30/2004).

Robust multiplexed model predictive control

This paper extends the recently developed multiplexed model predictive control (MMPC) concept to ensure satisfaction of hard constraints despite the action of persistent, unknown but bounded disturbances. MMPC uses asynchronous control moves on each input channel instead of synchronised moves on all channels. It offers reduced computation, by dividing the online optimisation into a smaller problem for each channel, and potential performance improvements, as the response to a disturbance is quicker, albeit via only one channel. Robustness to disturbances is introduced using the constraint tightening approach, tailored to suit the asynchronous updates of MMPC and the resulting time-varying optimisations. Numerical results are presented, involving a simple mechanical example and an aircraft control example, showing the potential computational and performance benefits of the new robust MMPC.