### Abstract

Original language | English |
---|---|

Pages (from-to) | 2175-2187 |

Journal | Physical review E: Statistical physics, plasmas, fluids, and related interdisciplinary topics |

Volume | 59 |

Issue number | 2 |

DOIs | |

Publication status | Published - 1998 |

### Fingerprint

### Keywords

- METIS-140433

### Cite this

}

*Physical review E: Statistical physics, plasmas, fluids, and related interdisciplinary topics*, vol. 59, no. 2, pp. 2175-2187. https://doi.org/10.1103/PhysRevE.59.2175

**Brownian Dynamics Simulation of a Hard-Sphere Suspension.** / Strating, P.

Research output: Contribution to journal › Article › Academic › peer-review

TY - JOUR

T1 - Brownian Dynamics Simulation of a Hard-Sphere Suspension

AU - Strating, P.

PY - 1998

Y1 - 1998

N2 - In this paper we discuss the nonequilibrium shear viscosity of a suspension of hard spheres that is modeled by neglecting hydrodynamic interactions in a consistent way. The aim is to establish the true capabilities of this model in predicting the properties of real suspensions. A Brownian dynamics algorithm is used to simulate the movements of hard spheres immersed in a Newtonian solvent in a nonequilibrium steady shear flow. A new development is the treatment of the overlap of spheres as elastic collisions, to simulate the no-flux boundary condition on the surfaces of rigid particles. This algorithm is compared with other algorithms suggested in the literature, and is shown to be simple and accurate even for two spheres at close distance. This provides an algorithm that is very suitable for calculating the pair distribution function and especially its hard-sphere contact value, both in equilibrium and nonequilibrium simulations. The algorithm is used to study the nonequilibrium stationary shear flow in the low shear limit. The simulations correctly reproduce the exact low-density limit of the perturbation of the pair distribution function. The perturbation of the pair distribution function in shear flow can be extracted from the simulation data and used to compute the stationary shear viscosity for a system of diffusing hard spheres without hydrodynamic interactions. This yields a flow curve for this model system including the low shear limit. It is found that the model shear viscosity fails at intermediate and high shear rates as can be expected from the neglect of hydrodynamic interactions, but also in the low shear limit at small and moderate volume fractions.

AB - In this paper we discuss the nonequilibrium shear viscosity of a suspension of hard spheres that is modeled by neglecting hydrodynamic interactions in a consistent way. The aim is to establish the true capabilities of this model in predicting the properties of real suspensions. A Brownian dynamics algorithm is used to simulate the movements of hard spheres immersed in a Newtonian solvent in a nonequilibrium steady shear flow. A new development is the treatment of the overlap of spheres as elastic collisions, to simulate the no-flux boundary condition on the surfaces of rigid particles. This algorithm is compared with other algorithms suggested in the literature, and is shown to be simple and accurate even for two spheres at close distance. This provides an algorithm that is very suitable for calculating the pair distribution function and especially its hard-sphere contact value, both in equilibrium and nonequilibrium simulations. The algorithm is used to study the nonequilibrium stationary shear flow in the low shear limit. The simulations correctly reproduce the exact low-density limit of the perturbation of the pair distribution function. The perturbation of the pair distribution function in shear flow can be extracted from the simulation data and used to compute the stationary shear viscosity for a system of diffusing hard spheres without hydrodynamic interactions. This yields a flow curve for this model system including the low shear limit. It is found that the model shear viscosity fails at intermediate and high shear rates as can be expected from the neglect of hydrodynamic interactions, but also in the low shear limit at small and moderate volume fractions.

KW - METIS-140433

U2 - 10.1103/PhysRevE.59.2175

DO - 10.1103/PhysRevE.59.2175

M3 - Article

VL - 59

SP - 2175

EP - 2187

JO - Physical review E: covering statistical, nonlinear, biological, and soft matter physics

JF - Physical review E: covering statistical, nonlinear, biological, and soft matter physics

SN - 2470-0045

IS - 2

ER -