Source author record

Brennan Sprinkle

Brennan Sprinkle appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.soft math.NA Numerical Analysis physics.flu-dyn

Catalog footprint

What is connected

3works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

libMobility: A Python library for hydrodynamics at the Smoluchowski level

Effective hydrodynamic modeling is crucial for accurately predicting fluid-particle interactions in diverse fields such as biophysics and materials science. Developing and implementing hydrodynamic algorithms is challenging due to the complexity of fluid dynamics, necessitating efficient management of large-scale computations and sophisticated boundary conditions. Furthermore, adapting these algorithms for use on massively parallel architectures like GPUs adds an additional layer of complexity. This paper presents the libMobility software library, which offers a suite of CUDA-enabled solvers for simulating hydrodynamic interactions in particulate systems at the Rotne-Prager-Yamakawa (RPY) level. The library facilitates precise simulations of particle displacements influenced by external forces and torques, including both the deterministic and stochastic components. Notable features of libMobility include its ability to handle linear and angular displacements, thermal fluctuations, and various domain geometries effectively. With an interface in Python, libMobility provides comprehensive tools for researchers in computational fluid dynamics and related fields to simulate particle mobility efficiently. This article details the technical architecture, functionality, and wide-ranging applications of libMobility. libMobility is available at https://github.com/stochasticHydroTools/libMobility.

preprint2022arXiv

The hydrodynamics of a twisting, bending, inextensible fiber in Stokes flow

In swimming microorganisms and the cell cytoskeleton, inextensible fibers resist bending and twisting, and interact with the surrounding fluid to cause or resist large-scale fluid motion. In this paper, we develop a novel numerical method for the simulation of cylindrical fibers by extending our previous work on inextensible bending fibers [Maxian et al., Phys. Rev. Fluids 6 (1), 014102] to fibers with twist elasticity. In our "Euler" model, twist is a scalar function that measures the deviation of the fiber cross section relative to a twist-free frame, the fiber exerts only torque parallel to the centerline on the fluid, and the perpendicular components of the rotational fluid velocity are discarded in favor of the translational velocity. In the first part of this paper, we justify this model by comparing it to another commonly-used "Kirchhoff" formulation where the fiber exerts both perpendicular and parallel torque on the fluid, and the perpendicular angular fluid velocity is required to be consistent with the translational fluid velocity. We then develop a spectral numerical method for the hydrodynamics of the Euler model. We define hydrodynamic mobility operators using integrals of the Rotne-Prager-Yamakawa tensor, and evaluate these integrals through a novel slender-body quadrature, which requires on the order of 10 points along the fiber to obtain several digits of accuracy. We demonstrate that this choice of mobility removes the unphysical negative eigenvalues in the translation-translation mobility associated with asymptotic slender body theories, and ensures strong convergence of the fiber velocity and weak convergence of the fiber constraint forces. We pair the spatial discretization with a semi-implicit temporal integrator to confirm the negligible contribution of twist elasticity to the relaxation dynamics of a bent fiber and study the instability of a twirling fiber.

preprint2020arXiv

Driven dynamics in dense suspensions of microrollers

We perform detailed computational and experimental measurements of the driven dynamics of a dense, uniform suspension of sedimented microrollers driven by a magnetic field rotating around an axis parallel to the floor. We develop a lubrication-corrected Brownian Dynamics method for dense suspensions of driven colloids sedimented above a bottom wall. The numerical method adds lubrication friction between nearby pairs of particles, as well as particles and the bottom wall, to a minimally-resolved model of the far-field hydrodynamic interactions. Our experiments combine fluorescent labeling with particle tracking to trace the trajectories of individual particles in a dense suspension, and to measure their propulsion velocities. Previous computational studies [B. Sprinkle et al., J. Chem. Phys., 147, 244103, 2017] predicted that at sufficiently high densities a uniform suspension of microrollers separates into two layers, a slow monolayer right above the wall, and a fast layer on top of the bottom layer. Here we verify this prediction, showing good quantitative agreement between the bimodal distribution of particle velocities predicted by the lubrication-corrected Brownian Dynamics and those measured in the experiments. The computational method accurately predicts the rate at which particles are observed to switch between the slow and fast layers in the experiments. We also use our numerical method to demonstrate the important role that pairwise lubrication plays in motility-induced phase separation in dense monolayers of colloidal microrollers, as recently suggested for suspensions of Quincke rollers [D. Geyer et al., Physical Review X, 9(3), 031043, 2019].