Concurrent Real-Time Estimation of State of Health and Maximum Available Power in Lithium-Sulfur Batteries

: Lithium-sulfur (Li-S) batteries are an emerging energy storage technology with higher performance than lithium-ion batteries in terms of speciﬁc capacity and energy density. However, several scientiﬁc and technological gaps need to be ﬁlled before Li-S batteries will penetrate the market at a large scale. One such gap, which is tackled in this paper, is represented by the estimation of state-of-health (SOH). Li-S batteries exhibit a complex behaviour due to their inherent mechanisms, which requires a special tailoring of the already literature-available state-of-charge (SOC) and SOH estimation algorithms. In this work, a model of SOH based on capacity fade and power fade has been proposed and incorporated in a state estimator using dual extended Kalman ﬁlters has been used to simultaneously estimate Li-S SOC and SOH. The dual extended Kalman ﬁlter’s internal estimates of equivalent circuit network parameters have also been used to the estimate maximum available power of the battery at any speciﬁed instant. The proposed estimators have been successfully applied to both fresh and aged Li-S pouch cells, showing that they can accurately track accurately the battery SOC, SOH, and power, providing that initial conditions are suitable. However, the estimation of the Li-S battery cells’ capacity fade is shown to be more complex, because the practical available capacity varies highly with the applied current rates and the dynamics of the mission proﬁle.


Introduction
Growing demands from original equipment manufacturers (OEMs) and users for lighter batteries with a higher capacity has lead to research into and development of new energy storage technologies such as Lithium-Sulfur (Li-S) batteries. This technology represents a suitable alternative to nowadays broadly used Lithium-ion (li-ion) battery cells as they are characterized by six times higher theoretical specific energy and two times higher theoretical energy density than traditional Li-ion batteries [1]. Furthermore, Li-S batteries are expected to be cheaper, as they use sulfur as a raw material, [1] and safer [2]. However, there are still obstacles to their successful large-scale commercialization. The first type of obstacles is related to the actual cell performance and lifetime, which are hindered problems such as the fast capacity fade during cycling [3] and the high self-discharge at high state-of-charge (SOC) levels [4]. While the chemical challenges are addressed at cell development or material level [5], there are also practical challenges that have to be addressed by system design, in order to operate the Li-S cells in a safe and optimal way due to their complex chemistry [6].

State-of-Charge and State-of-Health Estimation
The battery management system (BMS) is vital in the practical use of batteries, as it takes care of the safe and efficient operation of the battery [7]. Typically, advanced BMSs include state estimation functions, which provide information about state-of-charge (SOC) and state-of-health (SOH), which are crucial for effective and reliable battery function and also useful to the user.
Many different methods for the state estimation have been proposed in the literature. Frequently-mentioned methods for SOC estimation are: ampere-hour counting, open-circuit voltage (OCV) based, model based, impedance based, static battery characteristics based, fuzzy logic and machine learning based estimations [8]. Each of these methods has its advantages and drawbacks, which makes them suitable for different battery technologies. SOH estimation is usually based on various capacity and/or impedance estimation techniques and also correlations of those quantities to some observable effects [7]. These estimation techniques are usually tailored to specific battery chemistries, as they might have unique attributes and behaviour. It is possible to find many works which have been done on widely used Lead-acid [9][10][11], Nickel-metal Hydride (NiMH) [12][13][14] and Lithium-ion [11,15,16] batteries. However, these methods, successfully applied to other battery chemistries, are not directly applicable for Li-S batteries, due to the particularities of this technology. For example, the Coulomb counting method is not suitable due to the high self-discharge [4] and due to the variable shape and length of voltage charging curve, which is caused by the polysulfide shuttle and 'history' effect [17]. Moreover, the Li-S batteries have a unique shape of the open-circuit voltage [18], which is flat for the interval 0 to 70% SOC. Furthermore, a voltage 'dip' occurs during discharging between the high and the low voltage plateaus, which causes non-monotonic voltage trend. This reduces their observability and it prevents the use of any open-circuit voltage based technique.
To the present day, in the field of Li-S batteries, there has been proposed and demonstrated only a limited number of functional approaches to SOC estimation. A technique based on Kalman variant estimators (extended Kalman filter (EKF), unscented filter and particle filter) was presented in reference [19]; in that work, the filter was applied to the battery model, developed in reference [20], in a classical manner, meaning that the SOC was estimated according to a combination of Coulomb counting and voltage response of the battery. The SOC was estimated with a root square mean error (RSME) of 0.0114 for mixed pulse discharge profile scenario and 0.0217 RSME for New European Driving Cycle (NEDC) driving profile for the correct initial state. In a scenario, when the initial SOC was lower than the actual battery state (0.7 instead of 1), the estimated SOC by EKF converged quickly towards the reference state. However, when the initial SOC was set lower (i.e., 0.6), the estimated state did not converge over one full discharge. Moreover, when the estimation was applied to partially charged cell (initial reference SOC was 0.6), the estimator showed a poor accuracy and in the case of initial estimator SOC being set to 1, it also did not converge over one full discharge. Another approach for state estimation of Li-S batteries, discussed in reference [18], is based on the observability of the SOC from the battery parameters, such as open-circuit voltage and internal resistance. This method can provide advantages as these identified parameters might be easier to track and express during the laboratory degradation tests than the complex battery voltage response used in previously mentioned approach. The first step, following this approach, is to implement an online battery parameter identification algorithm. Three identification algorithms: gradient descent, genetic algorithm and prediction error minimization were investigated and compared in reference [21], where also the SOC estimation of a NiMH battery based on an adaptive neuro-fuzzy inference system (ANFIS) was demonstrated, which can be in a similar way applied to the Li-S battery. Furthermore, the parameter identification can be also based on the extended Kalman filter, which was presented for Li-S batteries in reference [22], and it was directly used in a dual Kalman filter manner for the SOC estimation of Li-S batteries. Accordingly, we can summarize that there exists already some base ground for the SOC estimation at Li-S batteries. Nevertheless, the robustness of estimation needs to be improved, especially in cases when the initially set SOC does not match the actual battery SOC, then the convergence has to be ensured. Moreover, the estimation needs to be validated not only at the fresh cell, but also at an aged cell, as battery parameters change during aging.
The area of SOH estimation for Li-S batteries has not been covered at all. Therefore, in this paper we introduce a combined SOC and SOH estimation algorithm, which is based on the parameter identification method introduced by Propp et al. [22]. The proposed SOH estimation method allows to estimate the capacity fade and resistance change caused by the battery aging without the use of any degradation model, which is especially beneficial at the current stage of reduced knowledge and practical experience with the Li-S batteries.

Maximum Available Power Estimation
Besides the SOC and SOH estimation, another function required of the BMS is estimation of the available power. Especially in electric vehicle applications, there are certain situations when the high power is required from the battery, for example during the acceleration. In these situations, the safety operation limits can be reached quickly, which will result into limitation of the provided power and its sudden decrease, causing an unexpected and possibly dangerous situation for a driver. Therefore, the estimation of the available power is applied to predict the maximum power which can be provided, typically in period of 1 to 20 seconds. In this way, the provided power can be limited from the beginning, causing the smooth operation and expected continuous decrease later on [7].
There are usually considered to be three main methods for the estimation of the available power. The first is based on a characteristic map, which is stored in a memory of the BMS. It is composed of static relationships of power towards the other battery states and quantities and it is obtained offline from the battery tests. During aging, some of the battery parameter changes, and thus the characteristic map has to be updated. The advantage of this method is its easy implementation and simplicity. However, it does not reflect the dynamic states of the battery accurately and it might be difficult to adjust it during the aging due to lack of the reference measurements (at maximum power). Moreover, it might required considerable space of the memory to store all the data with the inter-dependencies [7,23].
The second method is based on a dynamic battery model. The maximum power is then predicted according the accuracy of the battery model, reflecting the dynamic states. The obvious advantage of this method in our context is that the battery parameters can be identified online. Therefore, the battery parameters would be adjusted according to aging and also specific chemistry effects (a history effect of the Li-S batteries) to provide more precise estimated [7,24].
The third method, considered for this purpose, uses ANFIS. However, an additional filtering is required for training on power pulses. Moreover, the real-time applicability on a low cost hardware is questionable [7,25].
To the authors' best knowledge, there are no examples of maximum available power estimator for Li-S batteries in the literature. Thus we target this area in this work and present such estimator. This paper is structured as follows. The methodology for modelling and state estimation is described in Section 2. This includes definitions of battery states in terms of state of charge, health and maximum available power; and a description of the considered Li-S battery model. The section introduces a standard implementation of an Extended Kalman filter, an online parameter identification method, a maximum available power estimation approach and a metric for an estimation evaluation. The specific implementation of the SOC and SOH estimation algorithms is presented in Section 3. The results for SOC, SOH and maximum available power estimation are described and shown in Section 4. The paper is concluded in Section 5.

State Definitions
States of the battery receive a great amount of an interest. However, they are not always perceived in the same and unique way. There is more consensus about the state-of-charge definition, which is usually understood as an amount of charge, which is available to be extracted from the actual battery capacity [7,26]. Sometimes, the nominal capacity term is used instead of the actual capacity [27]. However, that might lead to an incorrect assumption of neglecting the capacity change during the battery life. Moreover, the amount of the extractable charge in the battery varies with the temperature or the applied current. The capacity dependence of the current is typically described by Peukert's law [28]. Extending the Peukert equation by a temperature dependency was described in [29].
As an indicator for the state of health of the battery, the most commonly used is a change in the actual capacity, a capacity fade [10,13]. Some methods estimate the capacity fade according to its relationship to the change of impedance, which is usually growing during the battery aging [30]. Furthermore, the impedance was accommodated as an another indicator for the SOH [15,26] and it is also related to the maximum available power of the battery, as the impedance growth contributes a lot to the power fade. However, other factors might also be considered for the SOH estimation [27], as for example the porosity of the electrodes [31] or the terminal voltage [32].
The maximum available power is the power actually achievable from the battery with respect to its safety limits and internal and external conditions [7].
To summarize, it is important to clearly state which definition is followed in the work. Definitions in this work are based on [7,26], where the initial capacityQ cap is the maximum extractable charge from the fully charged battery at the beginning of life under the specific conditions (i.e., temperature, C-rate, cut-off limits). The actual capacity Q cap represents the maximum extractable charge from the fully charged battery at the actual battery age under the specific conditions. The definitions for the initial internal resistanceR 0 and the actual internal resistance R 0 correspond to similar definitions to the capacity definitions, but to the internal resistance instead of the capacity.
State-of-charge: SOC represents a present amount of charge Q present , which is possible to extract from the battery, related to the actual capacity. It is noted as χ and expressed as: (1) State of health: SOH consists of two quantities representing capacity fade η Q and internal resistance change η R . The SOH can be written as: The capacity fade (η Q ) is computed as the actual capacity divided by the nominal capacity at the beginning of life: The internal resistance change (η R ) is computed as the internal resistance at the beginning of life divided by the actual internal resistance: Maximum available power: P dis.max /P ch.max is considered in this work to be limited only by current and voltage limits. Generally, it can be written as: where both current I and voltage V are limited to not exceed their own safety limits.

Li-S Battery Model and Its Parametrization
For the estimation of defined states, it is necessary to relate the behaviour of the Li-S cells to them. For the relationship of parameters over the SOC, the Thevenin electrical circuit network, illustrated in Figure 1a, and its behavioural representation, illustrated in Figure 1b, were applied. The details of the derivation of this model are presented in our previous work [20]. In this paper, we introduce the model briefly as the algebra of the following sections rely on it. Tests were based on discharging OXIS Energy 3.4 Ah long life Li-S cells in different temperatures with varying current rates shown in Figure 1c. The original idea of the application of different current pulses was to capture the open-circuit voltage (V OC ), the internal resistance (R 0 ) and the transient behaviour of the cell (R 1 , C 1 ) over its discharge range in addition to exploring the effects of the current rate. However, it was found that the rate dependencies of the model parameters are small and that they could be neglected without decreasing the model accuracy significantly. Since the model's intentional usage is also the state estimation, polynomial functions in the form of were fitted describing the parameters over the SOC. While fitted polynomial functions are generally a standard method, it is difficult to reproduce the harsh transitions between the high-and low-plateau of the Li-S cells with them. Therefore, two polynomial functions were combined for the V OC and R 0 and joined together at the transition point. The model implements combined functions for the V OC and R 0 Since the linearisation of the applied Kalman filter methods need differentiable functions, γ combines both polynomials smoothly where the conditions a, b and c stands for the different ranges of the function: where m is a scaling factor for the maximal gradient of the sinusoidal function, determining the transition range between both polynomials and c represents the point where both functions are equally represented. The resulting state space model has the forṁ which relies on the nonlinear functions of V OC , R 0 , C 1 , R 1 over the SOC. The dynamic states x = [x 1 x 2 ] T of the system are the voltage over the RC circuit V p and the SOC χ, calculated through Coulomb counting. The corresponding state space representation gives Alternatively, the core Thevenin battery model equations are written as: The behavioural battery model representations, described in [20], introduces parameters representing the dynamic bandwidth (Ω), the total steady-state resistance (R int ) and the dynamic fraction of the response (ρ), instead of R 0 , R 1 and C 1 . The parameters inter-relations are defined as: By substituting Equations (15), (16) and (17) into Equations (13) and (14), the following behavioural model equations are obtained:

Standard Extended Kalman Filter Implementation
The EKF implementation described in detail can be found for example in [33,34]. The EKF consists of two steps: time update and measurement update. During the time update, the state is estimated ahead in time (predicted) and, during the measurement update, the state is adjusted (corrected) according an actual measurement at that time. The representation of the EKF update equations with approximation of the state and measurement vector without noise values can be written as: Time update (prediction): Measurement update (correction): The process and noise measurement covariance matrices are formed as

Online Parameter Identification
In contrast to the standard state estimation, where a model identified offline is used to predict the system's output and the error between the prediction and measurements is used to correct the states of interest, an online parameter identification method can be used to identify the parameters directly. However, this means that the identification has to be fast enough to run on common BMS hardware and that the identification results have to be reliable with different current profiles. While the main concept of online identification was proposed in [35][36][37] for common Li-ion batteries, Propp et al. [22] adapted the identification for Li-S batteries. Since the principles of this work are used here as well, they are also introduced here briefly. The online parametrization fits the parameters with the derivative of the terminal battery voltage V L of the Thevenin model Generally, the OCV of the battery is dependent on the SOC, the operating temperature (T) and usage history (h). Therefore, the corresponding definition of V OC can be described as a function of SOC, T, h, whose derivation leads tȯ For the representation of the battery in a discrete manner, usually time steps of less than one second are used. Slow changing parameters therefore can be eliminated from the equation for simplification.
Thus, the SOC δSOC δt ≈ 0 , the change of the OCV over time δV OC δt ≈ 0 , the temperature changes δT δt ≈ 0 and the long-term history δh δt ≈ 0 can be approximated as 0. It is furthermore assumed that the internal resistance variation are negligible (Ṙ 0 ≈ 0) for small periods. The resulting simplified relationship for the terminal voltage over timė was then changed to a behavioural interpretation to improve the identification process. For the terminal voltage,V and for the voltage drop over the RC circuit, By re-writing the behavioural interpretation of the equation for the terminal battery voltage to bring V p on one side, substituting it into Equation (31) and including the result forV p in Equation (30), the behavioural state transition equation can be derived aṡ For the online identification with a Kalman filter based algorithm, the system is now expressed by its parameters as the state vector With the corresponding state transition functions from Equations (31) and (33), the input current I L = u and the previous considerations the state transition functions can be populated as where f 1 to f 6 are defined as: Since the measured terminal voltage of the battery is represented by the second state, the measurement equation is The Jacobians for the presented functions f 1 to f 6 , linearising the system around the current mean for the extended Kalman filter, are therefore: C can be obtained as This method has been proven to show robust identification results as long as there are some dynamics in the current profile. To apply the introduced online parametrization, the system has to be implemented in an extended Kalman filter algorithm, which is introduced in the following section.

Structure of the Maximum Available Power Estimation
The implemented maximum available power estimation is based on the dynamic battery model approach, when the battery parameters are identified online. The estimation is derived accordingly: The Laplace transform of equation (19) is which can be rearranged as and again toV Let us assume that the future current has a constant value, I L . Note that Substituting this in equations (18) and (44) gives Now, Ω s(s + Ω) We can take the inverse Laplace transform: or, equivalently, Let us assume that at some point T ≥ 0 seconds in the future, we want to know the minimum and maximum currents such that V L (T) ∈ [V min , V max ]. For convenience, we can write and In practical applications, V , R > 0 since V OC > V p0 and ρ < 1. We can then write To which becomes where which becomes where In practice, we are likely to also want to ensure that the current is constrained to a specified window, I L ∈ [I min , I max ]. The most conservative of these will be the most important. Effectively, the current that gives the maximum discharge power is given by with the corresponding power given by Similarly, the current that gives the maximum charge power is with the corresponding power given by

