Constructing accurate equivalent electrical circuit models of lithium iron phosphate and lead-acid battery cells for solar home system applications

The past few years have seen strong growth of solar-based off-grid energy solutions such as Solar Home Systems (SHS) as a means to ameliorate the grave problem of energy poverty. Battery storage is an essential component of SHS. An accurate battery model can play a vital role in SHS design. Knowing the dynamic behaviour of the battery is important for the battery sizing and estimating the battery behaviour for the chosen application at the system design stage. In this paper, an accurate cell level dynamic battery model based on the electrical equivalent circuit is constructed for two battery technologies: the valve regulated lead–acid (VRLA) battery and the LiFePO 4 (LFP) battery. Series of experiments were performed to obtain the relevant model parameters. This model is built for low C-rate applications (lower than 0.5 C-rate) as expected in SHS. The model considers the non-linear relation between the state of charge ( S O C ) and open circuit voltage ( V OC ) for both technologies. Additionally, the equivalent electrical circuit model for the VRLA battery was improved by including a 2nd order RC pair. The simulated model differs from the experimentally obtained result by less than 2%. This cell level battery model can be potentially scaled to battery pack level with flexible capacity, making the dynamic battery model a useful tool in SHS design.


Introduction
As of 2016, almost 1.1 billion people around the world lacked access to electricity [1].While gridextension is not always the most plausible solution, solar-based off-grid solutions such as solar home systems (SHS) have become prominent in use.Battery storage is an essential component of an SHS, as it is not only the most expensive component, but also the component with usually the least lifetime [2].Therefore, understanding battery behaviour is critical for designing and analysing SHS.
In this study, two battery technologies commonly used in SHS were explored: the valve regulated lead-acid (VRLA) battery and the LiFePO 4 (LFP) battery.Lead-acid is the cheapest technology in the market in comparison with the other common battery technologies [3].In particular, the valve regulated type of lead-acid battery is sealed, which reduces the risk of liquid leaking when compared to flooded lead-acid battery.Moreover, VRLA battery is entirely maintenance free [4].In the case of LFP battery, it has surged as a predominant technology due to its appropriate chemical and thermal stability together with its high energy capacity [5].Additionally, LFP batteries have experienced a reduction in price that is expected to continue [6].

Selecting the Battery Model
Depending on the application, accuracy, complexity, compatibility and universality, there are several technical methods to approach an expected battery model.Generally, there are three levels of battery models according to the degree of physical insight, viz., white box (e.g., electrochemical/physical model), grey box (e.g., equivalent electrical circuit model) and black box (e.g., artificial neural network model) [7].
The equivalent electrical circuit model (EECM) presents an equilibrium among computational time, precision, and complexity.As the EECM models the battery with only electrical components, this kind of models is highly compatible with other electrical models, which make them easy to implement with battery management systems and in other electrical engineering applications.Moreover, the construction of EECM does not have a very comprehensive co-relation with electrochemical processes inside the battery.As a result, there is no limitation on distinct battery types when modelling the battery by EECMs.Hence, the electrical equivalent circuit model (EECM) is selected in study.

Importance of Low Current Battery Models for SHS
Battery modelling is essential in SHS design, as they help to elucidate the relationship between power required by a predefined load and the operational conditions in which the batteries are used.Particularly, in SHS, the C-rates (related to the rate at which a battery is charged or discharged) are usually lower than 0.5 C, with the lowest C-rate being 1/120 C to 1/100 C [8,9].However, the existing battery cell level models do not focus on low C-rates; instead, they normally study a wide range of C-rates with high sampling resolution.For example, C-rates of 0.1 C, 0.5 C, 1 C, 2 C, 3 C, 4 C, 7 C, and 9 C have been proposed in [10], while other authors have explored C-rates of 0.04 C, 0.1 C, 1 C, 2 C, 3 C, 5 C, and 10 C [11], and 0.1 C, 0.5 C, 1 C, 2 C, 3 C, and 4 C [12].Therefore, the applicability of these models in wide C-rate ranges does not necessarily validate their suitability for SHS.
Consequently, we focus on the low battery C-rates reflecting typical SHS applications.To obtain an accurate battery model which suitable for low C-rate applications, experiments with narrower C-rate range but with high sample resolutions are proposed.The applicable C-rate range of this model is 0.05 C-0.5 C with the sample resolution of 0.05 C-0.1 C.

Contributions
This papers contributes to the current state of the art by

Battery Parameters
State of the Charge (SOC): SOC is usually defined as the proportion of the maximum possible electric charge that is present inside a rechargeable battery [13].It is calculated by the ratio between the difference of the battery capacity and the net electric charge discharged from a battery since the last full state of charge.The SOC can be expressed by Equation (1).
where Q t is the temperature dependent actual capacity, Q e is the capacity that be extracted from the battery starting from its full state, and Q r is the capacity that remaining in the battery.I m is the main branch current, which is the result of the load current applied to the battery minus the current consumed in side-reactions (I p ).I m is positive if the battery is discharging and I m is negative if the battery is charging.
Electromotive force (emf): The electromotive force of a cell is the algebraic sum of the potential difference of two half cells, which is the chemical potential of redox reactions happening on these two electrodes [14].
Open circuit voltage (OCV, or V OC ) and cte-OCV: The open circuit voltage (OCV) is the potential difference across its terminals when there is no current flow in or out of a reversible cell [15].If the battery rests enough time (depending on the battery technology) with no current flow, the OCV can be called Close-to-Equilibrium open circuit voltage (cte-OCV) [16].This value can be regarded as the estimated emf value.
Current rate (C-rate): The current rate, which is also called C-rate, is defined as the magnitude of charge/discharge current in Ampere to the nominal capacity of the battery in Ampere-hour.The C-rate represents the speed of charge/discharge, and it can be calculated using Equation (4) [17].
where I batt is the current applied on the battery in the unit of ampere (A), and Q n is the nominal battery capacity in the unit of ampere hour (Ah).

Construction of the Dynamic Battery Model
A good battery model should be capable of predicting both the battery storage capability information and the voltage response to the load.An equivalent electrical circuit model (EECM) can describe both characteristics.In some cases, modelling the side reactions is required in terms of battery losses, which is also possible to be realized with the EECM.
As shown in Figure 1, there are three main parts in the EECM to model the different features of the battery.First part is the storage circuit.Three parameters should be modelled in this sub-circuit: the battery capacity (which is the charge storage ability), the amount of charge is left in the battery (SOC), and the open circuit voltage V OC .
The second part is the voltage response circuit, which focuses on the internal structure of the battery.This makes it possible to imitate the electrical response of the battery under different loads.In this paper, a Thevenin model is applied, and it consists of a resistance in series with RC pairs.The last part is the parasitic branch, which is used to model the side reactions that induce loss of charge in the battery.There are many ways to model these battery losses.For example, Figure 1 shows the parasitic branch in the form of a resistance (R p ).For both battery technologies, the basic construction of EECM is applicable, as shown in Figure 1.

Storage Circuit
As discussed in Section 2.1, the cte-OCV can be regarded as the emf value.As the emf values are SOC dependent, the cte-OCV can be used as a SOC indicator.In this study, the SOC value is selected as the independent variable, while the coulomb counting method is used to determine the SOC.The SOC can be expressed by Equation (1).

