Search NASA⌕ Search

SEARCH · Search NASA

Results for “Calculation”

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 127 records · Page 7

Perturbative treatment of nonlocal chiral interactions in auxiliary-field diffusion Monte Carlo calculations

Nuclear many-body systems, ranging from nuclei to neutron stars, are some of the most interesting physical phenomena in our universe, and quantum Monte Carlo (QMC) approaches are among the most accurate many-body methods currently available to study them. In recent decades, interactions derived from chiral effective field theory (EFT) have been widely adopted in the study of nuclear many-body systems. One drawback of the QMC approach is the requirement that the nuclear interactions need to be local, whereas chiral EFT interactions usually contain nonlocalities. In this work, we leverage the capability of computing second-order perturbative corrections to the ground-state energy in order to develop a self-consistent approach to including nonlocal operators in QMC calculations. In conclusion, we investigate both the deuteron and the neutron-matter equation of state in order to show the robustness of our technique and pave the way for future QMC calculations at higher orders in the EFT, where nonlocal operators cannot be avoided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Constraining calculations of Pu 239 ( n , γ ) cross sections with multistep cascade γ -ray spectra measured at the DANCE facility

Measurements of multistep cascade gamma-ray spectra from several low-lying 1 + resonances in the Pu 239 ( n , γ ) reaction were made using the DANCE detector. The measured spectra are compared to spectra simulated using several models for the photon strength function and the nuclear level density, representing the first analysis of this nature on a fissile nucleus. In general, models of the photon strength function that included a strong M 1 “scissors mode” contribution coupled with a Gilbert-Cameron nuclear level density provided a reasonable description of the spectra. Simulations using photon strength functions and nuclear level densities specifically for Pu 240 , extracted by the Oslo method, resulted in a very good representation of the measured spectra. Neutron-capture cross section calculations that included the scissors mode provided a good estimate of the previously measured values, without additional normalization. However, the scissors-mode energy that produced the best description of the multistep cascade spectra was different from the value that best reproduced the capture cross section. In addition, the 〈 Γ γ 〉 inferred from the multistep cascade calculations was somewhat less than the tabulated value, which was determined from resonance analysis. Published by the American Physical Society 2025

Ullmann, J. L. (ORCID:0000000207784405)↗

Fast and accurate algorithm for calculating long-baseline neutrino oscillation probabilities with matter effects

Neutrino oscillation experiments will be entering the precision era in the next decade with the advent of high statistics experiments like DUNE, HK, and JUNO. Correctly estimating the confidence intervals from data for the oscillation parameters requires very large Monte Carlo datasets involving calculating the oscillation probabilities in matter many, many times. In this paper, we leverage past work to present a new, fast, precise technique for calculating neutrino oscillation probabilities in matter optimized for long-baseline neutrino oscillations in the Earth’s crust including both accelerator and reactor experiments. For ease of use by theorists and experimentalists, we provide fast ++ and codes . Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Range-separated thermal hybrid exchange-correlation density functional for accurate band-gap calculations of warm dense matter

Here, in this study, we present a range-separated thermal hybrid exchange-correlation (XC) functional, termed RS-KDT0, which provides highly accurate band-gap predictions at elevated temperature while reducing the computational cost relative to thermal global hybrid functionals. At zero temperature, RS-KDT0 recovers the widely used ground-state range-separated hybrid functional HSE06. To assess its performance across temperature regimes, we conducted static calculations of electronic pressure, band gaps, and band structures. RS-KDT0 exhibits the same level of accuracy as HSE06 in the low-temperature regime and demonstrates improved agreement with finite-temperature GW calculations at elevated temperatures when compared to other thermal XC functionals such as the generalized gradient approximation (GGA)-level KDT16 and the global hybrid KDT0.

density functional theory↗

Electromagnetic radii of light nuclei from variational Monte Carlo calculations

