The New Approximate Analytic Solution for Oxygen Diffusion Problem with Time-Fractional Derivative

Oxygen diffusion into the cells with simultaneous absorption is an important problem and it is of great importance in medical applications. The problem is mathematically formulated in two different stages. At the first stage, the stable case having no oxygen transition in the isolated cell is investigated, whereas at the second stage the moving boundary problem of oxygen absorbed by the tissues in the cell is investigated. In oxygen diffusion problem, a moving boundary is essential feature of the problem. This paper extends a homotopy perturbation method with time-fractional derivatives to obtain solution for oxygen diffusion problem. The method used in dealing with the solution is considered as a power series expansion that rapidly converges to the nonlinear problem. The new approximate analytical process is based on two-iterative levels. The modified method allows approximate solutions in the form of convergent series with simply computable components.


Introduction
The diffusion of oxygen into absorbing tissue was first studied in [1].First the oxygen is allowed to diffuse into a medium, some of the oxygen can be absorbed by the medium, and concentration of oxygen at the surface of the medium is maintained constant.This phase of the problem continues until a steady state is reached in which the oxygen does not penetrate any further and is sealed so that no oxygen passes in or out, the medium continues to absorb the available oxygen already in it, and, as a consequence, the boundary in the steady state starts to recede towards the sealed surface.
Crank and Gupta [2] also employed uniform space grid moving with the boundary and necessary interpolations are performed with either cube splines or polynomials.Noble [3] suggested repeated spatial subdivision, Reynolds and Dolton [4] also developed the heat balance integral method, and Liapis et al. [5] proposed an orthogonal collocation for solving the partial differential equation of the diffusion of oxygen in absorbing tissue.Gülkac ¸proposed two numerical methods for solving the oxygen diffusion problem [6].Mitchell studied the accurate application of the integral method [7].More references to this problem may be found in [8][9][10][11][12][13][14][15][16][17].
In recent years, fractional differential equations have drawn much attention.Many important phenomena in physics, engineering, mathematics, finance, transport dynamics, and hydrology are well characterized by differential equations of fractional order.Fractional differential equations play an important role in modelling the so-called anomalous transport phenomena and in the theory of complex systems.These fractional derivatives work more appropriately compared with the standard integer-order models.So, the fractional derivatives are regarded as very dominating and useful tool.For mathematical properties of fractional derivatives and integrals one can consult [18][19][20][21][22][23].
In the present work, we extend a homotopy perturbation method with time-fractional derivatives to obtain solution for oxygen diffusion problem.
We give some basic definitions of fractional derivatives as follows.
Definition 1.The Riemann-Liouville fractional integral of  ∈   of the order  ≥ 0 is defined as Mathematical Problems in Engineering where Γ denotes gamma function: Definition 2. The fractional derivatives of  ∈   of the order  ≥ 0, in Caputo sense, are defined as for  − 1 <  ≤ ,  ∈ ,  > 0,  ∈    ,  ≥ −1.

Analysis of Homotopy Perturbation Method with Time-Fractional Derivatives
Let us assume that nonlinear fractional differential equation is as follows: with the initial condition (, 0) = , where  is the operator,  is known functions, and (, ) is sough functions.Assume that operator  can be written as () = ()+(), where  is the linear operator and  is the nonlinear operator.Hence (7) can be written as follows: For solving (7) by homotopy perturbation method, we construct the homotopy or the equivalent one where  ∈ [0, 1] is an embedding or homotopy parameter, (, ; ) : Ω[0, 1] → , and  0 is the initial approximation for solution equation (8).
Clearly, the homotopy equations (, 0) = 0 and (, 1) = 1 are equivalent to the equations     −  0 = 0 and   −()−()−(, ) = 0, respectively.Thus, a monotonous change of parameter  from 0 to 1 corresponds to a continuous change of the trivial problem     −  0 = 0 to the original problem.Now, we assume that the solution of ( 8) can be written as a power series in embedding parameter , as follows: where  0 and  1 are functions which should be determined.Now, we can write (11) in the following form: ) Apply the inverse operator,    , which is the Riemann-Liouville fractional integral of order  > 0.
On both sides of ( 12), we have Suppose that the initial approximation of solutions equation ( 8) is in the following form: where   () for  = 1, 2, . . .are functions which must be computed.Substituting (11) and ( 14) into (13) we get Synchronizing the coefficients of the same powers leads to Mathematical Problems in Engineering 3 Now, we obtain the coefficients   (),  = 1, 2, . ... Therefore the exact solution can be obtained as follows: Efficiency and reliability of the method are shown.

Problem Description and Formulation
Crank and Gupta [1] were the first researchers to model oxygen diffusion problem mathematically.
The process includes two mathematical levels.At the first level, the stable condition occurs once the oxygen is injected into either the inside or outside of the cell; then the cell surface is isolated.
At the second level, tissues start to absorb the injected oxygen.The moving boundary problem is caused by this level.The aim of this process is to find a balance position and to determine the time-dependent moving boundary position.Writing down the time-fractional derivatives of oxygen diffusion problem in [1] is adopted, following [18].

Solution of Fractional Oxygen Diffusion Problem
We consider the following oxygen diffusion problem: with the following initial and boundary conditions: where 0 <  ≤ 1.To solve ( 18)-( 21) by present method, we construct the following homotopy: or where  ∈ [0, 1] and 0 <  ≤ 1.Consider Assume that the initial approximation of solutions equation ( 18) is in the following form: where   () for  = 1, 2, . . .are functions which must be computed.
Applying the inverse operator    of    on both sides of (24) we obtain  (, ) =  (, 0) +     0 (, ) Suppose the solution of (26) has the following form: where (, ) is functions which should be determined.Substituting ( 27) into (26), collecting the same powers of , and equating each coefficient of  to zero yield where then and from ( 29) and (30) we obtain (33) therefore, we obtain solution of (26) for  = 1: and we obtain the following solution: We can now obtain an expression for the location of the moving boundary ().We can write following [18]: and following initial and boundary conditions We construct the following homotopy for moving boundary as or where  ∈ [0, 1] and 0 <  ≤ 1.Consider applying the inverse operator    of    on both sides of (44), we obtain suppose the solution of (45) has the following form: Substituting ( 46) into (45), collecting the same powers of , and equating each coefficient of  to zero yield , where so Then we have from (46) if we let  = 1 Therefore, we obtain the solutions of moving boundary condition as or ] . (52)

Numerical Simulations
In this section numerical results for the solution of the oxygen diffusion problem using the constructed homotopy perturbation method with the time-fractional derivative are presented.These proposed homotopy perturbation methods are applied and figures present solutions are presented using different values for the derivative order .Figures 1(a

Conclusion
In this study, we extended homotopy perturbation method with time-fractional derivative to find the exact solution of oxygen diffusion problem with moving boundary.It is effortless and also easy to apply and we can say that the present method is an effective method and has appropriate technique to find the exact solution to many complex problems.