No abstract

Vesicles exposed to the human circulatory system experience a wide range of flows and Reynolds numbers. Previous investigations of vesicles in fluid flow have focused on the Stokes flow regime. In this work the influence of inertia on the dynamics of a vesicle in a shearing flow is investigated using a novel level-set computational method in two dimensions. A detailed analysis of the behaviour of a single vesicle at finite Reynolds number is presented. At low Reynolds numbers the results recover vesicle behaviour previously observed for Stokes flow. At moderate Reynolds numbers the classical tumbling behaviour of highly viscous vesicles is no longer observed. Instead, the vesicle is observed to tank-tread, with an equilibrium angle dependent on the Reynolds number and the reduced area of the vesicle. It is shown that a vesicle with an inner/outer fluid viscosity ratio as high as 200 will not tumble if the Reynolds number is as low as 10. A new damped tank-treading behaviour, where the vesicle will briefly oscillate about the equilibrium inclination angle, is also observed. This behaviour is explained by an investigation on the torque acting on a vesicle in shear flow. Scaling laws for vesicles in inertial flows have also been determined. It is observed that quantities such as vesicle tumbling period follow square-root scaling with respect to the Reynolds number. Finally, the maximum tension as a function of the Reynolds number is also determined. It is observed that, as the Reynolds number increases, the maximum tension on the vesicle membrane also increases. This could play a role in the creation of stable pores in vesicle membranes or for the premature destruction of vesicles exposed to the human circulatory system.

A three-dimensional numerical model of vesicle electrohydrodynamics in the presence of DC electric fields is presented. The vesicle membrane is modeled as a thin capacitive interface through the use of a semi-implicit level set Jet scheme. The enclosed volume and surface area are conserved both locally and globally by a new Navier-Stokes projection method. The electric field calculations explicitly take into account the capacitive interface by an implicit Immersed Interface Method formulation, which calculates the electric potential field and the trans-membrane potential simultaneously. The results match well with previously published experimental, analytic and two-dimensional computational works.Comment: This work was supported by NSF Award #125373

This paper proposes and implements a novel hybrid level set method which combines the numerical efficiency of the local level set approach with the temporal stability afforded by a semi-implicit technique. By introducing an extraction/insertion algorithm into the local level set approach, we can accurately capture complicated behaviors such as interface separation and coalescence. This technique solves a well known problem when treating a semi-implicit system with spectral methods, where spurious interface movements emerge when two interfaces are close to each other. Numerical experiments show that the proposed method is stable, efficient and scales up well into three dimensional problems.