Soliton Solutions of Mathematical Physics Models Using the Exponential Function Technique

Shumaila Javeed 1,*,† , Khurram Saleem Alimgeer 2,*,† , Sidra Nawaz 1,†, Asif Waheed 3 , Muhammad Suleman 1, Dumitru Baleanu 4,5,6 and M. Atif 7 1 Department of Mathematics, COMSATS University Islambad, Islamabad Campus, Park Road, Chak Shahzad Islamabad 45550, Pakistan; sidranawaz25@gmail.com (S.N.); muhammad.suleman@comsats.edu.pk (M.S.) 2 Department of Electrical computer Engineering, COMSATS University Islambad, Islamabad Campus, Park Road, Chak Shahzad Islamabad 45550, Pakistan 3 Department of Mathematics, COMSATS University Islamabad, Attock Campus 43600, Pakistan; asifwaheed8@gmail.com 4 Department of Mathematics, Cankaya University, Ankara 06790, Turkey; dumitru@cankaya.edu.tr 5 Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40447, Taiwan 6 Department of Mathematical Physics, Institute of Space Sciences, Magurele-Bucharest 077125, Romania 7 Department of Physics and Astronomy, College of Science, King Saud University, Riyadh 11451, Saudi Arabia; atifphy@gmail.com * Correspondence: shumaila_javeed@comsats.edu.pk (S.J.); khurram_saleem@comsats.edu.pk (K.S.A.) † These authors contributed equally to this work.


Introduction
When a physical system is modeled mathematically, nonlinear evolution equations are generated. Many physical phenomenon are inherently nonlinear. Consequently, modeling such phenomena mathematically gives the nonlinear evolution equation. Investigation of the wave solutions relating to the nonlinear PDEs in imperative for contemplating nonlinear physical occasions.
Recently, scientists and engineers have been attracted towards PDEs of non-integer orders. Ordinary Calculus does not satisfy the needs of solving such PDEs with non-integer problems. These needs extended the boundary of Ordinary calculus for solving these PDEs. In this paper, three nonlinear PDEs, i.e., Burger's equation, Kortewegde vries (Kdv) equation and Zakharov-Kuznetsov (ZK) equation are discussed. All of these nonlinear PDEs have symmetries in many ways. For example, Burger's equation is invariant, i.e., symmetric with respect to any space or time shift. Burger's equation also possesses scale symmetries. The said symmetry is based on the a priori premise that the division of variables into dependent and independent ones is a necessary part of a PDE. However, many arguments show that this point of view is limited. In particular, numerous tricks are available to find transformations which do not respect this division. Alexandre Vinogradov discussed four ideas to show the symmetries in nonlinear PDEs [1]. Readers are referred to these papers for understanding symmetry, its various variant of definitions, and application of symmetries on PDEs etc. [2][3][4][5][6][7].
Due to the applications of these PDEs in variety of practical problems like hydrodynamics, biomathematics, fluid dynamics, chemical physics, nonlinear optics, plasma physics, elastic media, optical fibers, geochemistry and chemical kinematics, there is a need to have general solutions to fit in for every kind of practical problem. We can explore hidden physical information extracted from nonlinear phenomenon by solving a nonlinear equation formulated against a complex event. In order to find the exact solutions of nonlinear PDEs, many powerful methods have been developed. A few of the methods are: the sub-equation method [8][9][10], the modified simple equation method [11], the exp(−φ(ξ)) method [12,13], homotopy analysis method [14,15] , first integral method [16,17], G/G expansion method [18][19][20], the functional variable method [21,22], the improved tan φ 2 -expansion method [23,24], etc.
In 2006, He and Wu presented an exp-function method to find the travelling wave solutions of nonlinear PDEs [25][26][27][28]. This method was reliable and effective to find a solution for such complex physical phenomena. The exp-function method is based on the homogeneous balancing principle. This is the reason that it can be applicable to nonlinear PDEs that satisfy the homogeneous balancing principle. The performance of the exponential function method is satisfactory. Moreover, it is not only reliable but also effective. A concise and promising tool for solving nonlinear PDEs could be built with the help of the Exp-function method. In this paper, the exp-function method is employed to extract the new exact solutions of the conformable ZK equation, conformable KdV equation and conformable Burger's equation.
This paper is divided into the following sections. The subsequent Section comprises of knowledge about preliminaries related to the conformable derivative, the same section deals with the properties of the derivatives, Section 3 comprises the exp-function method, Section 4 consists of implementation of exp-function for ZK equation, KdV equation and Burger's equation and finally conclusions and recommendations are presented in Section 5.

Groundwork for Conformable Derivative
Many scientists likes to have simple, effective and the most natural definitions; these definitions are preferred and have become popular due to ease of understanding. In this regard, khalil et al. described the definition of local derivative also known as conformable derivative of order ζ ∈ (0, 1] [29], which became popular soon after publishing. For any value of ζ, the definition can be generalized, but the most important case is ζ ∈ (0, 1]. Once it is demonstrated for particular case of ζ, the other cases become simple. In this section, conformable derivatives are explained in detail with examples. A conformable derivative of s of order ζ is defined for any generalized function S where s : [0, ∞) → R , as follows for all t > 0, ζ ∈ (0, 1). If s is ζ-differentiable for the range (0, a), a > 0 and lim t→0 + s (ζ) (t) exists, then define It is obvious that a constant function has zero conformable derivative. Also, if M ζ s(t) = 0 then s(x) = 0 for some x ∈ (c, d) using conformable fractional mean value theorem proved in [29].
The conformable integral of the same function is defined in following equation: Here, c ≥ 0 and ζ ∈ (0, 1]. The above integral is in the form of the famous Riemann integral in improper format. For ζ ∈ (0, 1] if r is a function which is differentiable with respect to ζat a point t > 0 then (M ζ r)(t) = t 1−ζ dr dt (t).
This definition has gained the attention of many scientist and a lot of work has been done on it. It was introduced in 2014 then developed in 2015 and is currently under intensive investigations [29][30][31]. Abdel Jawad used this derivative to demonstrate exponential functions, by parts integration, chain rule, Taylor power series expansion [30] and Laplace transform.

Presentation of Exp-Function Method for Nonlinear Conformable Pdes
The primitive work on the exponential function method was explored and performed by He and Wu in 2006 [25]. This work was further studied with detailed analysis in [32][33][34][35][36][37][38]. Homogeneous balancing principle was the base of this proposed method. The only nonlinear conformable PDEs that can be solved by this method are those which satisfy the homogeneous balancing principle. As a result, the exact solutions can be obtained for these equations. In mathematical physics there are many types of nonlinear evolution equations which appear to present any practical phenomenon, and requires effective, concise and easy method to solve these equations. In this way, the Exp-function is an excellent candidate to be used as solver for these equations. Symmetries are involved in the transformations of space and time. Spatial symmetries are related to spatial geometry of a physical system, while spatio-temporal symmetries involve changes in both space and time. The spatio-temporal symmetries are used to transform PDEs into ODEs (see Equation (6)). These symmetries play a vital role in finding the exact solution. Various methods are available to reduce PDEs into ODEs using some special ansatz for obtaining travelling wave solutions [2]. Almost all of these methods can be viewed as a special case of a symmetry reduction, potentially known as a non-classical symmetry. This approach is used (cf. Equation (6)) in this paper to obtain the exact solutions of nonlinear PDEs. In this paper, wave transformation is used to obtain nonlinear ODEs from nonlinear PDEs and hence making it easy to find a solution.
A general nonlinear conformable PDE [25] is considered in the form of equation: A new wave variable is introduced in the following equation where q 0, r 0,and p 0 are arbitrary constants, and we will evaluate these constants later. By using chain rule, we get Putting Equation (7) into Equation (5), it becomes a nonlinear ODE as Suppose the solution of the above ODE is where all the integers listed in the equation are positive integers and the values of these integers are determined by homogeneous balancing principle, a i and b j are unknown constants. The like linear term of the lowest order in Equation (8) is balanced with lowest order linear term, in order to have the values of c and m used in the same equation. Likewise, the highest order linear term of Equation (8) is balanced with the highest order nonlinear term, in order to find out the values of d and n variables present in the same equation. This procedure will lead to an exact solution of nonlinear PDEs.

Exact Solutions of Burger'S Equation, Zk Equation and Kdv Equation
In this section, three equations are selected for having the exact solutions. These equations are, conformable ZK, KdV and Burger's equations. The solutions of these equations are discussed in the following sections

Exact Solutions of Conformable Zk Equation
In this section, conformable ZK equations will be presented and then the solution of this equations will be discussed in details. The conformable ZK equation is given as where D ζ t shows the conformable derivative for ζ ∈ (0, 1]. Presenting a transformation ω = px + qy + r t ζ ζ , and then by chain rule we get an ODE in the following form can be obtained by putting Equation (11) into Equation (10), The solution of Equation (12) will be considered with unknown variables in the form of The method to have the values of different variables in the above equation is described in the previous section. To evaluate the value of variables, a linear term of highest order of Equation (12) is compared with the highest order nonlinear term u (ω) = ... + k 1 e (d+7n)ω + ... ... + k 2 e (8nω) + ... , and u 2 u = ... + k 3 e (n+3d)ω + ... ... + k 4 e (4nω) + ... = ... + k 3 e (5n+3d)ω + ... ... + k 4 e (8nω) + ... , where k i are coefficients, obtained by highest order term balancing in Equations (14) and (15), we get Comparing the linear term of lowest order of Equation (12) with the lowest order nonlinear term can provide the values of m and c, hence we have l i are coefficients. Lowest order terms in Equations (16) and (17) Switching Equation (18) into Equation (12), we have where , h i (i = −3, ..., 0, ..., 3) are unknown constants that are obtained by software Maple 13.0.
Following system of algebraic equations can be obtained by equating coefficients of e iω equals to zero, Solving algebraic equations for this system, following equation can be obtained Hence, ZK equation has the following exact solution where p, b 1 , a 0 and q are real numbers.

Exact Solution of Conformable Kdv Equation
Consider conformable kdV equation as Introducing a transformation ω = px + r t ζ ζ , and then by chain rule, Equation (23) is reduced into an ODE of the form assuming the solution of Equation (24) in the following form using exp-function method, where d, c, m and n are unknowns. These constants are obtained by comparing linear term of highest order of Equation (24) and balancing it with the highest order nonlinear term, i.e., u ω with u 2 . As a result, we obtain d = n, and c = m.
Like previous subsection, following system of the algebraic equations can be obtained by equating coefficients of e iω equals zero, following result can be obtained by solving this system of algebraic equations, Therefore, exact solution can be obtained for time-fractional kdv equation of the following form

Exact Solution of Conformable Burgers Equation
The nonlinear conformable Burgers equation is as follows where ζ ∈ (0, 1), u(x, t) is the function of velocity having spatial dimensions t and x, ν > 0.0. ν > 0.0 is the kinematic fluid viscosity, and ∂ ζ u ∂t ζ shows the derivative of function u(x, t) in conformable form. Introducing a transformation ω = x − r t ζ ζ , then using the chain rule we have Thus Equation (31) is transformed into an ODE, which is of the form By our proposed method, Equation (33) has a solution of the form u(ω) = a d e dω + ... + a −c e −cω b n e nω + ...
where d, c, m and n are unknowns. These constants are obtained through the method described previously in the same section. The linear term of highest order of Equation (33) following system of algebraic equations can be obtained by using software Maple 13.0 after Substituting Equation (35) into Equation (33), following solution can be obtained by solving this system of algebraic equations, Substituting Equation (38) into Equation (35), we have

Conclusions
The aim of this work was to find the exact solutions of nonlinear conformable PDEs. ZK equation, KdV equation and Burgers equation are selected to present the validity of the Exp-function method for obtaining the exact solutions of these equations. The proposed software assisted method allows us to perform complicated nonlinear conformable algebraic calculations. The solving procedure described in this paper for nonlinear conformable PDEs reveals that the exp-function method is a succinct tool for solving such equations, and we expect it to be successful in the future.