Estimation Evaluation
The estimation algorithms are evaluated according their absolute maximum error: and their absolute mean error: The subscript "Est" stands for estimated states and the subscript "Ref" represents reference values of the states. The reference states are obtained from an offline analysis of the measurements. Coulomb counting provides reference for χ. Moreover, obtained data from the parametrization (Section 2.2) are used as a reference for capacity fade η Q and the resistance change η R .

Test Procedure and Model Structure
The laboratory tests were performed as the base for tunning and evaluation of the estimators. The current profiles (in this work, specifically considered mixed pulse discharge, NEDC and Urban Dynamometer Driving Schedule (UDDS)) were applied to the battery, which was placed at the thermal chamber with the controlled temperature environment and it was connected to the battery test station.
The tests on the fresh cell were done in a MACCOR battery test system (supplied by MACCOR inc.,Tulsa, OK, USA) and the tests on the aged cell were performed in a Digatron battery test system (supplied by Digatron Power Electronics, Aachen, Germany). The measured quantities of the battery, such as current and voltage, were recorded and they were used as an input to the model with the estimators. The temperature for the simulations was considered constant of 20 • C. The layout of the test procedure is in Figure 2. The model layout with the estimators is shown in Figure 3. Current, voltage and temperature are the inputs. The first extended Kalman filter estimates online the circuit model parameters in terms of behaviour model, which are also translated to Thevenin model parameters. The estimated parameters V OC and R 0 are fed, together with current and temperature, into the SOC and SOH estimator, which returns the outputs of χ, η Q and η R . Furthermore, all the behaviour circuit model parameters are used as the input for the maximum available power estimator, which estimates the maximum available power in the specific time period ∆t for charging and discharging, together with maximum available current and voltage.

Modelling and Structure of the Filters
The states of the interest are SOC (χ) and SOH (η Q , η R ). Therefore, the state vector is constructed as: The SOC change is based on Coulomb counting. However, no process models for SOH are used, as the change of SOH is assumed to be very slow in comparison to SOC. η Q and η R are assumed to be observed from the increasing mean value of the error. The change of state then takes form with consideration of no noise as:ẋ The considered measurement vector has form: where V OC is open-circuit voltage and R 0 is internal resistance. They are obtained through online parameter identification. The function h(x − k , 0), relating the estimated measurements (estimated parameters V OC and R 0 ) to the estimated states, is based on the model and the fitted polynomials described in Section 2.2.
The full C matrix for EKF is: With the assumptions that the capacity fade does not have influence on V OC and R 0 , together with that resistance change does not influence V OC , the C matrix can be rewritten as: c 11 and c 21 are obtained by the derivative of h(x − k , 0) with respect to χ. c 23 takes the form: The full A matrix for EKF has the form: which, without use of any process model for η Q and η R , results in: The observability of such system can be determined from the 'observability grammian' having a full column rank. The grammian takes form as follows in (76), where c 11 and c 23 have to be non zero: The linearization of the state change is considered as: where T s is a simulation step size.

Numerical Values
The initial value for a posteriori error covariance is set as: The measurement noise covariance R represents a measurement noise of V OC and R 0 and is set to be emphasizing higher trusts in the more consistent open circuit voltage identification than the fluctuating internal resistance (Figure 4a). The process noise covariance Q, related to the uncertainty of the process models, is set as: Here, the arbitrary parameter tuning is based on slow changes, smaller variance values, for the χ and η Q states and faster parameter changes, or higher variances, for η R , that need to follow relatively quick changes at high states of charges (Figure 4c). Furthermore, the states are assumed to be uncorrelated as the diagonal Q matrix indicates. It shows that, when the cell ages, the inflection point between the high and low voltage plateaus is moving to the right (to lower SOC); (b) the mixed pulse discharge test profile and the estimated results for the fresh cell, when the discharge starts already at the low voltage plateau. If the SOC ini is also at the low voltage plateau, the state converges, otherwise they do not; (c) the mixed pulse discharge test profile and the estimated results for the aged cell. The estimates converge for all states. The accuracy is higher with the initial estimate closer to the reference state; (d) battery model identified parameters for the aged cell at various mission profiles. The absolute available capacity obtained during the discharge causes the shifting of the inflection point at the V OC and R 0 curves due to the SOC relativity to the available capacity. At the end of the discharge, there is a rapid increase of the resistance, dependent on the current rate and the dynamic of the mission profiles.

