Abstract

This paper is dedicated to proposing two numerical algorithms for solving the one- and two-dimensional heat partial differential equations (). In these algorithms, generalized Lucas polynomials () involving two parameters are utilized as basis functions. The two proposed numerical schemes in one and two- dimensions are based on solving the corresponding integral equation to the heat equation, and after that employing, respectively, the tau and collocation methods to convert the heat equations subject to their underlying conditions into systems of linear algebraic equations that can be treated efficiently via suitable numerical procedures. In this article, the convergence analysis is examined for the proposed generalized Lucas expansion. Five illustrative problems are numerically solved via the two proposed numerical schemes to show the applicability and accuracy of the presented algorithms. Our obtained results compare favourably with the exact solutions.

1. Introduction

Many mathematical models of real-world problems give rise to partial differential equations () of initial and boundary conditions. are frequently represented as mathematical equations that connect various amounts and their derivatives, e.g., heat transition, a particle’s movement in a straight line, the movement of a rocket, a molecule’s vibration, and a change in a substance’s molecular composition, etc. Every one of these issues is represented by hyperbolic, elliptic, or parabolic partial differential equation () and might be homogeneous, in one, two, or three dimensions, with non-local boundary conditions in addition to the initial conditions found in the prose. A parabolic is used to solve a variety of scientific problems, including ocean acoustic propagation as well as heat diffusion. The hyperbolic indicates the wave transformation and sound waves of an elastic string, whereas the elliptic describes the Laplace equation.

Fibonacci and Lucas polynomial sequences are crucial and they play vital roles in various disciplines. These sequences are employed to find approximate solutions of different types of DEs. For instance, Fibonacci polynomials were used to treat multi-term fractional DEs in [1]. In [2], Lucas polynomials are employed for the numerical treatment of sinh-Gordon equation. The authors in [3] developed a matrix method using Fibonacci polynomials for the treatment of the generalized pantograph equations with functional arguments. Another approach based on mixed Fibonacci and Lucas polynomials is followed in [4] to obtain numerical solutions of Sobolev equation in two dimensions. Lucas polynomials are employed in [5] to obtain numerical solutions of multidimensional Burgers-type equations. Lucas polynomials were also employed in [6] to solve the fractional-order electro-hydrodynamics flow model.

The Fibonacci and Lucas sequences can be generalized. For example, the authors in [7, 8] introduced two generalized families of Fibonacci and Lucas polynomials. In addition, they employed such generalized sequences to treat some fractional differential equations.

It is well-known that the heat equation is a parabolic that describes the distribution of heat. There are two types of heat equations: non-homogeneous and homogeneous. Non-homogeneous heat equations have source terms in the partial differential equations, whereas homogeneous heat equations do not have source terms. Many authors have researched theoretically and numerically the heat equations. For example, the authors in [9] obtained a numerical solution of the one-dimensional heat equation by using the Chebyshev wavelets method. In [10], the authors treated the same equation using a high-order compact boundary value method. The authors in [11] treated the heat equation using radial basis functions. In [12] a modified Crank-Nicolson scheme Richardson extrapolation is followed to treat the one-dimensional heat equation. Recently, the Chebyshev collocation algorithm is followed in [13] to treat the same equation.

A governs the temperature of a rod that is frequently defined as [14]: where is the temperature of a rod at position at time and is the thermal conductivity of the material, which measures the rod’s ability to conduct heat.

The solution’s domain is a semi-infinite wire of length that extends endlessly in time. In practice, the result is found only for a limited time. The solution with equation (1) necessitates the requirements of an initial condition at as well as boundary conditions at , and .

Initial condition:

Boundary conditions:

It is essential to refer here that (1) is called the homogeneous heat equation, whereas the non-homogeneous heat equation is given as: where is referred to as the heat source.

It is worth mentioning that the heat equation (1) governed by (2)-(4) can be extended to higher-dimensional heat equations. These types of equations were treated analytically and numerically by many authors. For example, the Adomian decomposition method was utilized for handling the two-dimensional heat equation in [14]. In addition, the collocation method together with the finite differences was employed to solve the same type of equations in [15]. Some analytical and numerical studies of a two-dimensional nonlinear heat equation with a source term were presented in [16]. Some other forms of the heat equations were handled in other contributions. For example, the authors in [17] applied the finite difference method of lines to treat the heat equation in three space variables. An Adomian decomposition method is applied to the treatment of a non-linear heat equation in [18]. For some other contributions relating to the heat equation, on can be referred to [1923].

There are numerous methods that have a significant impact on numerical analysis in general, see for example [2426]. Among these methods are the spectral methods, which play important roles in dealing with [27, 28], ordinary differential equations () [29, 30], and fractional differential equations () [3134]. The basic idea behind spectral methods is that the proposed approximate solution is written as linear combinations of many basic functions, which may be orthogonal or otherwise. The popular spectral approaches are Galerkin, collocation, and tau. In the context of numerical , each version has its own significance. Several authors have made extensive use of the latter methods. The Galerkin approach was followed to treat some types of differential equations. For example, the authors in [35] applied the Galerkin method to obtain spectral solutions of BVPs of even-orders, where the authors in [36] obtained approximate solutions of the fractional telegraph equation via implementing a spectral Legendre-Galerkin algorithm. Regarding the collocation method, it is an advantageous method from its capability for treating any type of differential equations governed by any underlying conditions. For example, it is followed in [37] to treat the initial value problems of any order with the aid of the operational matrices of some orthogonal polynomials. The tau method is different from the tau method in that no restrictions on choosing the basis functions. This of course makes its application to different types of is easier than the application of the Galerkin method. So, as a result, it is used for solving several types of differential equations.

The structure of this paper is as follows: Section 2 presents an overview of generalized Lucas polynomials and some of their fundamental properties. In Section 3, a numerical method based on the spectral tau method is applied to solve the one-dimensional partial differential heat equation. An extension to solve the two-dimensional heat equation is proposed in Section 4 based on the application of the collocation method. Section 5 examines the convergence and error analysis of the proposed expansion. Numerical outcomes and comparisons are presented in Section 6 to demonstrate the validity of our proposed methods. Section 7 is made up of a brief outline paper.

2. An Overview on Generalized Lucas Polynomials

The purpose of this section is to give an overview of the (). Furthermore, some of the basic formulas of these polynomials are presented.

The may be constructed with the aid of the following recursive formula:

They also may be generated by the following Binet’s formula:

The first few ones of the are given as follows:

It is important to point out that this kind of polynomials was employed in [8] to deal with some types of fractional .

Some celebrated polynomials can be obtained as special cases of the as a result of the existence of two parameters. In fact, the Lucas polynomials , Fermat-Lucas polynomials , Pell-Lucas polynomials , Chebyshev polynomials of the first kind , and Dickson polynomials of the first kind are special ones of the . Explicitly, we have

The have the following analytic formula ([8]): where [z] denotes the well-known floor function, which can also be written as: where

3. Numerical Treatment of the One-Dimensional Heat Equation

This section focuses on treating the one-dimensional partial differential heat equation We will analyze a numerical solution of the following one-dimensional linear non-homogeneous heat equation ([14]): governed by the non-homogeneous boundary conditions: and the initial conditions:

3.1. Integral Equation Corresponding to (13)-(15)

Our strategy to solve the one-dimensional heat equation (13) governed by the conditions (14) and (15) is to treat with its corresponding integral equation.

Now, integrating Eq. (13) with respect to the variable taking into the consideration the initial condition in (15), we get governed by the non-homogeneous boundary conditions:

3.2. Spectral Tau Treatment for the Heat Equation

The objective of the current section is to propose a spectral tau algorithm for numerically solving the corresponding integral form to the linear one-dimensional heat type equation. First, we consider the two families of basis functions and . Consider the next two spaces: where ; is the Sobolev space [38].

Now, the following approximation can be assumed for :

To use the spectral tau approach to Eq. (16) implies that we first compute the residual of Eq. (13). It is given by

The analytic form of in (11) allows us to express explicitly and in the following forms:

Based on the two Formulas (21) and (22), the residual in (20) can be rewritten as and therefore, the following system of equations can be acquired after the spectral tau technique is applied (see, [7]).

In addition, the boundary conditions (14) give:

Eqs. (24), (25), and (26) create a system of linear equations in the dimension with unknown expansion coefficients . The solution of this system can be found via the Gaussian elimination method.

4. Treatment of the Two-Dimensional Heat Equation

The distribution of heat flow in a two-dimensional space is governed by the following initial boundary value problem (see, [39, 40]) subject to the boundary conditions (BCs): and the initial condition (IC): where is the temperature of any point located at the position of a rectangular plate at any time , is the thermal diffusivity, and .

We suggest the following approximate spectral solution

By integrating (27) with respect to the variable t and making use of the IC (29), we get

Now, making use of (21) and (22), we can approximate in the form: where , can be expressed by (21), and can be expressed by (22).

Our strategy to solve numerically (27)-(29) is to utilize the spectral collocation method. For the residual of (31) is given by

We choose the following Riemann nodes , with

Hence, the application of the spectral collocation method implies that ([41]), and the use of the BCs leads to the following constraints:

Now, the above-mentioned equations build a system of algebraic equations of dimension , where

Thanks to the Gaussian elimination technique, we get the proposed approximate solution .

5. Error Analysis and Convergence of the Proposed GLPs Expansion

The goal of this section is to investigate the error analysis and convergence of the GLPs expansion that is used to solve the one-dimensional heat equation (13) governed by the underlying conditions (14) and (15). In the sequel, the next two lemmas are useful.

Lemma 1. Let and . For the , the following inequity is valid: where a and b are positive values.

Proof. We prove by mathematical induction. The inequality is satisfied for , since We now assume that (37) is satisfied for Finally, we demonstrate that validity of (37) for . Now, we have This ends the proof of Lemma 1.

Lemma 2. For all for every positive integer , and , the following inequity is valid for the :

Proof. By induction on , we will get started. Assume that the inequality (41) holds for and , and we have to prove that (41) itself holds. Now, our assumption implies that we have the following two inequalities: In virtue of the recurrence relation (6) and the Inequalities (42) and (43), we get Lemma 2 is now proved.

Theorem 3. let and belong to the space , and let , , . Let be expanded as We have the following: (1), which are positive constants.(2)The Series Comes to a Point of Absolute Convergence.

Proof. The first part of Theorem 3 can be demonstrated by following the same steps that were used in [8]). Now, we prove the remaining part of the theorem. Based on the first part, we have In virtue of Lemma 1, we get then the series comes to a point of absolute convergence.

Theorem 4. Let that belongs to the space satisfy the presumptions of Theorem 3, one obtains where the constants and are given as:

Proof. If we consider then, we have From Theorem 3, we get Based on Lemma 1, we can write and consequently, this leads to Because of simple inequity: ; .
Note:, , and denote, respectively, lower incomplete gamma, upper incomplete gamma, and gamma functions (see, [42]).

6. Numerical Outcomes and Comparisons

This section presents some examples to demonstrate the accuracy and performance of the following two proposed methods: (i)The generalized Lucas tau method () that employed for treating the one-dimensional heat equation.(ii)The generalized Lucas collocation method () that employed for treating the two-dimensional heat equation.

The error is represented by in the maximum norm, that is, in one dimension is computed by the formula:

We refer here that Mathematica software was used to perform all of the numerical data.

Example 5 (see [9, 12]). Consider the following heat equation:with the following initial boundary conditions:

The exact solution for Eq. (56) is: . In Table 1, the maximum absolute errors () obtained from the application of the are listed for and , while in Table 2, we compare the errors resulted from the application of the for the case corresponding to and with the best errors resulted from the application of the methods developed in [9, 12]. It is noticed from the obtained results in Table 2 that is more accurate than the two methods that developed in [9, 12]. Figure 1 displays the exact solution. Figure 2 shows the approximate solution for the case corresponding to a = b =1, whereas Figure 3 shows the resulting if the is applied.

Example 6 (see [9]). Consider the heat equation:governed by the following conditions:

The exact solution of Eq. (58) is: . A comparison of the of Example (55) resulting from the for the choices: , (), and (a = b =1) and the method applied in [9] is shown in Table 3. From the results in Table 3, it is evident that the is more accurate than the method developed in [9]. Furthermore, the exact solution, approximate solution (for the case corresponding to ), and absolute errors of the are displayed, respectively, in Figures 46.

Example 7 (see [9, 10]. The following homogeneous heat equation:governed by the following conditions:

The exact solution of (60) is: . For , and the three choices: , (a = b =1), and (a = 3, b =1), we compare the solutions behavior for GLTM and method in [9] as shown in Table 4. In Table 5, the for various values of M and a = b =1 is listed, which illustrates that the is more accurate than the method developed in [9]. Moreover, the approximate solution closes to the exact solution for and a = b =1 as shown in Table 6. For the three cases correspond to a = b =1, Figures 7 and 8 show the difference between the exact and the approximate solutions. Finally, Figure 9 plotted the absolute error when .

Example 8 (see [9]). Consider the following homogeneous heat equation:governed by the following conditions:

The exact solution of Eq. (62) is given by: . We present in Table 7 a comparison between the resulting error from the application of the for the case corresponding to , and for the three choices: , (a = b =1), and with those obtained by the application of the method presented in [9]. In Figures 1012, the exact solution, approximate solution, and absolute errors for the case corresponding to , resulted from the are, respectively, displayed.

Example 9 (see [39, 40]). Consider the following two-dimensional heat equation:with the boundary conditions: and the initial condition:

The exact solution of Eq. (64) is . For different times (t =1, t = 2, t = 3, and t = 4), we show the maximum pointwise error in Table 8 and the approximate solution in Figures 13, 14, 15, and 16, respectively.

7. Conclusions

In this paper, the generalized Lucas polynomials were utilized along with certain suitable spectral methods for obtaining numerical solutions of one- and two-dimensional heat equations. Two numerical approaches are followed for solving such equations. We showed that the proposed methods are superior if compared to some other methods. We have obtained more precise errors if the retained modes of the approximate expansions are small. Some estimations concerned with the generalized Lucas polynomials were proved and they served to investigate the convergence analysis of the suggested approximate expansion in one dimension. As future work, we plan to use the generalized Lucas polynomials to solve some other types of differential equations. In addition, we plan to use the generalized Lucas polynomials to solve some other types of heat equations.

Data Availability

No data is associated with this research.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Funding

No funding for this research.