High-order adaptive Gegenbauer integral spectral element method for solving non-linear optimal control problems

In this work, we propose an adaptive spectral element algorithm for solving non-linear optimal control problems. The method employs orthogonal collocation at the shifted Gegenbauer–Gauss points combined with very accurate and stable numerical quadratures to fully discretize the multiple-phase integral form of the optimal control problem. The proposed algorithm relies on exploiting the underlying smoothness properties of the solutions for computing approximate solutions efficiently. In particular, the method brackets discontinuities and ‘points of nonsmoothness’ through a novel local adaptive algorithm, which achieves a desired accuracy on the discrete dynamical system equations by adjusting both the mesh size and the degree of the approximating polynomials. A rigorous error analysis of the developed numerical quadratures is presented. Finally, the efficiency of the proposed method is demonstrated on three test examples from the open literature.     

High‐order numerical solution of second‐order one‐dimensional hyperbolic telegraph equation using a shifted Gegenbauer pseudospectral method

We present a high‐order shifted Gegenbauer pseudospectral method (SGPM) to solve numerically the second‐order one‐dimensional hyperbolic telegraph equation provided with some initial and Dirichlet boundary conditions. The framework of the numerical scheme involves the recast of the problem into its integral formulation followed by its discretization into a system of well‐conditioned linear algebraic equations. The integral operators are numerically approximated using some novel shifted Gegenbauer operational matrices of integration. We derive the error formula of the associated numerical quadratures. We also present a method to optimize the constructed operational matrix of integration by minimizing the associated quadrature error in some optimality sense. We study the error bounds and convergence of the optimal shifted Gegenbauer operational matrix of integration. Moreover, we construct the relation between the operational matrices of integration of the shifted Gegenbauer polynomials and standard Gegenbauer polynomials. We derive the global collocation matrix of the SGPM, and construct an efficient computational algorithm for the solution of the collocation equations. We present a study on the computational cost of the developed computational algorithm, and a rigorous convergence and error analysis of the introduced method. Four numerical test examples have been carried out to verify the effectiveness, the accuracy, and the exponential convergence of the method. The SGPM is a robust technique, which can be extended to solve a wide range of problems arising in numerous applications. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 307–349, 2016     

Solving a laminar boundary layer equation with the rational Gegenbauer functions

In this paper, a collocation method using a new weighted orthogonal system on the half-line, namely the rational Gegenbauer functions, is introduced to solve numerically the third-order nonlinear differential equation, af^?+ff^″=0${\mathit{af}}^{?}+{\mathit{ff}}^{″}=0$, where a   is a constant parameter. This method solves the problems on semi-infinite domain without truncating it to a finite domain and transforming the domain of the problems to a finite domain. For a=2$a=2$, the equation is the well-known Blasius equation, which is a laminar viscous flow over a semi-infinite flat plate. We solve this equation by considering 1?a?2$1?a?2$ and compare the new results with the established results to show the efficiency and accuracy of the new method.     

Block boundary value methods for solving Volterra integral and integro-differential equations

Reducible quadrature rules generated by boundary value methods are considered in block version and applied to solve the second kind Volterra integral equations and Volterra integro-differential equations. These extended block boundary value methods are shown to possess both excellent stability properties and high accuracy for Volterra-type equations. Numerical experiments are presented and the efficiency, accuracy and stability of the schemes are confirmed.     

On the numerical solution of linear and nonlinear volterra integral and integro-differential equations

Sinc bases are developed to approximate the solutions of linear and nonlinear Volterra integral and integro-differential equations. Properties of these sinc bases and some operational matrices are first presented. These properties are then used to reduce the integral and integro-differential equations to systems of linear or nonlinear algebraic equations. Numerical examples illustrate the pertinent features of the method and its applicability to a large variety of problems. The examples include convolution type, singular as well as singularly-perturbed problems.     

Existence and the number of solutions of nonresonant semilinear equations and applications to boundary value problems

In this paper, we study the existence and the finiteness of the number of solutions to nonlinear equations and to nonresonant semilinear equations involving nonlinear perturbations of Fredholm maps of index zero. Applications are given to semilinear elliptic boundary value problems and to semilinear parabolic and hyperbolic periodic-boundary value problems.     

Solving boundary value problems of ordinary differential equations with non-separated boundary conditions

This paper focuses on solving the two point boundary value problem, in which boundary conditions are systems of nonlinear equations. The shooting method was used together with a combination of Newton’s method and Broyden’s method, to update the initial values of the differential equations. The experiments showed that the proposed method performed well, in the sense that the overall amount of work was less than that of the Newton Shooting method.     

Boundary value methods for Volterra integral and integro-differential equations

New and effective quadrature rules generated by boundary value methods are introduced. We employ the introduced quadrature rules to construct quadrature methods for the second kind Volterra integral equations and Volterra integro-differential equations. These methods are shown to be effective and possess excellent convergence properties. The nonlinear multigrid method is applied to solve the discrete systems derived from the introduced numerical scheme. Numerical simulations are presented and confirm the efficiency and accuracy of the methods.     

Tension spline collocation methods for singularly perturbed Volterra integro-differential and Volterra integral equations

We consider the numerical discretization of singularly perturbed Volterra integro-differential equations (VIDE)

(*)

and Volterra integral equations (VIE)

(**)

by tension spline collocation methods in certain tension spline spaces, where is a small parameter satisfying 0<1, and q₁, q₂, g and K are functions sufficiently smooth on their domains to ensure that Eqs. (*) and (**) posses a unique solution.We give an analysis of the global convergence properties of a new tension spline collocation solution for 0<1 for singularly perturbed VIDE and VIE; thus, extending the existing theory for =1 to the singularly perturbed case.     

Numerical solution for the weakly singular Fredholm integro-differential equations using Legendre multiwavelets

An effective method based upon Legendre multiwavelets is proposed for the solution of Fredholm weakly singular integro-differential equations. The properties of Legendre multiwavelets are first given and their operational matrices of integral are constructed. These wavelets are utilized to reduce the solution of the given integro-differential equation to the solution of a sparse linear system of algebraic equations. In order to save memory requirement and computational time, a threshold procedure is applied to obtain the solution to this system of algebraic equations. Through numerical examples, performance of the present method is investigated concerning the convergence and the sparseness of the resulted matrix equation.     

Numerical solution of nonlinear Fredholm integral equations of the second kind using Haar wavelets

In this work, we present a computational method for solving nonlinear Fredholm integral equations of the second kind which is based on the use of Haar wavelets. Error analysis is worked out that shows efficiency of the method. Finally, we also give some numerical examples.     

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 rational approximation based on Bernstein polynomials for high order initial and boundary values problems

We introduce a new method to solve high order linear differential equations with initial and boundary conditions numerically. In this method, the approximate solution is based on rational interpolation and collocation method. Since controlling the occurrence of poles in rational interpolation is difficult, a construction which is found by Floater and Hormann [1] is used with no poles in real numbers. We use the Bernstein series solution instead of the interpolation polynomials in their construction. We find that our approximate solution has better convergence rate than the one found by using collocation method. The error of the approximate solution is given in the case of the exact solution f ∈ C^d+2[a, b].     

Eigenvalues of semi-positone Hammerstein integral equations and applications to boundary value problems

Eigenvalue intervals and the existence of finitely many positive eigenfunctions for semi-positone Hammerstein integral equations are obtained. The positive characteristic values and their upper and lower bounds of the corresponding linear Hammerstein integral operators are studied. Applications of the results are given to third-order differential equations with three-point boundary conditions.     

Convergence and stability of extended block boundary value methods for Volterra delay integro-differential equations

In this paper, we construct a class of extended block boundary value methods (B₂VMs) for Volterra delay integro-differential equations and analyze the convergence and stability of the methods. It is proven under the classical Lipschitz condition that an extended B₂VM is convergent of order p if the underlying boundary value methods (BVM) has consistent order p. The analysis shows that a B₂VM extended by an A-stable BVM can preserve the delay-independent stability of the underlying linear systems. Moreover, under some suitable conditions, the extended B₂VMs can also keep the delay-dependent stability of the underlying linear systems. In the end, we test the computational effectiveness by applying the introduced methods to the Volterra delay dynamical model of two interacting species, where the theoretical precision of the methods is further verified.     

Convergence of the collocation method and the mechanical quadrature method for systems of singular integro-differential equations in Lebesgue spaces

Computational schemes for the collocation method and the mechanical quadrature method for the approximate solution of systems of singular integro-differential equations with a Cauchy kernel are elaborated. The case where the systems of equations are defined on an arbitrary smooth closed contour of a complex plane is examined. The methods researched are based on Fejér points. Estimates of the rate of convergence in Lebesgue space are obtained.     

Solving singular second order three-point boundary value problems using reproducing kernel Hilbert space method

This paper investigates the numerical solutions of singular second order three-point boundary value problems using reproducing kernel Hilbert space method. It is a relatively new analytical technique. The solution obtained by using the method takes the form of a convergent series with easily computable components. However, the reproducing kernel Hilbert space method cannot be used directly to solve a singular second order three-point boundary value problem, so we convert it into an equivalent integro-differential equation, which can be solved using reproducing kernel Hilbert space method. Four numerical examples are given to demonstrate the efficiency of the present method. The numerical results demonstrate that the method is quite accurate and efficient for singular second order three-point boundary value problems.     

Wavelet analysis method for solving linear and nonlinear singular boundary value problems

In this paper, a robust and accurate algorithm for solving both linear and nonlinear singular boundary value problems is proposed. We introduce the Chebyshev wavelets operational matrix of derivative and product operation matrix. Chebyshev wavelets expansions together with operational matrix of derivative are employed to solve ordinary differential equations in which, at least, one of the coefficient functions or solution function is not analytic. Several examples are included to illustrate the efficiency and accuracy of the proposed method.