# A Conservation Law Treatment of Nonlinear KdV Hierarchies

Authors: Sameerah Jamal; Rivoningo Maphanga
Journal: Journal of Mathematics
 (2022)

---

## Abstract

We study the hierarchy commonly defined as an infinite sequence of partial differential equations which begins with the Korteweg–de Vries equation and its modified version. An important feature of the hierarchy is its highly nonlinear property. In this regard, obtaining solutions for the members of the hierarchy poses a great problem. In this paper, we propose a method to allow for the construction of solutions to the full hierarchy. Our approach involves a recursion operator in the conservation law of the hierarchy. The efficiency of the method is demonstrated by selected examples. In certain cases, we obtain snoidal solutions.


---

## Body

## 1. Introduction

The Korteweg–de Vries (KdV) equation,(1)ut=uxxx−6uux,attributed to the Dutch mathematicians Diederik Korteweg and Gustav de Vries (1895), is a mathematical model of waves on shallow water surfaces whose origin has a long history. To this day, the KdV equation is still considered as one of the most important nonlinear partial differential equations as it has a wide variety of applications [1]. It was originally derived to model the propagation of weakly dispersive nonlinear water waves and serve as a model equation for any physical system in consideration [2]. In most contexts, ux,t is the function that denotes the elongation of the wave at space x and time t [3].KdV equation (1) naturally extends to an infinite sequence of integrable nonlinear partial differential equations of solitonic characters [4] and can be considered as the initiator of the KdV hierarchy, denoted by(2)ut=Dx2−4u−2uxDx−1nux,n=1,2,…,where Dx denotes the total derivative and Dx−1 denotes the integral with respect to x. The second member of the hierarchy is given when n=2 in (2), viz.(3)ut=uxxxxx−20uxuxx−10uuxxx+30u2ux.Note that, for higher ordern, the equations become increasingly nonlinear and higher order in derivatives.If the nonlinear termuux in equation (1) is replaced by umux, then the most important case other than when m=1, is when m=2. This yields the so-called modified KdV (mKdV) equation, given by [5, 6](4)ut=uxxx−6u2ux,which may also be connected to the mKdV hierarchy:(5)vt=Dx2−4v2−4vxDx−1vnvx.Here, the case ofn=2 gives the second member of the mKdV hierarchy:(6)vt=30vxv4−10v2vxxx−40vvxvxx−10vx3+vxxxxx.Another famous hierarchy is the Burgers’ hierarchy and its solutions [7, 8].Completely integrable nonlinear equations, such as those of the KdV hierarchy, are endowed with many special mathematical properties. They are of interest due to their infinite conservation laws [9] and symmetries [10], bi- or tri-Hamiltonian structures, their Painlevé property [11], Lax pairs [12], etc. The solving of such equations is deeply connected to the inverse scattering transform [13, 14] and Hirota’s direct method [15]. The aforementioned importance of the KdV hierarchy has motivated this study.The mKdV equations are related to the KdV equation through the Miura transformation, which maps the solutions of the KdV equations to the solutions of the mKdV equations [16].The KdV hierarchies are an example of higher-order water wave models, which are of great significance, and they play crucial roles especially to study physical systems and the necessary material properties needed to manipulate waves in a desired manner [17].In this paper, we propose a method to solve the full KdV hierarchy and extend the method to include the mKdV hierarchy. This is a novel proposal as, to the best of our knowledge, no such endeavour has appeared in the literature. The entire hierarchy is highly nonlinear and of higher order in derivatives, thereby posing an extremely challenging problem to solve. We formulate our method based on a transformation of variables (derived from the point symmetries of the hierarchy), effectively reducing the partial differential hierarchy into an ordinary differential hierarchy. The latter is connected to a transformed conservation law of the hierarchy, and the knowledge of this transformed conservation law forms a general approach to solving the hierarchy for alln.This paper is organised as follows. In Section2, we briefly present some theoretical considerations, and Section 3 discusses some general properties of the hierarchies. Section 4 contains the main results and a description of our method, and Section 5 elaborates on some applications of our method. Section 6 provides some alternate solutions.


## 2. Preliminaries

The procedure for determining point symmetries for an arbitrary system of equations is well known [18]. Consider q unknown functions uα which depend on p independent variables xi, i.e., u=u1,…,uq and x=x1,…,xp, with indices α=1,…,q and i=1,…,p. Let(7)Gαx,uk=0be a system of nonlinear differential equations, where uk represents the kth derivative of u with respect to x. We consider the following symmetry:(8)X=ξi∂xi+ηα∂uα,given by(9)XGαx,uk=0,whenGαx,uk=0,where X is extended to all derivatives appearing in the equation through an appropriate prolongation. A current T=T1,…,Tn is conserved if it satisfies(10)DiTi=0,along the solutions of the given equation. Equation (10) is called a local conservation law.Suppose thatX is a symmetry of system (7) and T is a conserved vector of (7). Then, if X and T satisfy(11)XTi+TiDjξj−TjDjξi=0,the symmetry X is said to be associated with T [19]. Equation (11) is closely related to a Noether theorem, but in [19], it was proved that this result holds without the existence of Lagrangian. The transformation u=νx for the KdV equation enables the construction of a Lagrangian density [20] so that Noether’s theorem may be applied for its conservation laws. However, we have opted to study the KdV equation in the absence of Lagrangian.


## 3. Generalised Properties of the KdV and mKdV Hierarchy

A standard calculation of the symmetries of equation (1), using condition (9), reveals that it has the following four Lie point symmetries:(12)X1=∂t,X2=∂x,X3=6t∂x+∂u,X4=3t∂t−2u∂u+x∂x.The second member of the hierarchy, withn=2, or equation (3) has the following three symmetries:(13)X1,X2,X32=5t∂t+x∂x−2u∂u.If one repeats the Lie symmetry method for higher members, it is easy to see that hierarchy (2) possesses the Lie point symmetries(14)X1,X2,X1n=2n+1t∂t−2u∂u+x∂x,for n≥2 with Lie bracket relations A,B=AB−BA, given in Table 1.Table 1
Lie brackets of symmetries (14).
[,]X1X2X1nX1002n+1X1X200X2X1n−2n+1X1−X20A similar investigation of mKdV hierarchy (5) gives that n=1 has the following three symmetries:(15)X1,X2,X6=x3∂x+t∂t−v3∂v,and n=2 has the symmetries(16)X1,X2,X7=x5∂x+t∂t−v5∂v.As before, if one repeats the Lie symmetry method for higher members, it is easy to see that hierarchy (5) possesses the Lie point symmetries(17)X1=∂t,X2=∂x,X2n=t∂t−v2n+1∂v+x2n+1∂x,for n≥2 with Lie bracket relations in Table 2.Table 2
Lie brackets of symmetries (17).
[,]X1X2X2nX100X1/2n+1X200X2X2n−X1/2n+1−X20As for the conservation laws of the above hierarchies, we notice several interesting properties. There exist many ways to compute conservation laws, and we opt for the multiplier approach [21].Below are the cases of the KdV hierarchy whenn=1 and n=2, where the conservation laws are Ti=Tit,Tix where i=1,2,and3. Equation (1) has the following three conservation laws:(18)T1t=u,T1x=3u2−uxx,(19)T2t=u22,T2x=2u3+12ux2−uuxx,and finally,(20)T3t=12tu2−x6u,T3x=−ux6−tuuxx+2tu3−x2u2+t2ux2+x6uxx.These conservation laws can also be found via Noether’s theorem if the problem is reformulated to possess Lagrangian. The calculations are straightforward but tedious. Reports of these quantities or their equivalent appear in [20, 22], with more in [23].Equation (3), i.e., n=2, has the following two conservation laws:(21)T1t=u,T1x=−10u3+10uuxx+5ux2−uxxxx,(22)T2t=u22,T2x=−152u4+10u2uxx−uuxxxx+uxuxxx−12uxx2.As for the mKdV equation, equation (4) has the following two conservation laws:(23)T1t=12v2,T1x=−vvxx+12vx2+32v4,(24)T2t=v,T2x=2v3−vxx.Similarly, equation (6), n=2, has the conservation laws(25)T1t=12v2,T1x=−5v6+10vxxv3+5v2vx2−vvxxxx+vxvxxx−12vxx2,(26)T2t=v,T2x=−6v5+10v2vxx+10vvx2−vxxxx.We now establish a result that is the foundation of our approach and that is thenth conservation law of each hierarchy. As we shall show, such a conservation law can be manipulated to solve the entire hierarchy for all values of n. To begin, we establish the nth conserved vector of hierarchy (2) by the following theorem.Theorem 1.
The KdV hierarchy possesses the conserved vectorT=Tt,Tx=u,−Dx−1Dx2−4u−2uxDx−1nux along the solutions of equation (2), i.e., a component of the conserved vector admits a recursion operator.Proof.
Suppose the conservation law is(27)DtTt+DxTx=0,where Tt is the conserved density and Tx is the conserved flux. Then, from equation (2), we have(28)ut+Dx2−4u−2uxDx−1nux=Dtu+Dx−Dx−1Dx2−4u−2uxDx−1nux=0,along the solutions of equation (2), and the result follows.
Hence, the conserved density for every member of the hierarchy isTt=u, while the conserved flux is Tx=−Dx−1Dx2−4u−2uxDx−1nux. The fluxes for the first few members of hierarchy (18) for n=1 and (21) for n=2 are confirmed by the above theorem.
Similarly, we can prove a result for the mKdV hierarchy.Theorem 2.
The mKdV hierarchy possesses the conserved vectorT=Tt,Tx=v,−Dx−1Dx2−4v2−4vxDx−1vnvx along the solutions of equation (5), i.e., a component of the conserved vector admits a recursion operator.One can easily check that the conserved vectors (24) for n=1 and (26) for n=2 arise from this theorem.In the next section, we give a method to find solutions of the entire hierarchy.


