On the algorithmic and implementational aspects of a Discontinuous Galerkin method at finite strains

In this work, algorithmic modifications are proposed and analyzed for a recently developed stabilized finite strain Discontinuous Galerkin (DG) method. The distinguishing feature of the original method, referred to as VMDG, is a consistently derived expression for the numerical flux and stability tensor that account for evolving material and geometric nonlinearity in the vicinity of the interface. Herein, the proposed modifications involve simplifications to the residual force vector and tangent stiffness matrix of the VMDG method that lead to formulations similar to other existing DG methods but retain the enhanced definition for the stability parameters. The primary objective is to reduce the costs associated with implementing the method as well as executing simulations while retaining accuracy and flexibility, thereby making the formulation amenable to boarder material classes such as inelasticity. Each simplification has associated implications on the mathematical and algorithmic properties of the method, such as L2L2 convergence rate, solution accuracy, continuity enforcement, and stability of the nonlinear equation solver. These implications are carefully quantified and assessed through a comprehensive numerical performance study. The range of two and three dimensional problems under consideration involves both isotropic and anisotropic materials. Both triangular and quadrilateral element types are employed along with hh and pp refinement. The ability of the proposed methods to produce stable and accurate results for such a broad class of problems is highlighted.