Model-Based Condition Monitoring of a Vanadium Redox Flow Battery

: The safe, e ﬃ cient and durable utilization of a vanadium redox ﬂow battery (VRB) requires accurate monitoring of its state of charge (SOC) and capacity decay. This paper focuses on the unbiased model parameter identiﬁcation and model-based monitoring of both the SOC and capacity decay of a VRB. Speciﬁcally, a ﬁrst-order resistor-capacitance (RC) model was used to simulate the dynamics of the VRB. A recursive total least squares (RTLS) method was exploited to attenuate the impact of external disturbances and accurately track the change of model parameters in realtime. The RTLS-based identiﬁcation method was further integrated with an H-inﬁnity ﬁlter (HIF)-based state estimator to monitor the SOC and capacity decay of the VRB in real-time. Experiments were carried out to validate the proposed method. The results suggested that the proposed method can achieve unbiased model parameter identiﬁcation when unexpected noises corrupt the current and voltage measurements. SOC and capacity decay can also be estimated accurately in real-time without requiring additional open-circuit cells. the proposed model associated with the estimator has good potential for use in large VRB systems.


Introduction
Battery storage systems have been widely applied in many fields, such as renewables storage, peak shaving, emergency back-up, load-leveling and electric vehicles [1]. Specifically, the all-vanadium redox flow battery (VRB) proposed by Skyllas-Kazacos [2,3] and co-workers shows unique attributes, including cross contamination elimination and a long life cycle [4,5]. Many efforts have been made towards component optimization [6,7] and sizing [8]. However, the efficient and durable operation of VRB systems should be under the control of a high-fidelity management system that can accurately monitor the key operating parameters of the VRB.
Battery state of charge (SOC) is defined as the percentage of residual capacity to the maximum capacity. Effective monitoring of the SOC is critical for identifying the safety margin and protecting the VRB from fast degradation caused by unsuitable over-charge/discharge; thus, it underlies any further energy management strategies. Unlike other commonly-used batteries, such as lithium-ion batteries (LIBs), a VRB allows open-circuit voltage (OCV) monitoring using additional open-circuit cells to infer the SOC. The additional components, however, increase system complexity and cost accordingly. Meanwhile, the mapping from OCV to SOC is open-loop and sensitive to measurement noise. The coulomb counting (CC) method can be used to monitor the SOC of electrochemical storage systems. As an open-loop method, however, it is vulnerable to initialization and measurement errors. For VRBs, electrolyte conductivity and spectrophotometric properties can be measured to determine . z(t) = −ηI 0 (t)/Q (3) where I 0 denotes the load current, V pi the voltage across the i-th RC branch, V t0 the terminal voltage, z the SOC, η the coulomb efficiency, and Q the maximum capacity. OCV, as a function of the SOC, is calibrated as: where c i (i = 0, 1, 2, . . . , n p ) denotes the polynomial coefficients. and capacity decay with high accuracy, high adaptivity and strong disturbance tolerance is still lacking for VRBs. This paper aims to address the aforementioned challenges and proposes an online, adaptive, model-based method to jointly estimate the SOC and capacity decay of a VRB. An unbiased parameter identification method with high noise immunity is emphasized. To begin with, a firstorder RC model was used to simulate the dynamics of a VRB, while a total least squares (TLS) method and its recursive version (RTLS) were exploited to compensate for the noise effect and identify the unbiased model parameters either online or offline. Leveraging the identified battery model, the SOC and capacity of the VRB were estimated online by combining the TLS and an H-infinity filter (HIF). The proposed method was validated by both simulation and experimental results.
The remainder of the paper is organized as follows. Battery modeling and identification based on RTLS are detailed in Section 2. The HIF-based monitoring of SOC and capacity decay is described in Section 3. Experimental details and the major results are presented in Section 4, while conclusions are drawn in Section 5.

Battery Modeling
General n-th order RC models, such as the one shown in Figure 1, are widely used to describe the dynamics of a broad range of battery chemistries [35]. In Figure 1, the voltage source describes the SOC-dependent OCV, while Rs is the ohmic resistance. The RC parallel branches aim to simulate polarization effects, including charge transfer, diffusion, and passivation layer effect on electrodes. Governing equations of the N-th order RC model are given as: where I0 denotes the load current, Vpi the voltage across the i-th RC branch, Vt0 the terminal voltage, z the SOC, η the coulomb efficiency, and Q the maximum capacity. OCV, as a function of the SOC, is calibrated as: where ci (i = 0, 1, 2, …, np) denotes the polynomial coefficients.

Parameter Identification
The model parameters associated with Equations (1)-(3) should be identified accurately. Applying Laplace transform to Equation (1) and Equation (2) gives: where y 0 = V t0 − V oc . Applying bilinear transform s = 2 (q − 1)/t s /(q + 1) to Equation (5) gives: where q −1 is the backward shift operator. Rewriting Equation (6) in the discrete-time domain gives: where It is explicit that increasing the model order elevates the computing cost due to the cubic complexity of commonly used identification methods. In this paper, the model order was determined as 1, considering the extra computing burden brought about by the need of both online parameterization and state observation. For the chosen model order, the parameter vector and regressor simplify to where t s is the time interval of signal sampling in BMSs. Once the regression model (7) is solved, the model parameters can be determined by:

Total Least Squares
The regression problem characterized by Equation (7) is typically solved by using the least squares (LS) method considering its low computing cost. However, the measurements of both current and voltage are disturbed by noise in practice, i.e., where ∆I and ∆V t are noises on I 0 and V t0 . Noise corruption causes the solution of Equation (7) to become a typical error-in-variable (EIV) problem. The LS method assumes system inputs are free from disturbances; thus, its performance for solving EIV problems is questionable. In contrast, the TLS method aims to minimize deviations on both input and output, so as to achieve the best fitting for the regression problem.

Definition 1.
The auto-covariance of arbitrary vector p k is defined by Equation (11). The cross-covariance between an arbitrary vector p k and a scalar stochastic process r k is defined by Equation (12), whilst the auto-covariance function of scalar stochastic processes r k is defined by Equation (13).
where E[•] is the expected value operator.
To ease explaining the TLS, the input and output of Equation (7) can be re-written as: where ∆ϕ k is the noise contained in real measurements.
In terms of TLS, θ k is obtained by input-output fitting with minimum perturbations; i.e., solving minimum perturbation matrices (∆Φ k , ∆Y k ) and the estimated parameter vector to satisfy: where Y k = [y k , y k−1 , . . . , y 1 ] T and Φ k = [ϕ k , ϕ k−1 , . . . , ϕ 1 ] T . In this case, the TLS solution is given by: where ||•|| F denotes the Frobenius norm. The minimization problem shown in Equation (17) can be solved alternatively by:θ where N is the order of the model used and α k is the eigenvector associated with the least eigenvalue of the covariance matrix, given by: This completes the description of the TLS for parameter identification.

Recursive Version of TLS
Singular value decomposition is generally used to extract the eigenvector. However, it is difficult for direct application in a real-time system due to the high cubic complexity. This section aims to extend the TLS to a recursive version so that a real-time update of the model parameters can be achieved.
It has been reported that α k can be obtained alternatively by minimizing the Rayleigh Quotient (RQ) cost function, i.e.,α where the specific form of the cost function is given by: where β = σ v 2 /σ i 2 . Under the criterion-constrained RQ minimization, the eigenvector of interest is defined by α k = θ T k , −1 T , so that Equation (21) can be re-expressed by: The solution of the minimization problem suggested by Equation (20) gives rise to a couple of numerical methods in the literature, although in this paper the gradient search method is applied. In this framework, the parameter vector is updated recursively by: where µ k is the adaptive step size along the update direction, which can be determined by: It is explicit that augmented covariance defined by Equation (19) can be reorganized by: Combining Equation (22), Equation (24) and Equation (25) gives: where the step size can be easily updated by: The coefficients involved in Equation (26) can be calculated by: where The covariance terms are used frequently in the proposed method. It is critical to find a computationally efficient way to calculate them in real-time. Taking R ϕ,k as an example, the covariance can be propagated recursively with low cost by: , with λ-the forgetting factor-giving more weight to recently obtained information. Specifically, a large forgetting factor gives more weight to the old data so that the regression is more stable, although the estimation is lagging in real dynamics. In contrast, a small forgetting factor keeps the alertness to the new trend, but the estimation can suffer from inevitable fluctuation. Therefore, there is always a trade-off in the selection of the forgetting factor.

Joint Estimate of SOC and Capacity Fade
Leveraging the battery model parameterized by TLS/RTLS, a state joint estimator based on advanced filtering is built in this section to observe the SOC and capacity fade of a VRB simultaneously.

H-Infinity Filter
A variety of advanced filters have been studied for the purpose of state estimation in the literature, such as an EKF, unscented Kalman filter (UKF), particle filter, and slide-mode observer. Different filters exhibit distinct features and computing complexities. An HIF is applied in this paper as it is capable of tolerating uncertainty in battery modeling and works well without the knowledge of noise statistics [36]. Fundamentally, a set of discrete-time and nonlinear state-space models should be defined firstly, i.e., where the first two equations denote the state transition and output equation; u k , x k , and d k are the input, state and output, respectively; w k and v k are the process and measurement noises with covariance of R w and R v ; h k is a unity matrix if estimating x k directly; and ξ k is a linear combination of the states of interest.
The HIF aims to optimize the estimate of ξ k , i.e., minimizing ε k = ξ k −ξ k for arbitrary x 0 , w k , and v k . A typical description of the HIF-based estimate is given by: with the cost function defined by: where is the a priori state error covariance, and S k is a user-defined symmetric positive matrix. The optimized estimate of ξ k should fulfill sup J < 1/τ according to the worst-case performance principle, where "sup" denotes the supremum (least upper bound) and τ is a pre-defined threshold. Being aware of this, the HIF solution is given by: This completes the description of HIF. With the aforementioned knowledge, the completed algorithmic procedures of HIF are summarized in Table 1.
Priori error covariance update: Update of HIF gain: Posteriori error covariance update:

Co-Estimation of SOC and Capacity
The state-space model, as shown in Equation (31), should be formulated first to adopt the HIF. In this paper, the state to be observed is formulated by x = [V p z 1/Q] T , while I and V t are defined as the system input and measurement. Combining the discretized form of Equations (1)-(3) gives: where t s is the onboard sampling time interval. Note that the third entry ofĈ k , i.e., dV oc /d(1/Q), is approximated as 0 as the SOC-OCV correlation is not sensitive to capacity decay for the VRB [37]. Leveraging the derived reference matrices, the HIF can be applied to jointly estimate the SOC and capacity fade of a VRB. The proposed method is data-driven and requires only the measured onboard current and voltage. In this case, the commonly-used open-circuit cells installed for SOC estimation can be avoided.

Experiment
A 20 (4 × 5) cm 2 VRB cell was fabricated and used for experiments. Two reservoirs were used to store the half-cell electrolytes, each of which contained vanadium electrolyte (80 mL) comprising 1.6 M V +3.5 in 4.5 M sulfate solution (GFE, Nürnberg, Germany). An anion exchange membrane (FAP 450, Fumatech GmbH, Bietigheim-Bissingen, Germany) was used to separate the two half-cells. The requested current was loaded to the fabricated cell using a battery cycler. All the data, including the load current and terminal voltage, were sampled at 1 Hz. The half-cell solutions were cycled through the cell at 60 mL min −1 under the control of a peristaltic pump, while the ambient temperature was set to 22 • C.
Hybrid pulse tests were carried out to validate the proposed method. The VRB cell was first charged fully in constant current-constant voltage (CCCV) mode. The CC current, upper cut-off voltage, and stopping current were 1 A, 1.65 V and 0.1 A, respectively. The cell was then discharged with 1 A until it reached the pre-defined SOC. After 30 min rest, a hybrid pulse current profile was loaded until the cell reached its lower cut-off voltage of 0.9 V. The load current and terminal voltage were sampled and stored in the host computer. In order to simulate the adverse condition of noise corruption, random noises of 9 mA 2 and 100 mV 2 were added on the current and voltage measurements, respectively.
Using the aforementioned experimental schedule, the capacity of the VRB cell was obtained by leveraging the coulombic counting (CC) method. As the initial SOC was accurately known, the reference SOC could also be easily determined. The measured load current, terminal voltage and calibrated reference SOCs are shown in Figure 2. As such, the proposed method was able to be evaluated with the data obtained. The reference model parameters were also extracted to validate the proposed method. In this work, a user-defined number of time points were selected during the hybrid pulse experiment. Then, around each selected time point, several batches of current/voltage measurements were sampled. Each sampled data batch was leveraged to extract the model parameters offline. Specifically, Rs was calculated by dividing the voltage jump by the corresponding current change. As accurate knowledge of OCV was available from the reference values of the SOC, Rp and Cp were readily determined by fitting the predicted terminal voltages to real measurements.

Validation of Model Parameter Identification
Identification of model parameters is a prerequisite of accurate modeling and model-based estimation. In this section, validation of the proposed method is discussed for model identification by using the data obtained from the hybrid pulse experiment, as described in Section 4.1.
The online identified model parameters associated with the corresponding reference values are shown in Figure 3. In this figure, the reference values of the model parameters are highly dynamic during the entire experiment. This can be explained by the fact that the model parameters are sensitive to both SOC and current magnitude/direction. The real-time adaption of them is thereby paramount to ensure high modeling and estimation accuracy. By using the proposed method, the erroneously initialized model parameters converged quickly after a short period of initial offset correction. The proposed method was then able to follow the change of model parameters effectively. The proposed method hence manifests itself with the capability to learn the measured current and voltage and to online adapt the model parameters. The reference model parameters were also extracted to validate the proposed method. In this work, a user-defined number of time points were selected during the hybrid pulse experiment. Then, around each selected time point, several batches of current/voltage measurements were sampled. Each sampled data batch was leveraged to extract the model parameters offline. Specifically, R s was calculated by dividing the voltage jump by the corresponding current change. As accurate knowledge of OCV was available from the reference values of the SOC, R p and C p were readily determined by fitting the predicted terminal voltages to real measurements.

Validation of Model Parameter Identification
Identification of model parameters is a prerequisite of accurate modeling and model-based estimation. In this section, validation of the proposed method is discussed for model identification by using the data obtained from the hybrid pulse experiment, as described in Section 4.1.
The online identified model parameters associated with the corresponding reference values are shown in Figure 3. In this figure, the reference values of the model parameters are highly dynamic during the entire experiment. This can be explained by the fact that the model parameters are sensitive to both SOC and current magnitude/direction. The real-time adaption of them is thereby paramount to ensure high modeling and estimation accuracy. By using the proposed method, the erroneously initialized model parameters converged quickly after a short period of initial offset correction. The proposed method was then able to follow the change of model parameters effectively. The proposed method hence manifests itself with the capability to learn the measured current and voltage and to online adapt the model parameters.  With the online adapted model parameters, the modeling results are shown in Figure 4. In this figure, the terminal voltage of the VRB can be predicted accurately with the online adapted model. The predicting error on the terminal voltage is confined well to the 10 mV error bound most of the time, although the error enlarges to some extent in the low SOC region, but can still be kept within 20 mV. The mean absolute error (MAE) and rooted mean square error (RMSE) of the voltage prediction are listed in Table 2. The low modeling error further confirms the accuracy of the identified model parameters, and provides the basis for accurate model-based estimation.

Validation of SOC Estimation
Leveraging the online adapted model, the SOC of the VRB was estimated using the HIF. The results of the SOC estimation are shown in Figure 5. In spite of the large initialization error of the SOC, the proposed method was able to converge to the reference SOC rapidly and keep track of the evolving SOC with high precision. Throughout the experiment, it was observed that the maximum With the online adapted model parameters, the modeling results are shown in Figure 4. In this figure, the terminal voltage of the VRB can be predicted accurately with the online adapted model. The predicting error on the terminal voltage is confined well to the 10 mV error bound most of the time, although the error enlarges to some extent in the low SOC region, but can still be kept within 20 mV. The mean absolute error (MAE) and rooted mean square error (RMSE) of the voltage prediction are listed in Table 2. The low modeling error further confirms the accuracy of the identified model parameters, and provides the basis for accurate model-based estimation.  With the online adapted model parameters, the modeling results are shown in Figure 4. In this figure, the terminal voltage of the VRB can be predicted accurately with the online adapted model. The predicting error on the terminal voltage is confined well to the 10 mV error bound most of the time, although the error enlarges to some extent in the low SOC region, but can still be kept within 20 mV. The mean absolute error (MAE) and rooted mean square error (RMSE) of the voltage prediction are listed in Table 2. The low modeling error further confirms the accuracy of the identified model parameters, and provides the basis for accurate model-based estimation.

