Research Article | Open Access

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

**Academic Editor:**Evangelos J. Sapountzakis

#### Abstract

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.

#### 1. 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ülkaç 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–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–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 is defined as where denotes gamma function:

*Definition 2. *The fractional derivatives of of the order , in Caputo sense, are defined asfor , , , , .

*Definition 3. *The Caputo-time-fractional derivative operator of order is defined as

Lemma 4. *Let , , and , ; then*

Lemma 5. *If , , and then one has *

The Mittag-Leffler function plays a very important role in the fractional differential equations, in fact introduced by Mittag-Leffler in 1903 [24]. Mittag-Leffler function .

#### 2. Analysis of Homotopy Perturbation Method with Time-Fractional Derivatives

Let us assume that nonlinear fractional differential equation is as follows:with the initial condition , 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 onewhere is an embedding or homotopy parameter, , and is the initial approximation for solution equation (8).

Clearly, the homotopy equations and are equivalent to the equations and , respectively. Thus, a monotonous change of parameter from 0 to 1 corresponds to a continuous change of the trivial problem 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 and 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 .

On both sides of (12), we haveSuppose that the initial approximation of solutions equation (8) is in the following form:where for are functions which must be computed. Substituting (11) and (14) into (13) we get Synchronizing the coefficients of the same powers leads toNow, we obtain the coefficients , . Therefore the exact solution can be obtained as follows:Efficiency and reliability of the method are shown.

#### 3. 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].

#### 4. Solution of Fractional Oxygen Diffusion Problem

We consider the following oxygen diffusion problem:with the following initial and boundary conditions: where . To solve (18)–(21) by present method, we construct the following homotopy:orwhere and . ConsiderAssume that the initial approximation of solutions equation (18) is in the following form:where for are functions which must be computed.

Applying the inverse operator of on both sides of (24) we obtainSuppose 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 yieldwherethen and from (29) and (30) we obtaintherefore, we obtain solution of (26) for :orand we obtain the following solution:If thenIf thenWe can now obtain an expression for the location of the moving boundary . We can writefollowing [18]:and following initial and boundary conditions We construct the following homotopy for moving boundary asorwhere and . Considerapplying the inverse operator of on both sides of (44), we obtainsuppose the solution of (45) has the following form:Substituting (46) into (45), collecting the same powers of , and equating each coefficient of to zero yieldwhere soThen we have from (46) if we let Therefore, we obtain the solutions of moving boundary condition asor

#### 5. 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) and 1(b) show the surface concentration for , Figures 2(a) and 2(b) show the surface concentration for , and Figures 3(a) and 3(b) show the surface concentration for . Figures 4(a) and 4(b) show position of moving boundary for , Figures 5(a) and 5(b) show position of moving boundary for , and finally Figures 6(a) and 6(b) show position of moving boundary for .

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

#### 6. 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.

#### Competing Interests

The author declares that she has no competing interests.

#### References

- J. Crank and R. S. Gupta, “A moving boundary problem arising from the diffusion of oxygen in absorbing tissue,”
*IMA Journal of Applied Mathematics*, vol. 10, no. 1, pp. 19–33, 1972. View at: Publisher Site | Google Scholar - J. Crank and R. S. Gupta, “A method for solving moving boundary problems in heat flow using cubic splines or polynomials,”
*Journal of the Institute of Mathematics and Its Applications*, vol. 10, pp. 296–304, 1972. View at: Publisher Site | Google Scholar | MathSciNet - B. Noble, “Heat balance method in melting problems,” in
*Moving Boundary Problems in Heat Flow and Diffusion*, J. R. Ockendon and W. R. Hodkins, Eds., pp. 208–209, Clarendon Press, Oxford, UK, 1975. View at: Google Scholar - W. C. Reynolds and T. A. Dolton, “Use of integral methods in transient heat transfer analysis,”
*ASME Paper*58-A-248, 1958. View at: Google Scholar - A. I. Liapis, G. G. Lipscomb, O. K. Crosser, and E. Tsiroyianni-Liapis, “A model of oxygen diffusion in absorbing tissue,”
*Mathematical Modelling*, vol. 3, no. 1, pp. 83–92, 1982. View at: Publisher Site | Google Scholar | Zentralblatt MATH - V. Gülkaç, “Comparative study between two numerical methods for oxygen diffusion problem,”
*Communications in Numerical Methods in Engineering with Biomedical Applications*, vol. 25, no. 8, pp. 855–863, 2009. View at: Publisher Site | Google Scholar | MathSciNet - S. L. Mitchell, “An accurate application of the integral method applied to the diffusion of oxygen in absorbing tissue,”
*Applied Mathematical Modelling*, vol. 38, no. 17-18, pp. 4396–4408, 2014. View at: Publisher Site | Google Scholar | MathSciNet - E. Hansen and P. Hougaard, “On a moving boundary problem from biomechanics,”
*Journal of Institute of Mathematical and its Applications*, vol. 13, no. 3, pp. 385–398, 1974. View at: Publisher Site | Google Scholar - R. S. Gupta and A. Kumar, “Variable time-step method with coordinate transformation,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 44, no. 1, pp. 91–103, 1984. View at: Publisher Site | Google Scholar | Zentralblatt MATH - J. V. Miller, K. W. Morton, and M. J. Baines, “A finite element moving boundary computation with an adaptive mesh,”
*Journal of the Institute of Mathematics and Its Applications*, vol. 22, no. 4, pp. 467–477, 1978. View at: Publisher Site | Google Scholar | MathSciNet - S. G. Ahmed, “A numerical method for oxygen diffusion and absorption in a sike cell,”
*Applied Mathematics and Computation*, vol. 173, no. 1, pp. 668–682, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH - T. Özis and V. Gülkaç, “Application of variable interchange method for solution of two-dimensional fusion problem with convective boundary conditions,”
*Numerical Heat Transfer; Part A: Applications*, vol. 44, no. 1, pp. 85–95, 2003. View at: Publisher Site | Google Scholar - V. Gülkaç and T. Özis, “Erratum to ‘On a lod method for solution of two dimensional fusion problem with convective boundary conditions’,”
*International Communications in Heat and Mass Transfer*, vol. 31, no. 4, pp. 597–606, 2004. View at: Google Scholar - V. Gülkaç, “On the finite differences schemes for the numerical solution of two-dimensional moving boundary problem,”
*Applied Mathematics and Computation*, vol. 168, no. 1, pp. 549–556, 2005. View at: Publisher Site | Google Scholar | MathSciNet - V. Gülkaç, “A numerical solution of the two-dimensional fusion problem with convective boundary conditions,”
*International Journal for Computational Methods in Engineering Science and Mechanics*, vol. 11, no. 1, pp. 20–26, 2010. View at: Publisher Site | Google Scholar | MathSciNet - R. M. Furzeland, “A survey of the formulation and solution of free and moving boundary (Stefan) problems,”
*Brunel University Technical Report*TR/76, Brunel University, Uxbridge, UK, 1977. View at: Google Scholar - J. Crank,
*Free and Moving Boundary Problems*, Oxford Science Publications, The Clarendon Press, Oxford, UK, 1984. View at: MathSciNet - V. R. Voller, “Fractional stefan problems,”
*International Journal of Heat and Mass Transfer*, vol. 74, pp. 269–277, 2014. View at: Publisher Site | Google Scholar - X. Wang, F. Liu, and X. Chen, “Novel second-order accurate implicit numerical methods for the Riesz space distributed-order advection-dispersion equations,”
*Advances in Mathematical Physics*, vol. 2015, Article ID 590435, 14 pages, 2015. View at: Publisher Site | Google Scholar | MathSciNet - V. R. Voller, “An exact solution of a limit case Stefan problem governed by a fractional diffusion equation,”
*International Journal of Heat and Mass Transfer*, vol. 53, no. 23-24, pp. 5622–5625, 2010. View at: Publisher Site | Google Scholar - Y. Shao and W. Ma, “Finite difference approximations for the two-side space-time fractional advection-diffusion equations,”
*Journal Computational Analysis and Applications*, vol. 21, no. 2, pp. 369–379, 2016. View at: Google Scholar - Y. Zhou and L.-J. Xia, “Exact solution for Stefan problem with general power-type latent heat using Kummer function,”
*International Journal of Heat and Mass Transfer*, vol. 84, pp. 114–118, 2015. View at: Publisher Site | Google Scholar - A. J. Arenas, G. Gonzalas-Parra, and B. M. Chen-Charpentier, “Construction of nonstandard finite difference schemes for the SI and SIR epidemic models of fractional order,”
*Mathematics and Computers in Simulation*, vol. 121, pp. 48–63, 2016. View at: Publisher Site | Google Scholar | MathSciNet - M. G. Mittag-Leffler, “Sur la nouvelle fonction E
*α*(x), Comptes Rendus,”*Acad. Sci. Paris*, vol. 137, pp. 554–558, 1903. View at: Google Scholar

#### Copyright

Copyright © 2016 Vildan Gülkaç. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.