Fresh Cell-SOC and SOH Estimation
The mixed pulse discharge profile was applied to the fresh cell and the measurement data was fed to the estimator, as illustrated in Figure 2. The reference and estimated states are shown in Figure 5 and their error values are summarized in Table 1. When all the states were initialized with one (= fully charged, fresh cell), the average errors were 0.0080, 0.0125 and 0.0140; and the maximum errors were 0.0306, 0.0559 and 0.7714 for χ, η Q and η R , respectively. However, the maximum error value of η R is caused by the initialization of the parameter estimator and it is reached only at the very beginning. When the current is firstly applied, the resistance quickly converges close to the reference value already after one second. Therefore, the 'real' maximum error of the estimator performance is 0.2606 and the same metric is going to be used further on to evaluate the 'real' performance of the estimator and not the error caused by the parameter estimator initialization. In this ideal case, one can see that the state estimator provides relatively high accurate results with the average error below two percent and the maximum error below six percent. The least reliable state is the η Q . Moreover, the η R does not always match the reference, which is probably due to the settling times to reach convergence and also due to the parameter's dependence on the current. In Figure 4a, one can see that the parameters of V OC and R 0 used for the state estimation are identified very closely to the reference, while the remaining parameters of R 1 and C 1 vary a bit. The R 1 is identified to be slightly higher in the low voltage plateau and the C 1 is identified slightly lower in the high voltage plateau in comparison to the model. Table 1. The errors of the state estimation for the fresh cell with various initial conditions. It shows the capability of the estimator to converge to the states in the cases when the SOC ini is set to match the reference SOC value or to be mismatched, but still on the high voltage plateau as the SOC reference was. If the SOC ini was instead on the low voltage plateau, the estimator was not able to converge.

Initial Conditions Average Errors Maximum Errors
The next step was to see the convergence capability of the applied Kalman filter when the initial state values are wrongly set. In the first case, all the states are set to be 0.7 in the beginning. From Figure 5, it is apparent that the χ and the η R converged almost instantly. The χ estimation error was reduced already in the first step of the simulation, after one second, from 0.3 to 0.1264. The error of η R in the beginning of the simulation was 0.4714, which in 10 seconds of relaxation reduced to 0.4307 and, after the first second of the current being applied, it dropped to 0.1301. The η Q took a longer amount of time and it continuously settled in the middle of the high voltage plateau, after around 0.15 of SOC is discharged.
So far, the initial SOC (SOC ini ) was always set to be on the high voltage plateau. When the SOC ini was set to 0.6, which is on the low voltage plateau, the filter was not able to converge the SOC over the whole discharging period. Moreover, the originally well set η Q diverged in the second half of the discharge and the η R was highly off in the first half of the discharge.
To see what happens when the situation is the opposite, the reference SOC is on the low voltage plateau, while the SOC ini is at the high voltage plateau, the simulation was run for only a part of the mixed pulse discharge profile with the initial reference SOC to be 0.6764, which means that the discharge starts already at the low voltage plateau. The obtained results are shown in Figure 4b. When the SOC ini was one, all the states diverged. However, when the SOC ini was 0.6, the states were again estimated very well. Therefore, it is very important that the SOC ini is present at the same voltage plateau as is the reference or real SOC. In the case when the SOC ini was incorrectly set at the high voltage plateau, the η Q increased rapidly and steadily over the value of 2. Such a 'nonsense' value can be used as an indicator for the incorrect estimation and re-set the estimated SOC value.
This stability issue of the estimator based on the initial conditions of the SOC, whereas the cell is at a high or low voltage plateau can be solved in the following way. The initial conditions for the parameter identification, with a meaning of [V OC V L V p Ω ρ R int ], are adjusted to be: where V meas is the actually measured voltage at the cell. Subsequently, the initialization of the SOC/SOH Kalman filter contains the following condition: where the c stands for the SOC transition point between the high and low voltage plateau. Moreover, the V OC,max_lp is the maximum V OC value obtained from the Li-S battery model. In such a way, the estimator starts with the SOC ini value, which will lead to the convergence. Of course, if the more accurate guess of the initial condition exists, it shall be inputted.

