skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Shock waves simulated using the dual domain material point method combined with molecular dynamics

Abstract

Here in this work we combine the dual domain material point method with molecular dynamics in an attempt to create a multiscale numerical method to simulate materials undergoing large deformations with high strain rates. In these types of problems, the material is often in a thermodynamically nonequilibrium state, and conventional constitutive relations or equations of state are often not available. In this method, the closure quantities, such as stress, at each material point are calculated from a molecular dynamics simulation of a group of atoms surrounding the material point. Rather than restricting the multiscale simulation in a small spatial region, such as phase interfaces, or crack tips, this multiscale method can be used to consider nonequilibrium thermodynamic effects in a macroscopic domain. This method takes the advantage that the material points only communicate with mesh nodes, not among themselves; therefore molecular dynamics simulations for material points can be performed independently in parallel. The dual domain material point method is chosen for this multiscale method because it can be used in history dependent problems with large deformation without generating numerical noise as material points move across cells, and also because of its convergence and conservation properties. In conclusion, to demonstrate themore » feasibility and accuracy of this method, we compare the results of a shock wave propagation in a cerium crystal calculated using the direct molecular dynamics simulation with the results from this combined multiscale calculation.« less

Authors:
ORCiD logo [1];  [1]
  1. Los Alamos National Lab. (LANL), Los Alamos, NM (United States). Theoretical Division
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1340934
Alternate Identifier(s):
OSTI ID: 1413342
Report Number(s):
LA-UR-16-24306
Journal ID: ISSN 0021-9991; TRN: US1700968
Grant/Contract Number:  
AC52-06NA25396
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 334; Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
36 MATERIALS SCIENCE

Citation Formats

Zhang, Duan Z., and Dhakal, Tilak Raj. Shock waves simulated using the dual domain material point method combined with molecular dynamics. United States: N. p., 2017. Web. doi:10.1016/j.jcp.2017.01.003.
Zhang, Duan Z., & Dhakal, Tilak Raj. Shock waves simulated using the dual domain material point method combined with molecular dynamics. United States. doi:10.1016/j.jcp.2017.01.003.
Zhang, Duan Z., and Dhakal, Tilak Raj. Tue . "Shock waves simulated using the dual domain material point method combined with molecular dynamics". United States. doi:10.1016/j.jcp.2017.01.003. https://www.osti.gov/servlets/purl/1340934.
@article{osti_1340934,
title = {Shock waves simulated using the dual domain material point method combined with molecular dynamics},
author = {Zhang, Duan Z. and Dhakal, Tilak Raj},
abstractNote = {Here in this work we combine the dual domain material point method with molecular dynamics in an attempt to create a multiscale numerical method to simulate materials undergoing large deformations with high strain rates. In these types of problems, the material is often in a thermodynamically nonequilibrium state, and conventional constitutive relations or equations of state are often not available. In this method, the closure quantities, such as stress, at each material point are calculated from a molecular dynamics simulation of a group of atoms surrounding the material point. Rather than restricting the multiscale simulation in a small spatial region, such as phase interfaces, or crack tips, this multiscale method can be used to consider nonequilibrium thermodynamic effects in a macroscopic domain. This method takes the advantage that the material points only communicate with mesh nodes, not among themselves; therefore molecular dynamics simulations for material points can be performed independently in parallel. The dual domain material point method is chosen for this multiscale method because it can be used in history dependent problems with large deformation without generating numerical noise as material points move across cells, and also because of its convergence and conservation properties. In conclusion, to demonstrate the feasibility and accuracy of this method, we compare the results of a shock wave propagation in a cerium crystal calculated using the direct molecular dynamics simulation with the results from this combined multiscale calculation.},
doi = {10.1016/j.jcp.2017.01.003},
journal = {Journal of Computational Physics},
number = C,
volume = 334,
place = {United States},
year = {Tue Jan 17 00:00:00 EST 2017},
month = {Tue Jan 17 00:00:00 EST 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Save / Share: