Recent experimental studies have shown that confinement can profoundly affect selforganization in semi-dilute active suspensions, leading to striking features such as the formation of steady and spontaneous vortices in circular domains and the emergence of unidirectional pumping motions in periodic racetrack geometries. Motivated by these findings, we analyze the two-dimensional dynamics in confined suspensions of active self-propelled swimmers using a mean-field kinetic theory where conservation equations for the particle configurations are coupled to the forced Navier-Stokes equations for the self-generated fluid flow. In circular domains, a systematic exploration of the parameter space casts light on three distinct states: equilibrium with no flow, stable vortex, and chaotic motion, and the transitions between these are explained and predicted quantitatively using a linearized theory. In periodic racetracks, similar transitions from equilibrium to net pumping to traveling waves to chaos are observed in agreement with experimental observations and are also explained theoretically. Our results underscore the subtle effects of geometry on the morphology and dynamics of emerging patterns in active suspensions and pave the way for the control of active collective motion in microfluidic devices.

No abstract

No abstract

We present a solver for the Poisson-Boltzmann equation and demonstrate its applicability for biomolecular electrostatics computation. The solver uses a level set framework to represent sharp, complex interfaces in a simple and robust manner. It also uses non-graded, adaptive octree grids which, in comparison to uniform grids, drastically decrease memory usage and runtime without sacrificing accuracy. The basic solver was introduced in earlier works [16,27], and here is extended to address biomolecular systems. First, a novel approach of calculating the solvent excluded and the solvent accessible surfaces is explained; this allows to accurately represent the location of the molecule’s surface. Next, a hybrid finite difference/finite volume approach is presented for discretizing the nonlinear Poisson-Boltzmann equation and enforcing the jump boundary conditions at the interface. Since the interface is implicitly represented by a level set function, imposing the jump boundary conditions is straightforward and efficient.

No abstract