A Dynamical Analysis and New Traveling Wave Solution of the Fractional Coupled Konopelchenko–Dubrovsky Model

: The main object of this paper is to study the traveling wave solutions of the fractional coupled Konopelchenko–Dubrovsky model by using the complete discriminant system method of polynomials. Firstly, the fractional coupled Konopelchenko–Dubrovsky model is simplified into nonlinear ordinary differential equations by using the traveling wave transformation. Secondly, the trigonometric function solutions, rational function solutions, solitary wave solutions and the elliptic function solutions of the fractional coupled Konopelchenko–Dubrovsky model are derived by means of the polynomial complete discriminant system method. Moreover, a two-dimensional phase portrait is drawn. Finally, a 3D-diagram and a 2D-diagram of the fractional coupled Konopelchenko– Dubrovsky model are plotted in Maple 2022 software.


Introduction
In the present era, nonlinear evolution equations (NLEEs) [1][2][3][4] are employed in a number of areas like physics, chemistry, biology, fluid dynamics, engineering, optical fibers, plasma, and hydrodynamics.The analytical solutions of NLEEs can be applied to control complex behavior and difficult phenomena when the system displays [5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23].While there is no unified method to obtain the exact solution of nonlinear evolution equations, most of the time, NLEEs can be converted into an ordinary differential equation by taken the traveling wave transformation.Based on the efforts of many predecessors, various methods have been imposed to solve NLEEs precisely and analyze various wave phenomena.He and Wu [24] proposed the first-time Exp-function method to seek solitary solutions, periodic solutions and compacton-like solutions of the KdV equation and Dodd-Bullough-Mikhailov equation.By using Hirota's bilinear transformation method, Ma proved the existence of N-soliton solutions of the (2+1)-dimensional KdV equation, the Kadomtsev-Petviashvili equation, the (2+1)-dimensional Hirota-Satsuma-Ito equation, and a combined pKP-BKP equation [25,26], respectively.Li et al. presented the (w/g)-expansion method [27].Later, Arafat and his collaborators applied the customized (w/g)-expansion method to found the optical soliton solutions of the paraxial nonlinear Schrödinger equation and fractional Biswas-Arshed model [28,29], respectively.Wazwaz has derived the solitons and periodic solution of the Dodd-Bullough-Mikhailov and the Tzitzeica-Dodd-Bullough equations via the tanh-function approach [30].Seadawy and Iqbal analyzed the nonlinear damped Korteweg-de Vries equation in an unmagnetized collisional dusty plasma via the direct algebraic approach [31].Arafat also investigated scores of broad-spectral soliton solutions to the stated system via the auxiliary equation technique [32].However, research on the traveling wave solutions of more complex fractional order NLEEs is still ongoing, and there are still a large number of open issues that need to be addressed by researchers.
In this paper, we consider the fractional coupled Konopelchenko-Dubrovsky model [33] where D κ t u is the conformable fractional derivative.u = u(x, y, t) and v = v(x, y, t) represent the velocity components along the horizontal and vertical axes, respectively.a and b stand for the amplitude of the wave.When κ = 1, Equation (1) become the wellknown integer-order Konopelchenko-Dubrovsky model [34].The main object of this paper is to study the traveling wave solution of the fractional coupled Konopelchenko-Dubrovsky model by using the complete discriminant system method of polynomials.On the one hand, the main effort of this article is to focus on constructing the traveling wave solution of Equation (1).On the other hand, without solving Equation (1), its dynamic branch will be analyzed.
The conformable fractional derivative was first proposed by Khalil et al. [35].Compared with traditional fractional derivatives, the conformable fractional derivative has a more intuitive physical meaning.At present, it has been widely used in the construction of infectious disease dynamics models, nonlinear system modeling, and thermal science fields.Its definition is usually described as follows.
Definition 1 ([36]).Let f : [0, ∞) → R .Then, the conformable derivative of f of order κ is defined as and the function f is κ-conformable differentiable at a point t if the limit in Equation ( 2) exists.
The remaining sections of this article are arranged as follows: In Section 2, the traveling wave solutions of Equation (1) are constructed by using the complete discriminant system method.Moreover, a two-dimensional phase portrait is drawn.In Section 3, the threedimensional, two-dimensional, and density plots to some obtained solutions of Equation (1) are plotted.Finally, a brief summary is presented.

Traveling Wave Transformation
In this section, we first consider the wave transformation Substituting Equation (3) into Equation (1), we have Integrating the second equation of Equation ( 4), we obtain Substituting Equation (5) into the first equation of Equation (4), we have where d 1 is the integral constant.

Dynamical Analysis
Here, we consider the planar dynamic system of Equation ( 6) when where ).The first integration of Equation ( 7) is By setting the parameter values of fixed Equation ( 7), we draw the planar phase portrait of Equation ( 7), as shown in Figure 1.

Dynamical Analysis
Here, we consider the planar dynamic system of Equation ( 6) when d 1 " 0 where The first integration of Equation ( 7) is By setting the parameter values of fixed Equation ( 7), we draw the planar phase portrait of Equation ( 7), as shown in Figure 1.

Traveling Wave Solutions of Equation (1)
Multiplying both sides of Equation ( 6) by U 1 simultaneously and integrating it yields where here, d 2 is the integral constant.Here, we make a transformation: Substituting Equation (10) into Equation ( 9), we obtain: where ´b1 b 3 4b 4 `b0 .By integrating Equation (11) once, one has: where χ 0 denotes an integrating constant.Denote that Gpwq " w 4 `c2 w 2 `c1 w `c0 .We derive its complete discrimination system in the following form:

Traveling Wave Solutions of Equation (1)
Multiplying both sides of Equation ( 6) by U ′ simultaneously and integrating it yields where Here, we make a transformation: Substituting Equation (10) into Equation ( 9), we obtain: where (11) once, one has: where χ 0 denotes an integrating constant.Denote that G(w) = w 4 + c 2 w 2 + c 1 w + c 0 .We derive its complete discrimination system in the following form: According to the root-classifications of Equation ( 11), we will discuss the traveling wave solutions of Equation ( 1) under nine cases.
Case 1: By combining Equation (9) with Equation (3), the traveling wave solutions of Equation ( 1) can be expressed as below: Case 2: In this situation, the traveling wave solutions of Equation ( 1) take the form: Case 3: 2 , where α, β are real number, and α > β.
When w > α and w > β, or w < α and w < β, the traveling wave solutions of Equation ( 1) take the form: Case 6: The traveling wave solutions of Equation ( 1) take the form: where Case 7: When w > α 1 or w < α 4 , the traveling wave solutions of Equation ( 1) take the form: in which The traveling wave solutions of Equation (1) take the form: in which Case 9: , where l 1 , l 2 , s 1 , s 2 are real and s 1 ⩾ s 2 > 0. The traveling wave solutions of Equation ( 1) take the form: in which e 1 = l 1 e 3 + s 1 e 4 , e 2 = l 1 e 4 − s 1 e 3 , 4 ) e 2 3 +e 2 4 .

Numerical Simulation
In order to understand the dynamical processes and mechanisms of complex phenomena of the fractional coupled Konopelchenko-Dubrovsky model, numerical simulations of the obtained soliton solutions are given in this section.As is vividly shown in Figures 2a, 3a and 4a, u 1 (x, y, t), u 2 (x, y, t) and u 3 (x, y, t) stand for the tangent function solutions, the rational function solutions and the hyperbolic function solutions, respectively.Figures 2b, 3b and 4b denote the level curve at time t = 1.Furthermore, Figures 2c, 3c

Conclusions
In this article, we have presented the traveling wave solutions of Equation (1) via the complete discriminant system method, which is one of the most useful tools in solving NLEEs.The trigonometric function solutions, the rational function solutions, the hyperbolic function solutions, the exponential function solutions, and the Jacobian elliptic function solutions are obtained.Lastly, in order to understand the mechanisms of physical phenomena for Equation (1), we have also depicted two-dimensional and three-dimensional diagrams.In future work, we will focus on the traveling wave solutions and dynamic behavior of more complex NLEEs.Furthermore, we will also use the Darboux transformations to discuss the N-soliton solutions of more complex NLEEs.
Author Contributions: Software, Z.L.; Writing-original draft, J.W. All authors have read and agreed to the published version of the manuscript.

Conclusions
In this article, we have presented the traveling wave solutions of Equation (1) via the complete discriminant system method, which is one of the most useful tools in solving NLEEs.The trigonometric function solutions, the rational function solutions, the hyperbolic function solutions, the exponential function solutions, and the Jacobian elliptic function solutions are obtained.Lastly, in order to understand the mechanisms of physical phenomena for Equation (1), we have also depicted two-dimensional and three-dimensional diagrams.In future work, we will focus on the traveling wave solutions and dynamic behavior of more complex NLEEs.Furthermore, we will also use the Darboux transformations to discuss the N-soliton solutions of more complex NLEEs.
and 4c represent the density plots.Figures2d, 3d and 4dstand for the contour plots.

Fractal
Three-dimension graphic (b) Two-dimension graphic (c) Density plot (d) Contour plot
The authors declare no conflicts of interest.