Next Article in Journal
Performance Assessment of a nZEB Carbon Neutral Living/Office Space and Its Integration into a District Energy-Hub
Next Article in Special Issue
Recycling of Lithium Batteries—A Review
Previous Article in Journal
External Environmental Analysis for Sustainable Bike-Sharing System Development
Previous Article in Special Issue
Fabrication of Silicon Nanowire Metal-Oxide-Semiconductor Capacitors with Al2O3/TiO2/Al2O3 Stacked Dielectric Films for the Application to Energy Storage Devices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Constant Phase Element in the Time Domain: The Problem of Initialization

by
Juan Antonio López-Villanueva
* and
Salvador Rodríguez Bolívar
Department of Electronics and Computer Technology, CITIC, University of Granada, 18071 Granada, Spain
*
Author to whom correspondence should be addressed.
Energies 2022, 15(3), 792; https://doi.org/10.3390/en15030792
Submission received: 15 December 2021 / Revised: 13 January 2022 / Accepted: 17 January 2022 / Published: 21 January 2022
(This article belongs to the Special Issue Advances in Devices for Energy Generation and Storage)

Abstract

:
The constant phase element (CPE) is found in most battery and supercapacitor equivalent circuit models proposed to interpret data in the frequency domain. When these models are used in the time domain, the initial conditions in the fractional differential equations must be correctly imposed. The initial state problem remains controversial and has been analyzed by various authors in the last two decades. This article attempts to clarify this problem by proposing a procedure to prepare the initial state and defining a decay function that reveals the effect of the initial state in several illustrative examples. This decay function depends on the previous history, which is reflected in the time needed to prepare the initial state and on the current profile assumed for this purpose. This effect of the initial state is difficult to separate and can lead to the misinterpretation of the CPE parameter values.

1. Introduction

Fractional equivalent circuit models are widely used to represent the behavior of batteries and supercapacitors, both in the frequency and time domains [1,2,3,4,5], as they offer better accuracy with a lesser number of elements. These equivalent circuits include constant phase elements (CPEs), which are unequivocally observed in the experimental results obtained with electrochemical impedance spectroscopy (EIS) [6]. The phase of the impedance of a CPE is constant since it depends on the fractional exponent α , but does not depend on the frequency, as shown in (1) [1,7]:
Z C P E ( j ω ) = 1 Q · ( j ω ) α = 1 Q ω α e j α π 2
In addition, the phase is negative, so that this element is sometimes referred to as a fractional capacitor. The impedance of the CPE in the Laplace domain is:
Z C P E ( s ) = 1 Q s α .
There is interest in applying these models in the time domain [2,3,8,9,10] where the impedance (2) corresponds to a fractional derivative, as shown in (3), which is discussed in the next section.
i ( t ) = Q d α v ( t ) d t α
To solve differential equations that include fractional derivatives, the initial conditions must be handled appropriately. However, the initial state problem remains controversial and has been discussed by various authors in the last two decades. Lorenzo and Hartley first recognized that a time-varying initialization is actually required for the fractional-order derivative [11]. Some other authors also indicated that initial histories rather than initial values at a point should be considered [12,13,14]. The term “aberration” has been introduced in relation to this phenomenon [15] and has been used by other authors [16]. One procedure to deal with this problem is to define an initialized operator, related to the history, to be imposed in the fractional derivative at the initial point [17,18]. An alternative proposal is the diffusive representation introduced by Sabatier [19] and the continuous-frequency-distributed model proposed by Trigeassou and Maamri [20], which were merged in the construction of a partial differential equation, where the initial state is a vector of infinite dimensions [20,21,22]. The equivalence and compatibility of the above two initialization theories have been demonstrated [18,23].
We focus the analysis of this work on a fractional element with impedance given by (2), with 0 < α < 1 , which is the particular case of greatest interest in battery and supercapacitor models. One example of it corresponds to α = 0.5 , which is often called Warburg impedance [24,25].
In this paper, we bring this problem to the attention of the battery and supercapacitor modeling community. Our main goal is to provide a detailed study of how the previous history, i.e., how the initial voltage - v ( 0 ) - is reached, affects the subsequent behavior of the system. To do this, we hypothesized that the initial voltage of a CPE with 0 < α < 1 may be written as v ( 0 ) g ( t ) , where g ( t ) is a decay function between zero and one, with intermediate values and behavior between that of the resistor and that of the non-fractional capacitor. This decay function, closely related to the general initial function defined by other authors [17,18,26], is derived and analyzed for the CPE, ZARC (composed of the CPE in parallel with a resistor), and RC networks. In the paper, we also highlight the error tow which some definitions of the fractional derivative can lead.
A brief review of the basic definitions of fractional derivatives is provided in Section 2. Next, the decay function is calculated and discussed in Section 3 for a CPE, analyzing how the previous history of the CPE influences the time dependence of the function. The case of a ZARC element is discussed in Section 4, and the approximation of a CPE by a series RC network is considered in Section 5. Finally, a discussion of the results is presented in Section 6, and some conclusions and possible future work are provided in Section 7.

