Finite-difference methods for simulation models incorporating non-conservative forces

We discuss algorithms applicable to the numerical solution of second-order ordinary differential equations by finite-differences. We make particular reference to the solution of the dissipative particle dynamics fluid model, and present extensive results comparing one of the algorithms discussed with the standard method of solution. These results show the successful modeling of phase separation and surface tension in a binary immiscible fluid mixture.Comment: 27 pages RevTeX, 9 figures, J. Chem. Phys. (in press

Spinodal decomposition of off-critical quenches with a viscous phase using dissipative particle dynamics in two and three spatial dimensions

We investigate the domain growth and phase separation of hydrodynamically-correct binary immiscible fluids of differing viscosity as a function of minority phase concentration in both two and three spatial dimensions using dissipative particle dynamics. We also examine the behavior of equal-viscosity fluids and compare our results to similar lattice-gas simulations in two dimensions.Comment: 34 pages (11 figures); accepted for publication in Phys. Rev.

Computer simulations of domain growth and phase separation in two-dimensional binary immiscible fluids using dissipative particle dynamics

We investigate the dynamical behavior of binary fluid systems in two dimensions using dissipative particle dynamics. We find that following a symmetric quench the domain size R(t) grows with time t according to two distinct algebraic laws R(t) = t^n: at early times n = 1/2, while for later times n = 2/3. Following an asymmetric quench we observe only n = 1/2, and if momentum conservation is violated we see n = 1/3 at early times. Bubble simulations confirm the existence of a finite surface tension and the validity of Laplace's law. Our results are compared with similar simulations which have been performed previously using molecular dynamics, lattice-gas and lattice-Boltzmann automata, and Langevin dynamics. We conclude that dissipative particle dynamics is a promising method for simulating fluid properties in such systems.Comment: RevTeX; 22 pages, 5 low-resolution figures. For full-resolution figures, connect to http://www.tcm.phy.cam.ac.uk/~ken21/tension/tension.htm

Using Dissipative Particle Dynamics To Model Binary Immiscible Fluids

This paper looks at domain growth and phase separation in two-dimensional, binary, immiscible fluids using a new simulation technique called dissipative particle dynamics (DPD). We will find that the method can handle domain growth both qualitatively and quantitatively, yielding the correct scaling exponents and displaying a surface tension satisfying Laplace's law