Consistent Temperature Coupling with Thermal Fluctuations of Smooth Particle Hydrodynamics and Molecular Dynamics

Abstract

We propose a thermodynamically consistent and energy-conserving temperature coupling scheme between the atomistic and the continuum domain. The coupling scheme links the two domains using the DPDE (Dissipative Particle Dynamics at constant Energy) thermostat and is designed to handle strong temperature gradients across the atomistic/continuum domain interface. The fundamentally different definitions of temperature in the continuum and atomistic domain – internal energy and heat capacity versus particle velocity – are accounted for in a straightforward and conceptually intuitive way by the DPDE thermostat. We verify the here-proposed scheme using a fluid, which is simultaneously represented as a continuum using Smooth Particle Hydrodynamics, and as an atomistically resolved liquid using Molecular Dynamics. In the case of equilibrium contact between both domains, we show that the correct microscopic equilibrium properties of the atomistic fluid are obtained. As an example of a strong non-equilibrium situation, we consider the propagation of a steady shock-wave from the continuum domain into the atomistic domain, and show that the coupling scheme conserves both energy and shock-wave dynamics. To demonstrate the applicability of our scheme to real systems, we consider shock loading of a phospholipid bilayer immersed in water in a multi-scale simulation, an interesting topic of biological relevance.

Effect of Unsaturated Flow Modes on Partitioning Dynamics of Gravity‐Driven Flow at a Simple Fracture Intersection: Laboratory Study and Three‐Dimensional Smoothed Particle Hydrodynamics Simulations

“…To numerically study the given system, we employ a massively parallel SPH code implemented within the LAMMPS (large‐scale atomic/molecular massively parallel simulator; Ganzenmüller et al, ; Plimpton, ; Plimpton et al, ) framework. In its original form, LAMMPS has been developed for MD simulations.…”

Shock Wave-Induced Damage of a Protein by Void Collapse

“…MD simulations are routinely used to study atomic-level interactions in exquisite detail (15) and have proven to be an effective tool for aiding in the interpretation of experimental measurements. In particular, MD simulations have been previously used to study the effects of shock waves on membranes (9,10,(16)(17)(18)(19)(20)(21)(22). Many recent simulations (9,10,18) used the momentum-reflecting mirror method (see below) to imitate the motion of a piston to produce shock waves in the system under study.…”

Atomistic-Continuum Hybrid Simulation of Heat Transfer Between Argon Flow and Copper Plates

“…In fact, several effective hybrid simulation methods have been developed to study these particular phenomena caused by the size effect; these methods include atomic finite element method (AFEM) [5], atomistic-smooth particle method [6], and atomistic-finite volume method [7]. The AFEM has advantage of high computational efficiency for solid state problems, and the smooth-particle method is a simulation technique that is still under development [8], which also suffers issues from unclear physical meaning at boundary.…”