2. Fractional Derivatives

The most commonly used definitions of fractional derivatives start from the Riemann–Liouville fractional integral of order α , given by [27]:
0 I t α f ( t ) = 0 D t α f ( t ) = 1 Γ ( α ) 0 t f ( τ ) ( t τ ) 1 α d τ
where Γ ( α ) is the Gamma function [27].
Using Definition (4), the Riemann–Liouville fractional derivative is defined as [27]:
0 R L D t α f ( t ) = 0 D t m 0 D t α m f ( t ) ,
with m 1 α < m , m a positive integer, so that:
0 R L D t α f ( t ) = 1 Γ ( m α ) d m d t m 0 t f ( τ ) ( t τ ) 1 m + α d τ .
Applying the Laplace transformation to (5), the following equation is obtained [27]:
L [ 0 R L D t α f ( t ) ] = s α L [ f ( t ) ] k = 0 m 1 s k [ 0 R L D t α k 1 f ( t ) ] t = 0
This expression depends on fractional derivatives at the initial time. This is a drawback of the Riemann–Liouville definition and motivates the use of an alternative definition known as the Caputo fractional derivative:
0 C D t α f ( t ) = 0 D t α m 0 D t m f ( t ) ,
or:
0 C D t α f ( t ) = 1 Γ ( m α ) 0 t 1 ( t τ ) 1 m + α d m f ( τ ) d τ m d τ , m 1 α < m ,
whose Laplace transform is:
L [ 0 C D t α f ( t ) ] = s α L [ f ( t ) ] k = 0 m 1 s α k 1 f ( k ) ( 0 )
Equation (10) only uses integer derivatives at the initial time, so it is easier to apply in practice. However, the inability of the Caputo definition to provide a satisfactory solution to the initial condition problem has been pointed out [13,15], and we comment on it below.

3. Initialization of the CPE

The CPE, or fractional capacitor, with 0 < α < 1 , can be considered as an intermediate case between a resistor ( α = 0 ) and an ideal non-fractional capacitor ( α = 1 ). It is reasonable to think that it also behaves as an intermediate case in relation to the initial conditions. Allagui et al. [28] recently claimed that the weight of the voltage memory trace that results from the contribution of past voltage activity depends on the fractional exponent α , and therefore, the measured response of the device at any time is increasingly correlated with its past, while ideal capacitors do not exhibit such behavior and are memoryless devices. We interpret this same result in another way. If this two-terminal element supports an initial voltage drop, v ( 0 ) , at time t = 0 , and we apply a current i ( t ) for time t > 0 , the behavior in the two extreme cases is quite different. In the case of a resistor with resistance R, the voltage at time t > 0 is:
v ( t ) = R i ( t ) ,
and the initial voltage is immediately forgotten. In the opposite case of a non-fractional capacitor of capacitance C, the voltage is given by:
v ( t ) = v ( 0 ) + 1 C 0 t i ( τ ) d τ ,
and the initial voltage is always remembered, although its later effects do not depend on the previous history.
Let us apply (10) to (3), in the particular case of 0 < α < 1 , to obtain:
I ( s ) = Q { s α V ( s ) s α 1 v ( 0 ) } ,
or:
V ( s ) = I ( s ) Q s α + v ( 0 ) s ,
which in the time domain, for t > 0 , corresponds to:
v ( t ) = v ( 0 ) + 1 Q Γ ( α ) 0 t i ( τ ) ( t τ ) 1 α d τ .
Time-domain models for predicting the voltage–current characteristics of lithiumion batteries under arbitrary charging and discharging current profiles have been reported, showing a higher accuracy than the traditional models with multiple parallel RC units [29,30]. The initial voltage has been considered in an equation similar to (14). With this way of initialization, v ( 0 ) acts as an integration constant and is always remembered, as we would expect in the case of an ideal non-fractional capacitor, that it is not what we should obtain in the limit α 0 . Instead, we would expect to obtain:
v ( t ) = v ( 0 ) g ( t ) + 1 Q Γ ( α ) 0 t i ( τ ) ( t τ ) 1 α d τ
so that the Caputo definition does not seem to produce the correct results. The meaning of g ( t ) is discussed in the next section.
An alternative to the Riemann–Liouville and Caputo definitions of the fractional derivative is to use the discretized and truncated Grünwald–Letnikov derivative with a short and fixed memory length [5]. However, the voltage values during this memory interval, before the initial time, must also be known, so that the initialization problem is not solved by this approach.

4. The Decay Function

To obtain the value of g ( t ) , which we call the decay function, we propose to prepare an initial state starting from a state with zero initial conditions. The idea of preparing an initial state is not actually new [19]. It was applied in [13] and [14] to compute the internal force of an axially loaded viscoelastic bar, where an initial state was defined over an interval, and the initial condition was called the initial history [13]. The initialization function is a time-varying function and can be viewed as a generalization of the constant of integration required for the order-one integral [31]. Zhao et al. [26] also applied this procedure to the case of the internal force of an axially loaded viscoelastic bar and showed that the influence caused by the pre-initial process varies in time and reflects the memory of the whole pre-initial process, which explains the origin of the long memory property of fractional-order systems [26].
We distinguish the initialization interval, [ 0 , t 0 ] , and the observation interval, ( t 0 , ) , and assume that the initial state has been achieved at t 0 with an initialization current i 0 while the current applied during the observation phase is i 1 , according to the next equation.
i ( t ) = i 0 ( t ) 0 t < t 0 i 1 t t 0 t t 0
We also assumed zero initial conditions prior to the initialization phase. With zero initial conditions, the two definitions of fractional derivative (6) and (9) lead to the same Laplace transform and the ambiguity is avoided.
To obtain the initial voltage v 0 , we can apply a suitable constant current i 0 ( t ) = I 0 . With this current, the voltage can be obtained by:
v ( t ) = I 0 Q Γ ( α ) 0 t d τ ( t τ ) 1 α = I 0 t α Q Γ ( α + 1 ) ,
If I 0 is held for a time interval of length t 0 ,
v ( t 0 ) = I 0 t 0 α Q Γ ( α + 1 ) v 0 .
After t 0 , an arbitrary current i ( t ) flows through the CPE. The voltage is then obtained as:
v ( t ) = I 0 Q Γ ( α ) 0 t 0 d τ ( t τ ) 1 α + 1 Q Γ ( α ) t 0 t i ( τ ) ( t τ ) 1 α d τ .
The first term in (20) is:
I 0 Q Γ ( α ) 0 t 0 d τ ( t τ ) 1 α = I 0 Q Γ ( α + 1 ) [ t α ( t t 0 ) α ] ,
which can be rewritten, by using (19), as:
I 0 Q Γ ( α ) 0 t 0 d τ ( t τ ) 1 α = v 0 t t 0 α t t 0 t 0 α .
By defining the observation time, t , after t 0 ,
t t t 0 ,
the final result has the form expected in (16):
v ( t ) = v 0 g 1 ( t 0 ; t ) + 1 Q Γ ( α ) 0 t i ( τ + t 0 ) ( t τ ) 1 α d τ = v 0 g 1 ( t 0 ; t ) + 1 Q Γ ( α ) 0 t i 1 ( τ ) ( t τ ) 1 α d τ ,
with:
g 1 ( t 0 ; t ) = 1 + t t 0 α t t 0 α .
Function g 1 verifies g 1 ( t 0 ; 0 ) = 1 and l i m t g 1 ( t 0 ; t ) = 0 , so that the initial voltage, v 0 , is eventually forgotten, which justified the name “decay function” we used for it. However, although this procedure gives a method to prepare an initial state, the function g 1 depends on the time t 0 , which highlights the importance of taking into account the previous history since the same initial voltage can be obtained with different combinations of I 0 and t 0 values, and different decay rates of the initial voltage are therefore obtained.
The result also depends on the current profile used to prepare the initial state. To illustrate this, we can use a ramp current instead of a constant current for the same purpose. Let us assume a current i 0 ( t ) = a · t for 0 t < t 0 . The voltage for this time interval is then:
v ( t ) = a Q Γ ( α ) 0 t τ · d τ ( t τ ) 1 α = a t α + 1 Q Γ ( α + 2 )
and:
v ( t 0 ) = a t 0 α + 1 Q Γ ( α + 2 ) v 0 .
If an arbitrary current i ( t ) is applied after t 0 ,
v ( t ) = a Q Γ ( α ) 0 t 0 τ d τ ( t τ ) 1 α + 1 Q Γ ( α ) t 0 t i ( τ ) ( t τ ) 1 α d τ ,
where:
a Q Γ ( α ) 0 t 0 τ d τ ( t τ ) 1 α = a Q Γ ( α + 2 ) [ t α + 1 ( t t 0 ) α + 1 ( α + 1 ) t 0 ( t t 0 ) α ] .
Using (27),
a Q Γ ( α ) 0 t 0 τ d τ ( t τ ) 1 α = v 0 t t 0 α + 1 t t 0 t 0 α + 1 ( α + 1 ) t t 0 t 0 α
With t t t 0 , the final result is obtained as:
v ( t ) = v 0 g 2 ( t 0 ; t ) + 1 Q Γ ( α ) 0 t i ( τ ) ( t τ ) 1 α d τ ,
where:
g 2 ( t 0 ; t ) = 1 + t t 0 α + 1 1 + α + t t 0 t t 0 α
Again, g 2 ( t 0 ; 0 ) = 1 and l i m t g 2 ( t 0 ; t ) = 0 , so that g 2 is another form of the decay function obtained for a different initialization current profile. The decay functions g 1 and g 2 are shown in Figure 1 for α = 0.5 and two different values of t 0 . The decay rates are slow for long times and are different in the four cases, thus proving the influence of the previous history.

5. Initialization of ZARC

The CPE included in battery and supercapacitor models is often connected with a resistor in parallel, as shown in Figure 2. For this connection, the name ZARC was proposed earlier [32] and has been widely used ever since [1,2,10]. Its impedance in the Laplace domain is:
Z = R 1 + R Q s α
Instead of parameter Q, we can use a time constant, which can be used as a time reference, by defining τ R Q 1 α , so that:
Z = R 1 + ( τ s ) α
and:
1 + ( τ s ) α V ( s ) = I ( s ) R .
We apply now the same initialization procedure used above with a constant current, I 0 , to prepare the initial state. An arbitrary current i 1 is subsequently applied for t t 0 :
i ( t ) = I 0 0 t < t 0 i 1 t t 0 t t 0
Equation (36) can be rewritten as:
i ( t ) = I 0 [ u ( t ) u ( t t 0 ) ] + i 1 ( t t 0 ) u ( t t 0 ) ,
where u ( t ) is the unit step function. Its Laplace transform is [29]:
I ( s ) = I 0 1 e s t 0 s + I 1 ( s ) e s t 0 .
so that:
V ( s ) = I 0 R 1 e s t 0 s 1 + ( τ s ) α + I 1 ( s ) R e s t 0 1 + ( τ s ) α .
For t t 0 , or t 0 , where t = t t 0 , we obtain:
v ( t ) = v i n i ( t 0 ; t ) + v 1 ( t )
where:
v 1 ( t ) = L 1 I 1 ( s ) R 1 + ( τ s ) α
is the voltage that we would obtain in response to i 1 with zero initial conditions, and:
v i n i ( t 0 ; t ) = L 1 I 0 R ( 1 e s t 0 ) s 1 + ( τ s ) α .
is the response to the initial voltage in the observation phase. By using [27]:
L 1 1 s 1 + ( τ s ) α = 1 E α t τ α
where E α ( x ) is the one-parameter Mittag–Leffler function [27], defined according to:
E α ( x ) = n = 0 x n Γ ( n α + 1 )
the result, with t = t t 0 , is:
v i n i ( t 0 ; t ) = I 0 R 1 E α t + t 0 τ α 1 E α t τ α
The initial voltage can be written as:
v 0 v i n i ( t 0 ; 0 ) = I 0 R 1 E α t 0 τ α ,
so that the final result is:
v i n i ( t 0 ; t ) = v 0 g 3 ( t 0 ; t ) ,
where a new decay function g 3 is defined in this case according to:
g 3 ( t 0 ; t ) = E α t τ α E α t + t 0 τ α 1 E α t 0 τ α .
In the case of a non-fractional ZARC, which is obtained in the limit α = 1 , in which the CPE is substituted by a non-fractional capacitor, and taking into account that the Mittag–Leffler function reduces to the exponential function in this limit:
lim α 1 E α t τ α = e t τ ,
and the decay function is also an exponential function:
lim α 1 g 3 ( t 0 ; t ) = e t τ .
This is the well-known result in the case of the parallel RC network, in which the effect of the initial voltage decays exponentially.
Therefore, both in the fractional ZARC and in the parallel RC network, the effect of the initial voltage decays and eventually vanishes, but while in the case of the parallel RC network, the exponential decay only depends on the initial voltage, in the case of the fractional ZARC, the decay rate is smoother and depends on the previous history. This behavior is compared in Figure 3, for the case of α = 0.5 . The decay function g 3 decays rapidly at first, but the decay rate is quite slower than the exponential for long times.

6. Approximation of the CPR by an RC Network

To avoid the problems associated with the fractional derivatives, the fractional impedances are often approximated by functions with a finite number of poles and zeros. One possibility is that proposed by [33], and used by other authors [8], in which the fractional-order transfer function is approximated with a high-degree integer-order system. Another option is to approximate the CPE (or the ZARC) by a network composed by a series of RC branches [6,34], as shown in Figure 4. The approximation of the ZARC using finite and infinitely serially connected RC circuits was analyzed in detail in [35]. We consider the approximation of a CPE by a series RC network in this section.
We could think that by using non-fractional RC elements, in which the effect of the initial voltage does not depend on the previous history, the problem of the initial conditions that we discussed above can be avoided. However, this is not true, since given an initial voltage v 0 in the network, how this voltage is distributed among the initial voltages of the different elements actually depends on the previous history.
The voltage of the RC network, with N parallel RC elements connected in series, can be obtained as follows [24],
v ( t ) = v 0 g ( t ) + n = 1 N 1 C n 0 t i ( ξ ) e t ξ τ n d ξ ,
with:
g ( t ) = n = 1 N v n , 0 e t τ n n = 1 N v n , 0 ,
where v n , 0 is the initial voltage in the n-th element, and:
v 0 = n = 1 N v n , 0
We prepare the initial state again by starting from the zero initial condition and applying a constant current, I 0 , for a time t 0 . The initial voltage in the n-th element in the observation phase, at the end of this period, is obtained as:
v n , 0 = I 0 R n 1 e t 0 τ n
and the decay function is also obtained in this case:
g ( t ) = g 4 ( t 0 ; t ) = n = 1 N R n 1 e t 0 τ n e t τ n n = 1 N R n 1 e t 0 τ n ,
where time t is measured after the initialization period. The new decay function also has the expected limits l i m t g 4 ( t 0 ; t ) = 0 and l i m t 0 g 4 ( t 0 ; t ) = 1 .
The behavior of the decay function g 4 can be compared with the decay function g 1 of the CPE that the RC network approximates. To do so, we considered a Warburg-type CPE, with α = 0.5 , and approximated it by an RC network with five elements ( N = 5 ). We used typical parameter values of a Warburg impedance for a Li-ion battery [8,34]: (1) The values of R n and τ n = R n C n were obtained by a least-square fit to the real and imaginary parts of the CPE impedance in the range f m i n = 1 MHz to f m a x 1 Hz. (2) The real part of the CPE impedance at 1 MHz is 20 m Ω . The values of the RC network components are provided in Table 1.
Figure 5 shows a comparison of the real and imaginary parts of the impedances of the CPE (solid line) and the RC network (symbols). Since we chose α = 0.5 , the real and imaginary parts of the CPE match. A good fit was obtained within the frequency range used for the adjustment, but deviations were observed outside this range. The two decay functions g 1 and g 4 , given in (25) and (55), respectively, are compared in Figure 6 in logarithmic scales for two different values of t 0 . A good agreement was observed for times less than t m a x 1 2 π f m i n , but a deviation was produced for longer times, in which the initial voltage for the CPE decayed at a lower rate.

7. Discussion

Although the problem of the initialization of fractional differential equations has been studied in the last two decades, it has not received sufficient consideration in some practical applications. The case considered here is the use of constant-phase elements in battery and supercapacitor modeling. When equivalent circuits containing fractional elements are used to predict the evolution of voltage and current in the time domain, the correct treatment of the initial conditions poses considerable difficulties, since many times, the previous history is not well known and its influence has been clearly demonstrated. The response to an applied current has been shown to be superimposed on the initial voltage multiplied by a decay function that depends on the previous current profile and the time for which it has been applied. There are several possibilities. The simplest one is to be sure that we start from a state with zero initial conditions, but this requires long relaxation times and is not always applicable. Another option is to wait for a long period after which the effect of the initial conditions has decayed and has a negligible effect on the element voltage, but this is even more difficult with fractional elements since the decay rates can be quite slow. A third possibility is to try to adjust the decay function from the data actually measured by separating it from the response to the post-initial current. Du et al. mentioned the difficulties of obtaining the initialization function in practical situations [15], although an algorithm to accomplish this task has been recently proposed [36]. In fact, this is not an easy task, as we show in the example below.
Let us prepare an initial voltage of 50 mV in a CPE of the Warburg type according to (1) with α = 0.5 and Q = 446 Ω 1 s 1 2 (such that its real part is 20 m Ω for a frequency of 1 MHz). We then apply a constant current I 1 = 1 A. Measuring the time t after the initialization period and using (24) we find:
v ( t ) = v 0 g 1 ( t 0 ; t ) + I 1 t α Q Γ ( α + 1 ) ,
where function g 1 is given in (25). The result, for three different values of t 0 , namely 100 s, 500 s, and 3600 s, respectively, is shown in Figure 7. The curve corresponding to t 0 = 100 s initially decreases and then shows a turnaround behavior. This is a consequence of a high current in the initialization phase, since to obtain the initial value of 50 mV in t 0 = 100 s, a high initialization current, I 0 = 1.98 A, is needed. If the CPE is assumed to reflect, for example, a diffusion mechanism in a Li-ion battery [37], such a high current would have produced a highly nonuniform Li concentration in its particles [38], and the observed rebound is justified. However, in the other two cases, it is possible to misinterpret the data by assuming the wrong Equation (15) that corresponds to (56) with g 1 = 1 and different values of Q and α . A good fit, shown with symbols in Figure 7, was obtained with α = 0.77 and Q = 4730 Ω 1 s α for the t 0 = 500 s case and α = 0.53 and Q = 718 Ω 1 s α for the t 0 = 3600 s case, different from those used in (56). In the slower initialization case ( I 0 = 0.33 A during 3600 s), the error in Q and α is lower, but still important. The good fit obtained in some cases, such as those considered here with t 0 equal to 500 s and 3600 s, justifies the good result obtained by the authors, who used a non-decaying initial voltage [29,30]. However, the CPE parameter values obtained to represent the experimental results would have been different if a suitable decay function had been included. In other cases, such as the one shown here for t 0 = 100 s, it would have been much more difficult to obtain a good fit to the experimental results.

