Non-Linear Analytical Model for the Study of Double-Layer Supercapacitors in Different Industrial Uses

: It is generally considered that the representation of a double layer supercapacitor (DLSC) cannot be performed with the usual capacitance and resistance series connected, as it induces a relatively high level of inaccuracy in the results. In multiple previous studies, more advanced models have been developed with very different approaches: models with distributed parameter circuits, based on artiﬁcial neural networks (ANNs), fractional order, etc. A non-linear model, less complex than the previous ones and whose behavior adequately represents the DLSCs, is the one formed by a variable capacitance, dependent on its internal voltage. This paper presents a mathematical study to obtain analytical expressions of all the electrical variables of DLSCs, voltage, current, dissipated power and so on, by means of a previous model. This study is carried out considering that the DLSC is charged and discharged through a voltage source and also discharged through a resistor. In later sections, the operational conditions of the DLSC in numerous industrial applications are presented. Finally, a comparative analysis is made between the results produced by the conventional model, with constant capacitance, and the developed model. This analysis is ﬁnally followed by the conclusions.


Summary of the Available Models for DLSCs
There are multiple models in scientific literature to represent SCs [1]. The most widely used is the RC model, with constant capacitance C connected in series to its equivalent resistance (ESR) [2][3][4][5][6]. This simple model has several advantages; on the one hand, both capacitance and ESR are provided by the manufacturers in the cell data specifications. In addition, both variables are obtained by standardized tests. On the other hand, this model allows relatively simple analytical expressions to be obtained when using SCs in practically all their operating modes, even when charging or discharging at constant power, something that can be much more difficult when working with other more complex models [7]. This electrical model is also used when studying the thermal behavior of the SC [8,9] or when sizing the capacity of an SC bank for a particular application [10,11].
Despite all the above advantages, this is not the most suitable model for the analysis of DLSCs, hence the existence of studies that develop more sophisticated models, based on different mathematical approaches. Depending on the application in which the DLSC is used, it may be necessary to include additional parameters. For example, the model developed by V. Musolino et al. [12] allows for a complete electrical study when the DLSC works with currents containing different frequency components. Furthermore, in order to improve the dynamic response of the final model, it is possible to emulate DLSCs with

Main Modes of Operation of the DLSCs
When working with DLSCs, in most industrial applications, four modes of operation are usual: charging/discharging at a constant current, constant voltage, constant power, constant resistance or combinations of two or more of the above. A summary of the typical applications for each case, emphasizing the charging/discharging through a voltage source or discharging through a constant resistance, which is the most relevant to this study, are presented below.
Full charging of an SC, using only a voltage source, is not usually performed in industrial applications, since, if starting from zero voltage, the charging efficiency cannot exceed 50% [19]. However, the higher the initial voltage of the SC, the higher the efficiency when charging. In fact, if the SC starts charging from a voltage that is 50% of that which provides the source, the charging efficiency rises to 75%. It is therefore common to split the charge into two stages. The first stage is at a constant and high current, until the SC is close to full charge. In the second stage, a voltage source is used, whose unloaded voltage is the SC's rated one. This ensures a full load up to the rated voltage of the SC with very high efficiency. This procedure is known as the constant current-constant voltage method (CC/CV) [20][21][22]. Constant current charging and discharging is also used by cell manufacturers to obtain the electrical parameters of the constant capacity RC series model, by applying standardized methods, and also to obtain the thermal parameters (resistance and thermal capacity).
Recent studies have raised the possibility of integrating SCs into photovoltaic panels, so that the panel itself includes an energy storage system [23][24][25]. This new concept can have many applications for off-grid systems, such as uninterruptible power supply systems, as a backup in case of power failure in mobile phone towers or as off-grid lighting systems. In these applications, where DLSCs are directly connected to solar panels, over a wide range of voltages, during the charging process of the SCs, the panel can be considered as a real current source or as a real voltage source and, in both cases, it can be reduced to the latter.
In many power conversion applications, both the load and the source can be modelled as a system operating at a constant power [26][27][28], for example, in emergency lighting systems and DC micro-grids. In other cases, the load can be represented as a constant value resistor [29][30][31]. N. Kularatna et al. [29] proposed an instantaneous water heating system for domestic use, based on SCs, where a low voltage storage accumulator is needed, for safety reasons, but with the capacity to deliver the stored energy in a very short time. The aim is to heat the cold water remaining in the pipes very quickly to avoid waste. In this case, the SCs are discharged through a heating resistance that can be considered constant. The material used in their manufacturing is usually nichrome, an alloy whose resistivity is practically invariant in a wide range of temperatures. In addition, the heat produced will be quickly transferred to the water, for a very short time, causing the maximum temperature of the material to remain not too high. As the resistance of the heating system is much higher than the internal resistance of the SC bank itself, this application has a very high efficiency; practically all the energy extracted from the SC is dissipated in the heating resistor.
A usual case where SCs can be considered to be discharged through a constant resistance appears when studying the phenomenon of self-discharge. The cells discharge their stored energy through their own leakage resistance, which can be considered constant and of a high value [32]. Another industrial application of SCs in which they are used as a capacitive voltage source is in welding guns [33,34], where a bank of SCs replaces transformers. In this case, the discharge can be considered to take place through a resistor that can be roughly modelled with a constant value.
In the last decade, SCs have started to be used in a very specific application called Supercapacitor Assistant Surge Absorb, SCASA, [35][36][37][38], where they are used to absorb transient overvoltage. The advantage they present in this case is their ability to withstand transient voltage peaks of several kV. To analyze the ability of the SCs to endure these transient surges, the voltage pulse is usually modelled as a high voltage constant source and very short duration.
As has been shown, the number of industrial applications where DLSCs are charged by a voltage source or discharged through a constant resistance is significant and growing.

