Lithium-Ion Battery State-of-Charge Estimation Using Electrochemical Model with Sensitive Parameters Adjustment

: State-of-charge (SOC) estimation of lithium-ion (Li-ion) batteries with good accuracy is of critical importance for battery management systems. For the model-based methods, the electrochemical model has been widely used due to its accuracy and ability to describe the internal behaviors of the battery. However, the uncertainty of parameters and the lack of correction from voltage also induce errors during long-time calculation. This paper proposes a particle ﬁlter (PF) based method to estimate Li-ion batteries’ SOC using electrochemical model, with sensitive parameter identiﬁcation achieved using the particle swarm optimization (PSO) algorithm. First, a single particle model with electrolyte dynamics (SPME) is used in this work to reduce the computational burden of the battery electrochemical model, whose sensitive parameters are selected through the elementary effect test. Then, the representative sensitive parameters, which are difﬁcult to measure directly, are adjusted by PSO for a high efﬁciency. Finally, a model-based SOC estimation framework is constructed with PF to achieve accurate Li-ion battery SOC. Compared with extended Kalman ﬁlter and equivalent circuit model, the proposed method shows high accuracy under three different driving cycles.


Introduction
With the energy crisis and the deterioration of the ecological environment, electric vehicles (EVs) and smart grids have received widespread attention around the world.Among the existing energy storage devices, lithium-ion (Li-ion) batteries have been widely utilized for their distinctive advantages, such as high energy density, long service life, low self-discharging rate, and no memory effect [1][2][3].In order to ensure the safety and reliability and achieve the best working performance of the energy storage system with a proper management strategy, accurate state monitoring is needed for the battery management system (BMS).Battery state is a complicated concept that usually includes multitudes of branches, such as state of charge (SOC), state of health (SOH), state of power (SOP), etc.One of the key variables that must be estimated in the BMS is the SOC, which is the ratio of the battery's current remaining capacity to its maximum value [4][5][6].
However, it is inaccessible to accomplish the direct SOC measurement by BMS.Usually, some relevant measurable variables, such as current, voltage, and temperature, are needed for battery SOC estimation.The commonly used SOC estimation methods include the ampere-hour integration method [7,8], the open circuit voltage (OCV) method [9,10], and the model-based method.The ampere-hour integration method has obvious drawbacks of the measurement error accumulation and the uncertainty of the initial SOC, whereas it is the simplest method.Besides, the ampere-hour integration method only uses the current to estimate SOC without the correction and information from terminal voltage.
The OCV method can recompense the abovementioned shortcoming utilizing the determinate relationship between SOC and OCV.Nevertheless, the SOC estimation cannot be achieved online by the OCV method because battery OCV measurement requires a long relaxation period to reach its internal equilibrium.
The model-based methods commonly deploy the equivalent circuit model (ECM) or the electrochemical model (EM) to reach the goal of SOC estimation.Based on battery model, filtering algorithm, such as extended Kalman filter (EKF) [11], unscented Kalman filter (UKF) [12], and particle filter (PF) [13], can be applied to achieve SOC estimation.Meanwhile, a battery model can be combined with an observer, such as Luenberger observer (LO) [14], sliding-mode observer (SMO) [15], or nonlinear observer (NO) [16], to reach the same goal.ECM [17,18], including the Rint model, Thevenin model [19,20], firstorder RC model [21], partnership for new generation of vehicles (PNGV) model, etc., use resistors, inductors, and capacitors to simulate the external characteristics of a Li-ion battery, and have been widely applied to estimate SOC in practice thanks to their simple structure, high calculation efficiency, and few parameters to be identified.However, there remain some drawbacks for ECM, e.g., parameters of the ECM are not directly related to the electrochemical reactions inside the Li-ion battery.In this way, the parameters in ECM lose their specific physical meaning in reality, which further contributes to the incapability to reflect the internal states, such as the overpotential and the lithium-ions concentration, and also leading to limitations in terms of modeling accuracy.The SOC estimation errors can be divided into dynamic transient error and steady state error.The dynamic transient error can be eliminated if the gain of observer is selected to be a big value.However, the steady state error is independent of the observer gain and cannot be neglected.So, one consensus for the model-based SOC estimation method is that the estimation algorithm cannot well handle the battery modeling errors, which restricts the usage of ECMs for advanced BMS applications [22,23].Meanwhile, the EM can provide more internal details on the battery, which could support an accurate battery model, as well as the charging strategy and power prediction of the Li-ion battery.Thus, researchers focus on using EM to obtain a superior model-based SOC estimation method.
The original EM is a pseudo-two-dimensional porous electrode (P2D) model proposed by Doyle, Fuller, and Newman [24,25].The P2D model describes the battery internal reactions according to five subsections: porous positive electrode, porous negative electrode, separator, and the current collector in both positive and negative area.The P2D model is composed of a series of coupled partial differential equations (PDEs) that reflect the internal electrochemical reactions of the Li-ion battery, including the diffusion equation in solid and electrolyte phase, ohm equation in solid and electrolyte phase, charge conservation equation, and the Butler-Volmer equation.Although it can improve the modelling accuracy and widens the application range of the battery model to some extent, the computational burden of EM is aggravated owing to the difficulty of solving the coupled PDEs [26].An immense amount of research has emerged to reduce the complexity of EM.Mathematical approaches, such as finite element method (FEM) [27], method of linear (MOL) [28], finite difference method (FDM) [29], finite volume method (FVM) [30], and Padé approximation [31], have been widely studied.Meanwhile, B. Haran et al. [30] propose the concept of the single particle model (SPM) [32] using single spherical particles to represent the positive and negative electrodes.It is assumed the current is uniformly distributed in the electrode and the electrolyte dynamics are ignored.In this way, the mathematical calculation of the EM is simplified, and the SPM is widely used for the online estimation of SOC [33][34][35].Unfortunately, the simplification of electrolytes may decrease the calculation burden, and at the same time cause some unwanted errors which limit the application of the SPM in low C-rates circumstances.It is noted that SPM will lose fidelity once the C-rates is above C/2.Therefore, one challenge here is to properly choose a suitable EM for Li-ion battery SOC estimation.
Since the EM has a good interpretability and can precisely describe the state variations inside the battery, it has become a trend in the future development of the model-based SOC estimation method.The computational accuracy of EM depends remarkably on the accuracy of the parameters, so how to obtain parameters for such an intricate model is a crucial yet challenging task [36].Some parameters can be acquired by directly disassembling the commercial battery through laboratorial tests, such as geometric parameters and SOC-OCV curves in the electrodes.However, those microscopic parameters, such as diffusion coefficients and kinetic reaction rates, cannot be measured experimentally and need to be identified by specific algorithms.Li-ion battery parameter identification is generally performed using the nonlinear least squares (NLS) method [37,38], genetic algorithm (GA) [39,40], and particle swarm optimization (PSO) method [41,42].Due to the simplicity of the components and the fewness of parameters in ECM, the NLS method can well complete the identifying task.However, the parameters required to reflect the internal electrochemical reaction mechanisms are extensively increasing in EM.Meanwhile, EM has become highly nonlinear and due to a lack of analytical solutions, more advanced algorithms for parameter identification are in an urgent demand.Hence, the second challenge using EM for SOC estimation is the demand for a parameter adjustment strategy, which should be a trade-off of the accuracy and complexity.Although GA and PSO can well identify the parameters in EM, a large number of microscopic parameters will increase the computing burden.
This paper attempts to solve the abovementioned issues by utilizing a single particle model with electrolyte dynamics (SPME) based SOC online estimation framework using a sensitive parameter adjustment strategy and PF.In this thread, a proper selection procedure is designed to choose only the representative ones from SPME for parameter adjustment.There are three main contributions in this paper: (1) A simplification of the P2D model, SPME, is combined with PF for Li-ion battery SOC estimation, which involves a trade-off of both the modeling accuracy and simplicity.(2) The elementary effect test (EET) is used for the parameter sensitivity analysis (SA) of EM to improve the calculational efficiency, and only highly sensitive parameters are identified by PSO.(3) The performance of the proposed method is validated on three different driving cycles and compared with the EKF and the SPME without PF.
The remainder of this paper is organized as follows.Section 2 introduces the SPM with electrolyte dynamics and its simplification procedure.Section 3 describes representative parameter selection and the adjustment strategy with SA, EET, and PSO.Section 4 illustrates the state-space equation about the lithium-ions concentration in the solid phase, and then the SOC online estimator based on the PF and the SPME.The experimental results are exhibited in Section 5, and the conclusion is provided in Section 6.

