SEARCH · Search NASA

Results for “dependability”

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 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Direction-dependent elastic wave scattering and mode coupling in elastic plates

Asymmetric scattering is a phenomenon in which the field scattered from a discontinuity is dependent on the direction of incidence. In waveguide systems such as elastic plates, the existence of multiple propagating modes provides a platform to explore asymmetric scattering through direction-dependent mode coupling. This paper describes this concept in the context of reciprocal systems and how to utilize it in a general manner.

Haberman, Michael↗

DOE OSTI2026

The anaerobic fungus Caecomyces churrovis produces H2 via a non-3 bifurcating NADH-dependent enzyme complex

Anaerobic fungi (AF) decompose lignocellulose-based biomass into fermentable sugars through the production of powerful biomass-degrading enzymes. AF are unusual among fungi in that they generate energy via hydrogenosomes, which are also associated with the release of H2 though yet unknown metabolic mechanisms. In particular, it remains unclear how NAD(P)+ is regenerated within hydrogenosomes and how H2 is formed. Here, we reveal the molecular mechanism for hydrogenosomal H2 production in the AF strain C. churrovis by combining genomic search, proteomic analysis, and enzymology. Our enzyme assays on the large organelle fraction of C. churrovis revealed the activity of H2:NAD+ oxidoreductase but not pyruvate:ferredoxin oxidoreductase activity. We identified genes encoding [FeFe] hydrogenase (Hyd) and NADH dehydrogenase subunits E and F (NuoE, NuoF) in C. churrovis, and confirmed their expression in the isolated hydrogenosomal fractions by proteomic analysis. Combining the individually purified proteins, we found that the assay system consisting of Hyd-Strep and NuoEF-Strep reduced NAD+ with H2. Furthermore, this system formed H2 directly from NADH independent of ferredoxin, functioning as a non-bifurcating NADH-dependent enzyme rather than an electron-bifurcating enzyme. We identified homologs of hydrogenosomal NuoE, NuoF, and Hyd in many other AF, indicating this pathway is widely conserved among the early-branching AF. This work demonstrates the existence of a non-bifurcating NADH-dependent enzyme complex in eukaryotes. Moreover, this complex could be a target for controlling AF H2 production and altering fungal metabolism.

fungi↗

DOE OSTI2025

System size and energy dependence of the mean transverse momentum fluctuations at the LHC

Event-by-event fluctuations of the event-wise mean transverse momentum, $\langle$p T $\rangle$, of charged particles produced in proton–proton (pp) collisions at $\sqrt{s}$ = 5.02 TeV, Xe–Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV, and Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV are studied using the ALICE detector based on the integral correlator $\langle$$\langle$Δp T Δp T $\rangle$$\rangle$. The correlator strength is found to decrease monotonically with increasing produced charged-particle multiplicity measured at midrapidity in all three systems. In Xe–Xe and Pb–Pb collisions, the multiplicity dependence of the correlator deviates significantly from a simple power-law scaling as well as from the predictions of the HIJING and AMPT models. The observed deviation from power-law scaling is expected from transverse radial flow in semicentral to central Xe–Xe and Pb–Pb collisions. In pp collisions, the correlation strength is also studied by classifying the events based on the transverse spherocity, S 0 , of the particle production at midrapidity, used as a proxy for the presence of a pronounced back-to-back jet topology. Low-spherocity (jetty) events feature a larger correlation strength than those with high spherocity (isotropic). The strength and multiplicity dependence of jetty and isotropic events are well reproduced by calculations with the PYTHIA 8 and EPOS LHC models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Flavour-dependent chemical freeze-out of light nuclei in relativistic heavy-ion collisions

We study the production of light nuclei in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7–200 GeV and Pb+Pb collisions at $\sqrt{s_{NN}}$ = 2.76 and 5.02 TeV within a flavour-dependent chemical freeze-out scenario, assuming different flavoured hadrons undergo separate chemical freeze-out. Using the Thermal-FIST package, thermal parameters extracted from fits to various sets of hadron yields, including and excluding light nuclei, are used to calculate the ratios of the yields of light nuclei, namely, d/p, $\overline{d}$/$\overline{p}$, t/p, t/d, 4 He/ 3 He, and 3 H Λ / 3 He. A comparison with experimental data from the STAR and ALICE collaborations shows that a sequential freeze-out scenario provides a better description of light nuclei yield ratios than the traditional single freeze-out approach. These results suggest the flavour-dependent chemical freeze-out for final state light nuclei production persists in heavy-ion collisions at both RHIC and LHC energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Quantifying the Temperature Dependence of the Multi-Species, Multi-Reaction Model. Part 1: Parameterization for a Meso-Carbon Micro-Bead Graphite

In this work, the temperature impact on the Multi-Species, Multi-Reaction (MSMR) model is studied. This is accomplished by acquiring data from slow rate lithiation and delithiation of a meso-carbon micro-bead (MCMB) graphite. The MSMR model is used to simulate linear-sweep voltammetry data of a porous electrode composed of graphite, and because the electrode is close to a state of dynamic equilibrium, the peaks in the differential voltage spectroscopy plot can be analyzed. Through this analysis, the temperature impact on the total fraction of available host sites in a particular MSMR gallery (X j ), the impact on the reference potential (U° j ), and the impact on the parameter detailing the deviation from Nernstian behavior ( ω j ) can be found. This is the first time the temperature dependence of the MSMR parameters have been experimentally analyzed. In Part 2, the impact of the temperature dependence of the MSMR parameters on the entropy coefficient of an intercalation material will be studied.

Paul, Abigail (ORCID:0000000172892069)↗

DOE OSTI2024

Dependency Analysis Project

The Dependency Analysis Project creates datasets from package management ecosystems to determine what package a file was likely installed by, and provides tools to recognize implicit dependencies when source code is available for a variety of language ecosystems. The overarching goal is to help answer the question: what software does a system use?

Mast, Ryan↗

DOE OSTI2025

Dependence of Convective Cloud Microphysical Properties on Environmental Conditions during the TRACER and ESCAPE Field Campaigns: A Synergistic Approach of Observations, Machine Learning and Parcel Models

The sensitivity of convective clouds to aerosols and their interactions with environment, combined with limited observational constraints in parameterizations, introduces significant uncertainties in atmospheric models. Here, this study investigates the dependence of convective cloud microphysical properties on environmental conditions using a synergistic approach that combines unique observations from the TRACER and ESCAPE field campaigns, machine learning techniques, and parcel model simulations with a super-droplet microphysics scheme. A random forest algorithm identifies in-situ vertical velocity (w), temperature (T), and surface fine-mode aerosol mass concentration as the three most important environmental conditions influencing cloud properties including liquid water content (LWC), number concentration for particles with D max < 50 μm (N c ,<50), 50 μm ≤ D max ≤ 3000 μm (N c,50–3000 ), and droplet effective diameter (D e ). Results show that LWC, N c,<50 , and N c,50–3000 significantly increase with w in updrafts. Across w bins, as T decreases, LWC, D e , and N c,50–3000 increase, while N c,<50 decreases, which are closely linked to the distance above cloud bases. Warmer cloud bases yield higher LWC, greater N c,50–3000 , and smaller N c,<50 , while polluted environments produce greater N c,<50 . Parcel model simulations successfully replicate these observed dependencies. The simulation results indicate that warmer cloud bases enhance condensation generating larger droplets, and differences in droplet sizes are then amplified through collision-coalescence, resulting in a greater N c,50–3000 . Polluted conditions result in a greater N c,<50 primarily due to enhanced cloud condensation nuclei activation despite increased collision-coalescence rates compared to pristine conditions. This study provides observed quantitative patterns characterizing cloud microphysical properties as a function of key environmental parameters, offering valuable constraints for improving physics parameterizations and numerical models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Numerical eigen-spectrum slicing, accurate orthogonal eigen-basis, and mixed-precision eigenvalue refinement using OpenMP data-dependent tasks and accelerator offload

Performing a variety of numerical computations efficiently and, at the same time, in a portable fashion requires both an overarching design followed by a number of implementation strategies. All of these are exemplified below as we present transitioning the PLASMA numerical library from relying on dependence-driven large tasks to achieving utilization of fine grain tasking and offload to hardware accelerators while keeping its core dependence sets: OpenMP source code pragmas and runtime for most system-level functionality and basic low-level numerical kernels provided directly by hardware vendors or open source projects with vendor contributions. We also present new algorithmic methods and their efficient parallel implementations including fine grained tasking for eigen-spectrum slicing and offload for mixed-precision eigenvalue refinement. We provide performance, scaling, and numerical results showing sizable gains over the available solutions from either the open source and vendor-provided packages.

Luszczek, Piotr↗

DOE OSTI2026

Dynamic Transmission Line Switching Amid Wildfire-Prone Weather Under Decision-Dependent Uncertainty

During dry and windy seasons, environmental conditions significantly increase the risk of wildfires, exposing power grids to disruptions caused by transmission line failures. Wildfire propagation exacerbates grid vulnerability, potentially leading to prolonged power outages. To address this challenge, we propose a multistage optimization model that dynamically adjusts transmission grid topology in response to wildfire propagation, aiming to develop an optimal response policy. By accounting for decision-dependent uncertainty, where line survival probabilities depend on usage, we employ distributionally robust optimization to model uncertainty in line survival distributions. We adapt the stochastic nested decomposition algorithm and derive a deterministic upper bound for its finite convergence. To enhance computational efficiency, we exploit the Lagrangian dual problem structure for a faster generation of Lagrangian cuts. Using realistic data from the California transmission grid, we demonstrate the superior performance of dynamic response policies against two-stage alternatives through a comprehensive case study. In addition, after solving the multistage formulation, we construct easy-to-implement policies that significantly reduce computational burden while maintaining good performance in real-time deployment. History: Accepted by Russell Bent, Area Editor for Network Optimization: Algorithms and Applications. Funding: This work was supported by the U.S. Department of Energy, Office of Electricity [Grant DE-AC02-05CH11231]. The work of R. Jiang was supported in part by the U.S. National Science Foundation, Division of Electrical, Communications and Cyber Systems [Grant ECCS-1845980] and the U.S. Air Force Office of Scientific Research [Grant FA9550-23-1-0323]. Supplemental Material: The software that supports the findings of this study is available within the paper and its Supplemental Information ( https://pubsonline.informs.org/doi/suppl/10.1287/ijoc.2025.1210 ) as well as from the IJOC GitHub software repository ( https://github.com/INFORMSJoC/2025.1210 ). The complete IJOC Software and Data Repository is available at https://informsjoc.github.io/ .

Estrada-Garcia, Juan-Alberto↗

DOE OSTI2024

Modeling dynamics of acute HIV infection incorporating density-dependent cell death and multiplicity of infection

Understanding the dynamics of acute HIV infection can offer valuable insights into the early stages of viral behavior, potentially helping uncover various aspects of HIV pathogenesis. The standard viral dynamics model explains HIV viral dynamics during acute infection reasonably well. However, the model makes simplifying assumptions, neglecting some aspects of HIV infection. For instance, in the standard model, target cells are infected by a single HIV virion. Yet, cellular multiplicity of infection (MOI) may have considerable effects in pathogenesis and viral evolution. Further, when using the standard model, we take constant infected cell death rates, simplifying the dynamic immune responses. Here, we use four models—1) the standard viral dynamics model, 2) an alternate model incorporating cellular MOI, 3) a model assuming density-dependent death rate of infected cells and 4) a model combining (2) and (3)—to investigate acute infection dynamics in 43 people living with HIV very early after HIV exposure. We find that all models qualitatively describe the data, but none of the tested models is by itself the best to capture different kinds of heterogeneity. Instead, different models describe differing features of the dynamics more accurately. For example, while the standard viral dynamics model may be the most parsimonious across study participants by the corrected Akaike Information Criterion (AICc), we find that viral peaks are better explained by a model allowing for cellular MOI, using a linear regression analysis as analyzed by R 2 . These results suggest that heterogeneity in within-host viral dynamics cannot be captured by a single model. Depending on the specific aspect of interest, a corresponding model should be employed.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2026

Microwave-Assisted Hydrogen Generation from Hydrocarbon-Bearing Reservoir Rocks: Stage-Dependent Thermal Runaway and In-Situ Carbonate Engineering

Microwave-assisted hydrogen generation from hydrocarbon-bearing reservoir rocks is strongly influenced by mineralogy, methane activation, carbonate reactions, and thermal runaway behavior. This study investigates a new approach in which carbonate phases are generated in-situ through the reaction of internally produced CO2 with Ca(OH)2 under microwave heating conditions. The objective is to evaluate how rock mineralogy, methane injection, and Ca(OH)2 addition influence hydrogen generation, carbon redistribution, and stage-dependent reaction pathways during microwave exposure. Microwave heating experiments were conducted using Permian Basin reservoir rocks under three experimental conditions: rock-only experiments under Ar atmosphere, CH4–Ar experiments without additive, and CH4–Ar experiments containing 5 wt% Ca(OH)2. Methane-assisted experiments were performed under continuous injection of 30 standard cubic centimeters per minute (sccm) CH4 and 30 sccm Ar. Based on thermal runaway behavior, each experiment was divided into three operational stages: Before Thermal Runaway (BR), After Thermal Runaway–Decrease in Microwave Power (ARD), and After Thermal Runaway– Increase in Microwave Power (ARI). Temperature and gas composition were continuously monitored throughout the experiments. The rock-only experiments demonstrated that hydrogen generation can occur intrinsically from hydrocarbon-bearing rocks under microwave heating, even without externally injected methane. However, hydrogen production did not correlate solely with kerogen content, indicating that mineralogy strongly influences hydrogen-generation pathways. Correlation analyses suggested that kerogen decomposition initially generated CH4, CO, and CO2, followed by secondary hydrocarbon reactions associated with H2 and C2 hydrocarbon formation. Methane-assisted experiments substantially increased hydrogen production; however, identical methane injection rates produced significantly different hydrogen yields among the rock samples, confirming that mineralogical composition controls methaneconversion behavior under microwave heating conditions. The addition of Ca(OH)2 significantly altered carbon evolution behavior in a stage-dependent manner. During the BR stage, Ca(OH)2 reduced gas-phase CO2 production, particularly in carbonate-rich rocks, indicating favorable conditions for in-situ carbonation and carbonate deposition prior to extensive thermal decomposition. The suppression of CO2 during BR became more pronounced with increasing carbonate content of the rock system. After thermal runaway, carbonate-containing systems exhibited enhanced hydrogen generation behavior, suggesting that carbonatederived mineral transformations and carbonate-mediated reactions contribute to high-temperature hydrogen-generation pathways. In carbonate-poor rocks, Ca(OH)2 addition enabled simultaneous URTeC 4493775 2 enhancement of hydrogen production and partial suppression of CO2 release during the post-runaway stages. Overall, the results demonstrate that microwave-assisted hydrogen generation is governed by dynamically evolving interactions among kerogen decomposition, methane activation, mineral transformations, carbonate formation/decomposition, and thermal runaway behavior. This work introduces in-situ carbonate engineering with Ca(OH)2 as a strategy for coupling hydrogen generation with partial insitu carbon management under microwave heating conditions.

03 NATURAL GAS↗

DOE OSTI2025

Phase-field modeling of orientation-dependent crack growth in ductile single crystals with anisotropic elasticity

Crack growth in ductile single crystals (DuSCs) is orientation dependent due to the anisotropies of crystal plasticity and elastic tensor. This study develops a phase-field model incorporating both crystal plasticity and crack growth and proposes a general method to decompose the elastic energy into compressive and tensile parts to prevent crack growth under compression in the phase-field description. The phase-field model, in combination with three Euler angles, is employed to simulate orientation-dependent crack growth in DuSCs. The contributions from crystal plasticity and anisotropic elasticity are compared, and the former is found to dominate in the anisotropy of crack growth in copper single crystals. Furthermore, the simulation results demonstrate that crystal orientation strongly affects the heterogeneous distribution of plastic strain and the interaction between plastic strain and crack growth. High-throughput phase-field simulations are performed with exhaustive crystal orientations, and the results are explained based on the anisotropy of the Taylor factor.

Computational Solid Mechanics↗

DOE OSTI2024

Thickness-dependent transport in epitaxial CaIrO 3 perovskite thin films

The coexistence of strong electron–electron correlation and strong spin–orbit coupling makes CaIrO 3 a promising candidate for the exploration of emergent phenomena in thin film form. Here, we have realized epitaxial CaIrO 3 films on (La 0.18 Sr 0.82 )(Al 0.59 Ta 0.41 )O 3 single crystal substrates with excellent crystallinity and thicknesses ranging from 3 to 32 nm. Electronic properties range from semiconducting to insulating behavior as a function of decreasing thickness and correlate with increased tetragonal distortion due to substrate-imposed strain. Hall effect measurements indicate that conduction is dominated by electron-type carriers at room temperature but crosses over to transport dominated by hole-like carriers at lower temperatures. Longitudinal magnetoresistance measurements reveal both a quadratic Lorentz dependence and an additional positive component in the field dependence that cannot be explained purely by semi-classical orbital magnetoresistance effects.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Reducing the Parameter Dependency of Phase-Picking Neural Networks with Dice Loss

Training a neural network for picking seismic phase arrivals has been commonly posed as a segmentation problem. It is a highly imbalanced segmentation problem in the sense that the background vastly dominates the foreground because we are trying to pick the optimal single sample point that represents the arrival of a seismic phase in a many seconds long time window. Here, we test the Dice loss, which is a preferred loss function for highly imbalanced image segmentation problems. We show that phase-picking neural networks trained on the Dice loss behave in a binary fashion for which the prediction output is almost always either nearly 1 or nearly 0. This feature removes the strong dependence of data processing workflows on the prediction score threshold, which is an otherwise critical parameter to determine when using neural networks trained on the cross-entropy loss. When strategically used, models trained on the Dice loss can reduce the parameter dependency of machine learning-based seismic monitoring.

58 GEOSCIENCES↗

DOE OSTI2024

Hydropower Black Start: A Guidebook for Retrofitting Grid Dependent Hydropower

