# 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 »

- Authors:

- 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}

}