Abstract

In this paper, the subequation method and the sine-cosine method are improved to give a set of traveling wave solutions for the time-fractional generalized Fitzhugh–Nagumo equation with time-dependent coefficients involving the conformable fractional derivative. Various structures of solutions such as the hyperbolic function solutions, the trigonometric function solutions, and the rational solutions are constructed. These solutions may be useful to describe several physical applications. The results show that these methods are shown to be affective and easy to apply for this type of nonlinear fractional partial differential equations (NFPDEs) with time-dependent coefficients.

1. Introduction

Fractional calculus is considered as a generalization of classical concepts of integration and differentiation. The first appearance of this idea was in 1695. Then, its importance began to increase due to its many applications in several fields such as biology, plasma physics, solid state physic, engineering, economy, finance, liquid crystals, electrical network, numerical analysis, dynamical systems, and control systems [19].

Many mathematicians have proposed different kinds for the definition of the fractional derivative. The most popular of these are Riemann–Liouville, Caputo, Grunwald–Letnikov Erdelyi, Riesz, Hadamard, Marchaud, and Kober [1013].

There are many properties achieved in the classical derivative, but they are not satisfied with the definitions that are mentioned above. The most important of these properties are product rule, quotient rule, chain rule, Rolle theorem, and mean value theorem.

In 2014, Khalil et al. [14] introduced a new definition of fractional derivative, which is called the conformable fractional derivative. Unlike other definitions, this new definition satisfies the properties mentioned above [15].

A number of researchers presented analytical solutions to a large number of NFPDEs with constant coefficients by using the conformable fractional derivative, such as Al-Shawba et al. [16] solved the KdV-ZK equation with time-fractional derivative using the -expansion method. The space-time-fractional modified equal-width equation was solved by Zafar et al. [17] using a hyperbolic function method. A time-fractional biological population model was solved by Zhang and Zhang [18] using a fractional subequation method. Also, Çenesiz and Kurt [19] presented the solution of the space-fractional telegraph equation by introducing the conformable fractional complex transform.

Nowadays, the equations with time-dependent coefficients are more important than equations with constant coefficients as they describe cases that are more general. In 2020, Injrou [20] modified the subequation method to obtain a set of exact solutions for the space-time-fractional Zeldovich equation with time-dependent coefficients.

The Fitzhugh–Nagumo equation system has been derived by both Fitzhugh [21] and Nagumo et al. [22]. It is a simplified form of the Hodgkin–Huxley Model because it is too difficult to be solved analytically.

Recently, many researchers have been interested in the time-fractional Fitzhugh–Nagumo equation with different applications in the areas of neurophysiology, logistical population growth, flame spread, catalytic chemical reaction, and nuclear reactor theory where it combines diffusion and nonlinearity which are controlled by the term :

In 2012, Merdan [23] obtained analytical solutions to the time-fractional Fitzhugh–Nagumo equation (1) by a new application of fractional variational iteration method. At the same year, Pandir and Tandogan [24] presented analytical solution for equation (1) using the modified trial equation method. Ahmet Bekir et al. [25], in 2016, applied -expansion method to obtain exact solutions for equation (1). While in 2017, Bekir et al. [26] solved equation (1) by using the exp-function method. Taşbozan and Kurt [27], in 2020, introduced new exact solutions for equation (1) using the Sine–Gordon expansion method.

In this paper, we utilize the improved subequation method mentioned in [20] and improve the sine-cosine method to solve the time-fractional generalized Fitzhugh–Nagumo equation with time-dependent coefficients and linear dispersion termwhere , , and are arbitrary real-valued function of , is a constant, and is the unknown function depending on the temporal variable and the spatial variable . For , , equation (2) will be reduced to the standard fractional Fitzhugh–Nagumo equation. Equation (2) has never solved in the case where the coefficients are time-dependent in literature before.

The advantages of these methods are that they are more general because they are used in the predicted solution as time-dependent coefficients instead of using constant coefficients and they are effective and easy to apply to NFPDEs with time-dependent coefficients. However, we have noticed that some fractional differential equations are solved with time-dependent coefficients in a variety of methods but with constant coefficients in the predicted solution, for example [28, 29]. On the other hand, these methods require fewer calculations than other methods like the exp-function method or the -expansion method.

2. Preliminaries

2.1. Conformable Fractional Calculus
2.1.1. Definition of Conformable Fractional Derivative [14]

Given a function , the conformable fractional derivative of of order is defined byfor , If is -differentiable in some interval , , and exists, then

Furthermore, if , the definition is equivalent to the classical definition of the first-order derivative of the function.

2.1.2. Definition of Conformable Fractional Integral [14]

The conformable fractional integral of (function continuous) of order , on interval , is defined by function continuous:

2.1.3. Properties of Conformable Fractional Derivative [30]

Some important properties of the conformable fractional derivative and the conformable fractional integral are as follows.

Let and be -differentiable at a point . Then, there are the following properties

Property 1. (differential of a constant).

Property 2. (the linearity property).

Property 3. (product rule).

Property 4. (quotient rule).

Property 5. (chain rule).

Property 6.

Property7. Let and is be -differentiable at a point . If is differentiable, then

Property 8. Let and is any continuous function in a domain of , for , we have

2.1.4. Theorem [31]

Let , be a function such that is differentiable and also -differentiable. Let be a function define in the range of and also differentiable, then

2.2. Description of Two Methods
2.2.1. Description of the Improved Subequation Method

In this section, we summarize the main steps of the improved subequation method for finding exact solutions of (NFPDEs).

Suppose an NFPDEwhere is the conformable fractional derivative of , is an unknown function, and is a polynomial in and its conformable time-fractional partial derivatives and space partial derivatives.

The main steps of the improved subequation method are presented in [18, 20, 32] as follows:Step 1: suppose thatBy substituting (15), NFPDEs turn to the ordinary differential equation (ODE): where is a polynomial in and its derivatives, . Step 2: suppose that the solution of equation (15) can be expressed in the following form: where is the function of to be determined later, and satisfies the following Riccati equation: where is a constant. The following solutions of Riccati equation (19) are given by [33] where  Step 3: determine the positive integer by considering the homogeneous balance between the highest order derivatives and nonlinear terms appearing in equation (17), then substituting (18) with equation (19) into equation (17), then using the properties of the conformable fractional derivative (5)–(14), then collecting all terms with the same order of , and then setting each coefficient of to zero, one can get an overdetermined system of nonlinear differential equations for and . Step 4: assume that and can be obtained by solving the overdetermined system of Step 3, then substituting these results and the solutions, one can obtain the exact solutions of equation (15) immediately.

2.2.2. Description of the Improved Sine-Cosine Method

In this section, the main steps of the sine-cosine method [34,35] are given as follows: Step 1: suppose that where are all functions of to be determined later; equation (15) reduces to ODE where . Now, taking the target equation as considering the following: or Step 2: equating the highest order nonlinear term and highest order linear partial derivative in equation (22), it yields the value of . Step 3: substituting (21) and (23) into equation (15), and collecting all the terms with the same power of , then setting the coefficients to be zero, one can obtain an overdetermined system of nonlinear differential equations for and for . Step 4: assuming that , and can be obtained by solving the overdetermined system of Step 3, then substituting these results and the solutions, the exact solutions of equation (15) can be obtained immediately.

3. Applications of Two Methods

3.1. Applications of the Improved Subequation Method

Using the generalized traveling wave transformation,and Property 2.1.3 and Theorem 2.1.4, equation (2) can be reduced to the following nonlinear ordinary differential equation:where and . We suppose that equation (27) has a solution in the form of (18). Balancing the highest order derivative term and with nonlinear term in equation (27), one can get . So, we have

Substituting (28) into equation (2), collecting all the terms with the same power of, then setting the coefficients of to be zero, one can obtain the overdetermined system of nonlinear differential equations for ,, and as follows:

is a continuous funcation. Solving this system by Maple software, and for simplicity, we introduce the notation.

Finally, one can obtain these four cases: Case 1. If we assume that , , It can be obtained that When , the following hyperbolic function solution of equation (2) is obtained as follows: When , we obtain the following trigonometric function solution of equation (2), as follows: When , substituting into the abovementioned system and then solving the obtained system, one can find will not change but is changed and becomes equal to  When , we have the following rational solution of equation (2), as follows: Case 2. If we assume that , , it can be obtained that We construct the following hyperbolic function solution of equation (2), where is as follows: We construct the following hyperbolic function solution of equation (2), where is as follows: When, we have the following solution of equation (2), as follows: Case 3. If we assume that , , it can be obtained that where When , we obtain the following hyperbolic function solution of equation (2), as follows: When , we obtain the following trigonometric function solution of equation (2), as follows: When , we have the following rational solution of equation (2), as follows: Case 4. If we assume that , , it is obtained that When , we obtain the following hyperbolic function solution of equation (2), as follows:  When , we obtain the following trigonometric function solution of equation (2), as follows: When , we have the following rational solution of equation (2), as follows: where , , , , , and are arbitrary constants.

In Figure 1, we compare the space-time graph of solution (38) with the space-time graph of solution (39) at . In (a), we find the kink-shape wave soliton solution (38) in the plane, while (b) describes the propagation of the anti-kink-shape wave soliton solution in the plane, with the same time-dependent coefficients where .

In Figure 2. (a) clarifies the propagation of the kink-shape wave soliton solution in the plane. The wave moves with the same shape and amplitude over its the track in the domain, while (b) describes the propagation of the spike-shape wave soliton solution in the plane, with the same time-dependent coefficients where .

3.2. Applications of the Improved Sine-Cosine Method

Using the generalized traveling wave transformation,and Property 2.1.3 and Theorem 2.1.4, equation (2) can be reduced to the following nonlinear ordinary differential equation:

Balancing the highest order derivative term with nonlinear term , we can get . So, we have

Substituting (64) and (23) into equation (2) and collecting all the terms with the same power of for and , then setting the coefficients to be zero, we can obtain overdetermined system of nonlinear differential equations for , , and as follows:

Solving this system by Maple, for simplicity, we introduce the notation , where is a continuous function.

Finally, one finds these four cases as follows: Case 1. If we assume that , , it can be obtained that We substitute what mentioned before into (64) and using (24) then, we obtain the hyperbolic functions solution of equation (2), as follows: Case 2. If we assume that , , it can be obtained that where We substitute the solutions of the Case 2 into (64) and using (24), and then we give the hyperbolic functions solution of equation (2), as follows: Case 3. If we assume that , , it is obtained that We substitute the solutions of the Case 3 into (64) and using (24), and then we give the hyperbolic functions solution of equation (2), as follows: Case 4. If we assume that , , it is obtained that We substitute the solutions of the Case 4 into (64) and using (24), and then we give the hyperbolic functions solution of equation (2), as follows: where , , , , and are arbitrary constants.

In Figure 3, it is clearly shown that (a) the space-time graph of solution (73) when clarify the propagation of the kink-shape wave soliton solution and (b) the space-time graph of solution (79) when , represents the propagation of the anti-kink-shape wave soliton solution.

4. Conclusions

In this work, we have obtained three kinds of exact solutions that include the trigonometric function solutions, the hyperbolic function solutions, and the rational solutions which are successfully established by using the subequation method. Also, a kind of exact analytical solutions which are the hyperbolic function solutions are found by using the sine-cosine method and conformable fractional derivative of the time-fractional generalized Fitzhugh–Nagumo equation with time-dependent coefficients. We have found that the subequation method gives more general solutions than the solutions that are given by the sine-cosine method. Moreover, it is observed that the solutions obtained in this research may be important to describe certain nonlinear phenomena in mathematical physics, engineering, and biology. Remarkably, these solutions and the proposed traveling wave transformation have not been reported in other literature. Also, some 3D graphical representations are offered for the obtained results with the different time-dependent coefficients and different values of . We can conclude that the used method is a very powerful, convenient, and efficient technique and it can be used for many other partial fractional nonlinear differential equations.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by Tishreen University.