We present variational Monte Carlo calculations of charge and magnetic radii in 𝐴 ≤ 10 nuclei. The calculations are based on the Norfolk two- and three-nucleon interactions, and associated one- and two-nucleon electromagnetic charge and current operators derived up to next-to-next-to-next-to-leading order in the chiral expansion. The charge and magnetic radii are extracted from the respective form factors. We find that the charge radii are within 5% of the experimental values for the nuclei considered. For the magnetic radii, a comparison is available only with 3 H and 3 ⁢He electron scattering data that are affected by large error bars. We hope that our predictions foster an interest in precisely measuring magnetic radii of heavier systems.

lepton induced nuclear reactions↗

Connecting relativistic density functional theory to microscopic calculations

The development of systematic effective field theories (EFTs) for nuclear forces and advances in solving the nuclear many-body problem have greatly improved our understanding of dense nuclear matter and the structure of finite nuclei. For global nuclear calculations, density functional theories (DFTs) have been developed to reduce the complexity and computational cost required in describing nuclear systems. However, DFT often makes approximations and assumptions about terms included in the functional, which may introduce systematic uncertainties compared to microscopic calculations using EFTs. In this work, we investigate possible avenues of improving nuclear DFT using nonlinear relativistic mean-field (RMF) theory. We explore the impact of RMF model extensions by fitting the nonlinear RMF model to predictions of nuclear matter and selected closed-shell nuclei using four successful chiral EFT Hamiltonians. We find that these model extensions are impactful and important in capturing the physics present within chiral Hamiltonians, particularly for charge radii and neutron skins of closed-shell nuclei. However, there are additional effects that are not captured within the RMF model, particularly within the isoscalar sector of RMF theory. Additional model extensions and the reliability of the nonlinear RMF model are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development of a multiphase equation of state for gallium with experiments and ab initio free-energy calculations

Here, we present a five-phase equation of state (EOS) for elemental gallium (Ga) that is developed using both experimental data and new theoretical predictions. Four experimentally observed solid phases (Ga-I, Ga-II, Ga-III, and Ga-IV) and one liquid phase are included. To improve our understanding of the thermal behavior of Ga and its phase boundaries under compression, we have performed ab initio density functional theory (DFT) free-energy calculations for the Ga-III and liquid phases, which enables us to determine the melt temperature to be 2214 ± 100 K at 110 GPa, extending significantly beyond the existing experimental melt data, which are limited to 25 GPa. In order to best describe the electron-thermal contribution, which dominates the liquid free energy at high temperatures, we have carried out averaged-atom-in-jellium DFT calculations to cover the entire temperature and density range of the EOS. The resulting multiphase Ga EOS is able to accurately reproduce a diverse variety of data, including known phase boundaries, the principal Hugoniot, low-pressure liquid isobars, and diamond-anvil-cell isotherm measurements at high pressures. It agrees more closely with key experimentally measured properties than other Ga EOS models targeted for high-pressure applications.

