Spectral and Pseudospectral Approximations Using Hermite Functions: Application to the Dirac Equation

We consider in this paper spectral and pseudospectral approximations using Hermite functions for PDEs on the whole line. We first develop some basic approximation results associated with the projections and interpolations in the spaces spanned by Hermite functions. These results play important roles in the analysis of the related spectral and pseudospectral methods. We then consider, as an example of applications, spectral and pseudospectral approximations of the Dirac equation using Hermite functions. In particular, these schemes preserve the essential conservation property of the Dirac equation. We also present some numerical results which illustrate the effectiveness of these methods.     

The Linear Rational Pseudospectral Method with Preassigned Poles

We present a linear rational pseudospectral (collocation) method with preassigned poles for solving boundary value problems. It consists in attaching poles to the trial polynomial so as to make it a rational interpolant. Its convergence is proved by transforming the problem into an associated boundary value problem. Numerical examples demonstrate that the rational pseudospectral method is often more efficient than the polynomial method.     

Zhiber-Shabat方程的多辛几何和多辛隐式格式研究

Zhiber-Shabat方程,描述许多重要的物理现象,是一类重要的非线性方程,有着许多广泛的应用前景.本文给出Zhiber-Shabat方程的多辛几何结构和多辛Fourier拟谱方法.数值算例结果表明多辛离散格式具有较好的长时间的数值稳定性.     

FOURIER-CHEBYSHEV PSEUDOSPECTRAL METHOD FOR THREE-DIMENSIONAL VORTICITY EQUATION WITH UNILATERALLY PERIODIC BOUNDARY CONDITION

A Fourier-Chebyshev pseudospectral scheme is proposed for three-dimensionalvorticily equation with unilaterally periodic boundary condition. The generalized stability and convergence are analysed. The numerical results are presented.     

JACOBI PSEUDOSPECTRAL METHOD FOR FOURTH ORDER PROBLEMS

In this paper, we investigate Jacobi pseudospectral method for fourth order problems. We establish some basic results on the Jacobi-Gauss-type interpolations in non-uniformly weighted Sobolev spaces, which serve as important tools in analysis of numerical quadratures, and numerical methods of differential and integral equations. Then we propose Jacobi pseudospectral schemes for several singular problems and multiple-dimensional problems of fourth order. Numerical results demonstrate the spectral accuracy of these schemes, and coincide well with theoretical analysis.     

一类非线性伪抛物型方程的伪谱解法

本文考虑了一类非线性伪抛物型方程的Fourier伪谱方法,建立了该方程的Fourier伪谱方法的半离散格式和全离散格式．并利用Sobolev空间的正交映射理论,给出了这两种格式的误差估计．最后针对全离散格式给出了数值算例,数值结果表明Fourier伪谱格式能正确加解密,且计算误差较小,效率较高,具有较好的稳定性,可用于提高热流密码体制的加解密效率．     

Composite Legendre-Laguerre pseudospectral approximation in unbounded domains

A composite Legendre–Laguerre pseudospectral approximationin unbounded domains is developed. Some approximation resultsare obtained. As an application, a composite pseudospectralscheme is proposed for the Burgers equation on the half-line.The stability and convergence of the scheme are proved. By choosingappropriate base functions, the resulting system of this methodhas a sparse structure and can be solved in parallel. Numericalresults are given to show the efficiency of this new method.     

Chebyshev pseudospectral method for computing numerical solution of convection–diffusion equation

A method for computing highly accurate numerical solutions of 1D convection–diffusion equations is proposed. In this method, the equation is first discretized with respect to the spatial variable, transforming the original problem into a set of ordinary differential equations, and then the resulting system is integrated in time by the fourth-order Runge–Kutta method. Spatial discretization is done by using the Chebyshev pseudospectral collocation method. Before describing the method, we review a finite difference-based method by Salkuyeh [D. Khojasteh Salkuyeh, On the finite difference approximation to the convection–diffusion equation, Appl. Math. Comput. 179 (2006) 79–86], and, contrary to the proposal of the author, we show that this method is not suitable for problems involving time dependent boundary conditions, which calls for revision. Stability analysis based on pseudoeigenvalues to determine the maximum time step for the proposed method is also carried out. Superiority of the proposed method over a revised version of Salkuyeh’s method is verified by numerical examples.     

Numerical methods for computing ground states and dynamics of nonlinear relativistic Hartree equation for boson stars

Efficient and accurate numerical methods are presented for computing ground states and dynamics of the three-dimensional (3D) nonlinear relativistic Hartree equation both without and with an external potential. This equation was derived recently for describing the mean field dynamics of boson stars. In its numerics, due to the appearance of pseudodifferential operator which is defined in phase space via symbol, spectral method is more suitable for the discretization in space than other numerical methods such as finite difference method, etc. For computing ground states, a backward Euler sine pseudospectral (BESP) method is proposed based on a gradient flow with discrete normalization; and respectively, for computing dynamics, a time-splitting sine pseudospectral (TSSP) method is presented based on a splitting technique to decouple the nonlinearity. Both BESP and TSSP are efficient in computation via discrete sine transform, and are of spectral accuracy in spatial discretization. TSSP is of second-order accuracy in temporal discretization and conserves the normalization in discretized level. In addition, when the external potential and initial data for dynamics are spherically symmetric, the original 3D problem collapses to a quasi-1D problem, for which both BESP and TSSP methods are extended successfully with a proper change of variables. Finally, extensive numerical results are reported to demonstrate the spectral accuracy of the methods and to show very interesting and complicated phenomena in the mean field dynamics of boson stars.