No Arabic abstract
We present a generalization of the inertial coupling (IC) [Usabiaga et al. J. Comp. Phys. 2013] which permits the resolution of radiation forces on small particles with arbitrary acoustic contrast factor. The IC method is based on a Eulerian-Lagrangian approach: particles move in continuum space while the fluid equations are solved in a regular mesh (here we use the finite volume method). Thermal fluctuations in the fluid stress, important below the micron scale, are also taken into account following the Landau-Lifshitz fluid description. Each particle is described by a minimal cost resolution which consists on a single small kernel (bell-shaped function) concomitant to the particle. The main role of the particle kernel is to interpolate fluid properties and spread particle forces. Here, we extend the kernel functionality to allow for an arbitrary particle compressibility. The particle-fluid force is obtained from an imposed no-slip constraint which enforces similar particle and kernel fluid velocities. This coupling is instantaneous and permits to capture the fast, non-linear effects underlying the radiation forces on particles. Acoustic forces arise either because an excess in particle compressibility (monopolar term) or in mass (dipolar contribution) over the fluid values. Comparison with theoretical expressions show that the present generalization of the IC method correctly reproduces both contributions. Due to its low computational cost, the present method allows for simulations with many particles using a standard Graphical Processor Unit (GPU).
We present a numerical method that consistently implements thermal fluctuations and hydrodynamic interactions to the motion of Brownian particles dispersed in incompressible host fluids. In this method, the thermal fluctuations are introduced as random forces acting on the Brownian particles. The hydrodynamic interactions are introduced by directly resolving the fluid motions with the particle motion as a boundary condition to be satisfied. The validity of the method has been examined carefully by comparing the present numerical results with the fluctuation-dissipation theorem whose analytical form is known for dispersions of a single spherical particle. Simulations are then performed for more complicated systems, such as a dispersion composed of many spherical particles and a single polymeric chain in a solvent.
Frictional forces affect the rheology of hard-sphere colloids, at high shear rate. Here we demonstrate, via numerical simulations, that they also affect the dynamics of active Brownian particles, and their motility induced phase separation. Frictional forces increase the angular diffusivity of the particles, in the dilute phase, and prevent colliding particles from resolving their collision by sliding one past to the other. This leads to qualitatively changes of motility-induced phase diagram in the volume-fraction motility plane. While frictionless systems become unstable towards phase separation as the motility increases only if their volume fraction overcomes a threshold, frictional system become unstable regardless of their volume fraction. These results suggest the possibility of controlling the motility induced phase diagram by tuning the roughness of the particles.
Recent experimental studies have demonstrated that cellular motion can be directed by topographical gradients, such as those resulting from spatial variations in the features of a micropatterned substrate. This phenomenon, known as topotaxis, is especially prominent among cells persistently crawling within a spatially varying distribution of cell-sized obstacles. In this article we introduce a toy model of topotaxis based on active Brownian particles constrained to move in a lattice of obstacles, with space-dependent lattice spacing. Using numerical simulations and analytical arguments, we demonstrate that topographical gradients introduce a spatial modulation of the particles persistence, leading to directed motion toward regions of higher persistence. Our results demonstrate that persistent motion alone is sufficient to drive topotaxis and could serve as a starting point for more detailed studies on self-propelled particles and cells.
Generalized mode-coupling theory (GMCT) has recently emerged as a promising first-principles theory to study the poorly understood dynamics of glass-forming materials. Formulated as a hierarchical extension of standard mode-coupling theory (MCT), it is able to systematically improve its predictions by including the exact dynamics of higher-order correlation functions into its hierarchy. However, in contrast to Newtonian dynamics, a fully generalized version of the theory based on Brownian dynamics is still lacking. To close this gap, we provide a detailed derivation of GMCT for colloidal mixtures obeying a many-body Smoluchowski equation. We demonstrate that a hierarchy of coupled equations can again be established and show that these, consistent with standard MCT, are identical to the ones obtained from Newtonian GMCT when taking the overdamped limit. Consequently, the non-trivial similarity between Brownian and Newtonian MCT is maintained for our newly developed multi-component GMCT. As a proof of principle, we also solve the generalized mode-coupling equations for the binary Kob-Andersen Lennard-Jones mixture undergoing Brownian dynamics, and confirm the improved predictive power of the theory upon using more levels of the GMCT hierarchy of equations.
The short-time motion of Brownian particles in an incompressible Newtonian fluid under shear, in which the fluid inertia becomes important, was investigated by direct numerical simulation of particulate flows. Three-dimensional simulations were performed, wherein external forces were introduced to approximately form Couette flows throughout the entire system with periodic boundary conditions. In order to examine the validity of the method, the mean square displacement of a single spherical particle in a simple shear flow was calculated, and these results were compared with a hydrodynamic analytical solution that includes the effects of the fluid inertia. Finally, the dynamical behavior of a monodisperse dispersion composed of repulsive spherical particles was examined on short time scales, and the shear-induced diffusion coefficients were measured for several volume fractions up to 0.50.