A Novel Dynamic Li-Ion Battery Model for the Aggregated Charging of EVs

: Implementing successful aggregated charging strategies for electric vehicles to participate in the wholesale market requires an accurate battery model that can operate at scale while capturing critical battery dynamics. Existing models either lack precision or pose computational challenges for ﬂeet-level coordination. To our knowledge, most of the literature widely adopts battery models that neglect critical battery polarization dynamics favoring scalability over accuracy, donated as constant power models (CPMs). Thus, this paper proposes a novel linear battery model (LBM) intended speciﬁcally for use in aggregated charging strategies. The LBM considers battery dynamics through a linear representation, addressing the limitations of existing models while maintaining scalability. The model dynamic behavior is evaluated for the four commonly used lithium-ion chemistries in EVs: lithium iron phosphate (LFP), nickel manganese cobalt (NMC), lithium manganese oxide (LMO), and nickel cobalt aluminum (NCA). The results showed that the LBM closely matches the high-ﬁdelity Thevenin equivalent circuit model (Th-ECM) with substantially improved accuracy over the CPM, especially at higher charging rates. Finally, a case study was carried out for bidding in the wholesale energy market, which proves the ability of the model to scale.


Introduction
Electric vehicles have gained momentum as a sustainable and eco-friendly alternative to traditional internal combustion engine vehicles [1][2][3].Their widespread adoption has been driven by various factors, including environmental consciousness, government incentives, and advances in battery technology.Importantly, EVs possess a unique capability to address the issue of renewable energy intermittency [4,5].When connected to a grid, EVs can serve as mobile energy storage units, acting as a buffer to absorb excess energy during periods of high electricity production and release it during low-production or peak demand periods [6].This bidirectional energy flow, facilitated by intelligent vehicle-to-grid (V2G) systems, contributes to more balanced and stable grid operation [7].
The electric vehicle energy conversion chain starts with the grid charging power being conditioned through rectifiers and converters to charge the battery packs [8].The stored chemical energy is then converted to electrical energy through the battery's internal reactions during discharge.This electrical energy flows through inverters, converters, and motors to provide mechanical motion energy to the wheels for propelling the vehicle.Along this conversion process, the battery management system (BMS) provides monitoring and control functions, including estimating the State of Charge (SoC) and State of Health (SoH) parameters for the controlled battery pack [9].Among the various types of batteries available for EVs, Li-ion batteries have emerged as the predominant choice [10,11].They excel compared to other battery chemistries due to their high energy density, longer cycle life, fast charging capabilities, lower self-discharge rate, lightweight design, efficiency, suitability for regenerative braking, and comparatively reduced environmental impact [12,13].The ability of lithium-ion batteries to efficiently store and deliver energy renders them well-suited for V2G applications, enhancing EVs' potential to address renewable intermittency [14].
EV smart-scheduled charging can reduce peak demand, ultimately enhancing distribution network stability and lowering operating and upgrading costs for utility companies [15].Moreover, the coordinated charging and discharging of EVs can act as a valuable grid resource capable of smoothing out fluctuations in supply and demand, storing excess energy during periods of surplus, and then injecting it back into the grid during peak demand [16,17].Thus, collective smart fleet charging and discharging platforms (in short fleet charging platforms) have the potential to transform EVs from simple consumers of electricity to dynamic energy storage and grid support assets.
According to [18], fleet charging platforms face multiple complex challenges in managing a large number of EVs while interacting with the power grid.Recent research efforts have aimed to address each of these challenges through advanced algorithms, robust optimization techniques, pricing incentives, and stakeholder coordination.Specifically, the uncertainty in EV owner behavior, like arrival and departure times, poses difficulties in optimizing fleet charging schedules.Stochastic programming and robust optimization methods have been proposed to deal with these uncertain parameters [19].On the computational complexity side, decomposition techniques and distributed algorithms are being applied to break down large fleet optimization problems into more manageable sub-problems [20].To account for EV battery degradation, models have incorporated degradation costs into scheduling objectives and constraints [21].Approaches like stochastic programming and conditional value-at-risk methods have also emerged to manage the risks of market participation [22].For coordination needs, decentralized and agent-based solutions allow EV aggregators and stakeholders to interact productively [23].Infrastructure planning for charging stations now factors in aggregator operations and uses advanced metaheuristics [24].Regarding regulations, recent policy analyses have helped to value grid services properly and avoid issues, such as deviation penalties and market exclusion [25].Lastly, for customer engagement, pricing incentives and behavioral models are being developed to influence charging patterns [26].
The aforementioned research efforts tackle many challenges.However, many aspects still need further attention, refinements, and research [18,27].One of these aspects is to improve the accuracy of models used in fleet charging platforms to minimize market participation risks and the uncertainty of the whole process.Thus, it is fundamental for modern fleet charging platforms to implement advanced battery models that are accurate and can be scaled.Thevenin equivalent circuit models (Th-ECMs), which are commonly used in battery management systems (BMSs), are an accurate description of battery dynamics [28].However, utilizing Th-ECMs with large numbers of vehicles has two scalability challenges: (a) not being linear in terms of power and energy while fleet charging platforms deal with them, and (b) relying on lookup tables for monitoring the state of charge and open circuit voltage [29].These challenges lead to the wide adoption of what we call the battery constant power model (CPM) [30][31][32][33][34], which neglects all battery polarization dynamics that may lead to a mismatch between the expectations of fleet charging platforms and the reported values from the BMSs.The aim of this research was to develop a dynamic model for Li-ion batteries that narrows down this mismatch without adding many computational requirements and which improves the fleet charging platform's insights and decision accuracy.This paper is structured as follows: Section two summarizes Li-ion battery models, focusing on Th-ECM and CPM, while the proposed model is developed in Section three.In Section four, a detailed comparative analysis between the proposed model, Th-ECM, and CPM is conducted for the commonly used Li-ion battery chemistries in EVs.Finally, Section five provides a concise conclusion summarizing the key findings of this paper.For convenience, the definitions for the abbreviations and symbols used in this paper are provided in Abbreviations.

Li-Ion Battery Models
Modeling batteries using equivalent electrical circuits is widely used since it focuses on describing the battery's electrical behavior.The Randles circuit, Figure 1, accurately models the electrochemical behavior of Li-ion batteries [35].
This paper is structured as follows: Section two summarizes Li-ion battery models, focusing on Th-ECM and CPM, while the proposed model is developed in Section three.In Section four, a detailed comparative analysis between the proposed model, Th-ECM, and CPM is conducted for the commonly used Li-ion battery chemistries in EVs.Finally, Section five provides a concise conclusion summarizing the key findings of this paper.For convenience, the definitions for the abbreviations and symbols used in this paper are provided in Abbreviations.

Li-Ion Battery Models
Modeling batteries using equivalent electrical circuits is widely used since it focuses on describing the battery's electrical behavior.The Randles circuit, Figure 1, accurately models the electrochemical behavior of Li-ion batteries [35].R is the electrolyte resistance, R is the charge transfer resistance that considers the loading voltage drop voltage over the electrode-electrolyte interface, C is the doublelayer capacitance modeling the effect of charges building up in the electrolyte at the electrode surface, and the Warburg impedance, Z , models the diffusion of lithium ions in the electrodes.Modeling Z is challenging due to its frequency-dependent characteristics.Z is modeled as infinitely series-connected RC branches [36], as shown in Figure 2. The state of charge () dependent voltage source ( ) represents the battery open circuit voltage.This dependency is usually modeled as a lookup table, which depends on the battery chemistry [29].Figure 3 shows  against  for each of the commonly used types in EVs: lithium iron phosphate (LFP), nickel manganese cobalt (NMC), lithium manganese oxide (LMO), and nickel cobalt aluminum (NCA) [37].R e is the electrolyte resistance, R ct is the charge transfer resistance that considers the loading voltage drop voltage over the electrode-electrolyte interface, C dl is the double-layer capacitance modeling the effect of charges building up in the electrolyte at the electrode surface, and the Warburg impedance, Z W , models the diffusion of lithium ions in the electrodes.Modeling Z w is challenging due to its frequency-dependent characteristics.Z w is modeled as infinitely series-connected RC branches [36], as shown in Figure 2.
This paper is structured as follows: Section two summarizes Li-ion battery models, focusing on Th-ECM and CPM, while the proposed model is developed in Section three.In Section four, a detailed comparative analysis between the proposed model, Th-ECM, and CPM is conducted for the commonly used Li-ion battery chemistries in EVs.Finally, Section five provides a concise conclusion summarizing the key findings of this paper.For convenience, the definitions for the abbreviations and symbols used in this paper are provided in Abbreviations.

Li-Ion Battery Models
Modeling batteries using equivalent electrical circuits is widely used since it focuses on describing the battery's electrical behavior.The Randles circuit, Figure 1, accurately models the electrochemical behavior of Li-ion batteries [35].R is the electrolyte resistance, R is the charge transfer resistance that considers the loading voltage drop voltage over the electrode-electrolyte interface, C is the doublelayer capacitance modeling the effect of charges building up in the electrolyte at the electrode surface, and the Warburg impedance, Z , models the diffusion of lithium ions in the electrodes.Modeling Z is challenging due to its frequency-dependent characteristics.Z is modeled as infinitely series-connected RC branches [36], as shown in Figure 2. The state of charge () dependent voltage source ( ) represents the battery open circuit voltage.This dependency is usually modeled as a lookup table, which depends on the battery chemistry [29].Figure 3 shows  against  for each of the commonly used types in EVs: lithium iron phosphate (LFP), nickel manganese cobalt (NMC), lithium manganese oxide (LMO), and nickel cobalt aluminum (NCA) [37].The state of charge (ξ) dependent voltage source (v oc ) represents the battery open circuit voltage.This dependency is usually modeled as a lookup table, which depends on the battery chemistry [29].Figure 3 shows v oc against ξ for each of the commonly used types in EVs: lithium iron phosphate (LFP), nickel manganese cobalt (NMC), lithium manganese oxide (LMO), and nickel cobalt aluminum (NCA) [37].

Thevenin Equivalent Circuit Model (Th-ECM)
Figure 4 shows the Thevenin equivalent circuit model (Th-ECM), which is the most widely used approximation of Randles circuit in which the Warburg impedance is approximated by a finite number, N (usually one to three branches are used) of RC branches [36], the usually negligible double layer capacitance, C is omitted, and the charge transfer resistance, R , and the electrolyte resistance, R , are joined into a single resistance, R .According to Th-ECM, battery terminal voltage, , is polarized from  by the diffusion voltage,  = ∑  , and the loading voltage drop,  = R  as in: where  , is the  RC branch voltage, and  is the battery terminal current.Note that the charging current is adopted as positive.
Since  is  dependent, tracking  is required.The value of  can be updated at each time step by using the following equation:

Thevenin Equivalent Circuit Model (Th-ECM)
Figure 4 shows the Thevenin equivalent circuit model (Th-ECM), which is the most widely used approximation of Randles circuit in which the Warburg impedance is approximated by a finite number, N (usually one to three branches are used) of RC branches [36], the usually negligible double layer capacitance, C dl is omitted, and the charge transfer resistance, R ct , and the electrolyte resistance, R e , are joined into a single resistance, R 0 .