Electrochemical Model
An accurate and effective model is the key to precisely estimating the Li-ion batteries SOC.This section will start with the classic EM for Li-ion battery, the P2D model.One flaw of P2D is that it contains multitudes of PDEs and boundary conditions, which aggravates the computational burden.Thus, P2D model is simplified into SPME, which reduces the computational burden and retains accuracy at the same time.

P2D Model
The mechanism of P2D model is shown in Figure 1a, which geometrically divides the Li-ion battery into five regions: positive electrode, negative electrode, separator, and current collectors beside the electrodes.The negative electrode is generally made of graphite for its lattice can store the lithium ions.The material of positive electrode is usually the lithium metal oxide, such as lithium cobalt oxide, lithium manganate, lithium-ion phosphate, etc.Meanwhile, a physical isolation is formed between two electrodes by the separator, allowing the lithium ions to pass through and preventing the flow of electrons at the same time.During the charging process, in response to the external electric field, lithium ions in the positive electrode de-intercalate from the positive material into the electrolyte, passing through the separator, and finally intercalating into the negative electrode, or the graphite lattice.At the same time, electrons are released from the positive electrode, flowing to the negative electrode through the external electrical circuit.During the discharging process, lithium ions de-intercalate from the graphite lattice and intercalate into the positive material in an opposite way.The input and output of the P2D model are current and voltage, respectively.For the simulation, six sets of PDEs are usually coupled together, which is capable of reflecting the lithium-ions concentration in solid phase and electrolyte.The mathematical equations of the P2D model contain the following aspects: (1) Electrolyte lithium ions diffusion equations in positive electrode, negative electrode, and separator according to Fick's second law.(2) Solid phase lithium ions diffusion equations in electrodes due to Fick's second law.
(3) Electrolyte ohm equations in electrodes and separator.(4) Solid phase ohm equations in positive electrode and negative electrode.
(5) Charge conservation equations, including positive electrode, negative electrode, and separator.( 6) Butler-Volmer (BV) kinetic equations at the surface of particles in electrodes.( 7) More details and boundary conditions are list in Table 1 [43].

Charge Conservation Equation
In Table 1, r is the coordinate axis along the radius of the solid-phase particle, where r = 0 is the center of the particle. .The symbol 'brug' is the Bruggeman porosity, which equals to 1.5 under normal circumstances.However, this model contains dozens of PDEs and their boundary conditions, which will bring multitudes of computational burdens.To reduce such a burden, a simplified model SPME is introduced in Section 2.

SPM with Electrolyte Dynamics
Regarding the complicated coupling series of PDEs in the P2D model, it is easy to understand the P2D model suffers from the computing burden when used in a real application.Thus, the simplification of P2D model has become increasingly indispensable for EM in online battery SOC estimation.The P2D model is conventionally simplified into the SPM [44] assuming each electrode is a porous spherical particle and ignoring the electrolyte dynamics.In addition, there remaining leaves a primary limitation for the SPM, that is, once the C-rates exceeds 0.5 C, the model will have a steep increase in voltage error.Hence, the SPM with electrolyte dynamics (SPME) is used in this paper for balancing the simplicity and fidelity.The schematic diagram of SPME is shown in Figure 1b.
In order to clarify the SPME, the following assumptions are given for enabling the transformation from P2D to SPME [45]: (1) The solid phase lithium ion concentration is uniformly distributed in each electrode along the spatial x-axis, i.e., it is assumed that s c (x,r,t) and n j (x, t) are constant on the coordinate x.
(2) The exchange current density 0 i (x,t) is replaced by its average value along the xaxis 0 i (t) approximately, i.e., 0 i (x,t) is assumed to be independent of x.
( .Besides, the electrolyte ionic conductivity κ is assumed to be constant.By combining the above assumptions with the mathematical equations of the P2D model, the simplified SPME can be obtained.The molar flux of lithium ions at the surface of the active particles can be written as a proportional to current: Note that the subscripts n, p, and s represent the negative electrode, positive electrode, and separator, respectively; a is the specific surface area of the active particles, where  s s a 3ε / R .The simplified electrolyte diffusion equations can be described by Equations ( 2)-(4), using the same boundary conditions as the P2D model.The simplified solid diffusion equations are expressed by Equation (5).
Following the boundary conditions: The simplified electrolyte ohm equations are given by Equations ( 7) and (8).
According to the BV kinetic equation, the electrochemical reaction overpotential can be obtained by Equation (9).
The solid potential at the current collectors is described in Equation (10).
As we can see, SPME does solve the problem of the computational burden.However, electrochemical parameters are still hard to obtain, not to mention that some parameters, e.g., D, ε , and κ , can be changed with battery aging.So, it is important to chase down the sensitive parameters and implement the parameter identification to improve the accuracy of SPME.

Global Sensitivity Analysis and Parameter Identification
As previously described, the EM usually comprises a group of microscopic parameters related to battery design.Abovementioned parameters, such as solid phase lithium ions concentration s c (x,r,t) , exchange current density 0 i (x,t) , moles of lithium ions Li n , and electrolyte activity coefficient c /a e (d ln f / d ln c )(x, t) , are simplified in SPME, that is, assumed to be constant at any positions in electrodes and invariable along the x-coordinate.In order to achieve the goal of accurate SOC estimation and voltage prediction, more than twenty parameters still need to be identified in SPME.The parameters of the EM can be divided into the following three categories, as shown in Table 2.One is the geometric parameters which can be obtained through the battery disassembling test, the thickness of electrodes and separator p L , n L , s L , and the area of electrodes S belong to this kind of parameters.Another is the performance parameters, including diffusion coefficients, solid phase conductivity σ , electrolyte ionic conductivity κ , volume fraction ε , etc., which can be obtained from the relevant literature or parameter identification method.The third category is those fixed parameters, such as the molar gas constant R, Faraday's constant F, reaction transfer coefficients α , etc.It should be noted that the performance parameters are difficult to obtain directly from experimentation, while those parameters greatly affect the accuracy of the EM.Considering the computational burden of the BMS, the sensitivity of the parameters is analyzed and further identified by a special designed algorithm in this work.

Global Sensitivity Analysis
SA is used to verify the influence of the specific parameters on the output of the model.Through sensitivity analysis, the significant parameters of the EM can be selected from numerous variables to reduce the labor of parameter identification.One way to qualify the SA is changing the parameters within a specific range to evaluate the variation of the EM's output.This paper uses the EET to analyze the parameters' sensitivity [46].Assuming that the state-space equation of the system is:  y h(x,t,u,p) (13) where z is the state vector, y is the output vector, t represents the time, u is the system input, and p is the system parameters.Assuming each parameter i p (i = 1,2,….n)is independent in the system and varies in n-dimensional space, the elementary effect (EE) of j-th trajectory, i.e., the effect of i-th parameter on the output, can be defined as: Here, we set two metrics to evaluate the parameters: the average i μ and the stand- A high value of i μ means that the parameter has a significant effect on the output, in which indicates the parameter is sensitive.Similarly, a high i s means that the parameter interacts with other parameters, or that the effect of this parameter is nonlinear.Then, these interactions should be taken into consideration for parameter identification.After understanding the sensitive parameters, PSO can be used to identify these parameters.The principle and calculating steps are described as follows.

Parameter Identification Based on PSO
PSO is originally inspired by the behavior of birds foraging in nature.Through the cooperation and information sharing between individuals of the swarm, PSO can achieve the goal of finding the optimal solution.This algorithm treats each individual as a particle without weight and volume, flying with a certain speed in the n-dimensional searching space.The flying speed is synchronously affected by the flying experience of the individual and swarm.The flowchart of PSO is illustrated in Figure 2 and detailed steps are as follows: (1) Initialization: Randomly initialize the particle swarm, the position of i-th particle in n-dimensional space is expressed as i i 1 i 2 i n X (X , X ,..., X )  , the velocity can be described (2) Best particle selection: Record individual suitability, individual optimal suitability, and global optimal suitability; the particle with individual optimal suitability is the personal best particle i i 1 i 2 i n P (P , P ,..., P )  . Similarly, the particle with global optimal suitability is defined as the global best particle i i 1 i 2 i n G (G , G ,..., G )  .
(3) Update: Update the speed and position of particle flight.
where w is the inertial coefficient, 1 c , 2 c are the learning factors representing the influence of the best particles.If the suitability is more optimal, then update the best particle; otherwise, continue the algorithm.
(4) Recursion: Determine whether the termination criterion is met.The termination criterion could be the value of suitability or the cycle number.

Li-Ion Battery SOC Estimation
After all of the electrochemical parameters are obtained and adjusted, SPME with higher accuracy is built up.We proceed with the model-based SOC estimation through PF.Two SOC estimating methods are introduced in this section, the traditional amperehour integration method and electrochemical-model-based method.

Ampere-Hour Integration Method
SOC is defined as the ratio of the current stored to the available capacity, generally expressed as percentage.EM can achieve the purpose of SOC estimation by monitoring the average solid-phase lithium ion concentration of electrodes and the lithium ion concentration on solid particles' surface.For scenarios where there is no model-based method, SOC can only be estimated by the ampere-hour integration method [47] as expressed in Equation (18).
where amp SOC is the estimated SOC value using ampere-hour integration method, 0 SOC is the value of SOC at the initial time 0 t , C is battery capacity.However, the accu- mulated errors due to the long-time operation and the uncertainty of the initial SOC cannot be ignored.Besides, this method accomplishes the SOC estimation from only the current without the correction from voltage, which is why the accuracy of SOC estimation will deteriorate.

Model-Based Method
The battery's SOC is closely related to the solid phase lithium ions concentration, where the SOC of positive and negative electrodes can be described as Equation (19) and Equation (20), respectively.The SOC of negative electrode indicates the number of lithium ions that can be extracted from the active material in the electrode before depletion [3], and the SOC of positive electrode represents the number of lithium ions that can be absorbed by the active material in the electrode before saturation.
where s,p c , s,n c is the average solid phase lithium ions concentration.
Theoretically, the SOC of positive and negative electrodes should be equal, the average of them is taken as the battery SOC considering the calculation error.To achieve the acquisition of the above two state variables s,p c , s,n c , the dispersion of solid diffusion equations is needed.Refs.[43,45] use the FDM to simplify the solid diffusion equations, dividing the solid particles into n layers evenly along the radius direction with the spherical center as the original point.The state-space equation of the negative electrode is shown in Equations ( 22) and (23).
where s ,nk c is the lithium ions concentration of the k-th layer in the solid particle, In the FDM, if the value of n is too small, the calculation error will be increased.If the value of n is too large, a large number of state variables will be introduced, and the computational burden will also be increased.In this paper, a three-parameter parabola approximation method [16] is used to simplify the solid diffusion equations, as shown in Equations ( 24)- (26).
where s cf is the average concentration flux of lithium ions in active material.In this way, the state-space equation is shown in Equation (27), and the output function can be expressed in Equation (28).

Particle Filter
After the state-space equation is build up, it is combined with PF to achieve our final goal: SOC estimation.The flowchart of PF is showed in Figure 3. PF, using the sequential Monte Carlo and recursive Bayesian estimation, is a recently rising statistical filtering method with high potential.PF uses the Monte Carlo method to estimate the state value and uses random samples in the state space to estimate the posterior probability density.In general, the steps of PF are as follows [48].(1) Initialize: Assuming the initial value 0 z follows a normal distribution 2 N(μ , s ) , defining the state-space equation of a nonlinear system as Equations ( 29) and ( 30), generate random samples (i) 0 z , (i = 1, 2,... n) and weights (i) 0 w , where the weights (i) 0 w can be distributed proportionally, or define where  r f ( ) is the probability density function (PDF) of the observation noise, and (i) k w is the new weight after normalization.Generally, we assume that the observation noise k r follows a normal distribution.
(4) Resample: During the recursion of PF, the problem of particle degradation will arise in some circumstances, that is, a small number of particles have high weight values, while the remaining large number of particles have extremely low weight values.Even there is a possibility that after a few steps of recursion, only one particle may have a non-zero weight.Therefore, Equation ( 34) is defined to determine particle degradation and resampling is used to reduce the possible problem of the failures caused by particle degradation.
The smaller eff N is, the more particles are degraded and the more urgent resampling is needed.
(5) Recursion: Repeating steps ( 2)-( 4) while k = k + 1, the predicted value of k z can be obtained by calculating mathematical expectation as follows.
The work of this paper can be summed as follows.Firstly, the EE is applied in this paper to achieve SA and to find out the sensitive parameters.Then, PSO is used to identify those sensitive parameters.Afterwards, SPME is built up with more accurate parameters, which can improve the accuracy of SOC estimation drastically.Finally, through the PF, SOC estimation is achieved.The framework of the online SOC estimation using the EMbased method combined with PF is illustrated in Figure 4.

Validation and Discussion
The Li-ion battery used in this paper is a 2.6Ah cylindrical 18650 Li-ion battery of type ICR18650-26J manufactured by SAMSUNG.The specifications of the battery are summarized in Table 3.The battery experimental testing equipment are shown in Figure 5, and can be divided into temperature control box, charge and discharge tester, and host computer.During the charge and discharge test, the battery is put into the thermostat to ensure that the ambient temperature is unchanged.The charge and discharge tester can provide various complicated working conditions according to the test requirements.Values of current, cut-off voltage, cut-off current, charge time, etc., can be controlled by the host computer, and an automatic cycle test can be conducted by setting the cycle number in the host computer.

Parameter Identification
Considering the high accuracy of parameters as the foundation of accurate SOC estimation, parameter identification is first implemented with the SA.Before parameter identification, global sensitivity analysis is applied under 1 C and 0.5 C separately to reduce the calculational burden.Hence, 14 parameters are analyzed, as shown in Figure 6, and the initial parameters of the battery are shown in Table 4.
After the SA, we can find from the results that the radius of spherical active particles s R , thickness of the electrode L, and volume fraction ε are highly sensitive parameters.
Since the thickness of electrodes and separator L indicates geometric parameters, the radius of spherical active particles s R and the volume fraction are identified using PSO.
Meanwhile, some SOC-related parameters, such as lithium ions concentration when SOC=0 and SOC =1, are identified as well.The results of parameter identification are shown in Table 5. Squares in different colors represent elementary effect in different trajectory.

Model Validation
A precise model is indispensable for a model-based SOC online estimation.Hence, for better feasibility and applicability, the simulation experiments are carried out under  different working conditions.NEDC, FTP, and UDDS working conditions are applied in this paper.The errors under different conditions are analyzed as well, as shown in Figure 7.Despite the error rises, the current increases sharply in normal circumstances.The error of voltage varies in the range of 0-0.08 V under the NEDC conditions, and it is 0-0.09V under FTP conditions, 0-0.04 V under UDDS conditions.Such results imply the phenomenon that if the discharging current increases, the model error will increase as well.However, even in FTP conditions, which contain the greater error among the three conditions, the MAE = 0.04856 shows the great accuracy of the model.The results of the experiments imply that the SPME model has sufficient accuracy for SOC estimation.

SOC Estimation Results
According to the abovementioned works, the online SOC estimation is implemented using PF with SPME.It is worth mentioning that even though the SPME itself has the ability to monitor the lithium ions concentration change inside the solid particles, without the correctness from the voltage, the errors will be significant compared to the SPME with the assistance from PF.As shown in Figure 8, the proposed method is compared to the SPME model itself and the ECM with the EKF [49].It is obvious that the error of SPME has reduced dramatically compared to the ECM method.For example, while the maximum error of the ECM method is 0.078 under NEDC conditions, the maximum error of SPME has reduced to 0.045 under the same conditions.However, without the correction from voltage, the error of SPME shows an increasing trend.By combining SPME and PF, the maximum error has reduced to 0.00064 under NEDC conditions and the increasing trend is alleviated as well.The same trend can be seen in Figure 8b,c, both in FTP and UDDS conditions.The maximum error has reduced to 0.010 under FTP conditions and the error under UDDS is less than 0.011 in normal circumstances.
As shown in Figure 9, the evaluation criteria are mean absolute error (MAE) and root mean square error (RMSE).The method combined both SPME and PF can achieve the SOC estimation with MAE = 0.0025, RMSE = 0.0039 under NEDC condition, MAE = 0.0058, RMSE = 0.0074 under FTP condition, and MAE = 0.0076, RMSE = 0.0086 under UDDS condition.The working conditions applied in this experiment present distinctive differences in discharge rate, discharge time, and the waveform of current and voltage, which shows the inestimable potential presented by the higher accuracy and wider application range of the proposed method.

Conclusions
This paper developed SPME using PF to achieve the online accurate estimation of SOC under different conditions.Considering the accuracy and fidelity of EM, this paper uses SPME to find a compromise involving both accuracy and low computational burden.After completing the battery model, the elementary test was applied for the SA of parameters to reduce the calculational burden, and eight parameters were identified using PSO algorithm.The solid diffusion equations were simplified and transformed into state-space equations to adjust the PF and achieve the SOC estimation.
The battery tests were applied under NEDC, FTP, and UDDS conditions.The experimental results show that the method introduced in this paper has high accuracy.Among three different working conditions, the proposed method has maximum MAE = 0.0076, RMSE = 0.0086, while another model-based method has maximum MAE = 0.0641, RMSE = 0.0766.Therefore, the proposed method demonstrates remarkable robustness and reliability for different battery working conditions.
The proposed method is only experimentally validated on a regular lithium-ion batteries, namely the ICR18650-26J battery manufactured by SAMSUNG.With the development of batteries, solid-state and anode-free batteries are now widely used in BMS.With some adjustment to the model structure, the proposed method can be applied to these kinds of batteries.The adjustment of model structure also can be the future trend and our next research goal.

Figure 1 .
Figure 1.Schematics of (a) the P2D model and (b) the SPME.
a α and c α are the reaction transfer coefficients of anode and cathode, separately.Normally, we define   a c α α 0.5 .f ε . is the filler volume fraction, where    f e s ε 1 ε ε .The nomenclature of parameters is shown in Nomenclature.For the parameters with superscript eff, such as ) The moles of lithium ions in electrolyte and solid phases, Li ,e n , Li ,s n , are conserved in sum.The molar fluxes n j can be written as proportional to current density I thanks to the combination of this assumption and assumption (1).(4) The electrolyte activity coefficient c /a e (d ln f / d ln c )(x, t) is constant on the x-axis and can be approximated by c /a e (d ln f / d ln c )(t)

1 z when k = 1 .( 3 ) 1 z
sequences of states and observations, separately.(2) Predict: Predictions are made for each sample point(i) 0 z based on the statespace equation of the system, generating initial sample points (i) Update: Step (2) generates the sample points (i) while the output function is generated by the observation sample points (i) 1 y , which can be used to compare with the real observation value and update the weight values of sample points.After normalization, the new particle (i) 1 z and its weight (i) 1 w are obtained by Equations (31) and (32), and the normalization function is shown in Equation (33).

Figure 4 .
Figure 4.The framework of the online SOC estimation based on electrochemical model using PF4.

Figure 7 .
Figure 7. Model simulation results of voltage, current and the voltage error under (a) NEDC conditions, (b) FTP conditions and (c) UDDS conditions.

Table 1 .
Mathematical equations and boundary conditions of the P2D model.