SEARCH · Engineering Papers

Results for “Time Evolution”

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 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Situation awareness and dynamic ensemble forecasting of abnormal behavior in cyber-physical system

A plurality of monitoring nodes may each generate a time-series of current monitoring node values representing current operation of a cyber-physical system. A feature-based forecasting framework may receive the time-series of and generate a set of current feature vectors using feature discovery techniques. The feature behavior for each monitoring node may be characterized in the form of decision boundaries that separate normal and abnormal space based on operating data of the system. A set of ensemble state-space models may be constructed to represent feature evolution in the time-domain, wherein the forecasted outputs from the set of ensemble state-space models comprise anticipated time evolution of features. The framework may then obtain an overall features forecast through dynamic ensemble averaging and compare the overall features forecast to a threshold to generate an estimate associated with at least one feature vector crossing an associated decision boundary.

Abbaszadeh, Masoud↗

DOE OSTI2026

Roadmap on advanced and / real-time characterisation of solid state materials and devices for energy applications

A strong societal and political drive is motivating the development and optimization of novel energy conversion and storage systems for decarbonization. The successful implementation of solid state devices such as fuel cells and secondary batteries depends, however, on achieving ambitious targets in terms of performance, reliability and cost competitiveness. Research and technology are addressing these needs through a holistic approach including exploration of new materials and nanoarchitectures, as well as system engineering. These significant efforts require the support of appropriate characterization tools capable of assessing nanometer-scale phenomena such as concentration profiles of ionic and electronic charges, local chemical compositions and their evolution over time across interfaces. This roadmap provides an overview of selected advanced characterization techniques for energy materials and devices. Specific focus is put on in situ/operando methods for probing electrochemical phenomena in real-time under realistic working conditions. Experts in the field provide an extensive review of the current state of the art in 2025 and the current and future challenges for the characterization of local chemistry and kinetics in the bulk of the material, in nanoarchitectures (e.g. thin films) and at the interfaces (e.g. grain boundaries, phase contacts, solid/liquid and solid/gas interfaces) . The aim is to provide a detailed guide to the techniques, describing opportunities and bottlenecks for their practical deployment and examples of successful applications.

25 ENERGY STORAGE↗

DOE OSTI2022

Analytical solution for nonadiabatic quantum annealing to arbitrary Ising spin Hamiltonian

Abstract Ising spin Hamiltonians are often used to encode a computational problem in their ground states. Quantum Annealing (QA) computing searches for such a state by implementing a slow time-dependent evolution from an easy-to-prepare initial state to a low energy state of a target Ising Hamiltonian of quantum spins, H I . Here, we point to the existence of an analytical solution for such a problem for an arbitrary H I beyond the adiabatic limit for QA. This solution provides insights into the accuracy of nonadiabatic computations. Our QA protocol in the pseudo-adiabatic regime leads to a monotonic power-law suppression of nonadiabatic excitations with time T of QA, without any signature of a transition to a glass phase, which is usually characterized by a logarithmic energy relaxation. This behavior suggests that the energy relaxation can differ in classical and quantum spin glasses strongly, when it is assisted by external time-dependent fields. In specific cases of H I , the solution also shows a considerable quantum speedup in computations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Broadband Stability of the Habitable Zone Planet Finder Fabry–Pérot Etalon Calibration System: Evidence for Chromatic Variation

The comblike spectrum of a white light–illuminated Fabry–Pérot etalon can serve as a cost-effective and stable reference for precise Doppler measurements. Understanding the stability of these devices across their broad (hundreds of nanometers) spectral bandwidths is essential to realizing their full potential as Doppler calibrators. However, published descriptions remain limited to small bandwidths or short time spans. We present an ∼6 month broadband stability monitoring campaign of the Fabry–Pérot etalon system deployed with the near-infrared Habitable Zone Planet Finder (HPF) spectrograph. We monitor the wavelengths of each of ∼3500 resonant modes measured in HPF spectra of this Fabry–Pérot etalon (free spectral range = 30 GHz, bandwidth = 820–1280 nm), leveraging the accuracy and precision of an electro-optic frequency comb reference. These results reveal chromatic structure in the Fabry–Pérot mode locations and their evolution with time. We measure an average drift on the order of 2 cm s{sup –1} day{sup −1}, with local departures up to ±5 cm s{sup –1} day{sup −1}. We discuss these behaviors in the context of the Fabry–Pérot etalon mirror dispersion and other optical properties of the system and the implications for the use of similar systems for precise Doppler measurements. Our results show that this system supports the wavelength calibration of HPF at the ≲10 cm s{sup −1} level over a night and the ≲30 cm s{sup −1} level over ∼10 days. Our results also highlight the need for long-term and spectrally resolved study of similar systems that will be deployed to support Doppler measurement precision approaching ∼10 cm s{sup −1}.

42 ENGINEERING↗

DOE OSTI2022

Imaginary Time Propagation on a Quantum Chip

We report evolution in imaginary time is a prominent technique for finding the ground state of quantum many-body systems, and the heart of a number of numerical methods that have been used with great success in quantum chemistry, condensed matter, and nuclear physics. We propose an algorithm to implement imaginary time propagation on a quantum computer. Our algorithm is devised in the context of an efficient encoding into an optimized gate, drawing on the underlying characteristics of the quantum device of a unitary operation in an extended Hilbert space. However, we prove that for simple problems it can also be successfully applied to standard digital quantum machines. This work paves the way for porting quantum many-body methods based on imaginary-time propagation to near-term quantum devices, enabling the future quantum simulation of the ground states of a broad class of microscopic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Real-Time Alignment and Reorientation of Polymer Chains in Liquid Crystal Elastomers

Liquid crystal elastomers (LCEs) exhibit soft elasticity due to the alignment and reorientation of mesogens upon mechanical loading, which provides additional mechanisms to absorb and dissipate energy. This enhanced response makes LCEs potentially transformative materials for biomedical devices, tissue replacements, and protective equipment. However, there is a critical knowledge gap in understanding the highly rate-dependent dissipative behaviors of LCEs due to the lack of real-time characterization techniques that probe the microscale network structure and link it to the mechanical deformation of LCEs. In this work, we employ in situ optical measurements to evaluate the alignment and reorientation degree of mesogens in LCEs. The data are correlated to the quantitative physical analysis using polarized Fourier-transform infrared spectroscopy. The time scale of mesogen alignment is determined at different strain levels and loading rates. The mesogen reorientation kinetics is characterized to establish its relationship with the macroscale tensile strain, and compared to theoretical predictions. Overall, this work provides the first detailed study on the time-dependent evolution of mesogen alignment and reorientation in deformed LCEs. It also provides an effective and more accessible approach for other researchers to investigate the structural-property relationships of different types of polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Unraveling transition-metal-mediated stability of spinel oxide via in situ neutron scattering

The energy materials performance is intrinsically determined by structures from the average lattice structure to the atom arrangement, valence, and distribution of the containing transition metal (TM) elements. Understanding the mechanism of the structure transition and atom rearrangement via synthesis or processing is key to expediting the exploration of excellent energy materials. In this work, in situ neutron scattering is employed to reveal the real-time structure evolution, including the TM-O bonds, lattice, TM valence and the migration of the high-voltage spinel cathode LiNi 0.5 Mn 1.5 O 4 . The transition-metal-mediated spinel destabilization under the annealing at the oxygen-deficient atmosphere is pinpointed. Additionally, the formation of Mn 3+ is correlated to the TM migration activation, TM disordered rearrangement in the spinel, and the transition to a layered-rocksalt phase. The further TM interdiffusion and Mn 2+ reduction are also revealed with multi-stage thermodynamics and kinetics. The mechanisms of phase transition and atom migrations as functions of temperature, time and atmosphere present important guidance on the synthesis in various-valence element containing oxides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Modified gravity constraints from the full shape modeling of clustering measurements from DESI 2024

We present cosmological constraints on deviations from general relativity (GR) from the first-year of clustering observations from the Dark Energy Spectroscopic Instrument (DESI) in combination with other available datasets including the CMB data from Planck with CMB-lensing from Planck and ACT, BBN constraints on the physical baryon density, the galaxy weak lensing and clustering from DESY3 and supernova data from DESY5. We first consider the μ(a,k)–Σ(a,k) modified gravity (MG) parameterization (as well as η(a,k)) in a ΛCDM and a w 0 w a CDM cosmological backgrounds. Using a functional form for time-only evolution gives μ 0 = 0.11 +0.44 -0.54 from DESI(FS+BAO)+BBN and a wide prior on n s . Using DESI(FS+BAO)+CMB+DESY3+DESY5-SN, we obtain μ 0 = 0.05 ± 0.22 and Σ 0 = 0.008 ± 0.045 and similarly μ 0 = 0.02 +0.19 -0.24 and η 0 = 0.09 +0.36 -0.60 , in an ΛCDM background. In w 0 w a CDM we obtain μ 0 = -0.24 +0.32 -0.28 and Σ 0 = 0.006 ± 0.043, consistent with GR, and we still find a preference of the data for a dynamical dark energy with w 0 > -1 and w a < 0. Using functional dependencies in both time and scale gives μ 0 and Σ 0 with a same level of precision as above but other scale MG parameters remain hard to constrain. We then move to binned parameterizations in a ΛCDM background starting with two bins in redshift and obtain, μ 1 = 1.02 ± 0.13, μ 2 = 1.04 ± 0.11, Σ 1 = 1.021 ± 0.029 and Σ 2 = 1.022 +0.027 -0.023 , all consistent with the unity value of GR in the binning formalism. We then extend the analysis to combine two bins in redshift and two in scale giving 8 MG parameters that we find all consistent with GR. We note that we find here that the tension reported in previous studies about Σ 0 being inconsistent with GR when using Planck PR3 data goes away when we use the recent LoLLiPoP+HiLLiPoP likelihoods. As noted in previous studies, this seems to indicate that the tension is indeed related to the CMB lensing anomaly in PR3 which is also resolved when using the recent likelihoods. We then constrain the class of Horndeski theory in the effective field theory of dark energy approach. We consider both EFT-basis and α-basis in the analysis. Assuming a power law parameterization for the EFT function Ω, which controls non-minimal coupling, we obtain Ω 0 = 0.012 +0.001 -0.012 and s 0 = 0.996 +0.54 -0.20 from the combination of DESI(FS+BAO)+DESY5SN+CMB in a ΛCDM background, which are consistent with GR. Similar results are obtained when using the α-basis and assuming no-braiding (α B = 0) giving c M < 1.14 at 95% CL in a ΛCDM background, also in agreement with GR. However, we see a mild yet consistent indication for c B > 0 when α B is allowed to vary which will require further study to determine whether this is due to systematics or new physics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Numerical simulation of Li-pellet injection experiments for ELM-pacing in EAST

The injection of light impurity granules into fusion research discharges has shown a great potential for ELM pacing in future nuclear devices. The present work employed a model considering both pellet ablation and transport to systematically study for the first time the evolution of both impurity Li species and target plasma in the edge plasma during single and repeated Li pellet injections, with all initial information extracted from the EAST experiments. The evolutionary distributions of atomic and ionized Li species were presented and show good consistency with experimental CCD camera images. Simulation studies were carried out for single and repeated pellet injections at the outer mid-plane, and statistical uncertainty of ELM pacing rates triggering ELMs is discussed with regard to related injection conditions. It is found that the pressure locally caused by an ELM pacing pellet (with a diameter of 0.6 ~ 0.7 mm) was 1.8–2.4 times its unperturbed magnitude in the EAST shots of interest. Accumulation from repeated pellet injections results in the rise of pedestal height at a high injection frequency and this effect becomes stronger with increasing pellet size. Pressure changes are also evident in a scenario of multiple small pellets per stroke, which could lead to a nonnegligible, probabilistic ELM triggering rate observed experimentally. The effects of other possible factors are also discussed. With the assumption that the ELM trigger rate of a Li pellet of 0.7 mm diameter is 100%, the predicted trigger rates are in good agreement with those observed in experiments. This work can help understand the statistically probabilistic behavior of ELM triggering efficiency observed in previous studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Multi-fluid and kinetic models of partially ionized magnetic reconnection

Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomenon in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet–Parker scaling, where neutrals are unimportant, or a fully coupled Sweet–Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfvén speed, $v^*_A$. In this work, these differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $ν_{ni}L/v^*_A$ ≳1, where $ν_{ni}$ is the neutral–ion collision frequency and L is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $ν_{ni}L/v^*_A$ <1, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Multiscale Explanation of the Missing Gallium Vacancy in Gallium Arsenide

Irradiation of gallium arsenide (GaAs) produces immobile vacancies and mobile interstitials. Yet, after decades of experimental investigation, the immobile Ga vacancy continues to evade detection, raising the question: where is the Ga vacancy? Static first-principles calculations predict a Ga vacancy should be readily observed. We find that short-time dynamical evolution of primary defects is the key to explaining this conundrum. Using a dynamical multiscale atomistically informed device engineering (AIDE) method, we discover that during the initial displacement damage, the Ga vacancy (3-/2-) defect level pins the Fermi level near the midgap, producing oppositely charged vacancies and interstitials. Driven by Coulomb attraction, fast As interstitials preferentially annihilate Ga vacancies. The Ga vacancy population plummets below detectable limits—and the now unpinned Fermi level recovers—before being experimentally observed. This dynamical model solves the mystery of the missing Ga vacancy and reveals the importance of a multiscale approach to explore the dynamical chemical behavior in experimentally inaccessible short-time regimes.

Diaz, Leopoldo [Sandia National Laboratories (SNL-↗

DOE OSTI2023

Dynamic Rupture Simulations of the Mw7.2 1992 Landers,California, Earthquake

An important part of the simulation technique implemented in the Recipe is the earthquake rupture characterization based on established empirical scaling relationships between source parameters and seismic moment. The scaling relationships used in the rupture characterization are derived from recorded data and are tested against dynamic and kinematic rupture models of crustal earthquakes. Currently, a revised three-stage source scaling relationship between the source rupture area and seismic moment of inland crustal earthquakes, initially proposed by Irikura and Miyake, is being used. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, including the PI, Arben Pitarka, indicate that in order to reproduce the near-fault long-period ground motions and the permanent displacements, observed during surface-rupturing earthquakes, the recipe needs modifications to the kinematic rupture parametrization in the shallow crustal weak zone. In particular, improvements are required in characterizing the slip rate function at shallow depths, above the sesimogenic zone. The main objective of this research project was the development of a dynamic rupture model for the Mw7.2 1992 Landers, California earthquake. The dynamic rupture model was used to derive kinematic rupture parameters, such as slip, slip-rate, rupture time distributions, and slip rate functions in the weak (0-4km) and seismogenic crustal zones. In addition to the knowledge gained from the dynamic rupture modeling of Chi-chi, Ridgecrest, Kumamoto and other crustal earthquakes, studied by the PI, the rupture characterization of the Landers earthquake can be helpful in updating the conventional recipe for broadband ground motion simulations of crustal earthquakes with shallow rupture. The project emphasizes the need for performing physics-based earthquake rupture modeling which ultimately will inform further improvements of rupture characterization. In particular, the modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, used in adjustments of empirical source scaling relationships for large crustal earthquakes.

58 GEOSCIENCES↗

DOE OSTI2021

Multi-Fluid and Kinetic Models of Partially Ionized Magnetic Reconnection

Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomena in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet-Parker scaling, where neutrals are unimportant, or a fully coupled Sweet-Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfv\'en speed, $v_A^\star$. These differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $\nu_{ni}L/v_A^\star \gtrsim 1$, where $\nu_{ni}$ is the neutral-ion collision frequency and $L$ is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $\nu_{ni}L/v_A^\star < 1$, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.

Magnetic reconnection↗

DOE OSTI2021

StressNet - Deep learning to predict stress with fracture propagation in brittle materials

Abstract Catastrophic failure in brittle materials is often due to the rapid growth and coalescence of cracks aided by high internal stresses. Hence, accurate prediction of maximum internal stress is critical to predicting time to failure and improving the fracture resistance and reliability of materials. Existing high-fidelity methods, such as the Finite-Discrete Element Model (FDEM), are limited by their high computational cost. Therefore, to reduce computational cost while preserving accuracy, a deep learning model, StressNet, is proposed to predict the entire sequence of maximum internal stress based on fracture propagation and the initial stress data. More specifically, the Temporal Independent Convolutional Neural Network (TI-CNN) is designed to capture the spatial features of fractures like fracture path and spall regions, and the Bidirectional Long Short-term Memory (Bi-LSTM) Network is adapted to capture the temporal features. By fusing these features, the evolution in time of the maximum internal stress can be accurately predicted. Moreover, an adaptive loss function is designed by dynamically integrating the Mean Squared Error (MSE) and the Mean Absolute Percentage Error (MAPE), to reflect the fluctuations in maximum internal stress. After training, the proposed model is able to compute accurate multi-step predictions of maximum internal stress in approximately 20 seconds, as compared to the FDEM run time of 4 h, with an average MAPE of 2% relative to test data.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Fullerene-like defects in high-temperature neutron-irradiated nuclear graphite

Irradiation-induced defect evolution in graphite is particularly important for its application in graphite-moderated nuclear reactors. The evolution of defects directly influences macroscopically observed property changes in irradiated nuclear graphite which, in turn, can govern the lifetime of graphite components. This article reports novel defect structures and the irradiation response of microstructural features occurring in high-temperature irradiated nuclear graphite IG-110. High resolution transmission electron microscopy (HRTEM) was used to characterize specimens neutron-irradiated at a high temperature (≥ 800°C) at doses of 1.73 and 3.56 atomic displacements per atom (dpa). Concentric shelled and fullerene-like defects were found to result in swelling along the c-axis and contraction along the a/b-axis of crystallites. Furthermore, such defects are shown to occur within, and partially fill, Mrozowski cracks prior to turnaround dose. In addition, in situ TEM under similar irradiation conditions was used to capture the real-time dynamic evolution of defects, providing unambiguous analysis of the evolution of the graphite structures during irradiation. Here, the results suggest the mainstream theory for radiation damage in nuclear graphite (which assumes additional basal plane formation as the sole reason) to be an incorrect interpretation of defect evolution contributing to irradiation-induced property changes at higher temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

An introduction to Spent Nuclear Fuel decay heat for Light Water Reactors: a review from the NEA WPNCS

This paper summarized the efforts performed to understand decay heat estimation from existing spent nuclear fuel (SNF), under the auspices of the Working Party on Nuclear Criticality Safety (WPNCS) of the OECD Nuclear Energy Agency. Needs for precise estimations are related to safety, cost, and optimization of SNF handling, storage, and repository. The physical origins of decay heat (a more correct denomination would be decay power) are then introduced, to identify its main contributors (fission products and actinides) and time-dependent evolution. Due to limited absolute prediction capabilities, experimental information is crucial; measurement facilities and methods are then presented, highlighting both their relevance and our need for maintaining the unique current full-scale facility and developing new ones. The third part of this report is dedicated to the computational aspect of the decay heat estimation: calculation methods, codes, and validation. Different approaches and implementations currently exist for these three aspects, directly impacting our capabilities to predict decay heat and to inform decision-makers. Finally, recommendations from the expert community are proposed, potentially guiding future experimental and computational developments. One of the most important outcomes of this work is the consensus among participants on the need to reduce biases and uncertainties for the estimated SNF decay heat. If it is agreed that uncertainties (being one standard deviation) are on average small (less than a few percent), they still substantially impact various applications when one needs to consider up to three standard deviations, thus covering more than 95% of cases. The second main finding is the need of new decay heat measurements and validation for cases corresponding to more modern fuel characteristics: higher initial enrichment, higher average burnup, as well as shorter and longer cooling time. Similar needs exist for fuel types without public experimental data, such as MOX, VVER, or CANDU fuels. A third outcome is related to SNF assemblies for which no direct validation can be performed, representing the vast majority of cases (due to the large number of SNF assemblies currently stored, or too short or too long cooling periods of interest). A few solutions are possible, depending on the application. For the final repository, systematic measurements of quantities related to decay heat can be performed, such as neutron or gamma emission. This would provide indications of the SNF decay heat at the time of encapsulation. For other applications (short- or long-term cooling), the community would benefit from applying consistent and accepted recommendations on calculation methods, for both decay heat and uncertainties. This would improve the understanding of the results and make comparisons easier.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

Electrochemical Modeling of PID Leakage Current of PV Modules: Steady-State Current, Transient Current, and RC-Equivalent Circuit

Potential-induced degradation (PID) remains a significant reliability concern for photovoltaic (PV) modules, arising when a voltage difference between the module frame and the solar cells drives unintended leakage current through the glass-encapsulant stack. Although PID ultimately manifests as PID-s, PID-p, or PID-c, the underlying behavior of the leakage current-its magnitude and time dependence-requires clearer electrochemical interpretation. Traditional explanations attribute the initial transient current to bulk capacitive elements of the glass, encapsulant, and antireflection coatings, and the steady-state current according to their effective ohmic resistance. More recent studies, however, indicate that electrochemical charge-transfer processes at the encapsulant-metallization interface can play a dominant role in defining the leakage-current path. This paper develops a unified electrochemical framework for modeling PID leakage current. First, an RC-equivalent circuit is formulated by combining conventional RC elements with a Randles-type interface to capture transient leakage current through double-layer capacitance and faradaic processes at ionic-electronic boundaries. Second, the steady-state current-voltage behavior is explained using a linearized Butler-Volmer relationship, showing that the measured ohmic response corresponds to the low-overpotential limit of charge-transfer kinetics. Analytical results demonstrate that, for typical module materials-3.2-mm soda-lime glass and 0.45-mm encapsulant-the dominant modulators to PID leakage current are the glass surface resistance (under dry-surface conditions), the glass bulk capacitance, and the encapsulant resistance (under wet-surface conditions), with soda lime glass surface and EVA/POE encapsulant resistances primarily governing steady-state current. The proposed electrochemical model is validated against measured leakage-current data, showing good agreement in both the magnitude and the time-dependent evolution of PID leakage current.

14 SOLAR ENERGY↗