SEARCH · Search NASA

Results for “EM simulation”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗

DOE OSTI2024

Hub stability in the calcium calmodulin-dependent protein kinase II

The calcium calmodulin protein kinase II (CaMKII) is a multi-subunit ring assembly with a central hub formed by the association domains. There is evidence for hub polymorphism between and within CaMKII isoforms, but the link between polymorphism and subunit exchange has not been resolved. Here, we present near-atomic resolution cryogenic electron microscopy (cryo-EM) structures revealing that hubs from the α and β isoforms, either standalone or within an β holoenzyme, coexist as 12 and 14 subunit assemblies. Single-molecule fluorescence microscopy of Venus-tagged holoenzymes detects intermediate assemblies and progressive dimer loss due to intrinsic holoenzyme lability, and holoenzyme disassembly into dimers upon mutagenesis of a conserved inter-domain contact. Molecular dynamics (MD) simulations show the flexibility of 4-subunit precursors, extracted in-silico from the β hub polymorphs, encompassing the curvature of both polymorphs. The MD explains how an open hub structure also obtained from the β holoenzyme sample could be created by dimer loss and analysis of its cryo-EM dataset reveals how the gap could open further. An assembly model, considering dimer concentration dependence and strain differences between polymorphs, proposes a mechanism for intrinsic hub lability to fine-tune the stoichiometry of αβ heterooligomers for their dynamic localization within synapses in neurons.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

EM Physics

Geant4 provides a comprehensive set of electromagnetic (EM) processes and models for electron/positron, gamma and long-lived charged particles, spanning energies from 100 eV to 100 TeV. Covering diverse energy regions often requires multiple models, which can be constructed using pre-packaged or user-defined EM physics constructors. Geant4 also supports detailed low-energy EM physics through models like Livermore, Penelope, and ICRU73, offering extensive data for elements across a wide energy range (250 eV–100 GeV). Application domains include space, medical, and radiobiology simulations. Additionally, Geant4 provides robust options for simulating complex optical photon production and transportation processes, enhancing its versatility in physics research and applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

Rethinking materials simulations: Blending direct numerical simulations with neural operators

Abstract Materials simulations based on direct numerical solvers are accurate but computationally expensive for predicting materials evolution across length- and time-scales, due to the complexity of the underlying evolution equations, the nature of multiscale spatiotemporal interactions, and the need to reach long-time integration. We develop a method that blends direct numerical solvers with neural operators to accelerate such simulations. This methodology is based on the integration of a community numerical solver with a U-Net neural operator, enhanced by a temporal-conditioning mechanism to enable accurate extrapolation and efficient time-to-solution predictions of the dynamics. We demonstrate the effectiveness of this hybrid framework on simulations of microstructure evolution via the phase-field method. Such simulations exhibit high spatial gradients and the co-evolution of different material phases with simultaneous slow and fast materials dynamics. We establish accurate extrapolation of the coupled solver with large speed-up compared to DNS depending on the hybrid strategy utilized. This methodology is generalizable to a broad range of materials simulations, from solid mechanics to fluid dynamics, geophysics, climate, and more.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Dynamics and activation of membrane-bound B cell receptor assembly

B-cell receptor (BCR) complexes are expressed on the surface of a B-cell and are critical in antigen recognition and modulating the adaptive immune response. Even though the relevance of antibodies has been known for almost a hundred years, the antigen-dependent activation mechanism of B-cells has remained elusive. Several models have been proposed for BCR activation, including cross-linking, conformation-induced oligomerization, and dissociation activation models. Recently, the first cryo-EM structures of the human B-cell antigen receptor of the IgM and IgG isotypes have been published that validates the asymmetric organization of the BCR complex. Here, we carry out extensive molecular dynamics simulations to probe the conformational changes upon antigen binding and the influence of the membrane lipids. We identify two critical dynamical events that could be associated with antigen-dependent activation of BCR. First, antigen binding causes increased flexibility in regions distal to the antigen binding site. Second, antigen binding alters the rearrangement of IgM transmembrane helices, including the relative interaction of Igα/Igβ that mediates intracellular signaling. Furthermore, these transmembrane rearrangements lead to changes in localized lipid composition. Our work indirectly supports the conformational-change induced models of BCR activation and contributes to the understanding of the antigen-dependent activation mechanism of BCRs.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Multiphysics modeling of accelerators through code integration

