Investigation of the Electrochemical and Thermal Characteristics of NCM811-21700 Cylindrical Lithium-Ion Battery: A Numerical Study and Model Validation

: Recently, there has been growing recognition of the signiﬁcance of energy and environmental challenges. Utilization of lithium-ion batteries in electric vehicles has shown considerable potential and beneﬁts for tackling these issues. The effective management of battery temperature has become a crucial factor in the advancement and widespread adoption of lithium-ion batteries in electric vehicles. In this study, a thermo-coupled pseudo-two-dimensional (P2D) electrochemical model is employed to simulate the heat generation of the NCM811-21700 cylindrical battery cell at various discharge rates at an ambient temperature of 25 ◦ C, and is validated by experimental data. The validation results demonstrate that the thermo-coupled P2D model can effectively predict the battery voltage curve during the discharge process with less than 4% errors. Although there is a slightly larger error in the temperature prediction during the battery 2C and 3C discharge processes, the maximum error approaches 10%, which is still generally within an acceptable range. In addition, the battery’s electrochemical and thermal characteristics during discharge are presented. The suggested thermo-coupled electrochemical model can be used for applications in the thermal management system of the NCM811-21700 battery.


Introduction
Currently, a significant portion (approximately 80%) of the world's energy provision hinges on carbonaceous resources, such as anthracite, petroleum, and methane. The preponderance of energy production and use is often associated with ecological degradation and deleterious impacts on H. sapiens' wellness. In the face of our global community's mounting dependence on renewable energy pathways, vehicles propelled by electricity, or EVs, represent an auspicious means of creating a more pristine and sustainable environment. This is largely due to the batteries' superior capabilities to preserve a constant supply of energy derived from these greener alternatives [1][2][3]. Batteries are mechanisms that employ an electrochemical interaction between two or more cells to supply energy to electric appliances. They encompass numerous prevalent aspects, such as the self-discharge rate, specific energy, operational parameters, toxicity, recharging period, and C-rate. One essential trait is battery capacity, which is shaped by the reversible conditions of electrochemical activity. This attribute represents the maximum volume of electricity that a battery can accumulate and dispense. This capacity decreases with the deployment of electrodes or electrolytes. The C-rate refers to the charging or discharging rate and is used to measure the discharging of a battery relative to its maximum capacity [4].
Research on Battery Thermal Management Systems (BTMS) can be broadly classified into two approaches: direct measurements and indirect calculations. Direct measurement the thermal behavior of a 1S18P 21700 battery pack using a numerical method. The cell temperature was predicted with an error of less than 0.5 • C. Sahin et al. [20] simulated a BTMS for a module with 12 21700 cells. An electrochemical-thermal model was used to predict the thermal characteristics and was validated by experiments. The results indicate that the novel cooling method can effectively reduce the maximum cell temperature and the temperature difference between cells. A hexagonal cell module was developed for various discharge rates. These designs incorporated a phase change material (PCM) composed of a graphene-paraffin composite with high thermal conductivity. The results indicated that the PCM structure effectively regulated the thermal parameters, ensuring that the battery temperature remained stable and consistent [21]. Sheng et al. [22] developed a cellular liquid cooling jacket system for a 21700 lithium-ion battery module. The findings indicate that employing interlaced flow directions leads to a reduced standard deviation in the temperature and a more uniform thermal distribution.
From the literature review, it can be observed that there needs to be more research on NCM811-21700 batteries than on the more commonly studied 18,650 batteries, particularly in the field of numerical modeling. The P2D model, which can effectively represent the electrochemical and thermal characteristics of batteries, is frequently utilized in numerical studies of 18,650, 26,650 [7,9,10,23], and other types of batteries [8,11,24,25]. However, numerical studies utilizing P2D models to simulate the thermal performance of NCM811-21700 batteries are limited. Therefore, in this study, a numerical analysis of the thermal characteristics of an NCM811-21700 lithium-ion battery cell was performed using a thermally coupled Newman P2D electrochemical model. The heat generated by the battery cell was evaluated at various discharge rates and ambient temperatures. The results are validated using experimental data. Additionally, this study analyzes the electrochemical characteristics of the battery, and the developed model will be extensively applied in research related to the BTMS of the NCM811-21700 battery.

Numerical Methods
An electrochemical-thermal model was established for a 5 Ah NCM811-21700 cylindrical battery cell using the commercial CFD software ANSYS FLUENT 2022R1 [26]. The electrochemical model was based on the physical-based Newman P2D model, which evaluated the ion concentration and potential distribution of the battery cell [5]. The heat generation, heat dissipation, and temperature distribution of the battery cell were calculated by coupling P2D and thermal models. The cell structure consisted of several layers, including current collectors, porous electrodes, and a separator along the cell thickness, which was used for the thermally coupled Newman P2D electrochemical model, and is shown in Figure 1. level on a 21700 battery cell. The model was parameterized using an experimental method and exhibited high modeling and estimation accuracy. Research on 21700 batteries related to the module or pack level has also been conducted [18]. Kang et al. [19] examined the thermal behavior of a 1S18P 21700 battery pack using a numerical method. The cell temperature was predicted with an error of less than 0.5 °C. Sahin et al. [20] simulated a BTMS for a module with 12 21700 cells. An electrochemical-thermal model was used to predict the thermal characteristics and was validated by experiments. The results indicate that the novel cooling method can effectively reduce the maximum cell temperature and the temperature difference between cells. A hexagonal cell module was developed for various discharge rates. These designs incorporated a phase change material (PCM) composed of a graphene-paraffin composite with high thermal conductivity. The results indicated that the PCM structure effectively regulated the thermal parameters, ensuring that the battery temperature remained stable and consistent [21]. Sheng et al. [22] developed a cellular liquid cooling jacket system for a 21700 lithium-ion battery module. The findings indicate that employing interlaced flow directions leads to a reduced standard deviation in the temperature and a more uniform thermal distribution. From the literature review, it can be observed that there needs to be more research on NCM811-21700 batteries than on the more commonly studied 18,650 batteries, particularly in the field of numerical modeling. The P2D model, which can effectively represent the electrochemical and thermal characteristics of batteries, is frequently utilized in numerical studies of 18,650, 26,650 [7,9,10,23], and other types of batteries [8,11,24,25]. However, numerical studies utilizing P2D models to simulate the thermal performance of NCM811-21700 batteries are limited. Therefore, in this study, a numerical analysis of the thermal characteristics of an NCM811-21700 lithium-ion battery cell was performed using a thermally coupled Newman P2D electrochemical model. The heat generated by the battery cell was evaluated at various discharge rates and ambient temperatures. The results are validated using experimental data. Additionally, this study analyzes the electrochemical characteristics of the battery, and the developed model will be extensively applied in research related to the BTMS of the NCM811-21700 battery.

Numerical Methods
An electrochemical-thermal model was established for a 5 Ah NCM811-21700 cylindrical battery cell using the commercial CFD software ANSYS FLUENT 2022R1 [26]. The electrochemical model was based on the physical-based Newman P2D model, which evaluated the ion concentration and potential distribution of the battery cell [5]. The heat generation, heat dissipation, and temperature distribution of the battery cell were calculated by coupling P2D and thermal models. The cell structure consisted of several layers, including current collectors, porous electrodes, and a separator along the cell thickness, which was used for the thermally coupled Newman P2D electrochemical model, and is shown in Figure 1.

Model Description
The Newman P2D model is a physics-based model developed by Newman et al. that uses porous electrodes and concentrated solution theories [28]. The model is termed "pseudo" because its second dimension follows the electrode's radius rather than being a true second dimension. This model functions similarly to the 1D model but also incorporates the impact of diffusion in the electrolyte and solid (electrode) phases and the Butler-Volmer kinetics for simulating electrode kinetics. In this study, a thermally coupled electrochemical model was solved using commercial CFD software ANSYS FLUENT 2022R1. The conservation equations for the electrolyte and solid phases, effective properties, and energy conservation equations for the thermal model in the FLUENT solver are presented in Table 1. Some assumptions for the P2D modeling are listed as follows [24]: 1. The local reaction current density is described by the Butler-Volmer equation; 2. In both the positive and negative electrodes, the active material particles have a spherical shape and are evenly distributed; 3. Li+ ions move within the electrolyte mainly due to diffusion and electro-migration processes; 4. Side reactions are ignored (j side in Equation (4) is equal to 0); 5. The transportation phenomenon of the electrolyte is governed by the theory of concentration solutions. Table 1. Governing equations of thermally coupled Newman P2D model [26,27].

Electrochemical Model Governing Equations
Lithium conservation in the solid phase: Lithium conservation in the electrolyte phase: Charge conservation in solid phase: Charge conservation in the electrolyte phase: The Butler-Volmer equation: The current flux at the active material pore-wall surface due to lithium intercalation reaction as a function of over-potential η: where η is defined by: and i 0 is the exchange current density defined by: Effective properties: σ e e f f = κε β e σ )

Thermal model Governing equations
The energy conservation equation q = q re + q ohm + q act + q ex (16) q re = j int T ∂U ∂T (17)

Model Parameters
The dependability of the Newman P2D model is critically contingent on the precision of its constituent parameters. Each is meticulously fine-tuned to the specific design and chemistry of the individual cells. Consequently, the extrapolation of parameters from a single-cell design to a universal application is untenable. This underlines the primary challenge in battery modeling, which is to identify an optimal set of parameters tailored for a specific battery type. Certain researchers might endeavor to calibrate the model for congruence with empirical data, as suggested in [29]. However, given the intricate nature of the model and the sheer magnitude of its parameters, such a method frequently proves unfeasible without a robust initial approximation of these parameters. Therefore, the use of parameters sourced from the literature is another option. In this study, the physical and electrochemical parameters are obtained from identifying previous studies [30][31][32] and are listed in Table 2.
The effective thermal conductivity was calculated using the series-parallel connection method, which is similar to the calculation of series-parallel resistance. The effective thermal conductivity in the radial (u) and axial (z) directions is given by: An expression for the OCP curve was required in the simulation and was fitted as a function of stoichiometry for both the positive and negative electrodes from the experimental dataset [30]. In this study, the correlations in the literature are represented as polynomials in the stoichiometry by Equations (25) and (26) and are shown in Figure 2.   The electrolyte diffusivity for the electrode and conductivity of the electrolyte are concentration-dependent and are expressed by Equations (27) and (28) [30], as shown in Figure 3.  Equations (15)- (20) show that several mechanisms contribute to heat generation within battery systems. Notably, the irreversible and reversible heat components are perceived as crucial and primarily attributable to heat generation from active materials. In scenarios with high C-rates, more than half of the heat generation typically originates from irreversible heat (qohm), which is frequently referred to as ohmic heat loss. Conversely, at lower C-rates (1C or less), the significance of reversible heat generation (qre) was amplified. This form of heat generation is a consequence of entropy variations induced by intercalation reactions and is contingent upon the system's internal temperature and open-circuit voltage (OCV). At various lithiation stages, the positive electrode's entropic term was practically insignificant. The remaining lithiation stages did not demonstrate substantial value. This implies that the behavior of the whole cell is chiefly influenced by the negative Considering the cell composition, which includes numerous distinct thin layers, maintaining the same thermal model implies that each layer corresponds to a set of partial differential equations (PDEs). However, this approach, which is computationally demanding, offers only a marginal improvement in accuracy. Therefore, to simplify the model, the thermal and physical parameters of the active material were assumed to be constant and uniform. Then, the effective density and effective specific heat capacity can be evaluated [32]:  (25) and Figure 2 Equation (26)  The effective thermal conductivity was calculated using the series-parallel connection method, which is similar to the calculation of series-parallel resistance. The effective thermal conductivity in the radial (u) and axial (z) directions is given by: An expression for the OCP curve was required in the simulation and was fitted as a function of stoichiometry for both the positive and negative electrodes from the experimental dataset [30]. In this study, the correlations in the literature are represented as polynomials in the stoichiometry by Equations (25) and (26) and are shown in Figure 2. The electrolyte diffusivity for the electrode and conductivity of the electrolyte are concentration-dependent and are expressed by Equations (27) and (28) [30], as shown in Figure 3.
Equations (15)- (20) show that several mechanisms contribute to heat generation within battery systems. Notably, the irreversible and reversible heat components are perceived as crucial and primarily attributable to heat generation from active materials. In scenarios with high C-rates, more than half of the heat generation typically originates from irreversible heat (q ohm ), which is frequently referred to as ohmic heat loss. Conversely, at lower Crates (1C or less), the significance of reversible heat generation (q re ) was amplified. This form of heat generation is a consequence of entropy variations induced by intercalation reactions and is contingent upon the system's internal temperature and open-circuit voltage (OCV). At various lithiation stages, the positive electrode's entropic term was practically insignificant. The remaining lithiation stages did not demonstrate substantial value. This implies that the behavior of the whole cell is chiefly influenced by the negative electrode, following the same pattern but with an inverse magnitude. Heat generation in a cell that can be reversed is governed primarily by the negative electrode [31]. Therefore, in this study, the entropic term ∂U/∂T in Equation (17) is fixed as a function of State of Charge (SOC) (Figure 4) and can be expressed as [31]: Energies 2023, 16, x FOR PEER REVIEW 8 of 1 Figure 4. The fitted entropic term ∂U/∂T as a function of SOC.

Mesh, Boundary Conditions, and Initial Conditions
A sketch of the 5 Ah NCM811-21700 cylindrical cell is shown in Figure 5. The ce structure was simplified to an active zone, a positive tab, a negative tab, and a void. Thi study used the commercial computational fluid dynamics (CFD) software ANSYS FLU ENT to solve a thermally coupled P2D electrochemical model. In total, 317,372 tetrahedra

Mesh, Boundary Conditions, and Initial Conditions
A sketch of the 5 Ah NCM811-21700 cylindrical cell is shown in Figure 5. The cell structure was simplified to an active zone, a positive tab, a negative tab, and a void. This study used the commercial computational fluid dynamics (CFD) software ANSYS FLUENT to solve a thermally coupled P2D electrochemical model. In total, 317,372 tetrahedral and hexahedral elements were generated, as shown in Figure 5. Because most governing equations are solved in the solid domain, the mesh quality does not significantly impact the computational results; thus, a mesh-independent test is not required for this study. The shell of the battery was set to be 0.34 mm [31] with a constant convective heat transfer coefficient (15 W/m 2 K). The density, thermal conductivity, and specific heat capacity of the shell material are 7900 kg/m 3 , 15 W/mK, and 477 J/kgK, respectively [31]. The internal temperature of the cell was set as equal to the ambient temperature at initial time. The discharge process at various C-rates was simulated using the transient solution method with a time step size of 1s. Additionally, a User Define Function (UDF) was used to define the entropic term ∂U/∂T in the FLUENT environment.

Mesh, Boundary Conditions, and Initial Conditions
A sketch of the 5 Ah NCM811-21700 cylindrical cell is shown in Figure 5. The cell structure was simplified to an active zone, a positive tab, a negative tab, and a void. This study used the commercial computational fluid dynamics (CFD) software ANSYS FLU-ENT to solve a thermally coupled P2D electrochemical model. In total, 317,372 tetrahedral and hexahedral elements were generated, as shown in Figure 5. Because most governing equations are solved in the solid domain, the mesh quality does not significantly impact the computational results; thus, a mesh-independent test is not required for this study. The shell of the battery was set to be 0.34 mm [31] with a constant convective heat transfer coefficient (15 W/m 2 K). The density, thermal conductivity, and specific heat capacity of the shell material are 7900 kg/m 3 , 15 W/mK, and 477 J/kgK, respectively [31]. The internal temperature of the cell was set as equal to the ambient temperature at initial time. The discharge process at various C-rates was simulated using the transient solution method with a time step size of 1s. Additionally, a User Define Function (UDF) was used to define the entropic term ∂U/∂T in the FLUENT environment.

Results and Discussion
Numerical simulations were conducted on the 5 Ah NCM811-21700 cylindrical cell at discharge rates of 1C, 2C, and 3C under a temperature of 25 • C. This study aims to comprehend the cell thermal attributes under varying discharge rates during discharge and validate the thermally coupled Newman P2D model for thermal and electrochemical performance prediction. The outcomes are derived from the cell depth of discharge (DoD), which signifies the extent of battery discharge, where "1" indicates a fully charged state, and "0" represents a fully discharged state. The DoD of the discharge process can be evaluated as:

Model Validation
The developed thermally coupled electrochemical P2D model was validated using data obtained from an experimental study by Catenaro et al. [33]. Initially, the model was validated using experimental data of battery voltage versus time and DoD during the discharge process with the different discharge rates of 1C, 2C, and 3C at a steady ambient temperature of 25 • C. As shown in Figure 6, the simulation results of the cell voltage are consistent with the experimental data for every discharge rate. An error analysis was also performed. Given the discrepancy in the number of simulation results and experimental data, the error analysis was conducted according to the following steps: First, the simulation results were fitted into a polynomial about the DoD. The DoD values from the experimental data were substituted into this polynomial to obtain the corresponding fitted values. Finally, the fitted values were used for error calculations against the experimental data. Therefore, the error between simulation result and experimental can be defined as: where θ fitted and θ exp indicate the fitted data from the simulation result and the experiment data, respectively. The error in the cell voltage is shown in Figure 7. The maximum errors for 1C, 2C, and 3C were 2.8%, 2.7%, and 3.6%, respectively. Larger errors usually appear at the beginning or end of battery discharge. This may be due to unstable measurements of the battery at the beginning or end of discharge. The P2D model accurately reflected the voltage change during the battery discharge process. The thermal performance is also validated by the comparison of the cell surface temperature profile between simulation results and experimental data versus time and DoD during the discharge process with the different discharge rates of 1C, 2C, and 3C at a steady ambient temperature of 25 • C. As shown in Figures 8 and 9, the error in the temperature profile was relatively large, especially under the conditions of 2C and 3C. The same method was used to analyze the error in the temperature profile. The maximum errors at 1C, 2C, and 3C were 2.4%, 9.7%, and 9.8%, respectively. The temperature profile under the 1C discharge process was significantly smaller than those of 2C and 3C. This is because, in the P2D model, there are many parameters based on temperature variables, such as the specific heat of the battery and thermal conductivity. However, this study did not consider the effect of temperature on these parameters because of the lack of reference data. However, the measurement of battery temperature also has significant instability. This is affected by the measurement location, measurement accuracy, and external heat transfer conditions. Even in a constant temperature chamber, the heat transfer coefficient on the battery surface is uncertain. Therefore, when the battery was discharged at 2C and 3C, the temperature was much higher than the reference temperature, which led to large errors in the simulation results. In addition, the temperature error curve was not as smooth as the voltage error curve. This was likely because the temperature measurements were less stable than the voltage measurements during the experimental process. However, compared with other numerical analysis studies [8,13,14] on batteries of the same type, this error was not significant. Therefore, it can be said that the model's reliability has been verified under the ambient temperature of 25 • C.

Electrochemical and Thermal Analysis
The distribution of the current density inside the NCM811-21700 battery under the 1C, 2C, and 3C discharge processes is presented in Figure 10. The arrows are colored according to the magnitude of the current density. As shown in Figure 10, the direction of the arrows indicates that the electrical current is input via the negative current collector tab (bottom) and then output through the positive current collector tab (top). The magnitude of the current density correlates well with the discharge rate. In lithium-ion batteries, electrolyte solutions play the role of ionic conductors and unavoidably create concentration gradients during the operation of the battery. Recent studies [34] reveal that beyond these gradients, unique lithium-ion (Li+) concentration waves are generated in the electrolyte of graphite-based porous electrode/lithium cells. This particular phenomenon was also scrutinized in this simulation's approaches and is presented in Figure 11. The formation of these waves stems from the irregularities in the reaction distribution within the porous electrode, and they are influenced by both thermodynamic factors, such as open-circuit voltage (OCV), and kinetic factors, such as the unevenness of the charge transfer reaction. Figure 12 presents the distribution of the NCM811 21700 battery cell voltage of the 1C discharge rate at time intervals of 10, 30, and 50 min. During the discharging process, Energies 2023, 16, 6407 10 of 16 the cell voltage was not uniform across the entire cell. However, there was a slight drop in the overall battery working voltage from 2 to 50 min of discharge time, indicating a stable discharge voltage of the NCM811 battery. A significant potential gradient is observed near the positive end (top) of the cell. This can be attributed to the composition of the positive current collector and tab made of Al foil, in contrast to the negative current collector and tab (bottom tab) made of Cu foil, which has a higher electrical conductivity than Al foil. Moreover, the unevenness of the local cell voltage can be attributed to the heterogeneous electrochemical reaction rate in the battery. The local cell voltage was determined from the difference in potential between the positive and negative current collectors. The temperature distributions inside the battery cell under a 1C discharge rate at 10, 30, and 50 min are presented in Figure 13. The temperature gradient in the battery from the center to the outside surface is presented because of heat loss to the constanttemperature boundary condition. Although the temperature difference between the battery surface and the center is small, it can still be observed that this temperature difference increases with increasing discharge time.    The thermal performance is also validated by the comparison of the cell surface temperature profile between simulation results and experimental data versus time and DoD during the discharge process with the different discharge rates of 1C, 2C, and 3C at a steady ambient temperature of 25 °C. As shown in Figures 8 and 9, the error in the temperature profile was relatively large, especially under the conditions of 2C and 3C. The same method was used to analyze the error in the temperature profile. The maximum errors at 1C, 2C, and 3C were 2.4%, 9.7%, and 9.8%, respectively. The temperature profile under the 1C discharge process was significantly smaller than those of 2C and 3C. This is because, in the P2D model, there are many parameters based on temperature variables, such as the specific heat of the battery and thermal conductivity. However, this study did not consider the effect of temperature on these parameters because of the lack of reference data. However, the measurement of battery temperature also has significant instability. This is affected by the measurement location, measurement accuracy, and external heat transfer conditions. Even in a constant temperature chamber, the heat transfer coefficient on the battery surface is uncertain. Therefore, when the battery was discharged at 2C and 3C, the temperature was much higher than the reference temperature, which led to large errors in the simulation results. In addition, the temperature error curve was not as smooth as the voltage error curve. This was likely because the temperature measurements were less stable than the voltage measurements during the experimental process. However, compared with other numerical analysis studies [8,13,14] on batteries of the same type, this error was not significant. Therefore, it can be said that the model's reliability has been verified under the ambient temperature of 25 °C.

Electrochemical and Thermal Analysis
The distribution of the current density inside the NCM811-21700 battery under the 1C, 2C, and 3C discharge processes is presented in Figure 10. The arrows are colored according to the magnitude of the current density. As shown in Figure 10, the direction of the arrows indicates that the electrical current is input via the negative current collector tab (bottom) and then output through the positive current collector tab (top). The magnitude of the current density correlates well with the discharge rate. In lithium-ion batteries, electrolyte solutions play the role of ionic conductors and unavoidably create concentration gradients during the operation of the battery. Recent studies [34] reveal that beyond these gradients, unique lithium-ion (Li+) concentration waves are generated in the electrolyte of graphite-based porous electrode/lithium cells. This particular phenomenon was also scrutinized in this simulation's approaches and is presented in Figure 11. The formation of these waves stems from the irregularities in the reaction distribution within the porous electrode, and they are influenced by both thermodynamic factors, such as opencircuit voltage (OCV), and kinetic factors, such as the unevenness of the charge transfer reaction.

1C
2C 3C Figure 10. Current density inside of the cell at various discharge rates.

Electrochemical and Thermal Analysis
The distribution of the current density inside the NCM811-21700 battery under the 1C, 2C, and 3C discharge processes is presented in Figure 10. The arrows are colored according to the magnitude of the current density. As shown in Figure 10, the direction of the arrows indicates that the electrical current is input via the negative current collector tab (bottom) and then output through the positive current collector tab (top). The magnitude of the current density correlates well with the discharge rate. In lithium-ion batteries, electrolyte solutions play the role of ionic conductors and unavoidably create concentration gradients during the operation of the battery. Recent studies [34] reveal that beyond these gradients, unique lithium-ion (Li+) concentration waves are generated in the electrolyte of graphite-based porous electrode/lithium cells. This particular phenomenon was also scrutinized in this simulation's approaches and is presented in Figure 11. The formation of these waves stems from the irregularities in the reaction distribution within the porous electrode, and they are influenced by both thermodynamic factors, such as opencircuit voltage (OCV), and kinetic factors, such as the unevenness of the charge transfer reaction.

1C
2C 3C Figure 10. Current density inside of the cell at various discharge rates.   termined from the difference in potential between the positive and negative current collectors. The temperature distributions inside the battery cell under a 1C discharge rate at 10, 30, and 50 min are presented in Figure 13. The temperature gradient in the battery from the center to the outside surface is presented because of heat loss to the constant-temperature boundary condition. Although the temperature difference between the battery surface and the center is small, it can still be observed that this temperature difference increases with increasing discharge time.

Conclusions
This study employed a thermally coupled electrochemical P2D model to represent the electrochemical and thermal characteristics of an NCM811-21700 cylindrical lithiumion battery. The model parameters were obtained from other studies owing to the challenges and limitations of parameter identification. The commercial CFD software ANSYS FLUENT 2022R1 was used to solve the thermal and electrochemical submodels. Simulations were conducted to represent the discharge process under the discharge rates of 1C, 2C, and 3C, with a constant ambient temperature of 25 °C. The voltage curves and temperature profiles versus the discharge time and DoD were validated using experimental data. The validation results indicated that the model represented voltage curves with errors of less than 4% for every discharge rate. The maximum errors in the cell voltage at 1C, 2C, and 3C were 2.7%, 2.8%, and 3.6%, respectively. The errors of the temperature profiles were minor, and the maximum error at 1C, 2C, and 3C were 2.4%, 9.7%, and 9.8%, respectively. However, this error is still acceptable, considering the uncertainty in the selection of model parameters and the instability in measuring the battery temperature. Using the validated model, the electrochemical and thermal characteristics of an NCM811 battery were obtained from numerical simulations. The vector diagram of the current density, variation of lithium-ion concentration in the electrolyte phase, and the distribution of the internal voltage and temperature in the battery during the discharge process are discussed as being reasonable. This also indicates that the model developed in this study can be effectively applied to the thermal management research of the NCM811-21700 cylindrical battery.
A limitation of this study is that the effect of ambient temperature was not investigated owing to the lack of investigation of temperature-dependent parameters. Future experimental work is expected to involve performance testing, parameter identification, and numerical model validation for the NCM811-21700 battery under different ambient temperature conditions. Furthermore, the presented P2D model is only applied to inves-

Conclusions
This study employed a thermally coupled electrochemical P2D model to represent the electrochemical and thermal characteristics of an NCM811-21700 cylindrical lithium-ion battery. The model parameters were obtained from other studies owing to the challenges and limitations of parameter identification. The commercial CFD software ANSYS FLUENT 2022R1 was used to solve the thermal and electrochemical submodels. Simulations were conducted to represent the discharge process under the discharge rates of 1C, 2C, and 3C, with a constant ambient temperature of 25 • C. The voltage curves and temperature profiles versus the discharge time and DoD were validated using experimental data. The validation results indicated that the model represented voltage curves with errors of less than 4% for every discharge rate. The maximum errors in the cell voltage at 1C, 2C, and 3C were 2.7%, 2.8%, and 3.6%, respectively. The errors of the temperature profiles were minor, and the maximum error at 1C, 2C, and 3C were 2.4%, 9.7%, and 9.8%, respectively. However, this error is still acceptable, considering the uncertainty in the selection of model parameters and the instability in measuring the battery temperature. Using the validated model, the electrochemical and thermal characteristics of an NCM811 battery were obtained from numerical simulations. The vector diagram of the current density, variation of lithium-ion concentration in the electrolyte phase, and the distribution of the internal voltage and temperature in the battery during the discharge process are discussed as being reasonable. This also indicates that the model developed in this study can be effectively applied to the thermal management research of the NCM811-21700 cylindrical battery.
A limitation of this study is that the effect of ambient temperature was not investigated owing to the lack of investigation of temperature-dependent parameters. Future experimental work is expected to involve performance testing, parameter identification, and numerical model validation for the NCM811-21700 battery under different ambient temperature conditions. Furthermore, the presented P2D model is only applied to investigate the thermal and electrochemical performance during the discharge process at the cell level. Models for the charging process as well as at the module and pack levels still need to be developed and validated in future work.  Data Availability Statement: Data is unavailable due to privacy restrictions.