8. Conclusions and Future Research

To conclude, careful attention should be paid to the initialization of fractional elements, such as the CPE, in battery and supercapacitor models. The effect of the initial value decays with time and eventually vanishes, and it can be written as an initial function obtained as the product of the initial value and a decay function, but the initial function may be very difficult to separate from the rest of the response and may be misinterpreted by different values of the CPE parameters.
In many practical cases, the initial voltage can be measured or deduced if the open circuit voltage of the battery cell is known, but the previous history is usually not known, so it is not possible to obtain an accurate initialization function. However, in cases where the experiment cannot be designed starting from zero initial conditions, the decay functions proposed in this work can be used with an equivalent or effective t 0 parameter. This might not be totally accurate, but it would lead to better results than simply assuming an initial voltage that does not decay. We believe that this subject warrants further research. In those experimental works that seek to adjust the transient experimental results in the time domain with models that include the CPE, in addition to the CPE parameters, a new parameter defining the decay function can be included, and its effects on the accuracy or theoretical validity of the obtained results can be analyzed.

Author Contributions

Conceptualization, review, software, writing and editing, J.A.L.-V. and S.R.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the TIC105 research group of the University of Granada.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Andre, D.; Meiler, M.; Steiner, K.; Walz, H.; Soczka-Guth, T.; Sauer, D.U. A retrospective on lithium-ion batteries. J. Power Sources 2011, 196, 5349–5356. [Google Scholar] [CrossRef]
  2. Wang, B.; Li, S.E.; Peng, H.; Liu, Z. Fractional-order modeling and parameter identification for lithium-ion batteries. J. Power Sources 2015, 293, 151–161. [Google Scholar] [CrossRef]
  3. Bertrand, N.; Sabatier, J.; Briat, O.; Vinassa, J.M. Embedded Fractional Nonlinear Supercapacitor Model and Its Parametric Estimation Method. IEEE Trans. Ind. Electron. 2010, 57, 3991–4000. [Google Scholar] [CrossRef]
  4. Allagui, A.; Freeborn, T.J.; Elwakil, A.S.; Maundy, B.J. Reevaluation of Performance of Electric Double-layer Capacitors from Constant-current Charge/Discharge and Cyclic Voltammetry. Sci. Rep. 2016, 6, 38568. [Google Scholar] [CrossRef] [PubMed]
  5. Hu, X.; Yuan, H.; Zou, C.; Li, Z.; Zhang, L. Co-Estimation of State of Charge and State of Health for Lithium-Ion Batteries Based on Fractional-Order Calculus. IEEE Trans. Veh. Technol. 2018, 67, 10319–10329. [Google Scholar] [CrossRef]
  6. Gagneur, L.; Driemeyer-Franco, A.; Forgez, C.; Friedrich, G. Modeling of the diffusion phenomenon in a lithium-ion cell using frequency or time domain identification. Microelectron. Reliab. 2013, 53, 784–796. [Google Scholar] [CrossRef]
  7. Fouda, M.E.; Allagui, A.; Elwakil, A.S.; Das, S.; Psychalinos, C.; Radwan, A.G. Nonlinear charge-voltage relationship in constant phase element. AEU Int. J. Electron. Commun. 2020, 117, 153104. [Google Scholar] [CrossRef]
  8. Alavi, S.M.M.; Birkl, C.R.; Howey, D.A. Time-domain fitting of battery electrochemical impedance models. J. Power Sources 2015, 288, 345–352. [Google Scholar] [CrossRef]
  9. Allagui, A.; Freeborn, T.J.; Elwakil, A.S.; Fouda, M.E.; Maundy, B.J.; Radwan, A.G.; Said, Z.; Abdelkareem, M.A. Review of fractional-order electrical characterization of supercapacitors. J. Power Sources 2018, 400, 457–467. [Google Scholar] [CrossRef]
  10. Nasser-Eddine, A.; Huard, B.; Gabano, J.D.; Poinot, T. A two steps method for electrochemical impedance modeling using fractional order system in time and frequency domains. Control Eng. Pract. 2019, 86, 96–104. [Google Scholar] [CrossRef]
  11. Hartley, T.T.; Lorenzo, C. Dynamics and Control of Initialized Fractional-Order Systems. Nonlinear Dyn. 2002, 29, 201–233. [Google Scholar] [CrossRef]
  12. Fukunaga, M.; Shimizu, N. Role of prehistories in the initial value problems of fractional viscoelastic equations. Nonlinear Dyn. 2004, 38, 207–220. [Google Scholar] [CrossRef]
  13. Du, M.L.; Wang, Z.H. Initialized fractional differential equations with Riemann–Liouville fractional-order derivative. Eur. Phys. J. Spec. Top. 2011, 193, 49–60. [Google Scholar] [CrossRef]
  14. Du, B.; Wei, Y.; Liang, S.; Wang, Y. Estimation of exact initial states of fractional order systems. Nonlinear Dyn. 2016, 86, 2061–2070. [Google Scholar] [CrossRef]
  15. Du, M.; Wang, Z. Correcting the initialization of models with fractional derivatives via history-dependent conditions. Acta Mech. Sin. 2016, 32, 320–325. [Google Scholar] [CrossRef]
  16. Zhao, Y.; Wei, Y.; Chen, Y.; Wang, Y. A new look at the fractional initial value problem: The aberration phenomenon. J. Comput. Nonlinear Dyn. 2018, 13, 121004. [Google Scholar] [CrossRef]
  17. Lorenzo, C.F.; Hartley, T. Initialization of fractional-order operators and fractional differential equations. J. Comput. Nonlinear Dyn. 2008, 3, 021101. [Google Scholar] [CrossRef]
  18. Hartley, T.T.; Lorenzo, C.F.; Trigeassou, J.C.; Maamri, N. Equivalence of history-function based and infinite-dimensional-state initializations for fractional-order operators. J. Comput. Nonlinear Dyn. 2013, 8, 041014. [Google Scholar] [CrossRef]
  19. Sabatier, J.; Merveillaut, M.; Malti, R.; Oustaloup, A. How to impose physically coherent initial conditions to a fractional system? Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 1318–1326. [Google Scholar] [CrossRef]
  20. Trigeassou, J.C.; Maamri, N. Initial conditions and initialization of linear fractional differential equations. Signal Process. 2011, 91, 427–436. [Google Scholar] [CrossRef]
  21. Trigeassou, J.C.; Maamri, N.; Sabatier, J.; Oustaloup, A. State variables and transients of fractional order differential systems. Comput. Math. Appl. 2012, 64, 3117–3140. [Google Scholar] [CrossRef] [Green Version]
  22. Trigeassou, J.; Maamri, N.; Oustaloup, A. Lyapunov stability of noncommensurate fractional order systems: An energy balance approach. J. Comput. Nonlinear Dyn. 2016, 11, 041007. [Google Scholar] [CrossRef]
  23. Yuan, J.; Zhang, Y.; Liu, J.; Shi, B. Equivalence of initialized fractional integrals and the diffusive model. J. Comput. Nonlinear Dyn. 2018, 13, 034501. [Google Scholar] [CrossRef]
  24. Plett, G.L. Battery Management Systems, Volume I: Battery Modeling; Artech House: Noewood, MA, USA, 2015. [Google Scholar]
  25. Agudelo, B.O.; Zamboni, W.; Monmasson, E. A Comparison of Time-Domain Implementation Methods for Fractional-Order Battery Impedance Models. Energies 2021, 14, 4415. [Google Scholar] [CrossRef]
  26. Zhao, Y.; Wei, Y.; Shuai, J.; Wang, Y. Fitting of the initialization function of fractional order systems. Nonlinear Dyn. 2018, 93, 1589–1598. [Google Scholar] [CrossRef]
  27. Petras, I. Fractional-Order Nonlinear Systems. Modeling, Analysis and Simulation; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  28. Allagui, A.; Zhang, D.; Khakpour, I.; Elwakil, A.S.; Wang, C. Quantification of memory in fractional-order capacitors. J. Phys. D Appl. Phys. 2020, 53, 02LT03. [Google Scholar] [CrossRef]
  29. Cheng, C.S.; Chung, H.S.H.; Lau, R.W.H. Time-domain modeling of constant phase element for simulation of lithium batteries under arbitrary charging and discharging current profiles. In Proceedings of the 2017 IEEE Applied Power Electronics Conference and Exposition (APEC), Tampa, FL, USA, 26–30 March 2017; pp. 985–992. [Google Scholar] [CrossRef]
  30. Cheng, C.S.; Chung, H.S.H.; Lau, R.W.H.; Hong, K.Y.W. Time-Domain Modeling of Constant Phase Elements for Simulation of Lithium Battery Behavior. IEEE Trans. Power Electron. 2019, 34, 7573–7587. [Google Scholar] [CrossRef]
  31. Yuan, J.; Gao, S.; Xiu, G.; Shi, B. Equivalence of initialized Riemann–Liouville and Caputo derivatives. J. Appl. Anal. Comput. 2020, 10, 2008–2023. [Google Scholar] [CrossRef]
  32. Macdonald, J.R. New aspects of some small-signal ac frequency response functions. Solid State Ion. 1985, 15, 159–161. [Google Scholar] [CrossRef]
  33. Oustaloup, A.; Levron, F.; Mathieu, B.; Nanot, F. Frequency-band complex noninteger differentiator: Characterization and synthesis. IEEE Trans. Circuits Syst. I Fundam. Theory Appl. 2000, 47, 25–30. [Google Scholar] [CrossRef]
  34. Ren, H.; Zhao, Y.; Chen, S.; Yang, L.A. A comparative study of lumped equivalent circuit models of a lithium battery for state of charge prediction. Int. J. Energy Res. 2019, 43, 7306–7315. [Google Scholar] [CrossRef]
  35. Heil, T.; Hossen, A. Continuous approximation of the ZARC element with passive components. Meas. Sci. Technol. 2021, 32, 104011. [Google Scholar] [CrossRef]
  36. Krishnan, G.; Das, S.; Agarwal, V. An Online Identification Algorithm to Determine the Parameters of the Fractional-Order Model of a Supercapacitor. IEEE Trans. Ind. Appl. 2020, 56, 763–770. [Google Scholar] [CrossRef]
  37. Nasser-Eddine, A.; Huard, B.; Gabano, J.D.; Poinot, T.; Martemianov, S.; Thomas, A. Fast time domain identification of electrochemical systems at low frequencies using fractional modeling. J. Electroanal. Chem. 2020, 862, 113957. [Google Scholar] [CrossRef]
  38. Lopez-Villanueva, J.A.; Rodriguez-Bolivar, S.; Parrilla, L.; Finana, C. Simple Single Particle Model for Interpreting Fast Charge Results in Intercalation Batteries. In Proceedings of the XXXV Conference on Design of Circuits and Integrated Systems (DCIS), Segovia, Spain, 18–20 November 2020; pp. 1–5. [Google Scholar] [CrossRef]
