Concurrent coupling of atomistic simulation and mesoscopic hydrodynamics for flows over soft multi-functional surfaces

Abstract

We develop an efficient parallel multiscale method that bridges the atomistic and mesoscale regimes, from nanometer to micron and beyond, via concurrent coupling of atomistic simulation and mesoscopic dynamics. In particular, we combine an all-atom molecular dynamics (MD) description for specific atomistic details in the vicinity of the functional surface, with a dissipative particle dynamics (DPD) approach that captures mesoscopic hydrodynamics in the domain away from the functional surface. In order to achieve a seamless transition in dynamic properties we endow the MD simulation with a DPD thermostat, which is validated against experimental results by modeling water at different temperatures. We then validate the MD-DPD coupling method for transient Couette and Poiseuille flows, demonstrating that the concurrent MD-DPD coupling can resolve accurately the continuum-based analytical solutions. Subsequently, we simulate shear flows over polydimethylsiloxane (PDMS)-grafted surfaces (polymer brushes) for various grafting densities, and investigate the slip flow as a function of the shear stress. We verify that a "universal" power law exists for the sliplength, in agreement with published results. Having validated the MD-DPD coupling method, we simulate time-dependent flows past an endothelial glycocalyx layer (EGL) in a microchannel. Coupled simulation results elucidate the dynamics of EGL changing from an equilibrium state to a compressed state under shear by aligning the molecular structures along the shear direction. MD-DPD simulation results agree well with results of a single MD simulation, but with the former more than two orders of magnitude faster than the latter for system sizes above one micron. * The first two authors contribute equally to this work.

Understanding the Role of Endothelial Glycocalyx in Mechanotransduction via Computational Simulation: A Mini Review

“…A truly multiscale numerical method requires the dynamic coupling of different simulation methods in one simulation (Fedosov and Karniadakis, 2009). For example, a coupling between DPD and MD via a special multi-scale interface server (Wang et al, 2019) could solve the fluid problem and force transmission simultaneously. The challenge, however, is to find a suitable mapping between high-and low-resolution simulations in both spatial and temporal senses.…”

Particle–Continuum Coupling and its Scaling Regimes: Theory and Applications

“…In the first one, combining AdResS and a flux‐exchange hybrid MD/CFD method, [ 155 ] insertion of complex molecules to match the mass and momentum fluxes at the particle and continuum interface is made feasible by implementing there a low resolution model (blob‐molecules with soft effective interactions) and then using the AdResS to reintroduce the atomistic degrees of freedom and further couple with the bulk MD. [ 156,157 ] In the second one, coined the Triple‐decker, [ 107,172 ] the state‐variable coupling approach with the “handshaking” of the particle and continuum physical descriptions is used instead to couple atomistic and mesoscopic hydrodynamics.…”

CPL library — A minimal framework for coupled particle and continuum simulation

“…iii) Framework SediFOAM [45] CPL LIBRARY [49][50][51] MUSCLE [52] LIGGHTS/CFDEM [46] MUI [26] OASIS3 [53] OpenFOAM/LAMMPS [44] MaMiCo [47] MCT [54] Precice [55] Figure 2: [47] for Lattice Boltzmann-MD coupling and MUI [26] which allows general coupling of nonuniform grids and particle-particle coupling such as linking dissipative particle dynamics (DPD) and MD [48] . The Final classification considered are frameworks, codes which provide a series of tools and techniques as an environment for general coupling of two or more packages.…”