1 / 25

Error Analysis for Material Point Method and a case study from Gas Dynamics

Error Analysis for Material Point Method and a case study from Gas Dynamics. Le-Thuy Tran and Martin Berzins . Thanks to DOE for funding from 1997-2008. Material Point Method (MPM) Outline History See Brackbill IJNMF 2005 47 693-705. 1963 Harlow PIC methods then CiC ViC methods

janna
Télécharger la présentation

Error Analysis for Material Point Method and a case study from Gas Dynamics

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Error Analysis for Material Point Method and a case study from Gas Dynamics Le-Thuy Tran and Martin Berzins Thanks to DOE for funding from 1997-2008

  2. Material Point Method (MPM) Outline History See Brackbill IJNMF 2005 47 693-705 • 1963 Harlow PIC methods then CiC ViC methods • 1980s -90s Flip methods Brackbill et al. • 1990s Sulsky Brackbill MPM-PIC • 2000+ Sulsky et al. + GIMP Bardenhagen et al Since then proved effective on difficult problems involving large deformations fracture e.g CSAFE [Guilkey et al.] + [Sulsky et al ] + [Brydon] etc etc A KEY ISSUE NOW IS THE THEORETICAL UNDERPINNING OF THE METHOD.

  3. The Material Point Method (MPM) 1 2 3 4 6 5 Errors at each stage of this process Particles with properties (velocity, mass etc) defined on a mesh Particle properties mapped onto mesh points Forces, accelerations, velocities calculated on mesh points Mesh point motion calculated but only the particles moved by mapping velocities back to particles Handles deformation, contact, high strain

  4. Explicit MPM Algorithm Errors • Particle to grid projection Projection Error • Interpolation from grid back to particles Interpolation Error • Time integration of particles Time Error • Finite element + mass matrix lumping FE +Mass Error • Analysis of mesh/particle interplay from [Grigorev,Vshivkov and Fedoruk PIC 2002] - more particles may not be better on a given mesh. Both cell size and particles per cell matter. 1D : 4-8 particles per cell is optimal

  5. MPM Basis Overview grid basis fns particle basis fns S (x) . i Or delta fn .. … point mass point Velocity

  6. Forward Euler time integration n+1 V n.b. semi-implicit Local errors are second order in time Global error is first order in time Second derivatives of x and v are involved in error estimates These terms may be discontinuous

  7. Limits global accuracy To first order Acceleration continuous but its derivative is not Because of linear basis functions and moving points

  8. Particle Acceleration Derivative Potentially Discontinuous as Particle Crosses Boundary Limits global accuracy To first order

  9. Particle Methods applied to Euler Equations Energy Density Pressure Acceleration Positive density requires Note differing particle velocity gradients in different cells may lead to new extrema when particles cross over cells

  10. Arbitrary point Quadrature Error [Hickernell] Where

  11. Mass Projection Error First order in space for small mesh sizes, and may achieve a second order of convergence for large mesh.

  12. Momentum projection error Momentum projection error is also first order in space for small mesh sizes, and second order of convergence for large mesh.

  13. Force Projection Error Using Hickernel: Force Projection Error is first order in space for small mesh sizes, and may achieve up to second order of convergence for large meshes.

  14. Force Projection Error(cont.)

  15. Nodal Acceleration Projection Error Acceleration projection error is one order less in space than the order of convergence for Force projection error and mass projection error, since mass is first order in space.

  16. Nodal Velocity projection (particle to grid) error Let: is first order in space so the order of velocity projection error is determined by

  17. Nodal velocity projection error(cont) The velocity error for different mesh sizes So velocity projection error has constant order as the mesh size get smaller.

  18. Velocity gradient projection error

  19. Velocity gradient projection error(cont) Velocity projection error at nodes for h: Velocity projection error at nodes for 2h: ,We have: With:

  20. Velocity gradient projection error Therefore: So velocity gradient projection error is the same order as velocity projection error.

  21. Final Error Combination Region 1 and Region 3 capture the error at left and right rarefaction Region 4 captures the error at contact discontinuity Region 6 captures the error at shock-front Region 2 is the region of smooth solution of density and velocity Region 5 contains the error propagation from the shock-front

  22. Final Error Combination(cont.)

  23. Burgers’ Problem Solving with MPM: Follow the consequence steps of standard MPM method.

  24. Burgers’ Problem(cont.) Error at T=0.5 for different mesh sizes:

  25. Conclusions Nodal mass, momentum and force projection error converges with order Nodal acceleration projection error is order Nodal velocity projection error is order Nodal velocity gradient projection error is the same order with velocity projection error The overall error at particles as density, energy… will stop converging as the mesh size get smaller.

More Related