Brownian and Stokesian Dynamics

Author(s):  
G. Bossis ◽  
J. F. Brady
Keyword(s):  
Soft Matter ◽  
2021 ◽  
Author(s):  
Vincent Labalette ◽  
Alexis Praga ◽  
Florent Girard ◽  
Martine Meireles ◽  
Yannick Hallez ◽  
...  

A new numerical framework based on Stokesian dynamics is used to study a shear-induced glass-to-crystal transition in suspensions of clay-like anisotropically charged platelets.


2008 ◽  
Vol 99 (01) ◽  
pp. 108-115 ◽  
Author(s):  
Koichiro Yano ◽  
Ken-ichi Tsubota ◽  
Takuji Ishikawa ◽  
Shigeo Wada ◽  
Takami Yamaguchi ◽  
...  

SummaryWe propose a method to analyze platelet adhesion and aggregation computationally, taking into account the distinct properties of two plasma proteins, vonWillebrand factor (vWF) and fibrinogen (Fbg). In this method, the hydrodynamic interactions between platelet particles under simple shear flow were simulated using Stokesian dynamics based on the additivity of velocities. The binding force between particles mediated by vWF and Fbg was modeled using the Voigt model. Two Voigt models with different properties were introduced to consider the distinct behaviors of vWF and Fbg. Our results qualitatively agreed with the general observation of a previous in-vitro experiment, thus demonstrating that the significant development of thrombus formation in height requires not only vWF, but also Fbg. This agreement of simulation and experimental results qualitatively validates our model and suggests that consideration of the distinct roles of vWF and Fbg is essential to investigate the physiological and pathophysiological mechanisms of thrombus formation using a computational approach.


2015 ◽  
Vol 769 ◽  
pp. 369-386 ◽  
Author(s):  
A. Lefebvre-Lepot ◽  
B. Merlet ◽  
T. N. Nguyen

We address the problem of computing the hydrodynamic forces and torques among $N$ solid spherical particles moving with given rotational and translational velocities in Stokes flow. We consider the original fluid–particle model without introducing new hypotheses or models. Our method includes the singular lubrication interactions which may occur when some particles come close to one another. The main new feature is that short-range interactions are propagated to the whole flow, including accurately the many-body lubrication interactions. The method builds on a pre-existing fluid solver and is flexible with respect to the choice of this solver. The error is the error generated by the fluid solver when computing non-singular flows (i.e. with negligible short-range interactions). Therefore, only a small number of degrees of freedom are required and we obtain very accurate simulations within a reasonable computational cost. Our method is closely related to a method proposed by Sangani & Mo (Phys. Fluids, vol. 6, 1994, pp. 1653–1662) but, in contrast with the latter, it does not require parameter tuning. We compare our method with the Stokesian dynamics of Durlofsky et al. (J. Fluid Mech., vol. 180, 1987, pp. 21–49) and show the higher accuracy of the former (both by analysis and by numerical experiments).


2021 ◽  
Vol 925 ◽  
Author(s):  
Emma Gonzalez ◽  
Christian Aponte-Rivera ◽  
Roseanna N. Zia

We present a computational study of the equilibrium dynamics of a polydisperse hard-sphere colloidal dispersion confined in a spherical cavity. We account for many-body hydrodynamic and lubrication interactions between particles and with the confining cavity utilizing our confined Stokesian dynamics model, expanded here for size polydispersity. We find that, even though the tendency of polydispersity to homogenize structure in a suspension is still present in confinement, strong correlations induced by the cavity resist homogenization. Although seemingly opposite, these two effects have a common driver, which is to maximize configurational entropy of particles in the cavity interior. These structural effects couple with the hydrodynamics to change the particle dynamics: polydispersity weakens lubrication effects near the cavity wall, allowing small (large) particles to diffuse faster (slower) than in a monodisperse suspension. As a small (large) particle gets farther from the wall, polydispersity weakens many-body hydrodynamic couplings, driving diffusivity up (down). While the local cage dynamics dominates short-time self-diffusion, long-time dynamics is also affected. In the concentrated regime, polydispersity and confinement combine to induce radial de-mixing into size-segregated populations. The cavity becomes the most influential ‘nearest neighbour’, setting the length scale of and dynamics within these radial domains. This intermediate length-scale caging makes the angular dynamics insensitive to polydispersity but leads to radial long-time mean-square displacement that changes qualitatively with volume composition. These results hold promise for explaining colloidal-scale physics implicated in the functioning of biological cells, and the engineering of non-living confined colloids where size de-mixing could be useful in the design of encapsulated micro-reactors and therapeutic vesicles.


1994 ◽  
Vol 275 ◽  
pp. 157-199 ◽  
Author(s):  
Prabhu R. Nott ◽  
John F. Brady

