Research Article  Open Access
An Efficient Numerical Approach for Solving Nonlinear Coupled Hyperbolic Partial Differential Equations with Nonlocal Conditions
Abstract
One of the most important advantages of collocation method is the possibility of dealing with nonlinear partial differential equations (PDEs) as well as PDEs with variable coefficients. A numerical solution based on a Jacobi collocation method is extended to solve nonlinear coupled hyperbolic PDEs with variable coefficients subject to initialboundary nonlocal conservation conditions. This approach, based on Jacobi polynomials and GaussLobatto quadrature integration, reduces solving the nonlinear coupled hyperbolic PDEs with variable coefficients to a system of nonlinear ordinary differential equation which is far easier to solve. In fact, we deal with initialboundary coupled hyperbolic PDEs with variable coefficients as well as initialnonlocal conditions. Using triangular, soliton, and exponentialtriangular solutions as exact solutions, the obtained results show that the proposed numerical algorithm is efficient and very accurate.
1. Introduction
For several decades, numerical methods have been developed to obtain more accurate solutions of differential and integral equations. Spectral method [1â€“9] is one of the family of weighted residual numerical methods for solving various problems, including variable coefficient and nonlinear differential equations [10, 11], integral equations [12, 13], integrodifferential equations [14, 15], fractional orders differential equations [16â€“19], and function approximation and variational problems [20]. The collocation method [21â€“28] can be classified as a special type of spectral methods. In the last few years, the collocation method has been introduced as a powerful approximation method for numerical solutions of all kinds of initialboundary value problems.
Exact solutions for initial value problem for some nonconservative hyperbolic systems are presented in [29], while the analytical study of variable coefficient mixed hyperbolic partial differential problems is discussed in [30]. The solitary and periodic wave solutions have been studied for some kinds of hyperbolic KleinGordon equations in [31, 32]. Other numerical methods based on the boundary integral equation [33] and numerical integration techniques [34] are used to numerically solve different types of hyperbolic partial differential problems. In [35, 36], finite difference scheme is considered to numerically solve hyperbolic equations. Pseudospectral methods are used in [37â€“40] to solve KleinGordon equations. In [41], Dehghan and Shokri used the radial basis functions to solve a twodimensional SineGordon equation; moreover in [42] they developed numerical scheme to solve the onedimensional nonlinear KleinGordon equation with quadratic and cubic nonlinearity using collocation points and approximating the solution using Thin Plate Splines and RBFs.
There are no results on JacobiGaussLobatto collocation (JGLC) method for solving nonlinear coupled hyperbolic PDEs with variable coefficients subject to initialboundary and nonlocal conditions. Therefore, the objective of this work is to present this method to numerically solve four nonlinear coupled hyperbolic PDEs with variable coefficients. By using collocation method, exponential convergence for the spatial variables can be achieved to approximate the solution of PDE. The computerized mathematical algorithm is the main key to apply this method for solving the problem. Moreover, the nonlocal conservation conditions are efficiently treated by JacobiGaussLobatto quadrature rule at nodes to obtain a system of ODEs in time and then proper initial value software can be applied to solve this system of ODEs. Several illustrative problems with various kinds of exact solutions such as triangular, soliton, and exponentialtriangular solutions are presented for demonstrating the high accuracy of this scheme. Moreover, with the freedom of selecting the Jacobi indexes and , the scheme can be calibrated for a wide variety of problems. Finally, the accuracy of the proposed method is demonstrated by solving some test nonlinear problems.
A brief outline of this paper is as follows. We present some properties of Jacobi polynomials in the next section. The third section is divided into two subsections: the first one deals with coupled nonlinear hyperbolic PDE with initialboundary conditions. The numerical treatment of solve initialnonlocal conservation conditions is developed in Section 3.2. In Section 4 the proposed method is applied to four different test problems to show the accuracy of our method. In the last section, we present some observations and conclusions.
2. Jacobi Polynomials
Some basic properties of Jacobi polynomials have been recalled in this section. A basic property of the Jacobi polynomials is that they are the eigenfunctions to a singular SturmLiouville problem: We recall that the Jacobi polynomials satisfy the following recurrence relation: where Beside the following relations Moreover, the th derivative of can be obtained from Let ; then we define the weighted space as usual. The inner product and the norm of with respect to the weight function are defined as follows: The set of Jacobi polynomials forms a complete orthogonal system, and
3. The Problem and the Numerical Algorithm
In this section, we approximate the solution of coupled nonlinear hyperbolic types equations with two different kinds of boundary conditions for space variable by using the Jacobi collocation method.
3.1. InitialBoundary Conditions
In what follows, we propose an efficient numerical algorithm to solve the coupled nonlinear hyperbolic types equations in the following form: related to the initial conditions, and the boundary conditions, Starting with the transformations . Problem (8)â€“(10) will be a new problem in the spatial variable . This transformation enable us to use the Jacobi collocation method on , subject to a new set of initial and boundary conditions,
Now, we are interested in using the JGLC method to transform the previous coupled PDEs into system of ODEs. In order to do this, we approximate the spatial variable using JGLC method at some nodal points. The node points are the set of points in a specified domain where the dependent variable values are approximated. In general, the choice of the location of the node points is optional, but taking the roots of the Jacobi orthogonal polynomials referred to as Jacobi collocation points gives particularly accurate solutions for the spectral methods. Now, we outline the main step of the JGLC method for solving couples hyperbolic problem. Let us expand the dependent variable in a Jacobi series, And, in virtue of (6)(7), we evaluate and by The JacobiGaussLobatto quadrature has been used to evaluate the previous integrals accurately. For any , we have that For any positive integer , stands for the set of polynomials of degree at most , () and () are used as the nodes and the corresponding Christoffel numbers in the interval , respectively. Thanks to (6), the coefficients in terms of the solution at the collocation points can be approximated by Due to (17), the approximate solution can be written as Furthermore, if we differentiate (18) once and evaluate it at the first JacobiGaussLobatto collocation points, it is easy to compute the first spatial partial derivative of the numerical solution in terms of the values at these collocation points as where Accordingly, one can obtain the second spatial partial derivative as where In the proposed JGLC method the residual of (11) is set to zero at of JacobiGaussLobatto points; moreover, the boundary conditions (13) will be enforced at the two collocation points and . Therefore, the approximation of (11)â€“(13) is where This approach provides a system of second order ODEs in the expansion coefficients , , with the following initial conditions: or in matrix notation as with where The system of second order (27)(28) can be solved by using diagonally implicit RungeKuttaNystrÃ¶m (DIRKN).
3.2. InitialNonlocal Conservation Conditions
Here, we will implement the JGLC algorithm for the coupled nonlinear hyperbolic type equations with nonlocal conditions: subject to the initial conditions, and the boundary conditions, while the other two boundary conditions have the nonlocal conservation form Again, we used the change of variables , to reduce problem (30)â€“(33) into related to the new initial conditions, the boundary conditions, and the nonlocal conservation conditions, The problem now is how to deal with the nonlocal conditions (37). For this purpose, let us introduce a collocation treatment for the integral conservation conditions (37) as The above equations may be rearranged as or briefly where Consequently, and are expressed as the following expansion of and ,â€‰â€‰: Based on the information included in this subsection and the recent one, we obtain the following system of ODEs: with the following initial conditions: where , â€‰â€‰, , and are given in (36) and (42).
4. Test Problems
We test the numerical accuracy of the proposed method by introducing four test problems with different types of exact solutions.
4.1. Triangular Solution
As a first example, we consider the coupled nonlinear hyperbolic equation (8) with the following functions: subject to The exact solutions of this problem are The absolute errors in the given tables are where and are the exact and approximate solutions at the point , respectively. Moreover, the maximum absolute error is given by The root mean square (RMS) and errors may be given by Maximum absolute, root mean square, and errors of (45) are introduced in Table 1 using JGLC method with three different choices of , and in the interval . The approximate solutions and of problem (45) have been plotted in Figures 1 and 2, with values of parameters listed in their captions. Moreover, we plot the curves of approximate and exact solutions of at different values of and in Figures 3 and 4. Again, the curves of approximate and exact solutions of at different values of and are displayed in Figures 5 and 6.

4.2. Soliton Solution
Secondly, consider the coupled nonlinear hyperbolic equation (8) with the following functions: subject to The exact solutions are Table 2 shows the accurate results for maximum absolute, root mean square, and errors of (51) for various choices of , and in the interval . Figures 7 and 8 show that the absolute errors and are very small with values of parameters listed in their captions. We also plot the curves of approximate and exact solutions of at different values of and in Figures 9 and 10. Moreover, in Figures 11 and 12, the approximate and exact solutions of are plotted at different values of and .