Next Article in Journal
Using Lie Sphere Geometry to Study Dupin Hypersurfaces in Rn
Next Article in Special Issue
Stability and Convergence Analysis of the Discrete Dynamical System for Simulating a Moving Bed
Previous Article in Journal
Statistical Advancement of a Flexible Unitary Distribution and Its Applications
Previous Article in Special Issue
Optimal Fourth-Order Methods for Multiple Zeros: Design, Convergence Analysis and Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Homotopy Analysis Transform Method for a Singular Nonlinear Second-Order Hyperbolic Pseudo-Differential Equation

by
Said Mesloub
* and
Hassan Eltayeb Gadain
Department of Mathematics, College of Science, King Saud University, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(6), 398; https://doi.org/10.3390/axioms13060398
Submission received: 4 May 2024 / Revised: 7 June 2024 / Accepted: 8 June 2024 / Published: 14 June 2024
(This article belongs to the Special Issue Applied Mathematics and Numerical Analysis: Theory and Applications)

Abstract

:
In this study, we employed the homotopy analysis transform method (HATM) to derive an iterative scheme to numerically solve a singular second-order hyperbolic pseudo-differential equation. We evaluated the effectiveness of the derived scheme in solving both linear and nonlinear equations of similar nature through a series of illustrative examples. The stability of this scheme in handling the approximate solutions of these examples was studied graphically and numerically. A comparative analysis with existing methodologies from the literature was conducted to assess the performance of the proposed approach. Our findings demonstrate that the HATM-based method offers notable efficiency, accuracy, and ease of implementation when compared to the alternative technique considered in this study.

1. Introduction

Nonlinear partial differential equations serve as fundamental tools for describing a diverse array of phenomena across numerous scientific disciplines, including physics, biology, chemistry, economics, epidemiology, engineering, and beyond. Despite their ubiquity, exact solutions to the majority of these equations remain elusive. Consequently, various computational methodologies have been devised by researchers to address this challenge. Notable among these techniques are the Adomian decomposition method pioneered by G. Adomian [1,2,3,4], the variational iteration method introduced by J. He [5,6], the finite difference method developed by Crank–Nicolson in [7], the homotopy analysis method (HAM) proposed by Liao [8,9,10,11,12], the homotopy perturbation method advocated by J. He [13,14,15], and the differential transform method utilized by several researchers [16,17].
In recent years, a number of researchers have focused on investigating solutions to these partial differential equations through the integration of diverse techniques with the Laplace transform. Notable examples include the Laplace decomposition method detailed in references [18,19], the homotopy perturbation transform method explored in [20], and the introduction of the homotopy analysis transform method (HATM) as documented in [21]. Since then, the HATM has been widely implemented by many authors to solve various types of mathematical equations. For example, it was employed in [22] to solve a fractional equation modeling gas dynamics, and it was utilized in [23] to solve a fractional Fokker–Planck equation. In [24], it is applied to solve an integro-differential equation; in [25], the HATM was used to obtain exact solutions for fractional equations of type Korteweg–de Vries and Korteweg–de Vries–Burger’s; in [26], it was implemented to solve fractional Navier–Stokes equations; in [27], a modified version of the HATM was utilized to numerically solve a fractional-order smoking epidemic model; in [28], the HATM was employed to solve a nonlinear integro-differential equation; etc.
The goal of this study was to derive a numerical solution for a singular nonlinear second-order hyperbolic pseudo-differential equation based on the HATM. Moreover, the effectiveness of the proposed technique was explored through a set of illustrative example as well as a comparative analysis with existing methodologies in the literature. In particular, we employed the HATM to numerically solve the following pseudo-hyperbolic initial value problem:
2 t 2 U ( x , t ) α x x ( x U x ) β x 2 x t ( x U x ) γ x 2 x ( U U x ) + λ F ( U ) = G ( x , t ) ,
subject to the constraints
U ( x , 0 ) = g 1 ( x ) , x U ( x , 0 ) = g 2 ( x ) ,
where α , β , γ , and λ are constants, and F ,   G ,   g 1 , and g 2 are known functions.
In fact, at β = γ = λ = 0 , Equation (1) reduces to a singular one-dimensional wave equation, and at α = γ = λ = 0 , it reduces to a singular one-dimensional pseudo-wave equation.
The nonlinear singular one-dimensional pseudo-hyperbolic equation is a type of PDE that involves mixed partial derivatives in space and time, and it exhibits both hyperbolic and singular characteristics. In general, these equations are often used to model complex physical phenomena. In our model, the equation describes how waves propagate, attenuate, and interact in a medium with radial symmetry, such as sound waves in a cylindrical duct or stress waves in a cylindrical rod. The nonlinear terms in the equation represent complex interactions, like nonlinear elastic responses or external forces acting on the medium. These nonlinear terms can lead to complex behaviors such as wave steepening, shock formation, or solitons. They often represent interactions within the medium or external forces that depend on the wave’s state. More precisely, the equation models a radial wave propagation in an elastic medium, where the Bessel operator reflects the radial nature of wave propagation.
To achieve our goal, we need the Laplace transform for a second-order partial derivative of the function U with respect to t which is given by
L 2 t 2 U ( x , t ) = s 2 U ( x , s ) s   U ( x , 0 ) t U ( x , 0 ) .
Let us mention that linear and nonlinear versions of Equation (1) were considered in [29], where a combination of the double Laplace transform and the Adomian decomposition method was applied to solve these versions.
The rest of this paper is organized as follows. In Section 2, we provide the basic outlines of the HATM. In Section 3, the HATM is employed to derive the numerical method for solving model (1). Section 4 is devoted to the application of the derived iterative scheme to a set of test examples to evaluate its efficiency. Finally, in Section 5, we comment on our findings and conclusions.

