Journal of Applied Mathematics

Volume 2014 (2014), Article ID 726108, 7 pages

http://dx.doi.org/10.1155/2014/726108

Research Article

A New Model of Ultracapacitors Based on Fractal Fundamentals

Xiaodong Zhang,¹ Yi Sun,² Jia Song,¹ Aijun Wu,³ and Xiaoyan Cui⁴

¹School of Electrical Engineering, Beijing Jiaotong University, Beijing 100044, China

²China Electrotechnical Society, Beijing 100823, China

³Beijing Railway Bureau, Beijing 100860, China

⁴Automation School, Beijing University of Posts and Telecommunications, Beijing 100876, China

Received 24 January 2014; Accepted 25 June 2014; Published 23 July 2014

Academic Editor: Hongjie Jia

Copyright © 2014 Xiaodong Zhang 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.

Abstract.

An intuitive model is proposed in this paper to describe the electrical behavior of certain ultracapacitors. The model is based on a simple expression that can be fully characterized by five real numbers. In this paper, the measured impedances of three ultracapacitors as a function of frequency are compared to model results. There is good agreement between the model and measurements. Results presented in a previous study are also reviewed and the paper demonstrates that those results are also consistent with the newly described model.

1. Introduction

Smart grids and energy systems always need energy storage. Ultracapacitors represent an alternative to batteries for storing electrical energy and can help to compensate for the limited power density of batteries [1]. They resemble rechargeable batteries in terms of their ability to transport and store charge, but they employ a very different charge storage mechanism. Ultracapacitors store electric energy by accumulating and separating opposite charges physically, as opposed to batteries, which store energy chemically [2]. Opposing charges are separated by an electrode/electrolyte interface, which is referred to as an electrochemical double-layer. Compared to batteries, ultracapacitors have a much longer charge-discharge cycle life [2, 3]. The power density of ultracapacitors is considerably higher than that of batteries, and the energy density is higher than that of electrolytic capacitors for power applications. Ultracapacitors can store a high level of energy in a small volume and release this energy in a powerful burst [4]; so they are useful in power electronic systems and applications (e.g., power systems, automotive, telecommunication, and military) that need to provide or absorb sudden current surges [5]. The power output of ultracapacitors is limited by their internal impedance. The impedance needs to be identified and characterized in order to develop models for different applications. The development of these models requires measurements of their dynamic electrical behavior [2–4].

Often, ultracapacitors are modeled with simple circuits as shown in Figure 1. Models like this are sufficient for well-defined and stable electrical signals; however, they do not accurately describe the electrical behavior of these devices in dynamic and high-power situations.

Figure 1: circuit model for ultracapacitors.

Practical ultracapacitor models are more complicated. Conway [6] described the charge storage mechanism as a Faradic pseudocapacitance involving a redox reaction of microporous transition metal hydrous oxides. Some [3, 7] have modeled the behavior of ultracapacitors using transmission line equivalent circuits. The porous electrode is described by a line of and elements representing the elemental double layer capacitance and the respective electrolyte resistance at a particular pore depth. A more complete description of the porous electrode behavior was given by De Levie [8]. Gualous et al. [4] took into account the physics of the ultracapacitor and proposed an equivalent circuit that described the ultracapacitor electrical behavior with two branches as shown in Figure 2. This model considers the nonlinear relationship between the capacitance of activated carbon particles and their surface area that varies with the type of activated carbon used and the way it is treated. But the model has six parameters to be identified and does not accurately describe ultracapacitor behavior at low frequencies.

Figure 2: Equivalent circuit developed by Gualous et al.

Zubieta and Bonert [9] provided a model for the terminal behavior of an ultracapacitor based on physical reasoning. This model (see Figure 3) has three distinct time constants. , with the voltage-dependent capacitor () which reflects the voltage dependence of the capacitance dominates the behavior of the ultracapacitor in the time over a period of seconds in response to a charge action. The other branches, with and with , separately determine the terminal behavior in the range of minutes and the behavior for times longer than 10 mins. The resistor is a leakage resistor representing the self discharge property of the ultracapacitor. Others (e.g., [10]) have also developed equivalent circuit models based on variable time constants to fit the measured AC impedance data of ultracapacitors. In [10], the dependence of the resistance on frequency was divided into four distinct frequency zones and this model included a voltage-dependent capacitor as well.

Figure 3: Equivalent circuit developed by Zubieta and Bonert.

However, the pores of activated carbon in an ultracapacitor have a complex branch pore structure for which any impedance analysis model should account. Furthermore, the parameters in the above models, especially the voltage dependent capacitance, can be difficult to quantify.

Buller et al. [11, 12] presented a model shown in Figure 4, which can be used to describe the behavior of ultracapacitors over a wide range of frequencies. is the complex pore impedance related to the porosity of the ultracapacitor. The mathematical expression for is

Unfortunately, this function is not well suited for circuit simulation software because of the coth term [13]. Moreover, when the number is high, the series denominators of the formula in Figure 4

become very large and may not be accurately determined by measurements. Buller et al. [12] used 10 cells and 20 parameters to model a 1400 F ultracapacitor.

Figure 4: Equivalent circuit developed by Buller et al.

Qu and Shi [7] proposed an -ladder network model for ultracapacitors based on the pore structure of activated carbon which is shown in Figure 5. This model is particularly intuitive, because it illustrates how more capacitance becomes available as the time constant of the charging cycle is increased. and    can be treated as the resistance and capacitance of the pores with certain pore size. gives the unit of time and indicates how fast the pore of certain size is.

Figure 5: Equivalent circuit developed by Qu and Shi.

Itagaki et al. [14] proposed a model for ultracapacitors based on a fractal pore structure with three sizes of cylindrical pores. This model consisted of resistors and complex impedances connected in a tree-like structure.

This paper proposes a new electrical model for ultracapacitors that is also based on a fractal interpretation of their structure, but does not make any assumption about the size or shape of the pores. The new model is relatively simple and can be fully described by five parameters.

2. Model Description

Ultracapacitors are comprised of two highly porous activated carbon electrodes, which take surface area and charge separation distance to an extreme. The surface areas can be greater than 21,500 square feet per gram and the separation between the charged surfaces is reduced to distances on the order of nanometers [15]. The electrodes are immersed in a suitable electrolyte to facilitate the charge transfer and storage mechanism. Charges accumulate in the pores resulting in capacitance.

Consider the sample pore structure illustrated in Figure 6, a central conducting structure called a “post” is lined with many smaller structures referred to as “teeth.” Half of the posts in this structure are connected to one electrode of an ultracapacitor and the other half are connected to the other electrode. The ultracapacitor has posts, and every post has teeth. In this model, each tooth is a smaller version of the post that it is attached to.

Figure 6: Pore model of the ultracapacitor electrodes.

If we start by analyzing the capacitance of the structure without any teeth, the total capacitance is times to the post-to-post capacitance, because all of these capacitances are in parallel. The equivalent series resistance of the ultracapacitor is the resistance of one pair of posts divided by , since the post resistances are also in parallel, therefore, an equivalent circuit for an ultracapacitor consisting of interleaved posts of uniform size would be a simple circuit, where the value of is and the value of is . In this case, is the resistance associated with charge moving from one electrode into a single post and then to the second electrode. is the contribution to the capacitance of a single pair of posts. The time constant associated with charge populating the posts (i.e., the time constant) is

is proportional to the length of the post and inversely proportional to the cross-sectional area:

where is the electrical conductivity of the post, is the length of the post, and is the cross-sectional area of the post.

The capacitance, , can be expressed as the capacitance between any two posts times a constant that is determined by the number of posts in proximity to a given post. A uniform post distribution is illustrated in Figure 7. is proportional to the length and inversely proportional to the natural log of the ratio of the post separation , to the post radius , as indicated below:

where is the permittivity of the dielectric.

Figure 7: Capacitance scheme of porosity.

A circuit model for the post-only ultracapacitor is provided in Figure 8.

Figure 8: The post-only model.

Now, consider the structure of the intermeshed teeth between posts. In fractal geometry, each tooth would be a scaled down version of the post to which it was attached. If we assume that each post has teeth and the size of each tooth is times the size of the post, then the resistance of a single tooth is

The capacitance, , can be expressed as

For teeth per post in parallel, the overall capacitance is multiplied by and the resistance is divided by , yielding

From (8) and (9), it is clear that the time constant associated with charge moving from the posts out onto the teeth is the same as the time constant associated with charge moving out onto the posts and can be calculated with (3). The equivalent circuit for the capacitor with both posts and teeth is shown in Figure 9. Both resistors in this model have the same value. Both capacitors also have the same value.

Figure 9: Equivalent circuit for a capacitor with one set of posts and teeth.

Employing (6) through (9), it is relatively straightforward to show that adding smaller teeth to each tooth in the structure shown in Figure 6 would result in the same amount of additional capacitance provided through the same amount of additional resistance. The fractal geometry with many layers of repeating structures would yield the equivalent circuit in Figure 10, where all of the resistors would have the same value and so do all of the capacitors.

