General Decoupling and Sampling Technique for Reduced-Sensor Battery Management Systems in Modular Reconﬁgurable Batteries

: The capacity and voltage rating of battery packs for electric vehicles or stationary energy storages are increasing, which challenge battery management and monitoring. Breaking the larger pack into smaller modules and using power electronics to achieve dynamic reconﬁguration can be a solution. Reconﬁgurable batteries come with their own set of problems, including many sensors and complex monitoring systems, high-bandwidth communication interfaces, and additional costs. Online parameter estimation methods can simplify or omit many of these problems and reduce the cost and footprint of the system. However, most methods require many sensors or can only estimate a subset of the elements in the module’s equivalent circuit model (ECM). This paper proposes a simple decoupling technique to derive individual modules’ voltage and current proﬁles from the output measurements without direct measurement at the modules. The determined proﬁles can achieve a high sampling rate with minimum communication between the battery management system (BMS) and the modules. With accurate proﬁles, an estimation technique can easily determine the parameters of the modules. Provided simulations and experiments conﬁrm this claim by estimating the parameters of a ﬁrst-order ECM with a parallel capacitor. The proposed technique reduces the number of sensors from 2 N + 2 to only two at the pack’s output terminals.


Introduction
Battery electric vehicles are the preferred choice in comparison to combustion-engine vehicles to reduce greenhouse gas emissions of passenger cars and the transportation sector [1,2]. Concurrently, with the increased penetration of renewable energy sources and application of power electronics, interest in stationary storages for the purpose of grid stabilization is growing [3][4][5]. Moreover, most applications, particularly within the electromobility and stationary energy storages, observe an exponential growth in energy capacity and voltage rating of the battery system [6][7][8]. The high-energy and high-power battery packs in these applications typically consist of a large number (up to several thousand) of individual battery cells which are connected in series and parallel arrangements by soldering or welding joints [9][10][11]. However, large battery packs are less fault-tolerant and prone to parameter variations [12][13][14]. Furthermore, these require a more complex monitoring system which can make a cost-efficient design of the battery management system and balancing circuitry challenging [15][16][17]. In addition, the battery pack's lifetime and usable capacity are constrained by the weakest battery cell [18,19]; when one cell in the string is empty, the whole battery pack has to be considered empty; if one cell is hot, the power has to be limited for the whole string; and if a cell is dead, the string becomes unusable.
Alternatively, bulky battery packs can be divided into a collection of smaller modules where fully-controllable power electronics replace the fixed connections to offer dynamic reconfiguration capabilities [20][21][22]. Due to the modular and fully controllable architecture, such batteries are often referred to as (modular) reconfigurable or smart batteries [23][24][25][26]. The integration of the power electronics allows for rearranging the connection of each module within the battery cell string/array to balance the individual state of charges (SOC) or temperatures of the modules [27][28][29], perform fault prevention or mitigation [30,31], increase the system efficiency [32,33], and regulate the load flow of the individual modules to improve the reliability or lifetime of the system [34]. Therefore, reconfigurable batteries inherently combine many conventional functions of generic battery management and energy conversion systems [35,36]. Furthermore, dynamically reconfiguring the battery modules can create a multilevel-output waveform which improves the output current quality (total harmonic distortion) [37][38][39] which, in turn, increases the current control bandwidth [40], while reducing the common mode voltage [41,42]. In [41,43,44], a wide selection of modular battery topologies is described for both AC and DC applications. A higher modularity in the design of the pack can also be helpful with respect to the module arrangement and design of the physical shape of the pack. For example, it can allow more freedom to optimize the weight distribution for electromobility application. The added cost of wiring and complexity can also be improved with better integration of the electronics with the modular design of the system, to keep the complexity to a minimum.
Although modular reconfigurable batteries offer clear advantages in comparison with hard-wired packs, they still suffer some drawbacks, such as increased ohmic battery losses and intricate control requirements [45,46]. With many degrees of freedom and the need for monitoring many modules, a modular battery architecture can substantially increase the complexity and system costs due to the many sensors, gate drivers, and related electronics [47][48][49]. A monitoring system based on a conventional battery management systemm (BMS) requires voltage and current sensors for each module and a dedicated communication interface between the central controller and the modules [50,51]. Aside from costly measurement and isolation circuits, the high number of signals can overload the communication bus or decrease the update rate, diminishing the appeal of such methods [52][53][54].
As a simple solution, this work demonstrates that simple measurements at the output of a reconfigurable battery system combined with the knowledge of the modules' connections within the string contain the necessary information to determine each module's voltage and current profile without directly measuring them. The derived voltage and current profiles can help to estimate the parameters of each individual battery module [47,52]. In other words, since the arrangement of the modules within a battery string is constantly changing but known, this paper presents a decoupling method to derive the voltage and current profile of the individual modules from the voltage and current sensors at the output terminals of the battery system (requiring only two system-level-sensors). These two sensors already exist for controlling the output and the string current; hence, the method is even compatible in existing designs. This idea is inspired by the capacitor voltage estimation in modular multilevel converters (MMCs) [55]. For example, as demonstrated in [56], a simple sliding-mode observer to track the module voltages in the discrete-time domain minimizes the communication requirements. Similar other approaches exist in the literature [57][58][59][60][61] which focus on the capacitor voltage estimation in MMCs, with different levels of accuracy and complexity. Although these approaches are advantageous in reducing the system cost and complexity, they are only applicable for tracking the average terminal voltage of capacitor-based modules and mostly disregard other types of energy storages with more complex equivalent circuit models (ECM) [61][62][63][64]. Therefore, these methods do not scale well for batteries since these have a nonlinear voltage-current behavior, which can only be represented by intricate ECMs [47,52,65].
From the BMS perspective, the ECM parameters are crucial as most monitoring functionalities, including charge or health monitoring and power capability estimation, depend heavily on an accurate battery model [41,[66][67][68][69]. Studies have shown that these parameters vary heavily with the operating condition, SOC, and state of health (SOH), making developing accurate models even more challenging [70,71]. In a conventional BMS, an accurate battery model for specific SOC and currents is developed as a look-up table which is gradually updated as the battery ages [72,73]. However, methods that accurately update these parameters require constant monitoring of the voltage and current of the battery [74,75]. In recent years, there have been some attempts to improve the accuracy of the battery models for modular reconfigurable batteries that consider the peculiar behavior of batteries in response to a pulsating current [76,77]. Concurrently, attempts to improve some of the reduced sensor estimation methods considering the ECM of the batteries are made [47,52,78]. One of the simplest approaches is modeling the battery behavior with a variable voltage source in series with an equivalent circuit resistor. With the help of current variations, the equivalent series resistance of the module and the open-circuit voltages (OCVs) can be determined by tracking the voltage drop for different operating points [78,79]. While the method can appeal to some applications, the zero-order ECM needs to offer more information [80,81]. However, determining the parameters of higher-order ECMs for the batteries requires having the module's accurate voltage and current profile during the current variations, which goes beyond the capabilities of existing methods.
This paper proposes a decoupling technique that accurately determines every module's voltage and current profile without any direct measurement at the module level. The decoupling technique can achieve synchronized sampling frequencies as high as the module's switching rate, which is sufficiently large to estimate all the required parameters of the battery. The fixed and synchronized sampling rate simplifies the parameter estimation. Moreover, this paper develops a simplified second-order ECM for batteries with an optional parallel capacitor. The proposed model, combined with an estimator, demonstrates the applicability and accuracy of the developed decoupling technique. Compared to conventional BMSs, as shown in [78,82], the proposed technique can reduce the number of sensors from 2N + 2 (sensors for module current and module voltage) in a battery string with N modules to only two at the output terminals (output voltage and current of the pack). Additionally, the decoupling technique contains only simple scalar mathematics, does not require any computationally demanding calculations, and can be easily integrated into the switching function of the modules. As it is possible to track the exact voltage and current profile during the current variations, it is possible to estimate the ECM parameters of individual modules which capture the dynamics of the battery system more accurately [53,62,83]. Proposing a new estimation technique is not the main innovation of this paper. Consequently, a state-of-the-art estimator provides proof of concept and shows how the developed decoupling technique can be integrated with an estimation method. Similarly, other estimation methods with different levels of computation and complexity can also be combined with the proposed decoupling technique [84][85][86][87][88][89]. In addition to reducing the number of sensors and the size of the overall monitoring system, the necessary data communication can be significantly reduced. The effectiveness of the proposed technique is verified through simulations and experiments.
The remainder of this paper is organized as follows: Section 2 introduces the modularreconfigurable battery system and its working principles. Section 3 proposes an algorithm for estimating the electrical equivalent circuit parameters of the battery module. Section 4 discusses the simulations and the obtained experimental results. Finally, the last section concludes the paper with final remarks. Figure 1 depicts the modular reconfigurable battery where v b j and i b j are the battery voltage and current of the j th module, respectively [90]. The output voltage of the j th module is denoted by v t j . They form the vectors v b , i b , and v t where j ∈ {1, 2, . . . , N}. The modules have been connected in series and supply the load through a low-pass filter formed by L and C. v p and i p as the terminal voltage and current of the whole reconfigurable battery are independently measured for controlling the output voltage and current of the system. The battery management system or the central controller can also use these signals for monitoring purposes. Two switches and their anti-parallel diodes (in half-bridge topology) in parallel with the battery form one module which can offer bypass and series connections, as illustrated in Figure 2 for the j th module [47]. The upper switch in the j th modules is S 1 j , and the lower switch is S 2 j . The vector S denotes the state of all the modules. The j th element of S equals one if the j th module operates in the series mode and zero if it is bypassed. Therefore, considering switches and diode resistances, the general equation of the string voltage follows

System Description
where "×", ".", and "T" denote the matrix product, the normal product, and the transposition operator, respectively. r on is the on-state resistance of the switches or diodes in each module, R is the parasitic series resistance of the inductor L, Considering a single battery with an optional parallel capacitor (C p j ) to reduce the ripple as Figure 3a illustrates, the equivalent circuit based on the first-order ECM for the battery is shown in Figure 3b. Other parameters shown in Figure 3b are the open-circuit voltage (v oc ), and the ohmic resistance (R int ). The left RC network (R pa , C pa ) models the chemical reactions (e.g., polarization process) inside the battery that typically has a time constant of few-hundred milliseconds up to several seconds [91]. Recent research on the negative effect of high-frequency current ripples is inconclusive [92][93][94][95][96][97]. Still, in most cases, there is a parallel capacitor to reduce the high-frequency content of the pulsating load of the battery. The exact value of C p j depends on the switching frequency and operating conditions, but it is typically within the millifarad range for a few kilohertz switching rates [98].  As the size of C p j in comparison to C pa j is considerably smaller, the RC network modeling the chemical reactions acts as a constant voltage for C p j . Therefore, the ECM shown in Figure 3b can be approximated with Figure 3b. Similarly, the battery's opencircuit voltage can vary due to SOC or temperature, but within a few hundred milliseconds, it can be assumed constant. Hence, the ECM is simplified to Figure 3c which describes a second-order RC network. In the following, we consider this derived second-order RC network to present the decoupling technique. If there is no parallel capacitor with the batteries, C p j −→ 0 in all the derived equations.
Based on the the circuit depicted in Figure 3, the vector of the battery output voltage, where v oc is the vector of the open-circuit voltage, [v oc 1 v oc 2 . . . v oc N ] T . It generally changes according to the SOC, SOH, and some other predictable and unpredictable factors. v 1 is the vector of the voltages across the RC network formed by R int and C p . Similarly, v 2 is the vector of the voltage drops across both the R pa and C pa in Figure 3d. Replacing the differential equations for the RC networks of the j th module, The differential equations in (3) can be discretized using the backward Euler method for the vectors v 1 (k) and v 2 (k) as where k denotes the k th sample. Q 11 , Q 12 , Q 21 , and Q 22 are four vectors whose elements are C p j T s , Substituting (4) in (2) When the module is in series mode, the battery current is equal to the string current, and when in bypass mode, the current is zero. Therefore, the vector of the battery current can be expressed based on the vector of modules' connections (S) and the string current following Substituting (7) and (6) in (1) with some mathematical manipulations, the output terminal voltage for the k th sample is given by Equation (8) describes v p as a function of the modules' current and previous states (S, v b ) in addition to the string current. This relation can be leveraged to develop the method for decoupling the voltages of the modules from the string voltage.
Besides (8) as the key equation for the decoupling algorithm and the estimation algorithm, the average model of the half-bridge converter can provide a better understanding of the modules and the effect of the pulse-width modulation (PWM) on their behavior. The average of the output voltages of the modules v t is denoted by v t and follows v t = m T × v b (9) where m is the vector of the modulation indices, that is, m = [m 1 m 2 . . . m N ] T with amplitudes between one and zero [99]. The phase-shifted PWM is used for the modular configurable battery [100]. The average of the output voltage of the pack can be stated as i b denotes the average of the battery currents (i b ) during each switching cycle (T sw ). As long as the duration of the switching cycle is negligible compared to the time constants of the RC networks in the ECM shown in Figure 3, the effective current passing through the two RC networks equals the average module current per This approximation holds as long as the switching cycle follows

Decoupling Algorithm
Using (8), we can write the output voltage of the reconfigurable battery for the k th sample (v p (k)) as a function of The decoupling algorithm relies on the vector of the modules' states (S). If the j th module operates in the series mode, the j th battery will affect the output voltage as it can clearly be deduced from (8). However, the voltage of the j th module would not contribute to the overall pack voltage if it is in the bypass mode, i.e., S j = 0. Considering (8), if only the j th module changes its connection from the series to bypass (S j = 1 −→ 0) within the (k − 1) th and k th interval, the voltage variation observed between v p (k) and v p (k − 1) is due to the j th battery. In other words, the step in the output voltage of the pack exactly before and after the j th module changes its connection from series to bypass can be attributed to the terminal voltage of the j th module. Concurrently, using (11) and (7), the instantaneous and average current profile of the modules can be also determined.
Mathematically, this voltage difference can be stated as As long as multiple modules do not change states concurrently, which is the case for almost all available PWM-based modulations, the profiles of all modules can be updated using the corresponding negative edges of the voltage steps at the output. After each negative voltage step, the new data point is added to the corresponding vectors of the module per Performing the decoupling, particularly during load variations, can provide the necessary current and voltage profiles to estimate the parameters of the modules. Therefore, an online decoupling of voltage and current profiles during load variations can be beneficial, but the estimations can run slower during the idling intervals of the controller.
This simple yet effective sampling method forms the main principle of the decoupling technique. It enables the BMS or the central controller to determine the batteries' voltage and current profiles without using excessive sensors for each module. As the switching rate of the modules is relatively fixed in most modulation strategies (such as phase-shifted PWM) and the sampling rate of the voltage and current would be also synchronized to the switching actions and can be as high as the switching rate of the said module. This is a rather unique advantage of the proposed technique over the existing methods. Whereas in conventional methods, the sampling rate is usually restricted to tens of milliseconds due to the limited data bandwidth of the communication interface between the controller and the modules.
The decoupling process can be performed independent of the modulation process, or it can be integrated with it. If the estimator has only access to the measured output signals and the module states through a memory buffer, the sampling and decoupling process can be performed offline. This is usually the case if the monitoring and modulation routines are independently designed. However, it is possible to integrate both processes and perform an online synchronized sampling. Figure 4 shows the flowchart of the algorithm for an independent monitoring unit. First, the number of available samples, here, n, is determined depending on the buffer size. Then, the vector ∆S is calculated for each sample. It is the difference between S(k − 1) and S(k). Afterward, if only one of the elements inside ∆S (e.g., the j th element) equals one and the other elements equal zero, the voltage and current profile of the j th module will be updated per (14). The process is repeated for all of the n samples. Therefore, this process can even be added later on to a working system or even be performed in parallel. The sampling process can be integrated with the modulation using the interrupt routine to achieve synchronized sampling with minimum memory and computational requirement. This is achieved by two consecutive sampling before and after the negative edge of the lower switch (S 2 j control signal. The first sampling should be sufficiently before the start of the switching that the sampling is completed, and the second sampling should provide sufficient time for the completion of the switching interval and be short enough that the other modules' states can be assumed constant. For example, MOSFETs have a switching duration below a microsecond and IGBTs may need several microseconds, which both are much shorter than the time constants of the battery [101][102][103].

The Estimation Algorithm
The estimation process can be defined as the following optimization problem: where lb and ub are the lower and upper boundaries of the parameters which are supposed to be estimated. The function f is defined as where n is the number of times the falling edge has been detected for the j th module. ∆v mea j (k) is the determined voltage of the battery using the proposed decoupling technique, whereas ∆v est j (k) is the estimated voltage of the battery according to the estimated parameters.
As this is a simple optimization problem, many techniques can be used [104][105][106]. In this research study, the interior-point method is applied to solve it. As proposing a new optimization method is not the intent of this work, we refer to the literature for further detail [107].

Results
The simulation and experimental validation of the proposed method are shown in Figures 5-13. Matlab/Simulink serves to simulate a modular reconfigurable battery with seven modules. In practice, the sbRIO-9726 rapid-control prototyping board from National Instrument carries out the voltage modulation, records the modules' switch states (S) and measures the output current, i p , as well as the output voltage, v p . The switching rate of each module is set to 10 kHz, providing a maximum sampling rate of 10 kHz per module. The number of the modules (N), L, R, C are equal to 7, 1 mH, 50 mΩ, and 5 mF, respectively. To determine the true values in the experimental setup, battery modules are discharged with 0.025 C and the terminal voltage dynamics are constantly recorded. At extremely low-discharge rates, the terminal voltage of the battery is a good approximation of the open-circuit voltage after relaxation. Moreover, the terminal voltage is measured before and after the experiment when the batteries have reached equilibrium to determine the exact SOC of the modules. In this manner, the parameter sets of the internal resistance and the RC network parameters of the batteries, according to Figure 3, are separately derived through pulsed current at different SOCs and current amplitudes. Figure 5 shows an example of the output voltages of the modules for one complete switching cycle. It intuitively shows how they are combined to form the output voltage of the pack. Using the phase-shifted PWM, the output voltage of the pack is N-times the switching rate of each module. The negative and positive edges of the modules' voltages is clear on the profile of the pack's voltage. The voltages of the modules are approximately 8 V to 8.2 V, which are normalized for simplicity. The output voltage of the pack at each instant is the sum of all the modules' voltages. The time division in Figure 5 is 10 µs. The normalized voltage of the modules can approximate their corresponding state at each instance. As the module toggles from S j = 1 to S j = 0, the decoupling algorithm detects the falling edges in the vector S and records a new sample for the j th module. The sampling instances are marked with red arrows in Figure 5. The complete implemented process, from measuring the output voltage of the system to creating individual voltage profiles of each module, is shown in Figure 6. The estimation of the ECM conventionally is performed using pulse charging approaches, here, the load is changed to create an external disturbance for the system. After the decoupling procedure, as the voltage and current of the whole system (i p and v p ) contain a high level of noise, it is possible to reduce the noise level using a digital low-pass filter as long as the added delay due to the denoising procedure is negligible in comparison to the ECM time constants.
The following equation offers a typical example used to reduce the switching noises of the measured current i p at instant k For example, a q value of two means the output of the denoising function at the k th sample is the average value of the same sample along with the two samples before and after the main sample. In Figure 6, the actual and denoised profiles of the modules' voltages are depicted in the bottom left corner, respectively, in gray and blue, which clearly show the improved waveforms.
To study the accuracy of the estimations using the determined profiles, two systems with balanced and imbalanced modules are considered. Table 1 depicts the main parameters of the systems. In the case of imbalanced modules, each module has completely different parameters (±50%) to test the convergence of the system.  Figure 7 shows the estimation results for the ECM parameters of balanced modules (R int , C p , R pa , C pa , v oc ) using the resulted profiles through the decoupling method. Similarly, Figure 8 shows the results for the imbalanced scenario. The chart shown for each parameter contains the estimated value, the real value, and the absolute error for every module.
The absolute estimation error is calculated as:  The estimation error for all the parameters are shown in Figures 7 and 8 are below 4% which authenticate the usability of the derived voltage and current profiles for parameter estimation.
In both cases, the errors for estimation of the capacitors are well below 3%. Additionally, the absolute errors for the estimated resistances are less than 6 mΩ. Additionally, there is no change in the accuracy of estimations for the imbalanced system compared to the balanced one.
The estimation results for each iteration of the optimization process for the five parameters are presented in Figures 9 and 10. The initial estimations of the parameters are selected randomly with a ±50% uncertainty boundary. The estimated open-circuit voltage rapidly converges to the true value after a few iterations, but as the other two RC networks are relatively co-dependent, more iterations are required for the parameters to converge. As the depicted values are normalized in Figures 9 and 10, a result of one means the estimation is correctly converged. The optimization stops when there is no visible improvement in the cost value. It should be mentioned that the optimization process does not need to be online and can be completed during the idle intervals of the controller or even be off-loaded to the cloud services, but that falls outside the real focus of this paper [108].
Three scenarios help in evaluating the ability of the proposed estimator under different conditions: 1.
The first scenario investigates the behavior of a balanced system where all the modules have relatively similar states (their SOC and SOH as a case in point) and parameters.
There are only small inevitable inherent tolerances among their parameters due to manufacturing tolerances.

2.
In the second scenario, the uniform parameters of the modules are disturbed to simulate a possible imbalance in the system. The ohmic resistance of the sixth battery module (R int 6 ) is intentionally increased by externally connecting an extra 0.43 Ω resistance in series to the battery. The added ohmic resistance can emulate an aged or possibly faulty battery. Additionally, f sw is changed from 10 kHz to 1 kHz to study the effect of sampling frequency. Changing the battery resistance challenges the ability of the proposed technique to track the asymmetric voltage profiles and evaluate the estimation accuracy for parameters which are beyond normal ranges. 3.
The third scenario focuses on existence of imbalance in the modules' voltages. Therefore, battery modules 2 to 7 are charged to 8 V, whereas the first module is manually discharged to a lower open-circuit voltage (approximately 7.7 V). As a result of a higher charge depletion, other ECM parameters of Module One also change slightly. This scenario can emulate a SOC imbalance in the system, which can challenge the ability of the proposed decoupling technique to distinguish the profiles of different modules. Moreover, in the last two scenarios, estimation results of the other modules remain similar, showing that an existing imbalance in one module does not significantly impact the accuracy of the decoupling method for other modules.
The relevant experimental results for each scenario are shown in Figures 11-13, respectively. For each of them, the voltage profile of one of the modules (one, six, and one, respectively) is measured, denoised, and then plotted on the left side. Concurrently, the voltage profile of each module is derived and the parameters are estimated according to the derived profile. Afterwards, the corresponding output battery voltage are calculated based on the estimated ECM parameters. On the right side of each figure, the measured waveform is compared with the calculated one based on the estimated parameters. The matching results confirm the feasibility of the proposed decoupling technique.
The tables of Figures 11-13 show the measured parameters of the corresponding battery modules using the proposed decoupling technique. As seen, the estimated ohmic resistance in the second scenario is clearly increased. This confirms the authenticity of the results using the proposed method. Similarly, the open circuit voltage of the module in the third scenario is correctly estimated.

Conclusions
Proposing a feasible solution to determine the voltage and current profiles of modules for parameter estimation purposes in a reconfigurable battery is the main purpose of this paper. After detailed modeling and analysis of the reconfigurable batteries' operations, this paper proposes a decoupling technique that can easily be integrated into the monitoring routine. The resulting profiles for each module can be used to estimate the ECM parameters of the batteries with any state-of-the-art estimation technique. Conventionally, the estimation is accomplished based on the direct measurement at the module terminals, which requires a high number of sensors and a high-bandwidth communication interface, increasing the cost and complexity of the system. The proposed decoupling algorithm does not need any measurement at the module terminals and only uses the already available sensors (system-level measurements) at the output of the system. Therefore, the number of sensors compared to the conventional technique is decreased from 2N + 2 to only two. Additionally, no further communication between the controller and the modules is necessary, which reduces the constraints on the communication interface. The proposed concept can easily be implemented, reducing the system cost and improving the system's robustness. In addition to the advantages, the method can only estimate the time-constants that are relatively larger than the switching cycle of the modules. Therefore, it cannot be applied to estimate the sub-millisecond time constants of the battery. However, these time constants are rarely utilized in a BMS for estimating the module parameters. Increasing the band-width of the decoupling technique is an interesting goal that warrants further study as a future work. Additionally, switching noise can have a minor effect on the accuracy of the voltage and current profiles, but it can be mitigated through appropriate digital filters, e.g., Kalman filters. The practical results as well as the simulations confirm that the decoupling algorithm works properly. Combined with a suitable estimation technique, the decoupled voltage, and current profiles achieve very good accuracy; the estimation algorithm can accurately calculate the parameters and converge to the desired values.

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

Abbreviations
The following abbreviations are used in this manuscript: