Premium
An Eulerian particle level set method for compressible deforming solids with arbitrary EOS
Author(s) -
Mehmandoust B.,
Pishevar A.R.
Publication year - 2009
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.2607
Subject(s) - eulerian path , mathematics , finite volume method , riemann problem , courant–friedrichs–lewy condition , conservation law , boundary value problem , position (finance) , compressibility , boundary (topology) , mathematical analysis , compressible flow , level set method , conservation of mass , euler equations , mechanics , physics , discretization , computer science , riemann hypothesis , finance , lagrangian , segmentation , artificial intelligence , economics , image segmentation
In this study, an Eulerian, finite‐volume method is developed for the numerical simulation of elastic–plastic response of compressible solid materials with arbitrary equation of state (EOS) under impact loading. The governing equations of mass, momentum, and energy along with evolution equations for deviatoric stresses are solved in Eulerian conservation law form. Since the position of material boundaries is determined implicitly by Eulerian schemes, the solution procedure is split into two separate subproblems, which are solved sequentially at each time step. First, the conserved variables are evolved in time with appropriate boundary conditions at the material interfaces. In the present work a fourth‐order central weighted essentially non‐oscillatory shock‐capturing method that was developed for gas dynamics has been extended to high strain rate solids problems. In this method fluxes are determined on a staggered grid at places where solution is smooth. As a result, the method does not rely on the solution of Riemann problems but enjoys the flexibility of using any type of EOS. Boundary conditions at material interfaces are also treated by a special ghost cell approach. Then in the second subproblem, the position of material interfaces is advanced to the new time using a particle level set method. A fifth‐order Godunov‐type central scheme is used to solve the Hamilton–Jacobi equation of level sets in two space dimensions. The capabilities of the proposed method are evaluated at the end by comparing numerical results with the experimental results and the reported benchmark solutions for the Taylor rod impact, spherical groove jetting, and void collapse problems. Copyright © 2009 John Wiley & Sons, Ltd.