Density functional theory of phase coexistence in weakly polydisperse fluids

The recently proposed universal relations between the moments of the polydispersity distributions of a phase-separated weakly polydisperse system are analyzed in detail using the numerical results obtained by solving a simple density functional theory of a polydisperse fluid. It is shown that universal properties are the exception rather than the rule.Comment: 10 pages, 2 figures, to appear in PR

Relaxation of surface charge on rotating dielectric spheres: Implications on dynamic electrorheological effects

We have examined the effect of an oscillatory rotation of a polarized dielectric particle. The rotational motion leads to a re-distribution of the polarization charge on the surface of the particle. We show that the time averaged steady-state dipole moment is along the field direction, but its magnitude is reduced by a factor which depends on the angular velocity of rotation. As a result, the rotational motion of the particle reduces the electrorheological effect. We further assume that the relaxation of polarized charge is arised from a finite conductivity of the particle or host medium. We calculate the relaxation time based on the Maxwell-Wagner theory, suitably generalized to include the rotational motion. Analytic expressions for the reduction factor and the relaxation time are given and their dependence on the angular velocity of rotation will be discussed.Comment: Accepted for publications by Phys. Rev.

Generation of Porous Particle Structures using the Void Expansion Method

The newly developed "void expansion method" allows for an efficient generation of porous packings of spherical particles over a wide range of volume fractions using the discrete element method. Particles are randomly placed under addition of much smaller "void-particles". Then, the void-particle radius is increased repeatedly, thereby rearranging the structural particles until formation of a dense particle packing. The structural particles' mean coordination number was used to characterize the evolving microstructures. At some void radius, a transition from an initially low to a higher mean coordination number is found, which was used to characterize the influence of the various simulation parameters. For structural and void-particle stiffnesses of the same order of magnitude, the transition is found at constant total volume fraction slightly below the random close packing limit. For decreasing void-particle stiffness the transition is shifted towards a smaller void-particle radius and becomes smoother.Comment: 9 pages, 8 figure

Strong Phase Separation in a Model of Sedimenting Lattices

We study the steady state resulting from instabilities in crystals driven through a dissipative medium, for instance, a colloidal crystal which is steadily sedimenting through a viscous fluid. The problem involves two coupled fields, the density and the tilt; the latter describes the orientation of the mass tensor with respect to the driving field. We map the problem to a 1-d lattice model with two coupled species of spins evolving through conserved dynamics. In the steady state of this model each of the two species shows macroscopic phase separation. This phase separation is robust and survives at all temperatures or noise levels--- hence the term Strong Phase Separation. This sort of phase separation can be understood in terms of barriers to remixing which grow with system size and result in a logarithmically slow approach to the steady state. In a particular symmetric limit, it is shown that the condition of detailed balance holds with a Hamiltonian which has infinite-ranged interactions, even though the initial model has only local dynamics. The long-ranged character of the interactions is responsible for phase separation, and for the fact that it persists at all temperatures. Possible experimental tests of the phenomenon are discussed.Comment: To appear in Phys Rev E (1 January 2000), 16 pages, RevTex, uses epsf, three ps figure

The Influence of the Degree of Heterogeneity on the Elastic Properties of Random Sphere Packings

The macroscopic mechanical properties of colloidal particle gels strongly depend on the local arrangement of the powder particles. Experiments have shown that more heterogeneous microstructures exhibit up to one order of magnitude higher elastic properties than their more homogeneous counterparts at equal volume fraction. In this paper, packings of spherical particles are used as model structures to computationally investigate the elastic properties of coagulated particle gels as a function of their degree of heterogeneity. The discrete element model comprises a linear elastic contact law, particle bonding and damping. The simulation parameters were calibrated using a homogeneous and a heterogeneous microstructure originating from earlier Brownian dynamics simulations. A systematic study of the elastic properties as a function of the degree of heterogeneity was performed using two sets of microstructures obtained from Brownian dynamics simulation and from the void expansion method. Both sets cover a broad and to a large extent overlapping range of degrees of heterogeneity. The simulations have shown that the elastic properties as a function of the degree of heterogeneity are independent of the structure generation algorithm and that the relation between the shear modulus and the degree of heterogeneity can be well described by a power law. This suggests the presence of a critical degree of heterogeneity and, therefore, a phase transition between a phase with finite and one with zero elastic properties.Comment: 8 pages, 6 figures; Granular Matter (published online: 11. February 2012

Simulating (electro)hydrodynamic effects in colloidal dispersions: smoothed profile method

Previously, we have proposed a direct simulation scheme for colloidal dispersions in a Newtonian solvent [Phys.Rev.E 71,036707 (2005)]. An improved formulation called the ``Smoothed Profile (SP) method'' is presented here in which simultaneous time-marching is used for the host fluid and colloids. The SP method is a direct numerical simulation of particulate flows and provides a coupling scheme between the continuum fluid dynamics and rigid-body dynamics through utilization of a smoothed profile for the colloidal particles. Moreover, the improved formulation includes an extension to incorporate multi-component fluids, allowing systems such as charged colloids in electrolyte solutions to be studied. The dynamics of the colloidal dispersions are solved with the same computational cost as required for solving non-particulate flows. Numerical results which assess the hydrodynamic interactions of colloidal dispersions are presented to validate the SP method. The SP method is not restricted to particular constitutive models of the host fluids and can hence be applied to colloidal dispersions in complex fluids

A Dissipative-Particle-Dynamics Model for Simulating Dynamics of Charged Colloid

A mesoscopic colloid model is developed in which a spherical colloid is represented by many interacting sites on its surface. The hydrodynamic interactions with thermal fluctuations are taken accounts in full using Dissipative Particle Dynamics, and the electrostatic interactions are simulated using Particle-Particle-Particle Mesh method. This new model is applied to investigate the electrophoretic mobility of a charged colloid under an external electric field, and the influence of salt concentration and colloid charge are systematically studied. The simulation results show good agreement with predictions from the electrokinetic theory.Comment: 17 pages, 8 figures, submitted to the proceedings of High Performance Computing in Science & Engineering '1

Aspects of the dynamics of colloidal suspensions: Further results of the mode-coupling theory of structural relaxation

Results of the idealized mode-coupling theory for the structural relaxation in suspensions of hard-sphere colloidal particles are presented and discussed with regard to recent light scattering experiments. The structural relaxation becomes non-diffusive for long times, contrary to the expectation based on the de Gennes narrowing concept. A semi-quantitative connection of the wave vector dependences of the relaxation times and amplitudes of the final $\alpha$-relaxation explains the approximate scaling observed by Segr{\`e} and Pusey [Phys. Rev. Lett. {\bf 77}, 771 (1996)]. Asymptotic expansions lead to a qualitative understanding of density dependences in generalized Stokes-Einstein relations. This relation is also generalized to non-zero frequencies thereby yielding support for a reasoning by Mason and Weitz [Phys. Rev. Lett {\bf 74}, 1250 (1995)]. The dynamics transient to the structural relaxation is discussed with models incorporating short-time diffusion and hydrodynamic interactions for short times.Comment: 11 pages, 9 figures; to be published in Phys. Rev.