AN EXPONENTIAL WAVE INTEGRATOR PSEUDOSPECTRAL METHOD FOR THE SYMMETRIC REGULARIZED-LONG-WAVE EQUATION

An efficient and accurate exponential wave integrator Fourier pseudospectral （EWI-FP） method is proposed and analyzed for solving the symmetric regularized-long-wave （SRLW） equation, which is used for modeling the weakly nonlinear ion acoustic and space-charge waves. The numerical method here is based on a Gautschi-type exponential wave integrator for temporal approximation and the Fourier pseudospectral method for spatial discretization. The scheme is fully explicit and efficient due to the fast Fourier transform. Numerical analysis of the proposed EWI-FP method is carried out and rigorous error estimates are established without CFL-type condition by means of the mathematical induction. The error bound shows that EWI-FP has second order accuracy in time and spectral accuracy in space. Numerical results are reported to confirm the theoretical studies and indicate that the error bound here is optimal.

Output Feedback Stabilization of an Unstable Wave Equation with Observations Subject to Time Delay

This paper focuses on boundary stabilization of a one-dimensional wave equation with an unstable boundary condition,in which observations are subject to arbitrary fixed time delay.The observability inequality indicates that the open-loop system is observable,based on which the observer and predictor are designed:The state of system is estimated with available observation and then predicted without observation.After that equivalently the authors transform the original system to the well-posed and exponentially stable system by backstepping method.The equivalent system together with the design of observer and predictor give the estimated output feedback.It is shown that the closed-loop system is exponentially stable.Numerical simulations are presented to illustrate the effect of the stabilizing controller.

Error estimates of numerical methods for the nonlinear Dirac equation in the nonrelativistic limit regime

We present several numerical methods and establish their error estimates for the discretization of the nonlinear Dirac equation （NLDE） in the nonrelativistic limit regime, involving a small dimensionless parameter 0 〈 ε〈〈1 which is inversely proportional to the speed of light. In this limit regime, the solution is highly oscillatory in time, i.e., there are propagating waves with wavelength O（ ε^2） and O（1） in time and space, respectively. We begin with the conservative Crank-Nicolson finite difference （CNFD） method and establish rigorously its error estimate which depends explicitly on the mesh size h and time step τ- as well as the small parameter 0 〈 ε≤1 Based on the error bound, in order to obtain ＇correct＇ numerical solutions in the nonrelativistic limit regime, i.e., 0 〈 ε≤1 , the CNFD method requests the ε-scalability： τ- = O（ε3） and h = O（√ε）. Then we propose and analyze two numerical methods for the discretization of NLDE by using the Fourier spectral discretization for spatial derivatives combined with the exponential wave integrator and time- splitting technique for temporal derivatives, respectively. Rigorous error bounds for the two numerical methods show that their ε-scalability is improved to τ = O（ε2） and h = O（1） when 0 〈 ε 〈〈 1. Extensive numerical results are reported to confirm our error estimates.