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 »
- Authors:
-
- 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.
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
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.
@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}
}
Web of Science
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
Tied interface grid material point method for problems with localized extreme deformation
journal, August 2014
- Lian, Y. P.; Zhang, X.; Zhang, F.
- International Journal of Impact Engineering, Vol. 70
Simulation at Extreme-Scale: Co-Design Thinking and Practices
journal, February 2014
- Tian, Rong
- Archives of Computational Methods in Engineering, Vol. 21, Issue 1
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
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
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
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
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
Mass matrix formulation of the FLIP particle-in-cell method
journal, November 1992
- Burgess, D.; Sulsky, D.; Brackbill, J. U.
- Journal of Computational Physics, Vol. 103, Issue 1
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
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
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
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
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
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
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