Search

We present a novel computational framework for the simulation of rigid spherical Janus particle suspensions in Stokes flow. For a wide array of Janus particle types, we show long-range interactions may be resolved using fast, spectrally accurate boundary integral methods. We incorporate this to our rigid body Stokes platform, which resolves hydrodynamic interactions and contact. Our approach features the use of spherical harmonic expansions for spectrally accurate integral operator evaluation, complementarity-based collision resolution, and optimal O (n) scaling with the number of particles when accelerated via fast summation techniques. We demonstrate the versatility of our simulation platform through three test cases involving Janus particle systems prominent in applications: amphiphilic, bipolar electric and phoretic particles. For each test case, we formulate Janus particle interactions in boundary integral form and use our framework to demonstrate examples of self-assembly and complex collective behavior characteristic of these systems. [ABSTRACT FROM AUTHOR]

• Collision resolution by convex quadratic programming. • Scalable and efficient parallel algorithm and implementation. • Accurate second-kind boundary integral formulation. We describe a computational framework for simulating suspensions of rigid particles in Newtonian Stokes flow. One central building block is a collision-resolution algorithm that overcomes the numerical constraints arising from particle collisions. This algorithm extends the well-known complementarity method for non-smooth multi-body dynamics to resolve collisions in dense rigid body suspensions. This approach formulates the collision resolution problem as a linear complementarity problem with geometric 'non-overlapping' constraints imposed at each time-step. It is then reformulated as a constrained quadratic programming problem and the Barzilai-Borwein projected gradient descent method is applied for its solution. This framework is designed to be applicable for any convex particle shape, e.g., spheres and spherocylinders, and applicable to any Stokes mobility solver, including the Rotne-Prager-Yamakawa approximation, Stokesian Dynamics, and PDE solvers (e.g., boundary integral and immersed boundary methods). In particular, this method imposes Newton's Third Law and records the entire contact network. Further, we describe a fast, parallel, and spectrally-accurate boundary integral method tailored for spherical particles, capable of resolving lubrication effects. We show weak and strong parallel scalings up to 8 × 10 4 particles with approximately 4 × 10 7 degrees of freedom on 1792 cores. We demonstrate the versatility of this framework with several examples, including sedimentation of particle clusters, and active matter systems composed of ensembles of particles driven to rotate. [ABSTRACT FROM AUTHOR]