We propose a structure-preserving doubling algorithm for a quadratic eigenvalue problem arising from the stability analysis of time-delay systems. We are particularly interested in the eigenvalues on the unit circle, which are difficult to estimate. The convergence and backward error of the algorithm are analyzed and three numerical examples are presented. Our experience shows that our algorithm is efficient in comparison to the few existing approaches for small to medium size problems.

No abstract uploaded!

We present a uniformly first order accurate numerical method for solving the Klein-Gordon-Zakharov (KGZ) system with two dimensionless parameters 0 < ε ≤ 1 and 0 <γ ≤ 1, which are inversely proportional to the plasma frequency and the acoustic speed, respectively. In the simultaneous high-plasma-frequency and subsonic limit regime, i.e. ε < γ → 0+, the KGZ system collapses to a cubic Schrödinger equation, and the solution propagates waves with O ( ε2)-wavelength in time and meanwhile contains rapid outgoing initial layers with speed O (1/ γ ) in space due to the incompatibility of the initial data. By presenting a multiscale decomposition of the KGZ system, we propose a multiscale time integrator Fourier pseudospectral method which is explicit, efficient and uniformly accurate for solving the KGZ system for all 0 < ε < γ ≤ 1. Numerical results are reported to show the efficiency and accuracy of scheme. Finally, the method is applied to investigate the convergence rates of the KGZ system to its limiting models when ε < γ → 0+.

In this paper, we focus on error estimates to smooth solutions of semi-discrete discontinuous Galerkin (DG) methods with quadrature rules for scalar conservation laws. The main techniques we use are energy estimate and Taylor expansion first introduced by Zhang and Shu. We show that, with ${P}^k$ (piecewise polynomials of degree $k$) finite elements in 1D problems, if the quadrature over elements is exact for polynomials of degree $(2k)$, error estimates of $O(h^{k+1/2})$ are obtained for general monotone fluxes, and optimal estimates of $O(h^{k+1})$ are obtained for upwind fluxes. For multidimensional problems, if in addition quadrature over edges is exact for polynomials of degree $(2k+1)$, error estimates of $O(h^k)$ are obtained for general monotone fluxes, and $O(h^{k+1/2})$ are obtained for monotone and sufficiently smooth numerical fluxes. Numerical results validate our analysis.

For a class of convection-diffusion equations with variable diffusivity, we construct third order accurate discontinuous Galerkin (DG) schemes on both one and two dimensional rectangular meshes. The DG method with an explicit time stepping can well be applied to nonlinear convection–diffusion equations. It is shown that under suitable time step restrictions, the scaling limiter proposed in Liu and Yu (2014) [23]when coupled with the present DG schemes preserves the solution bounds indicated by the initial data, i.e., the maximum principle, while maintaining uniform third order accuracy. These schemes can be extended to rectangular meshes in three dimension. The crucial for all model scenarios is that an effective test set can be identified to verify the desired bounds of numerical solutions. This is achieved mainly by taking advantage of the flexible form of the diffusive flux and the adaptable decomposition of weighted cell averages. Numerical results are presented to validate the numerical methods and demonstrate their effectiveness.