ORIGINAL RESEARCH article
Sec. Statistical and Computational Physics
Volume 8 - 2022 | https://doi.org/10.3389/fams.2022.952668
Wave solutions of the DMBBM equation and the cKG equation using the simple equation method
- 1Department of Mathematics, Faculty of Science and Technology, Suratthani Rajabhat University, Suratthani, Thailand
- 2Education Program in Mathematics, Faculty of Education, Suratthani Rajabhat University, Suratthani, Thailand
In this article, we transform the (1 + 1)-dimensional non-linear dispersive modified Benjamin-Bona-Mahony (DMBBM) equation and the (2 + 1)-dimensional cubic Klein Gordon (cKG) equation, which are the non-linear partial differential equations, into the non-linear ordinary differential equations by using the traveling wave transformation and solve these solutions with the simple equation method (SEM) with the Bernoulli equation. Two classes of exact explicit solutions-hyperbolic and trigonometric solutions of the associated NLEEs are characterized with some free parameters; we obtain the kink waves and periodic waves.
Immediately, nonlinear evolution equations (NLEEs), i.e., partial differential equations with time derivatives and fractional differential equations (FDEs), have become a useful tool for describing the natural phenomena of science and engineering, such as optical fiber communications, atmospheric pollutant dispersion, solid state physics, signal processing, mechanical engineering, electric control theory, relativity, chemical reactions, etc. Exact solutions or numerical solutions have always played and still play an important role in properly understanding the qualitative features of many phenomena and processes in various fields of natural science. Various effective and powerful methods have been established to handle the NLEEs, such as the modified simple equation method [1, 2], the Jacobi elliptic function method , the (G′/G)-expansion method , the homotopy perturbation method , the variational iteration method , the Exp-function method [7, 8], the tanh function method , the F-expansion method , the Laplace-optimized decomposition method , the reproducing kernel algorithm [12–14] etc. The investigation of wave solutions of NLEEs plays a significant role in the study of nonlinear physical phenomena. The well-known wave solutions that are utilized to describe the exact solutions of NLEEs are lump solutions, lump-multi-kink solutions, and traveling wave solutions. Lump solutions and lump-multi-kink solutions to the (3 + 1) dimensional generalized Boiti–Leon–Manna–Pempinelli equation are studied in 2022 . Lump-soliton solutions to the KPI equation are shown in 2021 . Traveling wave solutions for the (2 + 1)-dimensional cubic nonlinear Klein–Gordon (cKG) equation and the (2 + 1)-dimensional nonlinear Zakharov–Kuznetsov modified equal width (ZK-MEW) equation were being investigated in 2019 . Additionally, the (1 + 1)-dimensional nonlinear dispersive modified Benjamin-Bona Mahony (DMBBM) equation and the (2 + 1)-dimensional cubic Klein Gordon (cKG) equation are non-linear partial differential equations that play a significant role in various scientific and engineering fields. The DMBBM equation was first derived to describe an approximation for surface-long waves in non-linear dispersive media. It can also characterize the hydromagnetic waves in cold plasma, acoustic waves in inharmonic crystals, and acoustic gravity. The DMBBM was studied in 2010 using the extended (G′/G)-expansion method , and, in 2014, using the modified simple equation method . The Klein-Gordon equation has been used to model a wide range of nonlinear phenomena, including the propagation of discrepancies in crystals and elementary particle behavior. This equation was studied in 2011  using the extended (G′/G)-expansion method, in 2013  using the modified simple equation method, and in 2019  using the Riccati Bernoulli sub ODE method.
The simple equation method (SEM) was presented by Nikolai Kudryashov . This method was created by two important ideas. One of them is to use the general solutions of the simplest non-linear differential equations. Another idea is to take into consideration all possible singularities in the equation studied. The simple equation method was used to investigate exact solutions of various equations, such as the Sharma Tasso Olver equation and the Burgers Huxley equation in 2008 , the Kodomtsev Petviashvili equation, the (2 + 1) dimensional breaking soliton equation, and the modified generalized Vakhnenko equation in 2016  and the fourth-order non-linear AKNS water equation in 2021 .
In this article, we use the traveling wave to transform the DMBBM equation and the cKG equation into the nonlinear ordinary differential equations. Then, using the simple equation method with the Bernoulli equation, we obtain the precise solution of these equations in terms of the exponential functions, and the physical wave solution is created in the form of kink waves and periodic waves. Moreover, we compare these results presented with other results in order to show that the simple equation method with the Bernoulli equation is more convenient and easier to understand.
Algorithm of simple equation method
In this section, we present a direct method, namely, the simple equation method, for finding the traveling wave solution to nonlinear equations. Suppose that the nonlinear partial equation, in two independent variablesx and t, is given by:
where u(x, t) is an unknown function, P is a polynomial of u(x, t) and its partial derivatives in which the highest-order derivatives and non–linear terms are involved. The main steps of the SE method  are as follows:
Step 1: Wave transformation
Combining the independent variables x and t into one variable ξ = x−ωt, we suppose that
where ω is the speed of a traveling wave. The traveling wave transformation Equation (2) permits us to convert Equation (1) into an ordinary differential equation (ODE) for u = u(ξ):
where Q is a polynomial of u(ξ) and its derivatives in which the prime indicates the derivative with respect to ξ.
Step 2: Solution assumption
Assume that the formal solution of Equation (3) is of the form:
where ai(i = 0, 1, 2, ..., n) are constants that will be determined later and F(ξ) are the functions that satisfy the simple equations (ordinary differential equations). The simple equation has two characteristics. First of all, it lacks the order of Equation (2); secondly, we are aware of the general solution to the simple equation. In this article, we shall use the Bernoulli and Riccati equations, which are well-known nonlinear ordinary differential equations, and their solutions can be described by simple functions. Regarding the Bernoulli equation,
Step 3: Finding the integer n
The positive integer n that occurs in Equation (4) can be estimated by taking into account the homogeneous balance between the highest-order derivative and the non-linear terms appearing in Equation (3).
Step 4: Solution attainment
We get the general solutions of the simple Equation (5) as follows:
Case I: if c > 0 and d < 0, we have
where ξ0 is a constant of the integration.
Case II: if c < 0 and d > 0, we have
where ξ0 is a constant of the integration.
In this section, we apply the proposed simple equation method to construct the traveling wave solution to solve the (1 + 1)-dimensional non-linear dispersive modified Benjamin-Bona-Mahony (DMBBM) and the (2 + 1)-dimensional cubic Klein Gordon equation.
Solution of the (1 + 1)-dimensional non-linear dispersive modified Benjamin-Bona- Mahony
The DMBBM equation is
where α is a non-zero positive constant. Using the traveling wave variable ξ = x − ωt is transformed (8) into the following ODE:
According to the balancing procedure that will be described, the balancing number n is a positive integer, which can be defined by balancing the highest-order derivative terms with the highest power non-linear terms in Equation (9), i.e., n+3 = 3n+1, thus n = 1. We have the solution of Equation (9) as follows:
where F satisfies Equation (5); consequently, u′, u‴ and u2 can be expressed as follows:
Substituting Equations (10) and (11) into Equation (9) then equating the coefficient of Fi to zero, where i ≥ 0, yields
Solving this system of algebraic equations, we obtain
Substituting Equations (13) and (14) into (10) yields the following two exact solutions. Then, we use the general solutions of Bernoulli equations (6) and (7). We get four exact solution (8) from the exponential term. Next, we set the constants c and d to obtain the exact solutions in hyperbolic form as follows:
Case I: if c = 1 and d = −1, we have
Case II: if c = −1 and d = 1, we have
where k is arbitrary constants, it might yield many new and more general exact solutions to the non-linear DMBBM equation.
Solution of the (2 + 1)-dimensional cubic Klein Gordon (cKG) equation
The (2 + 1)-dimensional cubic Klein Gordon (cKG) equation is
where α and β are non-zero positive constants. Using the traveling wave variable ξ = x + y − ωt is transformed (17) into the following ODE:
According to the balancing procedure that will be described, the balancing number n is a positive integer, which can be defined by balancing the highest-order derivative terms v″ with the highest power non-linear terms v3 in Equation (18), i.e., n + 2 = 3n; thus, n = 1. We have the solution of Equation (18) as follows:
where F satisfies Equation (5); consequently, and v3 can be expressed as follows:
Substituting Equations (19) and (20) into Equation (18) and then equating the coefficient of Fi to zero, here i ≥ 0, yields
Solving this system of algebraic equations, we obtain
Substituting Equations (22) and (23) into (19) yields the following two exact solutions. Then, we use the general solutions of Bernoulli equations (6) and (7). We get four exact solutions of (17) from the exponential term. Next, we set the constants c and d to obtain the exact solutions in hyperbolic and trigonometric form as follows:
Case I: if d = −1 and c > 0, we have
where and k is arbitrary constants.
Case II: if and c < 0, we have
where and k is arbitrary constants. It might yield many new and more general exact solutions to the non-linear cKG equation.
Graphical representation of some obtained solutions
In this section, we discuss the physical explanations and graphical representations of the solutions of the DMBBM equation and the cKG equation.
Graphical representation of the solution of the DMBBM equation
Upon utilizing the SE method, we achieve the traveling wave solutions of the DMBBM equation from Equations (15) and (16) when wave speed . The solution u1(x, t) represents by a singular kink-type wave for the parameter α = 6, k = 2 in the interval −10 ≤ x, t ≤ 10 and is shown in Figure 1 The kink wave rises or descends from one asymptotic state to another. The kink solution approaches a constant at infinity. The solutions u2(x, t), u3(x, t) and u4(x, t) are analogous to the figures of solution u1(x, t) and are presented in Figures 2–4, respectively.
Graphical representation of the cKG equation
In this subsection, we represent the shape of the solutions of the cKG equation in the form of Equations (24) and (25). The solution v1(x, t) is formed by a singular kink-type wave when the wave speed ω = 2 and the parametersβ = 1, k = 2and α = −1 in the interval −10 ≤ x, t ≤ 10 and is demonstrated in Figure 5. This wave rises from one state to another state when t decreases. The solution v2(x, t) also represents the same shape with v1(x, t), but this wave raises from one state to another state when increased and is shown in Figure 6.
The solutions v3(x, t) and v4(x, t) are a formed periodic traveling wave for the values β = 1, k = 2 and α = 1in the interval −10 ≤ x, t ≤ 10 and is presented in Figures 7, 8, respectively. A periodic traveling wave is a periodic function that moves with constant speed. Consequently, it is a special type of spatiotemporal oscillation that is a periodic function of both space and time.
In this section, we compare our solutions with some well-known methods.
Comparisons of the solution of the DMBBM equation
In 2014, Khan et al.  investigated the exact solutions of the DMBBM equation by using the MSE method, and they got
where c1 = −2c2(1−ω) are similar with our solutions, that is, c1 = −1 and c2 = 1, then u1, 2(x, t) is identical with u1(x, t) and u2(x, t) where k = 0 and u5, 6(x, t) are equal with u3(x, t) and u4(x, t) where k = 0.
Comparisons of the solution of the cKG equation
In 2019, Abdelrahman et al.  use the Riccati-Bernoulli sub-ODE method to investigate the exact solutions of the cKG equation.
where is an identical with our solutions v3(x, t) and v4(x, t).
where is an equal with our solutions v1(x, t) and v2(x, t).
The comparison of our exact solutions by using the SE method presented that the method is more convenient and the solutions have an easier format.
The simple equation method presented in this article has been successfully implemented to find the wave solutions for the DMBBM equation and the cKG equation. The method offers solutions with free parameters in the form of exponential, hyperbolic, and trigonometric functions that might be important to explain some intricate physical phenomena. Some special solutions, including the known traveling wave solution, are originated by setting appropriate values for the parameters. The wave effects of this equation were kink waves as shown in Figures 1–6 and periodic waves as displayed in Figures 7, 8. Compared to the currently proposed method with other methods, such as the MSE method and the Riccati-Bernoulli sub-ODE method, we assert that the SE method is a potent mathematical methodology, very simple, and has easy computations.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
5. Roozi A, Alibeiki E, Hosseini SS, Shafiof SM, Ebrahimi M. Homotopy perturbation method for special nonlinear partial differential equations. J King Saud Univ (Sci). (2011) 23:1. doi: 10.1016/j.jksus.2010.06.014
8. Naher H, Abdullah AF, Akbar MA. New traveling wave solutions of the higher dimensional nonlinear partial differential equation by the Exp-function method. J Appl Math. (2012). doi: 10.1155/2012/575387
11. Beghami W, Maayah B, Bushnaq S, Abu Arqub O. The Laplace optimized decomposition method for solving systems of partial differential equations of fractional order. Int J App Computational Math. (2022) 8:52. doi: 10.1007/s40819-022-01256-x
12. Djennadi S, Shawagfeh N, Abu Arqub O. A numerical algorithm in reproducing kernel-based approach for solving the inverse source problem of the time-space fractional diffusion equation. Partial Diff Equa Applied Math. (2021) 4. doi: 10.1016/j.padiff.2021.100164
13. Momani S, Abu Arqub O, Maayah B. Piecewise optimal fractional reproducing kernel solution and convergence analysis for the Atangana–Baleanu–Caputo model of the Lienard's equation. Fractals. (2020)28:8. doi: 10.1142/S0218348X20400071
14. Momani S, Maayah B, Abu Arqub O. The reproducing kernel algorithm for numerical solution of Van der Pol damping model in view of the Atangana–Baleanu fractional approach. Fractals. (2020) 28:08. doi: 10.1142/S0218348X20400101
16. Lü X, Chen SJ. New general interaction solutions to the KPI equation via an optional decoupling condition approach. Commun Non-Linear Sci Num Simulation. (2021) 103. doi: 10.1016/j.cnsns.2021.105939
17. Abdelrahman MAE, Sohaly MA, Alharbi A. The new exact solutions for the deterministic and stochastic (2+1)-dimensional equations in natural sciences. J Taibah Univ Sci. (2019) 13:1. doi: 10.1080/16583655.2019.1644832
18. Zayed EME, Joudi SA. Applications of an Extended (G'/G)-Expansion Method to Find Exact Solutions of Nonlinear PDEs in Mathematical Physics Mathematical Problems in Engineering. (2010). doi: 10.1155/2010/768573
19. Khan K, Akbar MA, Islam SMR. Exact solutions for (1+1)-dimensional nonlinear dispersive modified Benjamin-Bona-Mahony equation and coupled Klein-Gordon equations. Springerplus. (2014) 3:1. doi: 10.1186/2193-1801-3-724
21. Khan K, Akbar MA. Exact solutions of the (2+1)-dimensional cubic Klein–Gordon equation and the (3+1) dimensional Zakharov Kuznetsov equation using the modified simple equation method. J Assoc Arab Univ Basic Appl Sci. (2013) 15. doi: 10.1016/j.jaubas.2013.05.001
Keywords: traveling wave solution, simple equation method, non-linear evolution equations, dispersive modified Benjamin-Bona-Mahony equation, cubic Klein Gordon equation
Citation: Sanjun J and Chankaew A (2022) Wave solutions of the DMBBM equation and the cKG equation using the simple equation method. Front. Appl. Math. Stat. 8:952668. doi: 10.3389/fams.2022.952668
Received: 20 June 2022; Accepted: 22 July 2022;
Published: 23 August 2022.
Edited by:Jordan Yankov Hristov, University of Chemical Technology and Metallurgy, Bulgaria
Reviewed by:Omar Abu Arqub, Al-Balqa Applied University, Jordan
Xing Lu, Beijing Jiaotong University, China
Copyright © 2022 Sanjun and Chankaew. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Aungkanaporn Chankaew, firstname.lastname@example.org