Dynamic simulations of the pressure-driven flow in a channel of a non-Brownian suspension at zero Reynolds number were conducted using Stokesian Dynamics. The simulations are for a monolayer of identical particles as a function of the dimensionless channel width and the bulk particle concentration. Starting from a homogeneous dispersion, the particles gradually migrate towards the centre of the channel, resulting in an homogeneous concentration profile and a blunting of the particle velocity profile. The time for achieving steady state scales as (H/a)3a/〈u〉, where H is the channel width, a the radii of the particles, and 〈u〉 the average suspension velocity in the channel. The concentration and velocity profiles determined from the simulations are in qualitative agreement with experiment.A model for suspension flow has been proposed in which macroscopic mass, momentum and energy balances are constructed and solved simultaneously. It is shown that the requirement that the suspension pressure be constant in directions perpendicular to the mean motion leads to particle migration and concentration variations in inhomogeneous flow. The concept of the suspension ‘temperature’ – a measure of the particle velocity fluctuations – is introduced in order to provide a nonlocal description of suspension behaviour. The results of this model for channel flow are in good agreement with the simulations.


2016 ◽  
Vol 10 (6) ◽  
pp. 064117
Author(s):  
Yousef M. F. El Hasadi ◽  
Martin Crapper

1988 ◽  
Vol 20 (1) ◽  
pp. 111-157 ◽  
Author(s):  
J F Brady ◽  
G Bossis
Keyword(s):  

2000 ◽  
Vol 412 ◽  
pp. 279-301 ◽  
Author(s):  
ANUGRAH SINGH ◽  
PRABHU R. NOTT

We report the normal stresses in a non-Brownian suspension in plane Couette flow determined from Stokesian Dynamics simulations. The presence of normal stresses that are linear in the shear rate in a viscometric flow indicates a non-Newtonian character of the suspension, which is otherwise Newtonian. While in itself of interest, this phenomenon is also important because it is believed that normal stresses determine the migration of particles in flows with inhomogeneous shear fields. We simulate plane Couette flow by placing a layer of clear fluid adjacent to one wall in the master cell, which is then replicated periodically. From a combination of the traceless hydrodynamic stresslet on the suspended particles, the stresslet due to (non-hydrodynamic) inter-particle forces, and the total normal force on the walls, we determine the hydrodynamic and inter-particle force contributions to the isotropic ‘particle pressure’ and the first normal stress difference. We determine the stresses for a range of the particle concentration and the Couette gap. The particle pressure and the first normal stress difference exhibit a monotonic increase with the mean particle volume fraction ϕ. The ratio of normal to shear stresses on the walls also increases with ϕ, substantiating the result of Nott & Brady (1994) that this condition is required for stability to concentration fluctuations. We also study the microstructure by extracting the pair distribution function from our simulations; our results are in agreement with previous studies showing anisotropy in the pair distribution, which is the cause of normal stresses.


2002 ◽  
Vol 456 ◽  
pp. 239-275 ◽  
Author(s):  
J. BERGENHOLTZ ◽  
J. F. BRADY ◽  
M. VICIC

The non-Newtonian rheology is calculated numerically to second order in the volume fraction in steady simple shear flows for Brownian hard spheres in the presence of hydrodynamic and excluded volume interactions. Previous analytical and numerical results for the low-shear structure and rheology are confirmed, demonstrating that the viscosity shear thins proportional to Pe2, where Pe is the dimensionless shear rate or Péclet number, owing to the decreasing contribution of Brownian forces to the viscosity. In the large Pe limit, remnants of Brownian diffusion balance convection in a boundary-layer in the compressive region of the flow. In consequence, the viscosity shear thickens when this boundary-layer coincides with the near-contact lubrication regime of the hydrodynamic interaction. Wakes are formed at large Pe in the extensional zone downstream from the reference particle, leading to broken symmetry in the pair correlation function. As a result of this asymmetry and that in the boundary-layer, finite normal stress differences are obtained as well as positive departures in the generalized osmotic pressure from its equilibrium value. The first normal stress difference changes from positive to negative values as Pe is increased when the hard-sphere limit is approached. This unusual effect is caused by the hydrodynamic lubrication forces that maintain particles in close proximity well into the extensional quadrant of the flow. The study demonstrates that many of the non-Newtonian effects observed in concentrated suspensions by experiments and by Stokesian dynamics simulations are present also in dilute suspensions.


Author(s):  
Marco Vanni

The stresses acting on aggregates smaller than the Kolmogorov length scale in homogeneous isotropic turbulence were estimated by a two-scale numerical simulation. The fluid dynamics at the scales larger than the Kolmogorov length scale was calculated by a Direct Numerical Simulation of the turbulent flow, in which the aggregates were modeled as point particles. Then, we adopted Stokesian Dynamics to evaluate the phenomena governed by the smooth velocity field of the smallest scales. At this level the disordered structure of the aggregates was modeled in detail, in order to take into account the role that the primary particles have in generating and transferring the internal stress. From this result, it was possible to evaluate the internal forces acting at intermonomer contacts and determine the occurrence of breakup as a consequence of the failure of intermonomer bonds. The method was applied to disordered aggregates with isostatic and highly hyperstatic structures, respectively.


Sign in / Sign up

Export Citation Format

Share Document