Wavelet Methods for the Solutions of Partial and Fractional Differential Equations Arising in Physical Problems

The subject of fractional calculus has gained considerable popularity and importance during the past three decades or so, mainly due to its demonstrated applications in numerous seemingly diverse and widespread fields of science and engineering. It deals with derivatives and integrals of arbitrary orders. The fractional derivative has been occurring in many physical problems, such as frequency-dependent damping behavior of materials, motion of a large thin plate in a Newtonian fluid, creep and relaxation functions for viscoelastic materials, the PI D controller for the control of dynamical systems etc. Phenomena in electromagnetics, acoustics, viscoelasticity, electrochemistry, control theory, neutron point kinetic model, anomalous diffusion, Brownian motion, signal and image processing, fluid dynamics and material science are well described by differential equations of fractional order. Generally, nonlinear partial differential equations of fractional order are difficult to solve. So for the last few decades, a great deal of attention has been directed towards the solution (both exact and numerical) of these problems. The aim of this dissertation is to present an extensive study of different wavelet methods for obtaining numerical solutions of mathematical problems occurring in disciplines of science and engineering. This present work also provides a comprehensive foundation of different wavelet methods comprising Haar wavelet method, Legendre wavelet method, Legendre multi-wavelet methods, Chebyshev wavelet method, Hermite wavelet method and Petrov-Galerkin method. The intension is to examine the accuracy of various wavelet methods and their efficiency for solving nonlinear fractional differential equations. With the widespread applications of wavelet methods for solving difficult problems in diverse fields of science and engineering such as wave propagation, data compression, image processing, pattern recognition, computer graphics and in medical technology, these methods have been implemented to develop accurate and fast algorithms for solving integral, differential and integro-differential equations, especially those whose solutions are highly localized in position and scale. The main feature of wavelets is its ability to convert the given differential and integral equations to a system of linear or nonlinear algebraic equations, which can be solved by numerical methods. Therefore, our main focus in the present work is to analyze the application of wavelet based transform methods for solving the problem of fractional order partial differential equations. The introductory concept of wavelet, wavelet transform and multi-resolution analysis (MRA) have been discussed in the preliminary chapter. The basic idea of various analytical and numerical methods viz. Variational Iteration Method (VIM), Homotopy Perturbation Method (HPM), Homotopy Analysis Method (HAM), First Integral Method (FIM), Optimal Homotopy Asymptotic Method (OHAM), Haar Wavelet Method, Legendre Wavelet Method, Chebyshev Wavelet Method and Hermite Wavelet Method have been presented in chapter 1. In chapter 2, we have considered both analytical and numerical approach for solving some particular nonlinear partial differential equations like Burgers’ equation, modified Burgers’ equation, Huxley equation, Burgers-Huxley equation and modified KdV equation, which have a wide variety of applications in physical models. Variational Iteration Method and Haar wavelet Method are applied to obtain the analytical and numerical approximate solution of Huxley and Burgers-Huxley equations. Comparisons between analytical solution and numerical solution have been cited in tables and also graphically. The Haar wavelet method has also been applied to solve Burgers’, modified Burgers’, and modified KdV equations numerically. The results thus obtained are compared with exact solutions as well as solutions available in open literature. Error of collocation method has been presented in this chapter. Methods like Homotopy Perturbation Method (HPM) and Optimal Homotopy Asymptotic Method (OHAM) are very powerful and efficient techniques for solving nonlinear PDEs. Using these methods, many functional equations such as ordinary, partial differential equations and integral equations have been solved. We have implemented HPM and OHAM in chapter 3, in order to obtain the analytical approximate solutions of system of nonlinear partial differential equation viz. the Boussinesq-Burgers’ equations. Also, the Haar wavelet method has been applied to obtain the numerical solution of BoussinesqBurgers’ equations. Also, the convergence of HPM and OHAM has been discussed in this chapter. The mathematical modeling and simulation of systems and processes, based on the description of their properties in terms of fractional derivatives, naturally leads to differential equations of fractional order and the necessity to solve such equations. The mathematical preliminaries of fractional calculus, definitions and theorems have been presented in chapter 4. Next, in this chapter, the Haar wavelet method has been analyzed for solving fractional differential equations. The time-fractional Burgers-Fisher, generalized Fisher type equations, nonlinear time- and space-fractional Fokker-Planck equations have been solved by using two-dimensional Haar wavelet method. The obtained results are compared with the Optimal Homotopy Asymptotic Method (OHAM), the exact solutions and the results available in open literature. Comparison of obtained results with OHAM, Adomian Decomposition Method (ADM), VIM and Operational Tau Method (OTM) has been demonstrated in order to justify the accuracy and efficiency of the proposed schemes. The convergence of two-dimensional Haar wavelet technique has been provided at the end of this chapter. In chapter 5, the fractional differential equations such as KdV-Burger-Kuramoto (KBK) equation, seventh order KdV (sKdV) equation and Kaup-Kupershmidt (KK) equation have been solved by using two-dimensional Legendre wavelet and Legendre multi-wavelet methods. The main focus of this chapter is the application of two-dimensional Legendre wavelet technique for solving nonlinear fractional differential equations like timefractional KBK equation, time-fractional sKdV equation in order to demonstrate the efficiency and accuracy of the proposed wavelet method. Similarly in chapter 6, twodimensional Chebyshev wavelet method has been implemented to obtain the numerical solutions of the time-fractional Sawada-Kotera equation, fractional order Camassa-Holm equation and Riesz space-fractional sine-Gordon equations. The convergence analysis has been done for these wavelet methods. In chapter 7, the solitary wave solution of fractional modified Fornberg-Whitham equation has been attained by using first integral method and also the approximate solutions obtained by optimal homotopy asymptotic method (OHAM) are compared with the exact solutions acquired by first integral method. Also, the Hermite wavelet method has been implemented to obtain approximate solutions of fractional modified Fornberg-Whitham equation. The Hermite wavelet method is implemented to system of nonlinear fractional differential equations viz. the fractional Jaulent-Miodek equations. Convergence of this wavelet methods has been discussed in this chapter. Chapter 8 emphasizes on the application of Petrov-Galerkin method for solving the fractional differential equations such as the fractional KdV-Burgers’ (KdVB) equation and the fractional Sharma-TassoOlver equation with a view to exhibit the capabilities of this method in handling nonlinear equation. The main objective of this chapter is to establish the efficiency and accuracy of Petrov-Galerkin method in solving fractional differential equtaions numerically by implementing a linear hat function as the trial function and a quintic B-spline function as the test function. Various wavelet methods have been successfully employed to numerous partial and fractional differential equations in order to demonstrate the validity and accuracy of these procedures. Analyzing the numerical results, it can be concluded that the wavelet methods provide worthy numerical solutions for both classical and fractional order partial differential equations. Finally, it is worthwhile to mention that the proposed wavelet methods are promising and powerful methods for solving fractional differential equations in mathematical physics. This work also aimed at, to make this subject popular and acceptable to engineering and science community to appreciate the universe of wonderful mathematics, which is in between classical integer order differentiation and integration, which till now is not much acknowledged, and is hidden from scientists and engineers. Therefore, our goal is to encourage the reader to appreciate the beauty as well as the usefulness of these numerical wavelet based techniques in the study of nonlinear physical system

Numerical Approximate Methods for Solving Linear and Nonlinear Integral Equations

Integral equation has been one of the essential tools for various area of applied mathematics. In this work, we employed different numerical methods for solving both linear and nonlinear Fredholm integral equations. A goal is to categorize the selected methods and assess their accuracy and efficiency. We discuss challenges faced by researchers in this field, and we emphasize the importance of interdisciplinary effort for advancing the study on numerical methods for solving integral equations. Integral equations can be viewed as equations which are results of transformation of points in a given vector spaces of integrable functions by the use of certain specific integral operators to points in the same space. If, in particular, one is concerned with function spaces spanned by polynomials for which the kernel of the corresponding transforming integral operator is separable being comprised of polynomial functions only, then several approximate methods of solution of integral equations can be developed. This work, specially, deals with the development of different wavelet methods for solving integral and intgro-differential equations. Wavelets theory is a relatively new and emerging area in mathematical research. It has been applied in a wide range of engineering disciplines; particularly, wavelets are very successfully used in signal analysis for waveform representations and segmentations, time frequency analysis, and fast algorithms for easy implementation. Wavelets permit the accurate representation of a variety of functions and operators. Moreover, wavelets establish a connection with fast numerical algorithms. Wavelets can be separated into two distinct types, orthogonal and semi-orthogonal. The preliminary concept of integral equations and wavelets are first presented in Chapter 1. Classification of integral equations, construction of wavelets and multi-resolution analysis (MRA) have been briefly discussed and provided in this chapter. In Chapter 2, different wavelet methods are constructed and function approximation by these methods with convergence analysis have been presented. In Chapter 3, linear semi-orthogonal compactly supported B-spline wavelets together with their dual wavelets have been applied to approximate the solutions of Fredholm integral equations (both linear and nonlinear) of the second kind and their systems. Properties of these wavelets are first presented; these properties are then utilized to reduce the computation of integral equations to some algebraic equations. Convergence analysis of B-spline method has been discussed in this chapter. Again, in Chapter 4, system of nonlinear Fredholm integral equations have been solved by using hybrid Legendre Block-Pulse functions and xiii Bernstein collocation method. In Chapter 5, two practical problems arising from chemical phenomenon, have been modeled as Fredholm- Hammerstein integral equations and solved numerically by different numerical techniques. First, COSMO-RS model has been solved by Bernstein collocation method, Haar wavelet method and Sinc collocation method. Second, Hammerstein integral equation arising from chemical reactor theory has been solved by B-spline wavelet method. Comparison of results have been demonstrated through illustrative examples. In Chapter 6, Legendre wavelet method and Bernoulli wavelet method have been developed to solve system of integro-differential equations. Legendre wavelets along with their operational matrices are developed to approximate the solutions of system of nonlinear Volterra integro-differential equations. Also, nonlinear Volterra weakly singular integro-differential equations system has been solved by Bernoulli wavelet method. The properties of these wavelets are used to reduce the system of integral equations to a system of algebraic equations which can be solved numerically by Newton's method. Rigorous convergence analysis has been done for these wavelet methods. Illustrative examples have been included to demonstrate the validity and applicability of the proposed techniques. In Chapter 7, we have solved the second order Lane-Emden type singular differential equation. First, the second order differential equation is transformed into integro-differential equation and then solved by Legendre multi-wavelet method and Chebyshev wavelet method. Convergence of these wavelet methods have been discussed in this chapter. In Chapter 8, we have developed a efficient collocation technique called Legendre spectral collocation method to solve the Fredholm integro-differential-difference equations with variable coefficients and system of two nonlinear integro-differential equations which arise in biological model. The proposed method is based on the Gauss-Legendre points with the basis functions of Lagrange polynomials. The present method reduces this model to a system of nonlinear algebraic equations and again this algebraic system has been solved numerically by Newton's method. The study of fuzzy integral equations and fuzzy differential equations is an emerging area of research for many authors. In Chapter 9, we have proposed some numerical techniques for solving fuzzy integral equations and fuzzy integro-differential equations. Fundamentals of fuzzy calculus have been discussed in this chapter. Nonlinear fuzzy Hammerstein integral equation has been solved by Bernstein polynomials and Legendre wavelets, and then compared with homotopy analysis method. We have solved nonlinear fuzzy Hammerstein Volterra integral equations with constant delay by Bernoulli wavelet method and then compared with B-spline wavelet method. Finally, fuzzy integro-differential equation has been solved by Legendre wavelet method and compared with homotopy analysis method. In fuzzy case, we have applied two-dimensional numerical methods which are discussed in chapter 2. Convergence analysis and error estimate have been also provided for Bernoulli wavelet method. xiv The study of fractional calculus, fractional differential equations and fractional integral equations has a great importance in the field of science and engineering. Most of the physical phenomenon can be best modeled by using fractional calculus. Applications of fractional differential equations and fractional integral equations create a wide area of research for many researchers. This motivates to work on fractional integral equations, which results in the form of Chapter 10. First, the preliminary definitions and theorems of fractional calculus have been presented in this chapter. The nonlinear fractional mixed Volterra-Fredholm integro-differential equations along with mixed boundary conditions have been solved by Legendre wavelet method. A numerical scheme has been developed by using Petrov-Galerkin method where the trial and test functions are Legendre wavelets basis functions. Also, this method has been applied to solve fractional Volterra integro-differential equations. Uniqueness and existence of the problem have been discussed and the error estimate of the proposed method has been presented in this work. Sinc Galerkin method is developed to approximate the solution of fractional Volterra-Fredholm integro-differential equations with weakly singular kernels. The proposed method is based on the Sinc function approximation. Uniqueness and existence of the problem have been discussed and the error analysis of the proposed method have been presented in this chapte

Wavelets and Fast Numerical Algorithms

Wavelet based algorithms in numerical analysis are similar to other transform methods in that vectors and operators are expanded into a basis and the computations take place in this new system of coordinates. However, due to the recursive definition of wavelets, their controllable localization in both space and wave number (time and frequency) domains, and the vanishing moments property, wavelet based algorithms exhibit new and important properties. For example, the multiresolution structure of the wavelet expansions brings about an efficient organization of transformations on a given scale and of interactions between different neighbouring scales. Moreover, wide classes of operators which naively would require a full (dense) matrix for their numerical description, have sparse representations in wavelet bases. For these operators sparse representations lead to fast numerical algorithms, and thus address a critical numerical issue. We note that wavelet based algorithms provide a systematic generalization of the Fast Multipole Method (FMM) and its descendents. These topics will be the subject of the lecture. Starting from the notion of multiresolution analysis, we will consider the so-called non-standard form (which achieves decoupling among the scales) and the associated fast numerical algorithms. Examples of non-standard forms of several basic operators (e.g. derivatives) will be computed explicitly.Comment: 32 pages, uuencoded tar-compressed LaTeX file. Uses epsf.sty (see `macros'

A numerical solution for nonlinear heat transfer of fin problems using the Haar wavelet quasilinearization method

The aim of this paper is to study the new application of Haar wavelet quasilinearization method (HWQM) to solve one-dimensional nonlinear heat transfer of fin problems. Three different types of nonlinear problems are numerically treated and the HWQM solutions are compared with those of the other method. The effects of temperature distribution of a straight fin with temperature-dependent thermal conductivity in the presence of various parameters related to nonlinear boundary value problems are analyzed and discussed. Numerical results of HWQM gives excellent numerical results in terms of competitiveness and accuracy compared to other numerical methods. This method was proven to be stable, convergent and, easily coded

Numerical Solution of Fuzzy Fractional Differential Equation By Haar Wavelet

In this paper, we deal with a wavelet operational method based on Haar wavelet to solve the fuzzy fractional differential equation in the Caputo derivative sense. To this end, we derive the Haar wavelet operational matrix of the fractional order integration. The given approach provides an efficient method to find the solution and its upper bond error. To complete the discussion, the convergence theorem is subsequently expressed in detail. So far, no paper has used the Harr wavelet method using generalized difference and fuzzy derivatives, and this is the first time we have done so. Finally, the presented examples reflect the accuracy and efficiency of the proposed method