2. Basic Outline of HATM

The homotopy analysis transform method was first introduced in [21]. It is an elegant combination of the Laplace transform method and the homotopy analysis method. Its implementation does not require discretization or any other restrictive assumptions, which avoid round-off errors that may affect the reach of a closed-form solution of the problem. Moreover, it does not require decomposing the nonlinear terms in the equation under consideration, which usually leads to tedious algebraic work. Additionally, for some strongly nonlinear problems other techniques may fail to produce a convergent series solution, while the existence of the parameter in theh HATM can be used to adjust and control the convergence region to ensure the convergence of a series solution for such problems. This reveals that the HATM is a very effective and accurate method in handling the solution of various nonlinear types of partial differential equations. Thus, our motivation for this work was to determine a numerical solution for the singular problem (1) based on the HATM that is easy to implement and does not require any restrictive assumptions needed for the implementation of the aforementioned methods. To formulate the basic outlines of the HATM, we consider a general partial differential equation as follows:
t 2 U ( x , t ) + R U ( x , t ) + N U ( x , t ) = g ( x , t ) ,
subject to the initial conditions
U ( x , 0 ) = f 1 ( x )   and   U t ( x , 0 ) = f 2 ( x ) ,
where R represents a linear differential operator, N represents a nonlinear differential operator, and g ,   f 1 , and f 2 are known functions.
Taking the Laplace transform L for both sides of Equation (3) implies that
s 2 L U ( x , t ) s U ( x , 0 ) U t ( x , 0 ) + L R U ( x , t ) + N U ( x , t ) = L g ( x , t ) ,
that is,
L U ( x , t ) 1 s U ( x , 0 ) 1 s 2 U t ( x , 0 ) + 1 s 2 L R U ( x , t ) + N U ( x , t ) g ( x , t ) = 0 .
Thus, according to the theory of the HAM [8], we can define an operator N as follows:
N [ ϕ ( x , t ; q ) ] = L ϕ ( x , t ; q ) 1 s ϕ ( x , 0 ; q ) 1 s 2 ϕ t ( x , 0 ; q ) + 1 s 2 L { R ϕ ( x , t ; q ) + N ϕ ( x , t ; q ) g ( x , t ) } ,
where q [ 0 , 1 ] , and ϕ is a real-valued function in x , t , and q.
Then, we can define the zeroth-order deformation equation as follows:
( 1 q ) L ϕ ( x , t ; q ) U 0 ( x , t ) = q N [ ϕ ( x , t ; q ) ] ,
where q [ 0 , 1 ] is an embedding parameter; L is the Laplace transform operator; U 0 ( x , t ) is an initial approximation of the exact solution U ( x , t ) ; ϕ ( x , t ; q ) is an unknown function of x ,   t , and q; and is a nonzero auxiliary parameter, which plays a fundamental role in the convergence of the series solution.
It follows from Equation (4) that ϕ ( x , t ; 0 ) = U 0 ( x , t ) at q = 0 and ϕ ( x , t ; 1 ) = U ( x , t ) at q = 1 . Thus, as q moves continuously from 0 to 1, the function ϕ ( x , t ; q ) deforms from the initial guess U 0 ( x , t ) toward the exact solution U ( x , t ) .
Next, the Taylor series expansion of ϕ ( x , t ; q ) in powers of q yields
ϕ ( x , t ; q ) = U 0 ( x , t ) + i = 1 U i ( x , t )   q i ,
where
U i ( x , t ) = 1 i ! i ϕ ( x , t ; q ) q i q = 0 .
As Liao pointed out in [11], if the initial guess U 0 ( x , t ) and the parameter are chosen properly, then power series (5) will converge at q = 1 to a solution of the considered partial differential equation, and this solution is expressed in series form as follows:
U ( x , t ) = U 0 ( x , t ) + i = 1 q i U i ( x , t ) .
In fact, the success of the convergence of the series solution in the HATM depends on three elements, namely, the choice of the operator to be inverted, the value of the parameter , and the initial guess. The choice of the suitable operator takes into account that the operator is easily inverted. The existence of the parameter distinguishes the HATM from other semi-analytical methods, as its occurrence in the series solution enables one to adjust the convergence interval as required. But in the absence of easy theoretical techniques that can be implemented to determine its value, an alternative graphical technique based on the h-curve resulting from the plots of the approximate truncated series solutions is used to determine a range of values for this parameter, usually the interval on which the curve is almost horizontal. Finally, the choice of the initial guess depends on the operator to be inverted. It can be chosen from the initial or boundary conditions in the problem according to the inverted operator, and it should satisfy the differential equation under consideration.
Now, differentiating Equation (4) k times with respect to q, setting q = 0 , and dividing the result by k ! produce the following kth-order deformation equation:
L U k ( x , t ) χ k U k 1 ( x , t ) = R ( U k 1 ) ,
where
U k ( x , t ) = [ U 0 ( x , t ) , U 1 ( x , t ) , , U k ( x , t ) ] ,
and
R ( U k 1 ) = 1 ( k 1 ) ! k 1 q k 1 N ϕ ( x , t ; q ) q = 0 .
Finally, taking the inverse Laplace transform for both sides of Equation (6) implies that the components U k ( x , t ) ,   k = 1 , 2 , , can be determined by applying the following iterative formula:
U k ( x , t ) = χ k U k 1 ( x , t ) + L 1 R ( U k 1 ) ,
where
χ k = 0 ,     k 1 , 1 ,     k   >   1 .
Let us mention that the accuracy of the numerical solution will be improved as the number of terms in the truncated series solution is increased. On the other hand, the computational process can be terminated when either of the computed components lead to the closed-form solution or the terminated series solution provide the required accuracy of the numerical solution.

3. Application of the Method

To explore the applicability and efficiency of the HATM for solving problem (1) subject to initial conditions
U ( x , 0 ) = g 1 ( x ) , x U ( x , 0 ) = g 2 ( x ) ,
where g 1 , and g 2 are known functions, we apply the Laplace transform to both sides of this equation to obtain
L [ U ( x , t ) ] 1 s   g 1 ( x ) 1 s 2   g 2 ( x ) 1 s 2   L α x ( x U x ) x + β x ( x U x ) t x + γ x 2 ( U U x ) x λ F ( U ) + G ( x , t ) = 0 .
Thus, we define an operator N [ ϕ ( x , t ; q ) ] as follows:
N [ ϕ ( x , t ; q ) ] = L ϕ ( x , t ; q ) 1 s g 1 ( x ) 1 s 2 g 2 ( x ) 1 s 2 L [ α x ( x ϕ x ( x , t ; q ) ) x + β x ( x ϕ x ( x , t ; q ) ) t x + γ x 2 ( ϕ ( x , t ; q ) ϕ x ( x , t ; q ) ) x λ F ( ϕ ( x , t ; q ) ) + G ( x , t ) ] .
Therefore, the kth-order deformation equation can be defined as follows:
L [ U k ( x , t ) χ k U k 1 ( x , t ) ] = R ( U k 1 ) ,
where
R ( U k 1 ) = L U k 1 ( x , t ) 1 χ k 1 s g 1 ( x ) + 1 s 2 g 2 ( x ) 1 s 2 L [ α x x ( x x U k 1 ( x , t ) ) + β x 2 x t ( x x U k 1 ( x , t ) ) + γ x 2 x ( U k 1 ( x , t ) x U k 1 ( x , t ) ) λ F ( U k 1 ( x , t ) ) + 1 χ k G ( x , t ) ] .
Hence, the components of the series solution can be computed recursively from the iterative formula
U k ( x , t ) = χ k U k 1 ( x , t ) + L 1 R ( U k 1 ) ,     k = 1 , 2 , ,
and the solution is given by
U ( x , t ) = U 0 ( x , t ) + k = 1   U k ( x , t ) .

4. Numerical Results

In this section, we test the efficiency and accuracy of scheme (7) in handling the numerical solutions for problems of type (1). It is notable to mention that all computations in this section were conducted using the software Mathematica 8.04.0., and all the numerical results in the tables were computed using the default machine precision and accuracy of 17 significant digits rounded to 6 significant digits. Now, we apply this scheme to the following set of examples.
Example 1.
Consider the singular linear partial differential equation
2 t 2 U ( x , t ) 1 x x ( x U x ) 1 x 2 x t ( x U x ) = G ( x , t ) ,
subject to the following initial constraints:
U ( x , 0 ) = 0   a n d U t ( x , 0 ) = x 2 ,
where G ( x , t ) = x 2 sin t 4 sin t 4 cos t .
Solution.
Let U 0 ( x , t ) = U ( x , 0 ) = x 2 . Then, in view of (7), we obtain
U 1 ( x , t ) = χ 1 U 0 ( x , t ) +   L 1 R ( U 0 ) =   L 1 R ( U 0 ) = 4 ( 1 + t ) + 4 cos t + 4 + x 2 sin t ,
U 2 ( x , t ) = χ 2 U 1 ( x , t ) +   L 1 R ( U 1 ) = U 1 ( x , t ) +   L 1 R ( U 1 ) = 4 ( 1 + 2 ) ( 1 + t ) + ( 4 + 8 ) cos t + 4 + x 2 + 8 + x 2 sin t ,
U 3 ( x , t ) = χ 3 U 2 ( x , t ) +   L 1 R ( U 2 ) = U 2 ( x , t ) +   L 1 R ( U 2 ) = ( 1 + ) 4 ( 1 + 3 ) ( 1 + t ) + 4 ( 1 + 3 ) cos t + 4 + x 2 + 12 + x 2 sin t .
Thus, the series solution is given by
U ( x , t ) = U 0 ( x , t ) + U 1 ( x , t ) + U 3 ( x , t ) + = 12 ( 1 + ) 2 1 + t cos t 3 4 + x 2 + 3 8 + x 2 + 2 12 + x 2 sin t + .
It is noticed that at = 1 , the components of the series solution (8) are reduced to
U 1 ( x , t ) = 4 ( 1 + t ) + 4 cos t + 4 + x 2 sin t .
U 2 ( x , t ) = 4 ( 1 t + cos t + sin t ) .
U 3 ( x , t ) = 0 .
Continuing in this manner, we obtain U j ( x , t ) = 0 for all j 3 . Hence, the series solution (8) reduces to
U ( x , t ) = U 0 ( x , t ) + U 1 ( x , t ) + U 3 ( x , t ) + = x 2 sin t ,
which is the exact solution for this equation given in [29].
Figure 1 shows the h-curve corresponding to the 13th-order series solution. This figure shows that the values of that lead to the convergence of the series solution (9) should be selected from the range 1.9 < < 0.2 .
Figure 2 shows plots of the truncated series solutions for Example 1 with a different number of terms at x = 4 and = 1.2 . This figure shows the rapid convergence of these approximate solutions toward the exact solution. It illustrates that the curves of the truncated series solutions are getting closer to the exact solution as the number of terms in the truncated series solution is increased.
Table 1 and Table 2 show the absolute error E r [ i ] = | U ( x , t ) U [ i ] ( x , t ) | resulting from taking the difference between the exact solution and the approximate solution of several orders of Example 1 at different values of , t, and x. These data illustrate the rapid convergence of the approximate solution just after a few terms. Moreover, it appears from an error analysis of these tables that the convergence of the approximate solutions is more rapid at = 1 than at the other values of in the interval of convergence.
Example 2.
Consider the singular nonlinear partial differential equation
2 t 2 U ( x , t ) 1 x x ( x U x ) 1 x 2 x t ( x U x ) x 2   U U x + U 2 = x 2   e t ,
subject to the following initial conditions:
U ( x , 0 ) = x 2   a n d U t ( x , 0 ) = x 2 .
Solution.
Let U 0 ( x , t ) = x 2   e t . Then, in view of (7), we obtain
U 1 ( x , t ) = χ 1 U 0 ( x , t ) +   L 1 R ( U 0 ) =   L 1 R ( U 0 ) = 0 ,
U 2 ( x , t ) = χ 2 U 1 ( x , t ) +   L 1 R ( U 1 ) = U 1 ( x , t ) +   L 1 R ( U 1 ) = 0 ,
Continuing in this manner, we obtain U j ( x , t ) = 0 for all j 1 . Thus, the series solution is given by
U ( x , t ) = U 0 ( x , t ) + U 1 ( x , t ) + U 3 ( x , t ) + = x 2   e t ,
which is the exact solution for this equation presented in [29], without using Adomain polynomials.
Example 3.
Consider the singular nonlinear partial differential equation
2 t 2 U ( x , t ) 1 x x ( x U x ) 1 x 2 x t ( x U x ) x 2   x U U x = ln x t 1 ,
subject to the following initial conditions:
U ( x , 0 ) = ln x   a n d U t ( x , 0 ) = 1 .
Solution.
Let U 0 ( x , t ) = U ( x , 0 ) = ln x . Then, in view of (7), we obtain
U 1 ( x , t ) = 1 6 t 6 + t 2 ,
U 2 ( x , t ) = 1 6 t 6 + t 2 1 120 2 t 120 40 t 2 + t 4 ,
U 3 ( x , t ) = 1 6 t 6 + t 2 1 60 2 t 120 40 t 2 + t 4 + 3 t 5040 + 2520 t 2 126 t 4 + t 6 5040 ,
Thus, the series solution is given by
U ( x , t ) = U 0 ( x , t ) + U 1 ( x , t ) + U 3 ( x , t ) + = ln x + t 2 6 + t 2 2 t 40 120 40 t 2 + t 4 + 3 t 5040 5040 + 2520 t 2 126 t 4 + t 6 + .
In fact, at = 1 , the components of the series solution (9) are reduced to
U 1 ( x , t ) = t 1 t 2 6 ,
U 2 ( x , t ) = 1 120 t 3 20 t 2 ,
U 3 ( x , t ) = t 5 42 t 2 5040 ,
In general, we have
U j ( x , t ) = t 2 j 1 ( 2 j + 1 ) ! 2 j ( 2 j + 1 ) t 2 .
Hence, series solution (9) reduces to
U ( x , t ) = ln x + j = 1 t 2 j 1 ( 2 j + 1 ) ! 2 j ( 2 j + 1 ) t 2 .
Figure 3 shows the h-curve corresponding to the 10th-order series solution. This figure shows that the values of that lead to the convergence of series solution (9) should be selected from the interval 1.7 < < 0.4 .
Figure 4 shows plots of the truncated series solutions for Example 3 with a different number of terms at x = 0.1 and = 0.9 . This figure shows the rapid convergence of these approximate solutions toward the exact solution. It illustrates that the curves of the truncated series solutions are getting closer to the exact solution as the number of terms in the truncated series solution is increased.
The exact solution of Example 3 is U ( x , t ) = t ln x . Table 3 and Table 4 show the absolute error E r [ i ] = | U ( x , t ) U [ i ] ( x , t ) | resulting from taking the difference between the exact solution and the approximate solution of several orders at different values of , t, and x. These data demonstrate the rapid convergence of the approximate solution just after a few terms. Moreover, it appears from an error analysis of these tables that the convergence of the approximate solutions is more rapid at = 1 than at the other values of the parameter in the interval of convergence.
It is seen from the tables and figures in the above examples that the error is decreased and the accuracy of the approximate solution is refined as the number of iterations is increased, which reflects the reliability of the approximate solutions obtained using the HATM. On the other hand, these examples illustrate the stability of the derived iterative scheme (7) numerically and graphically. In both cases, this scheme generates iterations used to construct the approximate series solutions. The plots in Figure 2 and Figure 4 show the behavior of the first few terms in the sequence of partial sums of these series solutions, U [ n ] ( x , t ) = i = 0 n U i ( x , t ) . These plots display that the graph of U [ n ] ( x , t ) is getting closer and closer to the graph of the exact solution as n increases, which means that the error in each step is decreased and the global error is bounded, i.e., the scheme is stable, and the sequence of partial sums converges to the analytical solution. In addition to this, the numerical results in Table 1, Table 2, Table 3 and Table 4 show that the absolute error in the computed values of these partial sums at selected values of x and t decreases rapidly as n increases. This occurs for distinct values of the parameter in the interval of convergence. Again, this means that the scheme is stable; otherwise, the global error increases without a bound as n increases.

