Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parallel processing”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Assessment of ESM Readiness Level for Exascale HPC

Advancement of Earth System Models (ESMs) is becoming increasingly challenging due to a confluence of factors including increasing model complexity – to more fully represent the earth system, increasing spatial resolution - to achieve higher accuracy by resolving fine-scale dynamical to physical, biological, and chemical processes and their interaction, increasing ensemble size - to more accurately represent predictive uncertainty, and increased computing requirements – to enable more accurate and timely weather predictions and climate projections for societal benefit. The belief by many that computing will take care of itself is no longer valid given the disruptive changes in HPC that are driving up the cost of computing, increasing the difficulty of using emerging HPC effectively, and exposing limits in parallelism, portability and scalability of the ESM applications themselves.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Machine learning-based optimization of air-cooled heat sinks

Machine learning-based models using Artificial Neural Network (ANN) and greedy search algorithm are used to optimize air-cooled parallel plate-finned heat sinks (PPFHSs) subjected to laminar flow over an extensive range of design parameters. Here, the thermal and hydraulic performances of PPFHSs are represented by heat transfer coefficient (h) and pressure drop (ΔP), respectively. Optimization objectives for PPFHS designs can vary from industry to industry depending on their design priorities. The present study proposes a novel and generalized optimization method that defines practical optimization objectives and provides an accurate optimization process to design effective PPFHSs for a wide range of industrial applications with different design requirements. Three optimization objectives are presented in this study: (i) the largest h PΔ, (ii) the largest h within a specified maximum allowed flow rate, and (iii) the lowest weight that maximizes h for operation within the maximum allowed flow rate. While the shortcoming of the first objective is demonstrated, the other two objectives are found to be suitable for designing effective heat sinks (HSs) across different applications. Results suggest a promising trend from the third objective to develop HSs with ~ 37-68% lower weight, 80-85% reduced ΔP, and negligible penalty in h compared with optimized HSs obtained from the second objective. However, since the third objective leads to HSs with thinner fins, structural analysis should be performed to ensure reliable operation of the HSs.

42 ENGINEERING↗
DOE OSTI2024

Study of magnetic reconnection at low- β using laser-powered capacitor coils

Magnetic reconnection is a ubiquitous fundamental process in space and astrophysical plasmas that rapidly converts magnetic energy into some combination of flow energy, thermal energy, and non-thermal energetic particles. Over the past decade, a new experimental platform has been developed to study magnetic reconnection using strong coil currents powered by high-power lasers at low plasma beta, typical conditions under which reconnection is energetically important in space and astrophysics. KJ-class lasers were used to drive parallel currents to reconnect MG-level magnetic fields in a quasi-axisymmetric geometry, similar to the magnetic reconnection experiment or MRX, and thus this platform is named micro-MRX. This presentation summarizes two major findings from micro-MRX: direct measurement of accelerated electrons and observation of ion acoustic waves during anti-parallel reconnection. The angular dependence of the measured electron energy spectrum and the resulting accelerated energies, supported by particle-in-cell simulations, indicate that direct acceleration by the out-of-plane reconnection electric field is at work. Furthermore, a sudden onset of ion acoustic bursts has been measured by collective Thomson scattering in the exhaust of magnetic reconnection, followed by electron acoustic bursts with electron heating and bulk acceleration. These results demonstrate that the micro-MRX platform offers a novel and unique approach to study magnetic reconnection in the laboratory in addition to the capabilities provided by traditional magnetized plasma experiments such as MRX and the upcoming Facility for Laboratory Reconnection experiments (FLARE). Future prospects to study other particle acceleration mechanisms and ion acoustic waves from magnetic reconnection are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Online multimedia retrieval on CPU–GPU platforms with adaptive work partition

Nearest neighbors search is a core operation found in several online multimedia services. These services have to handle very large databases, while, at the same time, they must minimize the query response times observed by users. This is specially complex because those services deal with fluctuating query workloads (rates). Consequently, they must adapt at run-time to minimize the response times as the load varies. In this paper, we address the aforementioned challenges with a distributed memory parallelization of the product quantization nearest neighbor search, also known as IVFADC, for hybrid CPU–GPU machines. Overall, our parallel IVFADC implements an out-of-GPU memory execution scheme to use the GPU for databases in which the index does not fit in its memory, which is crucial for searching in very large databases. The careful use of CPU and GPU with work stealing led to an average response time reduction of 2.4 as compared to using the GPU only. Also, our approach to adapt the system to fluctuating loads, called Dynamic Query Processing Policy (DQPP), attained a response time reduction of up to 5 vs. the best static (BS) policy for moderate loads. The system has attained high query processing rates and near-linear scalability in all experiments. We have evaluated our system on a machine with up to 256 NVIDIA V100 GPUs processing a database of 256 billion SIFT features vectors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Deformation behavior of nanoscale Al–Al 2 Cu eutectics studied by in situ micropillar compression

