Wavelets collocation methods for the numerical solution of elliptic BV problems

Based on collocation with Haar and Legendre wavelets, two efficient and new numerical methods are being proposed for the numerical solution of elliptic partial differential equations having oscillatory and non-oscillatory behavior. The present methods are developed in two stages. In the initial stage, they are developed for Haar wavelets. In order to obtain higher accuracy, Haar wavelets are replaced by Legendre wavelets at the second stage. A comparative analysis of the performance of Haar wavelets collocation method and Legendre wavelets collocation method is carried out. In addition to this, comparative studies of performance of Legendre wavelets collocation method and quadratic spline collocation method, and meshless methods and Sinc–Galerkin method are also done. The analysis indicates that there is a higher accuracy obtained by Legendre wavelets decomposition, which is in the form of a multi-resolution analysis of the function. The solution is first found on the coarse grid points, and then it is refined by obtaining higher accuracy with help of increasing the level of wavelets. The accurate implementation of the classical numerical methods on Neumann’s boundary conditions has been found to involve some difficulty. It has been shown here that the present methods can be easily implemented on Neumann’s boundary conditions and the results obtained are accurate; the present methods, thus, have a clear advantage over the classical numerical methods. A distinct feature of the proposed methods is their simple applicability for a variety of boundary conditions. Numerical order of convergence of the proposed methods is calculated. The results of numerical tests show better accuracy of the proposed method based on Legendre wavelets for a variety of benchmark problems.     

A Haar wavelet collocation approach for solving one and two-dimensional second-order linear and nonlinear hyperbolic telegraph equations

We have developed a new numerical method based on Haar wavelet (HW) in this article for the numerical solution (NS) of one- and two-dimensional hyperbolic Telegraph equations (HTEs). The proposed technique is utilized for one- and two-dimensional linear and nonlinear problems, which shows its advantage over other existing numerical methods. In this technique, we approximated both space and temporal derivatives by the truncated Haar series. The algorithm of the method is simple and we can implement easily in any other programming language. The technique is tested on some linear and nonlinear examples from literature. The maximum absolute errors (MAEs), root mean square errors (RMSEs), and computational convergence rate are calculated for different number of collocation points (CPs) and also some 3D graphs are also drawn. The results show that the proposed technique is simply applicable and accurate.     

A Haar wavelets based approximation for nonlinear inverse problems influenced by unknown heat source

In this discussion, a new numerical algorithm focused on the Haar wavelet is used to solve linear and nonlinear inverse problems with unknown heat source. The heat source is dependent on time and space variables. These types of inverse problems are ill-posed and are challenging to solve accurately. The linearization technique converted the nonlinear problem into simple nonhomogeneous partial differential equation. In this Haar wavelet collocation method (HWCM), the time part is discretized by using finite difference approximation, and space variables are handled by Haar series approximation. The main contribution of the proposed method is transforming this ill-posed problem into well-conditioned algebraic equation with the help of Haar functions, and hence, there is no need to implement any sort of regularization technique. The results of numerical method are efficient and stable for this ill-posed problems containing different noisy levels. We have utilized the proposed method on several numerical examples and have valuable efficiency and accuracy.     

A numerical approach to solve the model of an electromechanical system

In this paper, the model of an electromechanical system, which is a system of linear differential equations, is studied. Haar wavelet collocation method (HWCM) is applied for finding the approximate solution of the model. HWCM reduces the system of the model into a matrix‐vector form that contains the unknown Haar coefficients, and these coefficients are easily calculated. To demonstrate the validity and applicability of HWCM, numerical solutions of the system for different parameter values in the system are presented. The obtained results demonstrate the efficiency and accuracy of the method. All of the computations are performed via a program written in Mathematica.     

An Accurate Space–Time Pseudospectral Method for Solving Nonlinear Multi-Dimensional Heat Transfer Problems

In this paper, we consider the numerical solution of the nonlinear one- and two-dimensional heat transfer problems subject to the given initial conditions and linear Robin boundary conditions. We propose a pseudospectral scheme in both time and spatial discretizations for these problems. The discretization processes are constructed through the multi-variate interpolation of the desired solutions in terms of Chebyshev Gauss Lobbato collocation points. Operational matrices of differentiation are constructed via the tensor products for speeding up of the proposed numerical algorithms’ implementation. Some test problems are provided and the numerical simulations are illustrated to show the spectral accuracy in both space and time of the suggested scheme.     

State Analysis and Optimal Control of Time-Varying Discrete Systems via Haar Wavelets

The state analysis and optimal control of time-varying discrete systems via Haar wavelets are the main tasks of this paper. First, we introduce the definition of discrete Haar wavelets. Then, a comparison between Haar wavelets and other orthogonal functions is given. Based upon some useful properties of the Haar wavelets, a special product matrix and a related coefficient matrix are proposed; also, a shift matrix and a summation matrix are derived. These matrices are very effective in solving our problems. The local property of the Haar wavelets is applied to shorten the calculation procedures.     

Computation of eigenvalues and solutions of regular Sturm–Liouville problems using Haar wavelets

The paper presents a novel method for the computation of eigenvalues and solutions of Sturm–Liouville eigenvalue problems (SLEPs) using truncated Haar wavelet series. This is an extension of the technique proposed by Hsiao to solve discretized version of variational problems via Haar wavelets. The proposed method aims to cover a wider class of problems, by applying it to historically important and a very useful class of boundary value problems, thereby enhancing its applicability. To demonstrate the effectiveness and efficiency of the method various celebrated Sturm–Liouville problems are analyzed for their eigenvalues and solutions. Also, eigensystems are investigated for their asymptotic and oscillatory behavior. The proposed scheme, unlike the conventional numerical schemes, such as Rayleigh quotient and Rayleigh–Ritz approximation, gives eigenpairs simultaneously and provides upper and lower estimates of the smallest eigenvalue, and it is found to have quadratic convergence with increase in resolution.     

Pseudospectral Legendre-based optimal computation of nonlinear constrained variational problems

It is well known that, spectrally accurate solution can be maintained if the grids on which a nonlinear physical problem is to be solved must be obtained by spectrally accurate techniques. In this paper, the pseudospectral Legendre method for general nonlinear smooth and nonsmooth constrained problems of the calculus of variations is studied. The technique is based on spectral collocation methods in which the trajectory, x(t), is approximated by the Nth degree interpolating polynomial, using Legendre-Gauss-Lobatto points as the collocation points, and Lagrange polynomials as trial functions. The integral involved in the formulation of the problem is approximated based on Legendre-Gauss-Lobatto integration rule, thereby reducing the problem to a nonlinear programming one to which existing well-developed algorithms may be applied. The method is easy to implement and yields very accurate results. Illustrative examples are included to confirm the convergence of the pseudospectral Legendre method. Moreover, a numerical experiment (on a nonsmooth problem) indicates that by applying a smoothing filter procedure to the pseudospectral Legendre approximation, one can recover the nonsmooth solution within spectral accuracy.     

Haar wavelet method for solving Fisher’s equation

In this paper, we develop an accurate and efficient Haar wavelet solution of Fisher’s equation, a prototypical reaction-diffusion equation. The solutions of Fisher’s equation are characterized by propagating fronts that can be very steep for large values of the reaction rate coefficient. There is an ongoing effort to better adapt Haar wavelet methods to the solution of differential equations with solutions that resemble shock waves or fronts typical of hyperbolic partial differential equations. Moreover the use of Haar wavelets is found to be accurate, simple, fast, flexible, convenient, small computation costs and computationally attractive.     

Numerical solutions of distributed order fractional differential equations in the time domain using the Müntz–Legendre wavelets approach

In this paper, a numerical method is presented to obtain and analyze the behavior of numerical solutions of distributed order fractional differential equations of the general form in the time domain with the Caputo fractional derivative. The suggested method is based on the Müntz–Legendre wavelet approximation. We derive a new operational vector for the Riemann–Liouville fractional integral of the Müntz–Legendre wavelets by using the Laplace transform method. Applying this operational vector and collocation method in our approach, the problem can be reduced to a system of linear and nonlinear algebraic equations. The arising system can be solved by the Newton method. Discussion on the error bound and convergence analysis for the proposed method is presented. Finally, seven test problems are considered to compare our results with other well‐known methods used for solving these problems. The results in the tabulated tables highlighted that the proposed method is an efficient mathematical tool for analyzing distributed order fractional differential equations of the general form.     

Optimal Control of Linear Time-Varying Systems via Haar Wavelets

This paper introduces the application of Haar wavelets to the optimal control synthesis for linear time-varying systems. Based upon some useful properties of Haar wavelets, a special product matrix, a related coefficient matrix, and an operational matrix of backward integration are proposed to solve the adjoint equation of optimization. The results obtained by the proposed Haar approach are almost the same as those obtained by the conventional Riccati method.     

Coupled Galerkin and parameterization methods for optimal control of discretely connected parallel beams

Galerkin and wavelet methods for optimal boundary control of a couple of discretely connected parallel beams are proposed. First, the problem with boundary controls is converted into a problem with distributed controls. The problem is, then, reduced by a Galerkin-based approach into determining the optimal control of a linear time-invariant lumped parameter system, which will be solved by a wavelet-based method using Legendre wavelets. The integration-operational matrix and Kronecker product are utilized to significantly simplify the optimization problem into a system of linear equations. A numerical example is presented to demonstrate the applicability and the efficiency of the proposed method.     

Numerical solution of nonlinear Volterra-Fredholm-Hammerstein integral equations via collocation method based on radial basis functions

A numerical technique based on the spectral method is presented for the solution of nonlinear Volterra-Fredholm-Hammerstein integral equations. This method is a combination of collocation method and radial basis functions (RBFs) with the differentiation process (DRBF), using zeros of the shifted Legendre polynomial as the collocation points. Different applications of RBFs are used for this purpose. The integral involved in the formulation of the problems are approximated based on Legendre-Gauss-Lobatto integration rule. The results of numerical experiments are compared with the analytical solution in illustrative examples to confirm the accuracy and efficiency of the presented scheme.     

Numerical solution of nonlinear two-dimensional integral equations using rationalized Haar functions

Two-dimensional rationalized Haar (RH) functions are applied to the numerical solution of nonlinear second kind two-dimensional integral equations. Using bivariate collocation method and Newton–Cotes nodes, the numerical solution of these equations is reduced to solving a nonlinear system of algebraic equations. Also, some numerical examples are presented to demonstrate the efficiency and accuracy of the proposed method.     

On the numerical treatment of the eigenparameter dependent boundary conditions

In this paper, we consider the numerical treatment of singular eigenvalue problems supplied with eigenparameter dependent boundary conditions using spectral methods. On the one hand, such boundary conditions hinder the construction of test and trial space functions which could incorporate them and thus providing well-conditioned Galerkin discretization matrices. On the other hand, they can generate surprising behavior of the eigenvectors hardly detected by analytic methods. These singular problems are often indirectly approximated by regular ones. We argue that spectral collocation as well as tau method offer remedies for the first two issues and provide direct and efficient treatment to such problems. On a finite domain, we consider the so-called Petterson-König’s rod eigenvalue problem and on the half line, we take into account the Charney’s baroclinic stability problem and the Fourier eigenvalue problem. One boundary condition in these problems depends on the eigenparameter and additionally, this also could depend on some physical parameters. The Chebyshev collocation based on both, square and rectangular differentiation and a Chebyshev tau method are used to discretize the first problem. All these schemes cast the problems into singular algebraic generalized eigenvalue ones which are solved by the QZ and/or Arnoldi algorithms as well as by some target oriented Jacobi-Davidson methods. Thus, the spurious eigenvalues are completely eliminated. The accuracy of square Chebyshev collocation is roughly estimated and its order of approximation with respect to the eigenvalue of interest is determined. For the problems defined on the half line, we make use of the Laguerre-Gauss-Radau collocation. The method proved to be reliable, accurate, and stable with respect to the order of approximation and the scaling parameter.     

Haar wavelet approximation for magnetohydrodynamic flow equations

This study proposes Haar wavelet (HW) approximation method for solving magnetohydrodynamic flow equations in a rectangular duct in presence of transverse external oblique magnetic field. The method is based on approximating the truncated double Haar wavelets series. Numerical solution of velocity and induced magnetic field is obtained for steady-state, fully developed, incompressible flow for a conducting fluid inside the duct. The calculations show that the accuracy of the Haar wavelet solutions is quite good even in the case of a small number of grid points. The HW approximation method may be used in a wide variety of high-order linear partial differential equations. Application of the HW approximation method showed that it is reliable, simple, fast, least computation at costs and flexible.     

Solving fractional integral equations by the Haar wavelet method

Haar wavelets for the solution of fractional integral equations are applied. Fractional Volterra and Fredholm integral equations are considered. The proposed method also is used for analysing fractional harmonic vibrations. The efficiency of the method is demonstrated by three numerical examples.     

A New Collocation Algorithm for Solving Even-Order Boundary Value Problems via a Novel Matrix Method

This paper is dedicated to presenting and analyzing a numerical algorithm for the solution of even-order boundary value problems. The proposed solutions are spectral and they depend on introducing a new matrix of derivatives of certain shifted Legendre polynomial basis, along with the application of the collocation method. The nonzero elements of the introduced matrix are expressed in terms of the well-known harmonic numbers. Numerical examples provide favorable comparisons with other existing methods and ascertain the efficiency and applicability of the proposed algorithm.     

Free vibration of non-uniform Euler–Bernoulli beam under various supporting conditions using Chebyshev wavelet collocation method

This paper proposes operational matrix of rth integration of Chebyshev wavelets. A general procedure of this matrix is given. Operational matrix of rth integration is taken as rth power of operational matrix of first integration in literature. But, this study removes this disadvantage of Chebyshev wavelets method. Free vibration problems of non-uniform Euler–Bernoulli beam under various supporting conditions are investigated by using Chebyshev Wavelet Collocation Method. The proposed method is based on the approximation by the truncated Chebyshev wavelet series. A homogeneous system of linear algebraic equations has been obtained by using the Chebyshev collocation points. The determinant of coefficients matrix is equated to the zero for nontrivial solution of homogeneous system of linear algebraic equations. Hence, we can obtain ith natural frequencies of the beam and the coefficients of the approximate solution of Chebyshev wavelet series that satisfied differential equation and boundary conditions. Mode shapes functions corresponding to the natural frequencies can be obtained by normalizing of approximate solutions. The computed results well fit with the analytical and numerical results as in the literature. These calculations demonstrate that the accuracy of the Chebyshev wavelet collocation method is quite good even for small number of grid points.     

On spline collocation methods for boundary integral equations in the plane

Nowadays boundary elemen; methods belong to the most popular numerical methods for solving elliptic boundary value problems. They consist in the reduction of the problem to equivalent integral equations (or certain generalizations) on the boundary Γ of the given domain and the approximate solution of these boundary equations. For the numerical treatment the boundary surface is decomposed into a finite number of segments and the unknown functions are approximated by corresponding finite elements and usually determined by collocation and Galerkin procedures. One finds the least difficulties in the theoretical foundation of the convergence of Galerkin methods for certain classes of equations, whereas the convergence of collocation methods, which are mostly used in numerical computations, has yet been proved only for special equations and methods. In the present paper we analyse spline collocation methods on uniform meshes with variable collocation points for one-dimensional pseudodifferential equations on a closed curve with convolutional principal parts, which encompass many classes of boundary integral equations in the plane. We give necessary and sufficient conditions for convergence and prove asymptotic error estimates. In particular we generalize some results on nodal and midpoint collocation obtained in [2], [7] and [8]. The paper is organized as follows. In Section 1 we formulate the problems and the results, Section 2 deals with spline interpolation in periodic Sobolev spaces, and in Section 3 we prove the convergence theorems for the considered collocation methods.