Laguerre Spectral Method for High Order Problems

In this paper, we propose the Laguerre spectral method for high order problems with mixed inhomogeneous boundary conditions. It is also available for approximated solutions growing fast at infinity. The spectral accuracy is proved. Numerical results demonstrate its high effectiveness.     

Legendre-Gauss Spectral Collocation Method for Second Order Nonlinear Delay Differential Equations

In this paper, we present and analyze a single interval Legendre-Gauss spectral collocation method for solving the second order nonlinear delay differential equations with variable delays. We also propose a novel algorithm for the single interval scheme and apply it to the multiple interval scheme for more efficient implementation. Numerical examples are provided to illustrate the high accuracy of the proposed methods.     

A Fourier spectral method for fractional-in-space Cahn–Hilliard equation

In this paper, a fractional extension of the Cahn–Hilliard (CH) phase field model is proposed, i.e. the fractional-in-space CH equation. The fractional order controls the thickness and the lifetime of the interface, which is typically diffusive in integer order case. An unconditionally energy stable Fourier spectral scheme is developed to solve the fractional equation with periodic or Neumann boundary conditions. This method is of spectral accuracy in space and of second-order accuracy in time. The main advantages of this method are that it yields high precision and high efficiency. Moreover, an extra stabilizing term is added to obey the energy decay property while maintaining accuracy and simplicity. Numerical experiments are presented to confirm the accuracy and effectiveness of the proposed method.     

Solving the driven flow by linearized elimination with differential quadrature

In the Navier‐Stokes equations the velocity and the pressure are coupled together by the incompressibility condition div u = 0( u = (u,v)^T) which makes the equations difficult to solve numerically. In this article, a new method named linearized elimination of unknowns with differential quadrature method is applied to the Navier‐Stokes equations. The method is of second‐order accuracy in time and of spectral accuracy in space. The numerical results show that our method is of high accuracy, of good convergence with little computational efforts. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

HIGH ORDER FINITE DIFFERENCE/SPECTRAL METHODS TO A WATER WAVE MODEL WITH NONLOCAL VISCOSITY

In this paper, efficient numerical schemes are proposed for solving the water wave model with nonlocal viscous term that describe the propagation of surface water wave. By using the Caputo fractional derivative definition to approximate the nonlocal fractional operator, finite difference method in time and spectral method in space are constructed for the considered model. The proposed method employs known 5/2 order scheme for fractional derivative and a mixed linearization for the nonlinear term. The analysis shows that the proposed numerical scheme is unconditionally stable and error estimates are provided to predict that the second order backward differentiation plus 5/2 order scheme converges with order 2 in time, and spectral accuracy in space. Several numerical results are provided to verify the efficiency and accuracy of our theoretical claims. Finally, the decay rate of solutions is investigated.     

Petrov-Galerkin Spectral Element Method for Mixed Inhomogeneous Boundary Value Problems on Polygons

The authors investigate Petrov-Galerkin spectral element method. Some results on Legendre irrational quasi-orthogonal approximations are established, which play important roles in Petrov-Galerkin spectral element method for mixed inhomogeneous boundary value problems of partial differential equations defined on polygons. As examples of applications, spectral element methods for two model problems, with the spectral accuracy in certain Jacobi weighted Sobolev spaces, are proposed. The techniques developed in this paper are also applicable to other higher order methods.     

Finite difference spectral approximation for the time‐space fractional telegraph equation and its parameter estimation

The object of this paper is to present the numerical solution of the time‐space fractional telegraph equation. The proposed method is based on the finite difference scheme in temporal direction and Fourier spectral method in spatial direction. The fast Fourier transform (FFT) technique is applied to practical computation. The stability and convergence analysis are strictly proven, which shows that this method is stable and convergent with (2?α) order accuracy in time and spectral accuracy in space. Moreover, the Levenberg‐Marquardt (L‐M) iterative method is employed for the parameter estimation. Finally, some numerical examples are given to confirm the theoretical analysis.     

A Petrov–Galerkin spectral method for fourth‐order problems

In this paper, we consider the Petrov–Galerkin spectral method for fourth‐order elliptic problems on rectangular domains subject to non‐homogeneous Dirichlet boundary conditions. We derive some sharp results on the orthogonal approximations in one and two dimensions, which play important roles in numerical solutions of higher‐order problems. By applying these results to a fourth‐order problem, we establish the H²‐error and L²‐error bounds of the Petrov–Galerkin spectral method. Numerical experiments are provided to illustrate the high accuracy of the proposed method and coincide well with the theoretical analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

Legendre spectral method for the fractional Bratu problem

In this paper, the Legendre spectral collocation method (LSCM) is applied for the solution of the fractional Bratu's equation. It shows the high accuracy and low computational cost of the LSCM compared with some other numerical methods. The fractional Bratu differential equation is transformed into a nonlinear system of algebraic equations for the unknown Legendre coefficients and solved with some spectral collocation methods. Some illustrative examples are also given to show the validity and applicability of this method, and the obtained results are compared with the existing studies to highlight its high efficiency and neglectable error.     

Spectral collocation and radial basis function methods for one-dimensional interface problems

Differential equations with singular sources or discontinuous coefficients yield non-smooth or even discontinuous solutions. This problem is known as the interface problem. High-order numerical solutions suffer from the Gibbs phenomenon in that the accuracy deteriorates if the discontinuity is not properly treated. In this work, we use the spectral and radial basis function methods and present a least squares collocation method to solve the interface problem for one-dimensional elliptic equations. The domain is decomposed into multiple sub-domains; in each sub-domain, the collocation solution is sought. The solution should satisfy more conditions than the given conditions associated with the differential equations, which makes the problem over-determined. To solve the over-determined system, the least squares method is adopted. For the spectral method, the weighted norm method with different scaling factors and the mixed formulation are used. For the radial basis function method, the weighted shape parameter method is presented. Numerical results show that the least squares collocation method provides an accurate solution with high efficacy and that better accuracy is obtained with the spectral method.     

THE NUMERICAL SOLUTION FOR A PARTIAL INTEGRO-DIFFERENTIAL EQUATION WITH A WEAKLY SINGULAR KERNEL

In this paper, a first order semi-discrete method of a partial integro-differential equation with a weakly singular kernel is considered. We apply Galerkin spectral method in one direction, and the inversion technique for the Laplace transform in another direction, the result of the numerical experiment proves the accuracy of this method.     

A pseudo‐spectral method that uses an overlapping multidomain technique for the numerical solution of sine‐Gordon equation in one and two spatial dimensions

In this article, we study an explicit scheme for the solution of sine‐Gordon equation when the space discretization is carried out by an overlapping multidomain pseudo‐spectral technique. By using differentiation matrices, the equation is reduced to a nonlinear system of ordinary differential equations in time that can be discretized with the explicit fourth‐order Runge–Kutta method. To achieve approximation with high accuracy in large domains, the number of space grid points must be large enough. This yields very large and full matrices in the pseudo‐spectral method that causes large memory requirements. The domain decomposition approach provides sparsity in the matrices obtained after the discretization, and this property reduces storage for large matrices and provides economical ways of performing matrix–vector multiplications. Therefore, we propose a multidomain pseudo‐spectral method for the numerical simulation of the sine‐Gordon equation in large domains. Test examples are given to demonstrate the accuracy and capability of the proposed method. Numerical experiments show that the multidomain scheme has an excellent long‐time numerical behavior for the sine‐Gordon equation in one and two dimensions. Copyright © 2013 John Wiley & Sons, Ltd.     

The alternative Legendre Tau method for solving nonlinear multi-order fractional differential equations

In this paper, the alternative Legendre polynomials (ALPs) are used to approximate the solution of a class of nonlinear multi-order fractional differential equations (FDEs). First, the operational matrix of fractional integration of an arbitrary order and the product operational matrix are derived for ALPs. These matrices together with the spectral Tau method are then utilized to reduce the solution of the mentioned equations into the one of solving a system of nonlinear algebraic equations with unknown ALP coefficients of the exact solution. The fractional derivatives are considered in the Caputo sense and the fractional integration is described in the Riemann-Liouville sense. Numerical examples illustrate that the present method is very effective for linear and nonlinear multi-order FDEs and high accuracy solutions can be obtained only using a small number of ALPs.     

Simple finite element method in vorticity formulation for incompressible flows

A very simple and efficient finite element method is introduced for two and three dimensional viscous incompressible flows using the vorticity formulation. This method relies on recasting the traditional finite element method in the spirit of the high order accurate finite difference methods introduced by the authors in another work. Optimal accuracy of arbitrary order can be achieved using standard finite element or spectral elements. The method is convectively stable and is particularly suited for moderate to high Reynolds number flows.

     

Space‐time spectral method for two‐dimensional semilinear parabolic equations

In this paper, a high‐order accurate numerical method for two‐dimensional semilinear parabolic equations is presented. We apply a Galerkin–Legendre spectral method for discretizing spatial derivatives and a spectral collocation method for the time integration of the resulting nonlinear system of ordinary differential equations. Our formulation can be made arbitrarily high‐order accurate in both space and time. Optimal a priori error bound is derived in the L²‐norm for the semidiscrete formulation. Extensive numerical results are presented to demonstrate the convergence property of the method, show our formulation have spectrally accurate in both space and time. John Wiley & Sons, Ltd.     

Numerical approximation of the conservative Allen–Cahn equation by operator splitting method

In this paper, a second‐order fast explicit operator splitting method is proposed to solve the mass‐conserving Allen–Cahn equation with a space–time‐dependent Lagrange multiplier. The space–time‐dependent Lagrange multiplier can preserve the volume of the system and keep small features. Moreover, we analyze the discrete maximum principle and the convergence rate of the fast explicit operator splitting method. The proposed numerical scheme is of spectral accuracy in space and of second‐order accuracy in time, which greatly improves the computational efficiency. Numerical experiments are presented to confirm the accuracy, efficiency, mass conservation, and stability of the proposed method. Copyright © 2017 John Wiley & Sons, Ltd.     

Stability analysis for a fully discrete spectral scheme for Boussinesq systems

In this paper we perform a stability analysis of a fully discrete numerical method for the solution of a family of Boussinesq systems, consisting of a Fourier collocation spectral method for the spatial discretization and a explicit fourth order Runge–Kutta (RK4) scheme for time integration. Our goal is to determine the influence of the parameters, associated to this family of systems, on the efficiency and accuracy of the numerical method. This analysis allows us to identify which regions in the parameter space are most appropriate for obtaining an efficient and accurate numerical solution. We show several numerical examples in order to validate the accuracy, stability and applicability of our MATLAB implementation of the numerical method.     

Spectral approximation for drainage of an Elastico‐viscous liquid and error analysis

A Fourier‐Galerkin spectral method is proposed and used to analyze a system of quasilinear partial differential equations governing the drainage of liquids of the Oldroyd four‐constant type. It is shown that, Fourier‐Galerkin approximations are convergent with spectral accuracy. An efficient and accurate algorithm based on the Fourier‐Galerkin approximations to the system of quasilinear partial differential equations are developed and implemented. Numerical results indicating the high accuracy and effectiveness of this algorithm are presented. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 492–505, 2012     

High‐order energy‐preserving schemes for the improved Boussinesq equation

This article proposes a class of high‐order energy‐preserving schemes for the improved Boussinesq equation. To derive the energy‐preserving schemes, we first discretize the improved Boussinesq equation by Fourier pseudospectral method, which leads to a finite‐dimensional Hamiltonian system. Then, the obtained semidiscrete system is solved by Hamiltonian boundary value methods, which is a newly developed class of energy‐preserving methods. The proposed schemes can reach spectral precision in space, and in time can reach second‐order, fourth‐order, and sixth‐order accuracy, respectively. Moreover, the proposed schemes can conserve the discrete mass and energy to within machine precision. Furthermore, to show the efficiency and accuracy of the proposed methods, the proposed methods are compared with the finite difference methods and the finite volume element method. The results of several numerical experiments are given for the propagation of the single solitary wave, the interaction of two solitary waves and the wave break‐up.     

A new numerical approach for single rational soliton solution of Chen‐Lee‐Liu equation with Riesz fractional derivative in optical fibers

In this paper, time‐splitting spectral approximation technique has been proposed for Chen‐Lee‐Liu (CLL) equation involving Riesz fractional derivative. The proposed numerical technique is efficient, unconditionally stable, and of second‐order accuracy in time and of spectral accuracy in space. Moreover, it conserves the total density in the discretized level. In order to examine the results, with the aid of weighted shifted Grünwald‐Letnikov formula for approximating Riesz fractional derivative, Crank‐Nicolson weighted and shifted Grünwald difference (CN‐WSGD) method has been applied for Riesz fractional CLL equation. The comparison of results reveals that the proposed time‐splitting spectral method is very effective and simple for obtaining single soliton numerical solution of Riesz fractional CLL equation.