Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical relativity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quantum stabilization of unexpected ordered phases on the honeycomb lattice

In this work, I discuss and showcase the utility of the method called minimally-augmented spin-wave theory (MAGSWT) as a relatively simple semi-analytical approach to the phase diagrams of quantum magnets. It complements numerical methods by providing physical insight into which states are competitive and by yielding approximate phase boundaries that agree well with much more numerically intensive calculations. This approach enabled me to construct the faithful phase diagrams of two paradigmatic honeycomb-lattice models in the quantum S=1/2 limit, the J 1 –J 3 FM-AF and J 1 –J 2 AF models, for the collinear quantum phases that replace or extend the classical ones. The results are in good qualitative and semi-quantitative agreement with state-of-the-art numerical studies for these models, correctly capturing the emergence of unexpected quantum phases and the suppression of classically favored spiral orders by quantum fluctuations. This study provides a much-needed important guidance to the ongoing theoretical and experimental searches of the unconventional quantum states. This work will be of significant and timely interest to both theorists and experimentalists in the field of quantum magnetism, broadly defined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Accident Progression and Source Term Demonstration Calculations for the Advanced Burner Test Reactor (ABTR)

As advanced reactor designers approach broad deployment of new reactor technologies and the associated applications, a critical milestone on this path includes licensing or authorization of these designs in the current regulatory environment. Commercial licensing and DOE authorization both consist of numerous aspects to not only certify the designed plant will operate as intended, but also related to confirmation of the safe operation of these facilities that comply with specific standards and requirements. To that end, domestic industry and regulators require toolsets and methods which can be used to confirm that advanced reactor designs not only align with expected performance during normal operation, but also a series of transient scenarios. These codes require a satisfactory validation basis, and in the absence of such a basis, code-to-code comparisons are often utilized as a surrogate to confirm model adequacy. To support this need for fast reactor technologies, a series of analyses using SAS4A/SASSYS- 1 and the Simplified Radionuclide Transport (SRT) code have been conducted for the Advanced Burner Test Reactor (ABTR). These analyses are expected to not only provide a baseline for transient and source term analyses of ABTR, but to also provide context for the key physics that determine evolution of specific transient scenarios in a sodium-cooled fast reactor (SFR). These codes are utilized for this work given their adoption in the SFR industry and by key SFR stakeholders.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Simulating a numerical UV completion of quartic Galileons

The Galileon theory is a prototypical effective field theory that incorporates the Vainshtein screening mechanism—a feature that arises in some extensions of general relativity, such as massive gravity. The Vainshtein effect requires that the theory contain higher order derivative interactions, which results in Galileons, and theories like them, failing to be technically well posed. While this is not a fundamental issue when the theory is correctly treated as an effective field theory, it nevertheless poses significant practical problems when numerically simulating this model. These problems can be tamed using a number of different approaches: introducing an active low-pass filter and/or constructing a UV completion at the level of the equations of motion, which controls the high momentum modes. These methods have been tested on cubic Galileon interactions, and have been shown to reproduce the correct low-energy behavior. Here we show how the numerical UV-completion method can be applied to quartic Galileon interactions, and present the first simulations of the quartic Galileon model using this technique. We demonstrate that our approach can probe physics in the regime of the effective field theory in which the quartic term dominates, while successfully reproducing the known results for cubic interactions. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2026

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗
DOE OSTI2025

Two-Step Procedure to Detect Cosmological Gravitational Wave Backgrounds with Next-Generation Terrestrial Gravitational-Wave Detectors

Cosmological gravitational-wave backgrounds are an exciting science target for next-generation ground-based detectors, as they encode invaluable information about the primordial Universe. However, any such background is expected to be obscured by the astrophysical foreground from compact-binary coalescences. We propose a novel framework to detect a cosmological gravitational-wave background in the presence of binary black holes and binary neutron star signals with next-generation ground-based detectors, including Cosmic Explorer and the Einstein Telescope. Our procedure involves first removing all the individually resolved binary black hole signals by notching them out in the time-frequency domain. Then, we perform joint Bayesian inference on the individually resolved binary neutron star signals, the unresolved binary neutron star foreground, and the cosmological background. For a flat cosmological background, we find that we can claim detection at 5⁢𝜎 level when Ω ref ≥ 2.7 × 10 −12 /$\sqrt{𝑇_{obs}/yr}$, where 𝑇 obs is the observation time (in years), which is within a factor of ≲ 2 from the sensitivity reached in the absence of these astrophysical foregrounds.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Fundamental absorption bandwidth to thickness limit for transparent homogeneous layers

Abstract Past work has considered the analytic properties of the reflection coefficient for a metal-backed slab. The primary result established a fundamental relationship for the minimal layer thickness to bandwidth ratio achievable for an absorber. There has yet to be establishment of a similar relationship for non-metal-backed layers, and here we present the universal result based on the Kramers–Kronig relations. Our theory is validated with transfer matrix calculations of homogeneous materials, and full-wave numerical simulations of electromagnetic metamaterials. Our results place more general fundamental limits on absorbers and thus will be important for both fundamental and applied studies.

Padilla, Willie J.↗
DOE OSTI2024

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Decoding the formation of hammerhead ion populations observed by Parker Solar Probe

Context. In situ observations by the Parker Solar Probe (PSP) have revealed new properties of the proton velocity distributions (VDs), including hammerhead features that suggest a non-isotropic broadening of the beams. Aims. The present work proposes a very plausible explanation for the formation of hammerhead proton populations through the action of a proton firehose-like instability triggered by the proton beam. Methods. We investigated a self-generated firehose-like instability driven by the relative drift of ion populations using a simplified moment-based quasi-linear (QL) theory. While simpler and faster than advanced numerical simulations, this toy model provided rapid insights and concisely highlighted the role of plasma micro-instabilities in relaxing the observed anisotropies of particle VDs in the solar wind and space plasmas. Results. The QL theory proposed here shows that the resulting transverse waves are right-hand polarized and have two consequences on the protons: (i) They reduce the relative drift between the beam and the core, but above all, (ii) they induce a strong perpendicular temperature anisotropy specific to the observed hammerhead ion beam. Moreover, the long-run QL results suggest that these hammerhead distributions are rather transitory states that are still subject to relaxation mechanisms, in which instabilities such as the one discussed here are very likely involved. This foundational work motivates future detailed studies using advanced methods.

Shaaban, Shaaban M. (ORCID:000000030465598X)↗
DOE OSTI2025

Game theoretic modeling and optimization of competition and collaboration in dual channel electronic waste supply chains

The rapid growth of electronic waste (e-waste) presents critical challenges for sustainable resource recovery and environmental protection. This study develops a dual-channel closed-loop supply chain (CLSC) model formulated as a hierarchical Stackelberg game, that integrates dynamic pricing and cost-sharing mechanisms to optimize both economic and environmental outcomes. The model explicitly captures strategic interactions between manufacturer-led and third-party recycling channels, accounting for consumer behavior, regulatory incentives, and market competition. Numerical simulations conducted (implemented over a four-iteration horizon using a commercial optimization solver) show that, relative to the baseline equilibrium, manufacturer profit increases from 11.6 thousand USD to 37.9 thousand USD (+226.8%), total recycled volume rises from 7,848 to 7,942 units (+1.2%), and collector profit nearly doubles under cost-sharing, enabling more equitable profit distribution. Furthermore, scenario-based simulations across Sub-Saharan Africa, high-income economies, and emerging Asian industrial countries reveal that infrastructure quality, policy intensity, and labor costs critically shape recycling efficiency and profit allocation. These findings demonstrate that subsidies alone are insufficient to ensure system efficiency. Instead, coordinated strategies that integrate internal incentive alignment with context-sensitive policy support are required. Overall, this study offers a robust framework for designing resilient, efficient, and regionally adaptable e-waste management systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Tailored ultrasound propagation in microscale metamaterials via inertia design

The quasi-static properties of micro-architected (meta)materials have been extensively studied over the past decade, but their dynamic responses, especially in acoustic metamaterials with engineered wave propagation behavior, represent a new frontier. However, challenges in miniaturizing and characterizing acoustic metamaterials in high-frequency (megahertz) regimes have hindered progress toward experimentally implementing ultrasonic-wave control. Here, we present an inertia design framework based on positioning microspheres to tune responses of 3D microscale metamaterials. We demonstrate tunable quasi-static stiffness by up to 75% and dynamic longitudinal-wave velocities by up to 25% while maintaining identical material density. Using noncontact laser-based dynamic experiments of tunable elastodynamic properties and numerical demonstrations of spatio-temporal ultrasound wave propagation, we explore the tunable static and elastodynamic property relation. This design framework expands the quasi-static and dynamic metamaterial property space through simple geometric changes, enabling facile design and fabrication of metamaterials for applications in medical ultrasound and analog computing.

42 ENGINEERING↗
DOE OSTI2024

A Global Increase in Nearshore Tropical Cyclone Intensification

Abstract Tropical Cyclones (TCs) inflict substantial coastal damages, making it pertinent to understand changing storm characteristics in the important nearshore region. Past work examined several aspects of TCs relevant for impacts in coastal regions. However, few studies explored nearshore storm intensification and its response to climate change at the global scale. Here, we address this using a suite of observations and numerical model simulations. Over the historical period 1979–2020, observations reveal a global mean TC intensification rate increase of about 3 kt per 24‐hr in regions close to the coast. Analysis of the observed large‐scale environment shows that stronger decreases in vertical wind shear and larger increases in relative humidity relative to the open oceans are responsible. Further, high‐resolution climate model simulations suggest that nearshore TC intensification will continue to rise under global warming. Idealized numerical experiments with an intermediate complexity model reveal that decreasing shear near coastlines, driven by amplified warming in the upper troposphere and changes in heating patterns, is the major pathway for these projected increases in nearshore TC intensification.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Extreme radiation emission regime for electron beams in strong focusing ion channels and undulators

A fundamental comparison between a magnetic undulator and an ion channel, or betatron, radiation from relativistic electrons is presented. While conventional theories nominally range from the undulator (𝐾 <1) to the wiggler (𝐾 >1) regime, they are only applicable for sufficiently large Lorentz factors (𝛾 0 ≫𝐾). They therefore do not account for high 𝐾/𝛾 0 cases, for which we show that particle trajectories and radiation characteristics strongly deviate from the linear predictions in both magnetic undulators and ion channels. This problem arises from the fundamental differences between a magnetostatically and electrostatically induced oscillation. A reformulation of both the ion channel betatron wavelength and amplitude, as well as the same parameters in a magnetic undulator, permits us to compare cases with equivalent oscillation period and amplitude in the two different scenarios. The notable differences in spectral features of the two radiation mechanisms can then be addressed via numerical simulations of single particle as well as full beam dynamics. Additionally, we identify and quantify a novel transverse orbit precession effect in ion channels for particles with initial angular momentum relative to the device axis. This effect, which is significant in cases of strong transverse kinetic energy oscillations, alters both the radiation divergence and the beam emittance. In this paper, we present this new theoretical framework and compare its results with numerical simulation applied to realizable experimental tests of such radiation sources.

Frazzitta, Andrea [University of Rome “La Sapienza↗
DOE OSTI2024

HDBind: encoding of molecular structure with hyperdimensional binary representations

Traditional methods for identifying “hit” molecules from a large collection of potential drug-like candidates rely on biophysical theory to compute approximations to the Gibbs free energy of the binding interaction between the drug and its protein target. These approaches have a significant limitation in that they require exceptional computing capabilities for even relatively small collections of molecules. Increasingly large and complex state-of-the-art deep learning approaches have gained popularity with the promise to improve the productivity of drug design, notorious for its numerous failures. However, as deep learning models increase in their size and complexity, their acceleration at the hardware level becomes more challenging. Hyperdimensional Computing (HDC) has recently gained attention in the computer hardware community due to its algorithmic simplicity relative to deep learning approaches. The HDC learning paradigm, which represents data with high-dimension binary vectors, allows the use of low-precision binary vector arithmetic to create models of the data that can be learned without the need for the gradient-based optimization required in many conventional machine learning and deep learning methods. This algorithmic simplicity allows for acceleration in hardware that has been previously demonstrated in a range of application areas (computer vision, bioinformatics, mass spectrometery, remote sensing, edge devices, etc.). To the best of our knowledge, our work is the first to consider HDC for the task of fast and efficient screening of modern drug-like compound libraries. We also propose the first HDC graph-based encoding methods for molecular data, demonstrating consistent and substantial improvement over previous work. We compare our approaches to alternative approaches on the well-studied MoleculeNet dataset and the recently proposed LIT-PCBA dataset derived from high quality PubChem assays. We demonstrate our methods on multiple target hardware platforms, including Graphics Processing Units (GPUs) and Field Programmable Gate Arrays (FPGAs), showing at least an order of magnitude improvement in energy efficiency versus even our smallest neural network baseline model with a single hidden layer. Our work thus motivates further investigation into molecular representation learning to develop ultra-efficient pre-screening tools. We make our code publicly available at https://github.com/LLNL/hdbind.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Accelerating magnonic simulations with the pseudospectral Landau-Lifshitz equation

The pseudospectral Landau-Lifshitz (PS-LL) model can describe atomic-scale magnetic exchange interactions within a continuum framework. This is achieved by employing a convolution kernel that models the nonlocal interaction in a grid-independent manner. Even though the PS-LL was originally introduced to address atomic exchange, any nonlocal kernel can be modeled. In the field of magnonics, the dipole field is fundamental to describe the dispersion relation of magnons, the quasiparticle representation of angular momentum. Because dipole-dipole interactions are long-range, numerical approaches typically rely on convolutions. Here, we demonstrate that the PS-LL model can be used to perform magnonic simulations with a single convolution kernel derived from analytical solutions. We demonstrate a twofold increase in computational speed compared with the full dipole calculation. This approach is valid insofar as the excitations are linear, which is typically the case for magnons. Our results have the potential to accelerate magnonic research, particularly for the inverse design method, where several simulations must be performed to achieve the desired outcome.

Mathematics and computing↗
DOE OSTI2024

A Parametric Study of the SASI Comparing General Relativistic and Nonrelativistic Treatments

We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics with GR gravity, and nonrelativistic (NR) hydrodynamics with Newtonian gravity, in an idealized model setting, we vary the initial radius of the shock, and by varying its mass and radius in concert, the proto-neutron star compactness. We investigate four compactnesses expected in a post-bounce core-collapse supernova (CCSN). We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the highest-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the highest-compactness suite. We discuss implications of our results for CCSN simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Plant Defense Proteins: Recent Discoveries and Applications

Proteins play pivotal roles in safeguarding plants against numerous biotic and abiotic stresses. Understanding their biological functions and mechanisms of action is essential for advancing plant biology, agriculture, and biotechnology. This review considers the diversity and potential applications of plant defense proteins including pathogenesis-related (PR) proteins, chitinases, glucanases, protease inhibitors, lectins, and antimicrobial peptides. Recent advances, such as the omics technologies, have enabled the discovery of new plant defense proteins and regulatory networks that govern plant defense responses and unveiled numerous roles of plant defense proteins in stress perception, signal transduction, and immune priming. The molecular affinities and enzymatic activities of plant defense proteins are essential for their defense functions. Applications of plant defense proteins span agriculture, biotechnology, and medicine, including the development of resistant crop varieties, bio-based products, biopharmaceuticals, and functional foods. Future research directions include elucidating the structural bases of defense protein functions, exploring protein interactions with ligands and other proteins, and engineering defense proteins for enhanced efficacy. Overall, this review illuminates the significance of plant defense proteins against biotic stresses in plant biology and biotechnology, emphasizing their potential for sustainable agriculture and environmental management.

G-proteins↗
DOE OSTI2025

Numerical design considerations for vapor transport deposition of metal-halide perovskite thin films

While metal-halide perovskites (MHPs) offer high efficiency and potential application in single junction and tandem solar cells, challenges remain in translating solution-based methods to commercial-scale production. Vapor transport deposition (VTD) offers advantages such as precise control over film composition and purity, as well as the possibility of achieving high thin film deposition rates for scale-up. In this work, a numerical model is developed to simulate the deposition process, analyzing the impact of parameters including source material temperature, carrier gas and dilution gas flow rates, and chamber pressure to optimize deposition rates and film uniformity. Geometrical parameters, such as throw distance, are also found to significantly affect rate and uniformity. Furthermore, the modeling results presented here offer valuable insight for improving the scalability of VTD for the processing of metal-halide perovskites.

Chemical vapor deposition↗