Search NASASearch

SEARCH · Search NASA

Results for “Temporal modeling”

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.

531 records · Page 11

Nuclear Physics Made Very, Very Easy

The fundamental approach to nuclear physics was prepared to introduce basic reactor principles to various groups of non-nuclear technical personnel associated with NERVA Test Operations. NERVA Test Operations functions as the field test group for the Nuclear Rocket Engine Program. Nuclear Engine for Rocket Vehicle Application (NERVA) program is the combined efforts of Aerojet-General Corporation as prime contractor, and Westinghouse Astronuclear Laboratory as the major subcontractor, for the assembly and testing of nuclear rocket engines. Development of the NERVA Program is under the direction of the Space Nuclear Propulsion Office, a joint agency of the U. S. Atomic Energy Commission and the National Aeronautics and Space Administration. This report is being reprinted for use in the U. S. Atomic Energy Commission and National Aeronautics and Space Administration educational and technology utilization programs.

22 GENERAL STUDIES OF NUCLEAR REACTORS

The Role of Nuclear Data Sensitivities in Prompt α-Eigenvalue Predictions of Delayed Critical Benchmarks

Alpha (α) eigenvalues, which describe the logarithmic time derivative of the neutron population in a multiplying system, are integral to time-dependent behavior and diagnostic applications. However, uncertainties in the evaluated nuclear data can significantly impact the accuracy of transport simulations for such quantities. This work explores the use of machine learning models to predict two key outputs, α-eigenvalues and keff bias, using input features derived from α-eigenvalue sensitivities to nuclear data. The criticality safety benchmark models used in this study come from the International Handbook of Evaluated Criticality Safety Benchmark Experiments. Three models, random forest, XGBoost, and NGBoost, are trained on both energy-resolved and energy-summed α sensitivities. For the α-eigenvalue bias prediction, NGBoost achieved the highest R 2 (0.9476) using energy-resolved features, while XGBoost performed best using summed sensitivities. In contrast, when predicting the keff bias, all the models showed moderate predictive capability (best R 2 ≈ 0.72), as the mapping from the static α-sensitivities to the static keff bias was less direct. SHAP (SHapley Additive exPlanations) analysis was used to interpret the model predictions. Across both prediction tasks, the features associated with neutron capture [H-1 (n, γ)], uranium scattering reactions (such as 235 U elastic/inelastic), and actinide capture/fission reactions (such as 239 Pu and 234 U) were consistently identified as the most impactful. This highlights the key role of specific nuclear reactions and energy ranges in shaping both time-dependent and steady-state criticality behavior. These results demonstrated that α-sensitivities, despite being computed for time-dependent metrics, can provide valuable insights for predicting both α-eigenvalues and the keff bias. Moreover, machine learning models offer a promising pathway for uncovering important nuclear data dependencies and guiding future data evaluation efforts.

Nuclear data

Ion Transport in Concentrated Crosslinked Solid Polymer Electrolytes

Crosslinking polymers is a common approach to create mechanically stable solid materials such as polymer electrolytes for lithium batteries. In conventional liquid electrolytes, the solvent molecules move freely to accommodate the field-induced motion of ions. However, in crosslinked polymer electrolytes, the rearrangement of polymer segments is constrained by the deformation limits of the network. Herein, we develop a new transport model that accounts for both the formation of concentration gradients and the elasticity of the electrolyte. The elasticity is incorporated by adding an additional term related to the entropy of crosslinked strands to the electrochemical potential of the salt. The resulting Crosslink Model contains two adjustable parameters: $\mathcal{N}$, the average number of monomers in a strand, and λ crit , the maximum strain the network can sustain. These solid-like constraints produce singularities in the governing transport equations, fundamentally altering the concentration profiles. Plateaus in salt concentrations emerge near the electrodes, and network elasticity introduces a strain overpotential. When compared to a Baseline Model ($\mathcal{N}$ → ∞, equivalent to concentrated solution theory), which predicts steepest gradients near the electrodes, both models yield similar current–voltage relationships. Model predictions are compared to electrochemical data for a poly(ethylene oxide)-based crosslinked polymer electrolyte.

Patel, Vivaan [University of California, Berkeley,

Entanglement Structure of Non-Gaussian States and How to Measure It

Rapidly growing capabilities of quantum simulators to probe quantum many-body phenomena require new methods to characterize increasingly complex states. Here, we present a protocol that constrains quantum states using experimentally measured correlation functions. This method enables measurement of a quantum state’s entanglement structure, opening a new route to study entanglement-related phenomena. Our approach extends Gaussian state parameterizations by systematically incorporating higher-order correlations. We show the protocol’s usefulness in conjunction with current and forthcoming experimental capabilities, focusing on weakly interacting fermions as a proof of concept. Here, the lowest nontrivial expansion quantitatively predicts early time thermalization dynamics, including signaling the onset of quantum chaos indicated by the entanglement Hamiltonian.

Fermi gases

Superdiffusion resilience in Heisenberg chains with two-dimensional interactions on a quantum processor

Superdiffusive spin transport in the one-dimensional (1D) Heisenberg model is a key theoretical discovery in nonequilibrium quantum many-body physics. Although extensively studied in 1D systems, the breakdown and sustenance of superdiffusion in two-dimensional (2D) lattices with integrability-breaking terms, as found in real materials, remains an open question. To address this, we develop a toy model that extends the 1D Heisenberg model with a representative set of 2D interaction types and tunable strengths. Our model exhibits varying degrees of superdiffusion breakdown depending on the interaction type, spanning ballistic to diffusive regimes. We establish and justify a hierarchy of 2D interactions based on their resilience against superdiffusion breakdown: Heisenberg >𝑋⁢𝑋 > Ising. This precise control over the superdiffusive behavior also enables rigorous benchmarking of quantum hardware, and our simulations on IBM's Heron devices confirm the hardware's ability to accurately capture these many-body nonequilibrium phenomena. Overall, our results are relevant not only to simulating superdiffusion in real materials, such as the 1D Heisenberg compound KCuF3, which contains modest nonintegrable 2D terms, but also to extending superdiffusive behavior to larger 2D qubit lattices and other 2D materials.

Alagarsamy Manikandan, Keerthi Kumaran [ORNL]

Stoichiometry dependent properties of cerium hydride: An active learning developed interatomic potential study

Cerium hydride has a variety of interesting properties, including a known lattice contraction and densification with increasing hydrogen content. However, precise stoichiometric control is not experimentally straightforward and ab initio approaches are not computationally feasible for many properties such as melting and low temperature diffusion. Therefore, we develop a machine-learned interatomic potential for cerium hydride that is valid for H to Ce ratios from 2.0 to 3.0. A query-by-committee active learning approach is used to develop the training set. Leveraging classical molecular dynamics simulations, we assess a range of properties and provide fundamental mechanisms for the trends with stoichiometry. Finally, a majority of the properties follow the trend of lattice contraction, being governed by the stronger lattice binding induced by adding octahedral atoms.

36 MATERIALS SCIENCE

Dark matter and baryon asymmetry from monopole-axion interactions

We introduce a novel mechanism where the kinetic energy of a rotating axion can be dissipated by the interactions with dark magnetic monopoles. This mechanism leads to a framework where the QCD axion and dark monopoles account for the dark matter density, and the observed baryon asymmetry is generated through the rotating QCD axion via axiogenesis. The monopoles acquire masses from a nonzero axion field, and they can transition between different quantized dyonic levels in the presence of a rotating axion field. The axion kinetic energy is dissipated by the transition, and thus the axion abundance is depleted to the observed dark matter abundance. We predict that the axion decay constant should be below 10 9 GeV to explain the observed dark matter and baryon densities.

Axions and ALPs

From zonal to nodal capacity expansion planning: Spatial aggregation impacts on a realistic test-case

Solving power system capacity expansion planning (CEP) problems at realistic spatial resolutions is computationally challenging. Thus, a common practice is to solve CEP over zonal models with low spatial resolution rather than over full-scale nodal power networks. Due to improvements in solving large-scale stochastic mixed integer programs, these computational limitations are becoming less relevant, and the assumption that zonal models are realistic and useful approximations of nodal CEP is worth revisiting. Here, this work is the first to conduct a systematic computational study on the assumption that spatial aggregation can reasonably be used for ISO-scale CEP. By considering a realistic, large-scale test network based on the state of California with over 8000 buses, we find that well-designed small spatial aggregations can yield good approximations but that coarser zonal models may result in large distortions of investment decisions, e.g., capacity under-investment of up to 41% for the lowest resolution model considered.

24 POWER TRANSMISSION AND DISTRIBUTION

Lepto-axiogenesis with light right-handed neutrinos

We study lepto-axiogenesis in theories where the right-handed neutrino is light enough that its dynamics affect the determination of the baryon asymmetry. When compared with theories of high-scale lepto-axiogenesis where the Majorana neutrino mass may be treated as an effective dimension-five operator, we find that the predicted saxion mass is lower. Two distinct scenarios emerge. In the first, processes that generate the baryon asymmetry are in equilibrium down to the mass of the right-handed neutrino. In the second, the relevant processes never reach equilibrium; the baryon number freezes in. We comment on implications for supersymmetric spectra and discuss constraints on late decays of supersymmetric relics and from dark radiation. In contrast to high-scale lepto-axiogenesis, which predicts superpartners with masses of 10–100 TeV or more, we find this scenario is consistent with a wider range of superpartner masses, all the way down to current direct search bounds.

Axions and ALPs

PowerModelsGAT-AI: Physics-Informed Graph Attention for Multi-System Power Flow With Continual Learning

Solving the alternating current power flow equations in real time is essential for secure grid operation, yet classical Newton–Raphson solvers can be slow under stressed conditions. Existing graph neural networks for power flow are typically trained on a single system and often degrade on different systems. We present PowerModelsGAT-AI, a physics-informed graph attention network that predicts bus voltages and generator injections. The model uses bus-type-aware masking to handle different bus types and balances multiple loss terms, including a power-mismatch penalty, using learned weights. We evaluate the model on 14 benchmark systems (4 to 6,470 buses) and train a unified model on 13 of these under contingency conditions with up to two branch outages, achieving an average normalized mean absolute error of 0.89% for voltage magnitudes and R 2 >0.99 for voltage angles. We also show continual learning: when adapting a base model to a new 1,354-bus system, standard fine-tuning causes severe forgetting with error increases exceeding 1000% on base systems, while our experience replay and elastic weight consolidation strategy keeps error increases below 2% and in some cases improves base-system performance. Interpretability analysis shows that learned attention weights correlate with physical branch parameters (susceptance: r=0.38 ; thermal limits: r=0.22 ), and feature importance analysis supports that the model captures established power flow relationships.

24 POWER TRANSMISSION AND DISTRIBUTION

Characterizing the γ-Ray Emission from FR0 Radio Galaxies

Abstract FR0 galaxies constitute the most abundant jet population in the local Universe. With their compact jet structure, they are broadband photon emitters and have been proposed as multimessenger sources. Recently, these sources have been detected for the first time inγrays. Using a revised FR0 catalog, we confirm that the FR0 population as a whole areγ-ray emitters, and we also identify two significant sources. For the first time, we find a correlation between the 5 GHz core radio luminosity andγ-ray luminosity in the 1–800 GeV band, having a 4.8σstatistical significance. This is clear evidence that the jet emission mechanism is similar in nature for FR0s and the well-studied canonical FR (FRI and FRII) radio galaxies. Furthermore, we perform broadband spectral energy distribution modeling for the significantly detected sources as well as the subthreshold source population using a one-zone synchrotron self-Compton model. Within the maximum jet power budget, our modeling shows that the detectedγrays from the jet can be explained as inverse Compton photons. To explain the multiwavelength observations for these galaxies, the modeling results stipulate a low bulk Lorentz factor and a jet composition far from equipartition, with the particle energy density dominating over the magnetic field energy density.

Astronomy & Astrophysics

Theoretical design and performance of three-dimensional, pillared FeS 2 cathodes

Three-dimensional (3D) electrode design can provide improved capacities and rate capabilities over conventional two-dimensional electrodes by enhancing electrical and ionic transport. Here, expanding upon our previous modeling efforts for conversion chemistry lithium-ion batteries, we develop a pseudo-four-dimensional (P4D) approach that is subsequently used to investigate the design of a pillared FeS 2 electrode. The model considers transport in three dimensions with an additional “fourth” dimension corresponding to the solid-state lithium transport within the active material particles. Additionally, we allow for expansion of the active material during the conversion reaction to understand how internal stresses impact the electrochemical performance of the cell. By optimizing the model with respect to areal capacity, we are able to predict areal capacities up to 16.8 mAh/cm 2 for an areal current density of 1.78 mA/cm 2 and a 103% improvement for the three-dimensional electrodes over planar electrodes of equal volume. Despite the promising results, our simulations suggest that 3D design may be difficult for conversion cathode materials due to the large internal stresses that arise during conversion. Nevertheless, the model is robust and adaptable to other materials that may be more suitable for 3D electrodes due to a lesser change in volume during discharge.

Conversion cathode materials

HTGR Multiphysics Application Drivers FY26 Updates

This report summarizes FY26 progress under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program's high-temperature gas-cooled reactor (HTGR) application driver work, covering a wide range of activities such as code validation and multi-physics code assessment. 1) A detailed SAM model of the High-Temperature Engineering Test Reactor (HTTR) was developed using a unique-block grouping approach, with an extended parallel thermal network method to capture block-to-block conduction and radiation heat transfer, and applied to steady-state simulations of the HTTR 30~MW and 9~MW cases. 2) In another activity, SAM's newly implemented multi-component gas flow model was validated against the Natural convection Shutdown heat removal Test Facility (NSTF) argon ingress experiment, correctly capturing the density-driven suppression and thermal recovery of natural circulation observed when argon is introduced into the air-cooled Reactor Cavity Cooling System (RCCS) loop. 3) For the OECD/NEA High Temperature Test Facility (HTTF) benchmark, we co-led the international benchmark activities as well as the OECD/NEA final benchmark report to be released at the end of this year. 4) Finally, the coupled Griffin-SAM modeling capability for pebble-bed HTGRs was advanced by verifying the Griffin neutronics solution against Serpent Monte Carlo for a realistic non-uniform temperature distribution, resolving several deficiencies in the SAM-to-Griffin temperature transfer scheme, and enabling distinct fuel kernel, moderator, and coolant temperatures for cross section feedback. These new features were demonstrated in a PBR load-following transient.

Lee, Alvin

Nonlocal Effect of Percolated Particle Networks on Viscoelasticity of Polymer–Filler Nanocomposites: A Mesoscale Simulation Study

With nanoparticles (NPs) as fillers, polymer nanocomposites (PNCs) usually exhibit enhanced mechanical properties. However, a direct connection between the microscopic structural relaxation and macroscopic mechanical properties of PNCs remains to be established. To investigate the micro-to-macro connection, we develop a mesoscale model, in which the NPbridging polymer chains are represented by a dynamic bonded interaction between NPs, and the bulk polymer matrix is implicitly modeled by overdamped Langevin dynamics. Extensive equilibrium simulations are performed to quantify the microscopic dynamics of model PNCs. Systematic analyses of modified Rouse dynamics, dynamic structure factor, and relaxation modulus uncover that the microscopic relaxation dynamics of PNCs are significantly decelerated across different length scales because of nonlocal effects of percolated particle networks a phenomenon that has not been adequately captured in prior simulation studies. We find that NPvolume- fraction and NP-bonding-energy barrier are the two critical variables that affect bulk viscoelasticity the most. The proposed mesoscale model is versatile and provides a powerful framework for studying structure−property relations of different PNCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Divergent evolution of slip banding in CrCoNi alloys

Abstract Metallic materials under high stress often exhibit deformation localization, manifesting as slip banding. Over seven decades ago, Frank and Read introduced the well-known model of dislocation multiplication at a source, explaining slip band formation. Here, we reveal two distinct types of slip bands (confined and extended) in compressed CrCoNi alloys through multi-scale testing and modeling from microscopic to atomic scales. The confined slip band, characterized by a thin glide zone, arises from the conventional process of repetitive full dislocation emissions at Frank–Read source. Contrary to the classical model, the extended band stems from slip-induced deactivation of dislocation sources, followed by consequent generation of new sources on adjacent planes, leading to rapid band thickening. Our findings provide insights into atomic-scale collective dislocation motion and microscopic deformation instability in advanced structural materials.

Science & Technology - Other Topics

Working Fluid Characterization and Performance Assessment of Subcritical Organic Rankine Cycles Based on the Lee–Kesler Approach for Energy Recovery

Here, a generalized model using the Lee–Kesler approach based on the corresponding states principle is developed to assess the performance of subcritical Organic Rankine Cycles operating with different working fluids. Each fluid is characterized by five parameters: the acentric factor, critical temperature, critical pressure, molar mass, and the ideal-gas ratio of specific heats at the critical temperature. The model was developed using the compressibility factor modified version of the Benedict–Webb–Rubin equation proposed by Lee and Kesler and the enthalpy and entropy functions to calculate thermodynamic state properties. The model was validated by comparing the results calculated with the model and working fluid thermodynamic properties obtained with the CoolProp database. This comparison was conducted for 91 working fluids, obtaining a relative error below 5% for 88 out of the 91 fluids (∼97%). A generalized parametric study was conducted to determine the influence of the pinch point and each fluid parameter on the performance of Organic Rankine Cycle (ORC) systems. It was found that efficiency increases with critical temperature, ideal-gas ratio of specific heats at the critical temperature, and acentric factor, reaching up to 13%. The developed model enables the evaluation of ORC system performance for existing working fluids. It also allows the formulation and evaluation of new fluids to enhance the performance of the ORC while retrieving energy from any kind of source; and likewise, the methodology can be applied to other power generation cycles.

24 POWER TRANSMISSION AND DISTRIBUTION

Heavy long-lived dark vector via a gluonic portal

We study a dark gauge boson Z′ that exclusively couples to the QCD gluons through higher dimensional operators. These operators are generated from integrating out of heavy ultraviolet resonances carrying both QCD and dark gauge charges. With SU(3)C gauge invariance, charge and parity symmetries preserved, we find that the leading effective operators are restricted to have the form of Z′GGG and Z′Z′GG at dimension-eight, which can naturally render the Z′ particle long-lived, and serve as a viable dark matter candidate. We investigate the phenomenology of these operators with both collider experiments and cosmological observation, without and with the assumption that this dark gauge boson plays the role of the dominant dark matter component. For an unstable Z′, we show that depending on its lifetime, it can be probed by various observables up to ultraviolet physics scale around 10 9 GeV. For Z′ being dark matter, we find that $m_{Z'}$ ≳ 1 TeV is consistent with the thermal freeze-out scenario. In contrast, in the freeze-in scenario, the extremely small couplings leave the relevant parameter space largely unconstrained by current experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Internally Catalyzed Hydrogen Atom Transfer (I-CHAT)—A New Class of Reactions in Combustion Chemistry

The current paradigm of low-T combustion and autoignition of hydrocarbons is based on the sequential two-step oxygenation of fuel radicals. The key chain-branching occurs when the second oxygenation adduct (OOQOOH) is isomerized releasing an OH radical and a key ketohydroperoxide (KHP) intermediate. The subsequent homolytic dissociation of relatively weak O–O bonds in KHP generates two more radicals in the oxidation chain leading to ignition. Based on the recently introduced intramolecular “catalytic hydrogen atom transfer” mechanism (J. Phys. Chem. 2024, 128, 2169), abbreviated here as I-CHAT, we have identified a novel unimolecular decomposition channel for KHPs to form their classical isomers—enol hydroperoxides (EHP). The uncertainty in the contribution of enols is typically due to the high computed barriers for conventional (“direct”) keto–enol tautomerization. Remarkably, the I-CHAT dramatically reduces such barriers. The novel mechanism can be regarded as an intramolecular version of the intermolecular relay transfer of H-atoms mediated by an external molecule following the general classification of such processes (Catal. Rev.-Sci. Eng. 2014, 56, 403). Here, we present a detailed mechanistic and kinetic analysis of the I-CHAT-facilitated pathways applied to n-hexane, n-heptane, and n-pentane models as prototype molecules for gasoline, diesel, and hybrid rocket fuels. We particularly examined the formation kinetics and subsequent dissociation of the γ-enol-hydroperoxide isomer of the most abundant pentane-derived isomer γ-C5-KHP observed experimentally. To gain molecular-level insight into the I-CHAT catalysis, we have also explored the role of the internal catalyst moieties using truncated models. All applied models demonstrated a significant reduction in the isomerization barriers, primarily due to the decreased ring strain in transition states. In addition, the longer-range and sequential H-migration processes were also identified and illustrated via a combined double keto–enol conversion of heptane-2,6-diketo-4-hydroperoxide as a potential chain-branching model. To assess the possible impact of the I-CHAT channels on global fuel combustion characteristics, we performed a detailed kinetic analysis of the isomerization and decomposition of γ-C5-KHP comparing I-CHAT with key alternative reactions—direct dissociation and Korcek channels. Calculated rate parameters were implemented into a modified version of the n-pentane kinetic model developed earlier using RMG automated model generation tools (ACS Omega, 2023, 8, 4908). Simulations of ignition delay times revealed the significant effect of the new pathways, suggesting an important role of the I-CHAT pathways in the low-T combustion of large alkanes.

Biochemistry & Molecular Biology