Aged Cell-SOC and SOH Estimation
The goal of the SOH estimation is to estimate and track the health of the cell. The cell under investigation was aged by 40 cycles of continuous charge (0.34 A) and discharge (0.68 A) at the elevated temperature of 50 • C; and stored at the shelf for approximately six months at room temperature. Afterwards, the continuous discharge, mixed pulse discharge and driving cycles of NEDC and UDDS were applied to it at 20 • C. The obtained capacity from the mixed pulse discharge test at 20 • C was 1.86 Ah, which results in η Q being 0.69.
The estimation for the mixed pulse discharge profile, shown in Figure 4c, works well, similarly as for the fresh cell, only with a lower accuracy. The error values are summarized in Table 2. The χ was slightly over estimated at the low voltage plateau with the highest misfit being 0.0550 for η Q ini = 1 and 0.0461 for η Q ini = 0.69. This χ over estimation is probably caused in the 'offset' present in the estimated value of the η Q at the low voltage plateau. In this region, the estimated value in an average was 0.74 for η Q ini = 1 and 0.71 for η Q ini = 0.69, while the reference was 0.69. It is also seen that, when the initial value of the η Q was closer to the reference, the accuracy was again improved. It is important to note that only the model of the fresh cell is implemented in the estimator and it has no knowledge, only assumptions, introduced in Section 3.1, of how the parameters will change with the change of states. In Figure 4a, the actual identified parameters are shown. The low voltage plateau has shrunk more than the high voltage plateau; and therefore the inflection point is moved to the right. Otherwise, the V OC is not changed. However, the R 0 seems to be lower than at the fresh cell, except the very low SOC, where it has increasing tendency.
Besides the mixed pulse discharge profile, four mission driving profiles were also applied to the cell and proceeded to the state estimation. The profiles were: NEDC with the maximum current of 1.2 A (=NEDC12) and of 2.9 A (=NEDC29), and the UDDS also with the maximum currents of 1.2 A (=UDDS12) and of 2.9 A (=UDDS29). The profiles with the estimation results are presented in Figure 6. The χ more or less follows the reference. However, the confidence is reduced and the error can be again related to the estimated η Q . From the reference η Q , one can see that it varies from 0.8110 to 1.1965, while the value obtained from the mixed pulse discharge profile was 0.69. The different numbers come from the obtained capacity from the specific tests, which was 1.8592, 2.8221, 2.2029, 3.2499 and 2.4172 Ah for the mixed pulse discharge, NEDC12, NEDC29, UDDS12 and UDDS29, consequently. The actual capacity varies due to its high dependence on the applied current and the dynamics of the profile, which was reported for the Li-S batteries in [38] and investigated and modelled in [39,40]. How it results into the identified parameters is shown in Figure 4d, where the inflection point at V OC and R 0 curves is moving according to the available capacity over the discharge because the SOC is relative to the actual capacity according to (1). Moreover, the dynamics related to the applied current rates and the relaxation periods are visible at the rapidly increasing resistance at the end of discharge. Therefore, the estimated η Q in our work is related only to the actual profile. If the applied profile remains similar, then the η Q will straightforwardly indicate the capacity fade ongoing in the cell. In order to have a referent, or so-called independent, actual capacity, which would be related to the capacity fade independent of the mission profile, the dynamics such as transport limitations have to be included in the model.
The estimated η R seems to have a similar character for all the driving cycles and also the mixed pulse discharge at the aged cell. It is shown in detail in Figure 7. Additionally, the average of the aged cell curves is plotted to highlight the observable trend. At the high voltage plateau, high SOC, the η R increases sharply and creates a peak around the value of 2. During the flat low voltage plateau in the middle SOC region, there is a smaller increase to around 1.2 and, at the end of discharge (SOC close to 0), there is a decrease down to around 0.6. If this behaviour would be observable over the other cells and mission profiles during aging, it could be also used to couple the resistance change to the capacity fade as it is often done at other battery chemistries. Figure 6. The driving cycle mission profiles and the estimated results for the aged cell. The χ roughly follows the reference. The accuracy of the χ is dependent on the estimated η Q , which varies with the total discharged capacity of the cell under different mission profiles, due to their different current rates and dynamics. The η R is in general followed well and is similar between the different driving cycles. Figure 7. The estimated η R at the aged cell compared to the fresh cell. The changes in the η R due to aging are very similar at any mission profile. In the high SOC region, the η R is rapidly growing. Over the middle SOC region, the η R slightly grows. Finally, in the low SOC region, the η R increases.

Maximum Available Power Estimation Validation
In order to validate the maximum available power estimation, the modified mixed pulse discharge profile was applied to the Li-S cell with considered maximum current pulses and it is shown in Figure 8. After the first set of the discharging pulses with currents of 0.29, 1.45 and 2.90 A, consequently, the discharging pulse with 6.8 A was applied and, after the second set of the discharging pulse, the charging pulse with 1.7 A was applied. When the voltage limits of 1.5 V for discharging and 2.45 V for charging were reached during these maximum current pulses, the operation mode switched from the constant current mode to constant voltage mode in order to obtain the maximum power in those conditions. We have considered two cases for the power estimation: first, the instantaneous power at T = 1 second, shown in Figure 8 and then for the 'accelerating' pulse for T = 10 seconds. The accuracy of the estimation is evaluated from the whole 30 second maximum current pulse for the instantaneous power. For the 'accelerating' pulse, only the last 20 seconds of the pulse are considered, as the estimated power is for 10 second application of the maximum current, so the estimated power can be compared just after the first 10 seconds of the applied pulse.
The absolute maximum and average error are computed similarly to equations (65) and (66), consequently, and they are presented in Table 3. Under the selected boundaries of maximum allowed voltage and current limits, the cell was able to provide 14.78 W for discharge and 4.14 W for charge during the applied pulses. The absolute maximum errors seem to reach high significant values related to the maximum cell performance. However, the absolute average errors are fairly acceptable. It is also important to keep in mind that the model was derived and accounted for maximum discharging current 2.9 A and no charging model or parameters are implemented. Therefore, the performance of the maximum available power estimator, while applying current of 6.8 A for discharge, which is more than two times that for the model that was parametrized, and 1.7 A charging current, is considered as more than acceptable.

Conclusions
In this work, a novel SOC and SOH estimation method for Li-S batteries has been introduced, together with a specifically tailored model for maximum available power estimation. The proposed estimation algorithm is based on two extended Kalman filters, where the first filter is used for online parameter identification of the battery electrical circuit model using voltage and current measurements. Subsequently, the estimated battery parameters are used in the second extended Kalman filter to estimate the SOC, the SOH, (in terms of both capacity fade and the resistance change) and the maximum available power for both charging and discharging pulses.
The robustness and the accuracy of the proposed SOC and SOH estimator has been analyzed using a mixed pulse discharge profile firstly for a fresh cell and then for an aged cell. The estimator does not include any process model of the Li-S battery degradation and it is solely based on the error correction. The estimator has shown the capability to track the SOC, the capacity fade and the resistance change for an aged cell with the worst case average errors being 0.0827, 0.1900 and 0.1796, respectively. However, the accuracy of the SOC is dependent on the estimated capacity fade, which varies depending on the current rates and dynamics of the considered mission profiles. Therefore, there is a need for the concept clarification of the actual capacity at the Li-S batteries according to their dynamics and history. Moreover, the model for the available capacity change would be helpful in order to provide information about the absolute capacity fade. In terms of the absolute capacity fade, the resistance change could be used as a supportive indicator, as it has shown a more consistent trend between the fresh and the aged cell. However, the question of how it evolves under different degradation conditions has to be still answered.
The implemented maximum available power estimator was presented and validated by comparing the obtained estimation results with the experimental tests. The proposed model proofed to be robust also when currents of 6.8 A were applied, which is more than two times the value of the current, which had been used for the model parametrization. Furthermore, the available power estimation algorithm provides relatively accurate results for the case when the available power should be predicted during battery charging, a case that was not considered either during the parametrization stage.

Acknowledgments:
The authors would like to thank OXIS Energy (Abingdon, Oxfordshire, UK) for supplying the Lithium-Sulfur battery cells.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations and symbols are used in this manuscript: