On Reduce Differential Transformation Method for Solving Damped Kawahara Equation

e damped Kawahara equation (KE) is nonintegrable equation and does not have analytical integration. In this work, the powerful numerical method, which is the reduce dierential transformation method (RDTM), is devoted to solve the damped KE. e accuracy of the method is proved. e results are compared with the dierent numerical methods. e numerical solution is axi-symmetric wave and shows the eect of damping term successfully. We conrmed that the RDTM is useful for solving nonintegrable equations.


Introduction
e partial deferential equations (PDEs) describe several important applications in many branches of science such as physics, engineering, medicine, and uid dynamic [1][2][3][4]. Mathematicians put forth high e orts to develop methods that are able to nd solutions of these PDEs [5][6][7][8]. Usually, as the PDEs describe a problem very well with taking all issues in account, there are some terms appear and make the PDEs are not solvable. erefore, the mathematicians improved the computational methods to nd di erent types of solutions such as exact, approximate, equivalent, numerical, and analytical.
One of the well-known PDEs is Korteweg-de Vries (KdV) equation and its family. e fth order of KDV is also known as Kawahara equation (KE). T. Ono and K. Ono [9] were the rst to discover this type of equation during the study of magneto-acoustic waves in a cool collision-free plasma. Kawahara numerically investigated this type of equation and discovered that it has both oscillatory and monotone solitary wave solutions [10]. In a uid medium like shallow water, the equation describes the propagation of soliton waves. e KE is governed by the following equation [11]: where α, β, and c are constants. e KE has been solved analytically and numerically in many researches [12][13][14]. e obtained solutions are N-soliton solutions [15], various solitons solutions [16], soliton and breathers [17], and different types of N-soliton and lump solutions [18]. e numerical solutions are obtained by using modi ed variational iteration algorithm-I and II [19,20], di erential quadrature [21], hybridizable discontinuous Galerkin (HDG) [22], and others. However, if a collisional e ect is taken into account in applications of KE equation, we obtain the damping term, and KE becomes damped KE with the following form: where C m/2 and m is the frequency of the ion-neutral collision. e damping term makes the (2) nonintegrable equation. In order to obtain the solutions, we aim to use a new improved technique. e di erential transformation method (DTM) is based on Taylor series expansion but di ers from the typical highorder Taylor series method, which takes a long time to calculate [23]. e DTM is one of the most powerful numerical methods. Pukhov was the rst who used the DTM to tackle linear and nonlinear initial value problems in electric circuit analysis [24]. Chen and Ho developed the DTM for solving PDEs and found closed form series solutions for a variety of linear and nonlinear initial value problems [25]. Abdel-Halim Hassan demonstrated that the DTM can be used on a wide range of PDES and easily obtain closed form solutions [26][27][28].
If the series of the solution has a closed form, then the numerical solution can be convergent to the exact solution, but this is not usually the case, especially in most realistic cases. us, the obtained solution is in series form. Since it is based on Taylor series, which is the local convergent [29], the DTM finds the solutions in small domain and about the initial point. It has been improved recently to reduce differential transformation method (RDTM) [30]. Keskin was the first who proposed the RDTM for finding exact solutions to PDEs [31,32]. Keskin and Oturanc created RDTM in recent years, in which the differential transformation is applied solely to one domain (time domain) [31]. e RDTM is a very effective and powerful tool for solving exact or approximate mathematical modeling solutions for a wide range of problems in technology, economics, engineering disciplines, and natural sciences such as biology, physics, chemistry, and earth science. It can solve both linear and nonlinear problems and provides results in the form of quick convergent successive approximations. e solutions by RDTM can also be classified as semiapproximate solution since the method applies the iteration only for the time domain. is technique is powerful compared to DTM and other methods. e novelty of this paper is proving that the RDTM is able to solve the class of nonintegrable equations, which does not have exact solutions. Such equations appear usually in physics applications when viscosity and ion-collisions are taken into account. We chose damped KE as an example of nonintegrable equations and devoted the RDTM to investigate the solution in long domain. e following is how the article is structured: Section 2 describes the used methods briefly, Section 3 presents the numerical solutions for KE and damped KE by RDTM, and Section 4 includes the conclusion of the work.

The Methodology
e DTM and its improved version (RDTM) are based on the following list of definitions. Definition 1 (differential transformation in two dimensions). e basic concept of the two-dimensional differential transform is as follows: let y(x, t) be analytic and continuously differentiable with respect to t and x, is the spectrum function [33]. e original function (lower case) y(x, t) is represented in this paper, whereas the converted function (upper case) Y(k, h) is represented. Using the twodimensional differential transformation (3), we present the differential transformation for several operators in Table 1.
Definition 2 (inverse differential transformation in two dimensions). e inverse differential transform of Y(k, h) is defined as follows [33]: Taking (3) and (4) together and assuming x 0 � t 0 � 0 yields to Definition 3 (reduce differential transformation and its inverse in two dimensions). If a(x, t) is analytical function in the domain of interest, then the spectrum function is used where a(x, t) is reduced transformed function. Lowercase a(x, t) refers to the original function, whereas uppercase A k (x) refers to the reduced transformed function. e differential inverse transformation of A k (x) is defined as [30] a( Combining (6) and (7) gives Table 2 shows the list of reduce differential transformation for several operators.

Kawahara Equation.
e first application is applying the DTM and RDTM into KE (1) in order to prove the accuracy of RDTM. In addition, we aim to prove the power of RDTM comparing to other methods in literature. Let's consider KE (1) with α � β � c � 1 and subjects to the initial condition [35].
e exact solution of this equation is given by where g � 1/2√13 and f � 36/169. We get the following scheme by using DTM in Definition 1 for k, h � 0, 1, 2, . . . , N, where N is the number of iterations: e initial condition is transformed into the following: e recursive equations deduced from (11) for different values of k, h are obtained as [36] k � 0, h � 0: We have noticed in Figure 1 that the numerical solution converges to exact solution in small interval about (− 4, 4) and diverges after that. Because of this disadvantage of DTM, the scheme is improved to RDTM as follows: e errors between the solution by RDTM and exact solution in defferent time are shown in Table 3. e solutions by DTM and RDTM are compared with the numerical solutions by optimal homotopy asymptotic method (OHAM) [35], homotopy perturbation and variational iteration method (VHPM) [37], homotopy perturbation method (HPM) [38], and Laplace homotopy perturbations method (LHPM) [39] in Table 4. e comparison reveals the accuracy of these methods. From Table 1: e fundamental operations by the two-dimensional differential transform method [33].  [30,34]. Table 4, we realized that the accuracy of RDTM and LHAM is better than that of the other methods, but RDTM is faster than LHPM. e speed of RDTM is 5.65 seconds, while for LHPM is 15.97 seconds for 6 iterations. erefore, RDTM is the optimal iteration method. Figure 2 shows the plot of the numerical solution of KE with IC [9]. e second example is, KE (1), where α � 3, β � 0.2, c � 0.4 and subjects to the IC [13].

Damped Kawahara Equation.
Because there is damping term in the Kawahara equation, the energy of the soliton is not conserved and decays with increasing both c and t, (2) is nonintegrable Hamiltonian system. We consider damped Kawahara (2) with α � 3, β � 0.2, c � 0.4 and subject to the IC [13]. Since we do not have exact solution, we can use the initial condition of Kawahara equation as initial condition of the damped Kawahara [13]. e scheme of the damped KE by RDTM is as follows: e numerical solution is shown in Figure 4. e amplitude of the wave decrease as the damping parameter increases.

Discussion and Conclusion
is paper studies the KdV-fifth order (Kawahara equation) within two cases: integrable KE and nonintegrable KE. e integrable KE has been solved in literature via different methods such as OHAM, VHPM, HPM, and LHAM. In this article, it is solved by DTM and RDTM to prove that RDTM converges to the solution faster than other methods with high accuracy. e new contribution in this work is solving nonintegrable KE, which includes damping term by RDTM. e two-dimensional DTM obtains the solutions in series form, but it is different from the traditional high-order Taylors series method, because it does not need symbolic computation of derivative for each term. Also, it does not require linearization, discretization, or other complected computation process. erefore, the DTM is faster than the Taylors series method. e DTM has been developed for solving ordinary and partial differential either linear or nonlinear equations. e improved version of the DTM is theRDTM, which is powerful to find numerical solutions for integrable equations as well as nonintegrable equations in several branches of science. MATLAB has been used for computations in this article. In future work, the RDTM can be applied to solve different new systems in physics and engineering that generate nonintegrable equations.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e author declares that there are no conflicts of interest.