Skip navigation

A Total Lagrangian upwind Smooth Particle Hydrodynamics algorithm for large strain explicit solid dynamics

A Total Lagrangian upwind Smooth Particle Hydrodynamics algorithm for large strain explicit solid dynamics

Lee, Chun Hean, Gil, Antonio J., Ghavamian, Ataollah and Bonet, Javier ORCID: 0000-0002-0430-5181 (2018) A Total Lagrangian upwind Smooth Particle Hydrodynamics algorithm for large strain explicit solid dynamics. Computer Methods in Applied Mechanics and Engineering, 344. pp. 209-250. ISSN 0045-7825 (Print), 1879-2138 (Online) (doi:https://doi.org/10.1016/j.cma.2018.09.033)

[img] PDF (Author Accepted Manuscript)
22018 BONET_Total Lagrangian_Upwind_Smooth_Particle_Hydrodynamics_Algorithm_2018.pdf - Accepted Version
Restricted to Registered users only until 14 October 2019.

Download (24MB) | Request a copy

Abstract

In previous work [1, 2], Lee et al. introduced a new Smooth Particle Hydrodynamics (SPH) computational framework for large strain explicit solid dynamics with special emphasis on the treatment of near incompressibility. A first order system of hyperbolic equations was presented expressed in terms of the linear momentum and the minors of the deformation, namely the deformation gradient, its co-factor and its Jacobian. Taking advantage of this representation, the suppression of numerical deficiencies (e.g. spurious pressure, long term instability and/or consistency issues) was addressed through well-established stabilisation procedures. In Reference [1], the adaptation of the very efficient Jameson-Schmidt-Turkel algorithm was presented. Reference [2] introduced an adapted variationally consistent Streamline Upwind Petrov Galerkin methodology. In this paper, we now introduce a third alternative stabilisation strategy, extremely competitive, and which does not require the selection of any user-defined artificial stabilisation parameter. Specifically, a characteristic-based Riemann solver in conjunction with a linear reconstruction procedure is used, with the aim to guarantee both consistency and conservation of the overall algorithm. We show that the proposed SPH formulation is very similar in nature to that of the upwind vertex centred Finite Volume Method presented in [3]. In order to extend the application range towards the incompressibility limit, an artificial compressibility algorithm is also developed. Finally, an extensive set of challenging numerical examples is analysed. The new SPH algorithm shows excellent behaviour in compressible, nearly incompressible and truly incompressible scenarios, yielding second order of convergence for velocities, deviatoric and volumetric components of the stress.

Item Type: Article
Uncontrolled Keywords: Conservation laws, SPH, Upwind, Riemann Solver, Explicit dynamics, Incompressibility
Subjects: T Technology > TA Engineering (General). Civil engineering (General)
Faculty / Department / Research Group: Vice-Chancellor's Group
Last Modified: 27 Nov 2018 12:21
Selected for GREAT 2016: None
Selected for GREAT 2017: None
Selected for GREAT 2018: None
Selected for GREAT 2019: None
URI: http://gala.gre.ac.uk/id/eprint/22018

Actions (login required)

View Item View Item

Downloads

Downloads per month over past year

View more statistics