Thevenin Equivalent Circuit Model (Th-ECM)
Figure 4 shows the Thevenin equivalent circuit model (Th-ECM), which is the most widely used approximation of Randles circuit in which the Warburg impedance is approximated by a finite number, N (usually one to three branches are used) of RC branches [36], the usually negligible double layer capacitance, C is omitted, and the charge transfer resistance, R , and the electrolyte resistance, R , are joined into a single resistance, R .According to Th-ECM, battery terminal voltage, , is polarized from  by the diffusion voltage,  = ∑  , and the loading voltage drop,  = R  as in: where  , is the  RC branch voltage, and  is the battery terminal current.Note that the charging current is adopted as positive.
Since  is  dependent, tracking  is required.The value of  can be updated at each time step by using the following equation: According to Th-ECM, battery terminal voltage, v, is polarized from v oc by the diffusion voltage, v d = ∑ n v n , and the loading voltage drop, v l = R 0 i as in: where v n , is the n th RC branch voltage, and i is the battery terminal current.Note that the charging current is adopted as positive.Since v oc is ξ dependent, tracking ξ is required.The value of ξ can be updated at each time step by using the following equation: where η c is the battery columbic efficiency, Q is the battery columbic capacity, and T s is the sampling period.
The use of multiple RC circuits was studied in [38], such that the voltage of each branch is described using the differential equation .
Using zero-order hold discretization, v n state transition could be expressed as: where R n , C n , and τ n are the n th RC branch resistance, capacitance, and time constant, respectively.Additionally, chargers must adhere to strict operational voltage and current limits to avoid rapid battery degradation [39,40].
where V min and V max are the battery operational voltage limits, and I chrg,max and I dis,max are the charging and discharging operational current limits, respectively.The solid region denotes the voltage limits in Figure 3.
To our knowledge, Th-ECM is the most widely used model in modern battery management systems (BMSs) [28].Th-ECM clearly states ( 2) and (3) as state transition equations with i as an input variable and ξ and v n as state variables.These state variables are used to inherently observe v oc via the lookup table while explicitly observing v using the output Equation (1).The integration of Th-ECM with fleet charging platforms requires observing, for every vehicle, stored energy (E) and terminal power (p) with: In terms of the E and p observables, the model is nonlinear with an increased number of variables and a lookup table for v oc [29].This computationally limits the scalability of the model to a large number of EVs.Thus, researchers used simplified battery models to integrate with their developed fleet charging platforms with wide dominance of what we call the constant power model (CPM).

Constant Power Model (CPM)
The CPM is a simple model that only considers the dynamics in stored energy with a single state transition equation: where η e is the battery energy efficiency.
Note that in this model, E is the only state variable with p as a single input variable.The major assumption used by the CPM is considering the E and p variables to have constant bounds [32] as: (10) where P chrg,max and P dis,max are the charging and discharging operational power limits, and E min and E max are the battery operational energy limits.Since power limits are not constant, P dis,max , and P chrg,max are average values that could be calculated as: Comparing the two models, the CPM ignores all polarization dynamics and the effect of open circuit voltage variation with the SoC and assumes constant energy efficiency independent of the charging current [32].These neglections of the CPM lead to a considerable mismatch between fleet charging platform expectations and the reported values using BMSs.For instance, by considering only the loading polarization effect for charging at the maximum rate of 1C, a 3.2 Ah NCA battery with the parameters V min = 3.3 V, V max = 4.2 V and a total internal resistance of 0.098 Ω exists [37].Then, if v oc = 4 V, the BMS would report a maximum charging power of 8.57 W according to (1) and ( 7) without breaking the operational voltage limits, while for the CPM according to (12), the maximum charging power is 12 W regardless of the current SoC or open circuit voltage.This shows the need for developing a battery model that narrows down this mismatch without adding many computational requirements, which improves the fleet charging platform's insights and decision accuracy.

Proposed Linear Battery Model (LBM)
Instead of adding more observables to the Th-ECM to track E and p, the developed model is a set of transformations and approximations applied to Th-ECM, mapping it to the power/energy domain.

