Abstract

In this paper, DSEK model with fractional derivatives of the Atangana-Baleanu Caputo (ABC) is proposed. This paper gives a brief overview of the ABC fractional derivative and its attributes. Fixed point theory has been used to establish the uniqueness and existence of solutions for the fractional DSEK model. According to this theory, we will define two operators based on Lipschitzian and prove that they are contraction mapping and relatively compact. Ulam-Hyers stability theorem is implemented to prove the fractional DSEK model’s stability in Banach space. Also, fractional Euler’s numerical method is derived for initial value problems with ABC fractional derivative and implemented on fractional DSEK model. The symmetric properties contribute to determining the appropriate method for finding the correct solution to fractional differential equations. The numerical solutions generated using fractional Euler’s method have been plotted for different values of where and different step sizes . Result discussion will be given, describing the changes that occur due to the step size .

1. Introduction

Fractional calculus is as historic as integer calculus but not until 1819, it was properly introduced in the form of definitions and functions. Many scientists, researchers, and mathematicians played their role in developing its theory such as [15]. Since fractional calculus was developed theoretically at first and had no practical application at the time, it was not as well known as integer calculus among other areas of science. However, after the contribution of Professor Mandelbrot’s fractal theory, fractional calculus theory developed rapidly and soon became the hot topic among all researchers around the globe.

The existing theory of nonlinear science is now seemed to be only focused on fractional order calculus theory and the theory of chaos and dissipative structure ([6, 7]). The fact that fractional calculus describes the heredity and memory of any physical phenomenon is fascinating [8, 9]. As a result, it is now used more than integer calculus in fluid dynamics, quantum mechanics, mathematical biology applications, chemistry, control and signal theory, economics, image processing, etc. Models developed in many areas of science and engineering are observed to be best explained by fractional differential equations. The symmetries can be found by solving a related set of partial fractional differential equations. Since integer-order models lack memory and heredity, they cannot adequately and sufficiently describe physical phenomena in many cases. These applications have also led to the rapid development of fractional calculus theory. The authors of [1013] have a great deal of literature on the subject describing applications and types of fractional derivatives. It is critical to note that every one of those fractional derivative order definitions has its own advantages and disadvantages.

Aside from the mathematical satisfactions of the fractional-order Atangana-Baleanu derivative, the new derivative is being studied due to the necessity of implementing a model depicting the behavior of orthodox viscoelastic materials, thermal medium, and other materials. The proposed mechanism can depict material heterogeneities as well as some structure or media at multiple scales.

The new kernel’s nonlocality enables the full description of memory inside of structure and media with multiple scales, which cannot be represented by classical fractional derivatives or those of the Caputo-Fabrizio type. Furthermore, we believe that Atangana-Baleanu derivatives can play an important role in the study of the microstructural behavior of some materials, particularly those involving nonlocal exchanges, which are important in defining the material’s properties states [14]. Atangana-Baleanu derivatives are thus extremely useful in describing a wide range of scientific, engineering, and technological problems.

Descemet’s stripping endothelial keratoplasty (DSEK) is the name given to eye surgery [1517] in which a damaged corneal layer is replaced with a healthy corneal layer from a donor or synthetic cornea. Cornea is a clear layer of the eye that is very important in the anterior part of the eye; if it is scratched or damaged, it affects vision. It itself is made up of five layers, see Figure 1. Keeping this in mind, the authors of [18, 19] created the DSEK model, which predicts the behavior of ocular parameters posttreatment. Because this procedure has never been studied mathematically as an ordinary system of differential equations, this work is extremely important. Since fractional calculus has been said to be the generalization of integer calculus, the fractional DSEK model is developed and studied theoretically and numerically in this paper.

In this paper, the first section gives an overview of the literature background of fractional calculus and the DSEK model. The second section gives the preliminary concepts of fractional calculus that will be used in this work. Section three is based on the explanation of the fractional DSEK model. Section four shows the existence of a solution by the fixed point theory of the fractional DSEK model. Section five describes the Ulam-Hyers stability analysis of fractional DSEK. Section 6 describes the computation of fractional Euler’s method for ABC fractional derivative and the application of fractional Euler’s method to fractional DSEK. The last section is the discussion of the results obtained and the conclusion of this paper.

2. Preliminaries of Fractional Calculus

2.1. Atangana-Baleanu Caputo Fractional Derivative

Definition 1. The authors of [20] introduced a new Caputo fractional derivative aswhere is the normalization function that follows the condition .

Definition 2. If the function does not follow the condition then it takes the form asThis equation can also take the form of the condition where , also .
This derivative was defined by [20] to involve an exponential kernel in fractional derivatives to represent the results of dynamic systems memory effects more accurately. With the passage of time, it occurred that this definition has a flaw in that it does not give the original function when . To overcome this problem, the authors of [21] presented the accurate kernel and modified this definition accordingly.

Definition 3. Let the new fractional derivative be defined aswhere also and has the same properties defined in [20]. Here, is the generalized Mittag-Leffler function defined as .
For the above definition, the constant function has a fractional derivative of zero. The above description would be helpful when solving real-world issues and will also provide a great benefit in utilizing the Laplace transform to solve any initial state physical problem. Nevertheless, if alpha is 0, we will not recover the initial function except when the function vanishes at the origin. We suggest the following definition, in order to avoid this problem.

Definition 4. Let the new fractional derivative be defined aswhere also and has the same properties defined in [20]. Here, is the generalized Mittag-Leffler function defined as .
Both definitions have a nonlocal kernel. For calculations in this paper, we will use definitions in (4) and (5).

2.2. Properties of Atangana-Baleanu Caputo Fractional Derivative

(i)Laplace transformation on equation (4)is(ii)Laplace transformation on equation (5)is(iii)Let then the following relation exists [20]:(iv)If is a continuous function on some closed interval . Then, the following inequality can be written on (v)Lipschitz condition Atangana-Baleanu Caputo fractional derivative satisfies the Lipschitz condition in Riemann and Caputo sense, and the following inequality exists:Similarly, for (5), the Lipschitz condition exists as(vi)AB fractional integral for the AB fractional integral for and nonlocal kernel is given as

When the ordinary integral is obtained, and for the initial function is obtained.

For proof of these, see [20].

Lemma 1. [21] Suggests that the proposed problem for has a solution; that is,its solution is given by

3. Fractional DSEK Model

As mentioned in Section 2, the definitions in (4) and (5) have nonlocal kernels, and therefore, Atangana-Baleanu Caputo fractional derivative operator’s performance in modeling eye surgery is better than any other definition. It inspired the valuable applications of several fractional operators in dynamic mathematical models; therefore, we are researching the dynamics of eye surgery derived in [18] by a system of nonlinear differential equations by involving fractional derivative.with initial conditionswhere is the Atangana-Baleanu Caputo fractional derivative of order . DSEK model is based on the same conditions given by [18]. Also, the defined parameters have the same description as given by [18]. Such as is the refractive index, is the axial length, is the corneal curvature, and is the central corneal thickness.

3.1. Preliminaries for Fractional DSEK Model

For fractional analysis of the DSEK model, let us define . To define the Banach space, let us say we have where . Then, the field can be written as under the norm supremum aswhere . Also, .

Definition 5. Let be a Banach space. Then, defined as will be a Lipschitzian if there exists a constant for which the inequality exists such thatfor all . Where is the Lipschitz constant for . If then is a contraction.

Theorem 1. Let B be a Banach space and be a contraction mapping. Then, there must exist a unique fixed point of .

Theorem 2. A subset of Banach space B is supposed to be . Let be convex, closed, and nonempty. Suppose that and map into , and the following relations exist:(i)(ii)F is continuous and compact(iii) is a contraction mappingThen, there exists such that .

4. Existence of Solutions for Fractional DSEK Model

By using the fixed-point theory, let us prove the uniqueness and existence of the DSEK model. To prove its uniqueness and existence, let us reformulate the DSEK model of (14).where

Let us consider system (14) aswith an initial condition where

In (21), the superscript represents the transpose. By using Lemma 1 and fractional integral, the (20) becomes the fractional integral equation as

Now, to prove the existence uniqueness, we consider two hypotheses based on Lipschitzian and some growth condition assumptions.

Hypothesis 1. For two constants , the inequality exists; that is,

Hypothesis 2. For a constant such thatfor each and .
Let us define two operators and aswhere .

Theorem 3. Consider a closed convex set where such that and prove thatfor This confirms that

Theorem 4. Prove that is a contraction.
To prove that is a contraction suppose . Then, by using Hypothesis 2, we haveAs we know that is a contraction mapping.

Theorem 5. Prove that is relatively compact.
We can prove that is relatively compact by showing that is continuous, uniformly bounded, and also equicontinuous.
As we know that is continuous, then is also continuous.
Let us assume that thenHence, proved that is uniformly bounded on . Now, we have to show that is equicontinuous. Assume and where . Then, we haveNow, the Arzelá-Ascoli theorem suggests that is relatively compact, and hence, it is completely continuous.

Theorem 6. If Hypothesis 1 and Hypothesis 2 hold, then the fractional integral equation that is equation (20) which is the solution of equation (12) has at least one solution only if where isBy using Theorems 25, it is proved that the integral equation given in (22) has at least one solution, and consequently, the DSEK model (14) under consideration also has at least one solution.

Theorem 7. Prove that integral equation (20) has a unique solution if under Hypothesis 2.
As we have defined asLet and . Then, we haveHence, suggests is a contraction. Hence, (22) has a unique solution which suggests that (14) also has a unique solution.

5. Ulam-Hyers Stability for DSEK Model

Stability analysis of nonlinear dynamical models is a must. So, in this work, we use Ulam-Hyers stability for DSEK model12 with some nonlinear functional analysis concepts. Ulam-Hyers stability was introduced in 1940 by [22, 23] as a stability study for functional equations. This acted as a motivator for various researchers, and then, this stability was discussed in many forms. Using the fixed-point technique, the authors in [24] investigated the Hyers-Ulam-Rassias and Hyers-Ulam stability of the fractional Volterra integral-differential equation. In a Banach space, some results on generalized Hyers-Ulam stability of the linear differential equation were introduced in [25]. In [26], the authors investigated the Hyers-Ulam stability of first-order linear differential equations and extended previous results using the integral factor approach. In [27], the Hyers-Ulam-Rassias stability of a certain fractional differential equation was discussed, as well as the Hyers-Ulam stability of a certain fractional differential equation. For a particular family of fractional integrodifferential equations, the stability of Ulam-Hyers, Ulam-Hyers-Rassias, and semi-Ulam-Hyers-Rassias on some intervals was studied in [28]. The Ulam-Hyers and generalized Ulam-Hyers-Rassias stabilities for the solution of a fractional-order pseudoparabolic partial differential equation were investigated using the Gronwall inequality [29]. The existence and uniqueness of solutions, as well as Ulam-Hyers-Rassias stability, of an impulsive certain fractional differential equation were investigated in [30]. Sometimes, it is the stability analysis of differential equation ordinary or partial, integral equations, functional equations, etc. Various types have been formed of Ulam-Hyers stability theory, namely, Ulam-Hyers-Rassias, semi-Ulam-Hyers-Rassias [28, 31], and Ulam stability [32].

Definition 6. For some ifthere must exist that satisfies the DSEK model (14) having an initial conditionwhere such thatwith this property, if there exists an , then it is said that the DSEK model (12) is UlamHyers stable.

Remark 1. Let be a small perturbation such that where has the properties given as follows:(i), where and (ii)For the model becomeswhere the superscript represents the transpose.

Lemma 2. Perturbed system (35)has a solution that satisfies the inequality represents the solution of the system (37),By using Remark 1 and Lemma 2, the solution of system (37) is given asAlso, we know thatNow, Remark 1 suggests that

Theorem 8. By using Theorem 4, it is proven that the DSEK system (35) is Ulam-Hyers stable in B. Let the DSEK system (12) with initial conditionshas a unique solution as and is the solution of inequality (34), and then,Since as suggested by an initial condition, hence, (45) becomesThen, by Lemma 2 and the hypothesis above, we havewhich implies thatSince hence from we obtain . Hence, proved that the DSEK system (35) is Ulam-Hyers stable.

6. Numerical Approximation of Fractional DSEK Model

In this section, the fractional DSEK model will be solved numerically by using fractional Euler’s method. There are several numerical techniques to compute the numerical results of a fractional system of differential equations, but in this case, even Euler’s method can analyze its solution. In order to do that, we first derive the fractional Euler’s method for Atangana-Baleanu Caputo fractional derivative.

6.1. Fractional Euler’s Method for Atangana-Baleanu Caputo Fractional Derivative

The authors of [33] proved that the generalized Taylor’s formula of Atangana-Baleanu Caputo fractional derivative is given as

Suppose that we have an initial value problem

Let be the interval on which we need to obtain the solution of our problem. For generalization instead of , we consider and use this for our approximation. Let the k subintervals of equal width be by using nodes for . Consider that etc. are continuous on then, by using (50), we expand about as

Upon neglecting higher-order terms because step size is considered as a smallest positive number and taking , (51) becomes

Equation (52) becomes the iterative equation for repeatedly calculating the points of that approximates the solution of . Hence, the general form of fractional Euler’s method for solving initial value problems with Atangana-Baleanu Caputo fractional derivative is

It can be observed easily that for , this becomes the classical Euler method.

Now, to solve the fractional DSEK model numerically, we use the parameter values and initial conditions given in [18]. According to that table, , and . Now, the fractional DSEK system in (14) with the iterative formula (54) becomes

By solving (55) with the help of software, we obtain the numerical solution for different values of in the form of Figures 25 and 6. Figure 2(a) represents the fractional solution of refractive index for different values of between and , whereas Figure 2(b) is obtained for the step size .

The only difference among the solutions presented in Figures 2(a) and 2(b) is the different values of . If we observe as , the fractional behavior is clearer to understand and gives us the accurate approximation for then .

Figures 3(a) and3(b) represent the numerical solution of axial length for different values of where and , respectively. By observing closely, the solutions depicted in Figure 3(a) show that as . Since this model represents a real-life case of eye surgery, hence, this result is unacceptable. As for , the graphical results in Figure 3(b) are more accurate because it suggests the or lies closer to as

Figures 4(a) and4(b) are the graphical illustration of corneal curvature for different values of and step sizes as and , respectively. Similar to the refractive index and axial length, the corneal curvature also depicts more realistic behavior when .

Figures 5(a) and 5(b) represent the numerical solution of corneal thickness for different values of where and , respectively. By observing closely, the solutions depicted in Figure 5(a) show that as . Since this model represents the real-life case of eye surgery, hence, this result is unacceptable.

As for , the graphical results in Figure 5(b) are more accurate because it suggests the or lies closer to as.

Figures 6(a) and 6(b) show the numerical solution of fractional DSEK by fractional Euler’s method. In Figure 6(a), results have been presented for and different values of between .

Graphical results are shown in Figures 25 and 6 described that as , the more accurate results we obtain. This is why the variables remain in for instead of .

Also, the fractional DSEK results have the hysteresis phenomenon, which means this system is influenced by the previous derivatives and values as well as the current conditions. The noninteger derivative given by different values of introduces the memory effect in the fractional DSEK model. As we explained in definitions (4) and (5), the exponential kernel when applied to the fractional DSEK model calculates the memory effect. This is why we can see the smoothness in Figures 25 and 6 as compared to graphical results in [18]. Results of refractive index, axial length, corneal curvature, and central corneal thickness are shown graphically of the ordinary system of differential equation in [18] showed huge oscillation whereas, in real life after surgery, the effect on vision is not that blurry or oscillated. The fractional DSEK model shows more realistic results of ocular parameters after Descemet’s stripping endothelial keratoplasty. It gives the same normal values but due to its fractal phenomenon, the oscillation among results is removed, and graphs are smoother giving the same normal values as the DSEK model in [18]. For more background about the numerical solutions of fractional-order differential equations, see [3437].

7. Conclusion

In this paper, we investigated the fractional DSEK model presented by fractional derivatives of the Atangana-Baleanu Caputo type. We proved the uniqueness and existence of its solutions by using fixed point theory. For this, we defined a hypothesis based on Lipschitzian and two operators . Then, we proved that and are contraction and relatively compact and hence proven the uniqueness and existence of those defined hypotheses. Furthermore, for the fractional DSEK model, proving its stability was a must so by UlamHyers stability in Banach space, we proved that fractional DSEK is Ulam-Hyers stable. Moreover, we have discussed the advantages of using the ABC fractional derivative instead of any other. In this paper, we presented and investigated the fractional behavior of the DSEK model and performed the numerical investigation using mathematical software. The numerical method “Euler” which is used to solve fractional DSEK is derived for initial value problems with ABC fractional derivatives, in this paper. Since eye surgery is a crucial process and with the passage of time, the results of surgery can be observed but with the help of the fractional DSEK model, a clearer picture of this surgery will be given.

Data Availability

The data used to support the findings of this study are available from the corresponding author on request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.