We report deformation behavior of nanoscale laser processed Al–Al 2 Cu eutectics at room temperature is characterized through in situ micro-pillar compression testing in a scanning transmission microscope. Interlamellar spacing of Al–Al 2 Cu eutectics varies from hundreds of nanometers to 20 nm. Three different sizes of micro-pillars are fabricated in order to study the deformation behaviors of single colony and multiple colonies, corresponding to the single crystal and polycrystal respectively. For single colonies, lamellar orientations parallel, normal or inclined to the loading direction were tested. The main findings are: 1) the plasticity mechanisms strongly depend on loading orientation: buckling and kinking in the parallel-loaded eutectics, planar sliding along Al–Al 2 Cu lamellar interfaces in the incline-loaded eutectics and localized shearing in the normal-loaded eutectics. 2) the incline-loaded eutectics exhibits the lowest compression flow strength, and the normal-loaded eutectics has the highest compression flow strength. 3) with decreasing inter-lamellar spacing, the strength increases and plasticity is uniformly distributed, as opposed to shear localization. Highest compressive plasticity is observed in polycrystalline eutectics with an ensemble of lamellar orientations with ~20 nm average spacing: 17.9% at flow stress of 1.63 GPa, and degenerate, bimodal morphology: 11.1% at flow stress of 1.36 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2020

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical sciences by providing tools to support and guide the experimental effort and for prediction of chemical and materials properties. In this regard, a special role has been played by electronic structure packages where complex chemical and materials processes can be modeled using first-principle-driven methodologies. Over the last few decades, the rapid development of computing technologies and a tremendous increase in computational power has offered a unique chance to study complex chemical transformations using sophisticated and predictive many-body techniques to describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, a critical role has been played by novel parallel algorithms capable of taking advantage of computational resources to address polynomial scaling of electronic structure methods. NWChem was among the first electronic structure codes that focused on delivering scalable performance for electronic structure simulations. Herein, we briefly review the NWChem suite of computational codes including its history, design principles, parallel tools, current capabilities, outreach and outlook.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

VWC-BERT: Scaling Vulnerability–Weakness–Exploit Mapping on Modern AI Accelerators

Defending cybersystems needs accurate mapping of software and hardware vulnerabilities to generalized descriptions of weaknesses, and weaknesses to exploits. These mappings enable cyber defenders to build plans for effective defense and assessment of potential risks to a cybersystem. With close to 170k vulnerabilities, manual mapping is not a feasible option. However, automated mapping is challenging due to limited training data, computational intractability, and limitations in computational natural language processing. Tools based on breakthroughs in Transformer-based language models have been demonstrated to classify vulnerabilities with high accuracy. We make three key contributions in this paper: (1) We present a new framework, \VWCBERT, that augments the Transformer-based hierarchical multi-class classification framework of Das et al. (\textsc{V2W-BERT}) with the ability to map weaknesses to exploits. (2) We implement \VWCBERT~ on modern AI accelerator platforms using two data parallel techniques for the pre-training phase and demonstrate nearly linear speedups across NVIDIA and Graphcore accelerator platforms. We observe nearly linear speedups for up to 16 V100 and 8 A100 GPUs, and about 3.4$\times$ speedup for A100 relative to V100 GPUs. We also observe excellent speedups on Graphcore, with $5.7\times$ speedup on 128 IPUs relative to 16 IPUs. Enabled by scaling, we also demonstrate higher accuracy using a larger language model, RoBERTa-Large. We show up to 87\% accuracy for strict and up to 98\% accuracy for relaxed classification. (3) We develop a novel parallel link manager for the link prediction phase and demonstrate up to 21$\times$ speedup with 16 V100 GPUs relative to one V100 GPU, and thus reducing the runtime from 2.5 hours to 10 minutes. We believe that generalizability and scalability of \VWCBERT~ will benefit both the theoretical development and practical deployment of novel cyberdefense solutions and vulnerability classification.

Das, Siddhartha Shankar↗
DOE OSTI2024

Evolution of the SLATE linear algebra library

SLATE (Software for Linear Algebra Targeting Exascale) is a distributed, dense linear algebra library targeting both CPU-only and GPU-accelerated systems, developed over the course of the Exascale Computing Project (ECP). While it began with several documents setting out its initial design, significant design changes occurred throughout its development. In some cases, these were anticipated: an early version used a simple consistency flag that was later replaced with a full-featured consistency protocol. In other cases, performance limitations and software and hardware changes prompted a redesign. Sequential communication tasks were parallelized; host-to-host MPI calls were replaced with GPU device-to-device MPI calls; more advanced algorithms such as Communication Avoiding LU and the Random Butterfly Transform (RBT) were introduced. Early choices that turned out to be cumbersome, error prone, or inflexible have been replaced with simpler, more intuitive, or more flexible designs. Applications have been a driving force, prompting a lighter weight queue class, nonuniform tile sizes, and more flexible MPI process grids. Of paramount importance has been building a portable library that works across several different GPU architectures – AMD, Intel, and NVIDIA – while keeping a clean and maintainable codebase. Here we explore the evolving design choices and their effects, both in terms of performance and software sustainability.

Gates, Mark↗
DOE OSTI2022

Cobalt-Free Cathodes for Next Generation Li-Ion Batteries

In this U.S. Department of Energy sponsored project Nexceris, in collaboration with project partners; The Ohio State University and Navitas Advanced Systems have advanced the technical maturity of a non-cobalt containing cathode for next-generation Li-ion batteries. The cathode is based on the lithium manganese nickel-titanium oxide, LiNi 0.5 Mn 1.5 TiO 4 (LNMTO) high voltage spinel. To address limitations with poor cycle and calendar life an microstucturally hierarchical LNMO/LNMTO core-shell cathode powder has been developed that enables the formation of a solid-electrolyte interface that effectively passivates the cathode surface. The microstructural enhancements of the cathode material focus on preferentially enriching the surface with titanium. In parallel, new, optimized binder and electrolyte chemistries have been incorporated to address degradation mechanisms associated with high-voltage systems. Single-layer pouch cell and large-format 2-Ah cell testing have shown that an optimized LNMO/LNMTO core-shell powder significantly improves initial cell capacity and cycle life compared to homogeneous LNMO powder. To support development and the fabrication of 2-Ah cells a novel Hybrid Alternative Wet-Chemical Synthesis (HAWCS) process have been developed. This low-cost, synthesis approach enables the excellent compositional and particle morphology control achieved with co-precipitation without the strict process controls and associated expensive process equipment.

25 ENERGY STORAGE↗
DOE OSTI2022

Heat conduction in an irregular magnetic field. Part 2. Heat transport as a measure of the effective non-integrable volume

Given the large anisotropy of transport processes in magnetized plasmas, the magnetic field structure can strongly impact heat diffusion: magnetic surfaces and cantori form barriers to transport while chaotic layers and island structures can degrade confinement. When a small but non-zero amount of perpendicular diffusion is included, the structure of the magnetic field becomes less important, allowing pressure gradients to be supported across chaotic regions and island chains. We introduce a metric for the effective volume over which the local parallel diffusion dominates based on the solution to the anisotropic heat diffusion equation. To validate this metric, we consider model fields with a single island chain and a strongly chaotic layer for which analytic predictions of the relative parallel and perpendicular transport can be made. We also analyse critically chaotic fields produced from different sets of perturbations, highlighting the impact of the mode number spectrum on the heat transport. Our results indicate that this metric coincides with the effective volume of non-integrability in the limit $\kappa _{\perp } \rightarrow 0$ , where $\kappa_{\perp}$ is the perpendicular diffusion coefficient. We propose that this metric be used to assess the impact of non-integrability on the heat transport in stellarator equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Porting Fragmentation Methods to Graphical Processing Units Using an OpenMP Application Programming Interface: Offloading the Fock Build for Low Angular Momentum Functions

Here, a framework to offload four-index two-electron repulsion integrals to graphical processing units (GPUs) using OpenMP is discussed. The method has been applied to the Fock build for low angular momentum s and p functions in both the restricted Hartree–Fock (RHF) and in the effective fragment molecular orbital (EFMO) framework. Benchmark calculations for the GPU code for the pure RHF method show an increasing speedup relative to the existing OpenMP CPU code in GAMESS from 1.04 to 52× for clusters of 70–569 water molecules. The parallel efficiency on 24 NVIDIA V100 GPU boards also increases when increasing the system size: from 75 to 94% for water clusters that contain 303–1120 molecules. In the EFMO framework, the GPU Fock build shows a high linear scalability up to 4608 V100s with a parallel efficiency of 96% for calculations on a solvated mesoporous silica nanoparticle system with ~67,000 basis functions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Overview of physics results from MAST upgrade towards core-pedestal-exhaust integration

Recent results from MAST Upgrade are presented, emphasising understanding the capabilities of this new device and deepening understanding of key physics issues for the operation of ITER and the design of future fusion power plants. The impact of MHD instabilities on fast ion confinement have been studied, including the first observation of fast ion losses correlated with Compressional and Global Alfvén Eigenmodes. High-performance plasma scenarios have been developed by tailoring the early plasma current ramp phase to avoid internal reconnection events, resulting in a more monotonic q profile with low central shear. The impact of m/n = 3/2, 2/1 and 1/1 modes on thermal plasma confinement and rotation profiles has been quantified, and scenarios optimised to avoid them have transiently reached values of normalised beta approaching 4.2. In pedestal and ELM physics, a maximum pedestal top temperature of ~350 eV has been achieved, exceeding the value achieved on MAST at similar heating power. Mitigation of type-I ELMs with n = 1 RMPs has been observed. Studies of plasma exhaust have concentrated on comparing conventional and Super-X divertor configurations, while X-point target, X-divertor and snowflake configurations have been developed and studied in parallel. In L-mode discharges, the separatrix density required to detach the outer divertors is approximately a factor 2 lower in the Super-X than the conventional configuration, in agreement with simulations. Detailed analysis of spectroscopy data from studies of the Super-X configuration reveal the importance of including plasma-molecule interactions and D 2 Fulcher band emission to properly quantify the rates of ionisation, plasma-molecule interactions and volumetric recombination processes governing divertor detachment. In H-mode with conventional and Super-X configurations, the outer divertors are attached in the former and detached in the latter with no impact on core or pedestal confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Nonlinear optical control of chiral charge pumping in a topological Weyl semimetal

Solids with topologically robust electronic states exhibit unusual electronic and optical properties that do not exist in other materials. A particularly interesting example is chiral charge pumping, the so-called chiral anomaly, in recently discovered topological Weyl semimetals, where simultaneous application of parallel DC electric and magnetic fields creates an imbalance in the number of carriers of opposite topological charge (chirality). Here, using time-resolved terahertz measurements on the Weyl semimetal TaAs in a magnetic field, we optically interrogate the chiral anomaly by dynamically pumping the chiral charges and monitoring their subsequent relaxation of the nonequilibrium state. Theory based on Boltzmann transport shows that the observed effects originate from an optical nonlinearity in the chiral charge pumping process. Our measurements reveal that the nonequilibrium chiral excitation relaxation time is much greater than 1 ns. The observation of terahertz-controlled chiral carriers with long coherence times and topological protection suggests the application of Weyl semimetals for quantum optoelectronic technology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Feasibility of Lead Fast Reactor Heat Exchanger Tube Online Monitoring

Within virtually all reactor systems with heat exchangers, the thinnest barriers between primary and secondary fluids are the heat exchanger tubes. Consistent with light water reactor experience, the advanced reactor heat exchangers are the most susceptible regions for corrosion and leakage (NRC, 2019). Moreover, traditional inspections and repairs using robotics will be expensive or impractical. On-line structural health corrosion monitoring (OLSHM) is therefore desirable to detect tube degradation prior to leaks that may allow mixing of heat exchanger fluids. Rhetorically many molten salt reactor (MSR) designers plan for on-line monitoring of critical reactor components including heat exchangers but there are no proven OLSHM sensors or systems for MSR heat exchangers. Guided wave ultrasound from the tube-end or ID using conventional low-temperature piezoelectric sensors has been shown to be sensitive to pits and cracks less than 50% through-wall in long pipes and tubes (Narayanan M.M., 2019) . Mode and frequency selection are important for practical sensor design decisions and ultimate performance. The L-0 wave-mode is easily excited by L-mode transducers but they have surface normal wave motion that will be sensitive to and more attenuative because of fluid interaction with the tube surface (Rose J.L., 1994). Shear Horizontal (SH) wave modes have material motion parallel to the tube surface therefore are virtually unaffected by fluid interaction. For SH-0 to be nondispersive in steel, the frequency-thickness product (MHz-mm) must be less than approximately 1.5. SH mode transducers must be coupled by a material that can support shear stress; a thin layer of honey works well for temporary sensors and an adhesive or a brazed metallurgical bond for permanent sensors. If sensors are incorporated into the fabrication process, the same kinds of SH-0 waves can be generated from the stagnant flow area of the tube OD without compromising the heat-exchange tube function. Piezoelectric materials are also available that can withstand harsh temperatures of molten salt heat exchangers (Stevenson T., 2015) (Machura et al., 2008). Though they may not have as high-performance characteristics, their performance only needs to be good enough to detect flaws of interest. This project attempts to demonstrate feasibility to develop a sensor that can be sensitive to flaws of interest and perform long-term monitoring of a molten salt reactor heat exchanger. Benefits to the advanced reactor system include: • Minimizing the expense of a reactor shut down for periodic inspection—only shutdown for cause. • Reduce risk of unexpected tube rupture by detecting flaws before they grow to through-wall leaks. • Cost reduction by minimizing need to design for manual/robotic inspectability. Transducer design, characterization of the beam profile, wedges vs. direct application of piezoelectric sensors, and ultrasonic pulse-echo responses from a representative heat exchange tube with fabricated flaw types are discussed.

