MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 807034 10.1155/2013/807034 807034 Research Article Modeling of an EDLC with Fractional Transfer Functions Using Mittag-Leffler Equations Quintana J. J. 1 Ramos A. 2 Nuez I. 1 Ahmad Bashir 1 Electronic and Automatic Engineering Department University of Las Palmas de G.C. 35017 Las Palmas Spain ulpgc.es 2 Process Engineering Department University of Las Palmas de G.C. 35017 Las Palmas Spain ulpgc.es 2013 3 12 2013 2013 31 07 2013 30 10 2013 31 10 2013 2013 Copyright © 2013 J. J. Quintana et al. 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.

Electrochemical double-layer capacitors (EDLC), also known as supercapacitors or ultracapacitors, are devices in which diffusion phenomena play an important role. For this reason, their modeling using integer-order differential equations does not yield satisfactory results. The higher the temporal intervals are, the more problems and errors there will be when using integer-order differential equations. In this paper, a simple model of a real capacitor formed by an ideal capacitor and two parasitic resistors, one in series and the second in parallel, is used. The proposed model is based on the ideal capacitor, adding a fractional behavior to its capacity. The transfer function obtained is simple but contains elements in fractional derivatives, which makes its resolution in the time domain difficult. The temporal response has been obtained through the Mittag-Leffler equations being adapted to any EDLC input signal. Different charge and discharge signals have been tested on the EDLC allowing modeling of this device in the charge, rest, and discharge stages. The obtained parameters are few but identify with high precision the charge, rest, and discharge processes in these devices.

1. Introduction

In recent years, the growing demand for new electrical energy storing systems has led to a remarkable development of electrochemical double-layer capacitors (EDLC). The EDLC are devices capable of storing energy and are characterized by their very rapid response during charge and discharge cycles, which allows them to provide high power and to hold a high number of charge and discharge cycles.

The complementary qualities of EDLC and batteries have allowed the generation of numerous hybrid applications for energy recovery or storage systems. Energy storage in EDLC is not supported by chemical processes. Moreover, supercapacitors have a high life cycle and do not need any maintenance. From the above it can be concluded that because of their characteristics these devices have raised great interesting expectations [1, 2].

It is a usual procedure that the dynamics of real systems are modeled by differential equations. In most cases, the differential equations are based on conventional derivatives, yielding sufficiently accurate mathematical models. However, there are a variety of systems or phenomena in which mathematical models based on ordinary differential equations do not provide satisfactory solutions. One of these elements is the EDLC. It is this case, when the application of models based on fractional derivatives, studied in the field of fractional calculus, is very useful. These equations are called fractional differential equations.

In multiple fields of physics, there are numerous examples in which a differential equation is used as fractional modeling tool . This is very common when considering geometry conditions with fractal dimension and distributed parameter systems . In the field of electric power, the current through a capacitor is proportional to the noninteger-order integral for electric current [9, 10]. An electric network composed of infinite RC elements can be modeled through fractional differential equations [11, 12].

In this paper, a dynamic fractional model for EDLC is developed, based on a differential equation with fractional derivatives. This equation has been solved by means of the Laplace transform, and in order to obtain the time domain solution used the Mittag-Leffler function Eα,β(z). The solution has been obtained for a unitary step input and for a generic input. This model has been applied to an EDLC for its identification in long time periods, obtaining very satisfactory results.

This paper is organized as follows: Section 2 shows a brief review of fractional calculus; in Section 3 the fractional model of the EDLC is deduced; Section 4 expounds the experimental data and the discussion; finally in Section 5 the conclusions are presented.

2. Brief Review of Fractional Calculus

In this section is exposed a brief mathematical background of the fractional calculus , for understanding methods, results, and conclusions presented in this paper. The following are some of the most popular functions, definitions, and properties of the fractional calculus.

2.1. Riemann-Liouville Derivative Definition