Voltage Transformation
When a Li-ion battery is fully charged, its OCV is at its highest, indicating substantial stored energy, and when it is empty, its OCV is at its lowest, indicating depleted stored energy.This provides a glimpse of an intrinsic relation (transformation) between v oc and E.
For instance, a small charge injection of dq changes the battery state of charge by dξ = η c dq/Q and the battery stored energy by dE = v oc (ξ)η c dq.Then, dE could be stated in terms of dξ as: By integrating dE, the value of stored energy, E, at any ξ could be expressed as: Equation ( 14) indirectly proves that there is one-to-one transformation between E and v oc .For instance, consider two states of charge, ξ 1 and ξ 2 , with corresponding stored energies, E 1 and E 2 , and open circuit voltages, v oc1 and v oc2 , respectively.Then, if E(ξ 1 ) = E(ξ 2 ), we have: Both integrals start from zero, and v oc (ξ) is monotonic with ξ, as shown in Figure 3, therefore ξ 1 must equal ξ 2 .If two definite integrals of an increasing function are equal f (x)dx is greater than zero due to the function's monotonicity.Figure 3 also shows that if ξ 1 = ξ 2 , then v oc1 should equal v oc2 .This proves that ξ, v oc , and E have bijective relationships; knowing one of them provides all the information needed to observe the others.
To visualize the bijective relationship between v oc and E, the stored energy is calculated for each ξ according to (14), then plotted against the corresponding v oc as shown in Figure 5, with the solid portions of the curves representing the practical operational region.Surprisingly, these figures depict curves that very closely resemble linearity within the practical operational range for the four commonly used Li-ion battery types.This suggests the approximation of the relationship between  and energy as a straight line, expressed as: Achieving a linear one-to-one correspondence between  and  that is independent of the history or the present and past currents greatly simplify the model and removes the model's dependency on lookup tables.Equation ( 16) could be interpreted as a transformation between the voltage and energy.If the same transformation is applied to the terminal voltage, , as in (17), a new energy term is introduced, donated as  , which stands for the apparent energy.
Then,  operational limits could be described in terms of  using the following equation: where E = α + α V and E = α + α V are the battery operational energy limits.
Equations ( 16) and ( 17) describe the analogy between voltages and energies.For instance, during the battery charging process,  deviates from  and it is crucial to ensure that  should not exceed V .Once  reaches V , the battery appears as if it fully charged.Similarly,  deviates from  and it is crucial to ensure that  should not exceed E .Hitting E causes the battery to appear as if it is fully charged.This comparison highlights the analogy between  and .
To describe the relation between  and , Equations ( 16) and ( 17) are substituted in (1), and results in  =  + ∑ α  + R α .This shows that  is deviated from  by the loading energy deviation R α  and diffusion energy deviation ∑ α  .Then, by introducing  = α  as the deviation in energy due to  RC branch, Equations ( 1) and (3) could be restated in terms of energies and current by the following equation: Surprisingly, these figures depict curves that very closely resemble linearity within the practical operational range for the four commonly used Li-ion battery types.This suggests the approximation of the relationship between v oc and energy as a straight line, expressed as: Achieving a linear one-to-one correspondence between v oc and E that is independent of the history or the present and past currents greatly simplify the model and removes the model's dependency on lookup tables.
Equation ( 16) could be interpreted as a transformation between the voltage and energy.If the same transformation is applied to the terminal voltage, v, as in (17), a new energy term is introduced, donated as E a , which stands for the apparent energy.
Then, v operational limits could be described in terms of E a using the following equation: where E min = α 0 + α 1 V min and E max = α 0 + α 1 V min are the battery operational energy limits.Equations ( 16) and ( 17) describe the analogy between voltages and energies.For instance, during the battery charging process, v deviates from v oc and it is crucial to ensure that v should not exceed V max .Once v reaches V max , the battery appears as if it fully charged.Similarly, E a deviates from E and it is crucial to ensure that E a should not exceed E max .Hitting E max causes the battery to appear as if it is fully charged.This comparison highlights the analogy between E a and v.
To describe the relation between E a and E, Equations ( 16) and ( 17) are substituted in (1), and results in E a = E + ∑ n α 1 v n + R 0 α 1 i.This shows that E a is deviated from E by the loading energy deviation R 0 α 1 i and diffusion energy deviation ∑ n α 1 v n .Then, by introducing E n = α 1 v n as the deviation in energy due to n th RC branch, Equations ( 1) and (3) could be restated in terms of energies and current by the following equation:

Current Transformation
Up to this point, the model uses i as the input variable, while the fleet charging platforms control terminal power instead.Fortunately, the operational current limits could be stated in terms of p and E a without approximation.This can be achieved by multiplying all sides of (5) by v, then substituting vi with p and v with (E a − α 0 )/α 1 , resulting in: Since E and ξ are proven to be bijective, then (6) could be used as a state transition equation instead of ( 2), which could be stated in terms of E and P as: To linearize this equation, the ratio v oc /v could be approximated using a constant value α r .Then, the linear form of ( 18) is: Notably, v oc /v is less than the unity during charging, while it exceeds the unity during discharging.In steady-state conditions, capacitors are open circuit, thus, this ratio could be expressed as: α r is defined as the average of v oc /v.Then, to find the precise value for α r , the joint probability distribution between v oc and i should be studied.However, in this work α r is simply estimated as: These maximum and minimum values of v oc /v could be obtained by optimizing (24) while subjected to (5) and: (26) where the term Equations ( 19) and ( 20) could be rewritten in terms of power as: To linearize both equations, the term 1/v could be approximated to the constant value α v .Then, the ( 27) and ( 28) linear representations are: α v is defined as the average of 1/v.To find the precise value for α v , the probability distribution of v should be studied.However, in this work α v is simply estimated as: The validity of using the average values for the ratio v/v oc and 1/v, could be proven by comparing the overall performance of the proposed model with Th-ECM, which is conducted in the next section.

