An acoustic Riemann solver for large strain computational contact dynamics

Runcie, C. J., Lee, C. H. , Haider, J., Gil, A. J. and Bonet, J. (2022) An acoustic Riemann solver for large strain computational contact dynamics. International Journal for Numerical Methods in Engineering, 123(23), pp. 5700-5748. (doi: 10.1002/nme.7085)

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



This paper presents a vertex-centred finite volume algorithm for the explicit dynamic analysis of large strain contact problems. The methodology exploits the use of a system of first order conservation equations written in terms of the linear momentum and a triplet of geometric deformation measures (comprising the deformation gradient tensor, its co-factor and its Jacobian) together with their associated jump conditions. The latter can be used to derive several dynamic contact models ensuring the preservation of hyperbolic characteristic structure across solution discontinuities at the contact interface, a clear advantage over the standard quasi-static contact models where the influence of inertial effects at the contact interface is completely neglected. Taking advantage of the conservative nature of the formalism, both kinetic (traction) and kinematic (velocity) contact interface conditions are explicitly enforced at the fluxes through the use of appropriate jump conditions. Specifically, the kinetic condition is enforced in the usual linear momentum equation, whereas the kinematic condition can now be easily enforced in the geometric conservation equations without requiring a computationally demanding iterative algorithm. Additionally, a Total Variation Diminishing shock capturing technique can be suitably incorporated in order to improve dramatically the performance of the algorithm at the vicinity of shocks. Moreover, and to guarantee stability from the spatial discretisation standpoint, global entropy production is demonstrated through the satisfaction of semi-discrete version of the classical Coleman–Noll procedure expressed in terms of the time rate of the so-called Hamiltonian energy of the system. Finally, a series of numerical examples is examined in order to assess the performance and applicability of the algorithm suitably implemented in OpenFOAM. The knowledge of the potential contact loci between contact interfaces is assumed to be known a priori.

Item Type:Articles
Glasgow Author(s) Enlighten ID:Runcie, Mr Callum and Lee, Dr Chun Hean
Authors: Runcie, C. J., Lee, C. H., Haider, J., Gil, A. J., and Bonet, J.
College/School:College of Science and Engineering > School of Engineering > Infrastructure and Environment
Journal Name:International Journal for Numerical Methods in Engineering
ISSN (Online):1097-0207
Published Online:05 August 2022
Copyright Holders:Copyright © 2022 The Authors
First Published:First published in International Journal for Numerical Methods in Engineering 123(23): 5700-5748
Publisher Policy:Reproduced under a Creative Commons License

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