Limitations of the Present Models
The existing models can be classified into three groups, depending on the possibility to calculate from them simple analytical expressions of the different electrical quantities (current, internal voltage, dissipated energy, etc.), for the four modes of operation described above. The first group would consist of the constant capacitance RC series model. This is the only model for which valid analytical expressions exist to study the four modes of operation. This is the reason why it is the most widely used, even if the results are less accurate than those obtained with more sophisticated approaches. In the second group would be the series RC model, with capacitance linearly varying with voltage and a single branch. For this model, only analytical studies of its operation at a constant current can be found. In fact, the methods for the estimation of its parameters are based on the analytical expressions for DLSC operation when discharged at a constant current. For this model, there are no published analytical expressions when operating at a constant power, nor is it possible to find studies when charging or discharging the cell through a voltage source or discharging at constant resistance. The third group would be formed by the rest of the models (models with several parallel RC branches, distributed parameter models, etc.), which are more accurate but much more complex, making it very difficult, if not impossible, to obtain simple analytical expressions that allow the calculation of the main electrical variables when operating in any of the four operational modes. Therefore, their use is usually restricted to numerical calculation or simulation.
Although numerical methods and the use of simulators are useful and, in certain cases, they are the only way to analyze the behavior of DLSCs, having analytical expressions is a great advantage, since they allow for a simple, fast and accurate analysis. Moreover, conclusions can be drawn much faster from them than with other methods, which imply multiple simulations and laboratory tests.

Objectives of the Study
Due to the interest of analyzing the behavior of a DLSC when it is charged or discharged through a voltage source (or discharged through a constant resistance) and the importance of having analytical expressions, the following section will develop the mathematical formulae needed to calculate all the electrical quantities, using the variable capacitance model. This will fill a small gap in the scientific literature on SCs.
So far, it has been most common to use the constant capacitance RC series model, due to its great simplicity. In this study, the analysis has been extended to a more accurate and complex case, by using the variable capacitance model, which, in addition, is transformed into the constant capacitance model by simply modifying only one parameter. It can be said that this study is an extension of those carried out to date. In addition, a comparison of the results produced by both models will be shown, so that researchers working in this field can select the option that most interests them, depending on what they intend to analyze and the precision they need.
The rest of this paper is arranged as follows: Section 2 will show the mathematical study to obtain the electrical variables of a DLSC connected to a voltage source. The equations obtained will be valid for both charging and discharging. In Section 3, several examples will be presented, where the DLSC will be charged by a voltage source and the solutions obtained will be compared with those of the RC series model and the variable capacitance model. Finally, Section 4 will summarize the conclusions of the whole study.