Developed Battery Model Evaluation
A comparative analysis was conducted between the dynamics of the LBM and CPM against Th-ECM, for the LFP, NMC, LMO, and NCA batteries.Battery data were obtained from [37] and listed in Table 1.For the LBM, the linear fit parameters α 0 and α 1 for each battery were estimated as described in Section 2.1.Then, α r and α v were estimated according to (25) and (31).For the CPM, the values of P dis,max , and P chrg,max were calculated by using ( 11) and ( 12), respectively.

Accuracy Evaluation of the Charging and Discharging Dynamics of the LBM
The charging and discharging simulations were performed using Google Colab by optimizing the AMPL models with a BONMIN solver.All batteries were simulated for full charge-discharge cycles.The charging dynamics, in terms of both the energy and power at a rate of 1C, are visualized in Figure 6, while the discharging dynamics are presented in Figure 7.The curves show that the LBM closely approximates the Th-ECM, whereas the CPM has a notable mismatch.
World Electr.Veh.J. 2023, 14, x FOR PEER REVIEW 9 of 15 α is defined as the average of 1/.To find the precise value for α , the probability distribution of  should be studied.However, in this work α is simply estimated as: The validity of using the average values for the ratio / and 1/, could be proven by comparing the overall performance of the proposed model with Th-ECM, which is conducted in the next section.

Developed Battery Model Evaluation
A comparative analysis was conducted between the dynamics of the LBM and CPM against Th-ECM, for the LFP, NMC, LMO, and NCA batteries.Battery data were obtained from [37] and listed in Table 1.For the LBM, the linear fit parameters α and α for each battery were estimated as described in Section 2.1.Then, α and α were estimated according to (25) and (31).For the CPM, the values of P , , and P , were calculated by using ( 11) and ( 12), respectively.

Accuracy Evaluation of the Charging and Discharging Dynamics of the LBM
The charging and discharging simulations were performed using Google Colab by optimizing the AMPL models with a BONMIN solver.All batteries were simulated for full charge-discharge cycles.The charging dynamics, in terms of both the energy and power at a rate of 1C, are visualized in Figure 6, while the discharging dynamics are presented in Figure 7.The curves show that the LBM closely approximates the Th-ECM, whereas the CPM has a notable mismatch.For the power dynamics depicted in Figures 6a and 7a, both the Th-ECM and LBM demonstrate a noticeable inflection point, signifying the transition from the constant current (CC) to the constant voltage (CV) phase.In contrast, in the CPM, the power transition indicates the end of the charging or discharging process.These deviations are more pronounced in batteries that have shorter periods of constant current (CC), such as those using lithium iron phosphate (LFP) chemistry and nickel cobalt aluminum oxide (NCA) batteries.These shorter CC phases occur because LFP has a limited voltage range, and NCA batteries have a high  value of 0.1082 Ω.
The noticeable difference in the power behavior between the CPM and Th-ECM leads to significant variations in energy patterns, as shown in Figures 6b and 7b.These differences, especially near full charge, lead to a huge mismatch between the expectations of the fleet charging platforms and the reported values from the BMS.This difference may cause the fleet management system to continuously readopt to the newly reported values, adding more computational burden.
One of the most significant trends in EV charging is the move towards faster charging speeds, which in turn requires batteries to handle higher C rates [41].For instance, Figure 8, shows the power and energy dynamics during batteries being charged at an accelerated rate of 1.5C.It is worth noting that, in comparison to the 1C, the duration of the constant current (CC) phases notably contracts, leading to a higher discrepancy between the CPM and ECM responses.For the power dynamics depicted in Figures 6a and 7a, both the Th-ECM and LBM demonstrate a noticeable inflection point, signifying the transition from the constant current (CC) to the constant voltage (CV) phase.In contrast, in the CPM, the power transition indicates the end of the charging or discharging process.These deviations are more pronounced in batteries that have shorter periods of constant current (CC), such as those using lithium iron phosphate (LFP) chemistry and nickel cobalt aluminum oxide (NCA) batteries.These shorter CC phases occur because LFP has a limited voltage range, and NCA batteries have a high R 0 value of 0.1082 Ω.
The noticeable difference in the power behavior between the CPM and Th-ECM leads to significant variations in energy patterns, as shown in Figures 6b and 7b.These differences, especially near full charge, lead to a huge mismatch between the expectations of the fleet charging platforms and the reported values from the BMS.This difference may cause the fleet management system to continuously readopt to the newly reported values, adding more computational burden.
One of the most significant trends in EV charging is the move towards faster charging speeds, which in turn requires batteries to handle higher C rates [41].For instance, Figure 8, shows the power and energy dynamics during batteries being charged at an accelerated rate of 1.5C.It is worth noting that, in comparison to the 1C, the duration of the constant current (CC) phases notably contracts, leading to a higher discrepancy between the CPM and ECM responses.