5. Conclusions

In this study, an iterative scheme was obtained to find a numerical solution for a singular nonlinear second-order hyperbolic pseudo-differential equation based on the HATM. A series of three examples was given to demonstrate the accuracy and efficiency of the obtained iterative scheme. The error analysis conducted for the computed numerical solutions of these examples was discussed and illustrated both graphically and numerically. The plots in Figure 2 and Figure 4 show the behavior of the first few truncated series solutions. These plots are getting closer and closer to the graph of the exact solution as the number of terms in these series increases, which means that the error is reduced in each step, i.e., the derived scheme is stable, and the approximate series converges to the analytical solution. Also, the results in Table 1, Table 2, Table 3 and Table 4 show that the computed absolute error is inversely proportional to the number of terms in the series solution, which shows that the numerical solutions obtained using the HATM are accurate and reliable. In addition, it is noticed from these tables that the best performance of this method is obtained at = 1 . On the other hand, the numerical solutions obtained using this method were compared with those obtained using the double Laplace transform given in [29]. It was found that the results obtained using the two methods coincide. But the implementation of the HATM avoids the use of Adomian polynomials required for the application of the other method, as these polynomials require tedious computational work.
Thus, we conclude that the HATM is an efficient method and can be implemented easily for solving this type of problems without any restrictive assumptions. Motivated by these results, we look forward to considering a numerical solution for a two-dimensional singular nonlinear hyperbolic equation based on the HATM in our future work. Finally, we would like to mention that the authors contributed equally to this work.

Author Contributions

Conceptualization, S.M.; methodology, S.M., Validation, H.E.G., investigation, H.E.G.; Writing original draft preparation, S.M.; Writing—review and editing, S.M. All authors have read and agreed to the published version of the manuscript.

Funding

Researchers Supporting Project (RSPD2024R975) King Saud University, Riyadh, Saudi Arabia.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors extend their appreciation to the Researchers Supporting Project (RSPD2024R975) King Saudi Arabia, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Adomian, G. Nonlinear Stochastic Operator Equations; Kluwer Academic Publishers: London, UK, 1986; ISBN 978-0-12-044375-8. [Google Scholar]
  2. Adomian, G. A Review of the Decomposition Method in Applied Mathematics. J. Math. Anal. Appl. 1988, 135, 501–544. [Google Scholar] [CrossRef]
  3. Adomian, G. A Review of The Decomposition Method and Some Recent Results for Nonlinear Equations. Comput. Math. Appl. 1991, 21, 101–127. [Google Scholar] [CrossRef]
  4. Abbaoui, K.; Cherruault, Y. Convergence of Adomian’s Method Applied to Nonlinear Equations. Math. Comput. Model. 1994, 20, 69–73. [Google Scholar] [CrossRef]
  5. He, J.H. Variational iteration method a kind of non-linear analytical technique: Some examples. Int. J. Non-Linear Mech. 1999, 34, 699–708. [Google Scholar] [CrossRef]
  6. He, J.H. Variational iteration method for autonomous ordinary differential systems. Appl. Math. Comput. 2000, 114, 115–123. [Google Scholar] [CrossRef]
  7. Crank, J.; Nicolson, P. A practical method for numerical evaluation of solutions of partial differential equations of the heat conduction type. Proc. Camb. Philos. Soc. 1947, 43, 50–67. [Google Scholar] [CrossRef]
  8. Liao, S.J. The Proposed Homotopy Analysis Technique for the Solution of Nonlinear Problems. Ph.D. Thesis, Shanghai Jiao Tong University, Shanghai, China, 1992. [Google Scholar]
  9. Liao, S.J. Homotopy analysis method a new analytical technique for nonlinear problems. Commun. Nonlinear Sci. Numer. Simul. 1997, 2, 95–100. [Google Scholar] [CrossRef]
  10. Liao, S.J. An explicit, totally analytic approximation of Blasius’ viscous flow problems. Int. J. Non-Linear Mech. 1999, 34, 759–778. [Google Scholar] [CrossRef]
  11. Liao, S.J. Beyond Perturbation: Introduction to the Homotopy Analysis Method; Chapman & Hall/CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  12. Liao, S.J. On the homotopy analysis method for nonlinear problems. Appl. Math. Comput. 2004, 147, 499–513. [Google Scholar] [CrossRef]
  13. He, J.H. Homotopy perturbation technique. Comptut. Methods Appl. Mech. Engrg. 1999, 178, 257–262. [Google Scholar] [CrossRef]
  14. He, J.H. A coupling method of a homotopy technique and a perturbation technique for non-linear problems. Int. J. Nonlinear Mech. 2000, 35, 37–43. [Google Scholar] [CrossRef]
  15. He, J.H. New interpretation of homotopy perturbation method. Int. J. Mod. Phys. B 2006, 20, 2561–2568. [Google Scholar] [CrossRef]
  16. Chang, S.H.; Chang, I.L. A new algorithm for calculating one-dimensional differential transform of nonlinear functions. Appl. Math. Comput. 2008, 195, 799–808. [Google Scholar] [CrossRef]
  17. Arikoglu, A.; Ozkol, I. Solution of fractional differential equations by using differential transform method. Chaos Solitons Fractals 2007, 34, 1473–1481. [Google Scholar] [CrossRef]
  18. Khuri, S.A. A Laplace decomposition algorithm applied to a class of nonlinear differential equations. J. Appl. Math. 2001, 1, 141–155. [Google Scholar] [CrossRef]
  19. Khuri, S.A. A new approach to Bratu’s problem. Appl. Math. Comput. 2004, 147, 131–136. [Google Scholar] [CrossRef]
  20. Khan, Y.; Wu, Q. Homotopy perturbation transform method for nonlinear equation using He’s polynomials. Comput. Math. Appl. 2011, 61, 1963–1967. [Google Scholar] [CrossRef]
  21. Khana, M.; Gondala, M.A.; Hussain, I.; Vanani, S.K. A new comparative study between homotopy analysis transform method and homotopy perturbation transform method on a semi-infinite domain. Math. Comput. Model. 2012, 55, 1143–1150. [Google Scholar] [CrossRef]
  22. Mohamed, M.S.; Al-Malki, F.; Al-humyani, M. Homotopy Analysis Transform Method for Time-Space Fractional Gas Dynamics Equation. Gen. Math. Notes 2014, 24, 1–16. [Google Scholar]
  23. Singh, J.; Kumar, R. An Analytical Study on Fractional Fokker-Planck Equation by Homotopy Analysis Transform Method. Recent Adv. Math. Stat. Comput. Sci. 2016, 313–321. [Google Scholar]
  24. Mohamed, M.S.; Gepreel, K.A.; Alharthi, M.R.; Alotabi, R.A. Homotopy Analysis Transform Method for Integro-Differential Equations. Gen. Math. Notes 2016, 32, 32–48. [Google Scholar]
  25. Saad, K.M.; AL-Shareef, E.H.F.; Alomari, A.K.; Baleanu, D.; Gómez-Aguilar, J.F. On exact solutions for time-fractional Korteweg-de Vries and Korteweg-de Vries-Burger’s equations using homotopy analysis transform method. Chin. J. Phys. 2020, 63, 149–162. [Google Scholar] [CrossRef]
  26. Nemah, E.M. Homotopy Transforms Analysis Method for Solving Fractional Navier Stokes Equations with Applications. Iraqi J. Sci. 2020, 61, 2048–2054. [Google Scholar] [CrossRef]
  27. Veeresha, P.; Prakasha1, D.G.; Baskonus, H.M. Solving smoking epidemic model of fractional order using a modified homotopy analysis transform method. Math. Sci. 2019, 13, 115–128. [Google Scholar] [CrossRef]
  28. Khanlari, N.; Paripour, M. Solving Nonlinear Integro-Differential Equations Using the Combined Homotopy Analysis Transform Method with Adomian Polynomials. Commun. Math. Appl. 2018, 9, 637–650. [Google Scholar]
  29. Eltayeb, H.; Mesloub, S.; Kiliçman, A. Application of double Laplace decomposition method to solve a singular one-dimensional pseudohyperbolic equation. Adv. Mech. Eng. 2017, 9, 1687814017716638. [Google Scholar] [CrossRef]
