Control of density fluctuations in atomistic-continuum simulations of dense liquids
Kotsalis E.M., Walther J.H., Koumoutsakos P., Physical Review E, 76, 016709, 2007
We present a control algorithm to eliminate spurious density fluctuations associated with the coupling of atomistic and continuum descriptions for dense liquids. A Schwartz domain decomposition algorithm is employed to couple molecular dynamics for the simulation of the atomistic system with a continuum solver for the simulation of the Navier-Stokes equations. The lack of periodic boundary conditions in the molecular dynamics simulations hinders the proper accounting for the virial pressure leading to spurious density fluctuations at the continuum-atomistic interface. An ad hoc boundary force is usually employed to remedy this situation. We propose the calculation of this boundary force using a control algorithm that explicitly cancels the density fluctuations. The results demonstrate that the present approach outperforms state-of-the-art algorithms. The conceptual and algorithmic simplicity of the method makes it suitable for any type of coupling between atomistic and continuum descriptions of dense fluids.