The Scalability of the LBM Compared to the ECM and CPM
To compare the scalability of the three models, a low computational overhead fleet participation model in the wholesale energy market was utilized.This simple model was formulated as: Subjected to: where j is the electric vehicle index and ranges from 1 to s, s is the fleet size, C k is the forecasted day-ahead cost of energy at instant k, k arr j and k arr j are the arrival and departure instants, and E arr j and E

The Scalability of the LBM Compared to the ECM and CPM
To compare the scalability of the three models, a low computational overhead fleet participation model in the wholesale energy market was utilized.This simple model was formulated as: Subjected to: where  is the electric vehicle index and ranges from 1 to ,  is the fleet size, C is the forecasted day-ahead cost of energy at instant , k and k are the arrival and departure instants, and E and E are charged to  To understand how well the ECM, LPM, and CPM scale, we optimized the market participation model for different scenarios of fleet sizes and battery models, employing the BONMIN solver with consistent Google Colab CPU and RAM configurations.For every scenario, the time to an optimal solution was averaged over five runs.This approach enhances the precision of our measurements.The collective findings are presented in Figure 9, showing the relationship between fleet size () and convergence time for the ECM, LBM, and CPM across 4 chemistries: LFP, NMC, LMO, and NCA.The x axis shows the size of the fleet and the y axis shows the time in seconds.
The graphs show that, as the number of vehicles increases, the time taken to reach the optimal solution also increases.The results indicate that convergence times are also dependent on battery chemistry for the Th-ECM and LBM.Notably, NMC and LMO converge faster than NCA and LFP, which may indicate a correlation between the CC period duration and the computational complexity of the problem.The data suggest that the CPM model is the most scalable, followed by the LBM and Th-ECM, with a noticeable gap.The results of this study proved our claim that the LBM has a better ability to scale than the Th-ECM.To understand how well the ECM, LPM, and CPM scale, we optimized the market participation model for different scenarios of fleet sizes and battery models, employing the BONMIN solver with consistent Google Colab CPU and RAM configurations.For every scenario, the time to an optimal solution was averaged over five runs.This approach enhances the precision of our measurements.The collective findings are presented in Figure 9, showing the relationship between fleet size (s) and convergence time for the ECM, LBM, and CPM across 4 chemistries: LFP, NMC, LMO, and NCA.The x axis shows the size of the fleet and the y axis shows the time in seconds.

Conclusions
In conclusion, this research successfully developed a Li-ion linear battery model (LBM) as a valuable tool for smart fleet charging platforms.The LBM addresses the critical battery polarization dynamics while maintaining scalability through linearity in the power and energy domain.The comparative analysis against the constant power model (CPM) and Thevenin equivalent circuit model (Th-ECM) across various Li-ion battery The graphs show that, as the number of vehicles increases, the time taken to reach the optimal solution also increases.The results indicate that convergence times are also dependent on battery chemistry for the Th-ECM and LBM.Notably, NMC and LMO converge faster than NCA and LFP, which may indicate a correlation between the CC period duration and the computational complexity of the problem.The data suggest that the CPM model is the most scalable, followed by the LBM and Th-ECM, with a noticeable gap.The results of this study proved our claim that the LBM has a better ability to scale than the Th-ECM.

