Shrinking microbubbles with microfluidics: mathematical modelling to control microbubble sizes.

Microbubbles have applications in industry and life-sciences. In medicine, small encapsulated bubbles (<10 µm) are desirable because of their utility in drug/oxygen delivery, sonoporation, and ultrasound diagnostics. While there are various techniques for generating microbubbles, microfluidic methods are distinguished due to their precise control and ease-of-fabrication. Nevertheless, sub-10 µm diameter bubble generation using microfluidics remains challenging, and typically requires expensive equipment and cumbersome setups. Recently, our group reported a microfluidic platform that shrinks microbubbles to sub-10 µm diameters. The microfluidic platform utilizes a simple microbubble-generating flow-focusing geometry, integrated with a vacuum shrinkage system, to achieve microbubble sizes that are desirable in medicine, and pave the way to eventual clinical uptake of microfluidically generated microbubbles. A theoretical framework is now needed to relate the size of the microbubbles produced and the system's input parameters. In this manuscript, we characterize microbubbles made with various lipid concentrations flowing in solutions that have different interfacial tensions, and monitor the changes in bubble size along the microfluidic channel under various vacuum pressures. We use the physics governing the shrinkage mechanism to develop a mathematical model that predicts the resulting bubble sizes and elucidates the dominant parameters controlling bubble sizes. The model shows a good agreement with the experimental data, predicting the resulting microbubble sizes under different experimental input conditions. We anticipate that the model will find utility in enabling users of the microfluidic platform to engineer bubbles of specific sizes.

Kinetics of surfactant desorption at an air-solution interface.

The kinetics of re-equilibration of the anionic surfactant sodium dodecylbenzene sulfonate at the air-solution interface have been studied using neutron reflectivity. The experimental arrangement incorporates a novel flow cell in which the subphase can be exchanged (diluted) using a laminar flow while the surface region remains unaltered. The rate of the re-equilibration is relatively slow and occurs over many tens of minutes, which is comparable with the dilution time scale of approximately 10-30 min. 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 time-dependent adsorption data. A key parameter of the model is the ratio of the depth of the diffusion layer, H(c), to the depth of the fluid, H(f), and we find that this is related to the reduced Péclet number, Pe*, for the system, via H(c)/H(f) = C/Pe*(1/2). Although from a highly idealized experimental arrangement, the results provide an important insight into the "rinse mechanism", which is applicable to a wide variety of domestic and industrial circumstances.

A comparison of the legal frameworks supporting water management in Europe and China.

This paper has compared the legal frameworks supporting water management in Europe and China, with special focus on integrated river basin management (IRBM) to identify synergies and opportunities in policymaking and implementation. The research shows that China has committed to the efficient management of water resources through various policy tools during the current period. This commitment, however, has often been interrupted and distorted by politics, resulting in the neglect of socioeconomic and environmental priorities. The European legal framework supporting water management underwent a complex and lengthy development, but with the adoption of the Water Framework Directive provides a policy model on which to develop an integrated and sustainable approach to river basin management, elements of which may help to meet the demands of the emerging 21st century Chinese society on these critical natural resources.

Curvature- and fluid-stress-driven tissue growth in a tissue-engineering scaffold pore.

Cell proliferation within a fluid-filled porous tissue-engineering scaffold depends on a sensitive choice of pore geometry and flow rates: regions of high curvature encourage cell proliferation, while a critical flow rate is required to promote growth for certain cell types. When the flow rate is too slow, the nutrient supply is limited; when it is too fast, cells may be damaged by the high fluid shear stress. As a result, determining appropriate tissue-engineering-construct geometries and operating regimes poses a significant challenge that cannot be addressed by experimentation alone. In this paper, we present a mathematical theory for the fluid flow within a pore of a tissue-engineering scaffold, which is coupled to the growth of cells on the pore walls. We exploit the slenderness of a pore that is typical in such a scenario, to derive a reduced model that enables a comprehensive analysis of the system to be performed. We derive analytical solutions in a particular case of a nearly piecewise constant growth law and compare these with numerical solutions of the reduced model. Qualitative comparisons of tissue morphologies predicted by our model, with those observed experimentally, are also made. We demonstrate how the simplified system may be used to make predictions on the design of a tissue-engineering scaffold and the appropriate operating regime that ensures a desired level of tissue growth.

Stochastic modelling of membrane filtration.

Membrane fouling during particle filtration occurs through a variety of mechanisms, including internal pore clogging by contaminants, coverage of pore entrances and deposition on the membrane surface. In this paper, we present an efficient method for modelling the behaviour of a filter, which accounts for different retention mechanisms, particle sizes and membrane geometries. The membrane is assumed to be composed of a series of, possibly interconnected, pores. The central feature is a conductivity function, which describes the blockage of each individual pore as particles arrive, which is coupled with a mechanism to account for the stochastic nature of the arrival times of particles at the pore. The result is a system of ordinary differential equations based on the pore-level interactions. We demonstrate how our model can accurately describe a wide range of filtration scenarios. Specifically, we consider a case where blocking via multiple mechanisms can occur simultaneously, which have previously required the study through individual models; the filtration of a combination of small and large particles by a track-etched membrane and particle separation using interconnected pore networks. The model is significantly faster than comparable stochastic simulations for small networks, enabling its use as a tool for efficient future simulations.

Dissecting the self-assembly kinetics of multimeric pore-forming toxins.

Pore-forming toxins are ubiquitous cytotoxins that are exploited by both bacteria and the immune response of eukaryotes. These toxins kill cells by assembling large multimeric pores on the cell membrane. However, a quantitative understanding of the mechanism and kinetics of this self-assembly process is lacking. We propose an analytically solvable kinetic model for stepwise, reversible oligomerization. In biologically relevant limits, we obtain simple algebraic expressions for the rate of pore formation, as well as for the concentration of pores as a function of time. Quantitative agreement is obtained between our model and time-resolved kinetic experiments of Bacillus thuringiensis Cry1Ac (tetrameric pore), aerolysin, Staphylococcus aureus α-haemolysin (heptameric pores) and Escherichia coli cytolysin A (dodecameric pore). Furthermore, our model explains how rapid self-assembly can take place with low concentrations of oligomeric intermediates, as observed in recent single-molecule fluorescence experiments of α-haemolysin self-assembly. We propose that suppressing the concentration of oligomeric intermediates may be the key to reliable, error-free, self-assembly of pores.

A combined network model for membrane fouling.

Membrane fouling during particle filtration occurs through a variety of mechanisms, including internal pore clogging by contaminants, coverage of pore entrances, and deposition on the membrane surface. Each of these fouling mechanisms results in a decline in the observed flow rate over time, and the decrease in filtration efficiency can be characterized by a unique signature formed by plotting the volumetric flux, Q^, as a function of the total volume of fluid processed, V^. When membrane fouling takes place via any one of these mechanisms independently the Q^V^ signature is always convex downwards for filtration under a constant transmembrane pressure. However, in many such filtration scenarios, the fouling mechanisms are inherently coupled and the resulting signature is more difficult to interpret. For instance, blocking of a pore entrance will be exacerbated by the internal clogging of a pore, while the deposition of a layer of contaminants is more likely once the pores have been covered by particulates. As a result, the experimentally observed Q^V^ signature can vary dramatically from the canonical convex-downwards graph, revealing features that are not captured by existing continuum models. In a range of industrially relevant cases we observe a concave-downwardsQ^V^ signature, indicative of a fouling rate that becomes more severe with time. We derive a network model for membrane fouling that accounts for the inter-relation between fouling mechanisms and demonstrate the impact on the Q^V^ signature. Our formulation recovers the behaviour of existing models when the mechanisms are treated independently, but also elucidates the concave-downward Q^V^ signature for multiple interactive fouling mechanisms. The resulting model enables post-experiment analysis to identify the dominant fouling modality at each stage, and is able to provide insight into selecting appropriate operating regimes.

On the predictions and limitations of the Becker-Döring 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 results from 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 timescales, in agreement with experiment and with previous theories. We conclude by exposing a limitation in the Becker-Döring theory for re-equilibration of surfactant solutions.