Figure 10: Electrical model of a fractal electrode structure.

Although the example in Figure 6 employs cylindrical pore structures, the resistance and capacitance will scale similarly with nearly any branch-like pore structure. Generally, the resistance of a branch will be proportional to its length and inversely proportional to its cross-sectional area. The capacitance between closely spaced branches will be proportional to its length and independent of cross-sectional area if the spacing between branches is also scaled. Therefore, the model in Figure 10 does not assume a particular pore structure as long as the successively smaller branches of the pores form fractal geometry.

As the different number of time between charging and discharging, the deterioration of ultracapacitors is also not the same. Few electrolytes will be decomposed to form the free insoluble product particles, which increase the resistance. Therefore, a resistor is paralleled in the model in Figure 10. The greater the loss is, the higher the resistor value will be.

Considering its ohmic resistance characteristics and ohmic capacitive characteristics, a resistor and a capacitor are in series; thus, the model in Figure 11 is proposed, which has 5 independent component values.

Figure 11: Ultracapacitor model described by only 5 real numbers.

Although the model in Figure 11 has only 5 independent component values, it has an infinite number of elements. In order to develop a closed-form expression for the input impedance, it is convenient to use a Fourier transform technique. For the impedance ladder network in Figure 12, is the equivalent input complex impedance and is the equivalent input complex impedance of the -order ladder network (see Figure 13). So

The recurrence formula for is

The solution of (11) can be obtained using a method presented in [16, 17]. In this case, the coefficient matrix of which is expressed by is shown as follows:

in (13) is the coefficient matrix of which is expressed by . If is expressed by its characteristic values and characteristic vectors, can be calculated as

The characteristic values of and its characteristic vectors are

Combining (16) into (15) results in

where

If , (18) can be simplified using (19) and as

It is assumed that

Therefore,

If , then the ,

Or else, if , then the ,

For the porous electrode model, , .

Figure 12: Complex impedance ladder network.

Figure 13: The equivalent circuit of ladder network in Figure 12.

In this case, , and (24) applies. Therefore, the formula for the electrode impedance in Figure 10 is

and the input impedance of the ultracapacitor model in Figure 11 is

3. Validation of Model

Electrochemical impedance spectroscopy (EIS) is one of the most frequently used analytical tools for the characterization of ultracapacitors. This method was used to determine input resistance and reactance of the ultracapacitor as a function of frequency at a given excitation voltage.

Two ultracapacitors with nominal values of 3.3 F and 2 F were measured at room temperature. Potentiostat cycling tests were performed in the capacitor designing lab of St. Jude Medical’s Cardiac Rhythm Manufacture, the Liberty, S.C., with a Gamry potentiostat (Reference 3000). The impedance was measured by applying a sinusoidal 5 mV excitation superimposed on a 1.2 V DC bias to the ultracapacitor and measuring the magnitude and phase of the current. The frequency range was from 100 mHz to 1 kHz.

In this study, a Panasonic 3.3 F/2.3 V (EEC-HW0D335) ultracapacitor which is a new one and a Taiyo Yuden 2 F/2.3 V (PAS1016LR2R3205) which had been charged and discharged for many times were measured. Figures 14(a) and 14(b) show the magnitude of the input impedance plotted as a function of the real part of the input impedance for the two ultracapacitors. At frequencies above 100 Hz, the electrical behavior of the ultracapacitors is more like a simple resistor than a capacitor. At low frequencies , the imaginary part of the impedance dominates.

(a) Nyquist plot of a 3.3 F ultracapacitor (values for the model calculation: , , , , and )

(b) Nyquist plot of a 2 F ultracapacitor (values for the model calculation: , , , , and )

Figure 14: The Comparison between measured numbers and calculated numbers in Nyquist plot.

The values of , , , , and for each capacitor model were obtained by curve-fitting the measured data to the model input impedance in (26). It seems that the value of old one (2 F) is larger than the new one (3.3 F). The experimental data curves of the real part and the imaginary part of were, respectively, fit using the MATLAB/CF Tool software. Figure 14 also shows the impedance obtained from the model in Figure 11. The plots show excellent agreement between the model and the measured values. Figure 15 shows the measured impedance of a new KAMCAP 400 F/2.7 V (HP-2R7-J407UY LL) ultracapacitor with a DC bias voltage of 1.2 V and 10 mV disturbance variable from 100 mHz to 100 Hz at +25°C. This experiment was carried out by Zahner IM6eX. The model data is also shown in Figure 15. In this model  Ohm,  F,  Ohm,  Ohm, and  F. It also shows excellent agreement between the measured and modeled values.

Figure 15: Complex-plane representation of a 400 F ultracapacitor (values for the model calculation: , , , , and ).

4. Conclusion

A new model for describing the electrical behavior of ultracapacitors is introduced based on an intuitive representation of the electrode pores as a fractal structure. The new model has an infinite number of elements but is fully represented by only five real numbers. A closed-form expression for the input impedance was derived making it relatively easy to fit measured results to the model.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work is supported by the Fundamental Research Funds for the Central Universities (no. 2014JBZ017).

References

P. Barrade, “Series connection of supercapacitors: comparative study of solutions for the active equalization of the voltages,” in Proceedings of the 7th International Conference on Modeling and Simulation of Electric Machines, Converters and Systems Electrimacs, Montréal, Canada, August 2002.

J. Van Mierlo, P. Van den Bossche, and G. Maggetto, “Models of energy sources for EV and HEV: Fuel cells, batteries, ultracapacitors, flywheels and engine-generators,” Journal of Power Sources, vol. 128, no. 1, pp. 76–89, 2004.

R. Kötz and M. Carlen, “Principles and applications of electrochemical capacitors,” Electrochimica Acta, vol. 45, no. 15-16, pp. 2483–2498, 2000.

H. Gualous, D. Bouquain, A. Berthon, and J. M. Kauffmann, “Experimental study of supercapacitor serial resistance and capacitance variations with temperature,” Journal of Power Sources, vol. 123, no. 1, pp. 86–93, 2003.

P. L. Taberna, P. Simon, and J. F. Fauvarque, “Electrochemical characteristics and impedance spectroscopy studies of carbon-carbon supercapacitors,” Journal of the Electrochemical Society, vol. 150, no. 3, pp. A292–A300, 2003.

B. E. Conway, Electrochemical Supercapacitors, Kluwer Academic Publishers/Plenum Press, New York, NY, USA, 1999.

D. Qu and H. Shi, “Studies of activated carbons used in double-layer capacitors,” Journal of Power Sources, vol. 74, no. 1, pp. 99–107, 1998.

R. De Levie, in Advances in Electrochemistry and Electrochemical Engineering, P. Delahay and C. T. Tobias, Eds., vol. 6, pp. 329–397, John Wiley & Sons, New York, NY, USA, 1967.

L. Zubieta and R. Bonert, “Characterization of Double-Layer Capacitors (DLCs) for power electronics applications,” in Proceedings of the IEEE Industry Applications Conference, pp. 1149–1154, October 1998.

F. Rafik, H. Gualous, R. Gallay, A. Crausaz, and A. Berthon, “Frequency, thermal and voltage supercapacitor characterization and modeling,” Journal of Power Sources, vol. 165, no. 2, pp. 928–934, 2007.

E. Karden, S. Buller, and R. W. De Doncker, “A frequency-domain approach to dynamical modeling of electrochemical power sources,” Electrochimica Acta, vol. 47, no. 13-14, pp. 2347–2356, 2002.

S. Buller, E. Karden, D. Kok, and R. W. de Doncker, “Modeling the dynamic behavior of supercapacitors using impedance spectroscopy,” IEEE Transactions on Industry Applications, vol. 38, no. 6, pp. 1622–1626, 2002.

D. Riu, N. Retière, and D. Linzen, “Half-order modelling of supercapacitors,” in Proceedings of the 39th IAS Annual Meeting on Record of the Industry Applications Conference, vol. 4, pp. 2550–2554, IEEE, October 2004.

M. Itagaki, S. Suzuki, I. Shitanda, K. Watanabe, and H. Nakazawa, “Impedance analysis on electric double layer capacitor with transmission line model,” Journal of Power Sources, vol. 164, no. 1, pp. 415–424, 2007.

M. Spandana and A. Jagruthi, “Energy Storage Ultracapacitors,” Yuva Engineers, 2nd BTech, EEE, 2010

X. Yuan, “Extraction of the general term formula of fractional recursive sequence of number— by applying matrix method,” Journal of China West Normal University (Natural Sciences), vol. 26, no. 1, pp. 102–104, 2005.

Y. Jin and Q. Wang, “The universal solution of ladder network equivalent impedance,” Journal of Higher Correspondence Education, vol. 23, no. 4, 2010.