Figure 1. Decay functions of the initial voltage obtained for a CPE with α = 0.5 by applying a constant current (lines) or a ramp current (symbols) for a time t 0 .
Figure 1. Decay functions of the initial voltage obtained for a CPE with α = 0.5 by applying a constant current (lines) or a ramp current (symbols) for a time t 0 .
Energies 15 00792 g001
Figure 2. ZARC composed by a resistor and a CPE connected in parallel.
Figure 2. ZARC composed by a resistor and a CPE connected in parallel.
Energies 15 00792 g002
Figure 3. Decay function of the initial voltage obtained for a ZARC with α = 0.5 by applying a constant current (lines), for different values of t 0 , compared to the exponential decay function of a parallel RC element (dots).
Figure 3. Decay function of the initial voltage obtained for a ZARC with α = 0.5 by applying a constant current (lines), for different values of t 0 , compared to the exponential decay function of a parallel RC element (dots).
Energies 15 00792 g003
Figure 4. Approximation of a CPE by the series connection of N RC elements.
Figure 4. Approximation of a CPE by the series connection of N RC elements.
Energies 15 00792 g004
Figure 5. Absolute values of the real part (resistance) and imaginary part (reactance) of a CPE with α = 0.5 of an approximation with a five-element RC network.
Figure 5. Absolute values of the real part (resistance) and imaginary part (reactance) of a CPE with α = 0.5 of an approximation with a five-element RC network.
Energies 15 00792 g005
Figure 6. Decay function of the initial voltage obtained by applying a constant current to a CPE with α = 0.5 and to an approximation with a five-element RC network.
Figure 6. Decay function of the initial voltage obtained by applying a constant current to a CPE with α = 0.5 and to an approximation with a five-element RC network.
Energies 15 00792 g006
Figure 7. Voltage for a CPE with α = 0.5 and Q = 446 Ω 1 Hz 1 2 , measured after initializing at 50 mV. Approximation with g 1 = 1 by fitting the values of α = 0.5 and Q with the voltage data.
Figure 7. Voltage for a CPE with α = 0.5 and Q = 446 Ω 1 Hz 1 2 , measured after initializing at 50 mV. Approximation with g 1 = 1 by fitting the values of α = 0.5 and Q with the voltage data.
Energies 15 00792 g007
Table 1. Values of the parameters for the RC elements.
Table 1. Values of the parameters for the RC elements.
n12345
R n (m Ω )1.04071.99914.524010.523971.3679
τ n (s)0.13521.801212.395170.7794686.9286
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

López-Villanueva, J.A.; Rodríguez Bolívar, S. Constant Phase Element in the Time Domain: The Problem of Initialization. Energies 2022, 15, 792. https://doi.org/10.3390/en15030792

AMA Style

López-Villanueva JA, Rodríguez Bolívar S. Constant Phase Element in the Time Domain: The Problem of Initialization. Energies. 2022; 15(3):792. https://doi.org/10.3390/en15030792

Chicago/Turabian Style

López-Villanueva, Juan Antonio, and Salvador Rodríguez Bolívar. 2022. "Constant Phase Element in the Time Domain: The Problem of Initialization" Energies 15, no. 3: 792. https://doi.org/10.3390/en15030792

APA Style

López-Villanueva, J. A., & Rodríguez Bolívar, S. (2022). Constant Phase Element in the Time Domain: The Problem of Initialization. Energies, 15(3), 792. https://doi.org/10.3390/en15030792

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop