A Study of (2+1)-Dimensional Konopelchenko-Dubrovsky (KD) System: Closed-Form Solutions, Solitary Waves, Bifurcation Analysis and Quasi-Periodic Solution


 Nonlinear evolution equations (NLEEs) are extensively used to establish the elementary propositions of natural circumstances. In this work, we study the Konopelchenko-Dubrovsky (KD) equation which depicts non-linear waves in mathematical physics with weak dispersion. The considered model is investigated using the combination of generalized exponential rational function (GERF) method and dynamical system method. The GERF method is utilized to generate closedform invariant solutions to the (2+1)-dimensional KD model in terms of trigonometric, hyperbolic, and exponential forms with the assistance of symbolic computations. Moreover, three-dimensional graphics are displayed to depict the behavior of obtained solitary wave solutions. The model is observed to have single and multiple soliton profiles, kink-wave profiles, and periodic oscillating nonlinear waves. These generated solutions have never been published in the literature. All the newly generated soliton solutions are checked by putting them back into the associated system with the soft computation via Wolfram Mathematica. Moreover, the system is converted into a planer dynamical system using a certain transformation and the analysis of bifurcation is examined. Furthermore, the quasi-periodic solution is investigated numerically for the perturbed system by inserting definite periodic forces into the considered model. With regard to the parameter of the perturbed model, two-dimensional and three-dimensional phase portraits are plotted.


Introduction
Nonlinear evolution equations (NLEEs) have extensive significance in the area of applied mathematics and physics. Finding the exact solutions for NLEE is an essential task as NLEE describes numerous phenomenon in nonlinear dynamics, engineering, optical fibre, plasma physics, fluid mechanics, natural sciences, etc. A large number of researchers and mathematicians have developed various effective techniques for computing exact solutions of NLPDEs (nonlinear partial differential equations), for instance, tanh function method [1], Hirota's bilinear method [2,3], the Jacobi elliptic function expansion method [4], Kudryashov method [5], the G ′ G -expansion method [6], Darboux transformation method [7], the Backlund transformation method [8], the inverse scattering method [9], Lie-symmetry analysis [10], multiple exp-function method, and many others. Among these techniques, GERF method [11][12][13][14] is very effective, robust and straightforward approach for finding the abundant exact soliton-form solutions of various NLPDEs.
The dynamics of NLPDEs grants us to understand and predict the acceptable structures of the associated complex nonlinear systems. A soliton or solitary wave is the particle-like object with the finite energy and amplitude, which save its form during propagation and restore it after the collision with another solitons. Nowadays, as a consequence, it is a very hot subject matter to derive the exact closed form solutions of NLPDEs. The soliton-form solutions of such type of NLPDEs are extensively favorable in the various areas such as nonlinear sciences, mathematical physics, plasma physics, applied mathematics, engineering, applied sciences and nonlinear dynamics.
In the recent study, we will analyze the following (2+1)-dimensional Konopelchenko-Dubrovsky (KD) model where u and v are the differentiable functions with respect to x, y and t variables. Here, a and b are arbitrary constants. Many reseachers have used some productive techniques to investigate the exact analytical solutions of the KD system. Konopelchenko and Dubrovsky [15] derived (2+1)-dimensional KD (Konopelchenko-Dubrovsky) equation in 1984. In 2014 Kumar et. al [16] used Lie symmetry analysis as well as travelling wave hypothesis to extract solutions to KD equations. In 2019 Alfalqi et al. [17] applied the modified simplest equation method and B-spline method to KD-equation. Recently, Khater et al. [18] implement modified auxiliary equation technique to this system to find analytical traveling wave solutions. Ren et al. [19] in 2016 obtained the non-local symmetries for the K-D equation with the truncated Painleve method and the Mobius (conformal invariant form.
Recently, the investigation of differential equations through the bifurcation analysis has become an important topic in the field of research. Bifurcation is a rapid quantitative shift in the model with a gentle change in the values of parameters. The exploration of the dynamics of nonlinear periodic forms is a significant part for the investigation of the physical propositions in detail. For example, the occurrence of homoclinic orbits, smooth heteroclinic orbits and periodic orbits for travelling wave models describes the periodic wave solutions, oscillatory travelling wave solutions, smooth kink wave solutions and smooth solitary solutions for considered PDEs, respectively. One can refer to [20][21][22][23][24] for the study of recent work in this field.
In this research article, the study of (2 + 1)-dimensional Konopelchenko-Dubrovsky (KD) model (1) is investigated using GERF (generalized exponential rational function) technique. The primary goal of this research study is to construct abundant exact analytical closed-form solutions for the system (1). Moreover, we exhibit the dynamics of solitary wave profiles of some soliton solutions in three dimensional and two dimensional graphics in numerical simulations, and hence we believe that the evolutionary profile dynamics of generated exact closed-form solutions are very impressive and advantageous for physical phenomena. The obtained soliton solutions are entirely new, novel and have never been reported in the previous findings. Also, we have investigated the dynamical system (1) using the bifurcation analysis. Numerically, the quasi-periodic solution is investigated for the perturbed system by including definite forces to the considered model, and two-dimensional and three-dimensional phase portraits are also plotted for the perturbed system.
The strategy of this article is organised as follows: Section 2 deals with the introduction and methodology of GERF approach is presented. In section 3, we find exact travelling wave solution of KD equations. This section also includes some particular three-dimensional shapes which provide more explanation to the behavior of these generated solutions. Graphically, multi-soliton structures, kink-wave structures, and the interaction of multi-soliton and kink wave solution have been observed for some soliton solutions. Section 4 deals with the bifurcation analysis of the dynamical system of eq. (1), and relative phase portraits are plotted for the considered system. Section 5 is related to the investigation of quasi-periodic solution for the perturbed system by inserting perturbation term to the associated model (1). Finally, the conclusion is given at the end.

2
Methodology of GERF method This GERF method was introduced by Ghanbari and Inc [11]. We will provide an explanation of GERF technique step-wise in this section: • Let us consider the system of two nonlinear PDEs including three variables x, y and t given as Here, u(x, y, t) = U (ξ) and v(x, y, t) = V (ξ) are the wave transformations where ξ = αx+βy −µt is used to obtain a system of ODEs (ordinary differential equations) which further provides an ordinary differential as where U = U (ξ), U ′ = dU dξ and α, β and µ are constants, to be calculated later. • Assume that the solitary wave solution of (3) can be represented as where Ψ(ξ) = r 1 e s 1 ξ + r 2 e s 2 ξ r 3 e s 3 ξ + r 4 e s 4 ξ , and r i , s i (1 ≤ i ≤ 4), R 0 , R k and S k (1 ≤ k ≤ M ) are constants to be calculated and M (a positive integer) is evaluated by homogeneous balancing method.
Equalizing the coefficients of A to zero, yields a set of algebraic equations.
• With the assistance of Mathematica, we find the values of the coefficients R 0 , R k , S k , α, β and µ by solving above mentioned algebraic equations.
• Inserting all these values into (4) we get the solution of determining equation (3). Accordingly we can get other function and hence we achieve some new solitary wave solutions of (2).

GERF method for the (2+1)-D Konopelchenko-Dubrovsky(KD) model
Making use of wave transformation u(x, y, t) = U (ξ) and v(x, y, t) = V (ξ) where ξ = αx + βy − µt in eq.(1), we obtain the following system of ODEs Integrating the first equation we find, Putting it into the second equation of (7) and integrating the resulting equation, and neglecting the constant of integration, the following ODE is obtained as Using balancing principle on terms U (3) and U ′′ of the equation (9), we get M + 2 = 3M which yields M = 1. Employing M = 1 in (4), we obtain the trial solution given as By using equation (10) into (9) and according to the step of GERF method, the following cases are considered for finding the exact solitary wave solutions of (1) with the assistance of Mathematica.
With the assistance of soft computation via Mathematica, we solve algebraic equation for obtaining the values of parameters and hence following set of solutions can be achieved. Solution set 1.1: Substituting the values of above known-constant parameters into eq.(10) and plugging it with eq. (9), we obtain the expression for U as Hence, from eq.(8), we obtain the expression for V as Accordingly, the solution of (1) is obtained as Substituting the values of above known-constants into eq.(10) and plugging it with eq. (9), we obtain the expression for U as Hence, from eq.(8), we obtain the expression for V as Solution set 1.3: Substituting the values of known-constants mentioned above into eq. (10) and plugging it with eq. (9), we obtain the value for U as Hence, from eq. (8), we obtain the expression for V as Accordingly, the solution of (1) is obtained as With the assistance of soft computation via Mathematica, we solve algebraic equation for obtaining the values of parameters and hence following set of solutions can be achieved. Solution set 2.1: We substitute the values of known-constants mentioned above and expression (21) into eq.(10), and plugging it with eq. (9), we obtain the value for U as Hence, from eq.(8), we obtain the expression for V as Accordingly, the solution of eq. (1) is obtained as Solution set 2.2: We substitute the values of known-constants mentioned above and expression (21) into eq.(10), and plugging it with eq. (9), we obtain the value for U as Hence, from eq.(8), we obtain the expression for V as Accordingly, the solution of (1) is obtained as Solution set 2.3: We substitute the values of known-constants described above and expression (21) into eq.(10), and plugging it with eq. (9), we obtain the expression for U as Hence, from eq.(8), we obtain the expression for V as Accordingly, the solution of (1) is obtained as With the assistance of soft computation via Mathematica, we solve algebraic equation for obtaining the values of parameters and hence following set of solutions can be achieved. Solution set 3.1: We substitute the values of known-constants described above and expression (31) into eq.(10), and plugging it with eq. (9), we obtain the expression for U as Hence, from eq.(8), we obtain the expression for V as .
Accordingly, the solution of eq. (1) is obtained as Solution set 3.2: We substitute the values of known-constants mentioned above and expression (31) into eq.(10), and plugging it with eq. (9), we obtain the expression for U as Hence, from eq.(8), we obtain the expression for V as Accordingly, the solution of (1) is obtained as Solution set 3.3: We substitute the values of known-constants mentioned above and expression (31) into eq.(10), and plugging it with eq. (9), we obtain the expression for U as Hence, from eq.(8), we obtain the expression for V as We substitute the values of known-constants described above and expression (41) into eq.(10), and plugging it with eq. (9), we obtain the expression for U as Hence, from eq.(8), we obtain the expression for V as Accordingly, the solution of (1) is obtained as u(x, y, t) = 6α tan 3 12t(6a 2 α 3 +αb 2 ) Solution set 4.2: We substitute the values of known-constants mentioned above and expression (41) into eq.(10), and plugging it with eq. (9), we obtain the expression for U as Hence, from eq.(8), we obtain the expression for V as Accordingly, the solution of (1) is obtained as Solution set 4.3: We substitute the values of known-constants described above and expression (41) into eq.(10), and plugging it with eq. (9), we obtain the expression for U as Hence, from eq.(8), we obtain the expression for V as Accordingly, the solution of (1) is obtained as (50) We substitute the values of known-constants mentioned above and expression (51) into eq.(10), and plugging it with eq. (9), we obtain the expression for U as Hence, from eq.(8), we obtain the expression for V as Accordingly, the solution of (1) is obtained as Solution set 5.2: Hence, from eq.(8), we obtain the expression for V as Accordingly, the solution of (1) is obtained as Solution set 5.3: We substitute the values of known-constants mentioned above and expression (51) into eq. (10), and plugging it with eq. (9), we obtain the expression for U as Hence, from eq.(8), we obtain the expression for V as Accordingly, the solution of (1) is obtained as

Bifurcation analysis
We explore the new dynamics of KD system (1), in this section, by utilizing the concepts of bifurcation theory.
After simplification, equation (9) can be rewritten in the following form of planar dynamical system as where, The three equilibrium points for the above system of differential equations are computed as (0, 0), (U 1 , 0) and (U 2 , 0) on U-axis, where is the coefficient matrix of the linearized system of (61) at equilibrium point (U i , 0); J and T be the determinant and trace of the matrix M , respectively. Here, By using the theory of planar dynamical system [25], we can discuss following definitions for the critical points (U i , 0).
(i) When J < 0, then (U i , 0) is a saddle point.
(ii) When J > 0 and T 2 − 4J ≥ 0, then (U i , 0) is a node; which is stable if T < 0 and unstable if T > 0.
(v) When J = 0 and Poincare index of (U i , 0) is zero, then it is called the zero point.
For different choices of parameters A 1 , A 2 and A 3 , various cases in detail are explained as:

Quasi-periodic solution
The exploration of different dynamics of the perturbed system of the main system (1) by adding perturbation term ρ 1 cos(ση) is investigated in this section. Thus, equation (7) with the addition of perturbation term becomes where ρ 1 is constant representing the intensity of the perturbed term and σ represents the frequency.
To analyze the chaotic and periodic behavior of system (7) with the existence of a perturbation term, we will fix the influence of force and frequency of disruption and retain the physical parameters of the system under observation.

Conclusion
In summary, we investigated the (2+1)-dimensional Konopelchenko-Dubrovsky (KD) model and constructed numerous exact closed-form solutions using GERF (generalized exponential rational function) approach with soft symbolic computations via Mathematica. The established soliton solutions exhibit that the KD equations (1) admits abundant exact closed-form solutions having arbitrary constant parameters. The findings depict rich dynamical formations of the generated closed-form solutions in the forms of multi-solitons, kink-wave profiles, the interaction of multi-soliton and kink-wave solutions, and multiple solitons in figures 1, 2, 3 and 4. The obtained solutions will be beneficial in the theory of solitons, nonlinear dynamics, fluid mechanics, applied physics, optic fiber, natural sciences, physics, and many other areas. The technique we have used is one of the powerful techniques for obtaining the exact analytical solutions of NLPDEs as this GERF method constructs an extensive scale of established solitary wave solutions to the associated system. This method is very effective, trustworthy, and efficient. Moreover, the dynamics of KD equations (1) is also examined by using the bifurcation analysis in which we have found various equilibrium points, and different types of phase portraits have been discussed, see figures 5 and 6. Furthermore, we numerically studied the existence of quasi-periodic solutions for the perturbed model, which is obtained after inserting the periodic forces to the considered KD system (1). Also, 2D and 3D phase portraits shown by figures 7, 8 and 9 for the perturbed system were also plotted.

Availability of data and material
Data sharing is not applicable to this article as no data sets were created or analyzed in this study.