Advances in Mathematical Physics

Volume 2013, Article ID 479634, 8 pages

http://dx.doi.org/10.1155/2013/479634

## A Fractional Anomalous Diffusion Model and Numerical Simulation for Sodium Ion Transport in the Intestinal Wall

School of Mathematics, Shandong University, Jinan 250100, China

Received 17 May 2013; Accepted 1 July 2013

Academic Editor: Changpin Li

Copyright © 2013 Bo Yu and Xiaoyun Jiang. 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

The authors present a fractional anomalous diffusion model to describe the uptake of sodium ions across the epithelium of gastrointestinal mucosa and their subsequent diffusion in the underlying blood capillaries using fractional Fick’s law. A heterogeneous two-phase model of the gastrointestinal mucosa is considered, consisting of a continuous extracellular phase and a dispersed cellular phase. The main mode of uptake is considered to be a fractional anomalous diffusion under concentration gradient and potential gradient. Appropriate partial differential equations describing the variation with time of concentrations of sodium ions in both the two phases across the intestinal wall are obtained using Riemann-Liouville space-fractional derivative and are solved by finite difference methods. The concentrations of sodium ions in the interstitial space and in the cells have been studied as a function of time, and the mean concentration of sodium ions available for absorption by the blood capillaries has also been studied. Finally, numerical results are presented graphically for various values of different parameters. This study demonstrates that fractional anomalous diffusion model is appropriate for describing the uptake of sodium ions across the epithelium of gastrointestinal mucosa.

#### 1. Introduction

The intestinal wall represents a complex system which allows the passage of substances either through the cells or in between the cells. The luminal surface of the intestine is covered with a typically leaky epithelium which enables the passage of ions via the intercellular route. The substance to be absorbed either penetrates into the intercellular space directly through the tight junction or enters the cell cytoplasm through the apical plasma membrane from the lumen of the intestine and then penetrates through the lateral plasma membrane to enter the intercellular space. The latter route leads to the underlying lamina propria, which consists of connective tissue, blood vessels, and lymph capillaries, and thus the substance enters the circulation (Figure 1) [1]. The process in which the ions enter the cell is passive diffusion under concentration gradient and potential gradient. This is mainly because transmural electrical potential differences of 5–12 mV have been reported from a variety of species during recent years [2, 3]. Although the potential differences across the intestinal wall are relatively small, they cannot be ignored in the studies of the intestinal transport of charged species [4].

Numerous techniques involving both in vivo and in vitro preparations have been employed in the study of intestinal transport. But because the cells are too small to provide continuous sections large enough for steady-state determinations of their transmission properties in actual physical situations, the distribution of the ions in the cellular and extracellular phases cannot be determined experimentally. Therefore, the idea of analysing such physiological problems using a theoretical approach has arisen. Fadali et al. [5] proposed an analytical model for water absorption in the intestine based on an integration of mass balance equation for active contact area for absorption. The model gave a solution for the amount of water absorbed in the intestine as a function of time following water ingestion using data from the physiological literature. Hills [6] proposed a two-phase model to study linear bulk diffusion into a continuous fluid in which a less permeable phase was distributed as particles of irregular profile; the overall uptake of solute by a parallel-faced section of tissue could be expressed as the sum of an infinite number of exponential terms. The advantage of this model was that it represented the histology of cellular tissue in the most realistic manner. Karmakar and Jayaraman [1] presented a linear diffusion model of the intestinal wall to describe the uptake of lead ions across the epithelium of gastrointestinal mucosa and their subsequent diffusion in the underlying blood capillaries. The model studied the variation of concentration with time in the extracellular phase and the cellular phase and the mean concentration available for absorption by the blood capillaries as a function of time, and it also reported the determination of membrane permeability for lead through theoretical analysis. Varadharajan and Jayaraman [7] presented a theoretical approach to study the uptake of sodium ions across the gastrointestinal mucosa and the concentrations at which they were taken up into the underlying blood capillaries. The model took into account both the diffusion under concentration gradient and potential gradient and active transport which was ATPase enzyme mediated, and appropriate partial differential equations for the two mechanisms of transport had been derived and were solved by iterative methods.

Recently, fractional calculus has been a subject of worldwide attention due to its surprisingly broad range of applications in physics, chemistry, engineering, economics, biology, and so forth [8–11]. In particular, fractional calculus is a key tool to study anomalous diffusion in transport processes which implies a fractional Fick’s law for the flux that accounts for spatial and temporal nonlocality. Many literatures have shown that the power-law behavior is a hallmark of many biological phenomena observed at different scales and at various levels of organization and that a rheological behavior that conforms to the power law can be described by using methods of fractional calculus [12–16]. For example, Magin et al. [12] describe the formulation of the bioheat transfer in one dimension in terms of the fractional order differentiation with respect to time. His study demonstrates that fractional calculus can provide a unified approach to examine the periodic heat transfer in peripheral tissue regions. In this paper, according to Magin’s idea, based on the previous analysis and fractional calculus theory, we consider a fractional anomalous diffusion model to describe the uptake of sodium ions across the epithelium of gastrointestinal mucosa using fractional Fick’s law. In Section 2, we present the fractional anomalous diffusion model, and appropriate partial differential equations describing the variation with time of concentrations of sodium ions in both the interstitial phase and the intracellular phase across the intestinal wall are obtained using Riemann-Liouville space-fractional derivative and are solved by numerical computation. In Section 3, numerical results are presented graphically for various values of different parameters. In Section 4, we have presented our conclusions.

#### 2. Materials and Methods

The fractional anomalous diffusion model considers a two-phase structure of the intestinal wall in which the epithelium is treated as a thin layer. The apical plasma membrane is adjacent to the lumen of the intestine and at the origin of a one-dimensional coordinate system. The rest of the cellular elements form a uniformly distributed array of identical cells. In between are the intercellular spaces which correspond to interstitial phase (Figure 1).

##### 2.1. Fractional Fick’s Law

Fick’s law is extensively adopted as a model for standard diffusion processes. For example, the simplest reaction diffusion model in spherical coordinates can be expressed as where is the concentration of solute (with radial symmetry), represents reaction kinetics, and is dispersive flux. Generally, Fick’s law is used in normal diffusion for dispersive flux based on empirical observations:

where is the diffusion coefficient.

However, requiring separation of scales, it is not suitable for describing nonlocal transport process. In order to study the anomalous diffusion, the fractional Fick’s law has been proposed [17], where the gradient of the solute concentration in the empirical flux equation is replaced by a fractional-order derivative:

where , , and is the anomalous diffusion coefficient. and are Riemann-Liouville operators which are defined as follows: where , . We name it as the time-space fractional Fick’s law [17]. Some special cases of this equation are as follows when , , it gives the classical Fick’s law; when , it gives time fractional Fick’s law; when , it gives space fractional Fick’s law. Here, we only consider the case of , that is, the space fractional Fick’s law.

##### 2.2. Fractional Anomalous Diffusion Model

The diffusion of sodium ions is complicated because its flux is determined by both the concentration gradient and the electrical gradient. Considering the motion of sodium ions under all forces, Macey [18] proposes that the flux equation is written as follows:

where is the diffusion coefficient, is the electrical potential, is the concentration of the sodium ions, is the distance across the wall measured from the lumen, is the charge on the ion ( for the sodium ion), is Faraday’s constant (), is the universal gas constant, and is the absolute temperature.

Here, according to Magin’s idea [12], based on the space fractional Fick’s law, the flux equation is expressed in the following form:

where , is the anomalous diffusion coefficient. The first term on the right stands for the concentration gradient, and the second term on the right stands for the electrical gradient.

We consider a two-phase model consisting of the interstitial phase and the intracellular phase. The mass balance equation in the interstitial phase, which accounts for the molecular diffusion flux and a uniformly distributed continuum of point sinks whose strength is proportional to the local concentration differences between the two phases [1], is where and are the concentrations of sodium ions in the interstitial phase and in the intracellular phase, respectively, and is the membrane permeability coefficient for the molecular diffusion of sodium ions into the cellular phase. Substituting (6) into (7), we can get the following equation:

And based on the assumption that diffusion does not contribute significantly to the total molecular transport inside the cell [1], the mass balance equation in the cellular phase is

which is justified by the fact that the dimensions of the cells are small compared to the thickness of the intestinal wall; therefore, the flux through them is independent of distance.

Meanwhile, we assume that , where is a constant to be determined. A justification for the constant field assumption can be found in the observation that if a membrane contains a large number of dipolar ions close to their isoelectric point, these dipoles will tend to alter their orientation in such a way that they tend to smooth out any irregularities and maintain a constant field [7]. The validity of this assumption has also been discussed by Goldman [19] and Cole [20]. Hence, (8) can be reduced to that is,

where .

The value corresponds to the lumen of the intestine and corresponds to serosa. We are interested in finding the ion concentration at , which corresponds to the blood capillary at which it is absorbed (Figure 1). In rats, the mucosal epithelium is approximately of the total intestinal wall thickness. Equations (9) and (11) are solved to obtain and as functions of and , and the mean concentration of sodium ions at is calculated from

where and are the interstitial and intracellular volume fractions, respectively.

Then, we introduce dimensionless parameters

to reduce (11) and (9) to the nondimensional form (with the notation dropped for convenience):

where is the concentration of sodium ions in the lumen. The parameter could be considered as the ratio of the membrane diffusion flux into the cellular phase to the molecular diffusion flux in the interstitial phase.

Based on the assumption that the concentration of sodium ions in the intestinal lumina surface is equal to the concentration at its abluminal surface for the epithelial is specially thin, the boundary conditions are given by

which mean that the concentration of sodium ions in the lumen is set to 1, whereas at the serosa it is set to 0 at all times. Further, the initial concentration is taken to be 0, a condition justified in the case of in vitro experiments. Both and are set to 0 at any point within the tissue when time equals 0. Therefore, the initial conditions are given by

##### 2.3. Numerical Computation

For the numerical solution of the problem above, we introduce a uniform grid of mesh points , with , , and , , where and are two positive integers, and are the uniform spatial and temporal mesh size, respectively. The theoretical solution at the point is denoted by ; the solution of an approximating difference scheme at the point will be denoted by . Similarly, the theoretical solution at the point is denoted by ; the solution of an approximating difference scheme at the point will be denoted by .

Then, we start to introduce the discretization of the differential operators. The first-order derivatives with respect to the temporal variable and are approximated by the following Euler backward difference, respectively:

and the first-order derivative with respect to the spatial variable is approximated by Euler forward difference:

As for the Riemann-Liouville fractional derivative, using the relationship between the Grünwald-Letnikov formula and Riemann-Liouville fractional derivative, we can approximate the fractional derivative by [21, 22]

where , for . There, we have adopted the shifted Grünwald-Letnikov formula for .

Finally, the finite difference method for the above problem is given as follows:

The boundary and initial conditions can be discretized by

The concentrations of the sodium ions in the intercellular phase and intracellular phase are determined at different steps of time and space, and their weighted mean concentration at the blood capillaries can also be obtained.

#### 3. Results and Discussion

The thickness of the intestinal wall is taken to be [1], which is measured from a cross-section of the rat intestinal wall using an ocular micrometer fitted to a simple microscope. Muller [23] reported in morphometric studies of rat gastric mucosa that the epithelial cells occupied while the remaining was occupied by lamina propria. Hence, a choice of is made for and for , arbitrarily as their values are not available in the literatures. According to Xu and Zhao [24], the permeability coefficient of sodium ions is taken to be . Varadharajan and Jayaraman [7] had studied that the numerical value of depended on the potential difference between the serosa and the mucosa, and, comparing with the experimental results of Lauterbach [25], Varadharajan obtained the value of to be and the optimum value of to be around 0.4, so is taken to be 0.4 in our studies.

##### 3.1. The Anomalous Diffusion Coefficient of Sodium Ions

According to the Stokes-Einstein formula , where is Boltzmann’s constant (), is the temperature (~310 K), is the viscosity of intercellular fluid (~0.001 Pas), and is the radius of the water molecule (~0.45 nm), we can obtain the diffusivity in water [1]. Nevertheless, it can also be considered as a reasonable approximation for the diffusion coefficient of sodium ions. Here, we change this parameter to , , , , and times of the value of diffusivity .

Figure 2, individual graph, demonstrates the variation of concentration with respect to for different values of when . During the course of diffusion, we find that the concentrations decrease in an exponential way, and the curves tend to change very little when the values of vary not too much, which indicates that the distribution of sodium ions is more uniform; therefore, we can consider any of these values as a reasonable approximation for the diffusion coefficient of sodium ions. By amplifying the figures, we also observe that smaller increases the amplitude of , which indicates that larger increases the speed of the absorption especially at the blood capillaries, a possible explanation is that faster movement of sodium ions makes it easier for diffusion.

##### 3.2. The Order of Fractional Derivative

Here, based on the above analysis, we take to be , which is about 0.7 times of the diffusivity , just as used in the literature [18]. Substituting the respective values of and , we can obtain that the actual time is about .

Figure 3 demonstrates that the variation of concentration with respect to at min for different values of . Obviously, we can observe that the concentrations decrease in an exponential way and the curves become smoother as the value of becomes larger, which indicate that the diffusion in the intercellular phase and the absorption especially at the blood capillaries become quicker as the value of becomes smaller. This fact demonstrates that the diffusion of sodium ions is anomalous superdiffusion.

Figure 4 demonstrates that the variation of concentration with respect to at different times when , . During the course of diffusion, we find that the concentrations decrease in an exponential way. By amplifying the figures, we obviously observe that the decay tends to be smoother and smoother when time increases, which indicates the distribution of sodium ions is more uniform. Figure 4(a) shows that at earlier times most of the sodium ions are absorbed by the cells, and for later times they tend to pass towards the serosa. Meanwhile, we observe that most of the obsorption takes place at from Figure 4(b), which can be explained as the distance at which the blood capillaries lie. This is quite reasonable, since in rats the mucosal epithelium is about of the total wall thickness [1]. All these phenomena are connective with the results of Varadharajan and Jayaraman [7], which indicate that the anomalous diffusion is appropriate for describing the uptake of sodium ions.

Figure 5(a) is plotted against at different when , . We find that, at a particular distance, increases with time, but at a farther distance, the flux is lower because the major absorption is made available at the distance where the blood capillaries lie. Figure 5(b) is plotted against at for different values of . We find that, at a particular distance , is lower when is smaller, which indicates that the diffusion is quicker when the value of is smaller, and this leads to the lower value of .

Figure 6(a) is , plotted against at for different values of . We find that, at a particular distance , is lower when is smaller, which indicates that the diffusion is quicker when the value of is smaller, and this leads to the lower value of . Figure 6(b) is plotted against at different when , . We find that about sodium ion absorption is achieved at a distance of 0.1–0.15 when we choose , and it is in good agreement with the experimental results of Lauterbach [25], which indicate that the concentration of Na in the cell water approaches of the initial concentration of the incubation medium after the addition to the luminal side.

#### 4. Conclusions

In summary, in this paper we have derived a fractional anomalous diffusion model for sodium ion transport in the intestinal wall using space fractional Fick’s law. Appropriate partial differential equations describing the variation with time of concentrations of sodium ions in both the interstitial phase and the intracellular phase across the intestinal wall are obtained using Riemann-Liouville space-fractional derivatives and are solved by finite difference methods. The numerical simulations have been discussed, and numerical results are presented graphically for various values of different parameters. It demonstrates that fractional anomalous diffusion model is appropriate for describing the uptake of sodium ions across the epithelium of gastrointestinal mucosa. This research also provides some new points for studying ions transferring processes in biological systems.

#### Acknowledgment

The project was supported by the National Natural Science Foundation of China (Grants nos. 11072134, 11102102, and 91130017).

#### References

