Mathematical Problems in Engineering

Volume 2016 (2016), Article ID 8409839, 7 pages

http://dx.doi.org/10.1155/2016/8409839

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

Department of Mathematics, Faculty of Arts and Science, Kocaeli University, Umuttepe Campus, Izmit, 41380 Kocaeli, Turkey

Received 9 March 2016; Accepted 22 May 2016

Academic Editor: Evangelos J. Sapountzakis

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.

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