Generalised dissipative particle dynamics with energy conservation: Density- And temperature-dependent potentials

Identifier:  imarina:6389712

Handle: https://hdl.handle.net/20.500.11797/imarina6389712

Authors:  Bonet Avalos, Josep; Lisal, Martin; Larentzos, James P; Mackie, Allan D; Brennan, John K

Abstract:
© 2019 the Owner Societies. We present a generalised, energy-conserving dissipative particle dynamics (DPDE) method appropriate for the non-isothermal simulation of particle interaction force fields that are both density- and temperature-dependent. A detailed derivation is formulated in a bottom-up manner by considering the thermodynamics of small systems with the appropriate consideration of the fluctuations. Connected to the local volume is a local density and corresponding local pressure, which is determined from an equation-of-state based force field that depends also on a particle temperature. Compared to the original DPDE method, the formulation of the generalised DPDE method requires a change in the independent variable from the particle internal energy to the particle entropy. As part of the re-formulation, the terms dressed particle entropy and the corresponding dressed particle temperature are introduced, which depict the many-body contributions in the local volume. The generalised DPDE method has similarities to the energy form of the smoothed dissipative particle dynamics method, yet fundamental differences exist, which are described in the manuscript. The basic dynamic equations are presented along with practical considerations for implementing the generalised DPDE method, including a numerical integration scheme based on the Shardlow-like splitting algorithm. Demonstrations and validation tests are performed using analytical equation-of-states for the van der Waals and Lennard-Jones fluids. Particle probability distributions are analysed, where excellent agreement with theoretical estimates is demonstrated. As further validation of the generalised DPDE method, both equilibrium and non-equilibrium simulation scenarios are considered, including adiabatic flash heating response and vapour-liquid phase separation.