A Fractional-Order Kinetic Battery Model of Lithium-Ion Batteries Considering a Nonlinear Capacity

Accurate battery models are integral to the battery management system and safe operation of electric vehicles. Few investigations have been conducted on the influence of current rate (C-rate) on the available capacity of the battery, for example, the kinetic battery model (KiBaM). However, the nonlinear characteristics of lithium-ion batteries (LIBs) are closer to a fractional-order dynamic system because of their electrochemical materials and properties. The application of fractional-order models to represent physical systems is timely and interesting. In this paper, a novel fractional-order KiBaM (FO-KiBaM) is proposed. The available capacity of a ternary LIB module is tested at different C-rates, and its parameter identifications are achieved by the experimental data. The results showed that the estimated errors of available capacity in the proposed FO-KiBaM were low over a wide applied current range, specifically, the mean absolute error was only 1.91%.


Introduction
Electric vehicles have the advantages of high fuel economy and zero exhaust emissions [1,2].As clean and efficient energy sources, power batteries are core components and are critical to the comprehensive performance of vehicles.Lithium-ion batteries (LIBs) show strong overall advantages in the field of power batteries because of their high energy density, long life, and excellent cycle performance [2,3].
Battery states mainly include the state of charge (SOC), state of health (SOH), state of power (SOP), state of energy (SOE), and state of function (SOF) [4][5][6].However, they cannot be measured directly; they can only be estimated by testing battery voltage, current, and temperature, among other factors.Battery state estimation is extremely important in battery management systems (BMSs) to guarantee the safe and reliable operation of batteries, and a multitude of research has investigated the methods of state estimation of LIBs based on an accurate model [7][8][9][10][11][12][13][14].Commonly used battery models include electrochemical models (EchMs) [11,14], analytical models (AMs) [7, [12][13][14], stochastic models (SMs) [7, 12,14], neural network models (NNMs) [14], and equivalent circuit models (ECMs) [8][9][10][11]14].EchMs are accurate in describing the internal electrochemical reaction using complex, nonlinear differential equations, but they are difficult to understand.AMs model the major properties of the battery using only a few equations, and are much easier to use than EchMs.SMs mainly concern the battery recovery characteristics as a Markov process, in which the pulse discharge characteristics of the battery can be described, but they are not applicable for variable current.NNMs have fast parallel processing capabilities as well as strong self-learning and self-organizing abilities, but they require a large amount of training data, and errors can arise from the training data and training methods.ECMs are widely used for electrical design and modelling simulations [9,10] because they can accurately describe the battery voltage-current (U-I) performance.Many improved models have been proposed by scholars that not only describe U-I performance (external characteristics) but also the capacity performance (internal characteristics).For instance, an ECM with a variable effective capacity for LIBs is proposed in [15], and the model is further optimized using computational intelligence techniques [16].
However, an accurate and concise battery model is not readily achievable because the models are highly nonlinear and complex.Few investigations have been conducted on the influence of the current rate (C-rate) on the available capacity of batteries because the characteristics of LIBs will change significantly under different conditions in view of its sensitivity to C-rate, temperature, cycle life, etc.As described in detail in Figure 1, the battery capacity is not similar to water in a bucket.It is mainly manifested as "capacity nonlinear effect" and "recovery effect" [7, 11,12].The capacity nonlinear effect is that the available capacity will decline nonlinearly with C-rates; the greater the discharge current, the less the available capacity.The recovery effect is that the battery's available capacity will rise up when discharge is stopped.For instance, a released capacity of 50 A will be much less than that at 5 A, and the capacity will be restored if the battery rests for a while.
can accurately describe the battery voltage-current (U-I) performance.Many improved models have been proposed by scholars that not only describe U-I performance (external characteristics) but also the capacity performance (internal characteristics).For instance, an ECM with a variable effective capacity for LIBs is proposed in [15], and the model is further optimized using computational intelligence techniques [16].
However, an accurate and concise battery model is not readily achievable because the models are highly nonlinear and complex.Few investigations have been conducted on the influence of the current rate (C-rate) on the available capacity of batteries because the characteristics of LIBs will change significantly under different conditions in view of its sensitivity to C-rate, temperature, cycle life, etc.As described in detail in Figure 1, the battery capacity is not similar to water in a bucket.It is mainly manifested as "capacity nonlinear effect" and "recovery effect" [7, 11,12].The capacity nonlinear effect is that the available capacity will decline nonlinearly with C-rates; the greater the discharge current, the less the available capacity.The recovery effect is that the battery's available capacity will rise up when discharge is stopped.For instance, a released capacity of 50 A will be much less than that at 5 A, and the capacity will be restored if the battery rests for a while.The remaining available battery capacity is critical in electric vehicles, similar to the role of the remaining fuel in internal combustion vehicles.Thus, the available capacity estimation of a battery considering C-rate is very important.Overall, commonly used models for available capacity estimation of the battery mainly include two classic analytical models: Peukert's law and the kinetic battery model (KiBaM).Proposed by Wilhelm Peukert, Peukert's law was used to estimate nonlinear delivered capacity and predict the battery run time of a rechargeable lead-acid battery at different constant discharge C-rates from the fully charged state [17].The nonlinear properties between the available capacity and the C-rate are considered, and the battery's run time can be approximated.The model itself is relatively simple compared to KiBaM.However, it does not consider the recovery effect of the battery when discharge stops.Fortunately, Manwell and McGowan proposed a kinetic battery model (KiBaM) to model lead-acid storage batteries in 1993 [13].It is intuitionistic and easy to understand based on perceptual knowledge.Moreover, it can be used in modelling and simulation [7, [12][13][14].In [14], a widely-used KiBaM was used to capture nonlinear capacity effects for accurate SOC tracking and runtime predictions of the battery.In [18], KiBaM was extended to consider the temperature effect on battery capacity.The proposed temperature-dependent KiBaM (T-KiBaM) can handle operating temperatures, and it can provide better estimates for battery lifetimes and voltage behaviors.However, the classic KiBaM is described by regular calculus.A precise and concise battery The remaining available battery capacity is critical in electric vehicles, similar to the role of the remaining fuel in internal combustion vehicles.Thus, the available capacity estimation of a battery considering C-rate is very important.Overall, commonly used models for available capacity estimation of the battery mainly include two classic analytical models: Peukert's law and the kinetic battery model (KiBaM).Proposed by Wilhelm Peukert, Peukert's law was used to estimate nonlinear delivered capacity and predict the battery run time of a rechargeable lead-acid battery at different constant discharge C-rates from the fully charged state [17].The nonlinear properties between the available capacity and the C-rate are considered, and the battery's run time can be approximated.The model itself is relatively simple compared to KiBaM.However, it does not consider the recovery effect of the battery when discharge stops.Fortunately, Manwell and McGowan proposed a kinetic battery model (KiBaM) to model lead-acid storage batteries in 1993 [13].It is intuitionistic and easy to understand based on perceptual knowledge.Moreover, it can be used in modelling and simulation [7,[12][13][14].In [14], a widely-used KiBaM was used to capture nonlinear capacity effects for accurate SOC tracking and runtime predictions of the battery.In [18], KiBaM was extended to consider the temperature effect on battery capacity.The proposed temperature-dependent KiBaM (T-KiBaM) can handle operating temperatures, and it can provide better estimates for battery lifetimes and voltage behaviors.However, the classic KiBaM is described by regular calculus.A precise and concise battery model at various conditions has always been challenging for researchers to create.Using fractional calculus with impedance models is quite common in the modelling of energy storage and generation elements, including capacitors/super capacitors [19][20][21] and batteries [19,22,23]; what is more, the fractional calculus has also been used in the state estimation and prediction of batteries [21,23,24].Actually, the nonlinear characteristics of lithium-ion batteries are closer to a fractional-order dynamic system, because the material diffusion and electrochemical properties have been successfully described using fractional calculus.The application of fractional-order models to represent physical systems is timely and interesting.

Capacity
In this paper, a novel fractional-order KiBaM (FO-KiBaM) is proposed to describe the nonlinear capacity characteristics of LIBs.The research ideas and arrangement of the rest of the paper are as follows.Firstly, the capacity nonlinear effect and recovery effect of KiBaM are analyzed in Section 2. The basic principle of fractional calculus and its application in the proposed FO-KiBaM are introduced in Section 3. In Section 4 the charge and discharge experiments of a LIB module are designed and conducted under different C-rates.Finally, the results and model error analyses are compared and illustrated in Section 5, followed by the conclusion in Section 6.

Kinetic Battery Model (KiBaM)
As shown in Figure 2, KiBaM uses two wells of different sizes to describe the dynamic changes in battery capacity.Thus, it is also called a "two-well" model.The two wells represent the "directly available capacity" and the "temporary capacity" of the battery, respectively [13,14].The directly available capacity can be obtained directly at discharge, denoted as y 1 , and its height is denoted as h 1 .The "temporary capacity" cannot be directly obtained at discharge, denoted as y 2 , and its height is denoted as h 2 .It is easy to understand that the sum of y 1 and y 2 is the total capacity of the battery, and the sum of y 1 and a part of y 2 is the available capacity of the battery.The letter k represents the rate that the charge flow from y 2 into y 1 ; and the letter c represents the capacity proportion of the two wells.The k and c variables affect the nonlinear battery capacity characteristics; notably, they are closely related to the materials and composition of the battery.
Electronics 2019, 8, x FOR PEER REVIEW 3 of 17 calculus with impedance models is quite common in the modelling of energy storage and generation elements, including capacitors/super capacitors [19][20][21] and batteries [19,22,23]; what is more, the fractional calculus has also been used in the state estimation and prediction of batteries [21,23,24].Actually, the nonlinear characteristics of lithium-ion batteries are closer to a fractional-order dynamic system, because the material diffusion and electrochemical properties have been successfully described using fractional calculus.The application of fractional-order models to represent physical systems is timely and interesting.
In this paper, a novel fractional-order KiBaM (FO-KiBaM) is proposed to describe the nonlinear capacity characteristics of LIBs.The research ideas and arrangement of the rest of the paper are as follows.Firstly, the capacity nonlinear effect and recovery effect of KiBaM are analyzed in Section 2. The basic principle of fractional calculus and its application in the proposed FO-KiBaM are introduced in Section 3. In Section 4 the charge and discharge experiments of a LIB module are designed and conducted under different C-rates.Finally, the results and model error analyses are compared and illustrated in Section 5, followed by the conclusion in Section 6.

Kinetic Battery Model (KiBaM)
As shown in Figure 2, KiBaM uses two wells of different sizes to describe the dynamic changes in battery capacity.Thus, it is also called a "two-well" model.The two wells represent the "directly available capacity" and the "temporary capacity" of the battery, respectively [13,14].The directly available capacity can be obtained directly at discharge, denoted as y1, and its height is denoted as h1.
The "temporary capacity" cannot be directly obtained at discharge, denoted as y2, and its height is denoted as h2.It is easy to understand that the sum of y1 and y2 is the total capacity of the battery, and the sum of y1 and a part of y2 is the available capacity of the battery.The letter k represents the rate that the charge flow from y2 into y1; and the letter c represents the capacity proportion of the two wells.The k and c variables affect the nonlinear battery capacity characteristics; notably, they are closely related to the materials and composition of the battery.
It is clear that the directly available capacity y1, the temporary capacity y2, and the heights h1 and h2 in the KiBaM are satisfied by the following: It is clear that the directly available capacity y 1 , the temporary capacity y 2 , and the heights h 1 and h 2 in the KiBaM are satisfied by the following: where δ h represents the height difference between the two wells.The principle of the KiBaM is as follows.When the battery is discharged, the charge of y 1 flows out, simultaneously, the charge of y 2 flows into y 1 slowly with k.The charge flowing out of y 1 is faster than flowing in, so the height difference between y 1 and y 2 will increase.The larger the discharge current, the less capacity that is released, which reflects the nonlinear effect of the battery's capacity.Additionally, the unavailable capacity of the battery is the capacity represented by the height differences of y 1 and y 2 .Once discharging is stopped, the charge of y 2 will flow into y 1 slowly until the heights of y 1 and y 2 are equal, and the charge of y 1 will pick up, which reflects the battery's recovery effect.From the above intuitive graphical description of KiBaM, the changes of the charges y 1 and y 2 in the two wells can be expressed as follows: dy 2 (t) From here, the capacity variance of the battery with time can be expressed as follows: where C rem , C unav , and C av represent the remaining capacity of battery, the unavailable capacity of battery, and the available capacity of battery, respectively.Take the constant-current discharge for example.As shown in Figure 3, the time course can be described as follows: where t 0 , t d , t r represents the initial time, the discharge end time, and the recovery end time, respectively.
Electronics 2019, 8, x FOR PEER REVIEW 4 of 17 than flowing in, so the height difference between y1 and y2 will increase.The larger the discharge current, the less capacity that is released, which reflects the nonlinear effect of the battery's capacity.Additionally, the unavailable capacity of the battery is the capacity represented by the height differences of y1 and y2.Once discharging is stopped, the charge of y2 will flow into y1 slowly until the heights of y1 and y2 are equal, and the charge of y1 will pick up, which reflects the battery's recovery effect.From the above intuitive graphical description of KiBaM, the changes of the charges y1 and y2 in the two wells can be expressed as follows: From here, the capacity variance of the battery with time can be expressed as follows: where Crem, Cunav, and Cav represent the remaining capacity of battery, the unavailable capacity of battery, and the available capacity of battery, respectively.
Take the constant-current discharge for example.As shown in Figure 3, the time course can be described as follows: , where 0 d r t t t , , represents the initial time, the discharge end time, and the recovery end time, respectively.Actually, when the battery is discharged in a time-varying current, the entire process can be divided into multiple discharge segments in accordance with the time interval.In each segment, the battery can be considered discharged with a constant current and then set aside for a while.Actually, when the battery is discharged in a time-varying current, the entire process can be divided into multiple discharge segments in accordance with the time interval.In each segment, the battery can be considered discharged with a constant current and then set aside for a while.
By the Laplace transform and inverse Laplace transform of (2), the following expression can be obtained: where k = k c(1−c) .The unavailable capacity of battery in (3) can be derived as follows: It can be further expressed as follows: As can be seen from ( 7), in the discharge period t 0 ≤ t ≤ t d , the battery's unavailable capacity is impacted nonlinearly by the discharge time; the longer the time and the larger the current, the larger the unavailable capacity.In the stationary period t d < t ≤ t r , the unavailable capacity decreases with time, because the charge of y 2 flows into y 1 , which reflects the recovery effect of the battery in an open circuit state.
It can be seen that when y 1 = 0 (or h 1 = 0), the battery is discharged completely.At this point, all the remaining capacity is unavailable.Therefore, judging whether the battery is fully discharged is shown as follows: The remaining available capacity of the battery can be expressed as follows: If the initial conditions t 0 = 0, and the battery is discharged at a constant current I, the initial state of y 1 and y 2 of the battery are shown as below: where C 0 is the initial total capacity of battery.Then, Formula (3) can be simplified as follows: 3. The Proposed Fractional-Order KiBaM (FO-KiBaM)

Fractional Calculus Theory
Fractional calculus is not a new concept, in fact, it can be traced back to the discussion of Leibniz and Hospital's research work in 1695, but due to a variety of calculation difficulties in practical applications, it was only a purely theoretical exploration in the early days.However, fractional calculus attracts great attention in complex engineering applications because of the continuous developments in natural science and computer technology.In recent decades, theoretical and mathematical tools have been used in research in a multitude of disciplines, and they are especially successful in high-energy physics, fluid mechanics, viscoelastic material mechanics, anomalous diffusion, electronic components analysis, and system control [20,25].Fractional calculus has been a research hotspot for its unique and irreplaceable advantages.
Fractional derivative definitions (FDD) are defined in different ways; the most commonly used include the Grunwald-Letnikov definition (GL-FDD), the Riemann-Liouville definition (RL-FDD), and the Caputo definition (Caputo-FDD) [20,25].The GL-FDD is expressed as: where G D α t represents the GL-FDD type; f (t) is an arbitrary integrable function; α is an arbitrary real number; n = [ t−t 0 h ] represents the integer part; and In fact, the RL-FDD is obtained on the basis of GL-FDD by simplifying the calculation process.RL-FDD can be expressed as follows: where R D α t represents the RL-FDD type; n is an integer; and Γ(•) is the Gamma function, a commonly used basic functions in fractional calculus, defined as follows: The function Γ(•) has the following properties: And the Laplace transform will be established by: where H(x) denotes the unit step function, which implies that it only needs x ≥ 0. The formula and its inverse transformation are often used in fractional calculus.
The Caputo-FDD is expressed as follows: where C D α t represents the Caputo-FDD type.The Laplace transform of the Caputo-FDD is expressed as follows: Thus, the Laplace transform of the Caputo-FDD under the zero initial conditions is: The derivation of constants in the Caputo-FDD is bounded, while the derivation of constants in the RL-FDD is unbounded.The RL-FDD needs to solve an initial value problem that it is theoretically feasible but lacks physical meaning.Therefore, the Caputo-FDD is more suitable for solving the initial value problem of fractional calculus.Thus, it was adopted in this paper.

Fractional-Order KiBaM
The internal electrochemical reaction of a power battery is extremely complex.The strong nonlinear characteristics of LIBs shows a fractional-order dynamic behavior [21][22][23].Therefore, fractional calculus can be used to model a novel fractional-order KiBaM (FO-KiBaM) with a higher accuracy [26].Fractional derivatives can be replaced to describe the battery capacity change process in (2).
where α is the order of fractional derivative equations, and 0 < α < 1.
In the case when the initial time t 0 = 0, and the battery is discharged at a constant current I, the Laplace transform of (20) will be expressed as: Similarly, assuming that the initial state of y 1 and y 2 are the same as in (10), the following will be obtained: The above inverse Laplace transform of the fractional calculus transfer function can be obtained using the Mittag-Leffler function, a commonly used basic function in fractional calculus [27][28][29][30].The Mittag-Leffler function has two different definition forms: the single-parameter form and the two-parameter form.The two-parameter form of the Mittag-Leffler function is defined as shown below: If β = 1, the single-parameter form of the Mittag-Leffler function will be obtained: The exponential function e z is critical in calculus.Similarly, the Mittag-Leffler function equally plays an important role, and it appears frequently in solutions of fractional differential equations.Sometimes the Mittag-Leffler function with two parameters is also called a generalized exponential function.In fact, e z can be seen as a special case of the Mittag-Leffler function, because they are equal if α = 1.
To facilitate the inverse Laplace transform, we defined a new function as shown below: Its Laplace transform can be obtained as follows: According to the above properties of (22), it is clear that β = α + 1 in the transfer function (27), and its inverse Laplace transform will be obtained as shown below: Similarly, the height difference of two wells can be obtained: Therefore, the capacity of the battery with fractional calculus can be expressed as follows: The SOC of the battery is similar to the fuel gauge of conventional internal combustion vehicles, which can be used to estimate the distance the vehicle can travel.The definition of SOC is: where C rem , C max represents the remaining capacity and the maximum available capacity of battery.The unavailable capacity of the battery is not considered in this definition.Thus, it cannot tell the driver the actual available battery capacity at different C-rates.To predict the remaining mileage of electric vehicles more accurately, it is necessary to improve the SOC definition.An improved SOC definition is shown as below: This definition of SOC takes into account the unavailable capacity of battery.Thus, it is useful in estimating the run time of the battery, and it can predict the remaining mileage of electric vehicles more accurately to relieve the "range anxiety" for the drivers.At the same time, by a more accurate definition of SOC, it can help to determine effective management strategies to avoid overcharging and over discharging the battery.
Further, we found that the height of y 1 in the proposed FO-KiBaM can be expressed as follows: As can be seen from ( 32) and ( 33), the height h 1 of the left "well" reflects the change in the battery's remaining available capacity, and it also explains why the height h 1 is an intuitive representation of the battery SOC.

Experimental Platform and Test Results
As shown in Figure 4, the platform of the battery test system consisted of a battery charging and discharging cycler MKLtech MCT8-100-05Z, a programmable temperature chamber, a LiNiMnCoO ternary LIB module, and a computer with control software.The voltage and the current range of the battery cycler was 0-5 V and ±100 A, respectively.The temperature of the chamber was maintained at 30 • C. The sampling rate of voltage and current was set at 1 Hz.The tested ternary LIB module had a capacity of 32.50 Ah, which consisted of 13 battery cells connected in parallel, as shown in Figure 5. Table 1 shows the specifications of the ternary LIB cells.The available capacity tests at different C-rates of 0.2, 0.67, 1, 1.5, 2, and 3 C were carried out.Here, 1 C indicated that the rechargeable battery was continuously discharged for 1 h at a current equal to the battery's nominal capacity.
As can be seen from ( 32) and ( 33), the height h1 of the left "well" reflects the change in the battery's remaining available capacity, and it also explains why the height h1 is an intuitive representation of the battery SOC.

Experimental Platform and Test Results
As shown in Figure 4, the platform of the battery test system consisted of a battery charging and discharging cycler MKLtech MCT8-100-05Z, a programmable temperature chamber, a LiNiMnCoO ternary LIB module, and a computer with control software.The voltage and the current range of the battery cycler was 0-5 V and ±100 A, respectively.The temperature of the chamber was maintained at 30 °C.The sampling rate of voltage and current was set at 1 Hz.The tested ternary LIB module had a capacity of 32.50 Ah, which consisted of 13 battery cells connected in parallel, as shown in Figure 5. Table 1 shows the specifications of the ternary LIB cells.The available capacity tests at different Crates of 0.2, 0.67, 1, 1.5, 2, and 3 C were carried out.Here, 1 C indicated that the rechargeable battery was continuously discharged for 1 h at a current equal to the battery's nominal capacity.In the test, first the ternary battery module was fully charged by the constant current, constant voltage (CCCV) charging strategy.The constant current (CC) was 10 A, and the constant voltage (CV) was 4.2 V.The cutoff point for the CV charging stage was when the charging current was less than 1/50 C, or the CV charging time had reached 1 h.Then, the battery was left standing in the opencircuit state.In the discharging procedure, the battery was fully discharged to the cutoff voltage (3.0 V) at 0.2 C.Then, the above test process was repeated to complete the available capacity tests at different C-rates.Throughout the testing, the discharge time, voltage, current, released capacity, and energy were monitored and recorded simultaneously.The current of the battery module in the charge and discharge processes at 2 C is shown in Figure 6.In the test, first the ternary battery module was fully charged by the constant current, constant voltage (CCCV) charging strategy.The constant current (CC) was 10 A, and the constant voltage (CV) was 4.2 V.The cutoff point for the CV charging stage was when the charging current was less than 1/50 C, or the CV charging time had reached 1 h.Then, the battery was left standing in the open-circuit state.In the discharging procedure, the battery was fully discharged to the cutoff voltage (3.0 V) at 0.2 C.Then, the above test process was repeated to complete the available capacity tests at different C-rates.Throughout the testing, the discharge time, voltage, current, released capacity, and energy were monitored and recorded simultaneously.The current of the battery module in the charge and discharge processes at 2 C is shown in Figure 6.In the test, first the ternary battery module was fully charged by the constant current, constant voltage (CCCV) charging strategy.The constant current (CC) was 10 A, and the constant voltage (CV) was 4.2 V.The cutoff point for the CV charging stage was when the charging current was less than 1/50 C, or the CV charging time had reached 1 h.Then, the battery was left standing in the opencircuit state.In the discharging procedure, the battery was fully discharged to the cutoff voltage (3.0 V) at 0.2 C.Then, the above test process was repeated to complete the available capacity tests at different C-rates.Throughout the testing, the discharge time, voltage, current, released capacity, and energy were monitored and recorded simultaneously.The current of the battery module in the charge and discharge processes at 2 C is shown in Figure 6.The available capacity of the tested LIB module at different C-rates is shown in Table 2.It can be seen that the battery released 96.12% of the maximum capacity at 0.2 C; while the battery released 94.52% of the maximum capacity at 1 C.When the current was further increased to 3 C, the battery only released 84.89%.The result was consistent with the "capacity nonlinear effect", and the relationship of the available capacity and current was nonlinear.The available capacity of the tested LIB module at different C-rates is shown in Table 2.It can be seen that the battery released 96.12% of the maximum capacity at 0.2 C; while the battery released 94.52% of the maximum capacity at 1 C.When the current was further increased to 3 C, the battery only released 84.89%.The result was consistent with the "capacity nonlinear effect", and the relationship of the available capacity and current was nonlinear.

Parameter Identification
The identification parameters of KiBaM, including capacity distribution ratio c, rate coefficient k, and fractional order α, were the key in achieving satisfactory accuracy based on the experimental data [31][32][33].The identification parameters were mainly arranged as follows.
Distribution ratio (c): The battery was charged until full, SOC = 1, as the initial state; and then the battery was discharged at a large constant current, here it was 95.69A (3 C).The released capacity was denoted as C 1 , and as shown in Table 2, C 1 = 27.59Ah.The maximum available capacity of the battery was 32.50 Ah, denoted as C max ; then, c = C 1 /C max = 0.849.
Rate coefficient (k): The available capacity of the battery was tested at different C-rates.For the classic KiBaM, we only needed one set of experimental data, and k could be obtained.The data measured at 1 C were selected here, and the data at other C-rates were used to verify the model's accuracy later.
A new function F(k ) can be defined and modified as follows: The change of the function F(k ) with k is depicted in Figure 7, and the solution, k = 0.000836, could be confirmed when F(k ) is 0. Thus, the unavailable capacity of the tested battery in KiBaM was obtained as follows: classic KiBaM, we only needed one set of experimental data, and k' could be obtained.The data measured at 1 C were selected here, and the data at other C-rates were used to verify the model's accuracy later.
-3468.6 31.88 1 3600 (1 0.849) .0.849 (32.5 30.72) A new function ( ) F k′ can be defined and modified as follows: -3468.6 31.88 1 ( ) (1 0.849) .0.849 1.78 3600 The change of the function F(k') with '  k is depicted in Figure 7, and the solution, ' k =0.000836, could be confirmed when F(k') is 0. Thus, the unavailable capacity of the tested battery in KiBaM was obtained as follows: As for the FO-KiBaM, k' and α of the tested battery module could also be identified by using one set of experimental data at 1 C.
The change of the function k and α is shown in Figure 8, where α changed from 0.1 to 0.9.We can see that it tended towards 0 when α was 0.9.Further, its trend was depicted in As for the FO-KiBaM, k and α of the tested battery module could also be identified by using one set of experimental data at 1 C.
And a new function F(α, k ) could be defined and modified as follows: The change of the function F(α, k ) with k and α is shown in Figure 8, where α changed from 0.1 to 0.9.We can see that it tended towards 0 when α was 0.9.Further, its trend was depicted in Figure 9 when α changed from 0.91 to 0.99.It was revealed that the numerical solution of F(α, k ) was not unique.The values of k and α could be a combination within a suitable range.Here, the parameters k = 0.000689 and α = 0.99 were selected for the FO-KiBaM.Thus, the unavailable capacity of the tested battery in FO-KiBaM was obtained.
was not unique.The values of ' k and α could be a combination within a suitable range.Here, the parameters ' k =0.000689 and 0.99 α = were selected for the FO-KiBaM.Thus, the unavailable capacity of the tested battery in FO-KiBaM was obtained.

Experiment Verification
The effect of the C-rate on the available capacity of the battery was verified by experiment.The variation of the total capacity and unavailable capacity of the battery at different C-rates are depicted in Figures 10-11, in which the discharge was finished when the total capacity of the battery (black line) intercepted the unavailable capacity (red/blue lines).At this point, all the remaining capacity of the battery was the unavailable capacity, so the discharge process was over.The prediction results of the available capacity of the battery are listed in Table 3.It can be seen that compared with the classic KiBaM, the errors of available capacity in the proposed FO-KiBaM were less, and it performed well over a wid applied current range.Specifically, its mean absolute error (MAE) was only 1.91%, with an improvement of 0.44%.Although the accuracy of KiBaM was already high, the proposed FO-KiBaM still had a smaller fitting error over a wide applied current range.Further study on the state estimations of LIBs and the high-precision optimization control of energy management will be useful, especially with the FO-KiBaM.
It should be pointed out that if the estimation accuracy was higher in some applications, we could achieve parameter identification according to current interval (segmentation) by using several sets of the tested data instead of just one set, so that the accuracy of the model could be improved.If we thought the parameters of the model were related to the current, we could even identify the parameters under different currents.Thus, that the accuracy of the model can be further improved, was not unique.The values of k and α could be a combination within a suitable range.Here, the parameters ' k =0.000689 and 0.99 α = were selected for the FO-KiBaM.Thus, the unavailable capacity of the tested battery in FO-KiBaM was obtained.

Experiment Verification
The effect of the C-rate on the available capacity of the battery was verified by experiment.The variation of the total capacity and unavailable capacity of the battery at different C-rates are depicted in Figures 10-11, in which the discharge was finished when the total capacity of the battery (black line) intercepted the unavailable capacity (red/blue lines).At this point, all the remaining capacity of the battery was the unavailable capacity, so the discharge process was over.The prediction results of the available capacity of the battery are listed in Table 3.It can be seen that compared with the classic KiBaM, the errors of available capacity in the proposed FO-KiBaM were less, and it performed well over a wid applied current range.Specifically, its mean absolute error (MAE) was only 1.91%, with an improvement of 0.44%.Although the accuracy of KiBaM was already high, the proposed FO-KiBaM still had a smaller fitting error over a wide applied current range.Further study on the state estimations of LIBs and the high-precision optimization control of energy management will be useful, especially with the FO-KiBaM.
It should be pointed out that if the estimation accuracy was higher in some applications, we could achieve parameter identification according to current interval (segmentation) by using several sets of the tested data instead of just one set, so that the accuracy of the model could be improved.If we thought the parameters of the model were related to the current, we could even identify the parameters under different currents.Thus, that the accuracy of the model can be further improved,

Experiment Verification
The effect of the C-rate on the available capacity of the battery was verified by experiment.The variation of the total capacity and unavailable capacity of the battery at different C-rates are depicted in Figures 10 and 11, in which the discharge was finished when the total capacity of the battery (black line) intercepted the unavailable capacity (red/blue lines).At this point, all the remaining capacity of the battery was the unavailable capacity, so the discharge process was over.The prediction results of the available capacity of the battery are listed in Table 3.It can be seen that compared with the classic KiBaM, the errors of available capacity in the proposed FO-KiBaM were less, and it performed well over a wid applied current range.Specifically, its mean absolute error (MAE) was only 1.91%, with an improvement of 0.44%.Although the accuracy of KiBaM was already high, the proposed FO-KiBaM still had a smaller fitting error over a wide applied current range.Further study on the state estimations of LIBs and the high-precision optimization control of energy management will be useful, especially with the FO-KiBaM.
Electronics 2019, 8, x FOR PEER REVIEW 13 of 17 but the amount of calculation would also increase.The whole process of parameter identification is the same, whether it is how to select the data or how many sets of data will be weighed.However, if the accuracy of the model identified by only one set of data met the application requirements, it would significantly reduce the amount of calculation.For a fair comparison and reduction in the number of calculations, we only selected one set of data to identify the parameters of both the KiBaM and the FO-KiBaM.

Conclusions
In recent years, various families of fractional-order systems have been found to be remarkably important and fruitful.Fractional calculus plays an important role in complex systems and, therefore, allows us to better describe real-world phenomena.By obtaining more parameters and degrees of freedom, fractional-order models can describe nonlinear performances of complex systems more accurately.Due to specific material and chemical properties of batteries, fractional calculus is more reasonable to describe the nonlinear performance of a battery's capacity.The proposed FO-KiBaM can describe the battery's nonlinear characteristics more accurately, with greater flexibility and novelty compared to the classic KiBaM.The proposed model can be applied in engineering.The estimation of the available capacity of LIBs is meaningful with a wide applicable current range, and  It should be pointed out that if the estimation accuracy was higher in some applications, we could achieve parameter identification according to current interval (segmentation) by using several sets of the tested data instead of just one set, so that the accuracy of the model could be improved.If we thought the parameters of the model were related to the current, we could even identify the parameters under different currents.Thus, that the accuracy of the model can be further improved, but the amount of calculation would also increase.The whole process of parameter identification is the same, whether it is how to select the data or how many sets of data will be weighed.However, if the accuracy of the model identified by only one set of data met the application requirements, it would significantly reduce the amount of calculation.For a fair comparison and reduction in the number of calculations, we only selected one set of data to identify the parameters of both the KiBaM and the FO-KiBaM.

Conclusions
In recent years, various families of fractional-order systems have been found to be remarkably important and fruitful.Fractional calculus plays an important role in complex systems and, therefore, allows us to better describe real-world phenomena.By obtaining more parameters and degrees of freedom, fractional-order models can describe nonlinear performances of complex systems more accurately.Due to specific material and chemical properties of batteries, fractional calculus is more reasonable to describe the nonlinear performance of a battery's capacity.The proposed FO-KiBaM can describe the battery's nonlinear characteristics more accurately, with greater flexibility and novelty compared to the classic KiBaM.The proposed model can be applied in engineering.The estimation of the available capacity of LIBs is meaningful with a wide applicable current range, and only a set of data at one C-rate is needed to accurately estimate the available capacity at different rates, which greatly reduces the number of calculations.What is more, the proposed FO-KiBaM provides a basic battery model for further research on SOC estimation methods.

Patents
A patent, termed fractional order KiBaM (kinetic battery model) that considers nonlinear capacity characteristics and parameter identification methods, resulted from the work reported in this manuscript, which can be seen on academic websites such as Google Patents.The progress in the battery b.The remaining available volume of water changes linearly with flow rate.c.The available volume of water has nothing to do with the temperature within a certain range.d.No recovery effect when stop releasing water.a.The available volume of water is a constant, nothing to do with flow rate.b.The remaining available capacity of battery changes nonlinearly with C-rate.c.The available capacity of battery is closely related to the temperature.d.Recovery effect when stop discharge battery.a.The available capacity of battery is related to discharge C-rate, the larger the current, the less the available capacity .

Figure 1 .
Figure 1.Battery capacity is not similar to water in a bucket.

Figure 1 .
Figure 1.Battery capacity is not similar to water in a bucket.

Figure 3 .
Figure 3. Battery capacity changes with the discharge time in the discharge process.Figure 3. Battery capacity changes with the discharge time in the discharge process.

Figure 3 .
Figure 3. Battery capacity changes with the discharge time in the discharge process.Figure 3. Battery capacity changes with the discharge time in the discharge process.

Figure 4 .
Figure 4.The battery test experimental platform.

Figure 5 .Table 1 .
Figure 5. Battery module with 13 cells connected in parallel.Table 1. Specification of the ternary lithium-ion battery (LIB) cell.Parameter Value Parameter Value Type LR1865SZ Series-parallel 1S-13P Nominal voltage 3.6 V Rated capacity 2.5 Ah End-of-charge voltage 4.2 V End-of-discharge voltage 3.0 V

Figure 5 .
Figure 5. Battery module with 13 cells connected in parallel.

Figure 7 .
Figure 7.The function curve that changes as k' changes.

Figure 7 .
Figure 7.The function curve that changes as k changes.

Figure 8 .
Figure 8. Trend of the function as k' and α (0.1 to 0.9) changes.

Figure 9 .
Figure 9. Trend of the function as k' and α (0.91 to 0.99) changes.

Figure 8 .
Figure 8. Trend of the function as k and α (0.1 to 0.9) changes.

Figure 8 .
Figure 8. Trend of the function as k' and α (0.1 to 0.9) changes.

Figure 9 .
Figure 9. Trend of the function as k' and α (0.91 to 0.99) changes.

Figure 9 .
Figure 9. Trend of the function as k and α (0.91 to 0.99) changes.

Figure 10 .
Figure 10.Battery capacity variance at different C-rates.

Table 2 .
Available capacity and discharge time at different C-rates.

Table 2 .
Available capacity and discharge time at different C-rates.

Table 3 .
Prediction results of the battery's available capacity at different C-rates.

Table 3 .
Prediction results of the battery's available capacity at different C-rates.