Technical Brief

An Efficient Dissipative Particle Dynamics-Based Algorithm for Simulating Ferromagnetic Colloidal Suspensions

[+] Author and Article Information
Wuming Li, Qingsheng Liu

Department of Applied Mathematics,
Northwestern Polytechnical University,
Xi'an 710129, China

Jie Ouyang

Department of Applied Mathematics,
Northwestern Polytechnical University,
Xi'an 710129, China
e-mail: jieouyang@nwpu.edu.cn

1Corresponding author.

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received December 27, 2014; final manuscript received March 9, 2015; published online August 26, 2015. Assoc. Editor: Paramsothy Jayakumar.

J. Comput. Nonlinear Dynam 11(2), 024501 (Aug 26, 2015) (6 pages) Paper No: CND-14-1333; doi: 10.1115/1.4030155 History: Received December 27, 2014

In this paper, the algorithm, Euler scheme-the modified velocity-verlet algorithm (ES-MVVA) based on dissipative particle dynamics (DPD) method, is applied to simulate a two-dimensional ferromagnetic colloidal suspension. The very desirable aggregate structures of magnetic particles are obtained by using the above-mentioned algorithm, which are in qualitatively good agreement with those in the literature obtained by other simulation methods for different magnetic particle–particle interaction strengths. At the same time, the radial distribution functions of magnetic particles and the mean equilibrium temperatures of the system are also calculated. Next, the mean equilibrium velocities of magnetic and dissipative particles are calculated, by comparing the results obtained by ES-MVVA with those obtained by other algorithm for different time step sizes, it shows the validity and good accuracy of the present algorithm. So, the DPD-based algorithm presented in this paper is a powerful tool for simulation of magnetic colloidal suspensions.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Fig. 1

Model of magnetic particle

Grahic Jump Location
Fig. 2

Model of interaction between magnetic and dissipative particles [26]

Grahic Jump Location
Fig. 3

Influence of dissipative particles mass m* and magnetic interaction strength λm on aggregate structures for dc*=0.4 based on algorithm of ES-MVVA: (a) for m*=0.05 and λm=10, (b) for m*=0.05 and λm=3, (c) for m*=0.01 and λm=10, (d) for m*=0.01 and λm=3, (e) for m*=0.005 and λm=10, and (f) for m*=0.005 and λm=3

Grahic Jump Location
Fig. 4

Radial distribution function for m*=0.01 by ES-MVVA (a) and MVVA–MVVA (b)




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In