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.

Sound transmission in strongly-curved slowly-varying cylindrical and annular lined ducts with flow

In this paper we consider the propagation of acoustic waves along a curved hollow or annular duct with lined walls. The curvature of the duct centreline and the wall radii vary slowly along the duct, allowing application of an asymptotic multiple scales analysis. This generalises Rienstra's analysis of a straight duct of varying cross-sectional radius. The result of the analysis is that the modal wavenumbers and mode shapes are determined locally as modes of a torus with the same local curvature, while the amplitude of the modes evolves as the mode propagates along the duct. The duct modes are found numerically at each axial location using a pseudo-spectral method. Unlike the case of a straight duct, there is a fundamental asymmetry between upstream and downstream propagating modes, with some mode shapes tending to be concentrated on either the inside or outside of the bend depending on the direction of propagation. The interaction between the presence of wall lining and curvature is investigated in particular; for instance, in a representative case it is found that the curvature causes the first few acoustic modes to be more heavily damped by the duct boundary than would be expected for a straight duct. Analytical progress can be made in the limit of very high mode order, in which case well-known 'whispering gallery' modes, localised close to the wall, can be identified.

Microgripper using a diamond like carbon/metal bilayer

A novel normally closed microcage has been fabricated and characterized. This device was made from a highly compressively stressed diamond like carbon (DLC) and electroplated Ni bimorph structure. The large stress in the DLC causes the bimorph layer to curve once it is released from the substrate. The radius of curvature is in the range of 18 - 50μm, and can be controlled by varying the DLC and the Ni thicknesses. The devices can be operated in a pulsed mode current with low operation temperature, and can be opened by ∼60μm laterally with a power consumption of only ∼16mW. © 2004 IEEE.

Electrical and field emission investigation of individual carbon nanotubes from plasma enhanced chemical vapour deposition

Plasma enhanced chemical vapour deposition (PECVD) is a controlled technique for the production of vertically aligned multiwall carbon nanotubes for field emission applications. In this paper, we investigate the electrical properties of individual carbon nanotubes which is important for designing field emission devices. PECVD nanotubes exhibit a room temperature resistance of 1-10 kΩ/μm length (resistivity 10-6 to 10-5 Ω m) and have a maximum current carrying capability of 0.2-2 mA (current density 107-108 A/cm2). The field emission characteristics show that the field enhancement of the structures is strongly related to the geometry (height/radius) of the structures and maximum emission currents of ∼ 10 μA were obtained. The failure of nanotubes under field emission is also discussed. © 2002 Elsevier Science B.V. All rights reserved.

Evaluating boundary dependent errors in QM/MM simulations.

Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations provide a powerful tool for studying chemical reactions, especially in complex biochemical systems. In most works to date, the quantum region is kept fixed throughout the simulation and is defined in an ad hoc way based on chemical intuition and available computational resources. The simulation errors associated with a given choice of the quantum region are, however, rarely assessed in a systematic manner. Here we study the dependence of two relevant quantities on the QM region size: the force error at the center of the QM region and the free energy of a proton transfer reaction. Taking lysozyme as our model system, we find that in an apolar region the average force error rapidly decreases with increasing QM region size. In contrast, the average force error at the polar active site is considerably higher, exhibits large oscillations and decreases more slowly, and may not fall below acceptable limits even for a quantum region radius of 9.0 A. Although computation of free energies could only be afforded until 6.0 A, results were found to change considerably within these limits. These errors demonstrate that the results of QM/MM calculations are heavily affected by the definition of the QM region (not only its size), and a convergence test is proposed to be a part of setting up QM/MM simulations.

Modelling the entrainment and motion of particles in a gap: Application to abrasive wear

The microscale abrasion or ball-cratering test is being increasingly applied to a wide range of bulk materials and coatings. The response of materials to this test depends critically on the nature of the motion of the abrasive particles in the contact zone: whether they roll and produce multiple indentations in the coating, or slide causing grooving abrasion. Similar phenomena also occur when hard contaminant particles enter a lubricated contact. This paper presents simple quantitative two-dimensional models which describe two aspects of the interaction between a hard abrasive particle and two sliding surfaces. The first model treats the conditions under which a spherical abrasive particle of size d can be entrained into the gap between a rotating sphere of radius R and a plane surface. These conditions are determined by the coefficients of friction between the particle and the sphere, and the particle and the plane, denoted by μs and μp respectively. This model predicts that the values of (μs + μp) and 2μs should both exceed √2d/R for the particles to be entrained into the contact. If either is less than this value, the particle will slide against the sphere and never enter the contact. The second model describes the mechanisms of abrasive wear in a contact when an idealized rhombus-sectioned prismatic particle is located between two parallel plane surfaces separated by a certain distance, which can represent either the thickness of a fluid film or the spacing due to the presence of other particles. It is shown that both the ratio of particle size to the separation of the surfaces and the ratio of the hardnesses of the two surfaces have important influences on the particle motion and hence on the mechanism of the resulting abrasive wear. Results from this model are compared with experimental observations, and the model is shown to lead to realistic predictions. © IMechE 2003.

Hydrodynamic entrapment of bacteria swimming near a solid surface.

The near-surface motility of bacteria is important in the initial formation of biofilms and in many biomedical applications. The swimming motion of Escherichia coli near a solid surface is investigated both numerically and experimentally. A boundary element method is used to predict the hydrodynamic entrapment of E. coli bacteria, their trajectories, and the minimum separation of the cell from the surface. The numerical results show the existence of a stable swimming distance from the boundary that depends only on the shape of the cell body and the flagellum. The experimental validation of the numerical approach allows one to use the numerical method as a predictive tool to estimate with reasonable accuracy the near-wall motility of swimming bacteria of known geometry. The analysis of the numerical database demonstrated the existence of a correlation between the radius of curvature of the near-wall circular trajectory and the separation gap. Such correlation allows an indirect estimation of either of the two quantities by a direct measure of the other without prior knowledge of the cell geometry. This result may prove extremely important in those biomedical and technical applications in which the near-wall behavior of bacteria is of fundamental importance.

DEVELOPMENT OF MULTI-ANCHOR EARTH RETENTION SYSTEMS.

This paper discusses two projects, the first relating to tests on 'reinforced earth' conducted on the 1. 5 m radius centrifuge at U. M. I. S. T. and the second concerning a short pilot study on the T. R. R. L concept of 'anchored earth' carried out on the 5 m radius Cambridge Geotechnical Centrifuge. The paper proposes a role for centrifugal model testing in design evaluation, whether by the designer himself, his centrification authority, a rival patent holder, or a research worker.

The effect of an upstream turbine on a low-aspect ratio vane

The interaction between a high-pressure rotor and a downstream vane is dominated by vortex-blade interaction. Each rotor blade passing period two co-rotating vortex pairs, the tip-leakage and upper passage vortex and the lower passage and trailing shed vortex, impinge on, and are cut by, the vane leading edge. In addition to the streamwise vortex the tip-leakage flow also contains a large velocity deficit. This causes the interaction of the tip-leakage flow with a downstream vane to differ from typical vortex blade interaction. This paper investigates the effect these interaction mechanisms have on a downstream vane. The test geometry considered was a low aspect ratio second stage vane located within a S-shaped diffuser with large radius change mounted downstream of a shroudless high-pressure turbine stage. Experimental measurements were conducted at engine-representative Mach and Reynolds numbers, and data was acquired using a fast-response aerodynamic probe upstream and downstream of the vane. Time-resolved numerical simulations were undertaken with and without a rotor tip gap in order to investigate the relative magnitude of the interaction mechanisms. The presence of the upstream stage is shown to significantly change the structure of the secondary flow in the vane and to cause a small drop in its performance.

Vortex filament simulation of trailing vortex merger

Aircraft in high-lift configuration shed multiple vortices. These generally merge to form a downstream wake consisting of two counter-rotating vortices of equal strength. The understanding of the merger of two co-rotating trailing vortices is important in evaluating the separation criteria for different aircraft to prevent wake vortex hazards during landing and take-off. There is no existing theoretical method on the basis of which such norms can be set. The present study is aimed at gaining a better understanding of the behaviour of wake vortices behind the aircraft. Two dimensional studies are carried out using the vortex blob method and compared with Bertenyi's experiment. It is shown that inviscid two dimensional effects are insufficient to explain the observations. Three dimensional studies, using the vortex filament method, are applied to the same test case. Two Lamb-Oseen profile vortices of the same dimensions and initial separation as the experiment are allowed to evolve from a straight starting condition until a converged steady flow is achieved. The results obtained show good agreement with the experimental distance to merger. Core radius and separation behaviour is qualitatively similar to experiment, with the exception of rapid increases at first. This may be partially attributable to the choice of filament-element length, and recommended further work includes a convergence study for this parameter. Copyright © 2005 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Friction of sliding surfaces carrying adsorbed lubricant layers

