Abstract
We study numerical methods for dissipative particle dynamics (DPD), which is a system of stochastic differential equations and a popular stochastic momentum-conserving thermostat for simulating complex hydrodynamic behavior at mesoscales. We propose a new splitting method that is able to substantially improve the accuracy and efficiency of DPD simulations in a wide range of the friction coefficients, particularly in the extremely large friction limit that corresponds to a fluid-like Schmidt number, a key issue in DPD. Various numerical experiments on both equilibrium and transport properties are performed to demonstrate the superiority of the newly proposed method over popular alternative schemes in the literature.
Original language | English |
---|---|
Article number | A1929–A1949 |
Pages (from-to) | A1929-A1949 |
Number of pages | 21 |
Journal | SIAM Journal on Scientific Computing |
Volume | 43 |
Issue number | 3 |
DOIs | |
Publication status | Published - 27 May 2021 |
Keywords
- Dissipative particle dynamics
- Invariant measure
- Order of convergence
- Splitting methods
- Stochastic differential equations
- Transport properties
ASJC Scopus subject areas
- Computational Mathematics
- Applied Mathematics