Small Collaboration: Advanced Numerical Methods for Nonlinear Hyperbolic Balance Laws and Their Applications (hybrid meeting)

This small collaborative workshop brought together experts from the Sino-German project working in the field of advanced numerical methods for hyperbolic balance laws. These are particularly important for compressible fluid flows and related systems of equations. The investigated numerical methods were finite volume/finite difference, discontinuous Galerkin methods, and kinetic-type schemes. We have discussed challenging open mathematical research problems in this field, such as multidimensional shock waves, interfaces with different phases or efficient and problem suited adaptive algorithms. Consequently, our main objective was to discuss novel high-order accurate schemes that reliably approximate underlying physical models and preserve important physically relevant properties. Theoretical questions concerning the convergence of numerical methods and proper solution concepts were addressed as well

Divided difference estimates and accuracy enhancement of discontinuous Galerkin methods for nonlinear symmetric systems of hyperbolic conservation laws

In this paper, we investigate the accuracy-enhancement for the discontinuous Galerkin (DG) method for solving one-dimensional nonlinear symmetric systems of hyperbolic conservation laws. For nonlinear equations, the divided difference estimate is an important tool that allows for superconvergence of the post-processed solutions in the local L2-norm. Therefore, we first prove that the L2-norm of the α-th order (1≤ α≤ k+1) divided difference of the DG error with upwind fluxes is of order k+(3-α)/2, provided that the flux Jacobian matrix, f'(u), is symmetric positive definite. Furthermore, using the duality argument, we are able to derive superconvergence estimates of order 2k+(3-α)/2 for the negative-order norm, indicating that some particular compact kernels can be used to extract at least (3k/2+1)-th order superconvergence for nonlinear systems of conservation laws. Numerical experiments are shown to demonstrate the theoretical results

A posteriori analysis of discontinuous galerkin schemes for systems of hyperbolic conservation laws

In this work we construct reliable a posteriori estimates for some semi- (spatially) discrete 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, which leads to error control in the case of smooth solutions. 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

A posteriori analysis of fully discrete method of lines DG schemes for systems of conservation laws

We present reliable a posteriori estimators for some fully discrete schemes applied to nonlinear systems of hyperbolic conservation laws in one space dimension with strictly convex entropy. The schemes are based on a method of lines approach combining discontinuous Galerkin spatial discretization with single- or multi-step methods in time. The construction of the estimators requires a reconstruction in time for which we present a very general framework first for odes and then apply the approach to conservation laws. The reconstruction does not depend on the actual method used for evolving the solution in time. Most importantly it covers in addition to implicit methods also the wide range of explicit methods typically used to solve conservation laws. For the spatial discretization, we allow for standard choices of numerical fluxes. We use reconstructions of the discrete solution together with the relative entropy stability framework, which leads to error control in the case of smooth solutions. We study under which conditions on the numerical flux the estimate is of optimal order pre-shock. While the estimator we derive is computable and valid post-shock for fixed meshsize, it will blow up as the meshsize tends to zero. This is due to a breakdown of the relative entropy framework when discontinuities develop. We conclude with some numerical benchmarking to test the robustness of the derived estimator

An Unsteady Entropy Adjoint Approach for Adaptive Solution of the Shallow-Water Equations

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/90693/1/AIAA-2011-3694-887.pd

Two types of spectral volume methods for 1-D linear hyperbolic equations with degenerate variable coefficients

In this paper, we analyze two classes of spectral volume (SV) methods for one-dimensional hyperbolic equations with degenerate variable coefficients. The two classes of SV methods are constructed by letting a piecewise $k$-th order ($k\ge 1$ is an arbitrary integer) polynomial function satisfy the local conservation law in each {\it control volume} obtained by dividing the interval element of the underlying mesh with $k$ Gauss-Legendre points (LSV) or Radaus points (RSV). The $L^2$-norm stability and optimal order convergence properties for both methods are rigorously proved for general non-uniform meshes. The superconvergence behaviors of the two SV schemes have been also investigated: it is proved that under the $L^2$ norm, the SV flux function approximates the exact flux with $(k+2)$-th order and the SV solution approximates the exact solution with $(k+\frac32)$-th order; some superconvergence behaviors at certain special points and for element averages have been also discovered and proved. Our theoretical findings are verified by several numerical experiments