Chebyshev‐Legendre spectral method for solving the two‐dimensional vorticity equations with homogeneous Dirichlet conditions

The Chebyshev‐Legendre spectral method for the two‐dimensional vorticity equations is considered. The Legendre Galerkin Chebyshev collocation method is used with the Chebyshev‐Gauss collocation points. The numerical analysis results under the L²‐norm for the Chebyshev‐Legendre method of one‐dimensional case are generalized into that of the two‐dimensional case. The stability and optimal order convergence of the method are proved. Numerical results are given. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Recovering a time‐dependent coefficient in a parabolic equation from overspecified boundary data using the pseudospectral Legendre method

The aim of this article is to discuss the problem of finding the unknown function u(x,t) and the time‐dependent coefficient a(t) in a parabolic partial differential equation. The pseudospectral Legendre method is employed to solve this problem. The results of numerical experiments are given. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

A numerical approach based on ln‐shifted Legendre polynomials for solving a fractional model of pollution

The model of pollution for a system of 3 lakes interconnected by channels is extended using Caputo‐Hadamard fractional derivatives of different orders α_i∈(0,1), i=1,2,3. A numerical approach based on ln‐shifted Legendre polynomials is proposed to solve the considered fractional model. No discretization is needed in our approach. Some numerical experiments are provided to illustrate the presented method.     

Two‐level discretization of the Navier‐Stokes equations with r‐Laplacian subgridscale viscosity

The r‐Laplacian has played an important role in the development of computationally efficient models for applications, such as numerical simulation of turbulent flows. In this article, we examine two‐level finite element approximation schemes applied to the Navier‐Stokes equations with r‐Laplacian subgridscale viscosity, where r is the order of the power‐law artificial viscosity term. In the two‐level algorithm, the solution to the fully nonlinear coarse mesh problem is utilized in a single‐step linear fine mesh problem. When modeling parameters are chosen appropriately, the error in the two‐level algorithm is comparable to the error in solving the fully nonlinear problem on the fine mesh. We provide rigorous numerical analysis of the two‐level approximation scheme and derive scalings which vary based on the coefficient r, coarse mesh size H, fine mesh size h, and filter radius δ. We also investigate the two‐level algorithm in several computational settings, including the 3D numerical simulation of flow past a backward‐facing step at Reynolds number Re = 5100. In all numerical tests, the two‐level algorithm was proven to achieve the same order of accuracy as the standard one‐level algorithm, at a fraction of the computational cost. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

A FEM–DtN formulation for a non‐linear exterior problem in incompressible elasticity

In this paper, we combine the usual finite element method with a Dirichlet‐to‐Neumann (DtN) mapping, derived in terms of an infinite Fourier series, to study the solvability and Galerkin approximations of an exterior transmission problem arising in non‐linear incompressible 2d‐elasticity. We show that the variational formulation can be written in a Stokes‐type mixed form with a linear constraint and a non‐linear main operator. Then, we provide the uniqueness of solution for the continuous and discrete formulations, and derive a Cea‐type estimate for the associated error. In particular, our error analysis considers the practical case in which the DtN mapping is approximated by the corresponding finite Fourier series. Finally, a reliable a posteriori error estimate, well suited for adaptive computations, is also given. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Gauss‐Lobatto‐Legendre‐Birkhoff pseudospectral approximations for the multi‐term time fractional diffusion‐wave equation with Neumann boundaryconditions

A second‐order finite difference/pseudospectral scheme is proposed for numerical approximation of multi‐term time fractional diffusion‐wave equation with Neumann boundary conditions. The scheme is based upon the weighted and shifted Grünwald difference operators approximation of the time fractional calculus and Gauss‐Lobatto‐Legendre‐Birkhoff (GLLB) pseudospectral method for spatial discretization. The unconditionally stability and convergence of the scheme are rigorously proved. Numerical examples are carried out to verify theoretical results.     

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     

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)     

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     

Web‐spline‐based finite element approximation of some quasi‐newtonian flows: Existence‐uniqueness and error bound

This article deals with the web‐spline‐based finite element approximation of quasi‐Newtonian flows. First, we consider the scalar elliptic p‐Laplace problem. Then, we consider quasi‐Newtonian flows where viscosity obeys power law or Carreau law. We prove well‐posedness at the continuous as well as the discrete level. We give some error bounds for the solution of quasi‐Newtonian flow problem based on the web‐spline method. Finally, we provide the numerical results for the p‐Laplace problem. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq31: 54–77, 2015     

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.     

On the Global Relation and the Dirichlet‐to‐Neumann Correspondence

The recently developed Fokas method for solving two‐dimensional Boundary Value Problems (BVP) via the use of global relations is utilized to solve axisymmetric problems in three dimensions. In particular, novel integral representations for the interior and exterior Dirichlet and Neumann problems for the sphere are derived, which recover and improve the already known solutions of these problems. The BVPs considered in this paper can be classically solved using either the finite Legendre transform or the Mellin‐sine transform (which can be derived from the classical Mellin transform in a way similar to the way that the sine transform can be derived from the Fourier transform). The Legendre transform representation is uniformly convergent at the boundary, but it involves a series that is not useful for many applications. The Mellin‐sine transform involves of course an integral but it is not uniformly convergent at the boundary. In this paper: (a) The Legendre transform representation is rederived in a simpler approach using algebraic manipulations instead of solving ODEs. (b) An integral representation, different that the Mellin‐sine transform representation is derived which is uniformly convergent at the boundary. Furthermore, the derivation of the Fokas approach involves only algebraic manipulations, instead of solving an ordinary differential equation.     

Non‐relativistic continuous states in arbitrary dimension for a ring‐shaped pseudo‐Coulomb and energy‐dependent potentials

In this research, the non‐relativistic particle scattering has been investigated for an alternative pseudo‐Coulomb potential plus ring‐shaped and an energy‐dependent potentials in D‐dimensional space. The normalized wave functions of continuous states on the k/2π scale are expressed in terms of the hyper‐geometric series, and formula of phase shifts is presented. Analytical properties of the scattering amplitude and thermodynamics properties are discussed. Some of the numerical results of energy levels have been calculated too. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical solution of multiterm variable‐order fractional differential equations via shifted Legendre polynomials

In this paper, shifted Legendre polynomials will be used for constructing the numerical solution for a class of multiterm variable‐order fractional differential equations. In the proposed method, the shifted Legendre operational matrix of the fractional variable‐order derivatives will be investigated. The fundamental problem is reduced to an algebraic system of equations using the constructed matrix and the collocation technique, which can be solved numerically. The error estimate of the proposed method is investigated. Some numerical examples are presented to prove the applicability, generality, and accuracy of the suggested method.     

Two‐grid methods for time‐harmonic Maxwell equations

In this paper, we develop several two‐grid methods for the Nédélec edge finite element approximation of the time‐harmonic Maxwell equations. We first present a two‐grid method that uses a coarse space to solve the original problem and then use a fine space to solve a corresponding symmetric positive definite problem. Then, we present two types of iterative two‐grid methods, one is to add the kernel of the curl ‐operator in the fine space to a coarse mesh space to solve the original problem and the other is to use an inner iterative method for dealing with the kernel of the curl ‐operator in the fine space and the coarse space, separately. We provide the error estimates for the first two methods and present numerical experiments to show the efficiency of our methods.Copyright © 2012 John Wiley & Sons, Ltd.     

Forbidden minors for wye‐delta‐wye reducibility

A graph is Y Δ Y reducible if it can be reduced to a single vertex by a sequence of series‐parallel reductions and Y Δ Y transformations. The class of Y Δ Y reducible graphs is minor closed. We found a large number of minor minimal Y Δ Y irreducible graphs: a family of 57578 31‐edge graphs and another 40‐edge graph. It is still an open problem to characterize Y Δ Y reducible graphs in terms of a finite set of forbidden minors. © 2004 Wiley Periodicals, Inc. J Graph Theory 47: 317–321, 2004     

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.     

Stability and convergence of two‐grid Crank‐Nicolson extrapolation scheme for the time‐dependent natural convection equations

In this paper, we consider the Crank‐Nicolson extrapolation scheme for the 2D/3D unsteady natural convection problem. Our numerical scheme includes the implicit Crank‐Nicolson scheme for linear terms and the recursive linear method for nonlinear terms. Standard Galerkin finite element method is used to approximate the spatial discretization. Stability and optimal error estimates are provided for the numerical solutions. Furthermore, a fully discrete two‐grid Crank‐Nicolson extrapolation scheme is developed, the corresponding stability and convergence results are derived for the approximate solutions. Comparison from aspects of the theoretical results and computational efficiency, the two‐grid Crank‐Nicolson extrapolation scheme has the same order as the one grid method for velocity and temperature in H¹‐norm and for pressure in L²‐norm. However, the two‐grid scheme involves much less work than one grid method. Finally, some numerical examples are provided to verify the established theoretical results and illustrate the performances of the developed numerical schemes.     

Mixed two‐grid finite difference methods for solving one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations

The aim of this paper is to propose mixed two‐grid finite difference methods to obtain the numerical solution of the one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations. The finite difference equations at all interior grid points form a large‐sparse linear system, which needs to be solved efficiently. The solution cost of this sparse linear system usually dominates the total cost of solving the discretized partial differential equation. The proposed method is based on applying a family of finite difference methods for discretizing the spatial and time derivatives. The obtained system has been solved by two‐grid method, where the two‐grid method is used for solving the large‐sparse linear systems. Also, in the proposed method, the spectral radius with local Fourier analysis is calculated for different values of h and Δt. The numerical examples show the efficiency of this algorithm for solving the one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations. Copyright © 2016 John Wiley & Sons, Ltd.