Not all United States (US) hydropower plants were designed to provide black start, but they are increasingly needed to uphold resilience in the evolving electric grid. This guidance is designed to help understand the minimal retrofits required for grid dependent hydropower (GDH) plants behind the point of interconnection (POI). For distribution connected hydropower plants or those with dedicated cranking paths, such upgrades can be sufficient for the plant to provide black start. For others, more coordination with the transmission system operator will be needed. This guidebook answers a number of questions relevant to retrofitting hydropower plants with black start capabilities. For example, the guidebook answers: • How flexible do the wicket gate controls need to be? • Who needs to do hydro governor model validation, why, and how? • How robust and flexible do the excitation and AVR controls need to be? • What protection settings need to be adjusted? • What relay(s) will need to be bypassed or overridden and at what risk? • What is the electrical energy demand of the station load or auxiliary power systems? • What should the strategy to energize transformer(s) along cranking path to address inrush currents be? • How should the critical load restoration be sequenced? In addition to outlining the specifications that hydropower plants need to meet for each component to be able to perform black start, this guidebook provides a set of case studies for specific upgrades needed at actual plants. Between the case studies of plants that have already performed black start retrofits and the examples of how this guidebook can be applied to scope future retrofits, five key themes have been identified for retrofit needs. 1. Protection needs “black start” mode: hydropower plants that are not designed with black start capabilities will have protections that prevent them from interconnecting to a “dead bus.” These protections will need to be overridden in every retrofit case and a separate black start mode should be established so that operators can safely switch between black start and grid connected modes, minimizing the risk to the plant. 2. Wicket gates need modern controls: digital governors accelerate the parameter tuning process and gate position sensors improve controllability, so plants with mechanical governors should be upgraded. Furthermore, a black start and islanding mode should be established for controls to maximize plant performance. 3. Robust excitation support: the DC system or excitation generator needs to be reliable enough to form and sustain the rotor electromagnetic field. These systems are typically undersized in plants that were not designed for black start, so they will need to be upgraded. 4. Turbine-governor model validation and operator training: validation of a standard hydro governor model is needed to characterize the dynamic response (i.e., inertial and primary frequency response) of the GDH. This is required for control development and old hydropower plants often have outdated or incorrect models. Operator training is also typically required to ensure the hardware retrofits are utilized correctly during the black start process. 5. Transformer and cranking path energization: any upgradation and control adjustment in front of the POI will depend upon the existing interconnection. Coordination with the transmission or distribution operator may be required.

13 HYDRO ENERGY↗

DOE OSTI2025

Impact of Time Dependent Reactor and Sensor Physics on Core Power Synthesis

Online synthesis of the power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins, and to provide essential knowledge associated with the burnup of the fuel. In light water reactors (LWRs), power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report aims to study the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a LWR. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative pressurized water reactor low-enriched uranium (LEU) core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used, as well, to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative (PBI) method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly under- or over-approximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumption leads to severe levels of error, approaching 25% at maximum. The authors also discovered that V SPNDs are extraordinarily robust in the low-enriched uranium fuel cycle considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Impact of Time-Dependent Reactor and Sensor Physics on Core Power Synthesis (Rev.1)

Online synthesis of power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins and to provide essential knowledge associated with the burnup of the fuel. In light-water reactors, power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report examines the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a pressurized water reactor (PWR), considering the typical low-enriched uranium (LEU, 3%-5% enrichment) fuel cycle as well as the higher enrichment LEU+ (5%-8% enrichment) fuel cycle. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative PWR LEU core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used as well to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori power distribution assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters and the extended nature of the LEU+ fuel cycle in comparison with LEU. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly underapproximated or overapproximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumptions leads to severe levels of error, approaching 25% at maximum (for LEU). The authors also discovered that V SPNDs are extraordinarily robust in both the LEU and LEU+ fuel cycles considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Tensor Network Space-Time Spectral Collocation Method for Time-Dependent Convection-Diffusion-Reaction Equations

Emerging tensor network techniques for solutions of partial differential equations (PDEs), known for their ability to break the curse of dimensionality, deliver new mathematical methods for ultra-fast numerical solutions of high-dimensional problems. Here, we introduce a Tensor Train (TT) Chebyshev spectral collocation method, in both space and time, for the solution of the time-dependent convection-diffusion-reaction (CDR) equation with inhomogeneous boundary conditions, in Cartesian geometry. Previous methods for numerical solution of time-dependent PDEs often used finite difference for time, and a spectral scheme for the spatial dimensions, which led to a slow linear convergence. Spectral collocation space-time methods show exponential convergence; however, for realistic problems they need to solve large four-dimensional systems. We overcome this difficulty by using a TT approach, as its complexity only grows linearly with the number of dimensions. We show that our TT space-time Chebyshev spectral collocation method converges exponentially, when the solution of the CDR is smooth, and demonstrate that it leads to a very high compression of linear operators from terabytes to kilobytes in TT-format, and a speedup of tens of thousands of times when compared to a full-grid space-time spectral method. These advantages allow us to obtain the solutions at much higher resolutions.

97 MATHEMATICS AND COMPUTING↗