A higher-order Lagrangian discontinuous Galerkin hydrodynamic method for solid dynamics

Abstract We present a new multidimensional high-order Lagrangian discontinuous Galerkin (DG) hydrodynamic method that supports hypoelastic and hyperelastic strength models for simulating solid dynamics with higher-order elements. We also present new one-dimensional test problems that have an analytic solution corresponding to a hyperelastic–plastic wave. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated high-order Taylor series polynomials. The stress fields are represented using nodal quantities. The constitutive models used to calculate the deviatoric stress are either a hypoelastic–plastic, infinitesimal strain hyperelastic–plastic, or finite strain hyperelastic–plastic model. These constitutive models require new methods for calculating high-order polynomials for the velocity gradient and deformation gradient in an element. The plasticity associated with the strength model is determined using a radial return method with a J 2 yield criterion and perfect plasticity. The temporal evolution of the governing equations is achieved with the total variation diminishing Runge–Kutta (TVD RK) time integration method. A diverse suite of 1D and 2D test problems are calculated. The new 1D piston test problems, which have analytic solutions for each elastic–plastic model, are presented and calculated to demonstrate the stability and formal accuracy of the various models with the new Lagrangian DG method. 2D test problems are calculated to demonstrate the stability and robustness of the new Lagrangian DG method on multidimensional problems with high-order elements, which have faces that can bend.

[1]  Nathaniel R. Morgan,et al.  An approach for treating contact surfaces in Lagrangian cell-centered hydrodynamics , 2013, J. Comput. Phys..

[2]  Michael Breu,et al.  Digital Libraries in Computer Science: The MeDoc Approach , 1998, Lecture Notes in Computer Science.

[3]  H. S. Udaykumar,et al.  An Eulerian method for computation of multimaterial impact with ENO shock-capturing and sharp interfaces , 2003 .

[4]  Christian Rohde,et al.  An Introduction to Recent Developments in Theory and Numerics for Conservation Laws: Proceedings of the International School on Theory and Numerics for Conservation Laws, Freiburg/Littenweiler, Germany, October 20-24, 1997 , 1999, Theory and Numerics for Conservation Laws.

[5]  Eleuterio F. Toro,et al.  A Second-Order Cell-Centered Lagrangian Method for Two-Dimensional Elastic-Plastic Flows , 2017 .

[6]  Nabil H. Abboud,et al.  Elastoplasticity with linear tetrahedral elements: A variational multiscale method , 2018 .

[7]  Tzanio V. Kolev,et al.  High-Order Curvilinear Finite Element Methods for Lagrangian Hydrodynamics , 2012, SIAM J. Sci. Comput..

[8]  Rémi Abgrall,et al.  A nominally second-order cell-centered Lagrangian scheme for simulating elastic-plastic flows on two-dimensional unstructured grids , 2013, J. Comput. Phys..

[9]  Michael Dumbser,et al.  A direct Arbitrary-Lagrangian-Eulerian ADER-WENO finite volume scheme on unstructured tetrahedral meshes for conservative and non-conservative hyperbolic systems in 3D , 2014, J. Comput. Phys..

[10]  Geoffrey Ingram Taylor,et al.  The use of flat-ended projectiles for determining dynamic yield stress I. Theoretical considerations , 1948, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[11]  G. J. Ball,et al.  A free-Lagrange augmented Godunov method for the simulation of elastic-plastic solids , 2002 .

[12]  M. Wilkins.  Calculation of Elastic-Plastic Flow , 1963 .

[13]  B. Carnes,et al.  A simple, stable, and accurate linear tetrahedral finite element for transient, nearly, and fully incompressible solid dynamics: a dynamic variational multiscale approach , 2016 .

[14]  E. Grüneisen,et al.  Theorie des festen Zustandes einatomiger Elemente , 1912 .

[15]  Veselin Dobrev,et al.  Curvilinear finite elements for Lagrangian hydrodynamics , 2011 .

[16]  Rémi Abgrall,et al.  A discontinuous Galerkin discretization for solving the two-dimensional gas dynamics equations written under total Lagrangian formulation on general unstructured grids , 2014, J. Comput. Phys..

[17]  Nathaniel R. Morgan,et al.  A Lagrangian staggered grid Godunov-like approach for hydrodynamics , 2014, J. Comput. Phys..

[18]  Guglielmo Scovazzi,et al.  Lagrangian shock hydrodynamics on tetrahedral meshes: A stable and accurate variational multiscale approach , 2012, J. Comput. Phys..

[19]  Xiaodong Liu,et al.  Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates , 2018, J. Comput. Phys..

[20]  A. Hunter,et al.  Implementation of a dislocation-density based single-crystal model into a continuum shock hydrodynamics code , 2018 .

[21]  M. Diehl,et al.  A spectral method solution to crystal elasto-viscoplasticity at finite strains , 2013 .

[22]  Michael Ortiz,et al.  A unified approach to finite deformation elastoplastic analysis based on the use of hyperelastic constitutive equations , 1985 .

[23]  Antonio J. Gil,et al.  An upwind vertex centred Finite Volume solver for Lagrangian solid dynamics , 2015, J. Comput. Phys..

[24]  Tzanio V. Kolev,et al.  High order curvilinear finite elements for elastic-plastic Lagrangian dynamics , 2014, J. Comput. Phys..

[25]  Nathaniel R. Morgan,et al.  Compatible, energy conserving, bounds preserving remap of hydrodynamic fields for an extended ALE scheme , 2018, J. Comput. Phys..

[26]  Nathaniel R. Morgan,et al.  A higher-order Lagrangian discontinuous Galerkin hydrodynamic method for elastic-plastic flows , 2019, Comput. Math. Appl..

[27]  C. Bronkhorst,et al.  A model for finite-deformation nonlinear thermomechanical response of single crystal copper under shock conditions , 2013 .

[28]  Rémi Abgrall,et al.  A Cell-Centered Lagrangian Scheme for Two-Dimensional Compressible Flow Problems , 2007, SIAM J. Sci. Comput..

[29]  Nabil H. Abboud,et al.  Implicit finite incompressible elastodynamics with linear finite elements: A stabilized method in rate form , 2016 .

[30]  Shudao Zhang,et al.  A new high-order discontinuous Galerkin spectral finite element method for Lagrangian gas dynamics in two-dimensions , 2011, J. Comput. Phys..

[31]  Mark A. Kenamond,et al.  Coupling continuum dislocation transport with crystal plasticity for application to shock loading conditions , 2016 .

[32]  Bernardo Cockburn,et al.  A hybridizable discontinuous Galerkin formulation for non-linear elasticity , 2015 .

[33]  Nabil H. Abboud,et al.  A dynamic variational multiscale method for viscoelasticity using linear tetrahedral elements , 2017 .

[34]  Dmitri Kuzmin,et al.  Slope limiting for discontinuous Galerkin approximations with a possibly non‐orthogonal Taylor basis , 2013 .

[35]  T. Bieler,et al.  Overview of constitutive laws, kinematics, homogenization and multiscale methods in crystal plasticity finite-element modeling: Theory, experiments, applications , 2010 .

[36]  Nathaniel R. Morgan,et al.  Reduction of dissipation in Lagrange cell-centered hydrodynamics (CCH) through corner gradient reconstruction (CGR) , 2015, J. Comput. Phys..

[37]  J. C. Simo,et al.  Remarks on rate constitutive equations for finite deformation problems: computational implications , 1984 .

[38]  Nathaniel R. Morgan,et al.  A cell-centered Lagrangian Godunov-like method for solid dynamics , 2013 .

[39]  Philip Eisenlohr,et al.  An elasto-viscoplastic formulation based on fast Fourier transforms for the prediction of micromechanical fields in polycrystalline materials , 2012 .

[40]  Thomas J. R. Hughes,et al.  Stabilized shock hydrodynamics: I. A Lagrangian method , 2007 .

[41]  Nathaniel R. Morgan,et al.  A high-order Lagrangian discontinuous Galerkin hydrodynamic method for quadratic cells using a subcell mesh stabilization scheme , 2019, J. Comput. Phys..

[42]  Wing Kam Liu,et al.  Nonlinear Finite Elements for Continua and Structures , 2000 .

[43]  A. Lew,et al.  A locking-free and optimally convergent discontinuous-Galerkin-based extended finite element method for cracked nearly incompressible solids , 2014 .

[44]  A. Hunter,et al.  Numerical implementation of a crystal plasticity model with dislocation transport for high strain rate applications , 2016 .

[45]  Xijun Yu,et al.  The cell-centered discontinuous Galerkin method for Lagrangian compressible Euler equations in two-dimensions , 2014 .

[46]  Nathaniel R. Morgan,et al.  A Lagrangian discontinuous Galerkin hydrodynamic method , 2018 .

[47]  Nathaniel R. Morgan,et al.  A 3D Lagrangian cell-centered hydrodynamic method with higher-order reconstructions for gas and solid dynamics , 2019, Comput. Math. Appl..