Wu, Christine J. [Lawrence Livermore National Labo↗

Calculating topological properties of artificial graphene in B-field

SAND2025-03258O Calculating topological properties of artificial graphene in B-field is a user-friendly tool designed to analyze artificial graphene systems influenced by magnetic fields. It helps researchers understand the unique properties of these materials by calculating the local Chern marker, a key indicator of their behavior. With just one file, this software runs easily on any device with Matlab, making it accessible for scientists and engineers. It provides valuable insights into the electronic characteristics of artificial graphene, which supports advancements in material science and technology, paving the way for innovative applications in electronics and beyond. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Spataru, Dan [Sandia National Lab. (SNL-CA), Liver↗

Proton pulse charge calculation algorithm in Beam Power Limiting System at Spallation Neutron Source

A proton pulse charge calculation algorithm in the Beam Power Limiting System (BPLS) at the Spallation Neutron Source (SNS) was developed and implemented in an FPGA. The algorithm calculates one-minute running average of the pulse charges and issues a fault to the Personal Protection System (PPS) and the Machine Protection System (MPS) when a limit is reached. A bit-accurate model of the algorithm was first developed and tested in Matlab® and then implemented and simulated in VHDL using Vivado® design environment. Finally, the algorithm was verified on a µTCA-based hardware platform.

Bobrek, Miljko [ORNL] (ORCID:0000000332763451)↗

Calculating Reactivity and Fusion Energy Yield

During my time at Lawrence Livermore National Lab, I had the opportunity to intern at the HEDS (High Energy Density Science) Center and become an IFE (Inertial Fusion Energy) summer student for three weeks with Dr. Veronika Kruse as my mentor. Within this time frame, my research project’s goal was to develop an understanding of nuclear fusion while building technical skills in calculating and analyzing the required data. The following report outlines the process of calculating reactivity curves and energy yields of specific fusion reactions. For context, nuclear fusion is when two atomic nuclei combine, and energy is released. NIF (the National Ignition Facility), located at LLNL, uses laser-based inertial confinement fusion to achieve ignition. To simplify the process, NIF uses lasers to heat a small sphere or capsule, made up of primarily deuterium and tritium, to extremely high temperatures and densities to achieve fusion ignition, where more energy is released than the lasers introduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of Increased Monte Carlo Parameters on Sensitivity Calculations with SCALE [Slides]

For both models and tests, NPG and NSK parameters have only small effects on calculated sensitivity coefficients. Outside of NPG=100, only differences in CFP affected sensitivity coefficient values. Fission reactions ( 235 U) require more NPG values than scattering ( 238 U) reactions – more particles are needed to locate fission sources in the model more accurately. This work confirms the previous results with the IFP method where the CFP parameter has the greatest impact on calculated sensitivity coefficients. While immediate work focuses on fast systems, other model specifications may require a different set of MC parameters.

97 MATHEMATICS AND COMPUTING↗

Measured vs. Calculated Dose Rates During Ring Injection Dump Replacement

One of the most irradiated Spallation Neutron Source (SNS) accelerator components is the ring injection dump (RID), which is located downstream of the accumulator ring injection section. The unstripped and only partially stripped H - beam, about 5% of the full SNS beam current, is discarded in the injection into the accumulator ring and is redirected to the RID. According to the accelerator operation plan, the beam stop and the window assemblies of the existing RID are removed and replaced when they have reached their end-of-life owing to radiation exposure. This procedure, which includes multiple steps, took place during a facility maintenance period in March 2023. To support work planning and meet as low as reasonably achievable (ALARA) requirements during removal and exchange of the components, dose rates were calculated for each stage of the exchange operation. During each stage, dose rates were measured before and after work. This comparison shows that measured dose rates are within a factor of 2 of the predicted (calculated) dose rates.

43 PARTICLE ACCELERATORS↗

Design Heating and Cooling Load Calculation Versus Building Load Simulation for Cold Climate Heat Pumps: Understanding the "Gap"

This report explores the differences between Manual J-equivalent block load calculations and building HVAC energy simulation results using EnergyPlus™ calculations when designing cold climate heat pump systems for residential use. This study will help HVAC researchers and advanced designers understand the impacts of oversizing heat pumps on home energy use.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

High-Temperature Gas Sensor Materials with Properties Predicted via First-Principles Calculations with Machine Learning Modeling and Experimental Corroboration

Understanding the temperature dependence of functional properties of sensing materials is vital for their applications in combustion environments. The electron-phonon coupling that derives the electronic structure change with temperatures is a key property of interest as it affects other sensing responses. Herein, we first assess the temperature dependence of band gap renormalization in sensing materials by employing Allen-Heine-Cardona (AHC) theory with density functional theory (DFT) simulations corroborated with experimental observation. As the AHC calculations are impractical for high-throughput screening of materials, we employ data-driven Gaussian process regression to predict the parameters employed in the O’Donnell empirical model from a set of physical features. To mitigate the reliability issues arising from the small size of the dataset, we apply a Bayesian technique to improve the generalizability of the data-driven models as well as to quantify the uncertainty associated with theoretical predictions. These models capture well the overall trend of the O’Donnell parameters with respect to a reduced feature set obtained by transforming the available physical features. Quantifying the associated uncertainty helps us understand the reliability of the predictions and, therefore, the variation of bandgap as a function of temperature for other novel materials. The predicted candidates from machine learning models are further validated by experiments and DFT calculations.

bandgap renormalization↗

Real time heat load calculation software based on EPICS for Fermilab PIP-II CM tests

Fermilab has a project to improve the proton beam energy which is called PIP-II (the 2nd Proton Improvement Plan). There is a superconducting linear accelerator, LINAC, to improve the proton beam power and the LINAC consists of 5 types of cryomodules (CM), 1 HWR CM, 2 SSR1 CM, 4 SSR2 CM, LB650 CM, and HB650 CM. The prototypes of these cryomodules are being tested at Fermilab’s CryoModule Test Facility (CMTF). Heat load measurements are an important part of the prototype CM testing. The CMTF cryogenic control system was developed based on the ACNET (Accelerator Control NETwork) for CM testing for other projects, but the PIP-II cryogenic control system will be implemented using the Experimental Physics and Industrial Control System (EPICS). As part of the prototype CM testing campaign an EPICS based control system has been implemented at CMTF. This EPICS cryogenic control system includes real time heat load calculation software utilizing the Fortran implementation of Hepak. This paper details the real time heat load calculation software developed for the prototype CM testing including the first results from the HB 650 CM.

Yoon, S. [Fermilab]↗

Benchmark Calculation for Turkey Point Unit 3 Cycles 1-3 Using the SCALE 6.3/Polaris–PARCS v3.4.2 Code Package

Benchmark calculations were performed for Turkey Point Unit 3 cycles 1–3 to validate the SCALE 6.3/Polaris–PARCS v3.4.2 code with the ENDF/B–VII.1 56–group library by comparing the simulated results with the measured data. The benchmark results will be used in evaluating the SCALE/Polaris–PARCS code package’s uncertainties for pressurized water reactor physics analysis. That future analysis will include key nuclear parameters such as reactivity, control bank worth, temperature coefficients, and pin and assembly power peaking factors. The present document details plant and fuel design specifications and input data for SCALE/Polaris, GenPMAXS, and PARCS. Additional details are provided with respect to the input and output files produced for the benchmark calculations. The benchmark results are summarized such that they can be used in evaluating uncertainties with other benchmark results for key nuclear parameters.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Benchmark Calculation for Surry Unit 1 Cycles 1-3 Using the SCALE 6.3/Polaris–PARCS v3.4.2 Code Package

The benchmark calculations were performed for Surry Unit 1 cycles 1–3 to validate the SCALE 6.3/Polaris–Purdue Advanced Reactor Core Simulator (PARCS) v3.4.2 with the ENDF/B–VII.1 56–group library by comparing the simulated results with the measured data. The benchmark results will be used to evaluate uncertainties of the SCALE/Polaris–PARCS code package for pressurized water reactor physics analysis for key nuclear parameters such as reactivity, control bank worth, temperature coefficients, and pin and assembly power peaking factors. This report details plant and fuel design specifications and input data for SCALE/Polaris, GenPMAXS, and PARCS. Additional details are provided for the input and output files produced for the benchmark calculations. The benchmark results were summarized such that they can be used in evaluating uncertainties with other benchmark results for key nuclear parameters.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Benchmark Calculation for the Quad Cities Unit 1 Cycles 1-3 Using the SCALE 6.3/Polaris–PARCS v3.4.2 Code Package

In this study, benchmark calculations were performed for the Quad Cities Unit 1 cycles 1–3 to validate the SCALE 6.3/Polaris–PARCS v3.4.2 code package with the ENDF/B-VII.1 AMPX 56-group library by comparing the simulated results with the measured data. The benchmark results will be used in evaluating uncertainties of the SCALE/Polaris–PARCS code package for boiling water reactor physics analysis for key nuclear parameters such as reactivity and assembly power peaking factors. This report details plant and fuel design specifications and input data for SCALE/Polaris, GenPMAXS, and PARCS; additionally, detailed information is provided for all the input and output files produced for the benchmark calculations. The benchmark results are summarized herein so that they can be used to evaluate uncertainties with other benchmark results for key nuclear parameters.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