We report a projection-based reduced order model (pROM) methodology has been developed for transient heat transfer problems involving coupled conduction and enclosure radiation. The approach was demonstrated on two test problems of varying complexity. The reduced order models demonstrated substantial speedups (up to 185×) relative to the full order model with good accuracy (less than 3% L∞ error). An attractive feature of pROMs is that there is a natural error indicator for the ROM solution: the final residual norm at each time-step of the converged ROM solution. Using example test cases, we discuss how to interpret this error indicator to assess the accuracy of the ROM solution. The approach shows promise for many-query applications, such as uncertainty quantification and optimization. The reduced computational cost of the ROM relative to the full-order model (FOM) can enable the analysis of larger and more complex systems as well as the exploration of larger parameter spaces.
Battery electrodes are composed of polydisperse particles and a porous, composite binder domain. These materials are arranged into a complex mesostructure whose morphology impacts both electrochemical performance and mechanical response. We present image-based, particle-resolved, mesoscale finite element model simulations of coupled electrochemical-mechanical performance on a representative NMC electrode domain. Beyond predicting macroscale quantities such as half-cell voltage and evolving electrical conductivity, studying behaviors on a per-particle and per-surface basis enables performance and material design insights previously unachievable. Voltage losses are primarily attributable to a complex interplay between interfacial charge transfer kinetics, lithium diffusion, and, locally, electrical conductivity. Mesoscale heterogeneities arise from particle polydispersity and lead to material underutilization at high current densities. Particle-particle contacts, however, reduce heterogeneities by enabling lithium diffusion between connected particle groups. While the porous composite binder domain (CBD) may have slower ionic transport and less available area for electrochemical reactions, its high electrical conductivity makes it the preferred reaction site late in electrode discharge. Mesoscale results are favorably compared to both experimental data and macrohomogeneous models. This work enables improvements in materials design by providing a tool for optimization of particle sizes, CBD morphology, and manufacturing conditions.
Inverse problems arise in a wide range of applications, whenever unknown model parameters cannot be measured directly. Instead, the unknown parameters are estimated using experimental data and forward simulations. Thermal inverse problems, such as material characterization problems, are often large-scale and transient. Therefore, they require intrusive adjoint-based gradient implementations in order to be solved efficiently. The capability to solve large-scale transient thermal inverse problems using an adjoint-based approach was recently implemented in SNL Sierra Mechanics, a massively parallel capable multiphysics code suite. This report outlines the theory, optimization formulation, and path taken to implement thermal inverse capabilities in Sierra within a unit test framework. The capability utilizes Sierra/Aria and Sierra/Fuego data structures, the Rapid Optimization Library, and an interface to the Sierra/InverseOpt library. The existing Sierra/Aria time integrator is leveraged to implement a time-dependent adjoint solver.
As computing power rapidly increases, quickly creating a representative and accurate discretization of complex geometries arises as a major hurdle towards achieving a next generation simulation capability. Component definitions may be in the form of solid (CAD) models or derived from 3D computed tomography (CT) data, and creating a surface-conformal discretization may be required to resolve complex interfacial physics. The Conformal Decomposition Finite Element Methods (CDFEM) has been shown to be an efficient algorithm for creating conformal tetrahedral discretizations of these implicit geometries without manual mesh generation. In this work we describe an extension to CDFEM to accurately resolve the intersections of many materials within a simulation domain. This capability is demonstrated on both an analytical geometry and an image-based CT mesostructure representation consisting of hundreds of individual particles. Effective geometric and transport properties are the calculated quantities of interest. Solution verification is performed, showing CDFEM to be optimally convergent in nearly all cases. Representative volume element (RVE) size is also explored and per-sample variability quantified. Relatively large domains and small elements are required to reduce uncertainty, with recommended meshes of nearly 10 million elements still containing upwards of 30% uncertainty in certain effective properties. This work instills confidence in the applicability of CDFEM to provide insight into the behaviors of complex composite materials and provides recommendations on domain and mesh requirements.
The objective of this project is to improve the fidelity of battery-scale simulations of abuse scenarios through the creation and application of microscale (particle-scale) electrode simulations.
Polymer foam encapsulants provide mechanical, electrical, and thermal isolation in engineered systems. In fire environments, gas pressure from thermal decomposition of polymers can cause mechanical failure of sealed systems. A 3-D finite element conduction-radiation model with porous media flow and a chemistry model was created to investigate the heat transfer and pressurization in such scenarios. Experiments show that the rate of pressurization and the temperature of select thermocouples are dependent on orientation with respect to gravity, indicating buoyancy-driven flow. In this work, the gas velocity is solved by applying the Darcy approximation, and the heat transfer and pressurization are determined by solving the continuity, species, and enthalpy equations in the condensed and gas phases. This work will describe the porous media model, explore material parameters (e.g. phase, permeability, conductivity) for use with PMDI polyurethane, compare predictions to experimental data, and recommend values for material properties. It will use multiple heating rates to validate the data, and show that incorporating gas motion into the model captures the divergent nature of the results in different orientations.
Lithium-ion battery electrodes are composed of active material particles, binder, and conductive additives that form an electrolyte-filled porous particle composite. The mesoscale (particle-scale) interplay of electrochemistry, mechanical deformation, and transport through this tortuous multi-component network dictates the performance of a battery at the cell-level. Effective electrode properties connect mesoscale phenomena with computationally feasible battery-scale simulations. We utilize published tomography data to reconstruct a large subsection (1000+ particles) of an NMC333 cathode into a computational mesh and extract electrode-scale effective properties from finite element continuum-scale simulations. We present a novel method to preferentially place a composite binder phase throughout the mesostructure, a necessary approach due difficulty distinguishing between non-active phases in tomographic data. We compare stress generation and effective thermal, electrical, and ionic conductivities across several binder placement approaches. Isotropic lithiation-dependent mechanical swelling of the NMC particles and the consideration of strain-dependent composite binder conductivity significantly impact the resulting effective property trends and stresses generated. Our results suggest that composite binder location significantly affects mesoscale behavior, indicating that a binder coating on active particles is not sufficient and that more accurate approaches should be used when calculating effective properties that will inform battery-scale models in this inherently multi-scale battery simulation challenge.
Battery performance, while observed at the macroscale, is primarily governed by the bicontinuous mesoscale network of the active particles and a polymeric conductive binder in its electrodes. Manufacturing processes affect this mesostructure, and therefore battery performance, in ways that are not always clear outside of empirical relationships. Directly studying the role of the mesostructure is difficult due to the small particle sizes (a few microns) and large mesoscale structures. Mesoscale simulation, however, is an emerging technique that allows the investigation into how particle-scale phenomena affect electrode behavior. In this manuscript, we discuss our computational approach for modeling electrochemical, mechanical, and thermal phenomena of lithium-ion batteries at the mesoscale. We review our recent and ongoing simulation investigations and discuss a path forward for additional simulation insights.
As LiCoO2 cathodes are charged, delithiation of the LiCoO2 active material leads to an increase in the lattice spacing, causing swelling of the particles. When these particles are packed into a bicontinuous, percolated network, as is the case in a battery electrode, this swelling leads to the generation of significant mechanical stress. In this study we performed coupled electrochemical-mechanical simulations of the charging of a LiCoO2 cathode in order to elucidate the mechanisms of stress generation and the effect of charge rate and microstructure on these stresses. Energy dispersive spectroscopy combined with scanning electron microscopy imaging was used to create 3D reconstructions of a LiCoO2 cathode, and the Conformal Decomposition Finite Element Method is used to automatically generate computational meshes on this reconstructed microstructure. Replacement of the ideal solution Fickian diffusion model, typically used in battery simulations, with a more general non-ideal solution model shows substantially smaller gradients of lithium within particles than is typically observed in the literature. Using this more general model, lithium gradients only appear at states of charge where the open-circuit voltage is relatively constant. While lithium gradients do affect the mechanical stress state in the particles, the maximum stresses are always found in the fully-charged state and are strongly affected by the local details of the microstructure and particle-to-particle contacts. These coupled electrochemical-mechanical simulations begin to yield insight into the partitioning of volume change between reducing pore space and macroscopically swelling the electrode. Finally, preliminary studies that include the presence of the polymeric binder suggest that it can greatly impact stress generation and that it is an important area for future research.
We are studying PMDI polyurethane with a fast catalyst, such that filling and polymerization occur simultaneously. The foam is over-packed to tw ice or more of its free rise density to reach the density of interest. Our approach is to co mbine model development closely with experiments to discover new physics, to parameterize models and to validate the models once they have been developed. The model must be able to repres ent the expansion, filling, curing, and final foam properties. PMDI is chemically blown foam, wh ere carbon dioxide is pr oduced via the reaction of water and isocyanate. The isocyanate also re acts with polyol in a competing reaction, which produces the polymer. A new kinetic model is developed and implemented, which follows a simplified mathematical formalism that decouple s these two reactions. The model predicts the polymerization reaction via condensation chemis try, where vitrification and glass transition temperature evolution must be included to correctly predict this quantity. The foam gas generation kinetics are determined by tracking the molar concentration of both water and carbon dioxide. Understanding the therma l history and loads on the foam due to exothermicity and oven heating is very important to the results, since the kinetics and ma terial properties are all very sensitive to temperature. The conservation eq uations, including the e quations of motion, an energy balance, and thr ee rate equations are solved via a stabilized finite element method. We assume generalized-Newtonian rheology that is dependent on the cure, gas fraction, and temperature. The conservation equations are comb ined with a level set method to determine the location of the free surface over time. Results from the model are compared to experimental flow visualization data and post-te st CT data for the density. Seve ral geometries are investigated including a mock encapsulation part, two configur ations of a mock stru ctural part, and a bar geometry to specifically test the density model. We have found that the model predicts both average density and filling profiles well. However, it under predicts density gradients, especially in the gravity direction. Thoughts on m odel improvements are also discussed.
Lithium-ion battery electrodes rely on a percolated network of solid particles and binder that must maintain a high electronic conductivity in order to function. Coupled mechanical and electrochemical simulations may be able to elucidate the mechanisms for capacity fade. We present a framework for coupled simulations of electrode mechanics that includes swelling, deformation, and stress generation driven by lithium intercalation. These simulations are performed at the mesoscale, which requires 3D reconstruction of the electrode microstructure from experimental imaging or particle size distributions. We present a novel approach for utilizing these complex reconstructions within a finite element code. A mechanical model that involves anisotropic swelling in response to lithium intercalation drives the deformation. Stresses arise from small-scale particle features and lithium concentration gradients. However, we demonstrate, for the first time, that the largest stresses arise from particle-to-particle contacts, making it important to accurately represent the electrode microstructure on the multi-particle scale. Including anisotropy in the swelling mechanics adds considerably more complexity to the stresses and can significantly enhance peak particle stresses. Shear forces arise at contacts due to the misorientation of the lattice structure. These simulations will be used to study mechanical degradation of the electrode structure through charge/discharge cycles.
We develop a capability to simulate reduction-oxidation (redox) flow batteries in the Sierra Multi-Mechanics code base. Specifically, we focus on all-vanadium redox flow batteries; however, the capability is general in implementation and could be adopted to other chemistries. The electrochemical and porous flow models follow those developed in the recent publication by [28]. We review the model implemented in this work and its assumptions, and we show several verification cases including a binary electrolyte, and a battery half-cell. Then, we compare our model implementation with the experimental results shown in [28], with good agreement seen. Next, a sensitivity study is conducted for the major model parameters, which is beneficial in targeting specific features of the redox flow cell for improvement. Lastly, we simulate a three-dimensional version of the flow cell to determine the impact of plenum channels on the performance of the cell. Such channels are frequently seen in experimental designs where the current collector plates are borrowed from fuel cell designs. These designs use a serpentine channel etched into a solid collector plate.