Investigation of optical solitons to the nonlinear complex Kundu-Eckhaus and Zakharov-Kuznetsov-Benjamin-Bona-Mahony equations in conformable.

This research manuscript focuses on the extraction of dark-bright solitons and periodic wave distributions of two models, namely, the Zakharov-Kuznetsov-Benjamin-Bona-Mahony equation and complex Kundu-Eckhaus equation with conformable derivative. To reach these important properties, the generalized exponential rational function method is considered. To observe wave distributions in periodic and singular sense, dynamical behaviour modulus of solutions are also archived. Strain conditions for validity of results obtained in this paper are also reported.

Introduction

Soliton theory is one of the most studied areas in applied sciences, in particular, nonlinear dynamical media. In this sense, Divo-Matos et. al. have firstly investigated the new model for gas adsorption isotherm at high pressures (Divo-Matos et al. 2021). This observation comes from the hydrogen which has been considered as an ideal energy source in the future. They have observed the new isotherm model for supercritical gas adsorption by using the Redlich-Kwong’s equation of state at high pressures. Another a new model to estimate permeability has been formed to explain the carbonate and shale samples (Liu et al. 2020). They have calculated the curvatures from the mercury saturation data by using the Fermic-Dirac function. Gunasekeran and his team introduced the digital health during COVID-19 (Gunasekeran et al. 2021). As it is well known by almost everybody, COVID-19 has affected many peoples from all over the world. In these works, they presented the mathematical perspective to the Covid-19 disease by giving deeper properties such as dynamical structures (Gao et al. 2020), modeling the dynamics of novel coronavirus (2019-nCov) (Khan and Atangana 2020), stochastic COVID-19 Lvy jump model with isolation strategy (Danane et al. 2021), unreported cases of 2019-nCOV epidemic outbreaks (Gao et al. 2020), threshold condition and non pharmaceutical interventions control strategies (Zamir et al. 2021), the spread of COVID-19 with new fractal-fractional operators (Atangana 2020), the total digitalization in education during the Covid pandemi (Yuce 2022). Moreover, many real world problems have been also studied by using mathematical methods and tools such as sine-Gordon expansion method (Eskitascioglu et al. 2019), tanh and extended tanh methods (Bibi and Mohyud-Din 2014; Zayed and Tala-Tebue 2018; Fan 2000). Some important properties of reduced higher-order and perturbed nonlinear Schrödinger equation (Kudryashov 2021, 2020) and so on were presented (Wazwaz 2018; Kaplan et al. 2017; Cattani 2018; Zhang et al. 2015; Bridges and Ratliff 2017a, b). Therefore, such nonlinear evaluation equations (NLEs) have always been one of the most inspiring tools for researchers to study real phenomena and models arising in mathematics, physics, engineering and various fields of science (Chen and Ren 2022; Nabti and Ghanbari 2021; Wang et al. 2020; Ghanbari and Kumar 2021; Hu et al. 2022; Akbulut and Kaplan 2018; Durur and Yokus 2021; Sulaiman et al. 2021; Gao et al. 2020; Ghanbari and Djilali 2020; Rajesh Kanna et al. 2020; Ghanbari 2020; Kaplan and Akbulut 2021; Ghanbari and Atangana 2020; McCue et al. 2021; Djilali and Ghanbari 2021; Gao et al. 2020; Ghanbari 2020a, b; Srivastava et al. 2019; Erturk and Kumar 2020; Ghanbari 2021; Saouli 2020; Ghanbari 2021; Munusamy et al. 2020; Ghanbari et al. 2020; Kaplan 2017; Ghanbari et al. 2019; Kudryashov 2020; Djilali and Ghanbari 2021), (Ali Akbar et al. 2021; Eslami et al. 2021; Rezazadeh et al. 2019; Akinyemi et al. 2021; Hosseini et al. 2020; Nisar et al. 2022; Pinar et al. 2020; Khodadad et al. 2021; Ozkan et al. 2021; Halidou et al. 2022; Can et al. 2020).

The remaining parts of research manuscript are organized as follows. In Sect. 2, overview of the conformable derivative and its general properties are presented. Section 3 provides the general properties of the projected scheme, namely, generalized exponential rational function method GERFM. Main results showing that the system can exhibit some interesting properties of the governing models such as Zakharov–Kuznetsov–Benjamin–Bona–Mahony (ZKBBM) equation and complex Kundu–Eckhaus equation with conformable derivative are given in Sect. 4. The meaning of the parameters of the new findings is given there in terms of different simulations. We discuss and conclude our results in the last section of the paper.

Firstly, ZKBMM equation (lzaidy 2013; Khodadad et al. 2017; Aksoy et al. 2016) given aswhere a, b are real-valued constants is studied.

Secondly, nonlinear complex Kundu–Eckhaus equation defined by Khater et al. (2018), Khodadad et al. (2017), Mirzazadeh et al. (2018):where are real-valued constants is investigated by using projected method to extract some important properties.

Overview of the conformable derivative

Recently, an operator called conformable was formulated by Khalil in Khalil et al. (2014). The conformable calculus satisfies all the properties of the standard calculus. This operator may be considered to be a natural extension of the classical properties (Usta 2018; Abdeljawad 2015; Benkhettou et al. 2016; Ünal and Gökdogan 2017).

Definition 2.1

Let , the conformable derivative of a function f(t) of order , is defined as (Khalil et al. 2014)

This new definition satisfies the properties in the following theorem:

Theorem 1

(Khalil et al. 2014) Let , f, g be -differentiable at a point t, then

, for .

, for .

.

.

If, in addition, f is differentiable, then . (Abdeljawad 2015) established the chain rule for conformable derivatives as following theorem:

Theorem 2

Let , be a function such that f is differentiable and also -conformable differentiable. Let g be a differentiable function defined in the range of f. Thenwhere prime denotes the classical derivatives with respect to t.

Overview of the generalized exponential rational function method

In this section, we state the main steps of GERFM as follows (Biazar and Ghanbari 2012; Ghanbari and Kuo 2021; Kumar et al. 2021; Ghanbari and Akgül 2020; Ismael et al. 2020; Ghanbari et al. 2020)

Let us take into account the following nonlinear partial differential equation Via the transformations and , in nonlinear partial differential equation (4), we attain which is indeed an ordinary differential equation; where the values of and will be found later.

Consider equation (5) has the solution of the form where The values of constants , and are determined, in such a way that solution (6) always persuade equation (5). By considering the homogenous balance principle the value of M is determined.

Putting equation (6) into Eq. (5) we give the following algebraic equation , in terms of for . After setting each of the coefficients of variables in to zero, a system of nonlinear equations in terms of these parameters is constructed.

By solving the above system of equations using any symbolic computation software, the values of , , and are determined, replacing these values in Eq. (6) provides us the soliton solutions of Equation (4).

Main results

In this section, we apply GERFM into the governing models to find some important properties.

The ZKBMM equation with conformable derivative

In this subsection, we investigate exact wave solutions in studying the conformable version of the ZKBMM Eq. (1), giving byBy considering the wave transformationwhere c, k are non zero constants. Utilizing the wave transformation (33), we getWhen we integrate once and setting the integration constant to zero in (10), we obtainApplying the balance principle on the terms and in Eq. (11), we have . Consequently, we find . Using Eq. (7) together with , we haveThe following results are provided by processing the general steps required in the method.

Set 1:

One obtains and , so Eq. (7) turns toCase 1.1: We also obtainPutting values in Equations (12) and (13), yields the following solutionConsequently, we get the solution of Equation (8) asFigure 1 shows the dynamic properties of for , and for two different values of ’s.

Fig. 1

Dynamic behaviours modulus of solution for

Case 1.2: We also obtainPutting values in Equations (12) and (13), yields the following solutionConsequently, we get the solution of Equation (8) asFigure 2 shows the dynamic properties of for , and for two different values of ’s.

Fig. 2

Dynamic behaviours modulus of solution for

Set 2:

One obtains and , so Eq. (7) turns toCase 2.1: We also obtainPutting values in Equations (12) and (16), yields the following solutionConsequently, we get the solution of Equation (8) asFigure 3 shows the dynamic properties of for , and for two different values of ’s.

Fig. 3

Dynamic behaviours modulus of solution for

Case 2.2: We also obtainPutting values in Equations (12) and (16), yields the following solutionConsequently, we get the solution of Equation (8) asFigure 4 shows the dynamic properties of for , and for two different values of ’s.

Fig. 4

Dynamic behaviours modulus of solution for

Set 3:

One obtains and , so Eq. (7) turns toCase 3.1: We also obtainPutting values in Equations (12) and (19), yields the following solutionConsequently, we get the solution of Equation (8) asCase 3.2: We also obtainPutting values in Equations (12) and (19), yields the following solutionConsequently, we get the solution of Equation (8) asFigure 5 shows the dynamic properties of for , and for two different values of ’s.

Fig. 5

Dynamic behaviours modulus of solution for

Set 4:

One obtains and , so Eq. (7) turns toWe also obtainPutting values in Equations (12) and (22), yields the following solutionConsequently, we get the solution of Equation (8) asFigure 6 shows the dynamic properties of for , and for two different values of ’s.

Fig. 6

Dynamic behaviours modulus of solution for

Set 5:

One obtains and , so Eq. (7) turns toWe also obtainPutting values in Equations (12) and (24), yields the following solutionConsequently, we get the solution of Equation (8) asFigure 7 shows the dynamic properties of for , and for two different values of ’s.

Fig. 7

Dynamic behaviours modulus of solution for

Set 6:

One obtains and , so Eq. (7) turns toWe also obtainPutting values in Equations (12) and (26), yields the following solutionConsequently, we get the solution of Equation (8) asSet 7:

One obtains and , so Eq. (7) turns toWe also obtainPutting values in equations (12) and (28), yields the following solutionConsequently, we get the solution of Equation (8) asFigure 8 shows the dynamic properties of for , and for two different values of .

Fig. 8

Dynamic behaviours modulus of solution for

Set 8:

One obtains and , so Eq. (7) turns toWe also obtainPutting values in Equations (12) and (30), yields the following solutionConsequently, we get the solution of Equation (8) as

The conformable Kundu–Eckhaus equation

In this part, we aim to construct exact wave solutions in studying the conformable version of the Kundu–Eckhaus equation (32), giving byIn order to find solutions of Equation (32), following new variables are introducedwhere and are constants to be determined.

Taking Eq. (33) into account in Eq. (32) yieldsIf we balance the highest derivative term of and nonlinear term of in Eq. (34) as , we obtain .

So, we need to use a new transformation in Eq. (34) to getNow if we apply the balance principle on the terms and in Eq. (35), we have , so . Using Eq. (7) together with , we haveThe following results are provided by processing the general steps required in the method.

Set 1:

One obtains and , so Eq. (7) turns toCase 1.1: In this case we also obtainPutting values in Equations (36) and (37), yields the following solutionandFigure 9 shows the dynamic behavior of solution for , and .

Fig. 9

Dynamic behaviours modulus of solution , real part (left) and imaginary part (right) for , and

Case 1.2: In this case we also obtainPutting values in Equations (36) and (37), yields the following solutionandFigure 10 shows the dynamic behavior of solution for , and .

Fig. 10

Dynamic behaviours modulus of solution , real part (left) and imaginary part (right) for , and

Set 2:

One obtains and , so Eq. (7) turns toCase 2.1: In this case we also obtainPutting values in equations(36) and (40), yields the following solutionandCase 2.2: In this case we also obtainPutting values in equations (36) and (40), yields the following solutionandSet 3:

One obtains and , so Eq. (7) turns toWe also obtainPutting values in equations(36) and (43), yields the following solutionandFigure 11 shows the dynamic behavior of solution for , and .

Fig. 11

Dynamic behaviours modulus of solution , real part (left) and imaginary part (right) for , and

Set 4:

One obtains and , so Eq. (7) turns toWe also obtainPutting values in equations (36) and (45), yields the following solutionandSet 5:

One obtains and , so Eq. (7) turns toWe also obtainPutting values in equation (36) and (47), yields the following solutionandFigure 12 shows the dynamic behavior of solution for , and .

Fig. 12

Dynamic behaviours modulus of solution , real part (left) and imaginary part (right) for , and

Set 6:

One obtains and , so Eq. (7) turns toWe also obtainPutting values in equations (36) and (49), yields the following solutionandFigure 13 shows the dynamic behavior of solution for , and .

Fig. 13

Dynamic behaviours modulus of solution , real part (left) and imaginary part (right) for , and

Conclusions

In this research, the projected method has been effectively applied to the nonlinear complex Kundu–Eckhaus and Zakharov–Kuznetsov–Benjamin–Bona–Mahony equations in conformable domain used to explain the most fascinating problems of modern optics. Some important optical soliton solutions such as single (dark, bright and singular), complex solitons, as well as a hyperbolic, travelling wave and trigonometric function solutions have been successfully extracted. The graphical simulations of the reported solutions have been also presented in Figs. (1, 2, 3, 4, 5, 6, 7, 8, 9, 10 , 11, 12, 13). While some of these figures symbolize singular wave properties, others gives travelling wave distributions. The results are entirely new, interesting and play an important roles in the field of the nonlinear Schrödinger equation because the studied model, namely nonlinear complex Kundu–Eckhaus equation is one of the part of NLSE. When we consider the obtained results, it is clear that the method has less limitations than other methods in determining the exact solutions of the equations. Despite the simplicity and ease of use of this method, it has a very powerful performance and is able to introduce a wide range of various types of solutions to such mathematical models. The idea used in this paper is readily applicable to solving other partial differential equations in mathematical physics. Finally, we observed that the propagation dynamics of these solutions obtained in this paper via GERFM may be used to explain the general properties of the nonlinear optical wave distributions (Weisstein 2002).