Translational and rotational friction on a colloidal rod near a wall

We present particulate simulation results for translational and rotational friction components of a shish-kebab model of a colloidal rod with aspect ratio (length over diameter) $L/D = 10$ in the presence of a planar hard wall. Hydrodynamic interactions between rod and wall cause an overall enhancement of the friction tensor components. We find that the friction enhancements to reasonable approximation scale inversely linear with the closest distance $d$ between the rod surface and the wall, for $d$ in the range between $D/8$ and $L$. The dependence of the wall-induced friction on the angle $\theta$ between the long axis of the rod and the normal to the wall is studied and fitted with simple polynomials in $\cos \theta$.Comment: 8 pages, 8 figure

Mesoscale modeling of the rheology of pressure sensitive adhesives through inclusion of transient forces

For optimal application, pressure-sensitive adhesives must have rheological properties in between those of a viscoplastic solid and those of a viscoelastic liquid. Such adhesives can be produced by emulsion polymerisation, resulting in latex particles which are dispersed in water and contain long-chain acrylic polymers. When the emulsion is dried, the latex particles coalesce and an adhesive film is formed. The rheological properties of the dried samples are believed to be dominated by the interface regions between the original latex particles, but the relation between rheology and latex particle properties is poorly understood. In this paper we show that it is possible to describe the bulk rheology of a pressure-sensitive adhesive by means of a mesoscale simulation model. To reach experimental time and length scales, each latex particle is represented by just one simulated particle. The model is subjected to oscillatory shear flow and extensional flow. Simple order of magnitude estimates of the model parameters already lead to semi-quantitative agreement with experimental results. We show that inclusion of transient forces in the model, i.e. forces with memory of previous configurations, is essential to correctly predict the linear and nonlinear properties.Comment: 29 pages, 8 figure

Force calculation on walls and embedded particles in multiparticle collision dynamics simulations

Colloidal solutions posses a wide range of time and length scales, so that it is unfeasible to keep track of all of them within a single simulation. As a consequence some form of coarse-graining must be applied. In this work we use the Multi-Particle Collision Dynamics scheme. We describe a particular implementation of no-slip boundary conditions upon a solid surface, capable of providing correct force s on the solid bypassing the calculation of the velocity profile or the stre ss tensor in the fluid near the surface. As an application we measure the friction on a spherical particle, when it is placed in a bulk fluid and when it is confined in a slit. We show that the implementation of the no-slip boundary conditions leads to an enhanced Ensko g friction, which can be understood analytically. Because of the long-range nature of hydrodynamic interactions, the Stokes friction obtained from the simulations is sensitive of the simulation box size. We address this topic for the slit geometry, showing that that the dependence on the system size differs very much from what is expected in a 3D system, where periodic boundary conditions are used in all directions.Comment: To appear in Physical Review

Efficient simulation of non-crossing fibers and chains in a hydrodynamic solvent

An efficient simulation method is presented for Brownian fiber suspensions, which includes both uncrossability of the fibers and hydrodynamic interactions between the fibers mediated by a mesoscopic solvent. To conserve hydrodynamics, collisions between the fibers are treated such that momentum and energy are conserved locally. The choice of simulation parameters is rationalised on the basis of dimensionless numbers expressing the relative strength of different physical processes. The method is applied to suspensions of semiflexible fibers with a contour length equal to the persistence length, and a mesh size to contour length ratio ranging from 0.055 to 0.32. For such fibers the effects of hydrodynamic interactions are observable, but relatively small. The non-crossing constraint, on the other hand, is very important and leads to hindered displacements of the fibers, with an effective tube diameter in agreement with recent theoretical predictions. The simulation technique opens the way to study the effect of viscous effects and hydrodynamic interactions in microrheology experiments where the response of an actively driven probe bead in a fiber suspension is measured.Comment: 12 pages, 2 tables, 5 figure

Transport coefficients of multi-particle collision algorithms with velocity-dependent collision rules

Detailed calculations of the transport coefficients of a recently introduced particle-based model for fluid dynamics with a non-ideal equation of state are presented. Excluded volume interactions are modeled by means of biased stochastic multiparticle collisions which depend on the local velocities and densities. Momentum and energy are exactly conserved locally. A general scheme to derive transport coefficients for such biased, velocity dependent collision rules is developed. Analytic expressions for the self-diffusion coefficient and the shear viscosity are obtained, and very good agreement is found with numerical results at small and large mean free paths. The viscosity turns out to be proportional to the square root of temperature, as in a real gas. In addition, the theoretical framework is applied to a two-component version of the model, and expressions for the viscosity and the difference in diffusion of the two species are given.Comment: 31 pages, 8 figures, accepted by J. Phys. Cond. Matte

Mechanical Properties of End-crosslinked Entangled Polymer Networks using Sliplink Brownian Dynamics Simulations

The mechanical properties of a polymeric network containing both crosslinks and sliplinks (entanglements) are studied using a multi-chain Brownian dynamics simulation. We coarse-grain at the level of chain segments connecting consecutive nodes (cross- or sliplinks), with particular attention to the Gaussian statistics of the network. Affine displacement of nodes is not imposed: their displacement as well as sliding of monomers through sliplinks is governed by force balances. The simulation results of stress in uniaxial extension and the full stress tensor in simple shear including the (non-zero) second normal stress difference are presented for monodisperse chains with up to 18 entanglements between two crosslinks. The cases of two different force laws of the subchains (Gaussian chains and chains with finite extensibility) for two different numbers of monomers in a subchain (no = 50 and no = 100) are examined. It is shown that the additivity assumption of slip- and crosslink contribution holds for sufficiently long chains with two or more entanglements, and that it can be used to construct the strain response of a network of infinitely long chains. An important consequence is that the contribution of sliplinks to the small-strain shear modulus is about ⅔ of the contribution of a crosslink

Systematic coarse-graining of the dynamics of entangled polymer melts: the road from chemistry to rheology

For optimal processing and design of entangled polymeric materials it is important to establish a rigorous link between the detailed molecular composition of the polymer and the viscoelastic properties of the macroscopic melt. We review current and past computer simulation techniques and critically assess their ability to provide such a link between chemistry and rheology. We distinguish between two classes of coarse-graining levels, which we term coarse-grained molecular dynamics (CGMD) and coarse-grained stochastic dynamics (CGSD). In CGMD the coarse-grained beads are still relatively hard, thus automatically preventing bond crossing. This also implies an upper limit on the number of atoms that can be lumped together and therefore on the longest chain lengths that can be studied. To reach a higher degree of coarse-graining, in CGSD many more atoms are lumped together, leading to relatively soft beads. In that case friction and stochastic forces dominate the interactions, and actions must be undertaken to prevent bond crossing. We also review alternative methods that make use of the tube model of polymer dynamics, by obtaining the entanglement characteristics through a primitive path analysis and by simulation of a primitive chain network. We finally review super-coarse-grained methods in which an entire polymer is represented by a single particle, and comment on ways to include memory effects and transient forces.Comment: Topical review, 31 pages, 10 figure

A Numerical Model for Brownian Particles Fluctuating in Incompressible Fluids

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.Comment: 6 pages, 8 figure

Velocity autocorrelation function of a Brownian particle

In this article, we present molecular dynamics study of the velocity autocorrelation function (VACF) of a Brownian particle. We compare the results of the simulation with the exact analytic predictions for a compressible fluid from [6] and an approximate result combining the predictions from hydrodynamics at short and long times. The physical quantities which determine the decay were determined from separate bulk simulations of the Lennard-Jones fluid at the same thermodynamic state point.We observe that the long-time regime of the VACF compares well the predictions from the macroscopic hydrodynamics, but the intermediate decay is sensitive to the viscoelastic nature of the solvent.Comment: 7 pages, 6 figure