Search Menu
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information
Search Scholarly Publications

Title: Material point methods applied to one-dimensional shock waves and dual domain material point method with sub-points

Abstract

Here, using a simple one-dimensional shock problem as an example, the present paper investigates numerical properties of the original material point method (MPM), the generalized interpolation material point (GIMP) method, the convected particle domain interpolation (CPDI) method, and the dual domain material point (DDMP) method. For a weak isothermal shock of ideal gas, the MPM cannot be used with accuracy. With a small number of particles per cell, GIMP and CPDI produce reasonable results. However, as the number of particles increases the methods fail to converge and produce pressure spikes. The DDMP method behaves in an opposite way. With a small number of particles per cell, DDMP results are unsatisfactory. As the number of particles increases, the DDMP results converge to correct solutions, but the large number of particles needed for convergence makes the method very expensive to use in these types of shock wave problems in two- or three-dimensional cases. The cause for producing the unsatisfactory DDMP results is identified. A simple improvement to the method is introduced by using sub-points. With this improvement, the DDMP method produces high quality numerical solutions with a very small number of particles. Although in the present paper, the numerical examples are one-dimensional,more » all derivations are for multidimensional problems. With the technique of approximately tracking particle domains of CPDI, the extension of this sub-point method to multidimensional problems is straightforward. Finally, this new method preserves the conservation properties of the DDMP method, which conserves mass and momentum exactly and conserves energy to the second order in both spatial and temporal discretizations.« less

Authors:
 [1]; ORCiD logo [1]
+ Show Author Affiliations
  1. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1457270
Alternate Identifier(s):
OSTI ID: 1359307
Report Number(s):
LA-UR-15-29121
Journal ID: ISSN 0021-9991; TRN: US1901343
Grant/Contract Number:  
AC52-06NA25396
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 325; Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 97 MATHEMATICS AND COMPUTING; Material point methods; Shock waves

Citation Formats

Dhakal, Tilak Raj, and Zhang, Duan Zhong. Material point methods applied to one-dimensional shock waves and dual domain material point method with sub-points. United States: N. p., 2016. Web. doi:10.1016/j.jcp.2016.08.033.
Copy to clipboard
Dhakal, Tilak Raj, & Zhang, Duan Zhong. Material point methods applied to one-dimensional shock waves and dual domain material point method with sub-points. United States. https://doi.org/10.1016/j.jcp.2016.08.033
Copy to clipboard
Dhakal, Tilak Raj, and Zhang, Duan Zhong. Tue . "Material point methods applied to one-dimensional shock waves and dual domain material point method with sub-points". United States. https://doi.org/10.1016/j.jcp.2016.08.033. https://www.osti.gov/servlets/purl/1457270.
Copy to clipboard
@article{osti_1457270,
title = {Material point methods applied to one-dimensional shock waves and dual domain material point method with sub-points},
author = {Dhakal, Tilak Raj and Zhang, Duan Zhong},
abstractNote = {Here, using a simple one-dimensional shock problem as an example, the present paper investigates numerical properties of the original material point method (MPM), the generalized interpolation material point (GIMP) method, the convected particle domain interpolation (CPDI) method, and the dual domain material point (DDMP) method. For a weak isothermal shock of ideal gas, the MPM cannot be used with accuracy. With a small number of particles per cell, GIMP and CPDI produce reasonable results. However, as the number of particles increases the methods fail to converge and produce pressure spikes. The DDMP method behaves in an opposite way. With a small number of particles per cell, DDMP results are unsatisfactory. As the number of particles increases, the DDMP results converge to correct solutions, but the large number of particles needed for convergence makes the method very expensive to use in these types of shock wave problems in two- or three-dimensional cases. The cause for producing the unsatisfactory DDMP results is identified. A simple improvement to the method is introduced by using sub-points. With this improvement, the DDMP method produces high quality numerical solutions with a very small number of particles. Although in the present paper, the numerical examples are one-dimensional, all derivations are for multidimensional problems. With the technique of approximately tracking particle domains of CPDI, the extension of this sub-point method to multidimensional problems is straightforward. Finally, this new method preserves the conservation properties of the DDMP method, which conserves mass and momentum exactly and conserves energy to the second order in both spatial and temporal discretizations.},
doi = {10.1016/j.jcp.2016.08.033},
journal = {Journal of Computational Physics},
number = C,
volume = 325,
place = {United States},
year = {Tue Aug 30 00:00:00 EDT 2016},
month = {Tue Aug 30 00:00:00 EDT 2016}
}
Copy to clipboard
Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record
https://doi.org/10.1016/j.jcp.2016.08.033
Copyright Statement
Other availability
Search WorldCat Search WorldCat to find libraries that may hold this journal
Citation Metrics:
Cited by: 10 works
Citation information provided by
Web of Science
Save / Share:
Export Metadata
Save to My Library
You must Sign In or Create an Account in order to save documents to your library.

Works referenced in this record:

A particle method for history-dependent materials
journal, September 1994

  • Sulsky, D.; Chen, Z.; Schreyer, H. L.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 118, Issue 1-2
  • DOI: 10.1016/0045-7825(94)90112-0

Tied interface grid material point method for problems with localized extreme deformation
journal, August 2014

Simulation at Extreme-Scale: Co-Design Thinking and Practices
journal, February 2014

Using the material-point method to model sea ice dynamics
journal, January 2007

  • Sulsky, Deborah; Schreyer, Howard; Peterson, Kara
  • Journal of Geophysical Research, Vol. 112, Issue C2
  • DOI: 10.1029/2005JC003329

Modelling of membranes in the material point method with applications: MODELLING OF MEMBRANES IN THE MPM WITH APPLICATIONS
journal, December 2014

  • Hamad, Fursan; Stolle, Dieter; Vermeer, Pieter
  • International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 39, Issue 8
  • DOI: 10.1002/nag.2336

A weighted least squares particle-in-cell method for solid mechanics: A WEIGHTED LEAST SQUARES PIC METHOD FOR SOLID MECHANICS
journal, September 2010

  • Wallstedt, P. C.; Guilkey, J. E.
  • International Journal for Numerical Methods in Engineering, Vol. 85, Issue 13
  • DOI: 10.1002/nme.3041

A convected particle domain interpolation technique to extend applicability of the material point method for problems involving massive deformations
journal, March 2011

  • Sadeghirad, A.; Brannon, R. M.; Burghardt, J.
  • International Journal for Numerical Methods in Engineering, Vol. 86, Issue 12
  • DOI: 10.1002/nme.3110

Material point method applied to multiphase flows
journal, March 2008

  • Zhang, Duan Z.; Zou, Qisu; VanderHeyden, W. Brian
  • Journal of Computational Physics, Vol. 227, Issue 6, p. 3159-3173
  • DOI: 10.1016/j.jcp.2007.11.021

Mass matrix formulation of the FLIP particle-in-cell method
journal, November 1992

Material point method enhanced by modified gradient of shape function
journal, July 2011

  • Zhang, Duan Z.; Ma, Xia; Giguere, Paul T.
  • Journal of Computational Physics, Vol. 230, Issue 16
  • DOI: 10.1016/j.jcp.2011.04.032

Distribution coefficient algorithm for small mass nodes in material point method
journal, October 2010

  • Ma, Xia; Giguere, Paul T.; Jayaraman, Balaji
  • Journal of Computational Physics, Vol. 229, Issue 20
  • DOI: 10.1016/j.jcp.2010.06.041

Second-order convected particle domain interpolation (CPDI2) with enrichment for weak discontinuities at material interfaces: ENRICHED SECOND-ORDER CONVECTED PARTICLE DOMAIN INTERPOLATION (CPDI2)
journal, July 2013

  • Sadeghirad, A.; Brannon, R. M.; Guilkey, J. E.
  • International Journal for Numerical Methods in Engineering, Vol. 95, Issue 11
  • DOI: 10.1002/nme.4526
    Sort by:
    [ × clear filter / sort ]

    Works referencing / citing this record:

    Conservative Taylor least squares reconstruction with application to material point methods
    journal, October 2018

    • Wobbes, Elizaveta; Möller, Matthias; Galavi, Vahid
    • International Journal for Numerical Methods in Engineering, Vol. 117, Issue 3
    • DOI: 10.1002/nme.5956

    Enhancement of the material point method using B-spline basis functions
    journal, July 2017

    • Gan, Yong; Sun, Zheng; Chen, Zhen
    • International Journal for Numerical Methods in Engineering, Vol. 113, Issue 3
    • DOI: 10.1002/nme.5620

    Investigation of the Propagation of Stress Wave in Nickel-Titanium Shape Memory Alloys
    journal, July 2018

    • Cui, Yehui; Zeng, Xiangguo; Chen, Huayan
    • Materials, Vol. 11, Issue 7
    • DOI: 10.3390/ma11071215

    Investigation of the Propagation of Stress Wave in Nickel-Titanium Shape Memory Alloys
    journal, July 2018

    • Cui, Yehui; Zeng, Xiangguo; Chen, Huayan
    • Materials, Vol. 11, Issue 7
    • DOI: 10.3390/ma11071215
      Sort by:
      [ × clear filter / sort ]

      Similar Records in DOE PAGES and OSTI.GOV collections: