A Numerical Method for Solving a Class of Nonlinear Second Order Fractional Volterra Integro-Differntial Type of Singularly Perturbed Problems

In this paper, we study a class of fractional nonlinear second order Volterra integro-differential type of singularly perturbed problems with fractional order. We divide the problem into two subproblems. The first subproblems is the reduced problem when e = 0. The second subproblems is fractional Volterra integro-differential problem. We use the finite difference method to solve the first problem and the reproducing kernel method to solve the second problem. In addition, we use the pade’ approximation. The results show that the proposed analytical method can achieve excellent results in predicting the solutions of such problems. Theoretical results are presented. Numerical results are presented to show the efficiency of the proposed method.


Introduction
Volterra integral equations are considered as a type of integral equations.In 1913, Volterra published the first book that talked about Volterra integral equations.In 1884, Volterra began working on integral equations, but his important study was in 1896.However, the name Volterra integral equation was first called by Lalesco in 1908.Volterra integral equations have many applications in science and engineering such as elasticity, semi-conductors, scattering theory, seismology, heat conduction, metallurgy, fluid flow, chemical reactions, population dynamics, and spread of epidemics [1].
Volterra integral equations have growlingly been recognized as useful tools for problems in science and engineering.In [2], they proposed and examined a spectral Jacobi-collocation approximation for fractional order integro-differntial equations.Ray et al. [3], used the Legendre wavelet method to find the solutions for a system of nonlinear Volterra integro-differntial equations.In [4], they used Lagurre polynomials and the collocation method to solve the pantograph-type Volterra integro-differntial equations under some initial conditions.Yang et al. [5], discussed the blow-up of Volterra integro-differntial equations with a dissipative linear term to show the differences of the solutions.In [6], they solved a non-linear system of higher order Volterra integro-differntial equations using the Single Term Walsh Series (STWS) method.Also in [7], they solved the fractional Fredholem-Volterra integro-differntial equations by the fractional-order functions based on the Bernoulli polynomials.We also indicate the interested reader to [8][9][10][11][12][13][14][15][16] for more research works on Volterra integro-differntial equations.
In 1904, A German physicist called Ludwig Prandtl was revolutionized fluid dynamics.He noted that the influence of friction is experienced only very near an object moving through a fluid.In [16], he presented the idea of the boundary layer and its significance for drag and streamlining.In his paper, Ludwig Prandtl assumed that the impact of friction was to cause the fluid instantly adjacent to the surface to stick to the surface.This boundary-layer notion had been the basis stone for the new fluid dynamics.Schlichting was one of the most famous books on boundary layer theory [17].The scientific justification of boundary layer theory gave us more general hypothesis to determine asymptotic expansions of the solutions to the complete equations of the motion.Singular perturbation problem was the result of reduced the problem which was then solved by the method of matched asymptotic expansions.In 1946, Friedrichs and Wasow were the first time used the expression "singular perturbation" [18].
The differential equations of the singularly perturbed problem indicated the study of a group of differential equations including an asymptotically small parameter.The singularly perturbed problem is very important to both applied and pure mathematicians, physicists and engineers because of the fact that the solutions exhibit some interesting behavior.For example, the boundary layer, interior layer, and resonance phenomena [19].
There are a lot of applications of the singularly perturbed problem such as the nonlinear problems of plates and shells by means of the singular perturbation method [20].Petar discussed typical applications of singular perturbation techniques to control problems in the last fifteen years [21].Kokotovic et al., showed results on singular perturbations surveyed as a tool for model order reduction and separation of time scales in control system design [22], Ghorbel and Spong, reviewed results of integral manifolds of singularly perturbed non-linear differential equations and outlined the basic elements of the integral manifold method in the context of control system design [23].Fridman, studied the H ∞ control problem for an affine nonlinear singularly perturbed system [24].Fridman, studied the infinite horizon nonlinear quadratic optimal control problem for a singularly perturbed system [25].Several other techniques to solve such problems are presented in [26][27][28][29][30][31].
In this paper, we consider the following class of fractional nonlinear second order Volterra integro-differntial type of singularly perturbed problems of the form where > 0 is a small positive parameter, y 0 and y 1 are constants, and K(x, t) and f (x) are smooth functions.The derivative which we use in this paper is in the Caputo sense.We organize our paper as follows.In Section 2, we present some preliminaries and the reproducing kernel (RKM) method which we use in this paper.In Section 3, we present some analytical results.In Section 4, we present the proposed method.Some numerical results are presented in Section 5. We end this paper by conclusions which presented in Section 6.

Reproducing Kernel Method
In this section, we present some preliminaries and RKM which we use in this paper.
Definition 1.Let α > 0 and x be a positive real number.Then, the Riemann-Liouville fractional integral of order α is given by where Γ is the Euler gamma function.Some of its properties are given as follows: where a and b are constants.
Definition 2. Let α > 0 and x be a positive real number.Then, the Caputo derivative of order α is given by Some of its properties are given as follows: where a and b are constants,
The second condition is called the reproducing property and a Hilbert space which possesses a reproducing kernel is called a reproducing kernel Hilbert space.
Consider the second order nonlinear fractional equation of the form where θ and φ are constants.First, we study the linear case when g(y) = a(x).To homogenize the initial condition, we assume u = y − φx − θ(1 − x).Thus, Equations ( 3) and ( 4) can be rewritten as In order to solve the linear Equations ( 5) and ( 6), we construct the kernel Hilbert spaces W 1 2 [0, 1] and W 3 2 [0, 1] in which every function satisfies the initial condition (6).Let The inner product in W 1 2 [0, 1] is defined as and the norm u W 1 2 [0,1] is given by where Theorem 1.The space W 1 2 [0, 1] is a reproducing kernel Hilbert space, i.e.; there exists R(s, y) ∈ W 1 2 [0, 1] and its second partial derivative with respect to y exists such that for any u Using the integration by parts, one can get and Since the characteristic equation of − ∂ 2 R ∂y 2 (s, y) = δ(y − s) is λ 2 = 0 and its characteristic value is λ = 0 with 2 multiplicity roots, we write R(s, y) as Using Conditions ( 8)-( 11), we get the following system of equations Next, we study the space and the norm u W 3 2 [0,1] is given by Theorem 2. The space W 3 2 [0, 1] is a reproducing kernel Hilbert space, i.e.; there exists K(s, y) ∈ W 3 2 [0, 1] which has its six partial derivative with respect to y such that for any u ∈ W 3 2 [0, 1] and each fixed y, s ∈ (1 + s 5 ), Proof.Using integration by parts, one can get Since u(y) and K(s, y (u(y), K(s, y) where δ is the dirac-delta function and Since the characteristic equation of ∂ 6 K ∂y 6 (s, y) = δ(s − y) is λ 6 = 0 and its characteristic value is λ = 0 with 6 multiplicity roots, we write K(s, y) as On the other hand, integrating ∂ 5 K ∂y 5 (s, y) = δ(s − y) from s − to s + with respect to y and letting → 0 to get Using the Conditions 13 and 15-20, we get the following system of equations We solved the last system using Mathematica to get (1 + s 5 ), Now, we present how to solve Equations ( 5) and ( 6).Let where L(σ i (s)) = D α σ i (s) + a(s)σ i (s) and L * is the adjoint operator of L. Using Gram-Schmidt orthonormalization to generate orthonormal set of functions and α ij are coefficients of Gram-Schmidt orthonormalization.In the next theorem, we show the existence of the solution of Equations ( 5) and (6).
Proof.First, we want to prove that {ψ For each fixed u(s Second, we prove Equation (22).Simple calculations imply that Let the approximate solution of Equations ( 5) and ( 6) be given by In the next theorem, we show the uniformly convergence of the ds for m = 0, 1, 2. Proof.First, we prove the theorem for m = 0.For any s ∈ [0, 1], From Theorem (4), one can see that which implies that {u N (s)} ∞ N=1 converges uniformly to u(s).Second, we prove the uniformly convergence for m = 1, 2. Since d m K(s,y) Thus, for any s ∈ [0, 1], which implies that Therefore, converges uniformly to ds m for m = 1, 2. Now, we discuss how to solve Equations ( 3) and ( 4).Let L(y(x)) = D α y(x) and N(y(x)) = g(x, y)y are the linear and nonlinear parts of Problem (3), respectively.We construct the homotopy as follows: where λ ∈ [0, 1] is an embedding parameter.If λ = 0, we get a linear equation D α y(x) = 0 which implies that y(x) = θ.If λ = 1, we turn out to be Problem (3).Following the Homotopy Perturbation method [34], we expand the solution in term of the Homotopy parameter λ as Substitute Equation (25) into Equation ( 24) and equating the coefficients of the identical powers of λ to get the following system To solve the above equations, we use the RKM which described above and we obtain where From Equation ( 26), it is easy to see the solution to Equations ( 3) and ( 4) is given by We approximate the solution of Equations ( 3) and ( 4) by

Analytical Results
In this section, we present the maximum principle, the stability theorem, and the uniqueness theorem.Firstly, Equations ( 1) and ( 2) is transformed into an equivalent problem as follows: The following conditions are needed in order to guarantee that Equations ( 1) and ( 2) do not have the turning-point problem for all x ∈ [0, 1], where k 1 , k 2 , k 3 , and k 4 are positive constants and y ∈ C 2 (0, 1) ∪ C[0, 1].
In the next theorem, the stability result is presented.

Solution Method
Consider the following of class of nonlinear second order fractional nonlinear Volterra integro-differential type of singularly perturbed problems of the form where > 0 is a small positive parameter, y 0 and y 1 are constant, and K(x,t) and f (x) are smooth functions.To solve Equations ( 1) and ( 2), we use the following steps.
Step 1: A reduced subproblem is obtained by setting = 0 in Equation ( 1) to get On most of the interval, the solution of Equation ( 34) behaves like the solution of Equations ( 1) and (2).However, there is small interval around x = 0 in which the solution of problem (1) and ( 2) does not agree with the solution of problem (1) and (2).To handle this situation, the boundary layer correction subproblem is introduced in step 2.
Step 2: Choose x = 1 α−1 s to get Thus, Hence, Equation (1) becomes or Setting = 0 in Equation ( 36) implies that Since the solution of the reduced problem in step 1 does not satisfy the initial condition at x = 0, then the solution of the above equation should satisfy it.This means, its solution has the form y 1 (0) + y 2 (x).Substitute y(x) = y 1 (0) + y 2 (x) in Equation (38) to get the boundary layer correction equation The solution of Equation (1) will be expressed in the form as and the initial condition must be satisfied by expression (40).When x = 0, the condition will be The solution of Equations ( 1) and ( 2) can be produced using the RKM as described in the previous section.More details can be found in [36][37][38][39].

Numerical Results
In this section, we present two of our examples to show the efficiency of the proposed method.
subject to y(0) = 0, y(1) = 0. (43) When = 0, we get We discretized the interval [0, 1] by x i = ih, h = 1 n , n ∈ N. Let y k ≈ y(x k ) for k = 0 : n.Using the backward finite difference method to approximate y (x k ) and the trapezoidal quadrature to approximate the integral Thus, we get the following system .
Using Mathematica, one can see that the solution of the above system for n = 12 is given in Figure 1.Using the change of variable x = 2 s, we get Using the RKM, we get Using the Pade' approximation of order [2,2], we have θ = −0.693147.In Figures 2-4, we plot the approximate solution for = 0.0001, 0.00001, and 0.000001, respectively.Let In Table 1, we present the error for x = 0, 0.1, ..., 1 for = 0.0001, 0.00001, and 0.000001.In Table 2, we present the computational time when = 0.0001, 0.00001, and 0.000001.When = 0, we get We discretized the interval [0, 1] by Using the backward finite difference method to approximate y (x k ) and the trapezoidal quadrature to approximate the integral Using Mathematica, one can see that the solution of the above system for n = 12 is given in Figure 5.Using the change of variable x = 2 s, we get −D Using the Pade' approximation of order [2,2], we have θ = 0.0927388622769557.In Table 3, we present the error for x = 0, 0.1, ..., 1 for = 0.001, 0.0001, and 0.00001.In Table 4, we present the computational time when = 0.001, 0.0001, and 0.00001.In Figures 6-8, we plot the approximate solution for =0.001, 0.0001, and 0.00001, respectively.

Conclusions
In this paper, we study a class of fractional nonlinear second order Volterra integro-differential type of singularly perturbed problems with fractional order.We divide the problem into two subproblems.The first subproblem is the reduced problem when = 0.The second subproblem is the second order fractional Volterra integro-differential problem.We use the finite difference method to solve the first subproblem and the reproducing kernel method to solve the second subproblem.The results show that the proposed analytical method can achieve excellent results in predicting the solutions of such problems.Theoretical results are presented.Numerical results are presented to show the efficiency of the proposed method.Figures 1-8 show the efficiency of the proposed method.

Example 1 :
Consider the following problem
Figure 1.The approximate solution y 1 .