Coupled Navier Stokes-Molecular Dynamics Simulation using Iterative Operator-Splitting Methods

Abstract

In this paper, we contribute a multi-scale method based on an iterative operator splitting method, which takes into account the disparity of the macroand microscopic scales. We couple Navier Stokes and Molecular Dynamics equations, while taken into account their underlying scales. Combining relaxation methods and averaging techniques we can optimize the computational effort. The proposed method enables simulations in which the continuum flow aspects of the flow are described by the Navier-Stokes equations at timescales appropriate for this level of modeling, while the viscous stresses within the Navier-Stokes equations are the result of Molecular Dynamics Simulations, with much smaller time-scales. The main benefit of the proposed method is that time-dependent flows can then be modeled with a computational work which is significantly smaller than when the complete flow would be modeled at molecular level, as a result of the different time-scales at the continuum and molecular levels, enabled by the application of the iterative operator-splitting method. We discuss the convergence analysis to the splitting methods, see also [27]. Finally we present numerical results for the modified methods and applications to real-life flow problems.