In concentrated contacts the behaviour of lubricants is much modified by the high local pressures: changes can arise both from molecular ordering within the very thin film lubricant layers present at the interface as well as from the deposition on the component surfaces of more solid-like polymeric boundary layers. These 'third bodies' separating the solid surfaces may have rheological or mechanical properties very different from those observed in the bulk. Classical elasto-hydrodynamic theory considers the entrapped lubricant to exhibit a piezo-viscous behaviour while the conventional picture of more solid boundary lubricant layers views their shear strength r as being linearly dependent on local pressure p, so that T = TO + ap where TO and a are constants. If TO is relatively small, then the coefficient of friction \i = T Ip ~ a and so Amonton's laws are recovered. However, the properties of adsorbed or deposited surface films, or indeed other third bodies such as debris layers, may be more complex than this. A preliminary study has looked quantitatively at the influence of the pressure dependence of the shear strength of any surface layer on the overall friction coefficient of a contact which is made up of an array of asperities whose height varies in a Gaussian manner. Individual contact points may be elastic or plastic. The analysis results in plots of coefficient of friction versus the service or load parameter PIH&NRa where P is the nominal pressure on the contact, HS the hardness of the deforming surface, N the asperity density, R the mean radius of curvature of the asperities, and a is the standard deviation of their height distribution. In principle, any variation oft withp can be incorporated into the model; however, in this initial study we have used data on colloidal suspensions from the group at the Ecole Centrale de Lyon as well as examining the effect of functional relationships of somewhat greater complexity than a simple linear form. Results of the analysis indicate that variations in fj. are possible as the load is varied which depend on the statistical spread of behaviour at individual asperity contacts. The value of this analysis is that it attempts to combine the behaviour of films on the molecular scale with the topography of real engineering surfaces and so give an indication of the effects at the full-size or macro-scale that can be achieved by chemical or molecular surface engineering.

The elastic-plastic indentation response of a columnar thermal barrier coating

Thermal barrier coatings with a columnar microstructure are prone to erosion damage by a mechanism of surface cracking upon impact by small foreign particles. In order to explore this erosion mechanism, the elastic indentation and the elastic-plastic indentation responses of a columnar thermal barrier coating to a spherical indenter were determined by the finite element method and by analytical models. It was shown that the indentation response is intermediate between that of a homogeneous half-space and that given by an elastic-plastic mattress model (with the columns behaving as independent non-linear springs). The sensitivity of the indentation behaviour to geometry and to the material parameters was explored: the diameter of the columns, the gap width between columns, the coefficient of Coulomb friction between columns and the layer height of the thermal barrier coating. The calculations revealed that the level of induced tensile stress is sufficient to lead to cracking of the columns at a depth of about the column radius. It was also demonstrated that the underlying soft bond coat can undergo plastic indentation when the coating comprises parallel columns, but this is less likely for the more realistic case of a random arrangement of tapered columns. © 2009 Elsevier B.V.

Estimating material elasticity by spherical indentation load-relaxation tests on viscoelastic samples of finite thickness

A two-step viscoelastic spherical indentation method is proposed to compensate for 1) material relaxation and 2) sample thickness. In the first step, the indenter is moved at a constant speed and the reaction force is measured. In the second step, the indenter is held at a constant position and the relaxation response of the material is measured. Then the relaxation response is fit with a multi-exponential function which corresponds to a three-branch general Maxwell model. The relaxation modulus is derived by correcting the finite ramp time introduced in the first step. The proposed model takes into account the sample thickness, which is important for applications in which the sample thickness is less than ten times the indenter radius. The model is validated numerically by finite element simulations. Experiments are carried out on a 10% gelatin phantom and a chicken breast sample with the proposed method. The results for both the gelatin phantom and the chicken breast sample agree with the results obtained from a surface wave method. Both the finite element simulations and experimental results show improved elasticity estimations by incorporating the sample thickness into the model. The measured shear elasticities of the 10% gelatin sample are 6.79 and 6.93 kPa by the proposed finite indentation method at sample thickness of 40 and 20 mm, respectively. The elasticity of the same sample is estimated to be 6.53 kPa by the surface wave method. For the chicken breast sample, the shear elasticity is measured to be 4.51 and 5.17 kPa by the proposed indentation method at sample thickness of 40 and 20 mm, respectively. Its elasticity is measured by the surface wave method to be 4.14 kPa. © 2011 IEEE.

Plasmonic band gaps and waveguide effects in carbon nanotube arrays based metamaterials.

Highly dense periodic arrays of multiwalled carbon nanotubes behave like low-density plasma of very heavy charged particles, acting as metamaterials. These arrays with nanoscale lattice constants can be designed to display extended plasmonic band gaps within the optical regime, encompassing the crucial optical windows (850 and 1550 nm) simultaneously. We demonstrate an interesting metamaterial waveguide effect displayed by these nanotube arrays containing line defects. The nanotube arrays with lattice constants of 400 nm and radius of 50 nm were studied. Reflection experiments conducted on the nanoscale structures were in agreement with numerical calculations.