On-Line Tube Monitor, Guided Wave Ultrasound, SH-0↗
DOE OSTI2020

Progress Toward Molten Salt Reactor Heat Exchanger On-Line Monitoring

Within virtually all reactor systems with heat exchangers, the thinnest barriers between primary and secondary fluids are the heat exchanger tubes. Consistent with light water reactor experience, the advanced reactor heat exchangers are the most susceptible regions for corrosion and leakage (NRC, 2019). Moreover, traditional inspections and repairs using robotics will be expensive or impractical. On-line structural health corrosion monitoring (OLSHM) is therefore desirable to detect tube degradation prior to leaks that may allow mixing of heat exchanger fluids. Rhetorically many molten salt reactor (MSR) designers plan for on-line monitoring of critical reactor components including heat exchangers but there are no proven OLSHM sensors or systems for MSR heat exchangers. Guided wave ultrasound from the tube-end or ID using conventional low-temperature piezoelectric sensors has been shown to be sensitive to pits and cracks less than 50% through-wall in long pipes and tubes (Narayanan M.M., 2019) . Mode and frequency selection are important for practical sensor design decisions and ultimate performance. The L-0 wave-mode is easily excited by L-mode transducers but they have surface normal wave motion that will be sensitive to and more attenuative because of fluid interaction with the tube surface (Rose J.L., 1994). Shear Horizontal (SH) wave modes have material motion parallel to the tube surface therefore are virtually unaffected by fluid interaction. For SH-0 to be nondispersive in steel, the frequency-thickness product (MHz-mm) must be less than approximately 1.5. SH mode transducers must be coupled by a material that can support shear stress; a thin layer of honey works well for temporary sensors and an adhesive or a brazed metallurgical bond for permanent sensors. If sensors are incorporated into the fabrication process, the same kinds of SH-0 waves can be generated from the stagnant flow area of the tube OD without compromising the heat-exchange tube function. Piezoelectric materials are also available that can withstand harsh temperatures of molten salt heat exchangers (Stevenson T., 2015) (Machura et al., 2008). Though they may not have as high-performance characteristics, their performance only needs to be good enough to detect flaws of interest. This project attempts to demonstrate feasibility to develop a sensor that can be sensitive to flaws of interest and perform long-term monitoring of a molten salt reactor heat exchanger. Benefits to the advanced reactor system include: • Minimizing the expense of a reactor shut down for periodic inspection—only shutdown for cause. • Reduce risk of unexpected tube rupture by detecting flaws before they grow to through-wall leaks. • Cost reduction by minimizing need to design for manual/robotic inspectability. Transducer design, characterization of the beam profile, and ultrasonic pulse-echo responses from a representative heat exchange tube with fabricated flaw types are discussed.

On-Line Tube Monitor, Guided Wave Ultrasound, SH-0↗
DOE OSTI2023

Data and scripts associated with a manuscript investigating impacts of solid phase extraction on freshwater organic matter optical signatures and mass spectrometry pairing

