Higher Order Solitary Wave Solutions of the Standard KdV Equations ()
1. Introduction
Many natural phenomena regularly occur in the universe. These phenomena can in certain cases be destructive for the environment in which they take place. For example, the cases of the propagation of the nonlinear excitations in fluid dynamics which sometimes, are manifested by tsunamis, hurricanes, tidal bore, and so on, and which very often are devastating when surge is across on the continent without warning. These different phenomena as presented in this case (among so many others) constitute real and permanent threats to our existence on the earth planet. It is following this observation that, a great number of researchers in particular, physicists and mathematicians, deploy themselves in everyday life in order to design and make available mathematical models [1] [2] [3] [4] [5] in an attempt to analyze, understand and explain these phenomena that can be observed in various fields such as nonlinear optics, mechanics, biology, and so on. Most of these models are nonlinear partial differential equations (NPDEs). From all these equations, we are interested, in this manuscript, in the one whose history throughout the last two centuries is intimately linked from the outset to the remarkable scientific discovery made in August 1834 [6] by John Scott Russell and who shows up under the form
(1)
Equation (1) is the standard form of the well-known fifth-order KdV (fKdV) equations [7] - [13] where
are real and nonzero arbitrary parameters,
is a sufficiently differentiable function, and, x, t are both independent spatial and temporal coordinates respectively.
and
are the two dispersive terms involved in this equation. Equation (1) is susceptible to a radical change in its characteristics due to the arbitrary and varied characters of the values that the parameters
can take. This changing and arbitrary character of the values that these different parameters can take gives rise to a variability of the fKdV equations [7] which are of particular interest in the literature and of which the best known are the Lax equation, the Ito equation, the Caudrey-Dodd-Gibbon (CDG) equation, the Sawada-Kotera (SK) equation and the Kaup-Kuperschmidt (KP) equation. Basically, the KdV equation models the propagation of weakly nonlinear dispersive waves in various fields such as plasma physics [14], waves of small amplitudes and long wavelengths on the surface of shallow waters (incompressible irrotational inviscid fluid) [15] [16] [17], the laser beam propagation [18], and so on. It is fundamental to point out here that the NPDEs admit different types of solutions, the most robust of which are solitary waves (solitons). However, one of the major difficulties remains the determination of the exact, approximate or forced solutions of these equations. In recent decades, numerous works have proposed various methods [19] - [26] of resolution in order to understand and explain phenomena which occur in systems whose dynamics are governed by these NLPDEs. Despite all this progress, a lot of work remains to be done because many other phenomena at the origin of the new predictable or unpredictable behaviors of these different systems still have to be detected, understood and explained in order to guarantee all of humanity a future safely.
So, in this manuscript, One tracks down, using the Bogning-Djeumen Tchaho-Kofané method (BDKm) extended to the new implicit Bogning’ (iB) functions, new prototypes of solitary waves of the standard KdV equation while revealing the hybrid character of these waves. Section 2 will be devoted to the brief description of the BDKm including new implicit Bogning’ functions (iB-functions) while, Section 3 will focus on finding the analytical solutions coupled with an intense numerical simulation. Section 4 will provide the gist of discussions. A conclusion followed by a perspective will complete this work in Section 5.
2. The BDKm Theory
This section takes care of highlighting all the mechanisms necessary for the implementation of the BDKm including iB-functions which will allow in section 3, to unearth new prototypes of the solitary wave solutions of Equation (1).
2.1. iB-Functions
These iB-functions [23] [27] [28] [29] have been highlighted thanks to the multiple research works [22] [23] [28] - [44] produced for more than a decade. It is during the repeated constructions of solitary wave solutions of certain types of equations in wave mechanics presenting dispersion terms coupled with non-linear terms (which can be of different orders) via the BDKm, that the fascinating properties of these functions were detected by Bogning Jean Roger. These iB-functions are noted [23] [27] [28] [29]
(2)
The left member is the implicit form and the right member is the explicit form of the iB-functions, where
are the parameters associated to the independent variables
, m is the power of the numerator, n that of the denominator. In one dimensional, according to the choice of the parameter
, Equation (2), is reduces under the form
(3)
where m and n are keeping the same characteristics as in Equation (2),
is a constant associated to the independent variable x. We associate here some of the fundamental properties of these functions, whose a large majority of which will be useful in the rest of this manuscript under the respective forms
(4)
(5)
(6)
(7)
(8)
(9)
(10)
and
(11)
It is important to note here that this function in its trigonometric form is written as
(12)
although we should not use it for this work, but for the knowledge of the reader. For a better understanding of the properties of these functions, it is needful to refer to [21] [25] [26] [27] where they are widely explained.
2.2. Implementation of the BDKm
The BDKm which has been proposed by three Cameroonian researchers finds its implementation field in nonlinear physics, wave mechanics, mathematics physics, and others. It is better suited for solving certain types of NPDEs of the form [22] [23] [28] - [44]
(13)
where
is an unknown function to be determined, X is some function of
and its derivatives with respect to x and X includes the highest order derivatives and the nonlinear terms. Generally, the solution sought is of the form
. As in the logic of a only variable, we can set the change of variable
. But in the case where we have a function of x and t;
, we can pose the change of variable
. Thus,
becomes
where
is the speed of the wave and Equation (13) becomes in these conditions
(14)
Equation (14) is an ordinary differential equation (ODE), where
,
represent respectively the first and second derivatives of the envelope
with respect to
. According to Equation (3), the solution we are trying to construct can be expressed as
(15)
where
is a real constant and
are the unknown constants to be determined. So, the combination of Equations (14) and (15) gives the main equation
(16)
where
are positive natural integers and
the real numbers [23] [27] [28] [29]. It can be noted here that Equation (16) is the one from which all the possible analyzes result. The identification of coefficients
at zero makes it possible to obtain the ranges of equations whose the resolutions could allow to obtain the expressions of the unknown coefficients
. It is important to point out here that, the resolution of this series of equations often leads to exact solutions [30] [38] for certain models and according to the form of the considered ansatz while, for other models and according to the form of the chosen ansatz , it (resolution) leads to approximate or forced solutions. In the case of approximate or forced solutions, the priority in the order of resolution is given to those from the highest clues of
, then to those of the highest clues of
. But, otherwise we go to those from the coefficients of lowest clues of
and
. Here, the priority makes reference to the series that permits to obtain good results or merely that tends more to the sought exact solution. Very often, the series of equations obtained by identifying at zero the coefficient of
gives satisfaction. On the opposite, the last series of equations given by the coefficient of
is not very important because it is considered rather like a confused domain for the obtainable good solutions. So, this resolution permits to obtain the possible expressions of the coefficients
of Equation (15) as a function of the parameters
, and those supplied by Equation (13), and then of the constraints which may result therefrom. Thus, the ansatz given in Equation (15) can be supported by Equation (13) as a solution. Recently, this method was used in [36] [37] and in [42] [43] [44] to construct hybrid solitary waves for the generalized Kuramoto-Sivashinsky equation, Multi-form solitary wave solutions of the KdV-Burgers-Kuaramoto equation and solitary wave solutions which propagate through transmission media such as electrical lines, respectively. Now, we apply it again to the standard form of the well-known fifth-order KDV equation [7].
3. Results
This part of the work is grouped into two subsections. The first subsection deals with the construction of the analytical hybrid solitary wave solutions of Equation (1). The second sub-section, for its part, is working on an intense numerical simulation in order to reassure itself of the stability of the obtained solutions with a view to a probable future application and for a possible confirmation of the hybrid characters (planned when choosing of the ansatz given by Equation (18) below) of these obtained solutions.
3.1. Analytical Higher Order Solitary Wave Solutions
By considering the change of variable
where
and
the wave speed, Equation (1) becomes
(17)
where
. We are looking for the solutions under the form
(18)
where a,b and c are real constants to be determined later,
, the inverse of the width of each component of Equation (18). It should be noted here that Equation (18) has a peculiarity in the sense that it is a collection of solitons taken individually in the two main families of existing solitary waves, namely the bright family and the dark family to make it a package [23] [27] [28] [29]. To be more precise, the first term of Equation (18) is a bright represented by the sech function, the second term is a hybrid soliton obtained by performing the product of a bright by a dark represented by
and the third term is a dark represented by tanh2. One just has to represent each of the terms to realize it. In this context, we believe that this mixture can be at the origin of interactions (or competitions) between these different components and generate hybrid structures, this, in relation to the values taken by the coefficients
and the parameters
of the wave, as well as those taken by the parameters
of the considered system. Following all that has just been said, we continue with our investigations by inserting Equation (18) into Equation (17) to produce of ranges main equation and contracted under the form
(19)
where
;
. So, Equation (19) has delivered in its formulation, two ranges of equations in the terms of
, and
, thus constituting the most important ranges according to the implementation of the above BDKm. In this context, priority in solving is given to the equations from the first range that is to say, the equations from the coefficients of the terms in
due to the complexity of the equations of the second range. Thus, by identifying at zero the different coefficients
of the terms in
and the different coefficients
of the terms in
, we obtain the both series of algebraic equations of unknowns a,b and c as being
the term in
,
(20)
the term in
,
(21)
the term in
,
(22)
the term in
,
(23)
the term in
,
(24)
the term in
,
(25)
the term in
,
(26)
the term in
,
(27)
the term in
,
(28)
the term in
,
(29)
the term in
,
(30)
the term in
,
(31)
the term in
,
(32)
At look to the structuring of the equations unveiled by this range, one can easily make the choices of the different families of the solutions to be constructed. Thus, in order to obtain non-trivial solutions, we will distinguish four cases from families of solutions, notably the case
; the case
; the case
and the case
. Thus, these investigations plunge us straight into what is essential in Subsection 3.1.
3.1.1. First Family of Solutions: Case: a ≠ 0; b ≠ 0
From Equations (21), (23) and (25), one obtains, respectively
(33)
(34)
and
(35)
with
. By setting the equality
between Equations (33) and (35) on the one hand, then, between Equations (34) and (35) on the other hand, it comes the constraint verified by the parameters
and
below
(36)
The insertion of Equation (36) into one of Equations (33), (34) and (35) gives
(37)
It emerges from Equation (37) that, the coefficient c is a linear function of parameter
or that c is a hyperbolic function of parameter
or that c is a parabolic function of the inverse of the width at half height
of the solitons, that one is searches to constructing according to the fact that the choice of the variable be worn on one of the parameters
or
. These observations sufficiently show what will be the importance of the impact of the variations of the coefficient c on the formation of the wave structures that we will obtain. By continuing, the insertion of Equation (37) into Equation (20) and Equation (26) gives, successively
(38)
(39)
and
(40)
Equation (40) reveals that the speed
has almost the same variations as in the case of Equation (37) with respect to the parameters
and
. Taking into account Equations (37) and (38) in the Equation (22) gives, respectively
(41)
and
(42)
It is important to emphasize here that, the two conditions which validate Equations (39) and (42) impose on to take
and
. And also, Equations (36), (37), (38), (40) and (41) verify Equation (24). We thus obtain the first family of solutions under the form
(43)
with
and
. It is very significant to underline here that, this family of solutions given by Equation (43) through its condition of existence:
, opens the way to other new varieties of fifth-order KdV equations [5] including her will be solution.
3.1.2. Second Family of Solutions: Case:
Under this condition, Equations (21), (23) and (25) are verified. However, Equation (26) gives
(44)
while Equation (20) gives, successively
(45)
and
(46)
with
. Next, Equations (22) and (24) become respectively
(47)
and
(48)
On the one hand, by taking into account of Equation (45) in the Equation (47), and on the other hand, the insertion of Equations (44) and (45) in Equation (48) lead to the same quadratic equation in c which is written
(49)
with
and
. Equation (49) has as discriminant
(50)
Looking at the architecture of Equation (50), a discussion according to the constraints linking
and
becomes necessary. So, this discussion concerns the cases
and the case
. From this investigation, the second family of solutions gives rise to two subfamilies.
1) First subfamily of the second family of solutions: case:
Due to the fact that
, the discriminant becomes
because
. Therefore, Equation (49) admits as solution
(51)
Equation (51) shows that, for a given value of
, the coefficient c is a parabolic function of the width at half height of each component of Equation (18). Thenceforward, the first subfamily of the second family of solutions is written considering Equation (46) under the form
(52)
where c is given by Equation (51). This solution presages interactions between hybrid and dark components of Equation (18). Which, in view of the constitution of the hybrid component of Equation (18), could generate bright, bright-dark, dark-bright, dark type structures, and so on, provided that
.
2) Second subfamily of the second family of solutions: case:
Under these conditions, the discriminant is that given by Equation (50). So, for
, the expression of c solution of Equation (49) that we are looking for is in the form
(53)
Therefore, the second subfamily of the second family of solutions is written considering Equation (46) as being
(54)
where c is given by Equation (53). This solution presages interactions between hybrid and dark components of Equation (18). Which, in view of the constitution of the hybrid component of Equation (18), could generate bright, bright-dark, dark-bright, dark type structures, and so on, provided that
. We cannot forget to point out, in the case of this second family of solutions of Equation (1) that:
- the first subfamily given by Equation (52) through its condition of existence:
, opens the way to other new varieties of fifth-order KdV equations [5] including her will be solution;
- the second subfamily given by Equation (54) through its condition of existence:
and
, opens the way to other new varieties of fifth-order KdV equations [5] including her will be solution. Now, One looks what is about the third family of solutions of the Equation (1).
3.1.3. Third Family of Solutions: Case:
For
and
, all the equations from (20) to (26) are verified for any nonzero values of the real coefficients
and those of the parameters
and
. According to the implementation of the BDKm as explained in Section 2, the third family of solutions can be approximated by the following expression
(55)
with
and
. Equation (55) presages interactions between bright and hybrid components of Equation (18). Which, in view of the constitution of the hybrid component of Equation (18), could generate bright, bright-dark, dark-bright, dark type structures, and so on, for nonzero values of
and those of the parameters
and
. To fix the a priori ideas, a study of the numerical stability of the obtained solutions is necessary.
3.1.4. Fourth Family of Solutions: Case:
By considering this condition, Equations (20) to (26) are verified. However, according to the BDK method, it is important in this case to examine what information can we derive from the equations of the second range (terms in
). In this context, Equations (28), (30) and (32) are checked while, Equations (27), (29) and (31) are reduced to, respectively
(56)
(57)
and
(58)
Equation (56) is a quadratic equation of the unknown coefficientc and the discriminant of which is given by
(59)
From the architecture of Equation (59), the discussions according to the constraints between the parameters
and
is necessary and thus offering two subfamilies of solutions for this fourth family, in particular, the subfamily characterized by
and the subfamily characterized by
.
1) First subfamily of the fourth family of solutions: case:
Under the condition
, Equation (59) leads to
and then, Equation (56) has as solution
(60)
The insertion of Equation (60) in Equations (57) and (58), and, taking into account the constraint
give successively
(61)
and
(62)
So, we get this first subfamily of the fourth family of solutions of Equation (1) considering Equation (61) in the form
(63)
Equation (63) shows that the fKdV equation has dark solitary wave solutions with a amplitude which is a parabolic function of the width at half height of the soliton for a given value of
.
2) Second subfamily of the fourth family of solutions: case:
For
and
, the solution of Equation (56) is written
(64)
Then, Equation (58) gives
(65)
where c is given by Equation (64) and should verify Equation (57). Thus, the second subfamily of the second family of solutions of Equation (1) has as expression
(66)
where
. As in the case of Equation (63), the result given by Equation (66) also indicates that Equation (1) admits dark type solitons as solutions.
It should be said here that the solutions given by Equations (63) and (66) are exact solutions of Equation (1) of which one of the classes had been proposed in [12] through the tanh method. We can retain here that, for certain solutions constructed in this subsection by using the discriminant
, only the cases
are interesting because complex solutions are not part of the aims of this work.
3.2. Numerical Simulations
This section is dedicated to the numerical simulations aiming to make observable, reliable and applicable the obtained solutions. They (numerical simulations) also aim to confirm the theoretical predictions on the hybrid characters of the new prototypes of the obtained solitary waves within the framework of this work. To implement all of this, and knowing that the boundary conditions of the profiles of the solutions constructed in this manuscript are not identical at the two borders. We have used the MATLAB toolbox pdepe [45] which solves initial-boundary value problems for parabolic-elliptic PDEs in 1-D, with zero flux boundary conditions. We have also used spatially extended grids in others to minimize boundary reflections that could induce spurious effects. It is also necessary to point out here that, these boundary conditions are appropriate to the profiles of the solutions studied in this work, instead of the periodic boundary conditions which require that when a wave passes from one end (of the computational spatial grid) to the other which is opposite to it, it should keep the same properties.
We cannot continue this study without providing the reader with elements useful for a good understanding of the different spatiotemporal evolution curves obtained. For example, to obtain Figure 1(a), we have set
to infer
of Equations (36), (40) and (41) respectively. In this context, discussions become necessary.
4. Discussions
From all that has been formulated in the previous lines, two poles clearly emerge and on which our discussions will focus. The first pole will be on the analytical results, while, the second pole will revolve around results revealed by numerical simulations. Thus:
*Speaking of analytical results, we can see that the existence conditions of the solutions given by Equations (43); (52) and (54), respectively, offer possibilities
Figure 1. Robust dynamics for the hybrid solitary waves given by Equation (43). (a):
;
;
;
;
;
;
;
. (b):
;
;
;
;
;
;
;
. (c):
;
;
;
;
;
;
;
. (d):
;
;
;
;
;
;
;
.
to develop varieties of the fifth-order KdV equations by simply modifying [7] the real values of the parameters
and
. Note here that all these varieties obtained must admit the ansatz given by Equation (18) as an approximate solution. To be a little clearer, one of the well-known fifth-order KdV equations is that of Lax, and which is characterized by:
(67)
For example, if
, we obtain
and
, successively. Thus, Equation (1) becomes:
(68)
Equation (68) is indeed the Lax equation which is one of the varieties of equations which derive from Equation (1) with values of parameters
and
given above. If we go back to our results:
- the solution given by Equation (43) has the existence condition
. For
, we have
and
. These values taken by the coefficients
and
give rise to another equation which derives from Equation (1) and which is written:
(69)
and characterized by:
and whose Equation (18) is an approximate solution.
- next, the solution given by Equation (52) has for existence condition
. Let us take
in the form:
. So, for
, we have
. So we get another form of Equation (1) as
(70)
whose solution given by Equation (18) is an approximate solution.
- eventually, a similar reasoning can also be carried out with the existence condition
of the solutions given by Equation (54), in order to obtain the corresponding fifth-order KdV equations. It is necessary to point out here that the five well-known fifth order KdV equations in the literature can be tidy in this latter case.
It should be point out here that, Equations (69) and (70) are two new forms [7] among so many other fifth-order KdV equations that we have just revealed through this manuscript and which are characterized by
and
, respectively, and whose Equations
(43) and (52) are respective approximate solutions. From these existence conditions, we can propose new varieties of fifth-order KdV equations by simply altering the values of the different parameters
and
.
*As for the numerical results, it appears by browsing from Figures 1-5 that:
- Figure 1 displays a set of four hybrid structures that we think are generated by balanced interactions between the three constituent solitons of the soliton package [23] [27] [28] [29] given by Equation (18). Figure 1(a) and Figure 1(c) are unfinished intermediate forms of a double-bright soliton while, Figure 1(b) and Figure 1(d) are unfinished intermediate forms of a double-dark soliton respectively, obtained by simply altering the values of the coefficients a, b and c.
- Figure 2 shows the dominant interactions of each of the terms making up Equation (18). We can observe this by noting that Figure 2(a) translates the domination of the bright soliton structure (first term of Equation (18)) while, Figure 2(b) and Figure 2(c) translate the domination of the hybrid soliton structure (second term of Equation (18)). Figure 2(d) translates the domination of the dark-soliton structure (third term of Equation (18)).
- Figure 3 presents a set of four profiles which reveal through Figure 3(a) and
Figure 2. Stable dynamics of the bright-dark solitary waves given by Equation (43). (a):
;
;
;
;
;
;
;
. (b):
;
;
;
;
;
;
;
. (c):
;
;
;
;
;
;
;
. (d):
;
;
;
;
;
;
;
.
Figure 3(b) a dominance of hybrid structures with strong bright-soliton tendencies and given by Equation (52), while, Figure 3(c) and Figure 3(d) reveal a dominance of hybrid structures with strong dark-soliton tendencies and given by Equation (54). All these structures result from the interactions between the second term (hybrid soliton) and the third term (dark soliton) of Equation (18) in the absence of the first term (bright soliton).
- Figure 4 reveals four profiles which present improved and balanced forms of the structures obtained in Figure 1. Figure 4(a) is an intermediate form of a double-bright soliton with balanced vertices while Figure 4(b) is an intermediate form of a double-dark soliton with balanced bottoms. On the other hand, Figure 4(c) presents an intermediate form of the double-bright soliton very evolved towards a bright-soliton with flat peak, then, Figure 4(d) is an intermediate form very evolved towards a dark-soliton form with flat bottom. One also has to point out here that all these figures are symmetrical with respect to a plane of equation
and that these properties of symmetry
Figure 3. Stable spatiotemporal evolution of the hybrid bright-dark solitary waves given by Equation (52) and (54). (a):
;
;
;
;
;
;
;
; (b):
;
;
;
;
;
;
;
; (c):
;
;
;
;
;
;
;
; (d):
;
;
;
;
;
;
;
.
could be of capital importance in the understanding and the explanation of certain phenomena which occur in systems whose dynamics are described by Equation (1).
- Finally, Figure 5 completes the advanced forms presented in Figure 4. This is how Figure 5(a) shows a bright soliton structure with a flat peak and Figure 5(b) shows a dark-soliton structure with a flat bottom. Figure 4 and Figure 5 translate the interactions between the first term (bright soliton) and the third term (dark soliton) of Equation (18) in the absence of the second term (hybrid soliton).
In summary, we find that, the analytical solutions obtained and of which we foresee hybrid characters have given rise to new varieties of fifth-order KdV equations while numerical simulations have confirmed the hybrid characters of most of the structures revealed starting from Figures 1-5. These solutions may find their applications in fluid mechanics, quantum mechanics, nonlinear optics, and so on. Abdul-Majid Wazwaz in [8] had used the simplified Hirota method to derive single, singular and multiple solitons with the auxiliary function, of
Figure 4. Robust dynamics of the multi-form solitary waves given by Equation (55). (a):
;
;
;
;
;
;
;
. (b):
;
;
;
;
;
;
;
. (c):
;
;
;
;
;
;
;
. (d):
;
;
;
;
;
;
;
.
Figure 5. Stable spatiotemporal evolution of bright and dark solitons with a flat top and a flat bottom, respectively, given by Equation (55). (a):
;
;
;
;
;
;
;
. (b):
;
;
;
;
;
;
;
.
Equation (1). He had also derived periodic solutions and compacton solutions of Equation (1) in [13]. Again, Abdul-Majid Wazwaz had used in [12] the sine-cosine and the tanh methods to successively derive exact periodic and soliton solutions of Equation (1). Bell shape multi-solitons, dark solitons and dark periodic cusp solutions of Equation (1) were found in [10] [11] while, exponential function solutions were revealed in [9]. It becomes easy to realize that, all these proposed solutions of Equation (1) including solutions found in [19] [25], are almost entirely different from those obtained in this manuscript.
5. Conclusion
In the final analysis, we can say with enthusiasm that, the application of the BDKm extended to iB-functions, to the standard form of the fifth-order KdV equations, has enabled to reveal new prototypes of solitary waves of higher order, which are for some, approximate hybrid (or multi-form) solutions, and for others, exact solutions of this equation. This method is very suitable for the construction of solitary wave solutions of certain types of equations in wave mechanics and which have terms of dispersion and terms of nonlinearity. The existence conditions of some of these solutions made it possible to derive from Equation (1) (and depending on the values of the parameters and depending on the values of the parameters
and
); some new varieties of the standard fifth-order KdV equations. Equation (18) is one of their multiple approximated solutions. Theoretical forecasts on the hybrid characters of the solutions obtained have been confirmed by intense numerical simulations as evidenced by the figures in Subsection 3.2. These figures also made it possible to understand that, when solitons are put together in the same package [23] [27] [28] [29], interactions occur between the different components of this package of solitons and what subsequently generates structures with hybrid characters. One estimates that, new varieties of equations derived from that of standard fifth-order KdV, including the new solutions proposed in this manuscript, will allow, during the propagations tests in the laboratories, to analyze, understand and explain some old or new phenomena which occur or which would occur in the systems whose dynamics are governed by Equation (1), in particular during the propagation of waves of low amplitudes and long wavelength on the surface of shallow waters [15] [16], during the transport of data through non-linear optical fibers and many other phenomena in quantum mechanics, in order to better secure our existence. One also hopes that, the properties of symmetry which certain obtained structures have, could be of capital importance in the understanding and the explanation of certain phenomena which occur in systems whose dynamics are described by Equation (1). It should also be pointed out that, the possibilities of formulating new varieties of fKdV-type equations offered by the results of this manuscript will help in the future, to devise new shallow fluid media with improved properties. However, the universe as a whole is being perpetually dynamic. We must further explore these models in order to obtain every day new information which is still hidden behind this equation for the benefit of all humanity and which will guarantee at the same time our survival.
Acknowledgements
The authors would like to thank both ministries of higher education of Cameroon and Gabon for their assistance through their respective research support programs. The authors’ thanks go again to Keutcha Padji Aurélien whose assistance in reading this manuscript has contributed to improving the quality of the English used. Authors would also like to thank the reviewers whose relevance of their criticism has enabled us to more improve the quality of this work.