Asymptotic solutions of glass temperature profiles during steady optical fibre drawing

In this paper we derive realistic simplified models for the high-speed drawing of glass optical fibres via the downdraw method, that capture the fluid dynamics and heat transport in the fibre via conduction, convection and radiative heating. We exploit the small aspect ratio of the fibre and the relative orders of magnitude of the dimensionless parameters that characterize the heat transfer to reduce the problem to one- or two-dimensional systems via asymptotic analysis. The resulting equations may be readily solved numerically and in many cases admit exact analytic solutions. The systematic asymptotic breakdown presented is used to elucidate the relative importance of furnace temperature profile, convection, surface radiation and conduction in each portion of the furnace and the role of each in controlling the glass temperature.\ud \ud The models derived predict many of the qualitative features observed in the real industrial process, such as the glass temperature profile within the furnace and the sharp transition in fibre thickness. The models thus offer a desirable route to quick scenario testing, providing valuable practical information into the dependencies of the solution on the parameters and the dominant heat-transport mechanism

The dynamics of liquid slugs forced by a syringe pump

Microfluidic processes for chemical synthesis have become popular in recent years. The small scale of the chemical reactions promise greater control over reaction conditions and more timely creation of products. The small scale of microfluidics poses its own set of problems, however. At the microscale, the dominant fluid forces are viscous resistance and surface tension. The effects of viscosity and scale reduce the Reynolds number and make mixing difficult. Much work has been done to control mixing at the microscale. This problem is concerned with a different microfluidic problem: delivering reactants to the site of reaction. A common setup is to attach syringes full of reactant to a reaction chamber by narrow hydrophobic tubing. Using a stepper motor, a controlled dose of liquid may be injected into the tube. The hydrophobosity causes the dose to curve outward on the sides, becoming a "slug" of reactant with air in front and behind. The syringe at the rear is then switched for one full of air, and air pressure is used to drive the slug to the reaction site. If too much pressure is applied, the slug will arrive with a significant back pressure that will be relieved through bubbling in the reaction site. This causes the formation of a foam and is highly undesirable. We present a simple model based on Boyle’s law for the motion of a slug through a tube. We then extend this model for trains of slugs separated by air bubbles. Last, we consider the case of a flooded reaction site, where the forward air bubble must be pushed through the flooding liquid. In conclusion, we have determined the dynamics of a single slug moving towards an empty reaction chamber giving the final equilibrium position of the slug. A phase-plane analysis then determined a condition on the size of the slug needed to ensure that it comes to rest without oscillating about the equilibrium position. The effect of a flooded reaction chamber was then considered. In this case it is impossible to avoid bubbling due to the design of the device. We found that it is possible, however, to reduce the bubbling by minimising the back pressure behind the slug. Finally, the dynamics of multiple slugs with or without a flooded reaction chamber has been investigated

On the predictions and limitations of the BeckerDoring model for reaction kinetics in micellar surfactant solutions

We investigate the breakdown of a system of micellar aggregates in a surfactant solution following an order-one dilution. We derive a mathematical model based on the Becker–Döring system of equations, using realistic expressions for the reaction constants fit to Molecular Dynamics simulations. We exploit the largeness of typical aggregation numbers to derive a continuum model, substituting a large system of ordinary differential equations for a partial differential equation in two independent variables: time and aggregate size. Numerical solutions demonstrate that re-equilibration occurs in two distinct stages over well-separated time-scales, in agreement with experiment and with previous theories. We conclude by exposing a limitation in the Becker–Döring theory for re-equilibration and discuss potential resolutions

An asymptotic theory for the re-equilibration of a micellar surfactant solution

Micellar surfactant solutions are characterized by a distribution of aggregates comprised predominantly of pre-micellar aggregates (monomers, dimers, trimers, etc.) and a region of proper micelles close to the peak aggregation number, connected by an intermediate region containing a very low concentration of aggregates. Such a distribution gives rise to a distinct two-timescale re-equilibration following a system dilution, known as the 1 and 2 processes, whose dynamics may be described by the Becker–Döring equations. We use a continuum version of these equations to develop a reduced asymptotic description that elucidates the behavior during each of these processes

The kinetics of surfactant desorption at the air–solution interface

The kinetics of desorption of the anionic surfactant sodium dodecylbenzene sulfonate at the air–solution interface have been studied using neutron reflectivity (NR). The experimental arrangement incorporates a novel flow cell in which the subphase can be exchanged (diluted) using a laminar flow whilst the surface region remains unaltered. The kinetics of the desorption is relatively slow and occurs over many tens of minutes compared with the dilution timescale of approximately 10–30 minutes. A detailed mathematical model, in which the rate of the desorption is determined by transport through a near-surface diffusion layer into a diluted bulk solution below, is developed and provides a good description of the timedependent adsorption data.\ud \ud A key parameter of the model is the ratio of the depth of the diffusion layer, Hc , to the depth of the fluid, Hf, and we find that this is related to the reduced Péclet number, Pe*, for the system, via Hc/Hf, = C/Pe* 1/ 2 . Although from a highly idealised experimental arrangement, the results provide an important insight into the ‘rinse mechanism’, which is applicable to a wide variety of domestic and industrial circumstances

Estimating specific surface area of fine stream bed sediments from geochemistry

Specific surface area (SSA) of headwater stream bed sediments is a fundamental property which determines the nature of sediment surface reactions and influences ecosystem-level, biological processes. Measurements of SSA – commonly undertaken by BET nitrogen adsorption – are relatively costly in terms of instrumentation and operator time. A novel approach is presented for estimating fine (2.5 mg kg−1), four elements were identified as significant predictors of SSA (ordered by decreasing predictive power): V > Ca > Al > Rb. The optimum model from these four elements accounted for 73% of the variation in bed sediment SSA (range 6–46 m2 g−1) with a root mean squared error of prediction – based on leave-one-out cross-validation – of 6.3 m2 g−1. It is believed that V is the most significant predictor because its concentration is strongly correlated both with the quantity of Fe-oxides and clay minerals in the stream bed sediments, which dominate sediment SSA. Sample heterogeneity in SSA – based on triplicate measurements of sub-samples – was a substantial source of variation (standard error = 2.2 m2 g−1) which cannot be accounted for in the regression model. The model was used to estimate bed sediment SSA at the other 1792 sites and at 30 duplicate sites where an extra sediment sample had been collected, 25 m from the original site. By delineating sub-catchments for the headwater sediment sites only those sub-catchments were selected with a dominant (>50% of the sub-catchment area) bedrock formation and land use type; the bedrock and land use classes accounted for 39% and 7% of the variation in bed sediment SSA, respectively. Variation in estimated, fine bed sediment SSA from the paired, duplicate sediment sites was small (2.7 m2 g−1), showing that local variation in SSA at stream sites is modest when compared to that between catchments. How the approach might be applied in other environments and its potential limitations are discussed

The Mersey Estuary : sediment geochemistry

This report describes a study of the geochemistry of the Mersey estuary carried out between April 2000 and December 2002. The study was the first in a new programme of surveys of the geochemistry of major British estuaries aimed at enhancing our knowledge and understanding of the distribution of contaminants in estuarine sediments. The report first summarises the physical setting, historical development, geology, hydrography and bathymetry of the Mersey estuary and its catchment. Details of the sampling and analytical programmes are then given followed by a discussion of the sedimentology and geochemistry. The chemistry of the water column and suspended particulate matter have not been studied, the chief concern being with the geochemistry of the surface and near-surface sediments of the Mersey estuary and an examination of their likely sources and present state of contamination

The effect of polar lipids on tear film dynamics

In this paper we present a mathematical model describing the effect of polar lipids on the evolution of a precorneal tear film, with the aim of explaining the interesting experimentally observed phenomenon that the tear film continues to move upwards even after the upper eyelid has become stationary. The polar lipid is an insoluble surface species that locally alters the surface tension of the tear film. In the lubrication limit, the model reduces to two coupled nonlinear partial differential equations for the film thickness and the concentration of lipid. We solve the system numerically and observe that the presence of the lipid causes an increase in flow of liquid up the eye. We further exploit the size of the parameters in the problem to explain the initial evolution of the system

Modelling Li+ Ion Battery Electrode Properties

We formulated two detailed models for an electrolytic cell with particulate electrodes based on a lithium atom concentration dependent Butler-Volmer condition at the interface between electrode particles and the electrolyte. The first was based on a dilute-ion assumption for the electrolyte, while the second assumed that Li ions are present in excess. For the first, we used the method of multiple scales to homogenize this model over the microstructure, formed by the small lithium particles in the electrodes. For the second, we gave rigorous bounds for the effective electrochemical conductivity for a linearized case. We expect similar results and bounds for the "full nonlinear problem" because variational results are generally not adversely affected by a sinh term. Finally we used the asymptotic methods, based on parameters estimated from the literature, to attain a greatly simplified one-dimensional version of the original homogenized model. This simplified model accounts for the fact that diffusion of lithium atoms within individual electrode particles is relatively much faster than that of lithium ions across the whole cell so that lithium ion diffusion is what limits the performance of the battery. However, since most of the potential drop occurs across the Debye layers surrounding each electrode particle, lithium ion diffusion only significantly affects cell performance if there is more or less complete depletion of lithium ions in some region of the electrolyte which causes a break in the current flowing across the cell. This causes catastrophic failure. Providing such failure does not occur the potential drop across the cell is determined by the concentration of lithium atoms in the electrode particles. Within each electrode lithium atom concentration is, to leading order, a function of time only and not of position within the electrode. The depletion of electrode lithium atom concentration is directly proportional to the current being drawn off the cell. This leads one to expect that the potential of the cell gradually drops as current is drawn of it. We would like to emphasize that all the homogenization methods employed in this work give a systematic approach for investigating the effect that changes in the microstructure have on the behaviour of the battery. However, due to lack of time, we have not used this method to investigate particular particle geometries