Validation of SOC Estimation
Leveraging the online adapted model, the SOC of the VRB was estimated using the HIF. The results of the SOC estimation are shown in Figure 5. In spite of the large initialization error of the SOC, the proposed method was able to converge to the reference SOC rapidly and keep track of the evolving SOC with high precision. Throughout the experiment, it was observed that the maximum

Validation of SOC Estimation
Leveraging the online adapted model, the SOC of the VRB was estimated using the HIF. The results of the SOC estimation are shown in Figure 5. In spite of the large initialization error of the SOC, the proposed method was able to converge to the reference SOC rapidly and keep track of the evolving SOC with high precision. Throughout the experiment, it was observed that the maximum absolute error did not exceed 2.3%. The MAE and RMSE of the SOC estimation are summarized in Table 2.
Although the additional open-circuit cell was not used, Table 2 shows that the SOC of the VRB can still be estimated accurately via the proposed method. The proposed method is thereby shown to be highly robust against erroneous initialization and highly reliable in the estimation of VRB SOC.  Table 2.
Although the additional open-circuit cell was not used, Table 2 shows that the SOC of the VRB can still be estimated accurately via the proposed method. The proposed method is thereby shown to be highly robust against erroneous initialization and highly reliable in the estimation of VRB SOC.

Validation of Capacity Decay Monitoring
After the hybrid pulse experiment explained in Section 4.1 was performed, a capacity decay test-including 140 cycles of charge and discharge-was carried out to introduce a decay to the VRB's capacity. The test of capacity decay was performed under a constant current of 1 A for both charge and discharge, while the upper and lower cut-off voltage were defined as 1.65 V and 0.9 V, respectively. The capacity decay of the VRB cell is illustrated in Figure 6, where a significant capacity drop may be observed.
A potentiometric titration was also performed to measure imbalance in vanadium concentrations and oxidation states between the two half-cells. Specifically, 1 mL of electrolyte from each reservoir was sampled after the 140th cycle when the SOC was 23.7%. The results showed that the SOCs were 76.42% and 14.29% for the positive and negative electrolytes, respectively. Moreover, the concentrations of the vanadium species for the positive and negative electrolytes were 2.01 M and 1.02 M, respectively. This finding further confirms a substantial imbalance of vanadium species, which explains the capacity loss of the VRB.

Validation of Capacity Decay Monitoring
After the hybrid pulse experiment explained in Section 4.1 was performed, a capacity decay test-including 140 cycles of charge and discharge-was carried out to introduce a decay to the VRB's capacity. The test of capacity decay was performed under a constant current of 1 A for both charge and discharge, while the upper and lower cut-off voltage were defined as 1.65 V and 0.9 V, respectively. The capacity decay of the VRB cell is illustrated in Figure 6, where a significant capacity drop may be observed.
A potentiometric titration was also performed to measure imbalance in vanadium concentrations and oxidation states between the two half-cells. Specifically, 1 mL of electrolyte from each reservoir was sampled after the 140th cycle when the SOC was 23.7%. The results showed that the SOCs were 76.42% and 14.29% for the positive and negative electrolytes, respectively. Moreover, the concentrations of the vanadium species for the positive and negative electrolytes were 2.01 M and 1.02 M, respectively. This finding further confirms a substantial imbalance of vanadium species, which explains the capacity loss of the VRB. Online monitoring of the imbalance and capacity decay was important for determining a suitable time for rebalancing the two half-cells. In this paper, the hybrid pulse experiment described in Section 5.1 was performed again after the cycling so as to evaluate the performance of the proposed method for capacity decay monitoring.
The results of the capacity decay estimation of the VRB are illustrated in Figure 7. As can be observed in the figure, the estimated capacity is able to converge with the reference capacity following incorrect initialization for the experiment before cycling. Although the initial stage shows an overshoot, the estimation tends to stabilize gradually around the reference value. Meanwhile, the estimation shows a similar performance after the cycling experiment, although a large imbalance of both vanadium species and bulk solution occurred. Rooted in the accurate real-time estimation of capacity, the capacity decay of the studied VRB cell from 1248.5 mAh to 646.4 mAh was monitored accurately. This is insightful for real application, as knowledge of capacity decay helps the operator to determine when to start the rebalancing process to recover the capacity. The relative error of capacity estimation for the experiment before and after cycling is shown in Figure 8. As can be observed, the estimation error oscillates at the initial stage to correct the initialization error; as time evolves, the estimation converges to a steady state with reasonable accuracy. The MAE, RMSE and convergence time of the estimation are summarized in Table 3. Note that the estimation converges once the relative estimation error enters the 10% error bound. Importantly, the errors were calculated only after the algorithm converges, as the initial converging process introduces unreasonably high error. Online monitoring of the imbalance and capacity decay was important for determining a suitable time for rebalancing the two half-cells. In this paper, the hybrid pulse experiment described in Section 4.1 was performed again after the cycling so as to evaluate the performance of the proposed method for capacity decay monitoring.
The results of the capacity decay estimation of the VRB are illustrated in Figure 7. As can be observed in the figure, the estimated capacity is able to converge with the reference capacity following incorrect initialization for the experiment before cycling. Although the initial stage shows an overshoot, the estimation tends to stabilize gradually around the reference value. Meanwhile, the estimation shows a similar performance after the cycling experiment, although a large imbalance of both vanadium species and bulk solution occurred. Rooted in the accurate real-time estimation of capacity, the capacity decay of the studied VRB cell from 1248.5 mAh to 646.4 mAh was monitored accurately. This is insightful for real application, as knowledge of capacity decay helps the operator to determine when to start the rebalancing process to recover the capacity. Online monitoring of the imbalance and capacity decay was important for determining a suitable time for rebalancing the two half-cells. In this paper, the hybrid pulse experiment described in Section 5.1 was performed again after the cycling so as to evaluate the performance of the proposed method for capacity decay monitoring.
The results of the capacity decay estimation of the VRB are illustrated in Figure 7. As can be observed in the figure, the estimated capacity is able to converge with the reference capacity following incorrect initialization for the experiment before cycling. Although the initial stage shows an overshoot, the estimation tends to stabilize gradually around the reference value. Meanwhile, the estimation shows a similar performance after the cycling experiment, although a large imbalance of both vanadium species and bulk solution occurred. Rooted in the accurate real-time estimation of capacity, the capacity decay of the studied VRB cell from 1248.5 mAh to 646.4 mAh was monitored accurately. This is insightful for real application, as knowledge of capacity decay helps the operator to determine when to start the rebalancing process to recover the capacity. The relative error of capacity estimation for the experiment before and after cycling is shown in Figure 8. As can be observed, the estimation error oscillates at the initial stage to correct the initialization error; as time evolves, the estimation converges to a steady state with reasonable accuracy. The MAE, RMSE and convergence time of the estimation are summarized in Table 3. Note that the estimation converges once the relative estimation error enters the 10% error bound. Importantly, the errors were calculated only after the algorithm converges, as the initial converging process introduces unreasonably high error. The relative error of capacity estimation for the experiment before and after cycling is shown in Figure 8. As can be observed, the estimation error oscillates at the initial stage to correct the initialization error; as time evolves, the estimation converges to a steady state with reasonable accuracy. The MAE, RMSE and convergence time of the estimation are summarized in Table 3. Note that the estimation converges once the relative estimation error enters the 10% error bound. Importantly, the errors were calculated only after the algorithm converges, as the initial converging process introduces unreasonably high error.   Table 3 shows that the proposed method takes longer to converge compared to the SOC estimate. However, this will not affect real application of the method for two reasons. Firstly, the capacity change is quite slow in real-world VRB systems; thus, the convergence time can be neglected considering the long monitoring period. Secondly, the convergence is actually counted by the calculated steps. The frequency of signal sampling and calculation can be much higher in typical realtime embedded systems. In this case, convergence will likely occur within several seconds.

Discussion
In this study, we found that the capacity of the VRB dropped significantly within 140 cycles. Two reasons exist for the fast capacity loss. Firstly, the volume of electrolytes was quite small, so that the diffusion-induced imbalance across the membrane easily decreased the overall capacity. Secondly, the differential pressure drop across the membrane was not controlled during the lab-scale experiments. Hence, the rapid capacity loss observed was an extreme case for a lab-scale investigation, and is not typical for VRB systems in practical applications. Commercial VRB systems generally lower the pressure across the membrane to an expected level via effective control of the flow rate; thus, the capacity decay is much slower than that observed in this paper.
The proposed estimation method relies on a computationally affordable circuit model. By lumping together multiple electrical components to describe the dynamics with different time constants, the proposed model associated with the estimator has good potential for use in large VRB systems.

Future Work
Temperature has a large impact on the performance of VRBs. Thermal precipitation of V 5+ occurs if the temperature rises above 40 °C for typical electrolytes with 1.8-2 M vanadium sulfate in 2.5-3 M sulfuric acid, while temperature below 5 °C leads to the precipitation of V 2+ /V 3+ [38]. As our model has not been validated for describing such thermal effects, the performance of the proposed modelbased estimator will decline under such unsuitable temperatures. However, in real applications, the electrolytes are typically controlled within a specific temperature range in order for the VRB to work efficiently. In our future work, a refined model which describes the temperature effect will be   Table 3 shows that the proposed method takes longer to converge compared to the SOC estimate. However, this will not affect real application of the method for two reasons. Firstly, the capacity change is quite slow in real-world VRB systems; thus, the convergence time can be neglected considering the long monitoring period. Secondly, the convergence is actually counted by the calculated steps. The frequency of signal sampling and calculation can be much higher in typical real-time embedded systems. In this case, convergence will likely occur within several seconds.

Discussion
In this study, we found that the capacity of the VRB dropped significantly within 140 cycles. Two reasons exist for the fast capacity loss. Firstly, the volume of electrolytes was quite small, so that the diffusion-induced imbalance across the membrane easily decreased the overall capacity. Secondly, the differential pressure drop across the membrane was not controlled during the lab-scale experiments. Hence, the rapid capacity loss observed was an extreme case for a lab-scale investigation, and is not typical for VRB systems in practical applications. Commercial VRB systems generally lower the pressure across the membrane to an expected level via effective control of the flow rate; thus, the capacity decay is much slower than that observed in this paper.
The proposed estimation method relies on a computationally affordable circuit model. By lumping together multiple electrical components to describe the dynamics with different time constants, the proposed model associated with the estimator has good potential for use in large VRB systems.

Future Work
Temperature has a large impact on the performance of VRBs. Thermal precipitation of V 5+ occurs if the temperature rises above 40 • C for typical electrolytes with 1.8-2 M vanadium sulfate in 2.5-3 M sulfuric acid, while temperature below 5 • C leads to the precipitation of V 2+ /V 3+ [38]. As our model has not been validated for describing such thermal effects, the performance of the proposed model-based estimator will decline under such unsuitable temperatures. However, in real applications, the electrolytes are typically controlled within a specific temperature range in order for the VRB to work efficiently. In our future work, a refined model which describes the temperature effect will be established and combined with the proposed method to estimate the SOC and capacity decay. This will endeavor to maintain the fidelity of monitoring in the case of unsuitable operating temperatures.
The capacity decay outlined in this paper is mostly due to the imbalance of half-cell electrolytes, which is rooted in the diffusion of vanadium species and water across the separator. Nevertheless, additional irreversible capacity decay, including that caused by side effects, passivation, and air oxidation, was not analyzed in this paper. More experiments will need to be performed to further validate the proposed method for monitoring such irreversible capacity decay in the future.

Conclusions
This paper proposes a novel method for monitoring the SOC and capacity decay of a VRB, which was based on an online adapted battery model. RTLS was employed to compensate for the noise effect and online identify the unbiased model parameters. Leveraging the online identified model, an HIF-based joint observer was proposed to co-estimate the SOC and capacity decay of a VRB. Experiments were performed to validate the proposed method. Results showed that the proposed method was able to accurately track the change of model parameters under the impact of noise disturbances. The online unbiased model parameterization further ensured a high modeling accuracy. The SOC estimate was found to be within a 2.3% error bound, while the capacity decay from 1248.5 mAh to 646.4 mAh was also monitored accurately.
The proposed method observes the states of interest by making use of the measured onboard current and voltage data, as well as avoiding the need to install open-circuit cells, which increases the complexity and cost of the overall system. Hence, the method appears to be a good candidate for the management of VRB systems.