This paper presents an explicit vertex centred finite volume method for the solution of fast transient isothermal large strain solid dynamics via a system of first order hyperbolic conservation laws. Building upon previous work developed by the authors, in the context of alternative numerical discretisations, this paper explores the use of a series of enhancements (both from the formulation and numerical standpoints) in order to explore some limiting scenarios, such as the consideration of near and true incompressibility. Both Total and Updated Lagrangian formulations are presented and compared at the discrete level, where very small differences between both descriptions are observed due to the excellent discrete satisfaction of the involutions. In addition, a matrix-free tailor-made artificial compressibility algorithm is discussed and combined with an angular momentum projection algorithm. A wide spectrum of numerical examples is thoroughly examined. The scheme shows excellent (stable, consistent and accurate) behaviour, in comparison with other methodologies, in compressible, nearly incompressible and truly incompressible bending dominated scenarios, yielding equal second order of convergence for velocities, deviatoric and volumetric components of the stress.
An upwind vertex centred finite volume algorithm for nearly and truly incompressible explicit fast solid dynamic applications: Total and Updated Lagrangian formulations
HASSAN, OSAMA IBRAHIM ISMAIL;Auricchio F.
2019-01-01
Abstract
This paper presents an explicit vertex centred finite volume method for the solution of fast transient isothermal large strain solid dynamics via a system of first order hyperbolic conservation laws. Building upon previous work developed by the authors, in the context of alternative numerical discretisations, this paper explores the use of a series of enhancements (both from the formulation and numerical standpoints) in order to explore some limiting scenarios, such as the consideration of near and true incompressibility. Both Total and Updated Lagrangian formulations are presented and compared at the discrete level, where very small differences between both descriptions are observed due to the excellent discrete satisfaction of the involutions. In addition, a matrix-free tailor-made artificial compressibility algorithm is discussed and combined with an angular momentum projection algorithm. A wide spectrum of numerical examples is thoroughly examined. The scheme shows excellent (stable, consistent and accurate) behaviour, in comparison with other methodologies, in compressible, nearly incompressible and truly incompressible bending dominated scenarios, yielding equal second order of convergence for velocities, deviatoric and volumetric components of the stress.I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.