Research Article | Open Access
Wavelet-Galerkin Quasilinearization Method for Nonlinear Boundary Value Problems
A numerical method is proposed by wavelet-Galerkin and quasilinearization approach for nonlinear boundary value problems. Quasilinearization technique is applied to linearize the nonlinear differential equation and then wavelet-Galerkin method is implemented to linearized differential equations. In each iteration of quasilinearization technique, solution is updated by wavelet-Galerkin method. In order to demonstrate the applicability of proposed method, we consider the various nonlinear boundary value problems.
The Galerkin method  is a very well-known method for finding the numerical solutions of differential equations. According to wavelet-Galerkin method, connection coefficients are the inner products of Daubechies scaling functions and their derivatives, because we are taking Daubechies scaling functions as a Galerkin basis. The exact and explicit representations of the differential operators in orthonormal bases of compactly supported wavelets are described by Beylkin  and he also discussed the sparse representations of shift operators in orthonormal bases of compactly supported wavelets. Latto et al.  gave the connection coefficients for zeroth level of resolution; these connection coefficients are essentially based on an unbounded domain. Chen et al.  provided the way of calculating the connection coefficients on a bounded interval and these finite integrals play a vital role in the wavelet-Galerkin approximation of differential equations. Restrepo and Leaf  reviewed the inner products of Daubechies wavelets and their derivatives and used the connection coefficients for approximation of differential operators.
Amaratunga et al.  implemented the wavelet-Galerkin technique for solving the one-dimensional counterpart of Helmholtz’s equation. Mishra and Sabina  used the wavelet-Galerkin method for solving linear, homogeneous boundary value problems with constant coefficients and compared the obtained solution with the exact solution by using a family of Daubechies wavelets and at different levels of resolution. Daubechies scaling functions as the Galerkin bases were used by Jianhua et al.  and the authors implemented the wavelet-Galerkin method for differential equations with a boundary layer. In order to implement the wavelet-Galerkin method, they considered the linear, inhomogeneous second order boundary value problem with constant coefficients. In [9–13] several applications of wavelet-Galerkin method are done and these applications are for linear boundary value problems. Motivated by the work of authors [6–13], we extended the wavelet-Galerkin method for the solution of nonlinear boundary value problems.
The quasilinearization approach was introduced by Kalaba and Bellman [14, 15] as a generalization of the Newton-Raphson method  to solve individual or systems of nonlinear ordinary and partial differential equations. The quasilinearization approach is suitable for a general nonlinear ordinary or partial differential equations of any order.
Jiwari  used a uniform Haar wavelet method with quasilinearization technique for the approximate solution of Burgers’ equation and compared the results with the solutions obtained by the other numerical methods and the exact solution. The same approach was used by Kaur et al.  for the solutions of nonlinear boundary value problems in which they treated the quadratic nonlinearity of unknown function. In , we extend the Haar wavelet quasilinearization technique for fractional nonlinear initial and boundary value problems.
In this paper, we consider the compactly supported Daubechies scaling functions  as a Galerkin basis and propose a numerical method by combining wavelet-Galerkin method with quasilinearization technique for solving nonlinear boundary value problems. The method deals with not only quadratic nonlinearities but also various other forms of nonlinearities. To the best of our knowledge, the wavelet-Galerkin method, with Daubechies scaling functions as Galerkin basis, has not been implemented for numerical solutions of nonlinear differential equations. Illustrative problems show the advantage of the proposed method.
The paper is arranged as follows: in Section 2, we give a brief review of Daubechies’ wavelets, while in Section 3 we describe the two-term connection coefficients. In Section 4, we give the procedure of implementation of the wavelet-Galerkin method. In Section 5, the quasilinearization technique for dealing with ordinary differential equations is discussed, while in Section 6 we apply the wavelet-Galerkin method with quasilinearization technique to some nonlinear boundary value problems. Finally, in Section 7 we conclude our work.
2. Daubechies’ Wavelets
Daubechies [22, 23] constructed a family of compactly supported orthonormal wavelets. A wavelet system consists of a scaling function and a wavelet function . There are two important relations in wavelet theory, which we called two-scale relation. Consider and the equation where . Relations (1) and (2) are known as refinement relations. The coefficients are called the wavelet filter coefficients. Relations (1) and (2) are also called refinement relations and (an even integer) is the number of wavelet filter coefficients in the refinement relations. The supports of the scaling function and wavelet function are and , respectively. Daubechies  constructed wavelet filter coefficients to satisfy the certain conditions. These conditions are linked with certain properties of scaling and wavelet functions . Consider
Relation (3) shows that scaling functions have unit area and relations (4) and (5) indicate the orthonormality of and orthogonality of and , respectively. Relation (6) shows that mth moment of is zero; that is, it has vanishing moments, which implies that we can express the elements of the set as a linear combination of , integer translate of .
Daubechies wavelet has no explicit expressions for the scaling function and the wavelet function at arbitrary .
The simplest Daubechies wavelet  is the Haar wavelet, which has explicit expression for calculating the scaling function . It is also called D2, which means the Daubechies wavelet having two filter coefficients. It is also called db1, which means the Daubechies wavelet having one vanishing moment. db2 (D4) has four wavelet coefficients, that is, , and two vanishing moments. Similarly, DN has -coefficients and -vanishing moments.
3. Two-Term Connection Coefficients
In the present work, we are dealing with the second order nonlinear boundary value problems. We will be concerned with the 2-term connection coefficients. Two-term connection coefficients are defined as
Take times derivative of the Daubechies scaling function (1), by assuming that it is times differentiable, to obtain
Let be a column vector with components which are connection coefficients: . Equation (9) gives a system of linear equations with as unknown vector; we can write (9) in vector form as where is a square matrix of order ; that is, , where indices and vary from to . We use the substitution  , throughout our work. This substitution provides a way of calculating the connection coefficients. Here, and are integers and is scaling factor. It corresponds to either the expansion or the contraction of the scaling or wavelet function.
Define connection coefficients as
Similarly, we can obtain
It is homogeneous system and thus does not have a unique nonzero solution. In order to make the system inhomogeneous, one equation is added which is derived from the moment equation of the scaling function  as follows: where is the dth moment of and we can compute it by considering the orthonormality of that is
Considering the substitution in (13), we arrive at
Differentiate (15) times to get
Taking inner product on both sides of (16) with , or
Latto et al.  derive an explicit formula to compute the moments of . Consider where are Daubechies wavelet coefficients. Finally, we get the system for the calculation of connection coefficients. Consider where is a row vector with all the .
4. Implementation of Wavelet-Galerkin Method
Consider the following form of boundary value problem: where , , , and are real constants. A trial solution for (22) is
For simplicity, use substitution . Also, we have .
Now (24) implies
Multiplying , on both sides of (25), and integrating, we get where , is a polynomial of degree in .
The orthonormality of Daubechies wavelets implies
Treatment of the boundary conditions  is as follows.
Conditions and imply
The quasilinearization approach is a generalized Newton-Raphson technique for functional equations [26, 27]. It converges quadratically to the exact solution, if there is convergence at all, and it has monotone convergence.
Let us consider the nonlinear nth order differential equation as follows: with the boundary conditions
Here, is a function of the function . Choose an initial approximation of the function , let us say ; it may be , for . The function can now be expanded around the function by the use of the Taylor series. Consider where second and higher order terms are ignored. Using (34) in (32), we get
we obtain a third approximation for , call it . Assume that the problem converges and continue the procedure for obtaining desired accuracy. Recurrence relation is of the form where is known and can be used for obtaining . Equation (37) is always a linear differential equation and boundary conditions are
Now consider the nonlinear second order differential equation of the form 
Here, the first derivative can be considered as another function and (39) implies with the same boundary conditions
Similarly, one can follow the same procedure for higher order nonlinear differential equations to obtain the recurrence relation where is order of the differential equation. Equation (42) is always a linear differential equation and can be solved recursively, where is known and one can use it to get .
In order to test the wavelet-Galerkin method with quasilinearization technique, four different nonlinear problems are considered.
In this section, we solve some nonlinear ordinary differential equations by the wavelet-Galerkin method along with quasilinearization technique and compare the results with those obtained by other methods and exact solution.
Example 1. We consider the nonlinear two-point boundary value problem  as follows:
subject to the boundary conditions . The exact solution is . Apply quasilinearization to (43); we get
with the boundary conditions , .
Applying wavelet-Galerkin method to (44), we have and boundary conditions imply where , with the initial approximation .
We solved (43) by using and fixed the level of resolution . Figure 1 shows the exact solution and approximate solution by proposed method at first and second iteration; that is, is the initial approximation, and by using we get , that is, the solution at first iteration, and then is used to get , which is the solution of (43) at second iteration. Table 1 is used to compare the approximate solution by proposed method at second iteration for and exact solution. We may get more accurate results while increasing level of resolution at higher iteration.
Example 2. We consider the nonlinear Bratu’s boundary value problem as follows:
The quasilinear form of (47) is with the boundary conditions , .
Wavelet-Galerkin method for (48) implies and, from boundary conditions, we have with the initial approximation .
Bratu’s boundary value problem is solved by using as Galerkin bases and at level of resolution . Decomposition method  is already implemented on (47). We compared our results with the results obtained by decomposition method and exact solution. Our results are more accurate as compared to decomposition method  at and as shown in Tables 2 and 3, respectively. We used the MATLAB command of one-dimensional data interpolation using spline to get the values at and plot the exact and approximate solutions at these points for and as shown in Figure 2.
Example 3. Consider the nonlinear Troesch’s boundary value problem as follows:
The quasilinearized form of (51) is where , with the boundary conditions , .
Implementation of wavelet-Galerkin method to (52) implies and boundary conditions lead to with the initial approximation . Tables 4 and 5 and Figure 3 represent the solution of (51) at second iteration. We use as Galerkin bases to find the solution of (51) at different level of resolutions and at and as shown in Tables 4 and 5, respectively. Solutions by proposed method are compared with variational iteration method  and with exact solution. Our results are in high agreement with exact solution and better than variational iteration method .