Search

We numerically investigate the hydrodynamics and membrane dynamics of multicomponent vesicles in two strongly confined geometries. This serves as a simplified model for red blood cells undergoing large deformations while traversing narrow constrictions. We propose a new parameterization for the bending modulus that remains positive for all lipid phase parameter values. For a multicomponent vesicle passing through a stenosis, we establish connections between various properties: lipid phase coarsening, size and flow profile of the lubrication layers, excess pressure, and the tank-treading velocity of the membrane. For a multicomponent vesicle passing through a contracting channel, we find that the lipid always phase separates so that the vesicle is stiffer in the front as it passes through the constriction. For both cases of confinement, we find that lipid coarsening is arrested under strong confinement and proceeds at a high rate upon relief from extreme confinement. The results may be useful for efficient sorting lipid domains using microfluidic flows by controlled release of vesicles passing through strong confinement

Active colloidal systems with non-equilibrium self-organization is a long-standing, challenging area in biology. To understand how hydrodynamic flow may be used to actively control self-assembly of Janus particles (JPs), we use a model recently developed for the many-body hydrodynamics of amphiphilic JPs suspended in a viscous background flow (JFM, 941, 2022). We investigate how various morphologies arise from tuning the hydrophobic distribution of the JP-solvent interface. We find JPs assembled into uni-lamella, multi-lamella and striated structures. To introduce dynamics, we include a linear shear flow and a steady Taylor-Green mixing flow, and measure the collective dynamics of JP particles in terms of their (a) free energy from the hydrophobic interactions between the JPs, (b) order parameter for the ordering of JPs in terms of alignment of their directors, and (c) strain parameter that captures the deformation in the assembly. We characterize the effective material properties of the JP structures and find that the uni-lamellar structures increases orientation order under shear flow, the multilamellar structure behaves as a shear thinning fluid, and the striated structure possesses a yield stress. These numerical results provide insights into dynamic control of non-equilibrium active biological systems with similar self-organization.

Lipid bilayer membranes have a native (albeit small) permeability for water molecules. Under an external load, provided that the bilayer structure stays intact and does not suffer from poration or rupture, a lipid membrane deforms and its water influx/efflux is often assumed negligible in the absence of osmolarity. In this work we use boundary integral simulations to investigate the effects of water permeability on the vesicle hydrodynamics due to a mechanical load, such as the viscous stress from an external flow deforming a vesicle membrane in free space or pushing it through a confinement. Incorporating the membrane permeability into the framework of Helfrich free energy for an inextensible, elastic membrane as a model for a semipermeable vesicle, we illustrate that, in the absence of an osmotic stress gradient, the semipermeable vesicle is affected by water influx/efflux over a sufficiently long time or under a strong confinement. Our simulations quantify the conditions for water permeation to be negligible in terms of the time scales, flow strength, and confinement. These results shed light on how microfluidic confinement can be utilized to estimate membrane permeability.

The physico-chemistry of surfactants (amphiphilic surface active agents) is often used to control the dynamics of viscous drops and bubbles. Surfactant sorption kinetics has been shown to play a critical role in the deformation of drops in extensional and shear flows, yet to the best of our knowledge these kinetics effects on a viscous drop in an electric fieldhave not been accounted for. In this paper we numerically investigate the effects of sorption kinetics on a surfactant-covered viscous drop in an electric field. Over a range of electric conductivity and permittivity ratios between the interior and exterior fluids, we focus on the dependence of deformation and flow on the transfer parameter $J$, and Biot number $\text{Bi}$ that characterize the extent of surfactant exchange between the drop surface and the bulk. Our findings suggest solubility affects the electrohydrodynamics of a viscous drop in distinct ways as we identify parameter regions where (1) surfactant solubility may alter both the drop deformation and circulation of fluid around a drop, and (2) surfactant solubility affects mainly the flow and not the deformation.

