No abstract uploaded!

A Deuflhard-type exponential integrator sine pseudospectral (DEI-SP) method is proposed and analyzed for solving the generalized improved Boussinesq (GIBq) equation. The numerical scheme is based on a second-order exponential integrator for time integration and a sine pseudospectral discretization in space. Rigorous analysis and abundant experiments show that the method converges quadratically and spectrally in time and space, respectively. Finally the DEI-SP method is applied to investigate the complicated and interesting long-time dynamics of the GIBq equation.

In this paper we develop an adaptive edge finite element method based on a reliable and efficient recovery type a posteriori error estimator for time-harmonic Maxwell equations. The asymptotically exact a posteriori error estimator is based on the superconvergence result proved for the lowest-order edge element on triangular grids, where most pairs of triangles sharing a common edge form approximate parallelograms. The efficiency and robustness of the proposed method is demonstrated by extensive numerical experiments for electromagnetic cloaking problems with highly anisotropic permittivity and permeability.

In this paper, we apply semi-Lagrangian discontinuous Galerkin (SLDG) methods for linear hyperbolic equations in one space dimension and analyze the error between the numerical and exact solutions under the L2-norm. In all the previous works, the theoretical analysis of the SLDG method would suggest a suboptimal convergence rate due to the error accumulation over time steps. However, numerical experiments demonstrate an optimal convergence rate and, if the terminal time is large, a superconvergence rate. In this paper, we will prove optimal convergence and optimal superconvergence rates. There are three main difficulties: 1. The error analysis on overlapping meshes. Due to the nature of the semi-Lagrangian time discretization, we need to introduce the background Eulerian mesh and the shifted mesh. The two meshes are staggered, and it is not easy to construct local projections and to handle the error accumulation during time evolution. 2. The superconvergence of time dependent terms under the L2-norm. The error of the numerical and exact solutions can be divided into two parts, the projection error and the time dependent superconvergence term. The projection strongly depends on the superconvergence rates. Therefore, we need to construct a sequence of projections and improve the superconvergence rates gradually. 3. The stopping criterion of the sequence of projections. The sequence of projections are basically of the same form. We need to show that the projections exist up to some certain order since the superconvergence rate cannot be innity. Hence, we will seek some \hidden" condition for the existence of the projections. In this paper, we will solve all the three difficulties and construct several local projections to prove the optimal convergence and superconvergence rates. Numerical experiments verify the theoretical ndings.

In this paper, we propose a model reduction method for solving multiscale elliptic PDEs with random coefficients in the multiquery setting using an optimization approach. The optimization approach enables us to construct a set of localized multiscale data-driven stochastic basis functions that give an optimal approximation property of the solution operator. Our method consists of the offline and online stages. In the offline stage, we construct the localized multiscale data-driven stochastic basis functions by solving an optimization problem. In the online stage, using our basis functions, we can efficiently solve multiscale elliptic PDEs with random coefficients with relatively small computational costs. Therefore, our method is very efficient in solving target problems with many different force functions. The convergence analysis of the proposed method is also presented and has been verified by the numerical simulations.