Computing collision stress in assemblies of active spherocylinders: Applications of a fast and generic geometric method.

In this work, we provide a solution to the problem of computing collision stress in particle-tracking simulations. First, a formulation for the collision stress between particles is derived as an extension of the virial stress formula to general-shaped particles with uniform or non-uniform mass density. Second, we describe a collision-resolution algorithm based on geometric constraint minimization which eliminates the stiff pairwise potentials in traditional methods. The method is validated with a comparison to the equation of state of Brownian spherocylinders. Then we demonstrate the application of this method in several emerging problems of soft active matter.

[1]  Anna-Karin Tornberg,et al.  Gravity induced sedimentation of slender fibers , 2009 .

[2]  J. Brady,et al.  Constant Stress and Pressure Rheology of Colloidal Suspensions. , 2015, Physical review letters.

[3]  M. Shelley,et al.  Connecting macroscopic dynamics with microscopic properties in active microtubule network contraction , 2017, 1706.10235.

[4]  Markus Bär,et al.  Large-scale collective properties of self-propelled rods. , 2009, Physical review letters.

[5]  John F. Brady,et al.  Brownian Dynamics simulation of hard-sphere colloidal dispersions , 2000 .

[6]  M Cristina Marchetti,et al.  Enhanced diffusion and ordering of self-propelled rods. , 2008, Physical review letters.

[7]  Reinhard Lipowsky,et al.  Enhanced ordering of interacting filaments by molecular motors. , 2006, Physical review letters.

[8]  Leslie Greengard,et al.  An integral equation formulation for rigid bodies in Stokes flow in three dimensions , 2016, J. Comput. Phys..

[9]  Ailing Gong,et al.  The stress tensor in a two-dimensional granular shear flow , 1986, Journal of Fluid Mechanics.

[10]  D. Eberly.  Robust Computation of Distance Between Line Segments , 2015 .

[11]  S. Fang.  A linearization method for generalized complementarity problems , 1984 .

[12]  Hysteresis, reentrance, and glassy dynamics in systems of self-propelled rods. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  S. Prager,et al.  Variational Treatment of Hydrodynamic Interaction in Polymers , 1969 .

[14]  Density functional theory for hard spherocylinders: phase transitions in the bulk and in the presence of external fields , 1999 .

[15]  Yu-Guo Tao,et al.  Brownian dynamics simulations of the self- and collective rotational diffusion coefficients of rigid long thin rods. , 2005, The Journal of chemical physics.

[16]  M. Anitescu,et al.  Large-scale parallel multi-body dynamics with frictional contact on the graphical processing unit , 2008 .

[17]  Roger Fletcher,et al.  Projected Barzilai-Borwein methods for large-scale box-constrained quadratic programming , 2005, Numerische Mathematik.

[18]  Jason E. Butler,et al.  Normal stress differences in suspensions of rigid fibres , 2014, Journal of Fluid Mechanics.

[19]  L. Onsager.  THE EFFECTS OF SHAPE ON THE INTERACTION OF COLLOIDAL PARTICLES , 1949 .

[20]  Andreas Deutsch,et al.  Self-propelled rods exhibit a phase-separated state characterized by the presence of active stresses and the ejection of polar clusters. , 2015, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  M. Anitescu,et al.  A matrix-free cone complementarity approach for solving large-scale, nonsmooth, rigid body dynamics , 2011 .

[22]  Shravan K. Veerapaneni,et al.  Boundary integral equation analysis for suspension of spheres in Stokes flow , 2017, J. Comput. Phys..

[23]  D. W. Rebertus,et al.  Molecular dynamics simulation of a fluid of hard spherocylinders , 1977 .

[24]  G J Vroege,et al.  Phase transitions in lyotropic colloidal and polymer liquid crystals , 1992 .

[25]  F. Peruani.  Active Brownian rods , 2015, 1512.07567.

[26]  Katarina Gustavsson,et al.  A numerical method for simulations of rigid fiber suspensions , 2006, J. Comput. Phys..

[27]  Charles S. Campbell,et al.  The stress tensor for simple shear flows of a granular material , 1989, Journal of Fluid Mechanics.

[28]  Hammad Mazhar,et al.  Using Nesterov's Method to Accelerate Multibody Dynamics with Friction and Contact , 2015, ACM Trans. Graph..

[29]  Mu Wang,et al.  Spectral Ewald Acceleration of Stokesian Dynamics for polydisperse suspensions , 2015, J. Comput. Phys..

[30]  Markus Bär,et al.  Mesoscale pattern formation of self-propelled rods with velocity reversal. , 2016, Physical review. E.

[31]  Aleksandar Donev,et al.  Brownian dynamics of confined rigid bodies. , 2015, The Journal of chemical physics.

[32]  Löwen.  Brownian dynamics of hard spherocylinders. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[33]  D. Stewart,et al.  Time-stepping for three-dimensional rigid body dynamics , 1999 .

[34]  D. Boyer,et al.  Order, intermittency, and pressure fluctuations in a system of proliferating rods. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[35]  Peter G. Bolhuis,et al.  Tracing the phase boundaries of hard spherocylinders , 1997 .

[36]  J. Brady,et al.  Swim pressure: stress generation in active matter. , 2014, Physical review letters.

[37]  Bertrand Maury,et al.  A time-stepping scheme for inelastic collisions , 2006, Numerische Mathematik.

[38]  Sarah Niebe,et al.  Numerical methods for linear complementarity problems in physics-based animation , 2015, SIGGRAPH '13.