Abstract
A numerical method for onedimensional Bratu’s problem is presented in this work. The method is based on Chebyshev wavelets approximates. The operational matrix of derivative of Chebyshev wavelets is introduced. The matrix together with the collocation method are then utilized to transform the differential equation into a system of algebraic equations. Numerical examples are presented to verify the efficiency and accuracy of the proposed algorithm. The results reveal that the method is accurate and easy to implement.
1 Introduction
In this paper, we consider the boundaryvalue problem and initial value problem of Bratu’s problem. It is well known that Bratu’s boundary value problem in onedimensional planar coordinates is of the form
with the boundary conditions
where θ satisfies
The problem has zero, one or two solutions when
It was evaluated in [13] that the critical value
In addition, an initial value problem of Bratu’s problem
with the initial conditions
Bratu’s problem is also used in a large variety of applications such as the fuel ignition model of the thermal combustion theory, the model of thermal reaction process, the Chandrasekhar model of the expansion of the universe, questions in geometry and relativity about the Chandrasekhar model, chemical reaction theory, radiative heat transfer and nanotechnology [411].
A substantial amount of research work has been done for the study of Bratu’s problem. Boyd [2,12] employed Chebyshev polynomial expansions and the Gegenbauer as base functions. Syam and Hamdan [8] presented the Laplace decomposition method for solving Bratu’s problem. Also, Aksoy and Pakdemirli [13] developed a perturbation solution to Bratutype equations. Wazwaz [10] presented the Adomian decomposition method for solving Bratu’s problem. In addition, the applications of spline method, wavelet method and SincGalerkin method for solution of Bratu’s problem have been used by [1417].
In recent years, the wavelet applications in dealing with dynamic system problems, especially in solving differential equations with twopoint boundary value constraints have been discussed in many papers [4,16,18]. By transforming differential equations into algebraic equations, the solution may be found by determining the corresponding coefficients that satisfy the algebraic equations. Some efforts have been made to solve Bratu’s problem by using the wavelet collocation method [16].
In the present article, we apply the Chebyshev wavelets method to find the approximate solution of Bratu’s problem. The method is based on expanding the solution by Chebyshev wavelets with unknown coefficients. The properties of Chebyshev wavelets together with the collocation method are utilized to evaluate the unknown coefficients and then an approximate solution to (1) is identified.
2 Chebyshev wavelets and their properties
2.1 Wavelets and Chebyshev wavelets
In recent years, wavelets have been very successful in many science and engineering
fields. They constitute a family of functions constructed from dilation and translation
of a single function called the mother wavelet
Chebyshev wavelets
where
and
We should note that the set of Chebyshev wavelets is orthogonal with respect to the
weight function
The derivative of Chebyshev polynomials is a linear combination of lowerorder Chebyshev polynomials, in fact [20],
2.2 Function approximation
A function
where
where C and
and
3 Chebyshev wavelets operational matrix of derivative
In this section we first derive the operational matrix D of derivative which plays a great role in dealing with Bratu’s problem.
In the interval
Applying (6) the derivative of
The function
where
In fact we have shown that
where
From (10), it can be generalized for any
4 Solution of Bratu’s problem
Consider Bratu’s problem given in (1). In order to use Chebyshev wavelets, we first
approximate
Applying (11) we can get
Thus we have
We now collocate (12) at
Suitable collocation points are
Thus with the boundary conditions
Equations (13), (14) and (15) generate
5 Error analysis
Theorem 5.1A function
Since the truncated Chebyshev wavelets series is an approximate solution of Bratu’s
problem, so one has an error function
The error bound of the approximate solution by using Chebyshev wavelets series is given by the following theorem.
Theorem 5.2Suppose that
Proof Applying the definition of norm in the inner product space, we have
Because the interval
where
Therefore, using the above equation, we would get
□
6 Numerical examples
To illustrate the ability and reliability of the method for Bratu’s problem, some examples are provided. The results reveal that the method is very effective and simple.
Example 6.1 Consider the first case for Bratu’s equation as follows, when
We solve the equation by using the Chebyshev wavelets method with
Table 1. Computed absolute errors for Example 6.1
Example 6.2 Consider the initial value problem [10,1416,24]
The exact solution is
Applying (12) we get
Using the initial condition, we obtain
Equations (18) and (19) generate a system of nonlinear equations. These equations can be solved for unknown coefficients of the vector C. A comparison between the exact and the approximate solutions is demonstrated in Figure 1. From Figure 1, it can be found that the obtained approximate solutions are very close to the exact solution. In addition, Table 2 shows the exact and approximate solutions using the method presented in Section 3 and compares the results with the method presented in [16]. Also, by comparing the results of the table, we see that the results of the proposed method are more accurate.
7 Conclusions
The aim of present work is to develop an efficient and accurate method for solving Bratu’s problems. The Chebyshev wavelet operational matrix of derivative together with the collocation method are used to reduce the problem to the solution of nonlinear algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
CY completed the main study, carried out the results of this article and drafted the manuscript. JH checked the proofs and verified the calculation. All the authors read and approved the final manuscript.
Acknowledgements
Project is supported by the Huaihai Institute of Technology (No. Z2001151).
References

Ascher, UM, Matheij, R, Russell, RD: Numerical Solution of Boundary Value Problems for Ordinary Differential Equations, SIAM, Philadelphia (1995)

Boyd, JP: Chebyshev polynomial expansions for simultaneous approximation of two branches of a function with application to the onedimensional Bratu equation. Appl. Math. Comput.. 14, 189–200 (2003)

Buckmire, R: Investigations of nonstandard Mickenstype finitedifference schemes for singular boundary value problems in cylindrical or spherical coordinates. Numer. Methods Partial Differ. Equ.. 19(3), 380–398 (2003). Publisher Full Text

Hsiao, CH: Haar wavelet approach to linear stiff systems. Math. Comput. Simul.. 64, 561–567 (2004). Publisher Full Text

Buckmire, R: Application of a Mickens finitedifference scheme to the cylindrical BratuGelfand problem. Numer. Methods Partial Differ. Equ.. 20(3), 327–337 (2004). Publisher Full Text

McGough, JS: Numerical continuation and the Gelfand problem. Appl. Math. Comput.. 89, 225–239 (1998). Publisher Full Text

Mounim, AS, de Dormale, BM: From the fitting techniques to accurate schemes for the LiouvilleBratuGelfand problem. Numer. Methods Partial Differ. Equ.. 22(4), 761–775 (2006). Publisher Full Text

Syam, MI, Hamdan, A: An efficient method for solving Bratu equations. Appl. Math. Comput.. 176, 704–713 (2006). Publisher Full Text

Li, S, Liao, SJ: Analytic approach to solve multiple solutions of a strongly nonlinear problem. Appl. Math. Comput.. 169, 854–865 (2005). Publisher Full Text

Wazwaz, AM: Adomian decomposition method for a reliable treatment of the Bratutype equations. Appl. Math. Comput.. 166, 652–663 (2005). Publisher Full Text

He, JH: Some asymptotic methods for strongly nonlinear equations. Int. J. Mod. Phys. B. 20(10), 1141–1199 (2006). Publisher Full Text

Boyd, JP: Onepoint pseudo spectral collocation for the onedimensional Bratu equation. Appl. Math. Comput.. 217, 5553–5565 (2011). Publisher Full Text

Aksoy, Y, Pakdemirli, M: New perturbation iteration solutions for Bratutype equations. Comput. Math. Appl.. 59, 2802–2808 (2010). Publisher Full Text

Caglara, H, Caglarb, N, Özer, M: Bspline method for solving Bratu’s problem. Int. J. Comput. Math.. 87(8), 1885–1891 (2010). Publisher Full Text

Jalilian, R: Nonpolynomial spline method for solving Bratu’s problem. Comput. Phys. Commun.. 181, 1868–1872 (2010). Publisher Full Text

Venkatesh, SG, Ayyaswamy, SK, Raja Balachandar, S: The Legendre wavelet method for solving initial value problems of Bratutype. Comput. Math. Appl.. 63, 1287–1295 (2012). Publisher Full Text

Rashidinia, J, Maleknejad, K, Taheri, N: SincGalerkin method for numerical solution of the Bratu’s problem. Numer. Algorithms. 62, 1–11 (2013). Publisher Full Text

Lepik, U: Numerical solution of differential equations using Haar wavelets. Math. Comput. Simul.. 68(2), 127–143 (2005). Publisher Full Text

Daubechies, I: Ten Lectures on Wavelet, SIAM, Philadelphia (1992)

Sezer, M, Kaynak, M: Chebyshev polynomials solutions of linear differential equations. Int. J. Math. Educ. Sci. Technol.. 27(4), 607–618 (1996). Publisher Full Text

Adibi, H, Assari, P: Chebyshev wavelet method for numerical solution of Fredholm integral equations of the first kind. Math. Probl. Eng.. 2010, 1–17 (2010)

Dahlquist, G, Björck, A: Numerical Methods in Scientific Computing, SIAM, Philadelphia (2008)

Biazar, J, Ebrahimi, H: Chebyshev wavelets approach for nonlinear systems of Volterra integral equations. Comput. Math. Appl.. 63(3), 608–616 (2012). Publisher Full Text

Batiha, B: Nnmerical solution of Bratutype equations by the variational iteration method. Hacet. J. Math. Stat.. 39, 23–29 (2010)