Conclusions
In conclusion, this research successfully developed a Li-ion linear battery model (LBM) as a valuable tool for smart fleet charging platforms.The LBM addresses the critical battery polarization dynamics while maintaining scalability through linearity in the power and energy domain.The comparative analysis against the constant power model (CPM) and Thevenin equivalent circuit model (Th-ECM) across various Li-ion battery chemistries demonstrates the LBM's superior accuracy, especially for LFP and NCA chemistries.Notably, the accuracy advantage of the LBM over the CPM is more pronounced at higher charging rates.Moreover, measuring the conversion times for the three models tested on a simple energy market participation problem showed a noticeably large improvement in LBM scalability over Th-ECM.These findings highlight the potential of the LBM to significantly enhance the performance of fleet charging platforms, advancing the integration of electric vehicles into more efficient grid support systems.The n th RC branch resistance C n The n th RC branch capacitance τ n The n th RC branch time constant V min Minimum operational voltage limit V max Maximum operational voltage limit I dis,max Maximum operational discharging current limit I chrg,max Maximum operational charging current limit η e Energy efficiency P dis,max Maximum operational discharging power limit P chrg,max Maximum operational charging power limit E min Minimum operational energy limit E max Maximum operational energy limit α 0 The zero-order coefficient for linearly fitting E to v oc α 1 The first-order coefficient for linearly fitting E to v oc α r The mean of v oc /v under typical battery usage pattern α v The mean of Terminal power for a fleet vehicle with index j E j Stored energy for a fleet vehicle with index j

Figure 3 .
Figure 3. Open circuit voltage for various battery chemistries.The solid portions represent the practical operational regions.

Figure 3 .
Figure 3. Open circuit voltage for various battery chemistries.The solid portions represent the practical operational regions.

Figure 3 .
Figure 3. Open circuit voltage for various battery chemistries.The solid portions represent the practical operational regions.
x)dx, then it implies that b 1 must equal b 2 , as assuming otherwise leads to a contradiction, since the difference between the integrals b 2 b 1

Figure 5 ,
Figure 5, with the solid portions of the curves representing the practical operational region.

Figure 6 .
Figure 6.The charging dynamics at a rate of 1C for the LFP, NMC, LMO, and NCA batteries: (a) terminal power dynamics; (b) stored energy dynamics.Figure 6.The charging dynamics at a rate of 1C for the LFP, NMC, LMO, and NCA batteries: (a) terminal power dynamics; (b) stored energy dynamics.

Figure 6 .Figure 7 .
Figure 6.The charging dynamics at a rate of 1C for the LFP, NMC, LMO, and NCA batteries: (a) terminal power dynamics; (b) stored energy dynamics.Figure 6.The charging dynamics at a rate of 1C for the LFP, NMC, LMO, and NCA batteries: (a) terminal power dynamics; (b) stored energy dynamics.

Figure 7 .
Figure 7. Discharging dynamics at a rate of 1C rate for the LFP, NMC, LMO, and NCA batteries: (a) terminal power dynamics; (b) stored energy dynamics.

Figure 8 .
Figure 8. Discharging dynamics at a rate of 2C for LFP, NMC, LMO, and NCA batteries: (a) terminal power dynamics; (b) stored energy dynamics.

Figure 8 .
Figure 8. Discharging dynamics at a rate of 2C for LFP, NMC, LMO, and NCA batteries: (a) terminal power dynamics; (b) stored energy dynamics.

15 Figure 9 .
Figure 9. Convergence time in seconds for reaching optimal solutions across a various number of EVs for the LFP, NMC, LMO, and NCA fleets.

Figure 9 .
Figure 9. Convergence time in seconds for reaching optimal solutions across a various number of EVs for the LFP, NMC, LMO, and NCA fleets.
1/v under typical battery usage pattern C kThe forecasted day-ahead energy cost at instant k Polarization in energy due to n th RC branch p j