This work aims to improve the ability of particle accelerator researchers to develop high-performance accelerator cavity designs by creating an overall multiphysics framework that integrates and couples existing application codes. This framework will allow accelerator researchers to build multiphysics models that will optimize cavity design, improve understanding of whole-device performance, and reduce the development and fabrication costs of accelerator research. We utilize the open-source VizSchema data standard as an intermediate data structure interface layer to standardize interfaces between individual application codes. VizScema is extensively documented online, and plugins for VizSchema are available for popular visualization packages, including VisIt and ParaView. Currently, the work focuses on coupling the EM field solver COMSOL and the electron gun code MICHELLE to allow COMSOL field-solve results to be seamlessly used by MICHELLE for particle-solve. Later work will extend this integration to include other fields, particles, and thermodynamics simulation codes.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

NeuroSEM: A hybrid framework for simulating multiphysics problems by coupling PINNs and spectral elements

Multiphysics problems that are characterized by complex interactions among fluid dynamics, heat transfer, structural mechanics, and electromagnetics, are inherently challenging due to their coupled nature. While experimental data on certain state variables may be available, integrating these data with numerical solvers remains a significant challenge. Physics-informed neural networks (PINNs) have shown promising results in various engineering disciplines, particularly in handling noisy data and solving inverse problems in partial differential equations (PDEs). However, their effectiveness in forecasting nonlinear phenomena in multiphysics regimes, particularly involving turbulence, is yet to be fully established. Here, this study introduces NeuroSEM, a hybrid framework integrating PINNs with the highfidelity Spectral Element Method (SEM) solver, Nektar++. NeuroSEM leverages the strengths of both PINNs and SEM, providing robust solutions for multiphysics problems. PINNs are trained to assimilate data and model physical phenomena in specific subdomains, which are then integrated into the Nektar++ solver. We demonstrate the efficiency and accuracy of NeuroSEM for thermal convection in cavity flow and flow past a cylinder. The framework effectively handles data assimilation by addressing those subdomains and state variables where the data is available. We applied NeuroSEM to the Rayleigh-B´enard convection system, including cases with missing thermal boundary conditions and noisy datasets. Finally, we applied the proposed NeuroSEM framework to real particle image velocimetry (PIV) data to capture flow patterns characterized by horseshoe vortical structures. Our results indicate that NeuroSEM accurately models the physical phenomena and assimilates the data within the specified subdomains. The framework’s plug-and-play nature facilitates its extension to other multiphysics or multiscale problems. Furthermore, NeuroSEM is optimized for efficient execution on emerging integrated GPU-CPU architectures. This hybrid approach enhances the accuracy and efficiency of simulations, making it a powerful tool for tackling complex engineering challenges in various scientific domains.

42 ENGINEERING↗

DOE OSTI2025

Development and Evaluation of a General Drag Model for Gas-Solid Flows via Deep Learning

This project presents the development and evaluation of a general drag model for gas–solid multiphase flows using deep learning techniques. A comprehensive database of more than 4,000 experimental and numerical data points for spherical and non spherical particles was compiled, incorporating geometric features such as sphericity, aspect ratio, and orientation. Several predictive approaches—including traditional em pirical correlations, machine learning, and deep neural networks—were benchmarked, with the proposed Drag Coefficient Correlation-aided Deep Neural Network (DCC DNN) demonstrating superior accuracy. To account for particle–particle interactions, additional drag data were generated using CFD-based simulations of packed and flu idized beds, leading to the development of a retrained model capable of incorporat ing volume fraction effects. Integration of the trained model with the MFiX CFD solver was achieved using FTorch, enabling drag predictions during discrete element method (DEM) simulations. Validation against experimental data for single particles and fluidized beds confirmed the model’s improved predictive ability, particularly for non-spherical geometries. While the model performed strongly under fluidized con ditions, limitations remained in unfluidized regimes, suggesting a need for expanded datasets. Overall, this study demonstrates the feasibility of combining deep learning with physics-informed CFD to improve drag modeling for gas–solid flows, with promis ing implications for scaling multiphase simulations in industrial applications.

42 ENGINEERING↗

DOE OSTI2024

Molecular architecture and functional dynamics of the pre-incision complex in nucleotide excision repair

Nucleotide excision repair (NER) is vital for genome integrity. Yet, our understanding of the complex NER protein machinery remains incomplete. Combining cryo-EM and XL-MS data with AlphaFold2 predictions, we build an integrative model of the NER pre-incision complex(PInC). Here TFIIH serves as a molecular ruler, defining the DNA bubble size and precisely positioning the XPG and XPF nucleases for incision. Using simulations and graph theoretical analyses, we unveil PInC’s assembly, global motions, and partitioning into dynamic communities. Remarkably, XPG caps XPD’s DNA-binding groove and bridges both junctions of the DNA bubble, suggesting a novel coordination mechanism of PInC’s dual incision. XPA rigging interlaces XPF/ERCC1 with RPA, XPD, XPB, and 5' ssDNA, exposing XPA’s crucial role in licensing the XPF/ERCC1 incision. Mapping disease mutations onto our models reveals clustering into distinct mechanistic classes, elucidating xeroderma pigmentosum and Cockayne syndrome disease etiology.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Neutrons in Structural Biology: Challenges and Opportunities (Workshop Report)

Gaining a thorough understanding of biological systems requires building our knowledge about biological processes from the level of atoms and electrons, and up to whole organisms. Such comprehensive knowledge will allow for a predictive understanding of complex biological systems behavior. It will guide us in the design and development of novel therapeutics and vaccines to tackle existing health threats and to prepare for future pandemics, and it will provide information necessary to create new biomaterials and bio-inspired technologies through manipulation of biological macromolecules, their assemblies, single cells and even microorganisms. Reaching these goals will require a synergistic combination of multiple experimental techniques with molecular calculations and predictive simulations, and the design and development of new techniques and capabilities that bridge current knowledge and technology gaps. Neutron scattering provides unique information about the biomacromolecular structure and function and can play a major role in achieving these goals. A workshop was held to engage the scientific community in identifying pressing challenges in biochemistry, structural biology, enzymology and structure-guided drug design not solved with the current neutron scattering technologies or utilizing other structural biology techniques such as X-ray crystallography, NMR, and cryo-EM. The workshop brought together structural biology, biochemistry and computational experts, as well as early career researchers and students, creating a forum for discussing scientific advancement and collaboration. The workshop included a one-day satellite training workshop where graduate students and postdoctoral researchers were educated in the application of neutron crystallography and small-angle scattering in structural biology. Furthermore, the Instrument Scientific Advisory Board (ISAB) for the development of a macromolecular neutron diffractometer at ORNL’s Second Target Station was introduced at the workshop. The major outcome was that neutrons can provide atomic-level understanding of biomacromolecular structure, function and dynamics which is of paramount importance for addressing the identified challenges. Neutron crystallography, in particular, can resolve long-standing biochemical issues regarding enzyme function by delineating the underlying chemistry and can have a major impact on the design of small-molecule therapeutics, especially in combination with molecular computation (quantum chemistry and molecular dynamics simulations) and the emerging artificial intelligence (AI)-assisted drug design technologies. The unique properties of neutrons, including their high sensitivity to hydrogen and their non-destructive nature, make them ideal probes of biological matter. There is a palpable need in the scientific community to expand and enhance the impact of neutron sciences on biology. Neutron crystallography is the only structural biology method capable of determining positions of all hydrogen atoms in proteins, nucleic acids and their complexes at near-physiological temperatures and of unstable species at cryogenic temperatures. Moreover, neutron analysis is non-ionizing, non-destructive and does not perturb the structure or redox chemistry of active site metal centers and clusters in proteins, which can be invaluable for studying radiation-sensitive metalloprotein complexes. Further, neutron energies used in scattering applications are similar to atomic motions, permitting neutron spectroscopies to characterize the dynamics of biomacromolecules on the picosecond to microsecond timescales. The different sensitivities of neutrons to protium (H) and deuterium (D) isotopes of hydrogen allow enhanced visibility of specific parts of biological complexes through isotopic labeling. The impact of neutrons will be most powerful when neutron scattering is combined with complementary experimental techniques that use photons and electrons, and with high-performance computing. The interconnection and mutuality of the experimental and theoretical capabilities will drive discoveries in biological and health sciences to generate more complete picture of complex biological systems. The major limitation in the field of biological neutron crystallography has been signal-to-noise, demanding large samples that are difficult to produce for the majority of biomacromolecules and limiting the applicability of this technique in biological sciences. A neutron crystallography instrument at the Second Target Station will revolutionize biological science with neutrons by engaging a large scientific community of structural biologists, enabling successful neutron diffraction experiments from radically smaller biomacromolecular crystals, resolving unanswered biochemical questions, and meaningfully contributing to rational drug design. The meeting highlighted 10 grand challenges that will be addressed with this advanced capability over the next decade and beyond, and the recommendations required to help address them are given below.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

FY24 Laboratory Directed Research and Development Annual Report

The Laboratory Directed Research and Development (LDRD) program yields foundational scientific research and development (R&D) essential to growing SRNL’s core competencies, in alignment with SRNL’s Strategic Plan to provide long-term benefits to the Department of Energy (DOE), the National Nuclear Security Administration (NNSA), and other customers and stakeholders. Five strategic goals are outlined in SRNL’s strategic plan: 1) Provide applied science and engineering for EM’s active clean-up sites and LM’s post closure management sites 2) Provide science-based solutions for gaps identified in nonproliferation strategic vision and support the government in activities impacting national security 3) Lead Science, Technology & Engineering as the central technical authority for processing tritium loaded reservoirs and support production of plutonium pits 4) Align science and energy security programs by focusing modern modeling, simulation, and data analytics tools on materials engineering and performance applications 5) Build a workforce for the future

Clark, Sue [Savannah River National Laboratory (SR↗

DOE OSTI2025

Relativistic gas accretion onto supermassive black hole binaries from inspiral through merger

Accreting supermassive black hole binaries are powerful multimessenger sources emitting both gravitational and electromagnetic (EM) radiation. Understanding the accretion dynamics of these systems and predicting their distinctive EM signals is crucial to informing and guiding upcoming efforts aimed at detecting gravitational waves produced by these binaries. To this end, accurate numerical modeling is required to describe both the spacetime and the magnetized gas around the black holes. In this paper, we present two key advances in this field of research. First, we have developed a novel 3D general relativistic magnetohydrodynamics (GRMHD) framework that combines multiple numerical codes to simulate the inspiral and merger of supermassive black hole binaries starting from realistic initial data and running all the way through merger. Throughout the evolution, we adopt a simple but functional prescription to account for gas cooling through photon emission. Next, we have applied our new computational method to follow the time evolution of a circular, equal-mass, nonspinning black hole binary for ∼200 orbits, starting from a separation of 20⁢𝑟 𝑔 and reaching the postmerger evolutionary stage of the system. We have shown how mass continues to flow toward the binary even after the binary “decouples” from its surrounding disk, but the accretion rate onto the black holes diminishes. We have identified how the minidisks orbiting each black hole are slowly drained and eventually dissolve as the binary compresses. We confirm previous findings that the system’s luminosity decreases by a factor of a few during inspiral; however, we observe an abrupt increase by ∼50% in this quantity at the time of merger, likely accompanied by an equally abrupt change in spectrum. Lastly, we have demonstrated that during the inspiral, fluid ram pressure regulates the fraction of the magnetic flux transported to the binary that attaches to the black holes’ horizons.

Accretion disk & black-hole plasma↗

DOE OSTI2025

Development of a Robust Reference Electrode in Aggressive Chemical and Radiation Environments in the Hanford Waste Tanks

The Hanford site stores more than 200 million liters of radioactive and chemically hazardous wastes from the production of weapons materials. The wastes are stored in 177 underground carbon-steel storage tanks, separated between 149 single shell tanks (SSTs) and 28 double shell tanks (DSTs). The DSTs provide critical retrieval and interim storage before the waste is vitrified in the Waste Treatment and Isolation Plant (WTP). The tanks will need to remain in-service far beyond the initial 40-year design life, and effective corrosion control practices must remain in force to extend the tanks’ lifespans. This effort includes direct measurements of corrosion rate (e.g., ultrasonic measurements and corrosion coupons) and electrochemical processes (e.g., linear polarization measurements and open circuit potential measurements). The Hanford site began monitoring the corrosion potential in select DSTs in 2008. Of the 45 reference electrodes that have been installed, 29 have failed and 6 others provided unreliable results. DOE-EM is supporting a 3-year program to develop a chemical and radiation resistant reference electrode for application in the Hanford tanks. The first year of the program focused on understanding the failure mechanism for the reference electrodes and identification of candidate construction materials that would mitigate degradation of the electrodes in the waste environment. During the second year of the program, the objectives were to: 1) test candidate materials under simulated waste conditions, 2) design components that will extend the service life of the electrode, 3) fabricate materials for prototype reference electrodes, and 4) assemble prototype reference electrodes for accelerated testing. The reference electrode is constructed of four principal parts: 1) junction, 2) casing, 3) inner chamber backfill materials, and 4) the sensing wire. Principally, improvements of the junction, casing, and inner chamber backfill materials are being pursued. The junction material at the interface between the waste and the inner chamber of the reference electrode was identified as a critical component in the failure of the reference electrodes. Nine candidate replacement junction materials were tested under simulated waste conditions to evaluate permeation rate. These materials included a variety of polymeric and ceramic materials, some of which were 3-D printed. Thus far, porous polyvinylidene fluoride materials have performed satisfactorily and are being considered for prototype development. The commercial electrode casing materials in general have performed well. Additionally, 3-D printing of chemically and mechanically stable materials is being investigated as a means for further improvement in fabrication consistency. SRNL has also investigated altering the reference electrode design to extend the service life. The new design of the interior of the reference electrode casing creates a longer, more tortuous path between the junction material and the electrode sensing wire. A finite element model was used to optimize the design without adversely impacting the circuit resistance of the electrode during the measurements, thus preserving the measurement accuracy while enhancing the service life. The inner chamber back fill materials are also critical to the performance of the reference electrode. Materials that are resistant to intruding tank waste and provide a conductive path to the sensing wire were investigated. Gel and powder materials that are interspersed with a conductive chloride bearing material were tested for their influence on diffusion and electrode resistance. All the investigated materials and components will be assembled, with collaboration from commercial vendors, to fabricate the initial prototypes. Accelerated testing of the prototypes will be initiated in Year 2 of the program and will be completed in Year 3. A recommendation on the materials of construction and the design of the new robust reference electrode will be presented to the Hanford tank farm facility.

Sykes, Kiana [Savannah River National Laboratory (↗

DOE OSTI2024

Estimating Sparse Direct Effects in Multivariate Regression With the Spike-and-Slab LASSO

The multivariate regression interpretation of the Gaussian chain graph model simultaneously parametrizes (i) the direct effects of p predictors on q outcomes and (ii) the residual partial covariances between pairs of outcomes. We introduce a new method for fitting sparse versions of these models with spike-and-slab LASSO (SSL) priors. We develop an Expectation Conditional Maximization algorithm to obtain sparse estimates of the p × q matrix of direct effects and the q × q residual precision matrix. Our algorithm iteratively solves a sequence of penalized maximum likelihood problems with self-adaptive penalties that gradually filter out negligible regression coefficients and partial covariances. Because it adaptively penalizes individual model parameters, our method is seen to outperform fixed-penalty competitors on simulated data. We establish the posterior contraction rate for our model, buttressing our method’s excellent empirical performance with strong theoretical guarantees. Using our method, we estimated the direct effects of diet and residence type on the composition of the gut microbiome of elderly adults.

EM algorithm↗

DOE OSTI2025

Discrete-Event Model of WIPP Operations

The Waste Isolation Pilot Plant (WIPP) is the critical component of the Department of Energy's (DOE) Transuranic Radioactive Waste (TRU) disposition infrastructure. Quantifying the operations ofWIPP with a discrete-event model demonstrates the capability to assess efficiency, identity bottlenecks, and improve future operations. Such a model has been developed with the simulation software ExtendSim. This report outlines the structure of that model, summarizes the model's successful reproduction of the annual amount of emplaced waste, and demonstrates that WIPP is successfully receiving and emplacing waste at a rate consistent with the rate at which the waste arrives. The model serves as a first step toward future production enhancements at WIPP. Those enhancements will rely on close collaboration with the Carlsbad Field Office (CBFO), accurate interpretation and incorporation of the model results, and effective planning with other DOE Environmental Management (DOE-EM) entities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Nonsolvating Fluoroaromatic Cosolvent Enabled Long-Term Cycling of High-Voltage Lithium-Ion Batteries with Organosulfur Electrolytes

Here, the structure–activity relationships of nonsolvating cosolvents for organosulfur-based electrolyte systems were revealed. The performance of nonsolvating dilutant fluorobenzene (FB) was compared to various fluorinated ether dilutants in high-voltage electrolytes containing a concentration of 1.2 M LiPF 6 dissolved in fluoroethylene carbonate (FEC), ethyl methyl sulfone (EMS), and the dilutant. In a high-voltage and high-loading LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) full cell configuration, the organosulfur-based electrolyte containing FB dilutant enabled superior electrochemical performance compared to the electrolytes using other nonsolvating fluorinated ether formulations. Moreover, the FB-containing electrolyte exhibited the highest ionic conductivity and lowest viscosity among all organosulfur-based electrolytes containing nonsolvating dilutant. These improvements are attributed to the enhanced physical properties of electrolyte and lithium-ion mobility. Furthermore, by employing first-principles simulations, the observed suppression of side reactions at high voltage is linked to FB’s lower reactivity toward singlet dioxygen, which is likely produced at the NMC interface. Overall, FB is considered an excellent diluent that does not impede cell operation by mass decomposition at the cathode.

25 ENERGY STORAGE↗

DOE OSTI2025

Quantitative Modeling of High-Energy Electron Scattering in Thick Samples Using Monte Carlo Techniques

Cryo-electron microscopy (cryo-EM) is a powerful tool for imaging biological samples but is typically limited by sample thickness, which is restricted to a few hundred nanometers depending on the electron energy. However, there is a growing need for imaging techniques capable of studying biological samples up to 10 µm in thickness while maintaining nanoscale resolution. This need motivates the use of mega-electron-volt scanning transmission electron microscopy (MeV-STEM), which leverages the high penetration power of MeV electrons to generate high-resolution images of thicker samples. In this study, we employ Monte Carlo simulations to model electron–sample interactions and explore the signal decay of imaging electrons through thick specimens. By incorporating material properties, interaction cross-sections for energy loss, and experimental parameters, we investigate the relationship between the incident and transmitted beam intensities. Key factors such as detector collection angle, convergence semi-angle, and the material properties of samples were analyzed. Our results demonstrate that the relationship between incident and transmitted beam intensities follows the Beer–Lambert law over thicknesses ranging from a few microns to several tens of microns, depending on material composition, electron energy, and collection angles. The linear depth of silicon dioxide reaches 3.9 µm at 3 MeV, about 6 times higher than that at 300 keV. Meanwhile, the linear depth of amorphous ice reaches 17.9 µm at 3 MeV, approximately 11.5 times higher than that at 300 keV. These findings are crucial for advancing the study of thick biological and semiconductor samples using MeV-STEM.

36 MATERIALS SCIENCE↗