Stable schemes for dissipative particle dynamics with conserved energy
Abstract
This article presents a new numerical scheme for the discretization of dissipative particle dynamics with conserved energy. The key idea is to reduce elementary pairwise stochastic dynamics (either fluctuation/dissipation or thermal conduction) to effective singlevariable dynamics, and to approximate the solution of these dynamics with one step of a Metropolis–Hastings algorithm. This ensures by construction that no negative internal energies are encountered during the simulation, and hence allows to increase the admissible timesteps to integrate the dynamics, even for systems with small heat capacities. Stability is only limited by the Hamiltonian part of the dynamics, which suggests resorting to multiple timestep strategies where the stochastic part is integrated less frequently than the Hamiltonian one.
 Authors:
 Publication Date:
 OSTI Identifier:
 22622304
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 340; Other Information: Copyright (c) 2017 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ALGORITHMS; APPROXIMATIONS; COMPUTERIZED SIMULATION; DYNAMICS; FLUCTUATIONS; HAMILTONIANS; HEAT; MATHEMATICAL SOLUTIONS; PARTICLES; SPECIFIC HEAT; STABILITY; STOCHASTIC PROCESSES; THERMAL CONDUCTION
Citation Formats
Stoltz, Gabriel, Email: stoltz@cermics.enpc.fr. Stable schemes for dissipative particle dynamics with conserved energy. United States: N. p., 2017.
Web. doi:10.1016/J.JCP.2017.03.059.
Stoltz, Gabriel, Email: stoltz@cermics.enpc.fr. Stable schemes for dissipative particle dynamics with conserved energy. United States. doi:10.1016/J.JCP.2017.03.059.
Stoltz, Gabriel, Email: stoltz@cermics.enpc.fr. 2017.
"Stable schemes for dissipative particle dynamics with conserved energy". United States.
doi:10.1016/J.JCP.2017.03.059.
@article{osti_22622304,
title = {Stable schemes for dissipative particle dynamics with conserved energy},
author = {Stoltz, Gabriel, Email: stoltz@cermics.enpc.fr},
abstractNote = {This article presents a new numerical scheme for the discretization of dissipative particle dynamics with conserved energy. The key idea is to reduce elementary pairwise stochastic dynamics (either fluctuation/dissipation or thermal conduction) to effective singlevariable dynamics, and to approximate the solution of these dynamics with one step of a Metropolis–Hastings algorithm. This ensures by construction that no negative internal energies are encountered during the simulation, and hence allows to increase the admissible timesteps to integrate the dynamics, even for systems with small heat capacities. Stability is only limited by the Hamiltonian part of the dynamics, which suggests resorting to multiple timestep strategies where the stochastic part is integrated less frequently than the Hamiltonian one.},
doi = {10.1016/J.JCP.2017.03.059},
journal = {Journal of Computational Physics},
number = ,
volume = 340,
place = {United States},
year = 2017,
month = 7
}

The dynamic properties of fluid, including diffusivity and viscosity, are temperaturedependent and can significantly influence the flow dynamics of mesoscopic nonisothermal systems. To capture the correct temperaturedependence of a fluid, an energyconserving dissipative particle dynamics (eDPD) model is developed by expressing the weighting terms of the dissipative force and the random force as functions of temperature. The diffusivity and viscosity of liquid water at various temperatures ranging from 273 K to 373 K are used as examples for verifying the proposed model. Simulations of a Poiseuille flow and a steady case of heat conduction for reproducing the Fourier law aremore »

Dissipative particle dynamics with attractive and repulsive particleparticle interactions
In molecular dynamics simulations, a combination of shortrange repulsive and longrange attractive interactions allows the behavior of gases, liquids, solids, and multiphase systems to be simulated. We demonstrate that dissipative particle dynamics (DPD) simulations with similar pairwise particleparticle interactions can also be used to simulate the dynamics of multiphase fluids. In these simulations, the positive, shortrange, repulsive part of the interaction potentials were represented by polynomial spline functions such as those used as smoothing functions in smoothed particle hydrodynamics, and the negative longrange part of the interaction has the same form but a different range and amplitude. If a singlemore » 
Pseudorandom number generation for Brownian Dynamics and Dissipative Particle Dynamics simulations on GPU devices
Highlights: {yields} Molecular Dynamics codes implemented on GPUs have achieved twoorder of magnitude computational accelerations. {yields} Brownian Dynamics and Dissipative Particle Dynamics simulations require a large number of random numbers per time step. {yields} We introduce a method for generating small batches of pseudorandom numbers distributed over many threads of calculations. {yields} With this method, Dissipative Particle Dynamics is implemented on a GPU device without requiring threadtothread communication.  Abstract: Brownian Dynamics (BD), also known as Langevin Dynamics, and Dissipative Particle Dynamics (DPD) are implicit solvent methods commonly used in models of soft matter and biomolecular systems. The interaction ofmore » 
Dissipative Particle Dynamics and Other Particle Methods for Multiphase Fluid Flow in Fractured and Porous Media
Particle methods are less computationally efficient than grid based numerical solution of the Navier Stokes equation. However, they have important advantages including rigorous mass conservation, momentum conservation and isotropy. In addition, there is no need for explicit interface tracking/capturing and code development effort is relatively low. We describe applications of three particle methods: molecular dynamics, dissipative particle dynamics and smoothed particle hydrodynamics. The mesoscale (between the molecular and continuum scales) dissipative particle dynamics method can be used to simulate systems that are too large to simulate using molecular dynamics but small enough for thermal fluctuations to play an important role. 
A Phase Field Approach to No Slip Boundary Conditions in Dissipative Particle Dynamics and Other Particle Models for Fluid Flow in Geometrically Complex Confined Systems
Dissipative particle dynamics (DPD) is an effective mesoscopic particle model with a lower computational cost than molecular dynamics because of the soft potentials that it employs. However, the soft potential is not strong enough to prevent the fluid DPD particles from penetrating solid boundaries represented by stationary DPD particles. A phase field variable, _(x,t) , is used to indicate the phase at point x and time t, with a smooth transition from 1 (phase 1) to +1 (phase 2) across the interface. We describe an efficient implementation of noslip boundary conditions in DPD models that combine solidliquid particleparticle interactions withmore »