Influence of coexisting phases on the surface dilatational viscosity of Langmuir monolayers

Coupling Vortical Bulk Flows to the Air–Water Interface: From Putting Oil on Troubled Waters to Surfactants on Protein Solutions

The air–water interface in flowing systems remains a challenge to model, even in cases where the interface is essentially flat. This is because even though each side is governed by the Navier–Stokes equations, the stress balance which provides the boundary conditions for the equations involves properties associated with surfactants that are inevitably present at the air–water interface. Aside from challenges in measuring interfacial properties, either intrinsic or flow-dependent, the two-way coupling of bulk and interfacial flows is non-trivial, even for very simple flow geometries. Here, we present an overview of the physics associated with surfactant monolayers of flowing liquid and describe how the monolayer affects the bulk flow and how the monolayer is transported and deformed by the bulk flow. The emphasis is primarily on cylindrical flow geometries, and both Newtonian and non-Newtonian interfacial responses are considered. We consider interfacial flows that are solenoidal as well as those where the surface velocity is not divergence free.

A finite difference method with subsampling for immersed boundary simulations of the capsule dynamics with viscoelastic membranes

The membrane or interfacial viscosity is an important property in many multiphase and biofluidic situations, such as the red blood cell dynamics and emulsion stability. The immersed boundary method (IBM), which incorporates the dynamic flow-membrane interaction via force distribution and velocity interpolation, has been extensively employed in simulations of such systems. Unfortunately, direct implementation of membrane viscosity in IBM suffers severe numerical instability, which causes an IBM calculation to break down before generating any useful results. Few attempts have been recently reported; however, several concerns exist in these attempts, such as the inconsistency to the classical definition of membrane viscosity, the inability to model the shear and dilatational viscosities separately, the unjustified mathematical formulations, and the complicated algorithms and computation. To overcome these concerns, in this paper, we propose a finite difference approach for implementing membrane viscosity in immersed boundary simulations. The viscous stress is obtained via finite difference approximations to the differential strain-stress relationship, with the help of a subsampling scheme to reduce the numerical noise in the calculated strain rates. This simple method has also avoided the complicated matrix calculations in previous attempts, and hence, a better computational efficiency is expected. Detailed mathematical description of the method and key steps for its implementation in immersed boundary programs are provided. Validation and illustration calculations are performed, and our results are compared with analytical solutions and previous publications with satisfactory agreement. The influences of membrane mesh resolution and simulation time step are also examined; and the results show no indication that our finite difference method has downgraded the general IBM accuracy. Based on these simulations and analysis, we believe that our method would be a better choice for future IBM simulations of capsule dynamics with viscoelastic membranes.

Predicting Steady Shear Rheology of Condensed-Phase Monomolecular Films at the Air-Water Interface

Predicting the non-Newtonian shear response of soft interfaces in biophysical systems and engineered products has been compromised by the use of linear (Newtonian) constitutive equations. We present a generalized constitutive equation, with tractable material properties, governing the response of Newtonian and non-Newtonian interfaces subjected to a wide range of steady shear. With experiments spanning six decades of shear rate, we capture and unify divergent reports of shear-thinning behavior of monomolecular films of the lipid dipalmitoylphosphatidylcholine, the primary constituent of mammalian cell walls and lung surfactant, at near-physiological packing densities.

Bursting bubbles and bilayers

This paper discusses various interactions between ultrasound, phospholipid monolayer-coated gas bubbles, phospholipid bilayer vesicles, and cells. The paper begins with a review of microbubble physics models, developed to describe microbubble dynamic behavior in the presence of ultrasound, and follows this with a discussion of how such models can be used to predict inertial cavitation profiles. Predicted sensitivities of inertial cavitation to changes in the values of membrane properties, including surface tension, surface dilatational viscosity, and area expansion modulus, indicate that area expansion modulus exerts the greatest relative influence on inertial cavitation. Accordingly, the theoretical dependence of area expansion modulus on chemical composition - in particular, poly (ethylene glyclol) (PEG) - is reviewed, and predictions of inertial cavitation for different PEG molecular weights and compositions are compared with experiment. Noteworthy is the predicted dependence, or lack thereof, of inertial cavitation on PEG molecular weight and mole fraction. Specifically, inertial cavitation is predicted to be independent of PEG molecular weight and mole fraction in the so-called mushroom regime. In the “brush” regime, however, inertial cavitation is predicted to increase with PEG mole fraction but to decrease (to the inverse 3/5 power) with PEG molecular weight. While excellent agreement between experiment and theory can be achieved, it is shown that the calculated inertial cavitation profiles depend strongly on the criterion used to predict inertial cavitation. This is followed by a discussion of nesting microbubbles inside the aqueous core of microcapsules and how this significantly increases the inertial cavitation threshold. Nesting thus offers a means for avoiding unwanted inertial cavitation and cell death during imaging and other applications such as sonoporation. A review of putative sonoporation mechanisms is then presented, including those involving microbubbles to deliver cargo into a cell, and those - not necessarily involving microubbles - to release cargo from a phospholipid vesicle (or reverse sonoporation). It is shown that the rate of (reverse) sonoporation from liposomes correlates with phospholipid bilayer phase behavior, liquid-disordered phases giving appreciably faster release than liquid-ordered phases. Moreover, liquid-disordered phases exhibit evidence of two release mechanisms, which are described well mathematically by enhanced diffusion (possibly via dilation of membrane phospholipids) and irreversible membrane disruption, whereas liquid-ordered phases are described by a single mechanism, which has yet to be positively identified. The ability to tune release kinetics with bilayer composition makes reverse sonoporation of phospholipid vesicles a promising methodology for controlled drug delivery. Moreover, nesting of microbubbles inside vesicles constitutes a truly “theranostic” vehicle, one that can be used for both long-lasting, safe imaging and for controlled drug delivery.

Effects of shearing flow with inertia on monolayer mesoscale structure

In an experimental flow system capable of imparting a well-controlled shear-rate distribution with inertia to a monolayer consisting of coexisting phases, we have studied the resulting phase morphology and domain fragmentation. These evolve on distinct time scales: the viscous time associated with the viscosity in the bulk and the Marangoni stress and the fragmentation/relaxation time associated with the phase morphology. A relationship between the microstructure (line tension) and macroflow (shear rate) determining the meso length scale of the coexisting phase domains has been deduced from dimensional analysis and was found to correlate well with the quantitative experimental observations.