Riemann-Liouville definition expresses the fractional derivative as a time convolution integral (1)Datαf(t)=1Γ(n-α)dndtnatf(τ)(t-τ)α-n+1dτ(n-1<α<n), where  n  is an integer and  α  is a real number and the fractional derivative order.

2.2. Caputo Derivative Definition

Caputo definition of a fractional derivative of a function is (2)DaCtαf(t)=1Γ(n-α)atf(n)(τ)(t-τ)α+1-ndτ(n-1α<n), where n is an integer and α is a real number and the fractional derivative order.

2.3. Grunwald-Letnikov Derivative Definition

Grunwald-Letnikov definition is a numerical form of the fractional derivative (3)Datαf(t)=limh01hαj=0[(t-a)/h](-1)j(αj)f(t-jh).

2.4. Laplace Transform of the Fractional Derivative

An interesting property of the fractional derivative operator is its Laplace transform (4)L[D0tαf(t)]=0e-stDatαf(t)dt=sαF(s)-k=0n-1skDatα-k-1f(t)|t=0(n-1<αn), where, if initial conditions are null, it yields (5)L[D0tαf(t)]=sαF(s).

2.5. Mittag-Leffler Function

The two-parameter function of the Mittag-Leffler type Eα,β(z) plays a very important role in the fractional calculus and is defined as (6)Eα,β(z)=k=0zkΓ(αk+β)(α>0,β>0). The exponential function ez is a particular case of the Mittag-Leffler function (7)E1,1(z)=k=0zkΓ(k+1)=ez. One of the main uses of the Mittag-Leffler function comes from the following Laplace transform (8)L[tαk+β-1Eα,β(k)(±atα)]=0e-sttαk+β-1Eα,β(k)(±atα)dt=k!sα-β(sαa)k+1, where  k  indicates the order of the derivative of the Mittag-Leffler function (9)Eα,β(k)(z)=dkdzkEα,β(z). Considering the specific case for  k=0  yields (10)L[tβ-1Eα,β(±atα)]=sα-βsαa. And by substituting the Mittag-Leffler function, the following is obtained: (11)L[tβ-1k=0(±atα)kΓ(αk+β)]=sα-βsαa. Which will be the relation used for calculating the solution in the time domain.

3. Dynamic Model of EDLC

Electric modeling of supercapacitors is an active line of investigation with numerous contributions in the last few years . There are not yet totally satisfactory models applicable to any operating mode of the EDLC, allowing the EDLC behavior to be observed and simulated. That is the reason why models based on heuristic techniques are mainly used . Models based on the physics processes occurring in EDLC are also used . Most of the models obtained make use of a considerable number of variables, either of a high number of passive elements or of fractional variables models . Either way, testing these devices during long time periods has allowed observing their behavior and their applicability to energy storage systems. Charge and discharge are rapid processes occurring in a short time interval, since high currents are used.

Due to diffusion processes resulting after the charging phase. The evolution of the processes occurring inside the EDLC generates a nonreversing thermal dissipation which makes the assessment of losses in the device possible. This process would occur in both the EDLC charge and discharge phases.

The model of EDLC has been considered as single input (current), single output (voltage) system. The electric model used comes from the model applied to ideal capacitors using integer-order differential equations. Subsequently, a model based on fractional order differential equations is proposed.

3.1. Electric Model of a Conventional Capacitor

The mathematical equation, relating the current and the voltage at the terminals of a real capacitor, can be deduced from an electric circuit made of an ideal capacitor and two parasitic resistors [9, 19].

Naming the series resistor as  R1  and the parallel one as  R2  and operating the following equation is obtained: (12)R2·Cdv(t)dt+v(t)=(R1+R2)·i(t)+R1·R2·Cdi(t)dt with  v(t)  being the voltage at the terminals of the EDLC and  i(t)  its current. The transfer function or impedance in the Laplace domain will be calculated as the quotient between voltage and current (13)G(s)=V(s)I(s)=R1(ss+(1/(R2·C)))+R1+R2R2·C(1s+(1/(R2·C))). Gathering the constants of the transfer function the following is obtained: (14)G(s)=V(s)I(s)=(k1ss+a+k21s+a), where (15)k1=R1,k2=R1+R2R2·C,a=1R2·C. The differential equation resulting from the transfer function can be solved for simple and complex signals using the existing simulation programs.

3.2. Electric Model of EDLC

The model obtained for real capacitor (12) fits well to the behavior of real capacitors in a wide range of frequencies. However, if applied to an EDLC, the adjustment is unsatisfactory, especially when long time periods are analyzed.

Many papers [3, 10, 20] exposed the fractional behavior of EDLC. Thus, by changing the ideal capacitor in the model in Figure 1 for a fractional one, the voltage and the current across such fractional capacitor are given by (16)C·Dαv(t)=i(t). Applying this relation to the circuit of Figure 1, it gives a differential equation of the form (17)D0tαv(t)+a·v(t)=k1·D0tαi(t)+k2·i(t). Applying Laplace transform (5) and assuming that the initial conditions are zero, the following is obtained: (18)G(s)=V(s)I(s)=(k1sαsα+a+k21sα+a). This transfer function is similar to (14), differing only in the fractional aspect of the capacitor. In this case differential equation (17) is of fractional order, being defined by variable  α. This equation can be solved for simple and complex signals input, using Mittag-Leffler functions [3, 21].

Equivalent circuit of a real capacitor.

3.3. Solution to the Fractional Equation of EDLC

In this section the solution to (18) will be deduced in the time domain for a variable input over time.

3.3.1. Solution to a Unitary Step Input

Multiplying transfer function (18) by the unit step input  i(s)=1/s, the response of the system to this signal is obtained: (19)v(s)=k1sα-1sα+a+k2s-1sα+av(t)=k1L-1[sα-1sα+a]+k2L-1[s-1sα+a]. Applying relation (11), (20)v(t)=k1j=0(-a·tα)jΓ(α·j+1)+k2·tα·j=0(-a·tα)jΓ(α·j+α+1). The response in the time domain to a unit step input is obtained.

3.3.2. Solution to a Generic Input

The response can be extended to any input. For this purpose, the superposition of steps displaced by the sampling period of the input signal has been used. Thus, for the first input signal, the temporal response matches the response obtained by (20). So, for the first time interval, (21)v(t)=i(0)·(k1j=0(-a·tα)jΓ(α·j+1)+k2·tα·j=0(-a·tα)jΓ(α·j+α+1)),0tT. As the input signal changes in the second period, its response can be obtained by superposing to the previous response a second step, delayed in time, a period  T, resulting in (22)v(t)=i(0)·(k1j=0(-a·tα)jΓ(α·j+1)+k2·tα·j=0(-a·tα)jΓ(α·j+α+1))+(i(T)-i(0))·(k1j=0(-a(t-T)α)jΓ(α·j+1)+k2(t-T)αj=0(-a(t-T)α)jΓ(α·j+α+1))Tt2T. Generalizing (22) for a generic time interval  N, the following is obtained: (23)v(t)=k=0N(i(kT)-i((k-1)T))·[k1j=0(-a·(t-kT)α)jΓ(α·j+1)+k2(t-kT)α·j=0(-a·(t-kT)α)jΓ(α·j+α+1)],NTt(N+1)T, where the initial conditions are zero.

4. Experimental Results, Identification, and Discussion 4.1. Instruments

The experimental data have been obtained in the laboratory making a test circuit and recording the electrical variables of the EDLC. The analyzed EDLC, manufactured by ELNA, has a capacity of 4.7 F and a voltage of 2.5 V.

Figure 2 shows the electric circuit for the charge-discharge control. A Crouzet Millenium III PLC has controlled the charge and discharge switches. Voltage and current data have been recorded through a NI USB-6009 data logger by National Instruments.

Schematic diagram of the experimental circuit.

4.2. Experimental Results

Several tests have been performed to the EDLC in order to monitor voltage and current signals. Because EDLC have memory effect, just before each test they have been short-circuited for 24 hours to ensure zero initial conditions.

The tests were performed as follows. In a first phase lasting a few minutes, was charged at a constant current between 20 mA and 0.9 A until the voltage reached 2.4 V; in the second phase, the current is cut off and left to rest for 8 hours to allow the voltage stabilization, and finally in the third phase is discharged through a resistance of 5 ohm until it was fully discharged.

Although the sampling period of voltage and current has been recorded at 0.1 s; sfor identification purposes it has been considered a 5 s sampling time. Figure 3 shows the experimental data obtained. The evolution of voltage in the EDLC throughout the process can be appreciated in Figure 3(a), while the current supplied at Figure 3(b). They are observed perfectly the three stages described above.

Charge-rest-discharge at 20 mA.

4.3. Identification

The coefficients for the proposed model have been defined according to transfer function (18) and implemented using (23).

For identification the MATLAB simulation software was used. The function created has as input parameters  α, a,  k1, and  k2. Using (23) the voltage response is generated taking all the current samples as input. The voltage is compared with the experimental data through standard deviation (24), generating the output of the function. Using this function in the command fminsearch of MATLAB, sought parameters are obtained. For implementing (23), the terms  j  of the sum may be truncated without causing significant errors.

The index chosen in this paper is the standard deviation defined as (24)σd=i=1N(yi,exp-yi,cal)  2N-1, where  N  is the number of samples and subindexes which represent the values obtained experimentally and calculated according to the model.

The EDLC identification has been carried out calculating the parameters  α, a, k1, and  k2  in (18) and has consisted of two parts: first, from the beginning until just before discharged and second, from the discharge onwards. This is because the internal processes in the EDLC are different [22, 23].

The data obtained for the charge and rest phases are shown in Table 1.

Charge + rest identification parameters.

Test I (A) α a k 1 k 2 σ d
1 0.02 A 0.941 5.974 E - 08 0.239 0.202 2.079 E - 02
2 0.1 A 0.952 8.671 E - 08 0.261 0.187 1.297 E - 02
3 0.5 A 0.950 1.02 E - 07 0.267 0.170 1.354 E - 02
4 0.5 A 0.950 1.020 E - 07 0.250 0.169 1.338 E - 02
5 0.9 A 0.953 9.917 E - 08 0.228 0.164 1.096 E - 02

Considering relation (15) between the parameters calculated with the resistors and the capacity in Figure 1, the parameters obtained for the charge are in Table 2.

Charge + rest electrical parameters.

Test I (A) α C R 1 R 2
1 0.02 A 0.941 4.95 0.239 3.382 E + 06
2 0.1 A 0.952 5.35 0.261 2.155 E + 06
3 0.5 A 0.950 5.89 0.267 1.660 E + 06
4 0.5 A 0.950 5.91 0.250 1.660 E + 06
5 0.9 A 0.953 6.08 0.228 1.658 E + 06

In order to compare the adjustment with traditional model (14) in which the fractional index is an integer, Table 3 and Figure 4 show the best fit to the data of test number 5 setting  α=1.

Parameters calculated for test number 5 setting α=1.

Test α C R 1 R 2
1 1.000 9.00 1.262 2.191 E + 04

Fitting in the resting phase using integer derivatives.

Operating in a similar way, the parameters for the discharge through the 5 Ω resistor are in Table 4.

Discharge identification parameters.

Test α a k 1 k 2 σ d
1 0.975 9.600 E - 08 0.235 0.187 1.391 E - 02
2 0.981 9.282 E - 08 0.244 0.185 1.281 E - 02
3 0.978 8.427 E - 08 0.267 0.192 2.068 E - 02
4 0.976 1.556 E - 07 0.251 0.202 2.338 E - 02
5 0.982 1.032 E - 07 0.240 0.171 1.285 E - 02

Of which the following electric parameters are obtained in Table 5.

Discharge electrical parameters.

Test α C R 1 R 2
1 0.975 5.34 0.235 1.950 E + 06
2 0.981 5.40 0.244 1.996 E + 06
3 0.978 5.21 0.267 2.278 E + 06
4 0.976 4.95 0.251 1.297 E + 06
5 0.982 5.85 0.240 1.656 E + 06

The following figures show graphically the fit of the proposed model with experimental data. Figure 5 shows that the model has a good precision all over the time range tested. Figures 6 and 8 display the model and the experimental data in the charging and discharging phases. Figures 4 and 7 show the fit when using transfer functions with integer and fractional indices, showing that the fit with the traditional model gives unsatisfactory results. In summary, it can be observed that the voltage deduced through the fractional model adjusts very close to the real data even in long time periods.

Experimental data of test number 5 versus proposed model.

Detail of fitting in charging phase.

Fitting in the resting phase using the proposed model.

Fitting in the discharging phase.

4.4. Discussion

The self-discharge adjustment of EDLC after charging and after long time resting phases is a complex phenomenon which is not easily applicable to modeling through integer-order transfer functions. In this work the EDLC has been modeled making use of the fractional derivative. Taking the electric circuit in Figure 1 as starting point and assuming that the capacitor has a fractional index  α, fractional transfer function (18) has been deduced. This function has been solved in the time domain using Mittag-Leffler functions and implemented in a simulation program through (23).

First of all, considering Figure 5, it can be observed that the proposed transfer function adjusts well to the real data, modeling satisfactorily the self-discharge phenomenon.

Unlike the traditional capacitor model, which is defined by the value of two resistors and the capacity, the proposed model adds a new parameter, that is, the fractional index  α  of the fractional capacitor. This new parameter models to a greater extent the self-discharging phenomenon.

Although constants in Tables 1 and 4 have been taken for the identification of equations, constants in Tables 2 and 5 will provide the data related to the value of the electric elements.

In all cases the adjustment has been very precise, with the typical deviation being less than 0.03.

The constancy of  R1,R2, and  C  values for both charge and discharge is especially noteworthy.

At last, the most significant difference in terms of charge and discharge is given by the value of the fractional index  α, which is around 0.95 in the charge and around 0.98 in the discharge.

5. Conclusions

This work has consisted of EDLC modeling by implementing a simple fractional model allowing a very satisfactory response for long time periods. Traditional models based on transfer functions with integer coefficients yield good results in those cases when there is no self-discharge phenomenon. Applying fractional mathematics implies an additional parameter to the already necessary parameters in order to define the basic model shown in Figure 1, which is the fractional index  α. According to the tests this index has different values for the charge and the discharge, with those values being apparently independent of the EDLC charging current.

Thanks to the proposed transfer function and the Mittag-Leffler functions, it is possible to estimate the voltage at the EDLC terminals after a long time, taking into account the charging current and the way it has been charged.

Burke A. Ultracapacitors: why, how, and where is the technology Journal of Power Sources 2000 91 1 37 50 2-s2.0-0034324931 10.1016/S0378-7753(00)00485-7 Conway B. Electrochemical Supercapacitors. Scientific Fundamentals and Technological Applications 1999 New York, NY, USA Kluwer Academic Podlubny I. Fractional Differential Equations 1999 198 San Diego, Calif, USA Academic Press Mathematics in Science and Engineering MR1658022 Oldham K. B. Spanier J. The Fractional Calculus 1974 1047 New York, NY, USA Academic Press MR0361633 Miller K. S. Ross B. An Introduction to the Fractional Calculus and Fractional Differential Equations 1993 New York, NY, USA John Wiley & Sons MR1219954 Oustaloup A. Ballouk A. Bansard B. Partial differential equations and non integer derivation Proceedings of the International Conference on Systems, Man and Cybernetics October 1993 166 173 2-s2.0-0027850118 Mbodje B. Montseny G. Boundary fractional derivative control of the wave equation IEEE Transactions on Automatic Control 1995 40 2 378 382 10.1109/9.341815 MR1312917 ZBL0820.93034 Montseny G. Audounet J. Mbodje B. Optimal model of fractional integrators and application to systems with fading memory Proceedings of the International Conference on Systems, Man and Cybernetics October 1993 65 70 2-s2.0-0027851592 Westerlund S. Ekstam L. Capacitor theory IEEE Transactions on Dielectrics and Electrical Insulation 1994 1 5 826 839 2-s2.0-0028530783 10.1109/94.326654 Bano N. Hashmi R. A. Electrical studies of leaves over wide frequency range IEEE Transactions on Dielectrics and Electrical Insulation 1996 3 2 229 232 2-s2.0-0030128540 10.1109/94.486774 Carlson G. E. Halijak C. A. Simulation of the fractional derivate operator and the fractional integral operator s 45 Proceedings of the Central State Simulation Council Meeting May 1961 7 Kansas State University Haba T. C. Martos M. Ablart G. Bidan P. Composants électroniques à impédance fractionnaire 5 Proceedings of the ESAIM conferece on Fractional Differential Systems: Models, Methods and Applications 1998 99 109 Quintana J. J. Ramos A. Nuez I. Identificaction of the fractional impedance of ultracapacitors Proceedings of the IFAC Workshop (FDA’06) 2006 Oporto, Portugal Zubieta L. Bonert R. Characterization of double-layer capacitors for power electronics applications IEEE Transactions on Industry Applications 2000 36 1 199 205 2-s2.0-0033875953 10.1109/28.821816 Bertrand N. Sabatier J. Briat O. Vinassa J.-M. Embedded fractional nonlinear supercapacitor model and its parametric estimation method IEEE Transactions on Industrial Electronics 2010 57 12 3991 4000 2-s2.0-78449235304 10.1109/TIE.2010.2076307 Martín R. Quintana J. J. Ramos A. de la Nuez I. Modeling of electrochemical double layer capacitors by means of fractional impedance Journal of Computational and Nonlinear Dynamics 2008 3 2 21303 21400 2-s2.0-47249084353 10.1115/1.2833909 021303 Musolino V. Piegari L. Tironi E. New full-frequency-range supercapacitor model with easy identification procedure IEEE Transactions on Industrial Electronics 2013 60 1 112 120 Wang Y. Hartley T. Lorenzo C. Adams J. Modeling Ultracapacitors as Fractional-Order Systems 2010 New Trends in Nanotechnology and Fractional Calculus Applications Spyker R. L. Classical equivalent circuit parameters for a double-layer capacitor IEEE Transactions on Aerospace and Electronic Systems 2000 36 3 829 836 2-s2.0-0034230134 10.1109/7.869502 Dzieliński A. Sarwas G. Sierociuk D. Time domain validation of ultracapacitor fractional order model Proceedings of the 49th IEEE Conference on Decision and Control (CDC '10) December 2010 3730 3735 2-s2.0-79953138374 10.1109/CDC.2010.5717093 Freeborn T. J. Maundy B. Elwakil A. S. Measurement of supercapacitor fractional-order model parameters from voltage-excited step response IEEE Journal on Emerging and Selected Topics in Circuits and Systems 2013 3 3 416 424 de Levie R. On porous electrodes in electrolyte solutions. I. Capacitance effects Electrochimica Acta 1963 8 10 751 780 2-s2.0-25744461413 Kaus M. Kowal J. Sauer D. U. Modelling the effects of charge redistribution during self-discharge of supercapacitors Electrochimica Acta 2010 55 25 7516 7523 2-s2.0-77955766951 10.1016/j.electacta.2010.01.002