Low-resolution simulations of vesicle suspensions in 2D

Abstract Vesicle suspensions appear in many biological and industrial applications. These suspensions are characterized by rich and complex dynamics of vesicles due to their interaction with the bulk fluid, and their large deformations and nonlinear elastic properties. Many existing state-of-the-art numerical schemes can resolve such complex vesicle flows. However, even when using provably optimal algorithms, these simulations can be computationally expensive, especially for suspensions with a large number of vesicles. These high computational costs can limit the use of simulations for parameter exploration, optimization, or uncertainty quantification. One way to reduce the cost is to use low-resolution discretizations in space and time. However, it is well-known that simply reducing the resolution results in vesicle collisions, numerical instabilities, and often in erroneous results. In this paper, we investigate the effect of a number of algorithmic empirical fixes (which are commonly used by many groups) in an attempt to make low-resolution simulations more stable and more predictive. Based on our empirical studies for a number of flow configurations, we propose a scheme that attempts to integrate these fixes in a systematic way. This low-resolution scheme is an extension of our previous work [51] , [53] . Our low-resolution correction algorithms (LRCA) include anti-aliasing and membrane reparametrization for avoiding spurious oscillations in vesicles' membranes, adaptive time stepping and a repulsion force for handling vesicle collisions and, correction of vesicles' area and arc-length for maintaining physical vesicle shapes. We perform a systematic error analysis by comparing the low-resolution simulations of dilute and dense suspensions with their high-fidelity, fully resolved, counterparts. We observe that the LRCA enables both efficient and statistically accurate low-resolution simulations of vesicle suspensions, while it can be 10× to 100× faster.

[1]  E. Shaqfeh,et al.  The shape stability of a lipid vesicle in a uniaxial extensional flow , 2013, Journal of Fluid Mechanics.

[2]  P. Moin,et al.  On the Effect of Numerical Errors in Large Eddy Simulations of Turbulent Flows , 1997 .

[3]  George Biros,et al.  Author ' s personal copy Dynamic simulation of locally inextensible vesicles suspended in an arbitrary two-dimensional domain , a boundary integral method , 2010 .

[4]  C. Pozrikidis Boundary Integral and Singularity Methods for Linearized Viscous Flow: Index , 1992 .

[5]  Gerhard Gompper,et al.  Deformation and dynamics of red blood cells in flow through cylindrical microchannels. , 2014, Soft matter.

[6]  T. Biben,et al.  Steady to unsteady dynamics of a vesicle in a flow. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[7]  George Biros,et al.  Vesicle migration and spatial organization driven by flow line curvature. , 2011, Physical review letters.

[8]  Axel Voigt,et al.  Diffuse interface models of locally inextensible vesicles in a viscous fluid , 2013, J. Comput. Phys..

[9]  Peter V Coveney,et al.  Deformability-based red blood cell separation in deterministic lateral displacement devices-A simulation study. , 2014, Biomicrofluidics.

[10]  Frank Dubois,et al.  Dynamics of Vesicle Suspensions in Shear Flow Between Walls , 2011 .

[11]  Raymond Quek,et al.  Separation of deformable particles in deterministic lateral displacement devices. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  George Biros,et al.  High-volume fraction simulations of two-dimensional vesicle suspensions , 2013, J. Comput. Phys..

[13]  Denis Zorin,et al.  Contact-aware simulations of particulate Stokesian suspensions , 2016, J. Comput. Phys..

[14]  Hong Zhao,et al.  A spectral boundary integral method for flowing blood cells , 2010, J. Comput. Phys..

[15]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[16]  C. Kelley,et al.  Convergence Analysis of Pseudo-Transient Continuation , 1998 .

[17]  R. Glowinski,et al.  A fictitious domain approach to the direct numerical simulation of incompressible viscous flow past moving rigid bodies: application to particulate flow , 2001 .

[18]  George Biros,et al.  A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D , 2009, J. Comput. Phys..

[19]  G. S. Patterson,et al.  Spectral Calculations of Isotropic Turbulence: Efficient Removal of Aliasing Interactions , 1971 .

[20]  G. Gompper,et al.  Dynamics and rheology of vesicle suspensions in wall-bounded shear flow , 2013 .

[21]  Gerhard Gompper,et al.  Behavior of rigid and deformable particles in deterministic lateral displacement devices with different post shapes. , 2015, The Journal of chemical physics.

[22]  Seifert,et al.  Fluid Vesicles in Shear Flow. , 1996, Physical review letters.

[23]  A. Acrivos,et al.  The shear-induced migration of particles in concentrated suspensions , 1987, Journal of Fluid Mechanics.

[24]  G. Biros,et al.  Complexity of vesicle microcirculation. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[25]  Jens Harting,et al.  Interplay between microdynamics and macrorheology in vesicle suspensions. , 2014, Soft matter.

[26]  Steven A. Orszag,et al.  Numerical simulation of incompressible flows within simple boundaries: accuracy , 1971, Journal of Fluid Mechanics.

[27]  T. Krüger,et al.  Breakdown of deterministic lateral displacement efficiency for non-dilute suspensions: A numerical study. , 2015, Medical engineering & physics.

[28]  C. Pozrikidis,et al.  The axisymmetric deformation of a red blood cell in uniaxial straining Stokes flow , 1990, Journal of Fluid Mechanics.

[29]  Chaouqi Misbah,et al.  Vacillating breathing and tumbling of vesicles under shear flow. , 2006, Physical review letters.

[30]  George Biros,et al.  A fast algorithm for simulating vesicle flows in three dimensions , 2011, J. Comput. Phys..

[31]  Dhairya Malhotra,et al.  A parallel algorithm for long-timescale simulation of concentrated vesicle suspensions in three dimensions , 2017 .

[32]  Shangjun Ye,et al.  Effects of the particle deformability on the critical separation diameter in the deterministic lateral displacement device , 2014, Journal of Fluid Mechanics.

[33]  T. Biben,et al.  Tumbling of vesicles under shear flow within an advected-field approach. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[34]  Z. Feng,et al.  The immersed boundary-lattice Boltzmann method for solving fluid-particles interaction problems , 2004 .

[35]  E. Sackmann,et al.  Supported Membranes: Scientific and Practical Applications , 1996, Science.

[36]  Bradley K. Alpert,et al.  Hybrid Gauss-Trapezoidal Quadrature Rules , 1999, SIAM J. Sci. Comput..

[37]  Gerhard Gompper,et al.  Margination of micro- and nano-particles in blood flow and its effect on drug delivery , 2014, Scientific Reports.

[38]  T. Hou,et al.  Removing the stiffness from interfacial flows with surface tension , 1994 .

[39]  Eitan Grinspun,et al.  Asynchronous variational contact mechanics , 2010, 1007.3240.

[40]  George Biros,et al.  Adaptive time stepping for vesicle suspensions , 2014, J. Comput. Phys..

[41]  Hong Zhao,et al.  The dynamics of a vesicle in simple shear flow , 2011, Journal of Fluid Mechanics.

[42]  Jason P Beech,et al.  Sorting cells by size, shape and deformability. , 2012, Lab on a chip.

[43]  J. Freund Leukocyte Margination in a Model Microvessel , 2006 .

[44]  Hong Zhao,et al.  Coarse-grained theory to predict the concentration distribution of red blood cells in wall-bounded Couette flow at zero Reynolds number , 2013 .

[45]  Udo Seifert,et al.  Configurations of fluid membranes and vesicles , 1997 .

[46]  H. Noguchi,et al.  Shape transitions of fluid vesicles and red blood cells in capillary flows. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[47]  Steven J. Ruuth,et al.  Implicit-explicit methods for time-dependent partial differential equations , 1995 .

[48]  Gerhard Gompper,et al.  Margination of white blood cells in microcapillary flow. , 2012, Physical review letters.

[49]  C. Canuto Spectral methods in fluid dynamics , 1991 .

[50]  Z. Feng,et al.  Hydrodynamic force on spheres in cylindrical and prismatic enclosures , 2002 .

[51]  Hong Zhao,et al.  The dynamics of a non-dilute vesicle suspension in a simple shear flow , 2013, Journal of Fluid Mechanics.

[52]  George Biros,et al.  High-order Adaptive Time Stepping for Vesicle Suspensions with Viscosity Contrast☆ , 2014, 1409.0212.

[53]  Adrianna Gillman,et al.  A Fast Algorithm for Simulating Multiphase Flows Through Periodic Geometries of Arbitrary Shape , 2015, SIAM J. Sci. Comput..

[54]  Steven A. Orszag,et al.  Fourier Series on Spheres , 1974 .

[55]  Victor Steinberg,et al.  Dynamics of interacting vesicles and rheology of vesicle suspension in shear flow , 2008 .

[56]  George Biros,et al.  Boundary integral method for the flow of vesicles with viscosity contrast in three dimensions , 2015, J. Comput. Phys..

[57]  G. Biros,et al.  Quantification of mixing in vesicle suspensions using numerical simulations in two dimensions. , 2016, Physics of fluids.

[58]  Eugene C. Eckstein,et al.  Self-diffusion of particles in shear flow of a suspension , 1977, Journal of Fluid Mechanics.

[59]  Michael T. Laub,et al.  Continuous Particle Separation Through Deterministic Lateral Displacement , 2004 .

[60]  Christian Wagner,et al.  Vesicle dynamics in a confined Poiseuille flow: from steady state to chaos. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[61]  M. Shelley,et al.  On the connection between thin vortex layers and vortex sheets , 1990, Journal of Fluid Mechanics.

[62]  David J. Jeffrey,et al.  The Rheological Properties of Suspensions of Rigid Particles , 1976 .

[63]  Jonathan B. Freund,et al.  A high-resolution fast boundary-integral method for multiple interacting blood cells , 2010 .

[64]  David Harmon,et al.  Asynchronous contact mechanics , 2009, SIGGRAPH 2009.

[65]  Roland Glowinski,et al.  A DLM/FD/IB method for simulating cell/cell and cell/particle interaction in microchannels , 2010 .

[66]  Paul J. Besl,et al.  A Method for Registration of 3-D Shapes , 1992, IEEE Trans. Pattern Anal. Mach. Intell..

[67]  J. Freund,et al.  Cellular flow in a small blood vessel , 2010, Journal of Fluid Mechanics.