The dynamics of an adhesive two-dimensional vesicle doublet under various flow conditions is investigated numerically using a high-order, adaptive-in-time boundary integral method. In a quiescent flow, two nearby vesicles move slowly towards each other under the adhesive potential, pushing out fluid between them to form a vesicle doublet at equilibrium. A lubrication analysis on such draining of a thin film gives the dependencies of draining time on adhesion strength and separation distance that are in good agreement with numerical results. In a planar extensional flow we find a stable vesicle doublet forms only when two vesicles collide head-on around the stagnation point. In a microfluid trap where the stagnation of an extensional flow is dynamically placed in the middle of a vesicle doublet through an active control loop, novel dynamics of a vesicle doublet are observed. Numerical simulations show that there exists a critical extensional flow rate above which adhesive interaction is overcome by the diverging stream, thus providing a simple method to measure the adhesion strength between two vesicle membranes. In a planar shear flow, numerical simulations reveal that a vesicle doublet may form provided that the adhesion strength is sufficiently large at a given vesicle reduced area. Once a doublet is formed, its oscillatory dynamics is found to depend on the adhesion strength and their reduced area. Furthermore the effective shear viscosity of a dilute suspension of vesicle doublets is found to be a function of the reduced area. Results from these numerical studies and analysis shed light on the hydrodynamic and rheological consequences of adhesive interactions between vesicles in a viscous fluid.

The turbulent solar convection zone has remained one of the most challenging and important subjects in physics. Understanding the complex dynamics in the solar con- vection zone is crucial for gaining insight into the solar dynamo problem. Many solar observatories have generated revealing data with great details of large scale motions in the solar convection zone. For example, a strong di erential rotation is observed: the angular rotation is observed to be faster at the equator than near the poles not only near the solar surface, but also deep in the convection zone. On the other hand, due to the wide range of dynamical scales of turbulence in the solar convection zone, both theory and simulation have limited success. Thus, cutting edge solar models and numerical simulations of the solar convection zone have focused more narrowly on a few key features of the solar convection zone, such as the time-averaged di erential rotation. For example, Brun & Toomre (2002) report computational finding of differential rotation in an anelastic model for solar convection. A critical shortcoming in this model is that the viscous dissipation is based on application of mixing length theory to stellar dynamics with some ad hoc parameter tuning. The goal of our work is to implement the subgrid scale model developed at CTR into the solar simulation code and examine how the differential rotation will be a affected as a result. Specifically, we implement a Smagorinsky-Lilly subgrid scale model into the ASH (anelastic spherical harmonic) code developed over the years by various authors. This paper is organized as follows. In x2 we briefly formulate the anelastic system that describes the solar convection. In x3 we formulate the Smagorinsky-Lilly subgrid scale model for unstably stratifed convection. We then present some preliminary results in x4, where we also provide some conclusions and future directions.

Thus we conduct numerical simulations of multiphase fluids stirred by two-dimensional turbulence to assess the possibility of self-similar drop size distribution in turbulence. In our turbulence simulations, we also explore the non-diffusive limit, where molecular mobility for the interface is vanishing. Special care is needed to transport the non-diffusive interface. Numerically, we use the particle level set method to evolve the interface. Instead of using the usual methods to calculate the surface tension force from the level set function, we reconstruct the interface based on phase- field modeling, and calculate the continuum surface tension forcing from the reconstructed interface.

In the present contribution we develop a level set method based on local anisotropic Cartesian adaptation as described in Ham et al. (2002). Such an approach should allow for the smallest possible Cartesian grid capable of resolving a given flow. The remainder of the paper is organized as follows. In section 2 the level set formulation for free surface calculations is presented and its strengths and weaknesses relative to the other free surface methods reviewed. In section 3 the collocated numerical method is described. In section 4 the method is validated by solving the 2D and 3D drop oscilation problem. In section 5 we present some results from more complex cases including the 3D drop breakup in an impulsively accelerated free stream, and the 3D immiscible Rayleigh-Taylor instability. Conclusions are given in section 6.