Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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.

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

At least 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Modeling soil CO<sub>2</sub> production and transport with dynamic source and diffusion terms: testing the steady-state assumption using DETECT v1.0

Abstract. The flux of CO2 from the soil to the atmosphere (soil respiration, Rsoil) is a major component of the global carbon (C) cycle. Methods to measure and model Rsoil, or partition it into different components, often rely on the assumption that soil CO2 concentrations and fluxes are in steady state, implying that Rsoil is equal to the rate at which CO2 is produced by soil microbial and root respiration. Recent research, however, questions the validity of this assumption. Thus, the aim of this work was two-fold: (1) to describe a non-steady state (NSS) soil CO2 transport and production model, DETECT, and (2) to use this model to evaluate the environmental conditions under which Rsoil and CO2 production are likely in NSS. The backbone of DETECT is a non-homogeneous, partial differential equation (PDE) that describes production and transport of soil CO2, which we solve numerically at fine spatial and temporal resolution (e.g., 0.01 m increments down to 1 m, every 6 h). Production of soil CO2 is simulated for every depth and time increment as the sum of root respiration and microbial decomposition of soil organic matter. Both of these factors can be driven by current and antecedent soil water content and temperature, which can also vary by time and depth. We also analytically solved the ordinary differential equation (ODE) corresponding to the steady-state (SS) solution to the PDE model. We applied the DETECT NSS and SS models to the six-month growing season period representative of a native grassland in Wyoming. Simulation experiments were conducted with both model versions to evaluate factors that could affect departure from SS, such as (1) varying soil texture; (2) shifting the timing or frequency of precipitation; and (3) with and without the environmental antecedent drivers. For a coarse-textured soil, Rsoil from the SS model closely matched that of the NSS model. However, in a fine-textured (clay) soil, growing season Rsoil was ∼ 3 % higher under the assumption of NSS (versus SS). These differences were exaggerated in clay soil at daily time scales whereby Rsoil under the SS assumption deviated from NSS by up to 35 % on average in the 10 days following a major precipitation event. Incorporation of antecedent drivers increased the magnitude of Rsoil by 15 to 37 % for coarse- and fine-textured soils, respectively. However, the responses of Rsoil to the timing of precipitation and antecedent drivers did not differ between SS and NSS assumptions. In summary, the assumption of SS conditions can be violated depending on soil type and soil moisture status, as affected by precipitation inputs. The DETECT model provides a framework for accommodating NSS conditions to better predict Rsoil and associated soil carbon cycling processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Stability of Cu(In x Ga 1− x )Se 2 Solar Cells Utilizing RbF Postdeposition Treatment under a Sulfur Atmosphere

Alkali halide postdeposition treatments (PDTs) have become a key tool to maximize efficiency in Cu(In x Ga 1− x )Se 2 (CIGS) photovoltaics. RbF PDTs have emerged as an alternative to the more common Na‐ and K‐based techniques. This study utilizes temperature‐dependent current–voltage ( JVT ) measurements to study a unique RbF PDT performed in a S atmosphere. The samples are measured before and after 6 months in a desiccator to study device stability. Both samples contain Na and K which diffuse from the soda–lime glass substrate. A reference sample and a RbF + S PDT sample both show the development of a rear contact barrier after aging. The contact barrier is higher for the RbF + S PDT sample, leading to decreased current in forward bias. Series resistance is also higher in the RbF + S PDT device which leads to lower fill factor. However, after aging the reference sample has a larger decrease in open‐circuit voltage ( V OC ). Ideality factor measurements suggest Shockley–Read–Hall recombination dominates both samples. V OC versus temperature and a temperature‐dependent activation energy model are used to calculate diode activation energies for each sample condition. Both techniques produce similar values that indicate recombination primarily occurs within the bulk absorber.

14 SOLAR ENERGY↗
DOE OSTI2022

Experimental validation of the mechanistic scale-up methodology of gas–solid spouted beds using radioactive particle tracking (RPT)

The very high-temperature reactors (VHTRs) are highly ranked among candidates of Generation IV of nuclear reactors due to their high efficiency, safety, the resistance to proliferation, and reliability. The VHTRs are preferentially fueled by Tristructural-isometric (TRISO) coated fuel particles which has fuel kernels of fissile material coated by four coating layers: a porous buffer pyrolysis carbon layer (buffer PyC), an inner dense pyrolysis carbon layer (IPyC), a silicon carbide layer (SiC) and an outer dense pyrolysis carbon layer (OPyC). The heart of the operation and safety of the VHTRs significantly depends on the reliability of the coating layers of TRISO particles to retain metallic and gaseous fission products within the particles. The technique used for coating TRISO particles are gas-solids spouted beds via chemical vapor deposition (CVD). Fabrication of high-quality low-defect TRISO fuel particles fuel at larger scale spouted beds is required to support the commercialization of the VHTRs. In this work, our new developed mechanistic scale-up methodology of gas-solids spouted beds based on matching the radial profile of gas-holdup has been demonstrated and validated using radioactive particle tracking (RPT). Two spouted beds of small and larger scales were used in the study. Three sets of conditions were carried out which include the conditions of the reference case in the large scale, conditions that provide similar gas holdup radial profile to that of the reference case and conditions that provided dissimilar gas holdup radial in the small-scale spouted beds. The results confirm the validation of the scale-up methodology in terms of the dimensionless values of the spout diameter, cumulative probability distribution of the solids particles penetration into the spout, fraction of cycle time in each region of the bed, the radial profiles of the dimensionless values of the root-mean-square particle velocities and solids eddy diffusivity. Finally, the results further advance the knowledge and understanding of the gas-solids spouted beds provide deeper insight into their solids dynamics and presenting important benchmarking data for validating computational fluid dynamics codes and models. At last, procedures are established for the implementation of the new scale-up methodology.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Best of both worlds: Enforcing detailed balance in machine learning models of transition rates

The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Distinct Gas-Particle Partitioning and Viscosity Characteristics of Secondary Organic Aerosols Derived from α-Pinene versus Ocimene

Secondary organic aerosols (SOA) have complex, multicomponent composition that controls particle viscosity and gas-particle partitioning, key factors to their atmospheric evolution. This study investigates the chemical composition, volatility and viscosity of SOA formed by ozonolysis of cyclic α-pinene (PSOA) and acyclic ocimene (OSOA) monoterpenes. Using Temperature-Programmed Desorption combined with Direct Analysis in Real-Time ionization and High-Resolution Mass Spectrometry, we determined the molecular composition and saturation mass concentration of individual SOA constituents. These data enabled gas-particle partitioning and viscosity estimates under varied atmospheric conditions. PSOA, composed of higher molecular weight and less oxidized species, shows higher condensability and viscosity under high total organic mass (tOM) loadings. Here, in contrast, OSOA, consisting of more oxidized, lower molecular weight species, exhibits greater sensitivity to tOM, with viscosity increasing significantly upon dilution. Poke-flow experiments support this trend, indicating that OSOA undergoes more dynamic compositional and phase changes during atmospheric aging. These observations reveal distinct dynamic trends in the atmospheric transformations and reactivity of SOA from cyclic and acyclic monoterpenes, with the latter showing greater compositional changes during aging that alter viscosity and diffusion. This highlights the importance of incorporating such dynamic transformations into atmospheric models to improve predictions of SOA atmospheric loadings, lifetimes, and impacts.

cyclic and acyclic monoterpenes↗
DOE OSTI2016

High efficiency 4H-SiC betavoltaic power sources using tritium radioisotopes

Realization of an 18.6% efficient 4H-silicon carbide (4H-SiC) large area betavoltaic power source using the radioisotope tritium is reported. A 200 nm 4H-SiC P{sup +}N junction is used to collect high-energy electrons. The electron source is a titanium tritide (TiH{sup 3}{sub x}) foil, or an integrated titanium tritide region formed by the diffusion of tritium into titanium. The specific activity of the source is directly measured. Dark current measured under short circuit conditions was less than 6.1 pA/cm{sup 2}. Samples measured with an external tritium foil produced an open circuit voltage of 2.09 V, short circuit current of 75.47 nA/cm{sup 2}, fill factor of 0.86, and power efficiency of 18.6%. Samples measured with an integrated source produced power efficiencies of 12%. Simulations were done to determine the beta spectrum (modified by self absorption) exiting the source and the electron hole pair generation function in the 4H-SiC. The electron-hole pair generation function in 4H-SiC was modeled as a Gaussian distribution, and a closed form solution of the continuity equation was used to analyze the cell performance. The effective surface recombination velocity in our samples was found to be 10{sup 5}–10{sup 6 }cm/s. Our analysis demonstrated that the surface recombination dominates the performance of a tritium betavoltaic device but that using a thin P{sup +}N junction structure can mitigate some of the negative effects.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Evaluating the stability and activity of dilute Cu-based alloys for electrochemical CO 2 reduction

Cu-based catalysts currently offer the most promising route to actively and selectively produce value-added chemicals via electrochemical reduction of CO 2 (eCO 2 R); yet further improvements are required for their wide-scale deployment in carbon mitigation efforts. Here, we systematically investigate a family of dilute Cu-based alloys to explore their viability as active and selective catalysts for eCO 2 R through a combined theoretical–experimental approach. Using a quantum–classical modeling approach that accounts for dynamic solvation effects, we assess the stability and activity of model single-atom catalysts under eCO 2 R conditions. Our calculations identify that the presence of eCO2R intermediates, such as CO*, H*, and OH*, may dynamically influence the local catalyst surface composition. Additionally, we identify through binding energy descriptors of the CO*, CHO*, and OCCO* dimer intermediates that certain elements, such as group 13 elements (B, Al, and Ga), enhance the selectivity of C 2+ species relative to pure Cu by facilitating CO dimerization. The theoretical work is corroborated by preliminary testing of eCO 2 R activity and selectivity of candidate dilute Cu-based alloy catalyst films prepared by electron beam evaporation in a zero-gap gas diffusion electrode-based reactor. Of all studied alloys, dilute CuAl was found to be the most active and selective toward C 2+ products like ethylene, consistent with the theoretical predictions. Finally, we attribute the improved performance of dilute CuAl alloys to more favorable dimerization reaction energetics of bound CO species relative to that on pure Cu. In a broader context, the results presented here demonstrate the power of our simulation framework in terms of rational catalyst design.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2015

Ultrafaint dwarf galaxies—The lowest-mass relics from before reionization

New observations suggest that ultrafaint dwarf galaxies (UFDs)—the least luminous systems bound by dark matter halos (≲10{sup 5} L{sub ⊙})—may have formed before reionization. The extrapolated virial masses today are uncertain, with estimates ranging from 10{sup 8} M{sub ⊙} to as high as 10{sup 9} M{sub ⊙} depending on the assumed form of the underlying potential. We show that the progenitor halo masses of UFDs can, in principle, be as low as M{sub vir} ≈ 10{sup 7} M{sub ⊙}. Under the right conditions, such a halo can survive the energy input of a supernova (SN) and its radiative progenitor. A clumpy (fractal) medium is much less susceptible to both internal and external injections of energy. It is less prone to SN sweeping (particularly if it is off-centered) because the coupling efficiency of the explosive energy is much lower than for a diffuse interstellar medium. With the aid of the 3D hydro/ionization code Fyris, we show that sufficient baryons are retained to form stars following a single SN event in dark matter halos down to M{sub vir} ≈ 10{sup 7} M{sub ⊙} in the presence of radiative cooling. In these models, the gas survives the SN explosion, is enriched with the specific abundance yields of the discrete events, and reaches surface densities where low-mass stars can form. Our highest-resolution simulations reveal why cooling is so effective in retaining gas compared to any other factor. In the early stages, the super-hot metal-enriched SN ejecta exhibit strong cooling, leading to much of the explosive energy being lost. Consistent with earlier work, the baryons do not survive in smooth or adiabatic models in the event of an SN. The smallest galaxies may not contribute a large fraction of matter to the formation of galaxies, but they carry signatures of the earliest epochs of star formation, as we show. These signatures may allow us to distinguish a small primordial galaxy from one that was stripped down to its present size through tidal interaction. We discuss these results in the context of local UFDs and damped Lyα systems (z ∼ 2) at very low metallicity ([Fe/H] ∼ −3). We show that both classes of objects are consistent with primordial low-mass systems that have experienced only a few enrichment events.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Temperature-Dependent Short-Circuit Capability of Silicon Carbide Power MOSFETs

Our paper presents a comprehensive short-circuit ruggedness evaluation and numerical investigation of up-to-date commercial silicon carbide (SiC) MOSFETs. The short-circuit capability of three types of commercial 1200-V SiC MOSFETs is tested under various conditions, with case temperatures from 25 to 200 degrees C and dc bus voltages from 400 to 750 V. It is found that the commercial SiC MOSFETs can withstand short-circuit current for only several microseconds with a dc bus voltage of 750 V and case temperature of 200 degrees C. Moreover, the experimental short-circuit behaviors are compared, and analyzed through numerical thermal dynamic simulation. Specifically, an electrothermal model is built to estimate the device internal temperature distribution, considering the temperature-dependent thermal properties of SiC material. Based on the temperature information, a leakage current model is derived to calculate the main leakage current components (i.e., thermal, diffusion, and avalanche generation currents). Finally, numerical results show that the short-circuit failure mechanisms of SiC MOSFETs can be thermal generation current induced thermal runaway or high-temperature-related gate oxide damage.

4H↗
DOE OSTI2016

Ecosystem fluxes of hydrogen in a mid-latitude forest driven by soil microorganisms and plants

Molecular hydrogen (H 2 ) is an atmospheric trace gas with a large microbe-mediated soil sink, yet cycling of this compound throughout ecosystems is poorly understood. Measurements of the sources and sinks of H 2 in various ecosystems are sparse, resulting in large uncertainties in the global H 2 budget. Constraining the H 2 cycle is critical to understanding its role in atmospheric chemistry and climate. We measured H 2 fluxes at high frequency in a temperate mixed deciduous forest for 15 months using a tower-based flux-gradient approach to determine both the soil-atmosphere and the net ecosystem flux of H 2 . We also found that Harvard Forest is a net H 2 sink (-1.4 ± 1.1 kg H 2 ha -1 ) with soils as the dominant H 2 sink (-2.0 ± 1.0 kg H 2 ha -1 ) and aboveground canopy emissions as the dominant H 2 source (+0.6 ± 0.8 kg H 2 ha -1 ). Aboveground emissions of H 2 were an unexpected and substantial component of the ecosystem H 2 flux, reducing net ecosystem uptake by 30% of that calculated from soil uptake alone. Soil uptake was highly seasonal (July maximum, February minimum), positively correlated with soil temperature and negatively correlated with environmental variables relevant to diffusion into soils (i.e., soil moisture, snow depth, snow density). Soil microbial H 2 uptake was correlated with rhizosphere respiration rates (r = 0.8, P < 0.001), and H 2 metabolism yielded up to 2% of the energy gleaned by microbes from carbon substrate respiration. Here, we elucidate key processes controlling the biosphere–atmosphere exchange of H 2 and raise new questions regarding the role of aboveground biomass as a source of atmospheric H 2 and mechanisms linking soil H 2 and carbon cycling. Our results should be incorporated into modeling efforts to predict the response of the H 2 soil sink to changes in anthropogenic H 2 emissions and shifting soil conditions with climate and land-use change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

On the role of soil water retention characteristic on aerobic microbial respiration

Soil water status is one of the most important environmental factors that control microbial activity and rate of soil organic matter (SOM) decomposition. Its effect can be partitioned into effect of water energy status (water potential) on cellular activity, effect of water volume on cellular motility, and aqueous diffusion of substrate and nutrients, as well as the effect of air content and gas-diffusion pathways on concentration of dissolved oxygen. Nonetheless,moisture functions widely used in SOM decomposition models are often based on empirical functions rather than robust physical foundations that account for these disparate impacts of soil water. The contributions of soil water content and water potential varyfrom soil to soil according to the soil water characteristic (SWC), which in turn is strongly dependent on soil texture and structure. The overall goal of this study is to introduce a physically based modeling framework of aerobic microbial respiration that incorporates the role of SWC under arbitrary soil moisture status. The model was tested by contrasting it with published data sets of SOM decomposition under laboratory conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A model for droplet heating and evaporation of water-in-oil emulsified fuel

Water-in-oil (W/O) emulsified fuel is a promising alternative fuel by inducing flash boiling of water at high temperature that can improve the atomization of fuel spray. The heating and evaporation process of emulsified fuel droplet is affected by the diffusion and coalescence of small dispersed water droplets in oil during the heating. In this study, a model is developed for the droplet heating and evaporation of W/O emulsified fuel with these key physics considered. The diffusivity of dispersed water droplets in oil is calculated using the Stokes-Einstein equation. The deactivation temperature of surfactant is selected as a criterion for the water coalescence since the water coalescence is due to the deactivation of surfactant. The process of water coalescence is simplified that the dispersed water droplets coalesce into one single water sub-droplet at the center of the oil droplet instantly as the droplet temperature reaches the deactivation temperature of surfactant. The model is validated against experimental data of single droplet under different heating temperatures, surfactant concentrations, and sizes of dispersed water droplets. Based on the proposed model, the effects of fuel properties and heating conditions on the droplet heating and evaporation are analyzed.

Droplet coalescence↗
DOE OSTI2017

A finite-volume HLLC-based scheme for compressible interfacial flows with surface tension

Shock waves are often used in experiments to create a shear flow across liquid droplets to study secondary atomization. Similar behavior occurs inside of supersonic combustors (scramjets) under startup conditions, but it is challenging to study these conditions experimentally. In order to investigate this phenomenon further, a numerical approach is developed to simulate compressible multiphase flows under the effects of surface tension forces. The flow field is solved via the compressible multicomponent Euler equations (i.e., the five equation model) discretized with the finite volume method on a uniform Cartesian grid. The solver utilizes a total variation diminishing (TVD) third-order Runge–Kutta method for time-marching and second order TVD spatial reconstruction. Surface tension is incorporated using the Continuum Surface Force (CSF) model. Fluxes are upwinded with a modified Harten–Lax–van Leer Contact (HLLC) approximate Riemann solver. An interface compression scheme is employed to counter numerical diffusion of the interface. The present work includes modifications to both the HLLC solver and the interface compression scheme to account for capillary force terms and the associated pressure jump across the gas–liquid interface. A simple method for numerically computing the interface curvature is developed and an acoustic scaling of the surface tension coefficient is proposed for the non-dimensionalization of the model. The model captures the surface tension induced pressure jump exactly if the exact curvature is known and is further verified with an oscillating elliptical droplet and Mach 1.47 and 3 shock-droplet interaction problems. The general characteristics of secondary atomization at a range of Weber numbers are also captured in a series of simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Schottky barrier height measurements of Cu/Si(001), Ag/Si(001), and Au/Si(001) interfaces utilizing ballistic electron emission microscopy and ballistic hole emission microscopy

The Schottky barrier heights of both n and p doped Cu/Si(001), Ag/Si(001), and Au/Si(001) diodes were measured using ballistic electron emission microscopy and ballistic hole emission microscopy (BHEM), respectively. Measurements using both forward and reverse ballistic electron emission microscopy (BEEM) and (BHEM) injection conditions were performed. The Schottky barrier heights were found by fitting to a linearization of the power law form of the Bell-Kaiser BEEM model. The sum of the n-type and p-type barrier heights are in good agreement with the band gap of silicon and independent of the metal utilized. The Schottky barrier heights are found to be below the region of best fit for the power law form of the BK model, demonstrating its region of validity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Axial magnetic flux amplification in Hall-magnetohydrodynamic simulations of externally magnetized z-pinches

It is shown that externally magnetized gas-puff z-pinches under compression can exhibit robust axial flux amplification under suitable axial boundary conditions. This effect relies upon the Hall term in the generalized Ohm's law to generate azimuthal currents in the presence of a driving axial electric field. Under dynamic compression, the total current tends to flow in a mostly force-free boundary layer, separating the predominately azimuthal field outside the layer and the predominately axial field inside the layer. The effect only appears to occur if the axial boundaries allow for outflow or absorption of electromagnetic energy. The effect is mitigated by imposing either periodic or conducting axial boundary conditions. A semi-analytic equilibrium analysis agrees with steady-state solutions of the time-dependent electron-magnetohydrodynamic equations and provides an estimate of the scaling of the boundary layer as well as suggesting a scenario for the formation of the boundary layer. When operative, the effect can significantly impede plasma compression due to the increase in axial flux that diffuses into the pre-compressed plasma or through the presence of conductors that inhibit movement of the generated axial flux. Several facilities have noted unusual implosion behavior in z-pinch experiments with applied axial fields that does not appear to be explainable within the standard magnetohydrodynamic model. In conclusion, it is suggested that these experiments can be explained by the axial flux amplification and concentration effect

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Iron Ore Reduction by Hydrogen Using a Laboratory Scale Fluidized Bed Reactor: Kinetic Investigation—Experimental Setup and Method for Determination

The reduction kinetics of hematite iron ore fines to metallic iron by hydrogen using a laboratory fluidized bed reactor were investigated in a temperature range between 873 K and 1073 K, by measuring the weight change of the sample portion during reduction. The fluidization conditions were checked regarding plausibility within the Grace diagram and the measured pressure drop across the material during experiments. The apparent activation energy of the reduction was determined against the degree of reduction and varied along an estimated two-peak curve between 11 and 55 kJ mol{sup −1}. Conventional kinetic analysis for the reduction of FeO to metallic iron, using typical models to describe gas–solid reactions, does not show results with high accuracy. Multistep kinetic analysis, using the Johnson–Mehl–Avrami model, shows that the initial stage of reduction from Fe{sub 2}O{sub 3} to Fe{sub 3}O{sub 4}, and partly to FeO, is controlled by diffusion and chemical reaction, depending on the temperature. Further reduction can be described by a combination of nucleation and chemical reaction, whereby the influence of nucleation increases with an increasing reduction temperature. The results of the kinetical analysis were linked to the shape of the curve from apparent activation energy against the degree of reduction.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Scaling Relations for Intercalation Induced Damage in Electrodes

Mechanical degradation, owing to intercalation induced stress and microcrack formation, is a key contributor to the electrode performance decay in lithium-ion batteries (LIBs). The stress generation and formation of microcracks are caused by the solid state diffusion of lithium in the active particles. Here in this work, scaling relations are constructed for diffusion induced damage in intercalation electrodes based on an extensive set of numerical experiments with a particle-level description of microcrack formation under disparate operating and cycling conditions, such as temperature, particle size, C-rate, and drive cycle. The microcrack formation and evolution in active particles is simulated based on a stochastic methodology. A reduced order scaling law is constructed based on an extensive set of data from the numerical experiments. The scaling relations include combinatorial constructs of concentration gradient, cumulative strain energy, and microcrack formation. Lastly, the reduced order relations are further employed to study the influence of mechanical degradation on cell performance and validated against the high order model for the case of damage evolution during variable current vehicle drive cycle profiles.

25 ENERGY STORAGE↗