The paper proposes a scheme by combining the Runge-Kutta discontinuous Galerkin method with a {delta}-mapping algorithm for solving hyperbolic conservation laws with discontinuous fluxes. This hybrid scheme is particularly applied to nonlinear elasticity in heterogeneous media and multi-class traffic flow with inhomogeneous road conditions. Numerical examples indicate the schemes efficiency in resolving complex waves of the two systems. Moreover, the discussion implies that the so-called {delta}-mapping algorithm can also be combined with any other classical methods for solving similar problems in general.
A moving mesh discontinuous Galerkin method is presented for the numerical solution of hyperbolic conservation laws. The method is a combination of the discontinuous Galerkin method and the mesh movement strategy which is based on the moving mesh partial differential equation approach and moves the mesh continuously in time and orderly in space. It discretizes hyperbolic conservation laws on moving meshes in the quasi-Lagrangian fashion with which the mesh movement is treated continuously and no interpolation is needed for physical variables from the old mesh to the new one. Two convection terms are induced by the mesh movement and their discretization is incorporated naturally in the DG formulation. Numerical results for a selection of one- and two-dimensional scalar and system conservation laws are presented. It is shown that the moving mesh DG method achieves the theoretically predicted order of convergence for problems with smooth solutions and is able to capture shocks and concentrate mesh points in non-smooth regions. Its advantage over uniform meshes and its insensitiveness to mesh smoothness are also demonstrated.
Motivated by considering partial differential equations arising from conservation laws posed on evolving surfaces, a new numerical method for an advection problem is developed and simple numerical tests are performed. The method is based on an unfitted discontinuous Galerkin approach where the surface is not explicitly tracked by the mesh which means the method is extremely flexible with respect to geometry. Furthermore, the discontinuous Galerkin approach is well-suited to capture the advection driven by the evolution of the surface without the need for a space-time formulation, back-tracking trajectories or streamline diffusion. The method is illustrated by a one-dimensional example and numerical results are presented that show good convergence properties for a simple test problem.
In this work we construct reliable a posteriori estimates for some discontinuous Galerkin schemes applied to nonlinear systems of hyperbolic conservation laws. We make use of appropriate reconstructions of the discrete solution together with the relative entropy stability framework. The methodology we use is quite general and allows for a posteriori control of discontinuous Galerkin schemes with standard flux choices which appear in the approximation of conservation laws. In addition to the analysis, we conduct some numerical benchmarking to test the robustness of the resultant estimator.
We present unconditionally energy stable Runge-Kutta (RK) discontinuous Galerkin (DG) schemes for solving a class of fourth order gradient flows. Our algorithm is geared toward arbitrarily high order approximations in both space and time, while energy dissipation remains preserved without imposing any restriction on time steps and meshes. We achieve this in two steps. First, taking advantage of the penalty free DG method introduced by Liu and Yin [J Sci. Comput. 77:467--501, 2018] for spatial discretization, we reformulate an extended linearized ODE system by the energy quadratization (EQ) approach. Second, we apply an s-stage algebraically stable RK method for temporal discretization. The resulting fully discrete DG schemes are linear and unconditionally energy stable. In addition, we introduce a prediction-correction procedure to improve both the accuracy and stability of the scheme. We illustrate the effectiveness of the proposed schemes by numerical tests with benchmark problems.
We present a novel implementation of the modal discontinuous Galerkin (DG) method for hyperbolic conservation laws in two dimensions on graphics processing units (GPUs) using NVIDIAs Compute Unified Device Architecture (CUDA). Both flexible and highly accurate, DG methods accommodate parallel architectures well as their discontinuous nature produces element-local approximations. High performance scientific computing suits GPUs well, as these powerful, massively parallel, cost-effective devices have recently included support for double-precision floating point numbers. Computed examples for Euler equations over unstructured triangle meshes demonstrate the effectiveness of our implementation on an NVIDIA GTX 580 device. Profiling of our method reveals performance comparable to an existing nodal DG-GPU implementation for linear problems.