Figure 1. The -curve corresponding to the 13th-order truncated series solution at x = 0.1 .
Figure 1. The -curve corresponding to the 13th-order truncated series solution at x = 0.1 .
Axioms 13 00398 g001
Figure 2. Truncated series solution U [ i ] ( x , t ) with various values of i for Example 1.
Figure 2. Truncated series solution U [ i ] ( x , t ) with various values of i for Example 1.
Axioms 13 00398 g002
Figure 3. The -curve corresponding to the 10th-order truncated series solution at x = 0.1 .
Figure 3. The -curve corresponding to the 10th-order truncated series solution at x = 0.1 .
Axioms 13 00398 g003
Figure 4. Truncated series solution U [ i ] ( x , t ) with various values of i for Example 3.
Figure 4. Truncated series solution U [ i ] ( x , t ) with various values of i for Example 3.
Axioms 13 00398 g004
Table 1. Absolute error at x = 0.1 and different values of , i, and t.
Table 1. Absolute error at x = 0.1 and different values of , i, and t.
i t
0.10.512
Er [ i ] Er [ i ] Er [ i ] Er [ i ]
−0.91 0.0186845 0.515250 2.22646 9.02557
5 9.30234 × 10 6 2.57433 × 10 4 1.11289 × 10 3 4.51242 × 10 3
10 1.85947 × 10 10 5.14819 × 10 9 2.2257 × 10 8 9.02475 × 10 8
15 2.77642 × 10 15 7.6621 × 10 14 3.34507 × 10 13 1.35459 × 10 12
−11 0.0206497 0.571968 2.47291 10.0274
5 8.45678 × 10 18 1.96024 × 10 16 1.17267 × 10 15 1.12930 × 10 15
10 8.45678 × 10 18 1.96024 × 10 16 1.17267 × 10 15 1.12930 × 10 15
15 8.45678 × 10 18 1.96024 × 10 16 1.17267 × 10 15 1.12930 × 10 15
−1.11 0.0226148 0.628685 2.71936 11.0292
5 1.13473 × 10 5 3.14534 × 10 4 1.36001 × 10 3 5.51498 × 10 3
10 2.27046 × 10 10 6.29116 × 10 9 2.72011 × 10 8 1.10300 × 10 7
15 3.57505 × 10 15 9.37245 × 10 14 4.08033 × 10 13 1.65583 × 10 12
Table 2. Absolute error at x = 5 and different values of , i, and t.
Table 2. Absolute error at x = 5 and different values of , i, and t.
i t
0.10.512
Er [ i ] Er [ i ] Er [ i ] Er [ i ]
−0.91 0.268168 1.71333 4.32929 11.2979
5 3.42507 × 10 5 3.77242 × 10 4 1.32318 × 10 3 4.73965 × 10 3
10 4.35430 × 10 10 6.34627 × 10 9 2.43598 × 10 8 9.25198 × 10 8
15 5.32907 × 10 15 9.05942 × 10 14 3.51719 × 10 13 1.37490 × 10 12
−11 0.0206497 0.571968 2.47291 10.0274
5 4.44089 × 10 16 3.55271 × 10 15 3.55271 × 10 15 3.55271 × 10 15
10 4.44089 × 10 16 3.55271 × 10 15 3.55271 × 10 15 3.55271 × 10 15
15 4.44089 × 10 16 3.55271 × 10 15 3.55271 × 10 15 3.55271 × 10 15
−1.11 0.226869 0.569399 0.616520 8.75689
5 1.36010 × 10 5 1.94726 × 10 4 1.14973 × 10 3 5.28774 × 10 3
10 2.24376 × 10 11 5.09308 × 10 9 2.50983 × 10 8 1.08028 × 10 7
15 1.33227 × 10 15 8.17124 × 10 14 3.87246 × 10 13 1.63780 × 10 12
Table 3. Absolute error at x = 0.1 and different values of , i, and t.
Table 3. Absolute error at x = 0.1 and different values of , i, and t.
i t
0.10.512
Er [ i ] Er [ i ] Er [ i ] Er [ i ]
−0.91 0.010150 0.068750 0.250000 1.40000
5 1.07568 × 10 6 1.65992 × 10 5 1.67883 × 10 4 5.12458 × 10 3
10 1.1533 × 10 11 3.46751 × 10 10 6.79202 × 10 9 6.18511 × 10 7
15 4.35207 × 10 14 4.64073 × 10 13 1.23013 × 10 13 3.69571 × 10 11
−11 1.66667 × 10 4 2.08333 × 10 2 0.166667 1.33333
5 0.1 . 22324 × 10 11 5.38229 × 10 9 2.50521 × 10 8 5.13067 × 10 5
10 7.01661 × 10 14 8.05134 × 10 13 3.31735 × 10 13 5.12479 × 10 12
15 6.21725 × 10 15 2.22045 × 10 15 1.64313 × 10 14 1.79412 × 10 13
−1.11 9.81667 × 10 3 2.70833 × 10 2 8.33333 × 10 2 1.26667
5 9.09339 × 10 7 3.50969 × 10 6 5.26508 × 10 6 8.32607 × 10 5
10 8.21831 × 10 12 6.05289 × 10 11 3.86545 × 10 10 5.89687 × 10 9
15 2.14051 × 10 13 5.7776 × 10 13 4.64073 × 10 13 5.77316 × 10 13
Table 4. Absolute error at x = 5 and different values of , i, and t.
Table 4. Absolute error at x = 5 and different values of , i, and t.
i t
0.10.512
Er [ i ] Er [ i ] Er [ i ] Er [ i ]
−0.91 0.010150 0.068750 0.250000 1.40000
5 1.07568 × 10 6 1.65992 × 10 5 1.67883 × 10 4 5.12458 × 10 3
10 1.15328 × 10 11 3.46751 × 10 10 6.79202 × 10 9 6.18511 × 10 7
15 2.28706 × 10 14 3.8658 × 10 13 1.30784 × 10 13 3.6539 × 10 11
−11 1.66667 × 10 4 2.08333 × 10 2 0.166667 1.33333
5 2.22045 × 10 6 1.22322 × 10 11 2.50521 × 10 8 5.13067 × 10 5
10 3.39728 × 10 14 8.41327 × 10 13 2.95541 × 10 13 4.70668 × 10 12
15 6.88338 × 10 15 2.88658 × 10 15 1.70974 × 10 14 1.59206 × 10 13
−1.11 9.81667 × 10 3 2.70833 × 10 2 8.33333 × 10 2 1.26667
5 9.09339 × 10 7 3.50969 × 10 6 5.26508 × 10 6 8.32607 × 10 5
10 8.2172 × 10 12 6.05211 × 10 11 3.86545 × 10 10 5.89683 × 10 9
15 1.36557 × 10 13 1.59206 × 10 13 5.00266 × 10 13 1.0687 × 10 12
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mesloub, S.; Gadain, H.E. Homotopy Analysis Transform Method for a Singular Nonlinear Second-Order Hyperbolic Pseudo-Differential Equation. Axioms 2024, 13, 398. https://doi.org/10.3390/axioms13060398

AMA Style

Mesloub S, Gadain HE. Homotopy Analysis Transform Method for a Singular Nonlinear Second-Order Hyperbolic Pseudo-Differential Equation. Axioms. 2024; 13(6):398. https://doi.org/10.3390/axioms13060398

Chicago/Turabian Style

Mesloub, Said, and Hassan Eltayeb Gadain. 2024. "Homotopy Analysis Transform Method for a Singular Nonlinear Second-Order Hyperbolic Pseudo-Differential Equation" Axioms 13, no. 6: 398. https://doi.org/10.3390/axioms13060398

APA Style

Mesloub, S., & Gadain, H. E. (2024). Homotopy Analysis Transform Method for a Singular Nonlinear Second-Order Hyperbolic Pseudo-Differential Equation. Axioms, 13(6), 398. https://doi.org/10.3390/axioms13060398

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop