Chun Hean Lee

An upwind vertex centred Finite Volume solver for Lagrangian solid dynamics

A vertex centred Jameson-Schmidt-Turkel (JST) finite volume algorithm was recently introduced by the authors (Aguirre et al., 2014 1) in the context of fast solid isothermal dynamics. The spatial discretisation scheme was constructed upon a Lagrangian two-field mixed (linear momentum and the deformation gradient) formulation presented as a system of conservation laws 2-4. In this paper, the formulation is further enhanced by introducing a novel upwind vertex centred finite volume algorithm with three key novelties. First, a conservation law for the volume map is incorporated into the existing two-field system to extend the range of applications towards the incompressibility limit (Gil et al., 2014 5). Second, the use of a linearised Riemann solver and reconstruction limiters is derived for the stabilisation of the scheme together with an efficient edge-based implementation. Third, the treatment of thermo-mechanical processes through a Mie-Gruneisen equation of state is incorporated in the proposed formulation. For completeness, the study of the eigenvalue structure of the resulting system of conservation laws is carried out to demonstrate hyperbolicity and obtain the correct time step bounds for non-isothermal processes. A series of numerical examples are presented in order to assess the robustness of the proposed methodology. The overall scheme shows excellent behaviour in shock and bending dominated nearly incompressible scenarios without spurious pressure oscillations, yielding second order of convergence for both velocities and stresses.

A two-step Taylor-Galerkin formulation for fast dynamics

Purpose – The purpose of this paper is to present a new stabilised low-order finite element methodology for large strain fast dynamics. Design/methodology/approach – The numerical technique describing the motion is formulated upon the mixed set of first-order hyperbolic conservation laws already presented by Lee et al. (2013) where the main variables are the linear momentum, the deformation gradient tensor and the total energy. The mixed formulation is discretised using the standard explicit two-step Taylor-Galerkin (2TG) approach, which has been successfully employed in computational fluid dynamics (CFD). Unfortunately, the results display non-physical spurious (or hourglassing) modes, leading to the breakdown of the numerical scheme. For this reason, the 2TG methodology is further improved by means of two ingredients, namely a curl-free projection of the deformation gradient tensor and the inclusion of an additional stiffness stabilisation term. Findings – A series of numerical examples are carried out ...