A combined Event-Driven/Time-Driven molecular dynamics algorithm for the simulation of shock waves in rarefied gases

A novel combined Event-Driven/Time-Driven (ED/TD) algorithm to speed-up the Molecular Dynamics simulation of rarefied gases using realistic spherically symmetric soft potentials is presented. Due to the low density regime, the proposed method correctly identifies the time that must elapse before the next interaction occurs, similarly to Event-Driven Molecular Dynamics. However, each interaction is treated using Time-Driven Molecular Dynamics, thereby integrating Newton's Second Law using the sufficiently small time step needed to correctly resolve the atomic motion. Although infrequent, many-body interactions are also accounted for with a small approximation. The combined ED/TD method is shown to correctly reproduce translational relaxation in argon, described using the Lennard-Jones potential. For densities between @r=10^-^4kg/m^3 and @r=10^-^1kg/m^3, comparisons with kinetic theory, Direct Simulation Monte Carlo, and pure Time-Driven Molecular Dynamics demonstrate that the ED/TD algorithm correctly reproduces the proper collision rates and the evolution toward thermal equilibrium. Finally, the combined ED/TD algorithm is applied to the simulation of a Mach 9 shock wave in rarefied argon. Density and temperature profiles as well as molecular velocity distributions accurately match DSMC results, and the shock thickness is within the experimental uncertainty. For the problems considered, the ED/TD algorithm ranged from several hundred to several thousand times faster than conventional Time-Driven MD. Moreover, the force calculation to integrate the molecular trajectories is found to contribute a negligible amount to the overall ED/TD simulation time. Therefore, this method could pave the way for the application of much more refined and expensive interatomic potentials, either classical or first-principles, to Molecular Dynamics simulations of shock waves in rarefied gases, involving vibrational nonequilibrium and chemical reactivity.

[1]  Yoichiro Matsumoto,et al.  Dynamic molecular collision (DMC) model for rarefied gas flow simulations by the DSMC method , 1999 .

[2]  Aleksandar Donev,et al.  Neighbor list collision-driven molecular dynamics simulation for nonspherical hard particles. , 2005 .

[3]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[4]  Giovanni Russo,et al.  Numerical solutions of the Boltzmann equation: comparison of different algorithms , 2008 .

[5]  Gerald Paul.  A Complexity O(1) priority queue for event driven molecular dynamics simulations , 2007, J. Comput. Phys..

[6]  Salomons,et al.  Usefulness of the Burnett description of strong shock waves. , 1992, Physical review letters.

[7]  L. Talbot,et al.  Measurement of Shock Wave Thickness by the Electron Beam Fluorescence Method , 1966 .

[8]  Vladimir Lunev.  Nonequilibrium Gas Flows , 2009 .

[9]  William G. Hoover,et al.  Shock-wave structure via nonequilibrium molecular dynamics and Navier-Stokes continuum mechanics , 1980 .

[10]  E. P. Muntz,et al.  Molecular velocity distribution functions in an argon normal shock wave at Mach number 7 , 1983 .

[11]  Paolo Valentini,et al.  Large-scale molecular dynamics simulations of normal shock waves in dilute argon , 2009 .

[12]  T. G. Cowling,et al.  The mathematical theory of non-uniform gases : notes added in 1951 , 1951 .

[13]  Anna Walsh.  STUDIES IN MOLECULAR DYNAMICS , 1965 .

[14]  Stefan Dietrich,et al.  Scalar and Parallel Optimized Implementation of the Direct Simulation Monte Carlo Method , 1996 .

[15]  Michael N. Macrossan,et al.  v -DSMV: a fast simulation method for rarefield flow , 2001 .

[16]  B. Alder,et al.  Studies in Molecular Dynamics. I. General Method , 1959 .

[17]  G. Pham-Van-diep,et al.  Nonequilibrium Molecular Motion in a Hypersonic Shock Wave , 1989, Science.

[18]  C. F. Curtiss,et al.  Molecular Theory Of Gases And Liquids , 1954 .

[19]  Stefan Schlamp,et al.  Higher moments of the velocity distribution function in dense-gas shocks , 2007, J. Comput. Phys..

[20]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[21]  H. Mott-Smith,et al.  The Solution of the Boltzmann Equation for a Shock Wave , 1951 .

[22]  W. Steckelmacher.  Molecular gas dynamics and the direct simulation of gas flows , 1996 .

[23]  H. Alsmeyer,et al.  Density profiles in argon and nitrogen shock waves measured by the absorption of an electron beam , 1976, Journal of Fluid Mechanics.

[24]  E. P. Muntz,et al.  Molecular Velocity Distribution Function Measurements in a Normal Shock Wave , 1969 .

[25]  Bryan E. Schmidt,et al.  Electron beam density measurements in shock waves in argon , 1969, Journal of Fluid Mechanics.

[26]  S. Miller,et al.  Event-driven molecular dynamics in parallel , 2004 .

[27]  J. M. Haile,et al.  Molecular dynamics simulation : elementary methods / J.M. Haile , 1992 .

[28]  Daniel A. Erwin,et al.  Nonequilibrium gas flows. I: A detailed validation of Monte Carlo direct simulation for monatomic gases , 1991 .

[29]  Aleksandar Donev,et al.  Neighbor list collision-driven molecular dynamics simulation for nonspherical hard particles. I. Algorithmic details , 2005 .

[30]  Carlo Cercignani,et al.  The structure of an infinitely strong shock wave , 1999 .

[31]  B. Alder,et al.  Phase Transition for a Hard Sphere System , 1957 .

[32]  Graeme A. Bird,et al.  Approach to Translational Equilibrium in a Rigid Sphere Gas , 1963 .