A Jacobi Dual-Petrov-Galerkin Method for Solving Some Odd-Order Ordinary Differential Equations
E. H. Doha,1A. H. Bhrawy,2,3and R. M. Hafez4
Academic Editor: Simeon Reich
Received31 Oct 2010
Revised13 Feb 2011
Accepted16 Feb 2011
Published10 Apr 2011
Abstract
A Jacobi dual-Petrov-Galerkin (JDPG) method is introduced and used for solving fully integrated reformulations of third- and fifth-order ordinary differential equations (ODEs) with constant coefficients. The reformulated equation for the th order ODE involves -fold indefinite integrals for
. Extension of the JDPG for ODEs with polynomial coefficients is treated using the Jacobi-Gauss-Lobatto quadrature. Numerical results with comparisons are given to confirm the reliability of the proposed method for some constant and polynomial coefficients ODEs.
1. Introduction
A well-known advantage of spectral methods is high accuracy with relatively fewer unknowns when compared with low-order finite-difference methods [1, 2]. On the other hand, spectral methods typically give rise to full matrices, partially negating the gain in efficiency due to the fewer degrees of freedom. In general, the use of the Jacobi polynomials ( with and is the polynomial degree) has the advantage of obtaining solutions of ordinary differential equations (ODEs) in terms of the Jacobi indices (see for instance, [3–5]). Several such pairs have been used for approximate solutions of ODEs (see [6–10]). We avoid developing approximation results for each particular pair of indices and instead carry out a study with general indices. With this motivation, we introduce in this paper a family of the Jacobi polynomials with general indices.
Third-order differential equations have applications in many engineering models, see for instance [11–14]). Fifth-order differential equations generally arise in the mathematical modeling of viscoelastic flows and other branches of mathematical physics and engineering sciences, see [15–17]. Existence and uniqueness of solutions of such boundary value problems are discussed, for instance, in [18].
In this paper, the proposed differential equations are integrated times, where is the order of the equation, and we make use of the formulae relating the expansion coefficients of integration appearing in these integrated forms of the proposed differential equations to the Jacobi polynomials themselves (see, Doha [19]). An advantage of this approach is that the general equation in the algebraic system then contains a finite number of terms. We, therefore, motivated our interest in integrated forms of these differential equations. The interested reader is referred to Doha and Bhrawy [7, 20].
The main aim of this paper is to propose a suitable way to approximate some integrated forms of third- and fifth-order ODEs with constant coefficients using a spectral method, based on the Jacobi polynomials such that it can be implemented efficiently and at the same time has a good convergence property. It is worthy to note here that odd-order problems lack the symmetry of even-order ones, so we propose a Jacobi dual-Petrov-Galerkin (JDPG) method. The method leads to systems with specially structured matrices that can be efficiently inverted. We apply the method for solving the integrated forms of third- and fifth-order ODEs by using compact combinations of the Jacobi polynomials, which satisfy essentially all the underlying homogeneous boundary conditions. To be more precise, for the JDPG we choose the trial functions to satisfy the underlying boundary conditions of the differential equations, and we choose the test functions to satisfy the dual boundary conditions. Extension of the JDPG for polynomial coefficient ODEs is obtained by approximating the weighted inner products in the JDPG by using the Jacobi-Gauss-Lobatto quadrature. Finally, examples are given to illustrate the efficiency and implementation of the method. Comparisons are made to confirm the reliability of the method.
The remainder of this paper is organized as follows. In Section 2 we give an overview of the Jacobi polynomials and their relevant properties needed hereafter. Section 3 is devoted to the theoretical derivation of the JDPG method for third-order differential equations with constant and variable polynomial coefficients. Section 4 gives the corresponding results for those obtained in Section 3, but for the fifth-order differential equations with constant coefficient and two choices of boundary conditions. In Section 5, we present some numerical results exhibiting the accuracy and efficiency of our numerical algorithms. Some concluding remarks are given in the final section.
2. Preliminaries
Let be the space of polynomials of degree at most on the interval . We set
and let be the Jacobi polynomials orthogonal with the weight functions , where .
Let , , be the zeros of . Denote by , , the weights of the corresponding Gauss-Lobatto quadrature formula, which are arranged in decreasing order. We define the discrete inner product and norm of weighted space as follows:
Obviously (see, e.g., formula (2.25) of [21]),
Thus, for any , the norms and coincide.
For any , the Jacobi-Gauss-Lobatto interpolation operator , satisfying
We also denote by and the Chebyshev-Gauss-Lobatto and Legendre-Gauss-Lobatto interpolation operators, respectively.
For any real numbers , the set of the Jacobi polynomials forms a complete -orthogonal system, and
where is the Kronecker function and
The following special values will be of fundamental importance in what follows (see, [22–24])
where and .
If we define the times repeated differentiation and integration of by and , respectively, then (cf. Doha [19, 22])
where
with being a polynomial of degree at most . It is to be noted that may be obtained from by replacing with negative . In general, the hypergeometric series cannot be summed in explicit form, but it can be summed by Watsons identity [25], if . The following two lemmas will be of fundamental importance in what follows.
Proof. It is immediately obtained from relation (2.9).
3. Third-Order Differential Equation
We are interested in using the JDPG method to solve the third-order differential equation
subject to
where , , and are constants and is a given source function. In this paper, we consider the fully integrated form of the ODE, given by
where
In this work we assume that , the three-fold indefinite integral form of , can be evaluated analytically. We set
then the dual-Petrov-Galerkin approximation to (3.3) is to find such that
where and is the inner product in the weighted space . The norm in will be denoted by .
3.1. The Jacobi Dual-Petrov-Galerkin Method
We choose compact combinations of the Jacobi polynomials as basis functions to minimize the bandwidth hoping to improve the condition number of the coefficient matrix corresponding to (3.6). We choose the test basis and trial functions of expansions and to be of the form
where , , , , , and are the unique constants such that and . From the boundary conditions, and (2.7), hence , and can be uniquely determined by using mathematica to give (see, [27])
Using (2.7), and , one verifies readily that
Now it is clear that (3.6) is equivalent to
where is the discrete inner product associated with the Jacobi-Gauss-Lobatto quadrature. The constants , , and would not appear if we take in (3.11), therefore we get
If we take and as defined in (3.7) and (3.8), respectively, and if we denote
then
and the nonzero elements , , , and for are given as follows:
where
Hence by setting
then (3.12) is equivalent to the following matrix equation:
All the analytical formulae of the nonzero elements of matrices , , and can be obtained by direct computations using the properties of the Jacobi polynomials (for details see, [27, 28]).
In the case of , , we can form explicitly the LU factorization, that is, . In general, the expense of calculating LU factorization of an dense matrix is operations, and the expense of solving , provided that the factorization is known, is (see, [27]). However, in the case of banded matrix of bandwidth , we need just operations to factorize and operations to solve a linear system. In the case of , , the square matrix has bandwidth of 13. So we need just operations to solve the linear system (3.18). If , this represents a very substantial saving.
It is worthy to note that, for , the algebraic system (3.18) resulting from fully integrated reformulation of (3.1) is sparse and is therefore cheaper to solve than those obtained from the differentiated form (see [27, Theorem 3.1]). Moreover, the savings in computational effort increase as the size of the systems grow. Thus, we have demonstrated the advantage of using the integrated forms over the differentiated ones for constant coefficients ODEs.
3.2. A Quadrature JDPG Method
The JDPG can be extended for ODEs with polynomial coefficients because of analytical form of a product of an algebraic polynomial, and the Jacobi polynomials are known.
Now the formula of the Jacobi coefficients of the moments of one single Jacobi polynomial of any degree (see, Doha [22]) is
with , , where
This formula can be used to facilitate greatly the setting up of the algebraic systems to be obtained by applying the spectral methods for solving differential equations with polynomial coefficients of any order.
Let us consider the following integrated form of the third-order differential equation:
where , , and are the variable polynomial coefficients of the differential equation. The quadrature dual-Petrov-Galerkin method for (3.21) is to find such that
where is the discrete inner product of and associated with the Jacobi-Gauss-Lobatto quadrature.
Let us consider
and using Lemma 2.2 and formula (3.19), we can obtain
And by setting
then the linear system of (3.22) becomes
4. Fifth-Order Differential Equation
In this section, we consider the fifth-order differential equation of the form
but by considering its integrated form, namely,
4.1. First Choice of Boundary Conditions
Here, we apply the dual-Petrov-Galerkin approximation to (4.2) subject to the boundary conditions
We set
then the Jacobi dual-Petrov-Galerkin approximation to (4.2) is to find such that
We consider the following the Jacobi dual-Petrov-Galerkin procedure for (4.1): find such that
Now, we choose the basis and the dual basis functions and to be of the form
It is not difficult to show that the basis functions and the dual basis functions .
We choose the coefficients , , , , and such that verifies the boundary conditions (4.3). Making use of (2.7), then the boundary conditions (4.3) lead to linear system for these coefficients. The computation of the exact solution of such linear system for the unknown coefficients is extremely tedious by hand, and we have resorted to the symbolic computation software mathematica version 6, hence , , , , and can be uniquely determined to give
Since the dual basis functions satisfy the dual boundary conditions, and making use of (2.7) then the unknown coefficients , , , , and are determined by using Mathematica to give
It is clear that (4.6) is equivalent to
The constants , , , , and would not appear if we take in (4.10), therefore we get
If we take and as defined in (4.7) and if we denote
then
and the nonzero elements , , , and for are given as follows:
where
Then (4.11) is equivalent to the following matrix equation:
4.2. Second Choice of Boundary Conditions
In this subsection, we consider the fifth-order differential equation (4.1) with the following boundary conditions:
Equation (4.1) subject to the boundary conditions (4.18) has been considered in [29, 30]. Let us denote
then the Jacobi dual-Petrov-Galerkin approximation of (4.1) subject to (4.18) consists of finding such that
We consider the following choice of basis functions:
and dual basis functions:
where , , , , , , , , , and are chosen to be the unique constants such that and , for all .
Equation (4.20) is equivalent to the following matrix equation:
where the elements of the matrices , , , and can be obtained similarly as in the previous sections, but details are not given here.
5. Numerical Results
In this section some examples are considered aiming to illustrate how one can apply the proposed algorithms presented in the previous sections. Comparisons between JDPG method and other methods proposed in [29–32] are made.
Example 5.1. Consider the one-dimensional third-order equation
with an exact smooth solution
Table 1 lists the maximum pointwise error of , using the JDPG with various choices of , , and .
Example 5.2. Consider the one-dimensional fifth-order differential problem
with the exact solution .
Table 2 lists the maximum pointwise error, using the JDPG method with various choices of , , , , , and . Numerical results of this example show that the JDPG method converges exponentially.
Example 5.3. Consider the following fifth-order boundary value problem (see [29–32]):
The analytic solution of this problem is . Approximate solutions () are obtained by using our proposed method. Table 3 exhibits a comparison between the error obtained by using JDPG method and the sixth-degree B-spline [31], sextic spline [32], nonpolynomial sextic spline [29], and the computational method in [30]. The numerical results show that JDPG method is more accurate than the existing methods.
Example 5.4. Consider the following third-order ODE with polynomial coefficients:
subject to
where is selected such that exact solution is .
Equation (5.5) can be rearranged to take the form
and accordingly its fully integrated form is
Using the quadrature dual-Petrov-Galerkin method described in Section 3.2, we evaluate the maximum pointwise error of with various choices of , , and in Table 4. Numerical results show that there is a very good agreement between the approximate solution obtained by the quadrature JDPG method and the exact solution and at the same time ascertain that the JDPG method converges exponentially.
6. Concluding Remarks
In this paper, we described a JDPG method for fully integrated forms of third- and fifth-order ODEs with constant coefficients. Because of the constant coefficients, the matrix elements of the discrete operators are provided explicitly, and this in turn greatly simplifies the steps and the computational effort for obtaining solutions. However, the integrated form of the source function (involving severalfold indefinite integrals) should be known analytically, and the right hand side vector require quadrature approximations. This approach is also considered for ODEs with polynomial coefficients. Numerical results exhibit the high accuracy of the proposed numerical methods of solutions.
Acknowledgment
The authors are very grateful to the referees for carefully reading the paper and for their comments and suggestions which have improved the paper.
References
J. P. Boyd, Chebyshev and Fourier Spectral Methods, Dover, Mineola, NY, USA, 2nd edition, 2001.
C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Dynamics, Springer Series in Computational Physics, Springer, New York, NY, USA, 1988.
Y. Chen and T. Tang, “Convergence analysis of the Jacobi spectral-collocation methods for Volterra integral equations with a weakly singular kernel,” Mathematics of Computation, vol. 79, no. 269, pp. 147–167, 2010.
E. H. Doha and A. H. Bhrawy, “Efficient spectral-Galerkin algorithms for direct solution for second-order differential equations using Jacobi polynomials,” Numerical Algorithms, vol. 42, no. 2, pp. 137–164, 2006.
P. W. Livermore and G. R. Ierley, “Quasi- norm orthogonal Galerkin expansions in sums of Jacobi polynomials. Orthogonal expansions,” Numerical Algorithms, vol. 54, no. 4, pp. 533–569, 2010.
K. Aghigh, M. Masjed-Jamei, and M. Dehghan, “A survey on third and fourth kind of Chebyshev polynomials and their applications,” Applied Mathematics and Computation, vol. 199, no. 1, pp. 2–12, 2008.
E. H. Doha and A. H. Bhrawy, “Efficient spectral-Galerkin algorithms for direct solution of the integrated forms of second-order equations using ultraspherical polynomials,” The ANZIAM Journal, vol. 48, no. 3, pp. 361–386, 2007.
M. Fernandino, C. A. Dorao, and H. A. Jakobsen, “Jacobi galerkin spectral method for cylindrical and spherical geometries,” Chemical Engineering Science, vol. 62, no. 23, pp. 6777–6783, 2007.
A. R. Aftabizadeh, C. P. Gupta, and J.-M. Xu, “Existence and uniqueness theorems for three-point boundary value problems,” SIAM Journal on Mathematical Analysis, vol. 20, no. 3, pp. 716–726, 1989.
F. Bernis and L. A. Peletier, “Two problems from draining flows involving third-order ordinary differential equations,” SIAM Journal on Mathematical Analysis, vol. 27, no. 2, pp. 515–527, 1996.
A. Boucherif, S. M. Bouguima, N. Al-Malki, and Z. Benbouziane, “Third order differential equations with integral boundary conditions,” Nonlinear Analysis: Theory, Methods & Applications, vol. 71, no. 12, pp. e1736–e1743, 2009.
A. R. Davies, A. Karageorghis, and T. N. Phillips, “Spectral Glarkien methods for the primary two-point boundary-value problems in modeling viscoelastic flows, Internat,” International Journal for Numerical Methods in Engineering, vol. 26, pp. 647–662, 1988.
A. Karageorghis, T. N. Phillips, and A. R. Davies, “Spectral collocation methods for the primary two-point boundary-value problems in modeling viscoelastic flows,” International Journal for Numerical Methods in Engineering, vol. 26, no. 4, pp. 805–813, 1988.
G. L. Liu, “New research directions in singular perturbation theory: artificial parameter approach and inverse-perturbation technique,” in Proceedings of the 7th Conference on Modern Mathematics and Mechanics, Shanghai, China, 1997.
E. H. Doha, “Explicit formulae for the coefficients of integrated expansions of Jacobi polynomials and their integrals,” Integral Transforms and Special Functions, vol. 14, no. 1, pp. 69–86, 2003.
E. H. Doha and A. H. Bhrawy, “A Jacobi spectral Galerkin method for the integrated forms of fourth-order elliptic differential equations,” Numerical Methods for Partial Differential Equations, vol. 25, no. 3, pp. 712–739, 2009.
B.-Y. Guo and L.-l. Wang, “Jacobi interpolation approximations and their applications to singular differential equations,” Advances in Computational Mathematics, vol. 14, no. 3, pp. 227–276, 2001.
E. H. Doha, “On the construction of recurrence relations for the expansion and connection coefficients in series of Jacobi polynomials,” Journal of Physics A, vol. 37, no. 3, pp. 657–675, 2004.
E. H. Doha and A. H. Bhrawy, “Efficient spectral-Galerkin algorithms for direct solution of fourth-order differential equations using Jacobi polynomials,” Applied Numerical Mathematics, vol. 58, no. 8, pp. 1224–1244, 2008.
Y. L. Luke, Mathematical Functions and Their Approximations, Academic Press, New York, NY, USA, 1975.
E. H. Doha, A. H. Bhrawy, and R. M. Hafez, “A Jacobi-Jacobi dual-Petrov-Galerkin method for third- and fifth-order differential equations,” Mathematical and Computer Modelling, vol. 53, no. 9-10, pp. 1820–1832, 2011.
A. H. Bhrawy and S. I. El-Soubhy, “Jacobi spectral Galerkin method for the integrated forms of second-order differential equations,” Applied Mathematics and Computation, vol. 217, no. 6, pp. 2684–2697, 2010.
S. S. Siddiqi, G. Akram, and S. A. Malik, “Nonpolynomial sextic spline method for the solution along with convergence of linear special case fifth-order two-point value problems,” Applied Mathematics and Computation, vol. 190, no. 1, pp. 532–541, 2007.
X. Lv and M. Cui, “An efficient computational method for linear fifth-order two-point boundary value problems,” Journal of Computational and Applied Mathematics, vol. 234, no. 5, pp. 1551–1558, 2010.
H. N. Çaglar, S. H. Çaglar, and E. H. Twizell, “The numerical solution of fifth-order boundary value problems with sixth-degree -spline functions,” Applied Mathematics Letters, vol. 12, no. 5, pp. 25–30, 1999.
S. S. Siddiqi and G. Akram, “Sextic spline solutions of fifth order boundary value problems,” Applied Mathematics Letters, vol. 20, no. 5, pp. 591–597, 2007.