Novel Ordinary Differential Equation for State-of-Charge Simulation of Rechargeable Lithium-Ion Battery

: Lithium-ion battery energy storage systems are rapidly gaining widespread adoption in power systems across the globe. This trend is primarily driven by their recognition as a key enabler for reducing carbon emissions, advancing digitalization, and making electricity grids more accessible to a broader population. In the present study, we investigated the dynamic behavior of lithium-ion batteries during the charging and discharging processes, with a focus on the impact of terminal voltages and rate parameters on the state of charge (SOC). Through modeling and simulations, the results show that higher terminal charging voltages lead to a faster SOC increase, making them advantageous for applications requiring rapid charging. However, large values of voltage-sensitive coefficients and energy transfer coefficients were found to have drawbacks, including increased battery degradation, overheating, and wasted energy. Moreover, practical considerations highlighted the trade-off between fast charging and time efficiency, with charging times ranging from 8 to 16 min for different rates and SOC levels. On the discharging side, we found that varying the terminal discharging voltage allowed for controlled discharging rates and adjustments to SOC levels. Lower sensitivity coefficients resulted in more stable voltage during discharging, which is beneficial for applications requiring a steady power supply. However, high discharging rates and sensitivity coefficients led to over-discharging, reducing battery life and causing damage. These new findings could provide valuable insights for optimizing the performance of lithium-ion batteries in various applications.


Introduction
Rechargeable lithium-ion batteries are widely used as a power source in many industrial sectors ranging from portable electronic devices to electric vehicles and power grid systems [1][2][3].In the context of energy management and distribution, the rechargeable lithium-ion battery has increased the flexibility of power grid systems, because of their ability to provide optimal use of stable operation of intermittent renewable energy sources such as solar and wind energy [4][5][6].In addition to storing renewable energy for later use, lithium-ion batteries are the major component dictating the performance of electric vehicle (EV) performance.Furthermore, their integration into EVs offers a promising solution to reducing CO 2 emissions.
Rechargeable lithium-ion batteries offer many opportunities, such as high power density, long life cycle [7,8], higher cell voltages [9,10], lower maintenance requirements [11,12], higher charging speeds with lower self-discharging rates [13,14], and no memory effect compared with the previous generations of rechargeable batteries [15].However, even if the rechargeable lithium-ion battery seems to provide a suitable solution to the limitations of renewable energy production and utilization, there are also some existing issues related to its technology.More explicitly, one of the fundamental problems with this device is battery management systems (BMSs) [16], and among all the parameters of BMSs, state of charge was found to be the most significant parameter [17,18].The state of charge is one of the crucial evaluation parameters of BMSs that confirms the optimal charging or discharging profiles of batteries [2].Moreover, the state of charge describes the remaining capacity in a battery energy storage system and is indicated in the form of a percentage, [%].Meanwhile, one can identify the charging or discharging status of a battery by checking its SOC value [19,20].High precision in state-of-charge estimation is a must for power battery control strategies [21,22].Indeed, to ensure the optimal charging/discharging profiles of the battery, the SOC should satisfy the following condition: where SOC min is the minimum state of charge and SOC max is the maximum state of charge.These conditions play an important role when designing a controller to decide when to charge and when to discharge a battery [2].Furthermore, these box constraints on the state of charge also help the controller to cut off the power to prevent the battery from damage due to over-discharging and overcharging [23].In ref. [24], for instance, it is reported that the overcharging of lithium-ion batteries leads to the generation of dendrites due to the deposition of lithium ions, which induces the generation of a large amount of heat due to internal short circuit, and the vaporization of the electrolyte solution.Another research work published has shown experimentally that when the lithium-ion battery is discharged to the cut-off voltage and then forced to discharge more than 20% of the state of charge, it will cause irreversible internal short-circuit damage to the battery [25,26].In this respect, the authors concluded that the possibility of short circuit has a non-linear relationship with the depth of over-discharging.In general, determining the SOC of rechargeable lithium-ion batteries is a complex task, due to the non-linear reaction dynamics caused by the electrochemical characteristics of lithium-ion batteries [27,28].This means that the SOC can only indirectly be calculated by measuring some physical parameters, such as charge/discharge current, terminal voltage, and internal resistance [29], or by using a mathematical model and algorithm [30,31].However, the physical parameters of lithium-ion batteries are usually affected by some uncertain factors, such as the open-circuit voltage (U oc ) of the battery, temperature variation, battery aging, and overall efficiency of the battery systems [32].This induces strong non-linear behavior in battery systems, which significantly increases the difficulty of estimating the battery state of charge.Therefore, accurately estimating the state of charge is crucial to keeping the lithium-ion battery working under safe conditions and maximizing its performance.

Some Existing Methods of Estimating the Lithium-Ion Battery State of Charge
At present, battery state-of-charge estimation methods can be divided into three categories, namely, experimental methods, data-driven methods, and model-based methods [33][34][35].

Experimental Methods
Experimental methods refer to the methods to explore battery characteristics through battery experiments.The first experimental method is referred to as the Coulomb counting method [36].This technique estimates the SOC in the battery by integrating the current transferred to the battery over time, to determine the remaining charge relative to the maximum battery state of charge.Coulomb counting is often utilized in combination with model-based methods due to its relatively simple implementation.However, this method suffers from inaccuracies due to sensor error accumulation and initial value errors; moreover, it is only suitable for online applications [37,38].The second experimental method is the so-called Hybrid Pulse Power Characterization (HPPC) [39].This method consists of monitoring voltage variations during short-/high-current charge/discharge pulses with relaxation periods.In this method, the cell voltage response is affected by the electrochemical kinetics of the battery, such as ohmic losses, double-layer capacity, and lithium-ion diffusion [40,41].The third existing experimental method is Electrochemical Impedance Spectroscopy (EIS), which measures the current and voltage relationship over a frequency range and can be used to model kinetic battery processes through the impedance spectrum and estimate physical properties like diffusion coefficients and reaction rates.Impedance Spectroscopy (EIS) requires linearity, stability, and causality, which makes it challenging to implement as an online method [42].The experimental methods alone are more suitable for offline state-of-charge estimation due to the required experimental conditions.

Data-Driven Method
Data-driven approaches are model-free black-box methods that, for a given dataset, can provide rapid and accurate results given enough training data.These applications are the fruit of outstanding progress in machine learning in the recent past, which has motivated a renewed interest in applying data-driven algorithms to the field of the battery state-ofcharge estimation.Indeed, data-driven methods take temperature as an input characteristic, which avoids the establishment of a new thermal coupling model.Nevertheless, in the studies of ref. [43], the authors utilized a BP neural network to estimate the state of charge of lithium-ion batteries.In this regard, the estimation error of the state of charge was around 5%.On the other hand, the work performed in ref. [44] proposed another method to estimate the battery state of charge using Support Vector Machine (SVM).In this case, the authors established an average error of 3%.Furthermore, ref. [45] proposed the use of LS-SVM to correct the cumulative error of the ampere-hour integration method, which favors online state-of-charge estimation.However, these approaches do not take into consideration the issue of the number of training data.Thus, the work of ref. [46] suggests that using an LSTM method to estimate the state of charge can reduce the demand for training data for modeling to some extent.However, in modeling state-of-charge estimation under different conditions, data-driven methods, including RNNs, do not have great advantages.Generating sufficient training data is important but difficult because battery cell chemistry can change at any time, making existing training data non-representative.

Model-Based Methods
Model-based methods attempt to model battery behavior by incorporating physical and chemical principles into mathematical equations to accurately estimate the state of charge [47].The authors in ref. [2] reported that state-of-charge models that define battery capacity in units of energy (kWh) are classified as energy reservoir models (ERMs).Those that define it in units of charge (Ah) are classified as charge reservoir models (CRMs), and those that define it in units of concentration (mol/L) are classified as concentrationbased models.
The model-based approach involves the use of electrochemical, electrical, or a combination of both models to describe the battery dynamics from which the state of charge is estimated.The electrical model employs an approximate equivalent circuit to represent battery dynamics, and it incorporates adaptive filter algorithms like the Kalman Filter (KF), as referenced in [31].In this approach, constructing the circuit models involves arranging various circuit elements in series, parallel, or combinations thereof to mimic the battery's dynamic behavior.Several equivalent circuit models have been proposed, including the Rint model, the RC model, and the Thevenin model, as discussed in [1,2].These models include a critical component, a voltage source, which signifies the Open Circuit Voltage (OCV) of the battery.The OCV, in turn, depends on the battery SOC.It is essential to note that the accuracy of SOC estimation is heavily impacted by the specific equivalent circuit model chosen to represent the battery.For SOC estimation within the electrical model, the Extended Kalman Filter (EKF) is utilized.The EKF is an enhanced version of the Kalman Filter (KF) designed for estimating the internal state of non-linear dynamic systems, using a state-space model, as described in [48].Essentially, it predicts the future state of the system based on prior data.The EKF encompasses two key equations.One equation involves matrices constructed by using equivalent circuit model (ECM) parameters, along with the system SOC, measurable input matrices, and non-measurable process noise.The parameters are identified through standard testing procedures.The second equation is the measurement equation, which relates the output voltage to the system state vectors, measurable input matrices, and measurement noise.In the context of SOC estimation, the EKF leverages advanced battery cell models and demands relatively high computational resources.Furthermore, the electrical model's performance relies on capacitance and resistance values that are influenced by factors like state of charge, current, or temperature.This necessitates conducting numerous experiments to establish the relationships among these parameters, as mentioned in [49].To enhance the accuracy of the battery equivalent circuit model, researchers, as highlighted in the papers [50], often consider temperature effects and develop thermal coupling models.However, this approach increases the complexity of the model and introduces errors in the state-of-charge estimation.
The electrochemical model uses partial differential equations (PDEs) to describe battery dynamics [51].Indeed, one of the widely used electrochemical models is the pseudo-2-D (P2D) model (refs.[52,53]), which is a set of six partial differential equations, coupled with algebraic constraints.The P2D model takes current as its input and returns voltage as its output.The partial differential equations used in the P2D model incorporate the following parameters [53]: • Electrolyte lithium-ion diffusion equations in the positive electrode, negative electrode, and separator according to Fick's second law.• Solid-phase lithium-ion diffusion equations in the electrodes due to Fick's second law.

•
Electrolyte ohm equations in the electrodes and separator.• Solid-phase ohm equations in the positive electrode and negative electrode.

•
Charge conservation equations, including a positive electrode and a negative electrode, and a separator.• Butler-Volmer (BV) kinetic equations at the surface of the particles in the electrodes.
However, the P2D model suffers from a notable drawback stemming from its intricacy, primarily because it incorporates a multitude of partial differential equations (PDEs) and boundary conditions.This heightened complexity significantly escalates the computational workload involved.Furthermore, the substantial computational demands associated with the application of state and parameter estimation algorithms to the PDAEs (partial differential algebraic equations) forming the P2D model make it challenging to implement these electrochemical models in real-time applications (as referenced in [54,55]).Stateof-charge estimation from the full P2D model using a modified Extended Kalman Filter (EKF) has been reported in [56], where the electrochemical model is solved by the Chebyshev orthogonal collocation method to reduce the computational burden of the battery P2D electrochemical model.However, parameter estimation of the model remains a significant challenge [18].To address this challenge, previous efforts have been made to estimate the state of charge of lithium-ion batteries by using simplified PDE-based models, such as the Single-Particle Model (SPM) and the Single-Particle Model with Electrolyte Dynamics (SPME) (as mentioned in [57][58][59]).These endeavors aimed to reduce the computational complexity associated with the P2D electrochemical model's extensive set of PDEs.An important drawback of these simplified models is their limited suitability for scenarios involving higher C-rates.It has been observed, as highlighted in references [60], that these models begin to lose their accuracy when C-rates surpass 0.5 C. Additionally, it is worth noting that the Single-Particle Model with Electrolyte Dynamics (SPME) captures less cell dynamics for its performance assessment.However, the set of PDEs that form a P2D model requires a lot of development and computation time in online battery state-of-charge estimation [18].In this respect, it is crucial to develop a new mathematical model that is suitable for real-time implementation of lithium-ion battery state of charge.Recently, the critical review of three mathematical models, namely, the energy reservoir model, the reservoir of charge model, and the concentration-based model for determining the evolution of the lithium-ion battery state of charge was provided to the power systems research community by Rosewater et al. in their work [2], where they used these models to calculate the optimal timing of a battery energy storage system for a peak-shaving application.The authors highlighted the advantages and disadvantages of each model from a computational point of view but mainly looked at references outside the typical applications of grid-level battery energy storage systems.However, their models do not directly incorporate the initial state of charge into the dynamic state-of-charge equation of the battery, nor the nominal and terminal voltages of the battery.
Although, as we have already stressed above, batteries are essential in modern applications because of their power and energy density, incorrect charging and discharging can impair their performance and cause breakdowns.A battery management system (BMS) is, therefore, essential for their smooth operation.SOC estimation is a key element of any battery management system and is essential to maintaining the health and efficiency of the battery.In this respect, the contribution of this paper is to provide a new mathematical model to simulate the state of charge of a lithium-ion battery, taking into account the efficiency, charge transfer coefficient, nominal voltage, and terminal voltage during the charging and discharging processes.To do so, we consider certain approximations in each modeling stage and discuss their physical meaning and characteristic units.Our proposed model assumes a one-dimensional anisotropic energy reservoir box in the 12-volt battery system, taking into account various parameters, such as battery length, energy reservoir cross-sectional area, nominal voltage, and battery voltage, during the charging and discharging processes.We numerically solve our models by using a fourth-order Runge Kutta algorithm implemented in Lunix and generate a global state-of-charge map for different sets of characteristic values of the model parameters.We assume charging and discharging efficiencies of 80% and 95%, respectively.
The remainder of this paper is organized as follows: In Section 3, an energy reservoir model configuration of a rechargeable lithium-ion battery is presented.In Section 4, we present a detailed exposition of our models.Section 5 presents a case study and discusses the results obtained from the numerical implementation of these models.In Section 6, we summarize and highlight the importance of meticulously choosing operational parameters to balance fast charging and discharging time efficiency.

Model Configuration of a Rechargeable Lithium-Ion Battery
This work considers a rechargeable lithium-ion battery to be a power-energy reservoir model, and the detailed structure is illustrated in Figure 1.In this simplified model configuration, x ∈ R is the one-dimensional spatial coordinate, L ∈ R + is the total length of the energy storage reservoir (i.e., yellow box), H ∈ R + is the cross-sectional area of the energy reservoir, E ch is the inlet energy within the battery, and E dis is the outlet energy from the battery system.

Mathematical Modeling of Rechargeable Lithium-Ion Battery State of Charge
Mathematical modeling plays an important role in gaining in-depth knowledge of the dynamic behavior of rechargeable lithium-ion batteries, ultimately leading to improvements in their operational efficiency.By using the above power-energy reservoir model presented in Figure 1, the energy balance equation for the lithium-ion battery during charging/discharging can be modeled by Equation (2), which includes the entire process in the spatial domain: ∀x ∈ [0, L].
HL β ch P ch + P dis To be more explicit, Equation ( 2) can help to capture charging or discharging scenarios, depending on the specific conditions applied, which we will provide later.In this framework, key parameters such as P ch and P dis quantify the rate of transfer of electrical energy into and out of the battery, while β ch and β dis characterize the efficiency of charging and discharging, respectively.The voltage U b is the terminal voltage of the battery when it reaches its fully charging or discharging state, while q designates, under specific conditions, the state of charge added during charging or the state of charge removed during discharging so that (κdq/dt) indicates the contribution of the flow of energy stored or supplied per unit time, where κ signifies the energy capacity of the energy reservoir box.Furthermore, C 0 represents the standard capacity, γ is the constant parameter that governs the energy transfer, and µ stands for the initial state of charge when the battery is neither charging nor discharging.In the term (µ ± q), the plus sign (+) means that we are in a charging condition, while the minus sign (−) means that we are in a discharging condition.The reason why P ch is multiplied by β ch is that when a battery is charged, not all the energy supplied is stored; some is lost in the form of heat due to internal resistance and other inefficiencies.The charging efficiency β ch reflects the fraction of the energy supplied that is stored in the battery.If one supplies P ch to the battery, the amount of energy stored in the battery is This means that only a part of the input power is effectively converted into stored energy.The parameter P dis is divided by β dis because when discharging a battery, not all the stored energy can be converted back into useful output power; again, some of it is lost as heat due to inefficiencies.The discharging efficiency β dis indicates the fraction of the stored energy that can be converted into useful output power.If the battery discharges power at a rate of P dis , the amount of stored energy consumed from the battery is This formula considers that more energy needs to be drawn from the battery than is available as useful output power due to losses during the discharging process.The parameters and units in Equation ( 2) are summarized in Table 1.The charging and discharging power in Equation ( 2) can be expressed as the rates of change in energy with respect to time, as evidenced by Equations ( 5) and ( 6): The negative sign in Equation ( 6) indicates that discharging power (P dis ) is the rate at which energy (E dis ) is released from the battery over time.Substituting Equations ( 5) and ( 6) into Equation ( 2) yields To assess the efficiency of the charging and discharging processes, we correlate the rate at which energy is delivered to the battery during charging (see Equation ( 8)) and the rate at which energy is delivered from the battery during discharging (see Equation ( 9)), both as a function of energy.This correlation provides valuable insights into the overall performance and efficiency of the energy transfer processes involved in the battery's charging and discharging cycles.
where U ch and U dis are charging and discharging voltages, respectively.By replacing Equations ( 8) and ( 9) into Equation ( 7), we obtain where U ch is the voltage supplied to the battery during charging, and U dis is the voltage delivered by the battery to the load during discharging.As we have already stressed above, Equation ( 10) can, under specific conditions, describe either the dynamics of the state of charge of a lithium-ion battery during the charging or the discharging processes.
To satisfy this condition, we should have U dis = 0 during the charging process and U ch = 0 during the discharging process.Moreover, the form of U b should be well defined.To follow the ultimate goal of the present study, we conceived the battery terminal voltage as an empirical function of SOC(t).Therefore, depending on the battery operation conditions, U b will take the form where A and B represent the sensitivity coefficients of the variation in terminal voltage during the charging and discharging periods, respectively.

Charging Process Condition
Equation (10), to simulate the charging process, should take the following form: where q int is the inlet state of charge that should be added to increase the battery state-ofcharge level.The battery will reach the full charging state when To maintain the SOC in [%], we consider that 0.01 = 1% and the new expression of U ch b obtained from Equation ( 11) is introduced with Equation ( 14) into Equation ( 13) to provide Rearranging Equation (15) yields where , and In Formula ( 17), V n is the nominal voltage of the battery.Moreover, α 1 represents the charging rate of the battery.It indicates how quickly the battery is being charged.A higher value of α 1 means a faster charging rate, while a lower value means a slower charging rate.To solve Equation ( 16), the initial condition which represents the initial state of the battery should be given.In fact, at t = 0, when the charging process has not yet started, the battery is in the following state: However, when the charging process has not yet started, we assume that q int (0) = 0 and µ 0 = 20%.It follows from this that the initial condition of Equation ( 16) is

Discharging Process Condition
To simulate the discharging process, Equation (10) will now take the form where q out denotes that the battery state of charge is diminishing (i.e., q out < µ 2 ).Furthermore, the battery will reach the discharging state when Substituting Equation (20) into Equation ( 19) yields As µ 2 is a constant, Equation ( 21) can be written in the form It follows from this that where , and In Equation ( 24), α 2 represents the rate of discharging of the battery.It indicates the rate at which the battery is discharged.In addition, a higher value of α 2 means a faster rate of discharging, while a lower value means a slower rate of discharging.
As in the case of the equation that simulates the charging process, we also need the initial condition to solve Equation (23).In fact, when the battery has not yet started to discharge, we have However, as in the case of the charging condition when the discharging of the battery has not yet started, we assume that q out (0) = 0 and µ 100 = 100%.It follows from this that the initial condition to solve Equation (23) is where the subscript 100 means that the battery starts its discharging process from the state of charge SOC(0) = 100.

Comparison of Developed Model with Old Model
The work of [2] has provided a critical review of battery energy storage systems to the power systems research community.In their work, the authors present different types of energy reservoir models and highlighted the advantages and disadvantages of each model from a computational point of view.However, their models do not directly incorporate the initial state of charge or the nominal and terminal voltages of the battery into the dynamic state-of-charge equation of the battery which we have taken into account in our model.The simplicity of the model derived in the present study lies in the fact that the battery can be considered an energy reservoir in which energy is stored and then withdrawn for consumption in the electricity distribution system.In addition, the charging and discharging dynamic models we have provided in our Equations ( 16) and ( 23) are advantageous because they do not involve complex mathematical terms and there is no need to distinguish the elementary electrochemical units or the type of electrochemistry inside the battery as presented in [54].Another advantage of the models we propose is that the initial state of charge of the battery is directly integrated into the mathematical model, and so are the nominal and terminal voltages of the battery.The control variables of these models are the terminal charging voltage U ch b and the terminal discharging voltage U dis b , while the state of charge is the only state variable.It is important to note that our proposed model accounts for the presence of losses during the charging and discharging cycles, which are quantified by the separate energy efficiency factors, β ch and β dis , assigned to each operation.In particular, the model formulated in this study provides a new theoretical framework for understanding the complex interaction of battery terminal voltage in the dynamic equations governing the state of charge.Some important parameters of the model are listed in Table 2.

Case Study
As a case study, the model developed in the present work is applied to a 12-volts lithium-ion battery.In fact, for this type of battery, the terminal voltage varies during charging and discharging, and its sensitivity coefficient reflects how it changes in response to current flow.The charging terminal voltage increases from 13.2 V to 14.6 V as energy is stored, while the discharging terminal voltage decreases from 12.9 V to 12 V as energy is released.Furthermore, the energy transfer coefficient characterizes the efficiency of energy conversion from electrical to chemical forms during charging or from chemical to electrical forms during discharging processes.All these parameter variations influence the battery SOC.Therefore, we have chose to vary these parameters in the finite range in the case of charging and discharging, and the results are summarized in Tables 3 and 4. Figure 2 shows the evolution of the state of charge, SOC(t), of the lithium-ion battery as a function of time for different values of the rate of charging α 1 and for V n = 12.8 V, U ch = 14.6 V, β ch = 80%, and µ 0 = 20%.Each curve corresponds to α 1 = 0.365 (blue curve), α 1 = 0.57 (orange curve), and α 1 = 1.2 (green curve).In Figure 2, we observe that at higher terminal charging voltages U ch b , voltage-sensitive coefficients A and energy transfer coefficients γ 1 lead to faster increases in time, making them advantageous for applications requiring fast charging.However, it should be noted that large values of these parameters can lead to drawbacks, such as increased battery degradation, overheating, and wasted energy.In addition, practical consideration of charging time, which varies from 8 minutes at a charging rate of 1.2 for SOC = 100% to 16 minutes at rates of 0.365 and 0.57 for SOC = 70% and SOC = 99%, respectively, highlights the trade-off between fast charging and time efficiency.These results are summarized in Table 3. Figure 3 shows the discharging profile of the lithium-ion battery as a function of time for different values of discharging rate α 2 and for V n = 12.8 V, U dis = 12 V, β dis = 95%, and µ 100 = 100%.Indeed, Figure 3 suggests that varying the terminal discharging voltage U dis b allows the discharging process to be controlled.More explicitly, a particular feature of these results that emerges from the numerical simulations of Equation ( 23) is the observation that a higher terminal voltage (i.e., 12.9 V) results in a slower discharging rate α 2 = 1.15, which can be advantageous when a longer battery life is desired.This variation allows the SOC to be adjusted to the desired level (i.e., 20%).In addition, setting the sensitivity coefficient (B) affects how the voltage reacts to variations in SOC.A lower B value (see, for example, B = 3.395 × 10 −2 ) results in lower sensitivity, giving a more stable voltage during discharging.This can be advantageous for applications requiring a steady power supply.However, a high γ 2 leads to excessive discharging, as observed when α 2 = 2.85.This can significantly reduce battery life and damage the battery.In addition, a higher α 2 (for example, 2.85) leads to rapid discharging, resulting in over-discharging to 9% SOC after 4 h.These results are summarized in Table 4.

Conclusions
This study provided valuable insights into the dynamic behavior of lithium-ion batteries during both the charging and discharging processes.We observed that varying terminal voltages, sensitivity coefficients, and energy transfer coefficients significantly impact the SOC and discharging rates.Higher terminal charging voltages, as demonstrated by our results, can lead to faster SOC increases, making them advantageous for applications requiring rapid charging.However, it is crucial to strike a balance, as excessive values of these parameters can result in battery degradation, overheating, and energy wastage.On the other hand, in the discharging process, a higher terminal voltage leads to a slower discharging rate, which is beneficial for extending battery life.Adjusting sensitivity coefficients can provide a more stable voltage during discharging, essential for applications demanding a consistent power supply.Yet, a high discharging rate can significantly reduce battery life and lead to over-discharging.These findings underscore the necessity of carefully selecting the operating conditions to achieve a balance among fast charging, time efficiency, and battery longevity in practical applications.However, certain factors, such as internal resistance, temperature, charging/discharging behavior, and aging, which can have an impact on the SOC, were not taken into account in this study.Therefore, for greater accuracy of the model, the inclusion of these factors should be considered in the future.

Figure 1 .
Figure 1.Simplified structure of power-energy reservoir model of lithium-ion battery.

2 Figure 2 .
Figure 2. (Color online.)The charging profile of a lithium-ion battery as a function of time for different values of charging rate α 1 and for V n = 12.8 V, U ch = 14.6 V, β ch = 80%, and µ 0 = 20%.

Figure 3 .
Figure 3. (Color online.)The discharging profile of a lithium-ion battery as a function of time for different values of discharging rate α 2 for V n = 12.8 V, U dis = 12 V, β dis = 95%, and µ 100 = 100%.

Table 1 .
Parameters and units.

Table 3 .
Simulation parameters for the charging condition.

Table 4 .
Simulation parameters for the discharging condition.