Electro-deformation of a moving boundary: a drop interface and a lipid bilayer membrane

This dissertation focuses on the deformation of a viscous drop and a vesicle immersed in a (leaky) dielectric fluid under an electric field. A number of mathematical tools, both analytical and numerical, are developed for these investigations. The dissertation is divided into three parts. First, a large-deformation model is developed to capture the equilibrium deformation of a viscous spheroidal drop covered with non-diffusing insoluble surfactant under a uniform direct current (DC) electric field. The large- deformation model predicts the dependence of equilibrium spheroidal drop shape on the permittivity ratio, conductivity ratio, surfactant coverage, and the elasticity number. Results from the model are carefully compared against the small-deformation (quasispherical) analysis, experimental data and numerical simulation results in the literature. Moreover, surfactant effects, such as tip stretching and surface dilution effects, are greatly amplified at large surfactant coverage and high electric capillary number. These effects are well captured by the spheroidal model, but cannot be described in the second-order small-deformation theory. The large-deformation spheroidal model is then extended to study the equilibrium deformation of a giant unilamellar vesicle (GUV) under an alternating current (AC) electric field. The vesicle membrane is modeled as a thin capacitive spheroidal shell and the equilibrium vesicle shape is computed from balancing the mechanical forces between the fluid, the membrane and the imposed electric field. Detailed comparison against both experiments and small-deformation theory shows that the spheroidal model gives better agreement with experiments in terms of the dependence on fluid conductivity ratio, electric field strength and frequency, and vesicle size. Asymptotic analysis is conducted to compute the crossover frequency where a prolate vesicle crosses over to an oblate shape, and comparisons show the spheroidal model gives better agreement with experimental observations. Finally, a numerical scheme based on immersed interface method for two-phase fluids is developed to simulate the time-dependent dynamics of an axisymmetric drop in an electric field. The second-order immersed interface method is applied to solving both the fluid velocity field and the electric field. To date this has not been done before in the literature. Detailed numerical studies on this new numerical scheme shows numerical convergence and good agreement with the large-deformation model. Dynamics of an axisymmetric viscous drop under an electric field is being simulated using this novel numerical code

Electrohydrodynamic Simulations of Capsule Deformation Using a Dual Time-Stepping Lattice Boltzmann Scheme

Capsules are fluid-filled, elastic membranes that serve as a useful model for synthetic and biological membranes. One prominent application of capsules is their use in modeling the response of red blood cells to external forces. These models can be used to study the cell’s material properties and can also assist in the development of diagnostic equipment. In this work we develop a three dimensional model for numerical simulations of red blood cells under the combined influence of hydrodynamic and electrical forces. The red blood cell is modeled as a biconcave-shaped capsule suspended in an ambient fluid domain. Cell deformation occurs due to fluid motion and electrical forces that arise due to differences in the electrical properties between the internal fluid, external fluid, and cell membrane. The electrostatic equations are solved using the immersed interface method. A finite element method is used to compute the membrane’s elastic forces and the membrane’s bending resistance is described by the Helfrich bending energy functional. The membrane forces are coupled to the fluid equations through the immersed boundary method, where the elastic, bending, and electric forces appear as force densities in the Navier-Stokes equations. The fluid equations are solved using a novel dual time-stepping (DTS) lattice Boltzmann method (LBM), which decouples the fluid and capsule discretizations. The computational efficiency of the DTS scheme is studied for capsules in shear flow where it is found that the newly proposed scheme decreases computational time by a factor of 10 when compared to the standard LBM capsule model. The method is then used to study the dynamics of spherical and biconcave capsules in a combined shear flow and DC electric field. For spherical capsules the effect of field strength, shear rate, membrane capacitance, and membrane conductance are studied. For biconcave capsules the effect of the electric field on the tumbling and tank-treading modes of biconcave capsules is discussed

Electrohydrodynamic Simulations of the Deformation of Liquid-Filled Capsules

A comprehensive two- and three-dimensional framework for the electrohydrodynamic simulation of deformable capsules is provided. The role of a direct current (DC) electric field on the deformation and orientation of a liquid-filled capsule is thoroughly considered numerically. This framework is based on lattice Boltzmann method for the fluid, finite element method for the membrane structure of the capsule, fast immersed interface method for the electric field and immersed boundary method being used to consider the fluid-structure-electric interaction. Under the effect of electric field, two different types of equilibrium states, prolate or oblate are obtained. The numerical algorithm is also applied to study the interfacial tension droplet and red blood cell under shear flow. The capsules are more deformed and arrive at equilibrium status more quickly under stronger electric field. Bending stiffness will suppress the deformation and cause transition from tank-treading to tumbling for the red blood cell. However, the applied electric field will slow down the transition from tank-treading to the tumbling motion or even stay in the tank-treading motion with stronger electric field

Spectrally-Accurate Close Evaluation Schemes for Stokes Boundary Integral Operators

Dense particulate flow simulations using integral equation methods demand accurate evaluation of Stokes layer potentials on arbitrarily close interfaces. In this thesis, two spectrally-accurate integration schemes for close evaluation of 2D Stokes layer potentials are developed -- a global quadrature for the moving particles (e.g., blood cells, vesicles) represented as smooth closed curves, and an adaptive panel quadrature for the stationary boundaries (e.g., vessel walls, microfluidic channels) which are more complex curves that can be non-smooth. Both schemes rely on expressing Stokes layer potentials in terms of Laplace potentials and related complex contour integrals, which are then evaluated accurately either through a singularity cancellation technique or using analytic expressions. Numerical examples are presented to demonstrate the robustness and super-algebraic convergence of both schemes. Finally, as an application of the integration schemes, we investigate the electrohydrodynamic interactions between (possibly deflated) vesicles, where interesting behaviors unique to vesicles, such as circulatory and oscillatory motions, are observed and analyzed.PHDApplied and Interdisciplinary MathematicsUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/153404/1/boweiwu_1.pd