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.     

Legendre spectral Galerkin method for second-kind Volterra integral equations

The Legendre spectral Galerkin method for the Volterra integral equations of the second kind is proposed in this paper. We provide a rigorous error analysis for the proposed method, which indicates that the numerical errors (in the L ₂ norm) will decay exponentially provided that the kernel function and the source function are sufficiently smooth. Numerical examples are given to illustrate the theoretical results.      

The Legendre Galerkin-Chebyshev collocation method for Burgers-like equations

A Legendre Galerkin–Chebyshev collocation method for Burgers-likeequations is developed. This method is based on the Legendre–Galerkinvariational form, but the nonlinear term and the right-handterm are treated by Chebyshev–Gauss interpolation. Errorestimates of the semi-discrete scheme and the fully discretescheme are given in the L²-norm. Numerical results indicatethat our method is as stable and accurate as the standard Legendrecollocation method, and as efficient and easy to implement asthe standard Chebyshev collocation method.     

Error analysis for solving the Korteweg‐de Vries equation by a Legendre pseudo‐spectral method

A Legendre pseudo‐spectral method is proposed for the Korteweg‐de Vries equation with nonperiodic boundary conditions. Appropriate base functions are chosen to get an efficient algorithm. Error analysis is given for both semi‐discrete and fully discrete schemes. The numerical results confirm to the theoretical analysis. © (2000) John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 513–534, (2000)     

A multistep collocation method for solving advection equations with delay

Advection equations with delay are appeared in the modeling of the dynamics of structured cell populations. In this article, we construct an efficient two-dimensional multistep collocation method for the numerical solution of a class of advection equations with delay. Equations with aftereffect and equations with both aftereffect and retardation of a state variable are considered. Computability of the algorithm and convergence properties of the proposed numerical method are analyzed for solutions in appropriate Sobolev spaces, and it is shown that the proposed scheme enjoys the spectral accuracy. Numerical examples are given and comparison with other existing methods in the literature is made to demonstrate the efficiency, superiority and high accuracy of the presented method.     

Rational Chebyshev collocation method for solving higher‐order linear ordinary differential equations

A collocation method to find an approximate solution of higher‐order linear ordinary differential equation with variable coefficients under the mixed conditions is proposed. This method is based on the rational Chebyshev (RC) Tau method and Taylor‐Chebyshev collocation methods. The solution is obtained in terms of RC functions. Also, illustrative examples are included to demonstrate the validity and applicability of the technique, and performed on the computer using a program written in maple9. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1130–1142, 2011     

带非局部边界条件的热方程的LEGENDRE配置法

对带非局部边界条件的热方程的初边值问题提出了LEGENDRE配置法,并给出其半离散逼近和全离散逼近的稳定性和收敛性分析.数值试验验证了方法的有效性.     

Least‐squares spectral collocation method for the Stokes equations

First‐order system least‐squares spectral collocation methods are presented for the Stokes equations by adopting the first‐order system and modifying the least‐squares functionals in 2 . Then homogeneous Legendre and Chebyshev (continuous and discrete) functionals are shown to be elliptic and continuous with respect to appropriate product weighted norms. The spectral convergence is analyzed for the proposed methods with some numerical experiments. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 128–139, 2004     

Cardinal functions for Legendre pseudo-spectral method for solving the integro-differential equations

In this article, we present a new numerical method to solve the integro-differential equations (IDEs). The proposed method uses the Legendre cardinal functions to express the approximate solution as a finite series. In our method the operational matrix of derivatives is used to reduce IDEs to a system of algebraic equations. To demonstrate the validity and applicability of the proposed method, we present some numerical examples. We compare the obtained numerical results from the proposed method with some other methods. The results show that the proposed algorithm is of high accuracy, more simple and effective.     

A new approach by two‐dimensional wavelets operational matrix method for solving variable‐order fractional partial integro‐differential equations

In this paper, a new computational scheme based on operational matrices (OMs) of two‐dimensional wavelets is proposed for the solution of variable‐order (VO) fractional partial integro‐differential equations (PIDEs). To accomplish this method, first OMs of integration and VO fractional derivative (FD) have been derived using two‐dimensional Legendre wavelets. By implementing two‐dimensional wavelets approximations and the OMs of integration and variable‐order fractional derivative (VO‐FD) along with collocation points, the VO fractional partial PIDEs are reduced into the system of algebraic equations. In addition to this, some useful theorems are discussed to establish the convergence analysis and error estimate of the proposed numerical technique. Furthermore, computational efficiency and applicability are examined through some illustrative examples.     

Determination of a control function in three‐dimensional parabolic equations by Legendre pseudospectral method

A Legendre pseudospectral method is proposed for solving approximately an inverse problem of determining an unknown control parameter p(t) which is the coefficient of the solution u(x, y, z, t) in a diffusion equation in a three‐dimensional region. The diffusion equation is to be solved subject to suitably prescribed initial‐boundary conditions. The presence of the unknown coefficient p(t) requires an extra condition. This extra condition considered as the integral overspecification over the spacial domain. For discretizing the problem, after homogenization of the boundary conditions, we apply the Legendre pseudospectral method in a matrix based manner. As a results a system of nonlinear differential algebraic equations is generated. Then by using suitable transformation, the problem will be converted to a homogeneous time varying system of linear ordinary differential equations. Also a pseudospectral method for efficient solving of the resulted system of ordinary differential equations is proposed. The solution of this system gives the approximation to values of u and p. The matrix based structure of the present method makes it easy to implement. Numerical experiments are presented to demonstrate the accuracy and the efficiency of the proposed computational procedure. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 74‐93, 2012     