Electrical Analysis of a DLSC Charged or Discharged with a Voltage Source
This section will show how to express the different electrical variables of a DLSC when connected to a voltage source. The voltage source will be represented by a constant no load voltage, E, and a series connected resistor, R C . As mentioned above, the supercapacitor model will be that of a capacitance linearly varying with its internal voltage. All the obtained functions will be valid for charging and discharging with a voltage source or discharging with a constant resistor. In the second case, it is sufficient to cancel the no-load voltage, E = 0 V, and the discharge resistance will be R C . The internal resistance of the DLSC will be denoted as R and its capacitance will vary with the internal voltage called u. The external voltage of the DLSC is u co and will obviously coincide with the terminal voltage of the voltage source ( Figure 1).

Objectives of the Study
Due to the interest of analyzing the behavior of a DLSC when it is charged or discharged through a voltage source (or discharged through a constant resistance) and the importance of having analytical expressions, the following section will develop the mathematical formulae needed to calculate all the electrical quantities, using the variable capacitance model. This will fill a small gap in the scientific literature on SCs.
So far, it has been most common to use the constant capacitance RC series model, due to its great simplicity. In this study, the analysis has been extended to a more accurate and complex case, by using the variable capacitance model, which, in addition, is transformed into the constant capacitance model by simply modifying only one parameter. It can be said that this study is an extension of those carried out to date. In addition, a comparison of the results produced by both models will be shown, so that researchers working in this field can select the option that most interests them, depending on what they intend to analyze and the precision they need.
The rest of this paper is arranged as follows: Section 2 will show the mathematical study to obtain the electrical variables of a DLSC connected to a voltage source. The equations obtained will be valid for both charging and discharging. In Section 3, several examples will be presented, where the DLSC will be charged by a voltage source and the solutions obtained will be compared with those of the RC series model and the variable capacitance model. Finally, Section 4 will summarize the conclusions of the whole study.

