New Challenges in Fractional SystemsView this Special Issue
An Implementation Solution for Fractional Partial Differential Equations
The link between fractional differentiation and diffusion equation is used in this paper to propose a solution for the implementation of fractional diffusion equations. These equations permit us to take into account species anomalous diffusion at electrochemical interfaces, thus permitting an accurate modeling of batteries, ultracapacitors, and fuel cells. However, fractional diffusion equations are not addressed in most commercial software dedicated to partial differential equations simulation. The proposed solution is evaluated in an example.
For an electrochemical system, species diffusion is often modeled by Fick’s laws . However, in almost all systems, the transport mechanism is different from the classical diffusion characterized by Fick’s laws. This kind of diffusion is denoted by anomalous diffusion . Anomalous diffusion is characterized by a mean squared displacement of the diffusing particles that has a power law dependence on time with between 0 and 2 (for classical diffusion ).
The theoretical approach of this type of diffusion is strongly related to fractional calculus . It was indeed demonstrated that stochastic processes of random walks can be represented by fractional diffusion equations .
Among all existing anomalous diffusion equations, three diffusion modes, respectively, called “anomalous diffusion Ia” (ADIA), “anomalous diffusion Ib” (ADIB), and “anomalous diffusion II” (ADII) are characterized by a Fick’s equations adaptation for phenomena considered .
In this paper, the authors have only considered ADIB type diffusion equations but the proposed contribution can be extended to others classes of equation. In the sequel, the link between fractional differentiation and diffusion is used to propose a solution for the implementation of a fractional diffusion equation in software such as COMSOL Multiphysics. These software applications are now powerful tools for engineers to simulate complex systems combining several physical domains such as electrochemistry and thermal. However they are not adapted to take into account anomalous diffusion and thus to model diffusion interfaces as in batteries, ultracapacitors, or fuel cells.
The link between fractional differentiation and diffusion equation is reminded in the second section of the paper. This link should be used to implement fractional differentiation in software dedicated to numerical solving of partial differential equation such as COMSOL Multiphysics software. However, as shown in Section 2, the diffusion equation form of a fractional system requires the computation of an inverse Fourier transform that is in most cases impossible to get analytically. This is why this paper proposes alternative partial differential equations approximation that exhibits a fractional behavior in a given frequency band. These differential equations can be easily implemented to simulate a fractional differentiator and thus a fractional diffusion equation.
2. Link between Fractional Systems and Partial Differential Equations
2.1. Partial Differential Equation Representation and Approximation of a Fractional System
For presentation simplicity, the following fractional system (fractional integrator) is considered with . Its link with diffusion equation can be demonstrated using the system impulse response  defined by the Mellin-Fourier integral of (1): where is greater than the abscissa of the singular points of . Using poles definition that can be found in [6, 7], this system does not generate poles and its impulse response is thus given by
Let the following state space representation can be obtained for system (1):
Such a representation can be generalised to a large class of fractional systems as demonstrated in [8, 9]. In these works, second relation in (7) is rewritten as and representations (7) and (8) are denoted by diffusive representation. For a fractional transfer function defined by with and where , function is defined by 
Relation (10) shows that a fractional integrator can thus be seen as an infinite dimensional system described by a diffusion equation. This interpretation is represented by Figure 1  in which(i)the input applied at the abscissa ,(ii)the real distributed state ,(iii)the output (weighed sum on the state) appears. This remark can be generalized to a large number of fractional systems and thus demonstrates their link with diffusion equations.
Implementation of relation (10) requires the integral truncation that can be done as follows:
In (13) where denotes the inverse Fourier transform. This relation is in practice impossible to compute analytically in most cases. To solve this problem, another partial differential equation is now proposed.
2.2. Another Partial Differential Equation Approximation
Using change of variable , relation (15) becomes
Implementation of such a transfer function requires the integral truncation, namely,
Note that , are homogenous to the logarithm of a frequency. Now Let be a function of the space variable of finite dimension () and of the time variable . This function satisfies the class of partial differential equations Also, let the system output be given by
This partial differential equation class has been studied in . Transfer function that links the system input and output is defined by Now if
It turns out that the solution of the partial differential equation defined by (15) and (16) is an approximation of the response of the system whose transfer is defined by (14) (and thus by (1) for a fractional integrator). Such a representation can thus be used to approximate a fractional system or a fractional integrator. It is now used to implement a fractional diffusion equation.
3. Application to Fractional Partial Differential Equations Implementation
Implementation problem of a fractional partial differential equation using simulation software such as Comsol Multiphysics is now addressed and the following class of equation is considered: with the following boundary and initial conditions:
It is an ADIb type anomalous diffusion equation that can also be rewritten as with According to the demonstration in Section 2 and in , fractional derivative of can be approximated by where functions and are, respectively, defined by
Implementation of (24) using approximations from (26) to (29) thus requires 2 geometries. In a first 1D-type geometry, (23) is implemented. To take into account relation (24), a second 2D-type geometry is created.
Values and for the 2D geometry along the axis are used to define the range of frequency for which the approximation of the fractional differentiation is expected using the relations Note that frequencies and depend, respectively, on the simulation duration and the sampling time . These frequencies can be defined using the following rules: and .
Information produced in the two geometries are then exchanged as described in Figure 4.
The following diffusion system with is considered: with the following initial and boundaries conditions: For , m/s, and As shown in the appendix, system (31) and (32) solution is defined by For s and s, this solution is represented by Figure 5.
The analytical solution is compared with the results produced by COMSOL Multiphysics for function using the implementation and geometries described in Figures 2–4. For the implementation, rd/s and rd/s. Figure 6 is a representation of the absolute error between function and the solution computed with COMSOL Multiphysics. This really small error permits us to validate the method we used for the implementation of a fractional diffusion equation using a partial differential equation for the approximation of a fractional derivative.
This paper proposes a method for the implementation of a fractional diffusion equation into simulation softwares such as COMSOL Multiphysics. These software applications are now powerful tools for engineers to simulate complex systems combining several physical domains such as electrochemical and thermal. However they are not adapted to take into account anomalous diffusion and thus to model diffusion interfaces as in batteries, ultracapacitors, or fuel cells. To permit the implementation, the link between fractional systems and diffusion equation is used. The fractional diffusion equation considered is splitted into two parts and the remaining fractional equation is approximated by a partial differential equation. For the implementation of this partial differential equation, an additional geometry is created (a 1D system is transformed into a 2D system). The efficiency of the proposed method is evaluated in an example. The results obtained showed the efficiency of the proposed method.
This appendix demonstrates how, using material provided in , the analytical solution of the following fractional diffusion equation: with the following initial and boundaries conditions is obtained. To obtain homogenous conditions at and , the following change of variable is used: with System thus becomes with Separation variable method leads to writing as: with Now let , , and Solution of system is given, according to , , and with Combining relations and leads to and thus Laplace transform applied to leads to and thus If is written as relation becomes Within the interval , and thus using .
If , function is a ramp and thus Inverse Laplace transform then permits Within the interval , , and thus .
Conflict of Interests
The authors declare no conflict of interests.
R. B. Bird, W. E. Stewart, and E. N. Lightfoot, Transport Phenomena, John Wiley & Sons, New York, NY, USA, 2nd edition, 2002.
A. Compte, “Stochastic foundations of fractional dynamics,” Physical Review E, vol. 53, no. 4, pp. 4191–4193, 1996.View at: Google Scholar
A. Oustaloup, Systèmes Asservis Linéaires D’ordre Fractionnaire, Masson, Paris, France, 1983.
J. Sabatier, M. Merveillaut, R. Malti, and A. Oustaloup, “How to impose physically coherent initial conditions to a fractional system?” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 5, pp. 1318–1326, 2010.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet