A Fast Discontinuous Galerkin Method for a Bond-Based Linear Peridynamic Model Discretized on a Locally Refined Composite Mesh

摘要

We develop a family of fast discontinuous Galerkin (DG) finite element methods for a bond-based linear peridynamic (PD) model in one space dimension. More precisely, we develop a preconditioned fast piecewise-constant DG scheme on a geometrically graded locally refined composite mesh which is suited for the scenario in which the jump discontinuity of the displacement field occurs at the one of the nodes in the original uniform partition. We also develop a preconditioned fast piecewise-linear DG scheme on a uniform mesh that has a second-order convergence rate when the jump discontinuity occurs at one of the computational nodes or has a convergence rate of one-half order otherwise. Motivated by these results, we develop a preconditioned fast hybrid DG scheme that is discretized on a locally uniformly refined composite mesh to numerically simulate the PD model where the jump discontinuity of the displacement field occurs inside a computational cell. We use a piecewise-constant DG scheme on a uniform mesh and a piecewise-linear DG scheme on a locally uniformly refined mesh of mesh size O(h(2)), which has an overall convergence rate of O(h). Because of their nonlocal nature, numerical methods for PD models generate dense stiffness matrices which have O(N-2) memory requirement and O(N-3) computational complexity, where N is the number of computational nodes. In this paper, we explore the structure of the stiffness matrices to develop three preconditioned fast Krylov subspace iterative solvers for the DG method. Consequently, the methods have significantly reduced computational complexity and memory requirement. Numerical results show the utility of the numerical methods.