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.

No abstract uploaded!

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.

In this paper, we consider linear stability issues for one-dimensional hyperbolic conservation laws using a class of conservative high order upwind-biased finite difference schemes, which is a prototype for the weighted essentially non-oscillatory (WENO) schemes, for initial-boundary value problems (IBVP). The inflow boundary is treated by the so-called inverse Lax-Wendroff (ILW) or simplified inverse Lax-Wendroff (SILW) procedure, and the outflow boundary is treated by the classical high order extrapolation. A third order total variation diminishing (TVD) Runge-Kutta time discretization is used in the fully discrete case. Both GKS (Gustafsson, Kreiss and Sundstr\"om) and eigenvalue analysis are performed for both semi-discrete and fully discrete schemes. The two different analysis techniques yield consistent results. Numerical tests are performed to demonstrate the stability results predicted by the analysis.

No abstract uploaded!