SEARCH · Engineering Papers

Results for “Source and propagation modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Assessing Photovoltaic Capacity Factor Variability Using Long-Term Satellite Derived Solar Resource Data Under Brazilian Climate

Accurate estimation of photovoltaic (PV) energy yield and its variability is essential for reducing financial risk and supporting reliable system planning for rapidly expanding PV markets. In Brazil, high solar adoption and increasing levels of distributed energy resources are beginning to introduce operational challenges such as curtailment and evolving grid requirements. Understanding how natural variability in solar resource propagates into PV system performance is therefore increasingly important for both project design and grid integration. Modern PV yield assessments commonly rely on multi-year meteorological datasets and probabilistic exceedance metrics (e.g., P50/P90) to quantify energy yield uncertainty for project financing. However, the implications of long-term solar resource variability for PV system design choices and high-adoption grid conditions remain less well characterized for rapidly expanding markets such as Brazil. In particular, understanding how weather-driven variability propagates into PV production distributions and capacity factor expectations is important for evaluating curtailment exposure, deployment strategies, and storage requirements in regions experiencing rapid growth of distributed and utility-scale PV. Seasonal and interannual variability in atmospheric conditions can produce substantial fluctuations in monthly PV energy production, which propagate into uncertainty in annual energy yield and capacity factor expectations. Characterizing this variability using long-term meteorological datasets allows probabilistic estimation of PV system performance and provides improved insight into the range of expected PV energy outcomes. This study explores the use of long-term satellite-derived meteorological data from the National Solar Radiation Database (NSRDB) to evaluate the variability of photovoltaic system performance across multiple locations in Brazil. Using a 27-year dataset (1998-2024), PV system simulations are performed to characterize the distribution of annual and seasonal capacity factors and energy yield outcomes, while propagating key sources of meteorological variability and model uncertainty through the PV modeling chain. The analysis also investigates the sensitivity of PV performance outcomes to key system design assumptions within the PV modeling chain, including tracking configuration and system sizing parameters. The resulting probabilistic performance characterization provides insight into how weather-driven variability influences PV production expectations and capacity factor distributions. These results provide a foundation for evaluating how weather-driven variability interacts with high PV adoption and potential storage or curtailment mitigation strategies.

14 SOLAR ENERGY↗

DOE OSTI2024

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Collaborative Research: Understanding air‐sea feedbacks to the MJO through process evaluation of observations and E3SM experiments (Final Technical Report)

This project explored the role of tropical oceans in regulating cloud populations within the Madden-Julian oscillation (MJO)–a large-scale tropical weather disturbance that influences extreme weather events across the globe by affecting the position of the jet stream and patterns of moisture transport. Our project sought to 1) document ocean feedback strength to the MJO in climate models, 2) understand how air-sea interactions affect the intensity, organization, and propagation of MJO convection, 3) utilize existing observations to identify sources of model bias for tropical air-sea interactions, and 4) better understand the effects of slow and fast ocean variability on MJO cloudiness. We found that transfers of heat and moisture across the air-sea interface (i.e., surface fluxes) in most climate models are over-supportive in maintaining MJO convection and under-supportive in propagating MJO convection. This can lead to misrepresentations in the frequency of extreme weather events across the globe through erroneous MJO feedbacks to jet stream position and moisture transport. Furthermore, using observations collected by ocean moorings across the tropical oceans, we showed that most models greatly overestimate surface evaporation in the subtropics, and that this bias is likely one contributor to a pervasive bias among climate models in their ability to simulate the mean state of rainfall over the tropical oceans (i.e., the “double ITCZ bias”). Replacing the default surface flux algorithm in two climate models with a state-of-the-art surface flux algorithm modestly reduced the double ITCZ bias in those models. Our work produced two model diagnostics packages that are now publicly available through DOE and NOAA model diagnostics archives. The first, the Analysis of Scales of Precipitation (ASoP) diagnoses scales of precipitation variability across space and time, with temporal resolution as fine as the model integration timestep. The second, the air-sea flux diagnostic diagnoses model surface fluxes and allows direct comparison of model fluxes to observations without having to recompute fluxes using model output and a state-of-the-art flux algorithm. These diagnostics, which were developed using climate models, can be readily applied to regional and global forecast model output to improve forecasts for the benefit of society.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Combining translational and rotational seismic motions to invert local-scale seismic data for time-variable moment tensors: do rotational motions help for high-frequency seismic data produced by underground explosions?

SUMMARY We present an analysis of combining translational and rotational seismic data in an inversion for the time-variable source time functions corresponding to the components of the seismic moment tensor. We conduct a series of numerical experiments where the data are simulated by a combination of an underground explosion and a co-located double couple shear source and recorded on surface-mounted seismometers within 1–2 km of the source. The experiments are designed to mimic explosion seismology experiments, and thus the data are in the 1–10 Hz frequency range and contain very few surface waves. We use a Monte Carlo method to propagate Earth model uncertainty into the estimates of seismic source parameters. In our experiments, we find that the uncertainty of the estimated seismic source parameters increases when we add rotational seismic motions to the inversion when using a constant number of data channels. In this case, the increased degree of uncertainty in the final results is most likely due to the near-surface Earth model uncertainty that we introduce in our simulations. However, for a fixed number of seismic stations, adding rotational seismic motions to the inversion acts to decrease the uncertainty of the estimated seismic source parameters, most likely due to the increase in the number of data channels used in the inversion.

Poppeliers, Christian (ORCID:0000000159526849)↗

DOE OSTI2022

Advanced Finite-Volume Numerics and Source Term Assumptions for Kernel and G-Equation Modelling of Propane/Air Flames

Here G-Equation models represent propagating flame fronts with an implicit two-dimensional surface representation (level-set). Level-set methods are fast, as transport source terms for the implicit surface can be solved with finite-volume operators on the finite-volume domain, without having to build the actual surface. However, they include approximations whose practical effects are not properly understood. In this study, we improved the numerics of the FRESCO CFD code’s G-Equation solver and developed a new method to simulate kernel growth using signed distance functions and the analytical sphere-mesh overlap. We analyzed their role for simulating propane/air flames, using three well-established constant-volume configurations: a one-dimensional, freely propagating laminar flame; a disc-shaped, constant-volume swirl combustor; and torch-jet flame development through an orifice from a two-chamber device. We tested the explicit (sub-cycled) vs. implicit formulation for the standard transport operators (advection, diffusion, compressibility). In addition to the accurate flame swept-volume method for chemistry and species source term, we developed a more accurate estimator for the burnt/unburnt split cell composition. Then, we developed a signed-distance-function (SDF) based method which provides a more stable reinitialization of the level-set field at every time-step. We found that simplifying assumptions common to several G-Equation implementations, for straightforward terms such as compressibility and advection, lead to large errors in predicting the propagation of even laminar flames, with deviations up to ~300% in simulated vs. formulated flame speed. Conversely, the enhanced numerics enabled through the SDF field reinitialization and improved chemistry source term improve simulation stability and smooth flame propagation even with significantly larger solver time-steps.

42 ENGINEERING↗

DOE OSTI2021

2021 Summer SPE Project

SPE is a project in order to develop new, more physics-based, seismic models of explosions (see Nelson et al. 2013). One key component of this effort is numerical modeling enabled by modern state-of-art code/software. Accurate modeling of the shape and amplitude of seismic waves from their generation to their propagation to remote monitoring seismic stations is important to our ability to determine the origin and strength of the source from remote recording. In this project, the modeling is performed by coupling two codes HOSS and SPECFEM3D. HOSS models the dynamic nonlinear processes happening near the explosion. SPECFEM3D computes the propagation of seismic waves as they travel through 3D complex Earth models. SPECFEM3D was modified in order to be driven by a set of time-series calculated by HOSS in lieu of a seismic source. The goal of this summer project is the investigation of several questions pertinent to the establishment of a full end-to-end modeling capability from the high-rate strain regime area to remote distances where seismic station record seismic waves generated by explosions. The investigated questions are: (1) How to perform proper filtering? Direct-solution modeling only sustains a limited range of frequency depending on the grid size. As we go from one modeling domain to the other via coupling, the mesh size is getting coarser to allow modeling at large scale but also to account for the fact that high-frequency waves do not physically travel to large distances. So filtering of the time-series generated by the near-field hydrodynamic modeling is a current practice often employed but its effects on the modeled waveforms has to be investigated. (2) Quantitative assessment of the efficiency of attenuation to remove high-frequency content of the wavefield. This assessment will allow to create meshes with a grid size appropriate to the actual physics of wave propagation for a given explosion. (3) Checking that the coupling process and the two codes respect the cylindrical symmetry that is expected in the case of a pure explosion in a half-space. (4) The effect of the state-of-stress in the near-source area on the modeled seismic waveforms. These questions will be investigated through the modeling of SPE-4P, the fourth explosion of this series because it was designed to have little interaction with the geologic surrounding and the free surface of the Earth so that it is the most explosion-like experiment, with the most symmetries to be verified.

58 GEOSCIENCES↗

DOE OSTI2022

3D full wave fast wave modeling with realistic HHFW antenna geometry and SOL plasma in NSTX-U

This paper reports the significant advancement of our ability to model and to understand how RF waves interact with the SOL plasma, by developing for the full torus 3D scrape-off layer (SOL) plasma simulation together with the antenna and core plasma in the NSTX-U device in HHFW frequency regime for a cold plasma model. The present paper extends the previous modeling (Bertelli et al 2020 AIP Conf. Proc. 2254 030001) to a full 3D device geometry including realistic antenna geometry in order to capture a more realistic experimental situation including the fully 3D wave field effects and the antenna plasma interaction in the SOL plasma, and, at the same time, resulting core wave propagation. Here, the central tool of the present work is the Petra-M code, which is a newly developed state-of-the-art generic electromagnetic simulation tool for modeling RF wave propagation based on MFEM (http://mfem.org), an open source scalable C++ finite element method library.

3D RF modeling↗

DOE OSTI2025

An infrasound source analysis of the OSIRIS-REx sample return capsule hypersonic re-entry

The OSIRIS-REx sample return capsule's hypersonic re-entry into the atmosphere is a rare opportunity to test a variety of sonic boom source models since the projectile dimensions are well characterized. While the as-flown flight path is unknown, the predicted flight path enables a rough approximation of the source Mach number and location. Six infrasound microphones deployed in the boom carpet along the predicted flight path recorded impulsive signals from the OSIRIS-REx re-entry. Using a suite of atmosphere profiles and the geometric acoustics approximation, we estimate locations with uncertainty estimates along the flight path from which the signals were emitted. Acoustic overpressure and signal duration predictions from Whitham's far field theory, Carlson's simplified sonic boom prediction method, and a drag-dominated hypersonic model are analyzed with uncertainty estimates from the location estimate. While the Carlson simplified sonic boom prediction method could be accurate, our preference is for the drag-dominated source model. Using this source model with an inviscid Burgers's equation solver for propagation, we obtained an excellent match to the recorded data. In conclusion, these results will help better inform future sample return capsule re-entry observation campaigns as well as contribute to a better understanding of high altitude infrasonic sources.

58 GEOSCIENCES↗

DOE OSTI2020

Predicting cell-to-cell failure propagation and limits of propagation in lithium-ion cell stacks

Thermal runaway of lithium-ion batteries is a risk that is magnified when stacks of lithium-ion cells are used for large scale energy storage. When limits of propagation can be identified so that systems can be designed to prevent large scale cascading failure even if a failure does occur, these systems will be safer. This work addresses the prediction of cell-to-cell failure propagation and the propagation limits in lithium-ion cell stacks to better understand and identify safe designs. A thermal-runaway model is presented based on recent developments in thermochemical source terms. It is noted that propagating failure is characterized by temperatures above which calorimetry data is available. Results show high temperature propagating failure predictions are too rapid unless an intra-particle diffusion limit is included, introducing a Damköhler number limiter into the rate expression. This new model form is evaluated against cell-to-cell failure propagation where the end cell of a stack is forced into thermal runaway through a nail-induced short circuit. Limits of propagation for this configuration are identified. Results show cell-to-cell propagation predictions are consistent with measurements over a range of cell states of charge and with the introduction of metal plates between cells to add system heat capacity representative of structural members. This consistency extends from scenarios where propagation occurs through scenarios where propagation is prevented.

25 ENERGY STORAGE↗

DOE OSTI2024

Performance evaluation of the USGS velocity model for the San Francisco Bay Area

In this study, we evaluated the performance of the United States Geological Survey velocity model developed for the San Francisco Bay Area (SFBA), version 21.1. The evaluation was performed through high-resolution three-dimensional physics-based ground motion simulations of seven small-magnitude earthquakes (ranging from magnitude 3.8 to 4.4) that occurred on the eastern side of the San Francisco Bay. The simulations were performed in the frequency range from 0 to 5 Hz with a minimum shear-wave velocity of 250 m/s, which allowed the capture of wave propagation effects of the near-surface soft materials that characterize local basins. Based on the direct comparison of Fourier amplitude spectra between recorded and simulated ground motions for more than 250 stations, we found that the velocity model generally performs well in the frequency range of 0.2–5 Hz. The median value of the Fourier amplitude residuals was found to be near zero for all seven earthquakes. The slight over-prediction of 0.2 log-natural units at frequencies above 3 Hz in our simulations was attributed to the potentially inaccurate representation of the source radiation pattern by a double-couple point source model, and simple representation of shallow small-scale underground structural complexity in the velocity model. Maps of spectral amplitude differences between the simulated and recorded data were used to identify areas responsible for systematic ground motion over-predictions or under-predictions. For example, while some sub-domains over soft sediments show over-prediction patterns, the block east of the Hayward fault is prone to exhibit patterns of under-prediction. These maps can be used to guide future refinements of the SFBA velocity model. Since our simulation methodology allows for the decoupling of the source and wave propagation effects, the ground motion data generated by our simulations can also be used to quantify the epistemic uncertainty due to the velocity model, in empirically based ground motion estimates for the SFBA.

58 GEOSCIENCES↗

DOE OSTI2021

PastoralScape: An Environment-Driven Model of Vaccination Decision Making Within Pastoralist Groups in East Africa

Economic and cultural resilience among pastoralists in East Africa is threatened by the interconnected forces of climate change, contagious diseases spread and evolving national and international trade. A key factor in the resilience of livestock that communities depend on is human decision making regarding vaccination against prevalent diseases such as Rift Valley fever and Contagious Bovine Pleuropneumonia. This paper describes an agent-based model that couples models of disease propagation, animal health, human decision making, and external GIS data sources capturing measures of foraging condition. We describe the design of the sub-models, their coupling, and demonstrate the sensitivity of the model to parameters that relate to controllable factors such as government and NGO information sources that can influence human decision making patterns. This model is intended to form the basis upon which richer economic and human factor models can be built.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Prediction of regional infrasound produced by supersonic sources using a ray-based Mach cone source

The geometry of the Mach cone produced by a supersonic source is analyzed and mapped into initial conditions used in acoustic ray tracing. The resulting source model is combined with spherical geometry ray tracing methods to enable propagation simulations for infrasonic signals produced by bolides, space debris, rockets, aircraft, and other fast-than-sound sources out to typical infrasonic observation distances of hundreds or thousands of kilometers. Idealized linear and parabolic trajectories typical of bolides and rockets, respectively, are used to demonstrate the calculation of regional infrasonic signals produced by such sources and characteristics of the radiated infrasonic waves are found to vary strongly with the geometry of the trajectory and atmospheric structure. Predicted regional infrasonic signals are compared with those observed from a November 2020 bolide that passed over Scandinavia using a combination of institutionally maintained infrasound stations and “citizen scientist” data from the Raspberry Shake data repository.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Prediction of regional infrasound produced by supersonic sources using a ray-based Mach cone source

42 ENGINEERING↗

DOE OSTI2023

Massively Parallel Bayesian Model Calibration and Uncertainty Quantification with Applications to Nuclear Fuels and Materials

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. Fuel behavior in a nuclear reactor is governed by a complex network of mechanisms interacting with various other physics aspects in the reactor system. Any model developed to represent the fuel behavior will likely be idealized resulting in uncertainties in their predictions compared to the observed data. As such, this report was motivated by the need to identify the sources of uncertainties and quantify and propagate them through the fuel model outputs. Such quantification of uncertainties will establish a level of model trustworthiness, identify approaches to improve the model trustworthiness, and even guide optimal experiment design for maximal information gain. To accomplish the uncertainty quantification for computational models, this report has relied on the Bayesian framework which provides probabilistic treatment of models their inputs and outputs. The current state-of-the-art on performing Bayesian Uncertainty Quantification (UQ) for nuclear engineering models using High Performance Computing (HPC) resources have been reviewed. Implementation of capabilities for massively parallel Bayesian UQ in Multiphysics Object-Oriented Simulation Environment (MOOSE) is discussed. Several verification cases are discussed to verify the accuracy of the quantified uncertainties using the developed computational capabilities in MOOSE. Then, the problem of quantifying the uncertainties in TRI-Structural isOtropic (TRISO) fuel silver release is addressed. For the first time, the uncertainties arising from the TRISO Fission Gas Release (FGR) model due to model inadequacy and experimental noise are quantified. Also, the Bayesian capabilities are applied to the calibration of the MATPRO creep model, a widely used model in several fuel assessment cases. The impact of the prediction uncertainties in the MATPRO model on the fuel cladding behavior as part of the TRIBULATION assessment case (which is an integral effects case) is investigated. This report concludes with a discussion on the future work for the UQ for computational models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Illuminating the North Korean nuclear explosion test in 2017 using remote infrasound observations

North Korea conducted its sixth underground nuclear explosion test (mb 6.3) on 2017 September 3. The underground explosion produced substantial low-frequency atmospheric waves, which were detected by infrasound arrays located up to a distance of 566 km. These infrasound waves are formed by the conversion of seismic energy to acoustic energy across the lithosphere–atmosphere interface. While infrasound records at regional distances produce estimates of ground motion amplitude over spatially extended regions covering about 26 500 km 2 , 3-D full seismo-acoustic simulations within the lithosphere and atmosphere provide quantitative information about seismo-acoustic energy partitioning. Our results demonstrate the capability of remote infrasound observations combined with 3-D propagation modelling to further develop discrimination methods for underground sources. These results contribute to enhance the confidence of source identification and characterization in nuclear test monitoring research, which is essential for the enforcement of the Comprehensive Nuclear-Test-Ban Treaty.

58 GEOSCIENCES↗

DOE OSTI2026

Seismic Observations of the OSIRIS-REx Sample Return Capsule Reentry: Deployment, Signal Characteristics, and Wavefield Phenomenology

Controlled spacecraft reentries from interplanetary trajectories provide rare, well-characterized hypersonic sources for advancing seismoacoustic observation techniques. Here we present seismic observations of the OSIRIS-REx sample return capsule (SRC) reentry on 24 September 2023, recorded by 16 three-component nodal seismometers deployed near Eureka, Nevada, at ground distances of 7–20 km from the capsule trajectory. Air-to-ground coupled signals are detected at all stations, exhibiting impulsive onsets consistent with ballistic shock arrivals from the descending Mach cone. We characterize the seismic wavefield through signal amplitude, period, waveform cross-correlation, and array processing. Signal periods decrease systematically with increasing distance from the trajectory within the airport array, indicating that higher-frequency content becomes more prominent at greater offsets, opposite to expectations from geometric spreading and atmospheric absorption. Seismic array processing identifies frequency-dependent back-azimuth variations whose origin remains unresolved; possible contributing factors include source geometry, scattering by fine-scale layered structure in the stratosphere, and near-surface effects. These observations document a spatially complex seismic wavefield from a well-characterized hypersonic line source and provide constraints for future modeling of atmospheric propagation and air-to-ground coupling.

OSIRIS-REx↗

DOE OSTI2026

Hot electron preheat effects in two-photon polymerization plastic lattice laser-driven shock-tube experiments on the OMEGA laser

We present simulations of Omega shock tube experiments designed to investigate hot electron preheat effects in 3D-printed, two-photon polymerization (2PP) plastic lattices. Preheat is inferred in the experiments from the expansion of a plastic witness disk embedded in the lattice. Using the Eulerian radiation-hydrodynamics code xRAGE, we model shock propagation and preheat from both radiative and hot electron energy sources to evaluate their relative impact. To simulate the transport of laser-generated hot electrons, the nonlocal electron heat transport model proposed by Schurtz, Nicolaï, and Busquet (SNB) is extended with a hot electron source term and an energy cascade algorithm. We explore how variations in ablator, lattice geometry, and laser drive affect the shock velocity and witness disk expansion. Simulations show that the inclusion of a 5 μm gold layer reduces shock pressure by 60% and shock speeds by 30%–40% but does not significantly reduce the hot electron preheat, and that different lattice geometries lead to enhanced shock velocities—up to 40% faster than in homogeneous foams. However, radiative and conductive preheat from classical mechanisms alone fail to match experiment. By including a hot electron source term, we reproduce experimental observables such as disk expansion rates and spatial radiographic features. We find that a hot electron population corresponding to 4%–8% of the incident laser energy with T hot = 50 keV produces expansion which agrees with the experimental data, suggesting hot electron preheat is the most plausible explanation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Dynamics of Intra-Cell Thermal Front Propagation in Lithium-Ion Battery Safety Issues

Thermal runaway (TR), a critical failure mode in lithium-ion batteries (LIBs), poses significant safety risks and hinders wider application of LIBs. TR typically begins at a localized heat source and spreads across the cell. Understanding thermal front propagation (TFP) characteristics, such as front and velocity, is crucial for assessing energy release and temperature distribution for battery hazardous estimation. Recent studies assume that TR within cells propagates at a near-constant velocity, based on the reaction kinetics and thermal properties. Here, in this study, an intra-battery TR model is further proposed and it indicates that TFP velocity stabilizes when the front is distanced from the heat source. Theoretical estimates for propagation velocity and front are developed and validated through numerical simulations and experimental tests from the NREL Battery Failure Databank. The energy release rate during TFP and the impact of preheating based on a point heat source are explored. This work clarifies the long-standing clouds of the thermal font propagation behaviors within the single cell, highlights the power and beauty of mathematics modeling to describe the complicated thermal behaviors, and provides important guidelines for thermal hazardous understanding for next-generation batteries.

25 ENERGY STORAGE↗