In this paper, we study the superconvergence behavior of the semi-discrete discontinuous Galerkin (DG) method for {\color{red} scalar} nonlinear hyperbolic equations in one spatial dimension. Superconvergence results for problems with fixed and alternating wind directions are established. On the one hand, we prove that, if the wind direction is fixed (i.e., the derivative of the flux function is bounded away from zero), both the cell average error and numerical flux error at cell interfaces converge at a rate of $2k+1$ when upwind fluxes and piecewise polynomials of degree $k$ are used. Moreover, we also prove that the function value approximation of the DG solution is superconvergent at interior right Radau points, and the derivative value approximation is superconvergent at interior left Radau points, with an order of $k+2$ and $k+1$, respectively. As a byproduct, we show a $k+2$-th order superconvergence of the DG solution towards the Gauss-Radau projection of the exact solution. On the other hand, superconvergence results for problems with alternating wind directions (i.e., the derivative of the flux function either changes sign or otherwise achieves the value zero in the domain) are also established. To be more precise, we first prove that the DG flux function is superconvergent towards a particular flux function of the exact solution, with an order of $k+2$, when Godunov fluxes are used. We then prove that the highest superconvergence rate of the DG solution itself is $k+\frac32$ when sonic points (i.e., the derivative of the flux function achieves zero) appear in the computational domain. As byproducts, we obtain superconvergence properties for the DG solution and the DG flux function at special points and for cell average. Numerical experiments demonstrate that most of our results are optimal, i.e., the superconvergence rates are sharp.

This article focuses on solving the generalized eigenvalue problems (GEP) arising in the source-free Maxwell equation with magnetoelectric coupling effects that models three-dimensional complex media. The goal is to compute the smallest positive eigenvalues, and the main challenge is that the coefficient matrix in the discrete Maxwell equation is indefinite and degenerate. To overcome this difficulty, we derive a singular value decomposition (SVD) of the discrete single-curl operator and then explicitly express the basis of the invariant subspace corresponding to the nonzero eigenvalues of the GEP. Consequently, we reduce the GEP to a null space free standard eigenvalue problem (NFSEP) that contains only the nonzero (complex) eigenvalues of the GEP and can be solved by the shift-and-invert Arnoldi method without being disturbed by the null space. Furthermore, the basis of the eigendecomposition is chosen carefully so that we can apply fast Fourier transformation (FFT)-based matrix vector multiplication to solve the embedded linear systems efficiently by an iterative method. For chiral and pseudochiral complex media, which are of great interest in magnetoelectric applications, the NFSEP can be further transformed to a null space free generalized eigenvalue problem whose coefficient matrices are Hermitian and Hermitian positive definite (HHPD-NFGEP). This HHPD-NFGEP can be solved by using the invert Lanczos method without shifting. Furthermore, the embedded linear system can be solved efficiently by using the conjugate gradient method without preconditioning and the FFT-based matrix vector multiplications. Numerical results are presented to demonstrate the efficiency of the proposed methods.

For decades, the widely used finite difference method on staggered grids, also known as the marker and cell (MAC) method, has been one of the simplest and most effective numerical schemes for solving the Stokes equations and Navier–Stokes equations. Its superconvergence on uniform meshes has been observed by Nicolaides (SIAM J Numer Anal 29(6):1579–1591, 1992), but the rigorous proof is never given. Its behavior on non-uniform grids is not well studied, since most publications only consider uniform grids. In this work, we develop the MAC scheme on non-uniform rectangular meshes, and for the first time we theoretically prove that the superconvergence phenomenon (i.e., second order convergence in the L2 norm for both velocity and pressure) holds true for the MAC method on non-uniform rectangular meshes. With a careful and accurate analysis of various sources of errors, we observe that even though the local truncation errors are only first order in terms of mesh size, the global errors after summation are second order due to the amazing cancellation of local errors. This observation leads to the elegant superconvergence analysis even with non-uniform meshes. Numerical results are given to verify our theoretical analysis.

In this paper, we introduce a nonlocal model for linear steady Stokes system with physical no-slip boundary condition. We use the idea of volume constraint to enforce the no-slip boundary condition. We prove that the nonlocal model is well-posed and the solution of the nonlocal system converges to the solution of the original Stokes system when the nonlocality vanishes.

In this paper, we adapt a simple weighted essentially non-oscillatory (WENO) limiter, originally designed for discontinuous Galerkin (DG) schemes on two-dimensional unstructured triangular meshes, to the correction procedure via reconstruction (CPR) framework for solving nonlinear hyperbolic conservation laws on two-dimensional unstructured triangular meshes with straight or curved edges. This is an extension of our earlier work in which the WENO limiter was designed for the CPR framework on regular meshes. The objective of this simple WENO limiter is to simultaneously maintain uniform high order accuracy of the CPR framework in smooth regions and control spurious numerical oscillations near discontinuities. The WENO limiter we adopt in this paper uses information only from the target cell and its immediate neighbors. Hence, it is particularly simple to implement and will not harm the conservativeness and compactness of the CPR framework. Since the CPR framework with this WENO limiter does not in general satisfy the positivity preserving property, we also extend the positivity-preserving limiters to the CPR framework. Numerical results for both scalar equations and Euler systems of compressible gas dynamics are provided to illustrate the good behavior of this procedure.