## 4. A Method to Solve the Full Hierarchy-Type I Solutions

To proceed, we require a symmetry generatorX to be associated with a conserved vector T of a given equation. Based on the previous section, (14), and (17), we notice that the symmetries X1 and X2 are possessed by the respective hierarchies for all n. In particular, we observe that the point symmetry ∂t, when applied to condition (11), is(29)∂tu−Dx−1Dx2−4u−2uxDx−1nuxand for the second symmetry,(30)∂xu−Dx−1Dx2−4u−2uxDx−1nuxTherefore, both symmetries satisfy the association condition, and we conclude that they are associated with every conservation law of Theorem1, i.e., with any of the nth conservation law of the KdV hierarchy member. A similar result holds for the mKdV hierarchy, and here, we conclude that the same symmetries are associated with every conservation law of Theorem 2.Next, we recall the fundamental theorem on double reduction [24, 25], which states that there exist functions Tr such that(31)DrTr=0.The transformed conserved quantity may be expressed as(32)Tr=TtDtr+TxDxrDtrDxs−DxrDts,where r and s are similarity variables connected to an associated symmetry X.SinceX1 and X2 are associated with the conserved vector T, we consider the linear combination X=X2+cX1 (c is a constant) to obtain the similarity transformation(33)r=cx−t,s=x,ux,t=ur,for the KdV hierarchy, and similarly,(34)r=cx−t,s=x,vx,t=vr,for the mKdV hierarchy.Therefore, we may establish the following results forTr.Theorem 3.
The conserved quantity of KdV hierarchy equation (2) can be reduced to(35)u+cDr−1c2Dr2−4u−2urDr−1nur,where u=ur.Proof.
Application of (32) gives us(36)Tr=u−1−Dx−1Dx2−4u−2uxDx−1nuxc−1and in the new variables, by transformation (33), equation (36) transforms to(37)Tr=u+cDr−1c2Dr2−4u−2urDr−1nur.
As examples,Tr corresponding to T1 of (18) is given by(38)Tr=urr⋅c3−3u2c+u,for n=1, and Tr for (21) in the case of n=2 is(39)Tr=urrrr⋅c5−10uurr⋅c3−5ur2⋅c3+10u3⋅c+uTheorem 4.
The conserved quantity of mKdV hierarchy equation (5) can be reduced to(40)v+cDr−1c2Dr2−4v2−4vrDr−1vnvr,where v=vr, and the proof is similar to that of Theorem 3.
Also, to this end, examples ofTr for the mKdV hierarchy include T2 of (24) given by(41)Tr=v−2v3⋅c+vrr⋅c3,for n=1, and Tr for T2 of (26) is given by(42)Tr=v+6v5⋅c−10v2vrr⋅c3−10v2vr2⋅c3+vrrrr⋅c5,for n=2.
That is, the above results can be used to findTr for any value of n, for both KdV and mKdV hierarchies. Based on equation (31), we have that Tr=κ, κ is a constant. Therefore, we have reduced the entire partial differential KdV and mKdV hierarchies to ordinary differential hierarchies. These ordinary differential hierarchies may then be solved for any n.


## 5. Type I Solutions

In this section, we illustrate the applicability of the above method and theory in establishing solutions to members of the KdV and mKdV hierarchy. The solutions obtainable via our method in Section4 will be referred to as type I solutions. Below, we set κ=0 for simplicity.


### 5.1. The KdV Hierarchy

Let us considern=1 in Theorem 3 to get the reduced conserved component (38); that is, we solve(43)urr⋅c3−3u2c+u=0.We find that this equation has an implicit solution(44)±∫urc21cC1c3+2a3c−a2da−r−C2=0.Suppose we choose the free parameters to beC1=1,C2=0,andc=1/2; then, the integral in (44) is evaluated to be(45)±225−1−2u+1−4u+5+14u+5−1EllipticF25+1−2u+1,2i5−118u3−8u2+1,where EllipticF is the incomplete elliptic integral of the first kind, and the solution to (43) becomes(46)ur=−25−28JacobiSN25−2r2,i2+i252−54+−10+1058+12,or in the original independent variables, by reversing transformation (33),(47)ux,t=−25−28×JacobiSN25−2cx−t2,i2+i252−54+−10+1058+12,where JacobiSN is an inverse of elliptic integrals and doubly periodic elliptic functions. These solutions appear graphically in Figure 1.Figure 1
2D and 3D evolutions of solutions from (46) and (47), respectively.Next, we consider the second member of the hierarchy,n=2, in Theorem 3 to get the reduced conserved component (39). That is, we solve(48)urrrr⋅c5−10uurr⋅c3−5ur2⋅c3+10u3⋅c+u=0.In this case, we find two solutions of the second member of the hierarchy, viz.(49)ux,t=−isech2−C1+1/2+i/2cx−t/2c5/42c,or secondly,(50)ux,t=i55sech2−C1+−i/c5/2−i5/c5/2cx−t/224+5sech2−C1+−i/c5/2−i5/c5/2cx−t/224−25102c.


### 5.2. The mKdV Hierarchy

This time, let us considern=1 in Theorem 4 to get the reduced conserved component (41); that is, we solve(51)v−2v3⋅c+vrr⋅c3=0.The solution of (51) is of two cases, namely,(52)vr=±2i2c3C1−c1−1−4c4C1×JacobiSN2cC1r2/1−1−4c4C1+2cC1C2r/1−1−4c4C1+cC1C22/1−1−4c4C11−1−4c4C1/1−4c4C1+11+1−4c4C1,and secondly,(53)vr=C2cC22−c+1JacobiSNr/c2−cc−1+C1cC22−c+1,C2c−1−cc−1.The latter may be expressed in original variables as(54)vx,t=C2cC22−c+1JacobiSNcx−t/c2−cc−1+C1cC22−c+1,C2c−1−cc−1.This solution has 2D and 3D plots in Figure2.Figure 2
Evolution of solutions corresponding to (53) and (54), respectively.


## 5.1. The KdV Hierarchy

Let us considern=1 in Theorem 3 to get the reduced conserved component (38); that is, we solve(43)urr⋅c3−3u2c+u=0.We find that this equation has an implicit solution(44)±∫urc21cC1c3+2a3c−a2da−r−C2=0.Suppose we choose the free parameters to beC1=1,C2=0,andc=1/2; then, the integral in (44) is evaluated to be(45)±225−1−2u+1−4u+5+14u+5−1EllipticF25+1−2u+1,2i5−118u3−8u2+1,where EllipticF is the incomplete elliptic integral of the first kind, and the solution to (43) becomes(46)ur=−25−28JacobiSN25−2r2,i2+i252−54+−10+1058+12,or in the original independent variables, by reversing transformation (33),(47)ux,t=−25−28×JacobiSN25−2cx−t2,i2+i252−54+−10+1058+12,where JacobiSN is an inverse of elliptic integrals and doubly periodic elliptic functions. These solutions appear graphically in Figure 1.Figure 1
2D and 3D evolutions of solutions from (46) and (47), respectively.Next, we consider the second member of the hierarchy,n=2, in Theorem 3 to get the reduced conserved component (39). That is, we solve(48)urrrr⋅c5−10uurr⋅c3−5ur2⋅c3+10u3⋅c+u=0.In this case, we find two solutions of the second member of the hierarchy, viz.(49)ux,t=−isech2−C1+1/2+i/2cx−t/2c5/42c,or secondly,(50)ux,t=i55sech2−C1+−i/c5/2−i5/c5/2cx−t/224+5sech2−C1+−i/c5/2−i5/c5/2cx−t/224−25102c.


## 5.2. The mKdV Hierarchy

This time, let us considern=1 in Theorem 4 to get the reduced conserved component (41); that is, we solve(51)v−2v3⋅c+vrr⋅c3=0.The solution of (51) is of two cases, namely,(52)vr=±2i2c3C1−c1−1−4c4C1×JacobiSN2cC1r2/1−1−4c4C1+2cC1C2r/1−1−4c4C1+cC1C22/1−1−4c4C11−1−4c4C1/1−4c4C1+11+1−4c4C1,and secondly,(53)vr=C2cC22−c+1JacobiSNr/c2−cc−1+C1cC22−c+1,C2c−1−cc−1.The latter may be expressed in original variables as(54)vx,t=C2cC22−c+1JacobiSNcx−t/c2−cc−1+C1cC22−c+1,C2c−1−cc−1.This solution has 2D and 3D plots in Figure2.Figure 2
Evolution of solutions corresponding to (53) and (54), respectively.


## 6. Type II Solutions

As seen above, both hierarchies admit conservation laws, such as (19) or (22), independent of Theorems 1 and 2. Now, we cannot transcribe these conservation laws to theorems with a recursion operator as was done in Theorems 1 and 2. Nonetheless, a Tr function may still be obtained in such cases, using the same formula (32) and transformation (33) or (35). This will lead to other solutions, which we call type II.For example, in the KdV hierarchy,Tr for T2 of (19) is given by(55)Tr=uurr⋅c3−12ur2⋅c3−2u3c+u22,for n=1, and Tr for T2 of (22) is given by(56)Tr=152u4⋅c−10uurr⋅c3+uurrr⋅c4−ururrr⋅c5+12urr2⋅c5+u22,for n=2.As for the mKdV hierarchy, forn=1, we have Tr for T1 of (23) which is given by(57)Tr=12v2+vvrr⋅c3−12vr2⋅c3−32v4⋅c,and Tr for T1 of (25) is given by(58)Tr=12v2+5v6⋅c−10v3vrr⋅c3−5v2vrr2⋅c5+vvrrrr⋅c5−vrvrrr⋅c5+12vrr2⋅c5,for n=2. Below, we explore some solutions that arise out of these Tr functions.


### 6.1. Type II Solution to the KdV Hierarchy

A type II solution corresponding to solving (55) yields an implicit solution(59)±∫urc21caC1c3+2_a2c−ada−r−C2=0.Here, the above integral is equal to(60)12c−1+−8C1c4+12curC1c3+2ur2c−ur×4cur+−8C1c4+1−1−1+−8C1c4+1×−4cur+−8C1c4+1+1−8C1c4+1×−cur−1+−8C1c4+1×EllipticF4cur+−8C1c4+1−1−1+−8C1c4+1,122−1+−8C1c4+1−8C1c4+1.Suppose we letC1=1,C2=0,andc=1/2; then, the explicit solution to (59) is(61)ur=−2+24×JacobiSNr24,i22−2242−1.or in original variables,(62)ux,t=−2+24×JacobiSNcx−t24,i22−2242−1.The progression of these solutions appears in Figure3, and they are visibly periodic in nature.Figure 3
2D and 3D progression of solutions corresponding to (61) and (62), respectively.


### 6.2. Type II Solution to the mKdV Hierarchy

A type II solution corresponding to solving (57) has an implicit solution(63)±∫vrc21caC1c3+a3c−ada−r−C2=0.The integral is evaluated to be(64)12c−1+−8C1c4+12cvrC1c3+2vr2c−vr×4cvr+−8C1c4+1−1−1+−8C1c4+1×−4cvr+−8C1c4+1+1−8C1c4+1×−cvr−1+−8C1c4+1×EllipticF4cvr+−8C1c4+1−1−1+−8C1c4+1,122−1+−8C1c4+1−8C1c4+1.Suppose we letC1=1,C2=0,andc=1/2; then, the solution to (63) is(65)vr=±1tanrtanr2+1,or(66)vx,t=±csccx−t,whose graphical representation appears in Figure 4.Figure 4
2D and 3D plots corresponding to (65) and (66), respectively.


## 6.1. Type II Solution to the KdV Hierarchy

A type II solution corresponding to solving (55) yields an implicit solution(59)±∫urc21caC1c3+2_a2c−ada−r−C2=0.Here, the above integral is equal to(60)12c−1+−8C1c4+12curC1c3+2ur2c−ur×4cur+−8C1c4+1−1−1+−8C1c4+1×−4cur+−8C1c4+1+1−8C1c4+1×−cur−1+−8C1c4+1×EllipticF4cur+−8C1c4+1−1−1+−8C1c4+1,122−1+−8C1c4+1−8C1c4+1.Suppose we letC1=1,C2=0,andc=1/2; then, the explicit solution to (59) is(61)ur=−2+24×JacobiSNr24,i22−2242−1.or in original variables,(62)ux,t=−2+24×JacobiSNcx−t24,i22−2242−1.The progression of these solutions appears in Figure3, and they are visibly periodic in nature.Figure 3
2D and 3D progression of solutions corresponding to (61) and (62), respectively.


## 6.2. Type II Solution to the mKdV Hierarchy

A type II solution corresponding to solving (57) has an implicit solution(63)±∫vrc21caC1c3+a3c−ada−r−C2=0.The integral is evaluated to be(64)12c−1+−8C1c4+12cvrC1c3+2vr2c−vr×4cvr+−8C1c4+1−1−1+−8C1c4+1×−4cvr+−8C1c4+1+1−8C1c4+1×−cvr−1+−8C1c4+1×EllipticF4cvr+−8C1c4+1−1−1+−8C1c4+1,122−1+−8C1c4+1−8C1c4+1.Suppose we letC1=1,C2=0,andc=1/2; then, the solution to (63) is(65)vr=±1tanrtanr2+1,or(66)vx,t=±csccx−t,whose graphical representation appears in Figure 4.Figure 4
2D and 3D plots corresponding to (65) and (66), respectively.


## 7. Conclusions

In the study of differential equations, equations that are highly nonlinear and that possess higher-order derivatives are almost impossible to solve. We have proposed a scheme to overcome this problem and aid the solution of, in particular, the KdV and mKdV infinite hierarchy.The well-known (solitary wave) solution of the KdV equation involves the hyperbolic secant function [26], but Korteweg and de Vries were interested in cnoidal solutions, expressible in terms of Jacobi’s elliptic CN functions [27]. Given the mathematical relations between the Jacobi SN and CN solutions, sin,cos, sech, and tanh functions, our above solutions for KdV n=1 may be related to the known ones, but in that case, the recovery of known or related solutions validates our approach. As for our solutions for KdV n=2, we find no connection to any known results. It is possible to find many more solutions.In the analysis of our solutions, we divided our solutions to be of two types: type I and type II. Type I is the most interesting solution as it is derived from a recursion operator within the conservation law of the KdV and mKdV infinite hierarchy. In both solution types, the knowledge of association between symmetry and conserved components was exploited and formed the basis of our approach to reduce the order of the partial differential hierarchy to an ordinary differential hierarchy. Consequently, our method has many significant uses and can be extended to solve other infinite hierarchies. Specifically, it may be applied to any hierarchy in possession of a recursion operator, for example, the Kaup–Kupershmidt hierarchy. Furthermore, it would be interesting to attempt a study of systems of nonlinear equations with known recursion operators, such as the Hirota–Satsuma system or the nonlinear Schrödinger system of real equations.An advantage of our approach is that it can easily be implemented into computer algebra programs such as Maple or Mathematica. A disadvantage is that, at higher-order members of the hierarchy, one may struggle to solve the reduced conservation laws, simply because the computations are too involved and computer algebra programs may run out of memory to complete the necessary calculations.In the known literature, there are numerous methods to solve members of the KdV hierarchy, for example, the generalized Kudryashov method [28], the double Laplace transform [29, 30], the differential transform method [31], the tanh-expansion method [32], the exp-function method [33], and the G′/G-expansion method [34]. Our approach involves a recursion operator and conservation law to aid the analysis of the nonlinear partial differential hierarchy. To the best of our knowledge, this is the first time that a study has conceived an approach for dealing with the entire KdV hierarchy.



---

*Source: 5378853-2022-01-18.xml*