Premium
Stability properties of the Discontinuous Galerkin Material Point Method for hyperbolic problems in one and two space dimensions
Author(s) -
Renaud Adrien,
Heuzé Thomas,
Stainier Laurent
Publication year - 2019
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.6239
Subject(s) - mathematics , discretization , discontinuous galerkin method , galerkin method , mathematical analysis , courant–friedrichs–lewy condition , regular grid , scalar (mathematics) , euler equations , finite volume method , upwind scheme , finite element method , grid , geometry , physics , mechanics , thermodynamics
Summary In this paper, stability conditions are derived for the Discontinuous Galerkin Material Point Method (DGMPM) on the scalar linear advection equation for the sake of simplicity and without loss of generality for linear problems. The discrete systems resulting from the application of the DGMPM discretization in one and two space dimensions are first written. For these problems, a second‐order Runge‐Kutta and the forward Euler time discretizations are respectively considered. Moreover, the numerical fluxes are computed at cell faces by means of either the Donor‐Cell Upwind or the Corner Transport Upwind methods for multidimensional problems. Second, the discrete scheme equations are derived assuming that all cells of a background grid contain at least one particle. Although a Cartesian grid is considered in two space dimensions, the results can be extended to regular grids. The von Neumann linear stability analysis then allows the computation of the critical Courant number for a given space discretization. Although the DGMPM is equivalent to the first‐order finite volume method if one particle lies in each element, so that the Courant number can be set to unity, other distributions of particles may restrict the stability region of the scheme. The study of several configurations is then proposed.