Inelastic Interaction and Blowup New Solutions of Nonlinear and Dispersive Long Gravity Waves

School of Mathematics, Qilu Normal University, Jinan, 250200 Shandong Province, China School of Control Science and Engineering, Shandong University, Jinan, 250061 Shandong, China Department of Mathematics, Faculty of Science, Jiangsu University, China Department of Mathematics, El Obour Institutes, Cairo, Egypt Department of Basic Science, Higher Technological Institute 10th of Ramadan City, Egypt


Introduction
Studying the energetic atomic projectiles (atoms or ions) is one of the most exciting recent fields which has impacts on surface physics, plasma, and some related applications [1][2][3]. 13;3]. Moreover, it has an essential role which multicharged ions play as part of the solar wind in the space environment [4]. These applications include plasma-wall interaction in controlled thermonuclear fusion devices [5], surface analysis [6], single-particle detection [7], and electrical discharges [8]. The projectile energy deposition is connected with the intensity and nature of ion (atom) surface. Thus, both the nature and intensity of the ion depend on the potential and kinetic energies carried by a projectile toward the surface [9]. The energetic ion that can penetrate a solid surface is able to transfer its kinetic energy to target electrons, leading to atomic excitation, collective electron excitations, or ionization [10].
Partial differential equations (PDEs) have been playing an essential role in the energetic atomic projectiles where many nonlinear evolution equations have been derived to describe the dynamical behaviour of several phenomena in atomic and nuclear physics. Thus, partial differential equations (PDEs) have been playing an important role in the emerging technologies where many nonlinear evolution equations have been derived to describe the dynamical behaviour of some distinct phenomena, for example, nonlinear optics, fluid dynamics, Bose-Einstein condensates, and quantum mechanicss. However, the inadequacy of the PDEs with an integer order has been clarified because of the nonlocal property where this kind of equation does not explain that kind of properties. Therefore, several nature phenomena have been formulated with nonlinear PDEs with fractional order . Thus, many fractional operators have been derived such as conformable fractional derivative, fractional Riemann-Liouville derivatives, Caputo, Caputo-Fabrizio definition, and so on [36][37][38][39][40].
These definitions have been employed to convert the fractional nonlinear partial differential equations to a nonlinear integer-order ordinary differential equation. Then, the computational and numerical schemes can be applied to get various types of solutions for these models and the examples of these schemes.
The ABR fractional operator is considered one of the most general recent fractional operators that is derived from avoiding the deficiencies and defects of some other fractional operator. This operator is defined as follows.
where ½ð0 < α < 1Þ, S = Sðx, tÞ, W = W ðx, tÞ is unknown functions in time (t) and place (x). System (3) is used to simulate the bidirectional propagation of long waves in shallow water. Moreover, it studies the interaction between dispersive and nonlinear long gravity waves. System (3) is an integrable system and is related to the classical Boussinesq equation that is given by [50] where ½λ, α, β are arbitrary constants. Equation (5) is a model of nonlinear dispersive waves. Employing the next wave transformation (the ABR fractional operator) [Y = ððð1 − αÞðc t −n α Þ/BðαÞ∑ ∞ n=0 ð−ðα/ð1 − αÞÞÞ n Γð1 − αnÞÞ + xÞ, S = SðYÞ, W = W ðYÞ, where ðcÞ is an arbitrary constant] to the system ((3)); then, differentiate the results once and substitute the second equation of the system into the first equation leading to convert the system into the next equation Applying the homogeneous balance principle to Equation (6), yields ðm = 1Þ.
The rest of research paper is organized as follows: Section 2 applies the modified Khater method and sech-tanh functions expansion method to the suggested model to get novel solitary wave solutions of it. Section 3 explains the physical interpretation of the shown sketches in our paper. Section 4 shows the novelty of our paper by comparing our results with those obtained in previous research papers. Section 5 explains the conclusion of all the steps of our paper in detail.

Application
Here, in this section, the modified Khater method and sechtanh functions expansion method are applied to the BK equation to simulate the bi-directional propagation of long waves in sha llow water.

The Modified Khater Method.
Applying the modified Khater method to Equation (6) leads to formulate the general solution of this model in the following formula where ½a 0 , a 1 , b 1 are arbitrary constants to be determined later. Additionally, FðYÞ is the solution function of the following ordinary differential equation where ½δ, ϱ, χ are arbitrary constants. Substituting Equation (7) along (7) into Equation (6) and collecting all terms with the same power of ½K iFðYÞ , i = −5,−4, ⋯, 4, 5, give a system of the algebraic equation. Using the Mathematica 12 program for solving this system, yields Thus, the explicit wave solutions of Equation (4) are formulated in the following formulas 2 Journal of Function Spaces For ½χ 2 − 4δϱ > 0&δ ≠ 0, For ½δϱ < 0&ϱ ≠ 0&δ ≠ 0&χ = 0, For ½χ = δ = κ&ϱ = 0, For ½ϱ = 0&χ ≠ 0&δ ≠ 0 Family 3.

The Sech-Tanh Function Expansion Method.
Applying the sech-tanh function expansion method to Equation (6) leads to formulate the general solution of this model in the following formula: where ½a 0 , a 1 , b 1 are arbitrary constants to be determined later. Substituting Equation (23) into Equation (6) and collecting all terms with the same power of ½sech ðYÞ, tanh ðYÞ , give a system of the algebraic equation. Using the Mathematica 12 program for solving this system, yields Thus, the explicit wave solutions of Equation (4) are formulated in the following formulas: Family 5.
Journal of Function Spaces Thus, the explicit wave solutions of Equation (4) are formulated in the following formulas:

Figure Interpretation
This section gives the physical interpretation of the shown figures in our paper. All our obtained solutions are consid-ered traveling wave solutions. In the following lines, we give the physical interpretation of the shown figures: (1) Figure 1 shows the bell wave solution (8) in the twodimensional plot (a) to explain the wave propagation pattern of the wave along the x-axis and the contour plot (b) to explain the overhead view of the solution when ½α = ð1/2Þ, δ = 1, χ = 3, ϱ = 2: (2) Figure 2 shows the singular wave solution (13) in the two-dimensional plot (a) to explain the wave propagation pattern of the wave along the x-axis and the  (3) Figure 3 shows the W-wave solution (14) in the twodimensional plot (a) to explain the wave propagation pattern of the wave along the x-axis and the contour plot (b) to explain the overhead view of the solution when ½α = ð1/2Þ, δ = 1, χ = 5, ϱ = 6: (4) Figure 4 shows the rouge-wave solution (15) in the two-dimensional plot (a) to explain the wave propagation pattern of the wave along the x-axis and the contour plot (b) to explain the overhead view of the solution when ½α = ð1/2Þ, δ = 1, χ = 5, ϱ = 6: (5) Figure 5 shows the Breath-wave solution (21) in the two-dimensional plot (a) to explain the wave propagation pattern of the wave along the x-axis and the contour plot (b) to explain the overhead view of the solution when ½α = ð1/2Þ: (6) Figure 6 shows the kink-wave solution (22) in the two-dimensional plot (a) to explain the wave propagation pattern of the wave along the x-axis and the

Results and Discussion
This section shows the novelty of this research paper by explaining the comparison between our obtained solutions and those obtained in previous paper.

Computational
Schemes. This paper investigated the analytical solutions of the fractional BK system by using two recent computational schemes (modified Khater method and sech-tanh function expansion method). These methods are considered recent analytical schemes in this field, and they were not applied to this model yet. (ii) All our obtained solutions are different from those obtained in [47][48][49] where the authors of [47][48][49] applied different methods to solve the BK system with an integer order. On the other hand, we investigated the fractional form of this model

Conclusion
In this paper, we investigated new soliton wave solutions of the nonlinear fractional BK system by using two recent analytical schemes. A new fractional operator is used to convert the fractional form of this model to a nonlinear partial differential equation with an integer order. The modified Khater method and sech-tanh function expansion method were applied to this model. Some new soliton wave solutions were obtained, and some of them were explained by plotting them in two-dimensional and contour plots of these solutions. The novelty of our paper was shown by making the comparison between our obtained solutions and that was purchased in previously published articles.

Data Availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.

Ethical Approval
The authors state that this research complies with ethical standards. This research does not involve either human participants or animals.