Electrical Analysis of a DLSC Charged or Discharged with a Voltage Source
This section will show how to express the different electrical variables of a DLSC when connected to a voltage source. The voltage source will be represented by a constant no load voltage, E, and a series connected resistor, RC. As mentioned above, the supercapacitor model will be that of a capacitance linearly varying with its internal voltage. All the obtained functions will be valid for charging and discharging with a voltage source or discharging with a constant resistor. In the second case, it is sufficient to cancel the noload voltage, E = 0 V, and the discharge resistance will be RC. The internal resistance of the DLSC will be denoted as R and its capacitance will vary with the internal voltage called u. The external voltage of the DLSC is uco and will obviously coincide with the terminal voltage of the voltage source ( Figure 1). The capacitance of the DLSC shall be modelled as the parallel association of two capacitances, one of constant value, C0, called "Initial Capacitance", which is the capacitance when the DLSC is fully discharged, and the other which is proportional to the internal The capacitance of the DLSC shall be modelled as the parallel association of two capacitances, one of constant value, C 0 , called "Initial Capacitance", which is the capacitance when the DLSC is fully discharged, and the other which is proportional to the internal voltage, C 1 (u). The final capacitance of the cell, C(u), will be the parallel association of the two.
The values of C 0 and k c are not specified in the datasheets of the cells provided by the manufacturers but can be expressed as a function of the rated capacitance, C N , and the rated voltage, U N . In this way, the capacitance can be expressed as a function of the internal voltage as follows: The value of k 0 (dimensionless), which is the quotient of the initial capacitance and the nominal capacitance, is known as the "Normalized Initial Capacitance". In a conventional capacitor, with constant capacitance, it is true that: k 0 = 1, k c = 0 and C N = C 0 . The constant k 0 is also not usually specified in the datasheets, but in most families of SCs from different manufacturers, k 0 presents values between 0.7 and 0.8. Considering a discharge process, according to the definition of an electric current, this is equal to the time variation of charges experienced by the device, i.e., The negative sign indicates that the current represents a discharge. In both conventional SCs and DLSCs, the capacitance establishes the well-known relationship between electric charge and internal voltage: Considering the capacitance value in (1), the stored charge can be expressed as a function of the internal voltage and the values of C 0 and k c : By deriving the charge q with respect to the internal voltage: If the value of the derivative calculated in (6) is substituted in the expression for the current, the relationship between the current and the internal voltage can be obtained: The C V value is called "Virtual Capacitance" or "Dynamic Capacitance". When the DLSC is discharged through a voltage source, the following must also be fulfilled: Therefore, the value of the current will be the following: Substituting in (9) the value of the current obtained in (7) gives the following first-order differential equation of separate variables: Assuming that at t = 0, the DLSC starts from a known internal voltage of value u (t = 0) = U 0 (in the case of the discharge, it is satisfied that U 0 > E), and separating the two variables of (10), it is given that: By integrating (11), the following equality, where time is related to the value of internal voltage, is obtained: During a charging process (E > U 0 ), it will also be true that E ≥ u, since the internal voltage can never exceed the value of the voltage source, and (12) would be equally valid. In order to simplify the mathematical development, two new constants, k 1 and k 2 , will be defined: k 1 is measured in V −1 and k 2 in s −1 . Both constants shall be considered positive. The sign of the voltage source, E, cannot be negative since all SCs, such as electrolytic capacitors, only admit one polarity. If E were negative, the SC could reach an internal voltage of u < 0 which would irreversibly damage it.
By means of k 1 and k 2 , (12) can be simplified as follows: The equality obtained in (15) shows the relationship between the internal voltage, u, and the time, t. From this equation, it is possible to obtain the internal voltage as a function of time, u(t): where W 0 (x) represents the main or superior branch of the Lambert W function which will be used in both charging (E ≥ u) and discharging (E ≤ u), as will be latter explained. E, k 1 and U 0 make it possible to define a new constant, k 3 : k 3 is dimensionless, positive at discharge (U 0 > E) and negative at charge (U 0 < E). The internal voltage, expressed as a function of this new constant, is as follows: In order to further simplify the rest of the mathematical analysis, a new function called g(t), measured in (V), will be defined with the following expression: The internal voltage, u, expressed as a function of g, gives the following equation: The expression of the variable g given in (19) is valid only for k c > 0 (k 0 < 1). To model a conventional capacitor of constant capacity by using (20) and entering the value k 0 = 1 would lead to an indeterminacy. However, it can be shown that the same function would be obtained by using the classical formulas without more than taking into account that: By means of (9) and (20), the discharge current as a function of g is also obtained: In the case of charging (E > u), it follows from (20) that the value of g will be negative, while during discharging, (u > E), g will always be positive. As in (22), the expression for the current supplied by the DLSC was presented, if E > u is satisfied, g will be negative, therefore the current will be absorbed by the DLSC.
As mentioned above, only the main branch of the Lambert W function, W 0 (x), can be used. This restriction is caused by the fact that during a discharge, as time, t, increases, the argument x of the function W 0 (x), presented in (19), is positive (k 3 > 0) and decreasing; for positive values of x, only the main branch of the function is defined (Figure 2). In this case, if time tends to infinity, the value of x will tend to zero; therefore, the value of g also tends to zero, so the internal voltage, u, will converge to the value of the voltage source E. During a charging process, since u < E, the argument of W 0 (x) will be negative (k 3 < 0) and its absolute value will be decreasing. Consequently, as time increases, the absolute value |x| tends to zero, as do the functions W 0 (-|x|) and g, and again, the internal voltage, u, will converge to the value of E, which is what must happen. If the secondary or lower branch of the Lambert W function, W -1 (-|x|), were taken as a solution, as |x| tends to zero, it would tend to -∞, as would the function g and hence the internal voltage of the DLSC, which would not make physical sense. Figure 2 shows the evolution of the variable g over time, both in the charging and discharging processes. In both cases, for very large values of time, t, the variable g will tend to zero, which will be its limiting value. Appl By means of (9) and (20), the discharge current as a function of g is also obtained: In the case of charging (E > u), it follows from (20) that the value of g will be negative, while during discharging, (u > E), g will always be positive. As in (22), the expression for the current supplied by the DLSC was presented, if E > u is satisfied, g will be negative, therefore the current will be absorbed by the DLSC.
As mentioned above, only the main branch of the Lambert W function, W0(x), can be used. This restriction is caused by the fact that during a discharge, as time, t, increases, the argument x of the function W0(x), presented in (19), is positive (k3 > 0) and decreasing; for positive values of x, only the main branch of the function is defined (Figure 2). In this case, if time tends to infinity, the value of x will tend to zero; therefore, the value of g also tends to zero, so the internal voltage, u, will converge to the value of the voltage source E. During a charging process, since u < E, the argument of W0(x) will be negative (k3 < 0) and its absolute value will be decreasing. Consequently, as time increases, the absolute value |x| tends to zero, as do the functions W0(-|x|) and g, and again, the internal voltage, u, will converge to the value of E, which is what must happen. If the secondary or lower branch of the Lambert W function, W-1(-|x|), were taken as a solution, as |x| tends to zero, it would tend to -∞, as would the function g and hence the internal voltage of the DLSC, which would not make physical sense. Figure 2 shows the evolution of the variable g over time, both in the charging and discharging processes. In both cases, for very large values of time, t, the variable g will tend to zero, which will be its limiting value. On the other hand, from Equations (8) and (22), it is easy to deduce the value of the external voltage as a function of the variable g: The power dissipated as heat in the internal resistance R of the DLSC, pd, can also be On the other hand, from Equations (8) and (22), it is easy to deduce the value of the external voltage as a function of the variable g: The power dissipated as heat in the internal resistance R of the DLSC, p d , can also be expressed as a function of the variable g in the following way: Similarly, the power dissipated in the internal resistance of the voltage source, R c , presents the same function as p d , changing R to R c ; that is: The power absorbed by the voltage source, E, representing the no-load voltage of the voltage source is: The output power of the DLSC, i.e., the power consumed by the voltage source, p out , is obtained by simply multiplying the external voltage by the current; therefore: The energy dissipated in the internal resistance, R, of the DLSC, e d , is obtained by integrating the power loss, p d , obtained in (24). In doing so, it is considered that: Deriving g(t) as a function of time, as presented in (19), after some simplifications, it is obtained that: The expression for g' in (29) is valid for both charging and discharging, without the need to change the sign. If the value g' obtained in (29) and the value of p d calculated in (24) are substituted into (28), integrating both members of the obtained equation, the result is the value of the energy dissipated in the resistance R, e d , depending on the variable g: In the expression given in (30), g 0 is the value of g evaluated in t = 0 that, according to (20), results in g 0 = U 0 − E. The energy dissipated in R c , e dRc , follows the same evolution as e d , by only changing the value of R that appears in the numerators by R c , i.e., Another variable that is important to consider is the instantaneous energy stored in the DLSC. This energy depends on the constants C 0 , k c and the internal voltage, u, according to the following equation: The C E value is called "Energetic Capacitance", as it is used to obtain the energy stored in a variable capacitance DLSC, with an equation similar to that used in a conventional capacitor. This, as with the dynamic capacitance, also depends on the internal voltage. By substituting the former (20) in (32), the instantaneous stored energy as a function of g is finally obtained: If the DLSC starts from an initial internal voltage, U 0 , and after a time, t, reaches an internal voltage, u, the discharged energy, e dch , will be: On the other hand, the energy consumed by the voltage source, E, e E , can be also obtained by integration, bearing in mind that: Replacing in (35) the value of p E as a function of g, obtained in (26), and the value of the derivative of g(t) with respect to time, calculated in (29), and integrating the energy consumed or generated, depending on whether it is a charge or discharge of the DLSC can be obtained: Finally, it is necessary to develop an equation that expresses time as a function of g. The variable g as a function of time was calculated earlier in (19). From the definition of the Lambert W function itself, it must be satisfied that: Taking neperian logarithms on both sides of (37), the function t(g) is obtained: By replacing k 1 , k 2 and k 3 , previously and respectively calculated in (13), (14) and (17) in (38), the time, t, can be expressed as a function of g and other known variables: As can be deduced from (39), for g to reach a zero value, the time, t, should tend to infinity. As in first-order circuits with conventional capacitors that are charged or discharged through a voltage source, the time constant, τ, can be defined as the elapsed time for the following equality to be satisfied, valid for both charging and discharging (e being the Euler's number): Considering the value of t obtained in (39) and the relationship between g, u and E from (20), the "time constant" of the circuit, τ, can be determined by the following equation: In contrast to conventional capacitors, the time constant is not only dependent on the passive elements but is also a function of the initial voltage of the DLSC, U 0 and the no-load voltage, E, of the source. This result is congruent, since for this model, the capacitance of the DLSC is a function with linear dependence on its internal voltage, whereas in a conventional capacitor, the capacitance remains invariant.

Mathematical Expresión of g as a Function of the Electrical Variables
Once a set of equations has been obtained with the main electrical variables and time as a function of the variable g, the complementary equations of g, as a function of the current, voltage, power and so on, will be calculated. In this way, the development of analytical equations will be completed, obtaining a set of functions where any variable can be expressed as a function of the others. The variable g as a function of time has already been found in Equation (19). The calculation of g as a function of internal voltage is immediate if (20) is used: The calculation of g as a function of the current is straightforward from (22): Similarly, g can also be obtained as a function of the external voltage by simply using (23): By means of (24), g can be expressed as a function of the power dissipated in the internal resistance of the DLSC, p d .
As in (45), the value of g as a function of the power dissipated in the resistance of the voltage source R c , p dRc , has the same form but replaces the value of R by R c ; thus: The expression for g as a function of the power consumed by E, p E , is also directly obtained from (26): To express g as a function of the DLSC output power, p out , it is necessary to solve the second-degree equation shown in (27). Of the two possible solutions, the positive one must be chosen both in the charging and discharging cases: The equation of g as a function of the energy dissipated by the source, E, e E , is obtained in almost the same way, by solving the second-degree equation obtained in (36) and also choosing the positive solution: Obtaining expressions of g as a function of e dch , e stored , e d and e dRc is more complex, since all four are cubic functions of g. They have therefore been omitted from this study, since, in addition, the set of analytical expressions already shown and their combinations allow practically all possible alternatives to be covered. For example, if the value of the internal voltage of the DLSC, u, as a function of the power consumed by the voltage source E, p E , is to be calculated, it would be sufficient to combine (20), where u(g) is defined, with (47), which presents the equation of g(p E ). The equation of u(p E ), would be directly obtained: If the objective was to calculate the time it would take for the external voltage, u co , to reach a certain value, i.e., to obtain the function t(u co ), the combination of (39), which defines t(g), and (44) which shows g(u co ), would directly yield t(u co ): As in the previous two examples, the process can be repeated with any two electrical variables (it is possible to include time, t) and to obtain an analytical expression of the associated function.

Results and Discussion
This section will show a case study in which a DLSC, modelled as a capacitance linearly varying with the internal voltage, will be charged through a voltage source. The evolution of the electrical variables will be studied using the set of analytical expressions of Section 2, and simulations will be carried out to analyze their variation with different values of the constant k 0 . Subsequently, the results obtained will be compared with the same DLSC, considering it as a constant capacitance with its equivalent series resistance. Table 1 shows its rated data together with the no-load voltage of the voltage source, E.
The wiring diagram of the case study is presented in Figure 3. To study the time evolution of the different variables, three values will be considered for the variable k 0 = [0.65, 0.85, 1] and four different values will be chosen for the internal resistance of the voltage source, R c (Ω) = [0. 5,1,3,5]. From the obtained results, the constants C 0 , k c , k 1 , k 2 and k 3 as defined in Equations (2), (13), (14) and (17) will be, respectively, calculated together with the value of the time constant, τ, according to (41), which, in this study, is the time taken for the internal voltage to reach 1.7067 V. All the above values are shown in Table 2.  The wiring diagram of the case study is presented in Figure 3. To study the time evolution of the different variables, three values will be considered for the variable k0 = [0.65, 0.85, 1] and four different values will be chosen for the internal resistance of the voltage source, Rc (Ω) = [0. 5,1,3,5]. From the obtained results, the constants C0, kc, k1, k2 and k3 as defined in Equations (2), (13), (14) and (17) will be, respectively, calculated together with the value of the time constant, τ, according to (41), which, in this study, is the time taken for the internal voltage to reach 1.7067 V. All the above values are shown in Table 2. Figure 4 shows the evolution of the DLSC internal voltage, u(t), for different values of the constant k0 and various values of the internal resistance of the voltage source, Rc. As can be appreciated from it, for the three models of different k0, there is a coincidence in three values of the internal voltage, u(t). Two of them are obvious, since they are the initial value (they all start from the same U0) and the final value, which converge to the no-load voltage of the voltage source E = UN = 2.7 V. The third one corresponds to a specific value of the g variable. This value, which will be denoted as gs, ensures that whatever the value of k0 and the resistance, Rc, it will be the same in all models. The value of gs can be calculated as follows: (52) Figure 3. Wiring diagram of the case study. Table 2. Constants C 0 , k c , k 1 , k 2 , k 3 and τ for different values of k 0 and R c . As can be appreciated from it, for the three models of different k 0 , there is a coincidence in three values of the internal voltage, u(t). Two of them are obvious, since they are the initial value (they all start from the same U 0 ) and the final value, which converge to the no-load voltage of the voltage source E = U N = 2.7 V. The third one corresponds to a specific value of the g variable. This value, which will be denoted as g s , ensures that whatever the value of k 0 and the resistance, R c , it will be the same in all models. The value of g s can be calculated as follows: As shown in (52), g s only depends on E, the rated voltage of the cell, U N and its initial voltage, U 0 , but is independent of all other variables, including R c . In the case study presented, g s = −0.5486 V. If this value of g s is replaced in (20), where u(g) was calculated, the internal voltage at which all the models coincide is 2.1514 V, as can be appreciated from Figure 4. Below this voltage, the variable capacitance model (k 0 < 1) produces higher internal voltages than the conventional model, while beyond this point, the trend is reversed. By substituting this value of g s in (39), the result is the time instant at which, for any value of k 0 , all the internal voltages cross at a common point. This instant does depend on the value of R c . Figure 4. Below this voltage, the variable capacitance model (k0 < 1) produces higher internal voltages than the conventional model, while beyond this point, the trend is reversed. By substituting this value of gs in (39), the result is the time instant at which, for any value of k0, all the internal voltages cross at a common point. This instant does depend on the value of Rc.  Table 3 shows the times at which, with all the models, the same internal voltage is reached as a function of the values of Rc. The value of the current at those instants is also displayed.  Table 3 shows the times at which, with all the models, the same internal voltage is reached as a function of the values of R c . The value of the current at those instants is also displayed.  Figure 5 shows the difference between the calculation of the internal voltage with variable capacitance, k 0 = [0.65, 0.85], and the conventional model with constant capacitance, k 0 = 1, for different values of R c . As can be seen, the maximum difference is independent of the value taken by R c , which only affects the instant at which it appears. In this example, the maximum discrepancy between the conventional DLSC model and the one with k 0 = 0.65 is 0.1738 V. Figure 6 shows the evolution of the instantaneous current of the circuit, i(t), calculated with the three values of the constant k 0 and for different values of R c . In the same instants in which the internal voltage, u, coincides, the current also coincides, whatever the value of k 0 , although in this case, the values of the current are different for each Rc value, since i(t) depends on this resistance, as was earlier indicated by means of expression (9). It can be observed how, in the calculation of the current, the difference between modelling the DLSC with variable capacitance and with constant capacitance is minimal. It can therefore be concluded that the usual RC series model is perfectly valid for i(t) calculation. lated with the three values of the constant k0 and for different values of Rc. In the same instants in which the internal voltage, u, coincides, the current also coincides, whatever the value of k0, although in this case, the values of the current are different for each Rc value, since i(t) depends on this resistance, as was earlier indicated by means of expression (9). It can be observed how, in the calculation of the current, the difference between modelling the DLSC with variable capacitance and with constant capacitance is minimal. It can therefore be concluded that the usual RC series model is perfectly valid for i(t) calculation.     Figure 7 shows the difference between the current obtained with variable capacitance, k 0 = [0.65, 0.85], and the conventional model with constant capacitance, k 0 = 1, for different values of R c . As can be appreciated, the higher the value of R c , the less difference there is between considering the model with variable capacitance and the RC series model with constant capacitance. Figure 8 shows the time evolution of the energy dissipated in the internal resistance of the DLSC, e d (t), presented in (30) and calculated with the three values of the constant k 0 and different values of R c . The higher the value of the constant k 0 , at any instant of time, the higher the energy dissipated. Furthermore, the larger the value of R c the smaller the difference between models with different k 0 values. Finally, Figure 9 shows the difference in the energy dissipated in the DLSC, e d (t), displayed in absolute value, between the model with variable capacitance, k 0 = [0.65, 0.85], and conventional capacitor, k 0 = 1, considering different values of R c .  Figure 7 shows the difference between the current obtained with variable capacitance, k0 = [0.65, 0.85], and the conventional model with constant capacitance, k0 = 1, for different values of Rc. As can be appreciated, the higher the value of Rc, the less difference there is between considering the model with variable capacitance and the RC series model with constant capacitance.

Conclusions
The main electrical variables of a DLSC, modelled as a capacitance linearly varying with the internal voltage, when charged or discharged through a voltage source, have been obtained. The set of analytical expressions presented is equally useful for the case in which the DLSC is discharged through a constant resistance, as occurs in industrial applications such as instantaneous heaters, or if it is desired to analyze the self-discharge pro-

Conclusions
The main electrical variables of a DLSC, modelled as a capacitance linearly varying with the internal voltage, when charged or discharged through a voltage source, have been obtained. The set of analytical expressions presented is equally useful for the case in which the DLSC is discharged through a constant resistance, as occurs in industrial applications such as instantaneous heaters, or if it is desired to analyze the self-discharge pro-

Conclusions
The main electrical variables of a DLSC, modelled as a capacitance linearly varying with the internal voltage, when charged or discharged through a voltage source, have been obtained. The set of analytical expressions presented is equally useful for the case in which the DLSC is discharged through a constant resistance, as occurs in industrial applications such as instantaneous heaters, or if it is desired to analyze the self-discharge process. Furthermore, from the new set of equations, the time constant of the DLSC has been defined by resemblance to conventional, constant capacitance capacitors. Unlike other more sophisticated models, which are valuable and very accurate, even in the representation of the physical structure of the DLSC, the set of analytical equations presented here is not restricted to the use of numerical calculation tools. It also allows analytical solutions to be easily obtained without the complexity involved in the models mentioned, where complex tools, from distributed parameter circuits to artificial neural networks, are used.
It is often claimed that the series RC model does not accurately model the behavior of a DLSC. Therefore, in this paper, a comparative study between the model with variable capacitance (k 0 < 1) and the conventional model with constant capacitance (k 0 = 1) has been carried out and several conclusions have been drawn. It can be finally concluded that when a DLSC is connected to a voltage source, the differences between the models with k 0 < 1 and with k 0 = 1 are not too relevant, and the traditional RC series model, which is simpler and more widely used, provides, in most cases, sufficiently good values so that it is not necessary to treat it as a variable capacitance SC. In the case of needing to improve accuracy and/or perform direct analytical calculations, the new set of analytical equations developed here can be applied in a straightforward and very simple way.