This data package is associated with the publication “Investigating the impacts of solid phase extraction on dissolved organic matter optical signatures and the pairing with high-resolution mass spectrometry data in a freshwater system” submitted to “Limnology and Oceanography: Methods.” This data is an extension of the River Corridor and Watershed Biogeochemistry SFA’s Spatial Study 2021 (https://doi.org/10.15485/1898914). Other associated data and field metadata can be found at the link provided. The goal of this manuscript is to assess the impact of solid phase extraction (SPE) on the ability to pair ultra-high resolution mass spectrometry data collected from SPE extracts with optical properties collected on ambient stream samples. Forty-seven samples collected from within the Yakima River Basin, Washington were analyzed dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC), absorbance, and fluorescence. Samples were subsequently concentrated with SPE and reanalyzed for each measurement. The extraction efficiency for the DOC and common optical indices were calculated. In addition, SPE samples were subject to ultra-high resolution mass spectrometry and compared with the ambient and SPE generated optical data. Finally, in addition to this cross-platform inter-comparison, we further performed and intra-comparison among the high-resolution mass spectrometry data to determine the impact of sample preparation on the interpretability of results. Here, the SPE samples were prepared at 40 milligrams per liter (mg/L) based on the known DOC extraction efficiency of the samples (ranging from ~30 to ~75%) compared to the common practice of assuming the DOC extraction efficiency of freshwater samples at 60%. This data package folder consists of one main data folder with one subfolder (Data_Input). The main data folder contains (1) readme; (2) data dictionary (dd); (3) file-level metadata (flmd); (4) final data summary output from processing script; and (5) the processing script. The R-markdown processing script (SPE_Manuscript_Rmarkdown_Data_Package.rmd) contains all code needed to reproduce manuscript statistics and figures (with the exception of that stated below). The Data_Input folder has two subfolders: (1) FTICR and (2) Optics. Additionally, the Data_Input folder contains dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) data (SPS_NPOC_Summary.csv) and relevant supporting Solid Phase Extraction Volume information (SPS_SPE_Volumes.csv). Methods information for the optical and FTICR data is embedded in the header rows of SPS_EEMs_Methods.csv and SPS_FTICR_Methods.csv, respectively. In addition, the data dictionary (SPS_SPE_dd.csv), file level metadata (SPS_SPE_flmd.csv), and methods codes (SPS_SPE_Methods_codes.csv) are provided. The FTICR subfolder contains all raw FTICR data as well as instructions for processing. In addition, post processed FTICR molecular information (Processed_FTICRMS_Mol.csv) and sample data (Processed_FTICRMS_Data.csv) is provided that can be directly read into R with the associated R-markdown file. The Optics subfolder contains all Absorbance and Fluorescence Spectra. Fluorescence spectra have been blank corrected, inner filter corrected, and undergone scatter removal. In addition, this folder contains Matlab code used to make a portion of Figure 1 within the manuscript, derive various spectral parameters used within the manuscript, and used for parallel factor analysis (PARAFAC) modeling. Spectral indices (SPS_SpectralIndices.csv) and PARAFAC outputs (SPS_PARAFAC_Model_Loadings.csv and SPS_PARAFAC_Sample_Scores.csv) are directly read into the associated R-markdown file. We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected some of these data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Equilipy: a python package for calculating phase equilibria

The CALPHAD (CALculation of PHAse Diagram) approach (Nigel Saunders & Miodownik, 1998) provides predictions for thermodynamically stable phases in multicomponent-multiphase materials across a wide range of temperatures. Consequently, the CALPHAD calculations became an essential tool in materials and process design (Luo, 2015). Such design tasks frequently require navigating a high-dimensional space due to multiple components involved in the system. This increasing complexity demands high-throughput CALPHAD calculations, especially in the rapidly evolving field of alloy design. In response to the need, we developed Equilipy an open-source Python package designed for calculating phase equilibria of multicomponent-multiphase systems. Equilipy is specifically tailored for high-throughput CALPHAD calculations, offering parallel computations across multiple processors and nodes with the given NPT input conditions namely elemental compositions (N), pressure (P), and temperature (T). Equilipy utilizes the program structure and Gibbs energy functions from the Fortran-based program, Thermochimica (Piro et al., 2013), with incorporating a new Gibbs energy minimization algorithm. This algorithm, originally developed by Capitani and Brown in 1987 (Capitani & Brown, 1987), has been revised and implemented to enhance the stability and performance of calculations. The Fortran codes are precompiled and interfaced with Python via F2PY, ensuring high computation speed. Benchmark tests shown in Figure 1 demonstrate that Equilipy’s computation speed is comparable to those of established commercial software, TC-Python and PanPython. This result highlights its efficiency and potential applications in various scientific and industrial fields.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