- N. Karmakar and G. Jayaraman, “Linear diffusion of lead in the intestinal wall: A Theoretical Study,”
*IMA Journal of Mathematics Applied in Medicine and Biology*, vol. 5, no. 1, pp. 33–43, 1988. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - W. M. Armstrong, “Cellular mechanisms of ion transport in the small intestine,”
*Physiology of the Gastrointestinal Tract 2*, pp. 1251–1266, 1987. View at Google Scholar - S. G. Schultz, “Principles of electrophysiology and their application to epithelial tissues, gastrointestinal physiology,” in
*Gastrointestinal Physiology*, pp. 87–91, Butterworth, London, UK, 1974. View at Google Scholar - S. G. Schultz and P. F. Curran,
*Handbook of Physiology*, vol. 3, Section 6, 1968. - M. S. Fadali, J. W. Steadman, and R. G. Jacquot, “A model of water absorption in human intestine,”
*Biomedical Sciences Instrumentation*, vol. 16, pp. 49–53, 1980. View at Google Scholar · View at Scopus - B. A. Hills, “Linear bulk diffusion into heterogeneous tissue,”
*The Bulletin of Mathematical Biophysics*, vol. 30, no. 1, pp. 47–59, 1968. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - M. Varadharajan and G. Jayaraman, “Sodium ion transport in the intestinal wall: a mathematical model,”
*IMA Journal of Mathematical Medicine and Biology*, vol. 11, no. 3, pp. 193–205, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - D. Baleanu, A. K. Golmankhaneh, R. Nigmatullin, and A. K. Golmankhaneh, “Fractional Newtonian mechanics,”
*Central European Journal of Physics*, vol. 8, no. 1, pp. 120–125, 2010. View at Publisher · View at Google Scholar · View at Scopus - J. Singh, P. K. Gupta, and K. N. Rai, “Solution of fractional bioheat equations by finite difference method and HPM,”
*Mathematical and Computer Modelling*, vol. 54, no. 9-10, pp. 2316–2325, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Metzler and J. Klafter, “The random walk's guide to anomalous diffusion: a fractional dynamics approach,”
*Physics Reports*, vol. 339, no. 1, pp. 1–77, 2000. View at Publisher · View at Google Scholar · View at MathSciNet - R. Metzler and J. Klafter, “The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics,”
*Journal of Physics A*, vol. 37, no. 31, pp. R161–R208, 2004. View at Publisher · View at Google Scholar · View at MathSciNet - R. L. Magin, S. C. Boregowda, and C. Deodhar, “Modeling of pulsating peripheral bio-Heat transferusing fractional calculus and constructal theory,”
*International Journal of Design, Nature, and Ecodynamics*, vol. 1, no. 18, pp. 18–33, 2007. View at Google Scholar - X. Y. Jiang and H. T. Qi, “Thermal wave model of bioheat transfer with modified Riemann-Liouville fractional derivative,”
*Journal of Physics*, vol. 45, no. 48, 10 pages, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - W. C. Tan, C. Fu, C. Fu, W. Xie, and H. Cheng, “An anomalous subdiffusion model for calcium spark in cardiac myocytes,”
*Applied Physics Letters*, vol. 91, no. 18, 3 pages, 2007. View at Publisher · View at Google Scholar · View at Scopus - R. G. Larson,
*The Structure and Rheology of Complex Fluids*, vol. 2, Oxford University Press, New York, NY, USA, 1999. - R. L. Magin, “Fractional calculus models of complex dynamics in biological tissues,”
*Computers & Mathematics with Applications*, vol. 59, no. 5, pp. 1586–1593, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X. Y. Jiang, M. Y. Xu, and H. T. Qi, “The fractional diffusion model with an absorption term and modified Fick's law for non-local transport processes,”
*Nonlinear Analysis: Real World Applications*, vol. 11, no. 1, pp. 262–269, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. I. Macey,
*Membrane Physiology*, Plenum, New York, NY, USA, 2nd edition, 1980. - D. E. Goldman, “Potential, impedence, and rectification in membranes,”
*The Journal of General Physiology*, vol. 27, no. 1, pp. 37–60, 1943. View at Publisher · View at Google Scholar - K. S. Cole, “Electrodiffusion models for the membrane of squid giant axon,”
*Physiological Reviews*, vol. 45, pp. 340–379, 1965. View at Google Scholar · View at Scopus - M. M. Meerschaert and C. Tadjeran, “Finite difference approximations for fractional advection-dispersion flow equations,”
*Journal of Computational and Applied Mathematics*, vol. 172, no. 1, pp. 65–77, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. Li and F. Zeng, “Finite difference methods for fractional differential equations,”
*International Journal of Bifurcation and Chaos in Applied Sciences and Engineering*, vol. 22, no. 4, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - O. Muller, “The quantitation of the rat gastric mucosa by morphometric methods,”
*Scandinavian Journal of Gastroenterology, Supplement*, vol. 101, no. 1, pp. 1–6, 1984. View at Google Scholar · View at Scopus - M. Y. Xu and G. L. Zhao, “Nonisotonic reabsorption of water and steady state distribution of sodiumion in renal descending limb of Henle,”
*Acta Mechanica Sinica*, vol. 1, no. 3, 1986. View at Google Scholar - F. O. Lauterbach,
*Intestinal Transport*, University Park Press, 1976.