A first-order hyperbolic arbitrary Lagrangian Eulerian conservation formulation for non-linear solid dynamics

Di Giusto, T. B.J., Lee, C. H. , Gil, A.J., Bonet, J. and Giacomini, M. (2024) A first-order hyperbolic arbitrary Lagrangian Eulerian conservation formulation for non-linear solid dynamics. International Journal for Numerical Methods in Engineering, (doi: 10.1002/nme.7467) (Early Online Publication)

[img] Text
320635.pdf - Published Version
Available under License Creative Commons Attribution.

9MB

Abstract

The paper introduces a computational framework using a novel Arbitrary Lagrangian Eulerian (ALE) formalism in the form of a system of first-order conservation laws. In addition to the usual material and spatial configurations, an additional referential (intrinsic) configuration is introduced in order to disassociate material particles from mesh positions. Using isothermal hyperelasticity as a starting point, mass, linear momentum and total energy conservation equations are written and solved with respect to the reference configuration. In addition, with the purpose of guaranteeing equal order of convergence of strains/stresses and velocities/displacements, the computation of the standard deformation gradient tensor (measured from material to spatial configuration) is obtained via its multiplicative decomposition into two auxiliary deformation gradient tensors, both computed via additional first-order conservation laws. Crucially, the new ALE conservative formulation will be shown to degenerate elegantly into alternative mixed systems of conservation laws such as Total Lagrangian, Eulerian and Updated Reference Lagrangian. Hyperbolicity of the system of conservation laws will be shown and the accurate wave speed bounds will be presented, the latter critical to ensure stability of explicit time integrators. For spatial discretisation, a vertex-based Finite Volume method is employed and suitably adapted. To guarantee stability from both the continuum and the semi-discretisation standpoints, an appropriate numerical interface flux (by means of the Rankine–Hugoniot jump conditions) is carefully designed and presented. Stability is demonstrated via the use of the time variation of the Hamiltonian of the system, seeking to ensure the positive production of numerical entropy. A range of three dimensional benchmark problems will be presented in order to demonstrate the robustness and reliability of the framework. Examples will be restricted to the case of isothermal reversible elasticity to demonstrate the potential of the new formulation.

Item Type:Articles
Additional Information:Funding information: European Union Horizon 2020, Grant/Award Number: 764636; EPSRC, Grant/Award Number: EP/R008531/1; UK AWE, Grant/Award Number: PO40062030; MCIN, Grant/Award Numbers: PID2020-113463RB-C33, PID2022-141957OB-C21, CEX2018-000797-S.
Status:Early Online Publication
Refereed:Yes
Glasgow Author(s) Enlighten ID:Lee, Dr Chun Hean
Authors: Di Giusto, T. B.J., Lee, C. H., Gil, A.J., Bonet, J., and Giacomini, M.
College/School:College of Science and Engineering > School of Engineering > Infrastructure and Environment
Journal Name:International Journal for Numerical Methods in Engineering
Publisher:Wiley
ISSN:0029-5981
ISSN (Online):1097-0207
Published Online:24 April 2024
Copyright Holders:Copyright: © 2024 The Authors
First Published:First published in International Journal for Numerical Methods in Engineering 2024
Publisher Policy:Reproduced under a Creative Commons licence

University Staff: Request a correction | Enlighten Editors: Update this record

Project CodeAward NoProject NamePrincipal InvestigatorFunder's NameFunder RefLead Dept
300129Strategic Support Package: Engineering of Active Materials by Multiscale/Multiphysics Computational MechanicsChristopher PearceEngineering and Physical Sciences Research Council (EPSRC)EP/R008531/1ENG - Infrastructure & Environment