Electrical Response Circuit
There is a deviation of voltage from the equilibrium state, also called polarization, whenever the battery is loaded.Conversely, the voltage recovers back towards the equilibrium state whenever the load is disconnected.Hence, a voltage relaxation phenomenon is observable in a no-load interval (also called rest interval) when the battery is operating.The battery internal impedance can be potentially extracted from this voltage relaxation.
The non-load rest interval during a discharge stage is shown in Figure 2b.The voltage response can be divided into two parts, viz., an instantaneous voltage jump right after the load is paused, followed by a gradual increase of the voltage towards the cte-OCV.Those two voltage responses are correlated with two electrical parts: the ohmic resistance R o , and the resistance/capacitance pairs (R 1 C 1 , R 2 C 2 ), as introduced in the construction of EECM model in Figure 1.
Generally, there is an expected time window for each of the voltage response.An EECM containing two RC-pairs is assumed.It is suggested that, within a five seconds time-window, the voltage boost (V o ) right after the load disconnection is induced by the ohmic resistance R o [18].The time window for the rest of the gradual voltage increase caused by the RC pairs is called the time constant.The time constant reflects the rate of change of voltage; if the voltage changes faster, the time constant is smaller [19].The increase in voltage right after the immediate voltage increase is the "short-term voltage response".This short time interval (τ 1 ) is expected to be around 10 seconds, and this voltage increase (V 1 ) is related to the first RC-pair, which is the R 1 , C 1 in Figure 2b [20].The following increase in voltage (V 2 ) is then treated as the long-term voltage response after a long time interval (τ 2 ), which corresponds to the second RC-pair (R 2 , C 2 ), as shown in Figure 2b.It is possible to increase the number of RC-pairs and divide the voltage relaxation into more parts, with which the accuracy of the model is improved, such as the R n , C n in Figure 2b.Many authors suggest to model the lead-acid battery with first order RC circuit.However, for Li-ion batteries, all three orders of RC circuit model are commonly seen [12,[21][22][23][24].Eventually, the relaxation voltage will reach cte-OCV if the rest interval is long enough.However, it is necessary to strive an optimum balance between accuracy and complexity.Therefore, the study described in this paper has considered two RC pairs.Rest interval Curve fitting method was applied on each of the voltage relaxation results.The recognition of each part of the voltage and the corresponding time constant was done by nonlinear least squares curve fit in Matlab [20].The distinct voltage and time constant values were then fitted into the RC circuit complete response equations to obtain the values of electrical elements in the model.The equations are listed below as Equations ( 5)-( 9).

Parasitic Reaction Circuit
The electrical components interpreted in the previous subsections describe the behaviour of the battery based on its main reactions.During operation, not all the battery energy is involved in its main (charging/discharging) reactions; there are unwanted side reactions responsible for loss of useful energy.Therefore, additional elements are needed to model a parasitic branch considering the side reactions.
Especially when charging a lead-acid battery, significant losses can occur, which has a huge impact on the battery dynamic behaviour.This is particularly true when the battery is at a high SOC (>95%).Therefore, the charge loss of VRLA during the charging stage is considered as the most dominant loss mechanism.This loss is mainly due to gassing and oxygen recombination, but the transient gassing phenomenon can be quite complex and unpredictable especially for VRLA batteries [25].
There are different methods to model the parasitic branch.Some models use a self-discharge conductance to represent the side reactions, where the conductance is the function of battery voltage, or side reaction related voltage and the temperature [21,22].Some other models represent the side reaction by means of gassing current, as the function of water decomposition voltage and the temperature, instead of modelling the side branch circuit [26].
In this study, an equation is derived from the one proposed in [27].This method is based on the principle that the losses incurred in the parasitic branch can be viewed as the Coulombic inefficiency of the battery.Therefore, this method models the parasitic branch in the form of the Coulombic efficiency as a function of SOC and the battery current, as shown in Equation (10).
where I ref is the reference current, and I test is the test current.a, b, and c are the experimentally derived empirical constants.It must be noted that, as the current terms in Equation ( 10) occur in the form of a ratio, the Coulombic efficiency can be considered as a function of C-rate as well.

Choice of Operational Variables
For both battery technologies, each of the internal impedances is known to be, to some extent, influenced by the SOC, the temperature, and the current [28].To obtain a proper dynamic battery model, a battery should ideally be tested with all three variables.However, the influence of operating temperature was not investigated in this study.All experiments were performed in a laboratory with the temperature controlled in the range of 20-22 • C. Therefore, all experiments can be considered as under constant room temperature conditions.SOC and current (in the form of C-rate) were the main operational variables considered in this study.

Overall Methodology
The EECM model is assumed as shown in Figure 1, with two RC pairs.The choice of two RC pairs is justified in Section 4. For the task of fully constructing the EECM model for both VRLA and LFP batteries, the following methodology was followed involving experiments, data analysis, and model verification.

1.
Experimentation.Series of experiments were performed to extract the parameters of each electrical element pertaining to every sub-circuit.These measurements were applied on each featured SOC and with selected current rate, as further described in Table 1.

2.
Parameter extraction.The parameters of each component in the EEMC were extracted and analyzed based on the experiments.The values of those parameters were then summarized into equations, and those equations can represent each of the electrical element in the electrical circuit.

3.
EECM construction and model verification.All parts of the electrical elements were put together to construct the full EECM.Finally, the constructed EECM of each battery was simulated and compared with the experimental results.

Equipment and Materials
In this study, the battery tester used is MACCOR R Series 4000 Automated Test System [29].The batteries selected for tests are Cyclon R AGM D single cell and A123systems R ANR26650M1B, and both of the battery cells have 2.5 Ah capacity [30,31].
The basic charge and discharge methods applied in this study are constant current (CC) method and constant current-constant voltage (CC-CV) method.The CC method involves charging or discharging the battery with a constant current until the battery voltage reaches a specific value.The specified voltage to stop the discharging process is called the end of discharge voltage (EODV).The CC-CV method involves first charging/discharging the battery with constant current until the battery voltage reaches the pre-set value, and then switch to constant voltage stage, which involves charging/discharging the battery with a constant voltage.
The in-detailed charge and discharge method of each battery technology is summarized in Table 1.This follows the recommendations and best practices outlined in the data sheets and previous articles [20,[30][31][32].

Storage Circuit
The purpose of this set of experiments is to obtain the relation between the cte-OCV and the SOC.Therefore, the cte-OCV should be measured at different SOC values, and the battery is fully discharged before the measurement.During the experiments, the battery was monotonically charged to full and then discharged to empty, with both procedures being executed in steps of 5% SOC.Coulomb-counting method was used to approach the designated SOC demarcation.After each step, the battery was rested for a constant time to reach the close-to-equilibrium state, at which point the voltage was measured and noted as the cte-OCV value.The C-rate was 0.2, while the rest time was 3 h for VRLA battery and 4 h for LFP battery.Three battery cells were tested to get a result with better accuracy.
The experimental procedure is illustrated in the flowchart shown in Figure 3.

Voltage Response Circuit
The electrical elements in the model were tested while the battery underwent eight different C-rates at different SOC levels.The test procedure esd similar to the cte-OCV test, as shown in Figure 3, with the only difference being the amount of energy being charged/discharged in each step, and the resting time interval.In this test, the step size of charge/discharge was 0.15 Ah, and the relaxation time between each step was 3 min for both battery technologies.The resting time between charge/discharge was 1 h for both battery technologies.

Parasitic Branch
A simplified measurement was applied to estimate the Coulombic efficiency in this study.The Coulombic efficiency was approximately measured and calculated from the integral recharge efficiency.The test procedure is outlined as follows.
1. Fully charge and discharge the battery to measure the initial battery capacity and the overall Coulombic efficiency.2. Discharge the battery to a specific SOC value (x%) and note down the discharged capacity Q d .3. Fully recharge the battery and record the recharged capacity Q c . 4. Calculate the recharge efficiency Q d /Q c assuming it is the averaged integral recharge efficiency at the mean SOC between x% and 100%.5. Go back to Step 2 and repeat the steps until enough data points are obtained.6. Test the whole procedure under different C-rates.
One thing should be noted: after each charge/discharge step, the battery needed to rest until reaching its equilibrium state, for 1 h in this case.
Then, the averaged integral recharge efficiency was fitted to Equation (10) to get the roughly estimated Coulombic efficiency as a function of SOC and C-rates.
The experimental settings are summarised in Table 2.

EECM Construction
After achieving all the required values of the electrical elements, the battery model can be constructed according to the EECM as shown in Figure 4.It must be noted that the parasitic resistance R p shown in Figure 4 is not an actual circuit element in the EECM but only a notional lossy element that signifies side-reaction.In the proposed model in this study, the VRLA losses are actually quantified through the Coulombic efficiency, as shown in Equation (10).
The SOC of the battery can be represented by Equation (1), while the main branch current I m is calculated through Equation (11).The Coulombic efficiency η c is obtained from Equation (10).
The dynamic battery voltage of each part of the circuit can be calculated by using the achieved electrical elements through simple circuit analysis, as shown in Section 2.4.The battery voltage is shown in Equation (12).
where V batt is the battery voltage and the ± sign in Equation ( 12) is dependent on the current direction, i.e., + implies charging and − implies discharging the battery.

Cte-OCV as a Function of SOC
The cte-OCV and SOC are correlated; as the SOC increases, the cte-OCV also rises.Although many papers and manufacturer datasheets suggest that the cte-OCV of VRLA versus SOC shows a linear relation [33,34], our study found that a rational curve fit is better than a linear fit.This new fit is especially relevant when the battery is at low SOC, as plotted in Figure 5. Hence, the empirical-based rational Equation ( 13) was proposed and employed instead of a linear relationship.In the case of LFP battery (Figure 6), there is a sudden increase in V OC when the SOC is close to 100%; similarly, the V OC drops when the SOC is almost 0%.The cte-OCV and SOC behaviour can be perfectly fitted into a double exponential curve [35].This has been verified by this study, as shown in Figure 6, where the hysteresis phenomenon between the V OC values when charging or discharging is depicted.This hysteresis phenomenon was modelled in the electrical response circuit instead of V OC -SOC equation for both batteries.
The fitted cte-OCV-SOC result of LFP battery is written in Equation ( 14).

Internal Impedance
To choose an optimised EECM structure, the fitting comparison from first to third orders of RC circuits are plotted in Figure 7a From these curves, it can be concluded that the second order RC circuits have sufficient accuracy for both battery technologies.In addition, modelling the battery with a higher order of RC circuit requires more parameters.Thus, the second order RC circuits were selected for modelling both battery technologies owing to their lower complexity but sufficient accuracy.
The measured and fitted results of electrical elements in EECMs of LFP and VRLA batteries are shown in Figures 8-11, where R o is the ohmic resistance; R s , C s are the values of RC pair representing the voltage relaxation in the short time-window (τ 1 in Figure 2b); R l , C l are the values of RC pair used to represent the voltage relaxation in the long time-window (τ 2 in Figure 2b).In Figures 8 and 9, some features of VRLA batteries can be inferred with respect to the internal impedance circuit elements, as listed below.
1.All the electrical elements show a strong dependency on SOC.
2. All the elements also show an inverse trend versus SOC in charge stage and discharge stage.3.In terms of current dependency, elements during charge and discharge stages have diverse behaviours.While charging, only R l and C l values show a clear operational current dependence.
While discharging, all the elements except R o are evidently influenced by the operational current.
Therefore, the R l and C l in discharge stage and all the elements in charge stage were modelled as the function of operational current.The detailed equations are listed as below: Charge R s_VAc = 4.246e − 13 • e 29.05•SOC + 0.03787 • e 1.467•SOC ( 16) Discharge  From the results in Figure 10 and 11, some conclusions for LFP batteries can be drawn: 1. Except for R o , all other electrical elements show a strong dependency on SOC. 2. All electrical elements exhibit an inverse behaviour versus SOC in charge and discharge.
3. The operational current influences R s , R l and C l in both charge and discharge processes.
Therefore, all the elements could be expressed as the function of SOC for better accuracy, even though Ro does not vary too much with different SOC values.Elements R s , R l and C l were modelled as the function of both SOC and C-rate.
In these two figures, it can be observed that the fitting of all the electrical elements for battery model construction has a high accuracy.The equations of all these electrical elements for LFP battery are listed below. Charge R s_LPc = 0.01036 • e 0. The Coulombic efficiency is expected to decrease when current increases.However, the efficiency measured with 0.2C is higher than that measured with 0.1 C, and the result measured with 0.4 C is larger than the one measured with 0.3 C.This is because 0.1 C and 0.3 C measurements were operated on the same battery, while a different battery was tested with 0.2 C and 0.4 C. The differences between batteries may lead to this consequence.Nevertheless, the Coulombic efficiency tested with higher C-rate is lower than one tested with lower C-rate for the same battery.Therefore, the Coulombic efficiency is fitted into the modified Equation (10), which was introduced in Section 3.6.The result is shown in Equation (35).
where C r is the C-rate.

Simulation and Validation
In this section, the simulation results and the experimental data are compared.The comparison of measured voltage and simulated voltage for both battery technologies are plotted in Figures 13 and 14  As stated in Section 2.2, there are three essential features that should be reflected in the model: the storage feature, the voltage response and the parasitic side reaction.In this study, the SOC is the independent variable, and the voltage response is the main feature being modelled.The Coulombic efficiency is included in the storage circuit as well as the voltage response circuit.Hence, the voltage is selected to evaluate the accuracy of the model.To quantify the model accuracy, all the modelled and measured voltage values are calculated and evaluated using Equation (36).Then, for each battery, there is an averaged error value achieved among the whole charge/discharge process.
where V sim is the simulated voltage and V exp is the measured voltage.
The constant voltage (CV) charging stages were omitted from the calculation because, in the real measurement, the pause steps are accomplished by the tester, and it depends on the natural battery features.This step counting outcome accomplished by the tester in CV stage is unable to be repeated with the simulation.How it works is reflected in Figures 13 and 14.For example, the comparison results of VRLA batteries is shown in Figure 13.The discrepancy in this figure between the experimental and the simulation results (e.g., 600-1000, 250-700 and 150-650 min in Figure 13a-c, respectively) are due to the step counting method difference.Therefore, the CV stages for both battery technologies are eliminated.
The comparison results of VRLA batteries are in Figure 13.The simulated voltage curve increases slower than the measured one during the last charging stage before the battery reaches the CV stage, which is, e.g., around 600 min in Figure 13a.This slower increase trend of voltage is due to the modelled Coulombic efficiency loss.This can be attributed to the fact that the number of measured Coulombic efficiency data points when the SOC > 98% are insufficient.The fitted Coulombic efficiency result then cannot model the VRLA battery behaviour accurately when it reaches the high SOC stage, especially when the SOC higher than 98%.
For VRLA batteries, the maximum averaged error among all eight different C-rates is 1.7664%, while discounting the CV charging stage.The Coulombic efficiency fitting results might be one of the reason causing the error here.
Figure 14 shows the comparison results of LFP batteries.For LFP batteries, the maximum averaged error among all different C-rates, which is eight batteries in total, is 1.6708%.The discrepancy of the wedges shaped charge/discharge steps is mainly due to the differences in ohmic resistance.The error from ohmic resistance counts at least 0.5%.
The dynamic models of both battery technologies are therefore considered accurate.

Conclusions
In this study, two battery technologies (VRLA and LFP) were chosen for experimentally determining the dynamic behaviour of the battery, with the eventual goal of constructing an EECM model.The constructed model works well throughout the C-rate range from 0.05 C to 0.5 C, with an error lower than 2% for both LFP and VRLA batteries.The chosen C-rate range gave a better granularity for the intended application of SHS.The non-linear relation between SOC and V OC , the use of 2nd order RC circuit in EECM, and the quantification of the parasitic branch of VRLA cells in the form of SOC and C-rate based Coulombic efficiency η c were identified as methodological contributions for modelling the dynamic behaviour of VRLA batteries.Additionally, the model developed in this work is at the cell level, which is fully scalable to battery pack level.The proposed model can be used by SHS designers for understanding application-specific battery behaviour in the context of system design and analysis.

Recommendations and Future Work
Firstly, it is recommended to identify and develop a more accurate methodology to measure the Coulombic efficiency.The methodology presented in this paper is sufficient over most of the SOC range, except beyond 98%, where the achievable SOC granularity of the method is lower than desirable.Even though more accurate industrial test methods for Coulombic efficiency measurement are mostly chemical-based, it is still worthwhile developing a simple method with sufficient accuracy over the entire operational SOC range, especially for electrical models.Secondly, the batteries' dynamic behaviour was modelled only when the batteries were new.The dynamic performance evolution through ageing is currently ignored.This is proposed to be improved in the future work.Thirdly, the impact of temperature needs to studied as an independent and controllable variable.Finally, this EECM is proposed to be applied in actual SHS use-cases by using battery charging and discharging profiles from SHS applications.

Figure 1 .
Figure 1.Schematic of the overall EECM construction.

Figure 2 .
Figure 2.An example of voltage relaxation during rest interval in discharge with corresponding circuit elements.(a) The electrical circuit elements; (b) The voltage relaxation curve.
Storage circuit.The cte-OCV measurement method was performed.(b) Voltage response circuit.Voltage relaxation (also called step response) method was used.(c) Parasitic branch.The differential recharge efficiency measurement method was used.

Figure 3 .
Figure 3. Test procedure of cte-OCV-SOC measurement and the Internal impedance circuit measurement.

Figure 4 .
Figure 4.The construction of the EECM.The subscripts s and l stand for short and long time duration respectively with respect to voltage relaxation (τ 1 , τ 2 in Figure 2b).

Figure 5 .
Figure 5. Measured open circuit voltage versus SOC of VRLA battery: (a) V OC measured during charge and discharge separately; and (b) the measured data points and the fitted curves.

Figure 6 .
Figure 6.Measured open circuit voltage versus SOC of LFP battery: (a) the hysteresis of V OC when charge and discharge the battery separately; (b) the measured data points and the fitted curve.

Figure 7 .
Figure 7.Comparison of accuracy with modelling battery into different orders of RC circuits: (a) VRLA battery; (b) LFP battery.

Figure 8 .Figure 9 .
Figure 8.The values of RC elements and their fitted results: VRLA batteries during charge.

Figure 10 .Figure 11 .
Figure 10.The values of RC elements and their fitted results: LFP batteries during charge. .

Figure 13 .Figure 14 .
Figure 13.Comparison of experimental result and simulation result-VRLA battery: (a) test with 0.1 C; (b) test with 0.25 C; and (c) test with 0.5 C.

Table 1 .
Specification and summarized charge/discharge method for both battery technologies.

Table 2 .
Experimental setting for Coulombic efficiency measurement for VRLA battery.