Traveling Wave Solutions to the Nonlinear Evolution Equation Using Expansion Method and Addendum to Kudryashov's Method

The inspection of wave motion and propagation of diffusion, convection, dispersion, and dissipation is a key research area in mathematics, physics, engineering, and real-time application fields. This article addresses the generalized dimensional Hirota–Maccari equation by using two different methods: the exp(−φ(ζ)) expansion method and Addendum to Kudryashov’s method to obtain the optical traveling wave solutions. By utilizing suitable transformations, the nonlinear PDEs are transformed into ODEs. The traveling wave solutions are expressed in terms of rational functions. For certain parameter values, the obtained optical solutions are described graphically with the aid of Maple 15 software.


Introduction
The study of the traveling wave solutions plays a significant role in understanding and describing the characters of nonlinear problems in physical science. The elements of physical systems usually operate on multiple time scales [1]. Although closed form descriptions for nonlinear partial differential equations (NLPDEs) of physical significance exist, we cannot obtain these forms explicitly. Such problems occur especially in various realistic problems in physical systems. In this scenario, we aim to investigate exact physically significant solutions, which are so important and useful due to their very wide application of NLPDEs in fluid mechanics, fluid dynamics, quantum mechanics, bio-science, physics, chemistry, and other areas of engineering.
Modeling of wave motion has attracted many researchers due to its crucial role in ocean, coastal, naval and marine engineering. In addition, waves are considered to be the major source of environmental actions on beaches or floating structures for most geographical areas. Most problems of practical interest create physical phenomena in nature that are often nonlinear and can be described by fractional differential equations. For example, heat conduction systems, nonlinear chaotic systems, plasma waves, and diffusion processes are modeled by fractional mathematics and equations [2][3][4]. Therefore, numerous investigations have been conducted to develop new methods to solve such equations. For example, Choucha et al. [5] investigated problems of a nonlinear viscoelastic wave equation in the presence of distributed delay, strong damping and source terms. By considering suitable conditions, they obtained a blow-up result of solutions. The improved sub-equation method has been used by Zhong et al. [6] to obtain exact solutions for a wide range of nonlinear problems of fractional equations. By converting nonlinear problems of fractional and integer order equations with the modified Riemann-Liouville derivative into Riccati equations, they obtained exact solutions, including wave solutions, soliton solutions, and complex solutions. Simbanefayi et al. [7] obtained traveling wave solutions for the Korteweg-de Vries-Benjamin-Bona-Mahony equation using Lie symmetry and Jacobi elliptic function expansion methods. The main objective of the Lie symmetry method is to transform the governing equations to a simpler equation while maintaining the invariance of the original equation. Others have previously used different methods to find the traveling wave solutions for nonlinear evolution equations (NEEs) (see, for example, [8][9][10][11][12][13][14][15][16][17][18][19][20]).
Over the past few years, many researchers have explored several direct methods to solve nonlinear evolution equations. One of the popular examples of such methods is the Kudryashov method, which typically only performs the calculation without the need for the form of a specific function. Examples of Kudryashov methods include the extended Kudryashov method [21], the generalized Kudryashov method [22], and the new extended generalized Kudryashov method [23]. The generalized Kudryashov method by Kaplan et al. [22] does usefully apply to the nonlinear Jaulent-Miodek hierarchy and (2 + 1)-dimensional Calogero-Bogoyavlenskii-Schiff equation. This method has successfully provided exact solutions to the nonlinear evolution equations. The Addendum to Kudryashov's method is also among the Kudryashov methods and was introduced by [24]. This method is the general form of previous Kudryashov methods because the trial equation is proposed as the general form of the trial equations in other Kudryashov methods.
For the purposes of discussion, we consider the nonlinear (2 + 1) dimensional Hirota-Maccari equation [27]. Let x and y be the independent spatial variables and t be the time variable. Consider the complex and the real scalar fields u = u(x, y, t) and v = v(x, y, t), respectively, satisfying the (2 + 1)-dimensional Hirota-Maccari equation Assuming that x = y , then the system (1) reduces to the (1 + 1)-dimensional Hirota equation [28].
In the last two decades, scientists have explored the Hirota-Maccari system using several approaches to efficiently compute and predict exact solutions. The approach by Painlevé has been applied successfully to construct the general solutions of certain nonlinear second-order equations [27]. Other methods that have been introduced include a new unified algebraic method [29], Weierstrass elliptic function expansion method [30] and extended trial method [31].
Here, we focus attention on using the expansion method and AKM to predict and derive exact solutions of the Hirota-Maccari system. This paper is organized as follows. Section 2 presents the methodology of the exp(−ϕ(ζ))-expansion method. Section 3 discusses the Hirota-Maccari equation. In addition, numerical computations are carried out by Maple to present the behavior of the optical traveling wave solution through figures. Section 4 addresses the (2 + 1)-dimensional Hirota-Maccari equation by using Addendum to Kudryashov's method (AKM) to find the straddled solitary wave solution and the singular soliton solution and presents them in graphs.

Methodology of the Expansion Method
This section presents the procedure of the exp(−ϕ(ζ))-expansion method. In the first step, we consider a nonlinear partial differential equation in the following form: where u = u(x, t) represents an unknown function, and F is a polynomial function in u(x, t) and its partial derivatives. We define an appropriate traveling wave transformation as follows, where c is the velocity of a traveling wave. By utilizing this appropriate traveling wave transformation (3), the nonlinear partial differential Equation (2) reduces to the ordinary differential equation (ODE): in which Q is a polynomial in u(ζ) and its total derivatives such that = d dζ . We seek traveling wave solutions for a large class of Equation (4) in the following form: where a i are arbitrary constants, N is a positive integer, and ϕ(ζ) is the elementary function that satisfies the following ODE: where λ and µ are arbitrary constants. By taking the transformation Ψ = exp(ϕ), the nonlinear ODE (6) can be written in the following form The above Equation (7) is the generalized Riccati first ODE's. By taking this transformation η = Ψ + λ/(2µ), the generalized Riccati first ODE's (7) reduces to the following form which can be solved by applying a separation of variables to obtain the following solutions: • First, when λ 2 − 4µ > 0, and µ = 0, we have the hyperbolic function solutions: or where E is the integration constant. • Second, when λ > 0, and µ = 0, we have solutions as: • Third, when λ 2 − 4µ < 0, we obtain the periodic solutions as: or • Fourth, when λ 2 − 4µ = 0, we obtain the solution as: • Fifth, when µ = λ = 0, we obtain the solution as: Substituting Equation (5) into Equation (4) and using ODE (6), the left-hand side is transferred to a polynomial in exp(−ϕ(ζ)). By equating the coefficient of this polynomial to zero, we obtain a system of algebraic equations that describes the model equations for a N , . . . , c, λ, and µ. These algebraic system of equations are solved by means of the Maple software. Substituting the values of the parameters a N , . . . , c, λ, µ, and the general solution of Equation (6) into Equation (5), we obtain the desired traveling wave solutions, as shown in Equation (2).

Application
The proposed method in the previous section is implemented successfully to the Hirota-Maccari system (1) to obtain explicit and exact traveling wave solutions. The traveling wave solutions of the Hirota-Maccari system (1) are expressed in the form [32] u where U and V are real functions, ζ = x + y + ct and c is the wave velocity, p is the frequency of the wave, while q and k are arbitrary constants. Using the traveling wave transformation (16) empowers us to transform the Hirota-Maccari system (1) into the following system of ODEs: where p = 1/3. By explicitly taking advantage of the homogeneous balance method [33], we balance the nonlinear term U 3 and the highest order derivative U in Equation (17) to obtain the balance constant N = 1. Therefore, the solution of Equation (17) takes the following form: where a 0 , a 1 are non-zero constants. Elementary algebra gives U and U 3 as: and Substitute Equations (18)- (20) into Equation (17) and collect all terms with the same power of exp(−iϕ), for (i = 0, 1, 2, 3). Then, setting the coefficients of exp(−iϕ), for (i = 0, 1, 2, 3) to zero, yields: Solving the above algebraic equations gives the values of the constants a 0 , a 1 , and µ as: Therefore, substituting these constants into the solution (18), we obtain To obtain the optical traveling wave solutions for the Hirota-Maccari system, we substitute Equations (9)-(15) into Equation (23) to obtain the following five optical traveling wave solutions: • Case 1: when λ 2 − 4µ = 4(p 3 −pq−k) 6(1−3p) > 0, and µ = 0, the complex and the real scalar fields can be expressed respectively as: and For special values of parameters, Figure 1 represents the real part (1a) and imaginary part (1b) of (24a) and its projections (1c) and (1d), respectively, while the exact solution (24b) and its projections are shown in Figure 2. Figure 1. Periodic wave solution of the real part (a) and imaginary part (b) of (24a) and its projections (c,d), respectively, when λ = 3, µ = 26/12, y = 2, and c = p = q = k = E = 1. • Case 2: when λ 2 = 4(p 3 −pq−k) 6(1−3p) > 0, µ = 0, and λ = 0 the complex and the real scalar fields can be expressed, respectively, as: For special values of parameters, Figure 3 represents the real part (3a) and imaginary part (3b) of (26a) and its projections (3c) and (3d), respectively, while the exact solution (26b) and its projections are shown in Figure 4. Figure 3. The singular solitary wave solution of the real part (a) and imaginary part (b) to (26a) and its projections (c,d), respectively, when λ = 1/ √ 3, µ = 0, y = 2, and c = p = q = k = E = 1. • Case 3: when λ 2 − 4µ = 4(p 3 −pq−k) 6(1−3p) < 0, and µ = 0, the complex and the real scalar fields can be expressed respectively as: Figure 5 represents the real part (5a) and imaginary part (5b) of (27a) and its projections (5c) and (5d), respectively, while the exact solution (27b) and its projections are shown in Figure 6, for special values of parameters. Figure 5. The periodic wave solution of the real part (a) and imaginary part (b) to (27b) and its projections (c,d), respectively, when λ = 3, µ = 7/3, y = 2, c = p = q = E = 1, and k = −1.  (28b) • Case 5: when λ 2 = 4(p 3 −pq−k) 6(1−3p) = 0, µ = 0, and λ = 0, the complex and the real scalar fields can be expressed, respectively, as: (29b) Figure 7 represents the singular Kink wave solution of the real part (7a) and imaginary part (7b) of (29a) and its projections (7c) and (7d), respectively, while the singular Kink wave solution (29b) and its projections are shown in Figure 8, for special values of parameters. Figure 7. The singular Kink solution of the real part (a) and imaginary part (b) to (29a) and its projections (c,d), respectively, when λ = µ = 0, y = 2, c = p = E = 1, and q = k = 1/2.

Hirota-Maccari Equation
For nonlinear evolution equations (NLEEs), Kudryashov [34] introduced a new approach to find highly dispersive optical solitons. This method is intended to be used to perform the calculation without the need for the form of a specific function. Inspired by the work of Kudryashov [34], Zayed et al. [24] recently introduced the Addendum to Kudryashov's method AKM. One aim of this section is to apply this method to find the straddled solitary wave solutions and the singular soliton solutions of the (2 + 1)-dimensional Hirota-Maccari equation.
For convenience, we start by presenting the main steps of the AKM as follows: • Step 1: Assuming that (17) has a solution in the following form where σ g = 0 for (g = 0, 1, 2, . . . , N) are constants can be determined, and R(ξ) satisfies the NODE so that where χ is an arbitrary constant. We can verify that Equation (31) can be written in the form: where A represents a non-zero real number, T is a natural number, and exp K (Tξ) = K Tξ . • Step  N). To evaluate the values of σ g for (g = 0, 1, 2, . . . , N) and c, we must solve this system of equations. Consequently, we will obtain the analytical solutions of Equation (17).

Conclusions
The search for exact solutions for NEEs has attracted the attention of many scientists in physics and mathematics. In this paper, I investigate the optical traveling wave solutions to a nonlinear evolution equation in mathematical physics, namely for the (2 + 1)-dimensional Hirota-Maccari equation by means of the expansion method. Additionally, by following the method of the Addendum to Kudryashov's method AKM, introduced by Zayed et al. [24], Section 4 analogously derives the straddled solitary wave solution and the singular soliton solution of the (2 + 1)-dimensional Hirota-Maccari Equation (1), which is presented with the support of graphs. The resulting solutions indicate that the proposed approach promises to empower us to address a wide class of nonlinear evolution equations arising in mathematical physics. The optical wave solutions obtained in this paper is just one important example of solutions that can have the potential to empower systematic analysis and understanding of the nonlinear evolution equation in many fields used widely in engineering and science.