Single and multi‐interval Legendre spectral methods in time for parabolic equations

In this article, we take the parabolic equation with Dirichlet boundary conditions as a model to present the Legendre spectral methods both in spatial and in time. Error analysis for the single/multi‐interval schemes in time is given. For the single interval spectral method in time, we obtain a better error estimate in L²‐norm. For the multi‐interval spectral method in time, we obtain the L²‐optimal error estimate in spatial. By choosing approximate trial and test functions, the methods result in algebraic systems with sparse forms. A parallel algorithm is constructed for the multi‐interval scheme in time. Numerical results show the efficiency of the methods. The methods are also applied to parabolic equations with Neumann boundary conditions, Robin boundary conditions and some nonlinear PDEs. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Solution of nonlinear weakly singular Volterra integral equations using the fractional‐order Legendre functions and pseudospectral method

In this article, our main goal is to render an idea to convert a nonlinear weakly singular Volterra integral equation to a non‐singular one by new fractional‐order Legendre functions. The fractional‐order Legendre functions are generated by change of variable on well‐known shifted Legendre polynomials. We consider a general form of singular Volterra integral equation of the second kind. Then the fractional Legendre–Gauss–Lobatto quadratures formula eliminates the singularity of the kernel of the integral equation. Finally, the Legendre pseudospectral method reduces the solution of this problem to the solution of a system of algebraic equations. This method also can be utilized on fractional differential equations as well. The comparison of results of the presented method and other numerical solutions shows the efficiency and accuracy of this method. Also, the obtained maximum error between the results and exact solutions shows that using the present method leads to accurate results and fast convergence for solving nonlinear weakly singular Volterra integral equations. Copyright © 2015 John Wiley & Sons, Ltd.     

A Crank–Nicolson collocation spectral method for the two‐dimensional viscoelastic wave equation

In this paper, we first establish the Crank–Nicolson collocation spectral (CNCS) method for two‐dimensional (2D) viscoelastic wave equation by means of the Chebyshev polynomials. And then, we analyze the existence, uniqueness, stability, and convergence of the CNCS solutions. Finally, we use some numerical experiments to verify the correctness of theoretical analysis. This implies that the CNCS model is very effective for solving the 2D viscoelastic wave equations.     

Numerical solution of two‐dimensional parabolic equation subject to nonstandard boundary specifications using the pseudospectral Legendre method

In this research, the problem of solving the two‐dimensional parabolic equation subject to a given initial condition and nonlocal boundary specifications is considered. A technique based on the pseudospectral Legendre method is proposed for the numerical solution of the studied problem. Several examples are given and the numerical results are shown to demonstrate the efficiently of the newly proposed method. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Legendre multiscaling functions for solving the one‐dimensional parabolic inverse problem

An inverse problem concerning diffusion equation with a source control parameter is investigated. The approximation of the problem is based on the Legendre multiscaling basis. The properties of Legendre multiscaling functions are first presented. These properties together with Galerkin method are then utilized to reduce the inverse problem to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the new technique. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

三阶微分方程的多区域Legendre-Petrov-Galerkin谱方法

提出三阶微分方程初边值问题的多区域Legendre-Petrov-Galerkin谱方法.对于三阶线性微分方程,证明该方法全离散格式的稳定性,并给出L~2-误差估计.进而将该方法和Legendre配置方法相结合,应用于某些非线性问题.数值算例对单区域和多区域方法的结果进行比较.     

A two‐level method in space and time for the Navier‐Stokes equations

A two‐level method in space and time for the time‐dependent Navier‐Stokes equations is considered in this article. The approximate solution u_M∈H_M is decomposed into the large eddy component v∈H_m(m < M) and the small eddy component w∈H. We obtain the large eddy component v by solving a standard Galerkin equation in a coarse‐level subspace H_m with a time step length k, whereas the small eddy component w is derived by solving a linear equation in an orthogonal complement subspace H with a time step length pk, where p is a positive integer. The analysis shows that our two‐level scheme has long‐time stability and can reach the same accuracy as the standard Galerkin method in fine‐level subspace H_M for an appropriate configuration of p and m. Moreover, some numerical examples are provided to complement our theoretical analysis. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Space‐time spectral collocation method for the one‐dimensional sine‐Gordon equation

In this article, we introduce a new space‐time spectral collocation method for solving the one‐dimensional sine‐Gordon equation. We apply a spectral collocation method for discretizing spatial derivatives, and then use the spectral collocation method for the time integration of the resulting nonlinear second‐order system of ordinary differential equations (ODE). Our formulation has high‐order accurate in both space and time. Optimal a priori error bounds are derived in the L²‐norm for the semidiscrete formulation. Numerical experiments show that our formulation have exponential rates of convergence in both space and time. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 670–690, 2015     

A space‐time spectral method for one‐dimensional time fractional convection diffusion equations

In this paper, we propose a space‐time spectral method for solving a class of time fractional convection diffusion equations. Because both fractional derivative and spectral method have global characteristics in bounded domains, we propose a space‐time spectral‐Galerkin method. The convergence result of the method is proved by providing a priori error estimate. Numerical results further confirm the expected convergence rate and illustrate the versatility of our method. Copyright © 2016 John Wiley & Sons, Ltd.