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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

High-temperature ion-thermal behavior from average-atom calculations

Atom-in-jellium calculations of the Einstein frequency were used to calculate the mean displacement of an ion over a wide range of compression and temperature. Expressed as a fraction of the Wigner-Seitz radius, the displacement is a measure of the asymptotic freedom of the ion at high temperature, and thus of the change in heat capacity from six to three quadratic degrees of freedom per atom. A functional form for free energy was proposed based on the Maxwell-Boltzmann distribution as a correction to the Debye free energy, with a single free parameter representing the effective density of potential modes to be saturated. Finally, this parameter was investigated using molecular dynamics simulations, and found to be ~0.2 per atom. In this way, the ion-thermal contribution can be calculated for a wide-range equation of state (EOS) without requiring a large number of molecular dynamics simulations. Example calculations were performed for carbon, including the sensitivity of key EOS loci to ionic freedom.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

First-principles calculation of the configurational energy density of states for a solid-state ion conductor with a variant of the Wang and Landau algorithm

In this work, a variant of the Wang and Landau algorithm for calculation of the configurational energy density of states is proposed. The algorithm was developed for the purpose of using first-principles simulations, such as density functional theory, to calculate the partition function of disordered sublattices in crystal materials. The expensive calculations of first-principles methods make a parallel algorithm necessary for a practical computation of the configurational energy density of states within a supercell approximation of a solid-state material. The algorithm developed in this work is tested with the two-dimensional (2d) Ising model to bench mark the algorithm and to help provide insight for implementation to a materials science application. Tests with the 2d Ising model revealed that the algorithm has good performance compared to the original Wang and Landau algorithm and the 1/$\textit{t}$ algorithm, in particular the short iteration performance. Further, a proof of convergence is presented within an adiabatic assumption, and the analysis is able to correctly predict the time dependence of the modification factor to the density of states. The algorithm was then applied to the lithium and lanthanum sublattice of the solid-state lithium ion conductor Li 0.5 La 0.5 TiO 3 . This was done to help understand the disordered nature of the lithium and lanthanum. The results find, overall, that the algorithm performs very well for the 2d Ising model and that the results for Li 0.5 La 0.5 TiO 3 are consistent with experiment while providing additional insight into the lithium and lanthanum ordering in the material. The primary result is that the lithium and lanthanum become more mixed between layers along the c axis for increasing temperature. In part, the simulation of the disordered Li 0.5 La 0.5 TiO 3 system serves as a benchmark for what size systems are currently and in the near future practical to calculate with density functional theory methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Static Self-Energy and Effective Mass of the Homogeneous Electron Gas from Quantum Monte Carlo Calculations

We discuss the methodology of quantum Monte Carlo calculations of the effective mass based on the static self-energy Σ(k,0). We then use variational Monte Carlo calculations of Σ(k,0) of the homogeneous electron gas at various densities to obtain results very close to perturbative G 0 W 0 calculations for values of the density parameter 1 ≤ r s ≤ 10. Here, the obtained values for the effective mass are close to diagrammatic Monte Carlo results and disagree with previous quantum Monte Carlo calculations based on a heuristic mapping of excitation energies to those of an ideal gas.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Advanced calculations of X-ray spectroscopies with FEFF10 and Corvus

The real-space Green's function code FEFF has been extensively developed and used for calculations of X-ray and related spectra, including X-ray absorption (XAS), X-ray emission (XES), inelastic X-ray scattering, and electron energy-loss spectra. The code is particularly useful for the analysis and interpretation of the XAS fine-structure (EXAFS) and the near-edge structure (XANES) in materials throughout the periodic table. Nevertheless, many applications, such as non-equilibrium systems, and the analysis of ultra-fast pump–probe experiments, require extensions of the code including finite-temperature and auxiliary calculations of structure and vibrational properties. To enable these extensions, we have developed in tandem a new version FEFF10 and new FEFF -based workflows for the Corvus workflow manager, which allow users to easily augment the capabilities of FEFF10 via auxiliary codes. This coupling facilitates simplified input and automated calculations of spectra based on advanced theoretical techniques. The approach is illustrated with examples of high-temperature behavior, vibrational properties, many-body excitations in XAS, super-heavy materials, and fits of calculated spectra to experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Collection: TD-DFT and EOM-CCSD Calculations for the GDB-9-Ex Dataset

We present two datasets that contain quantum chemical electronic structure calculations for organic molecules from the GDB-9-Ex dataset. The “GDB-9-Ex_TD-DFT-PBE0” dataset contains calculations performed using the time-dependent density functional theory (TD-DFT) first principles method, and the “GDB-9-Ex_EOMCCSD” dataset contains calculations performed using the equation-of-motion coupled cluster (EOM-CCSD) method. Both types of calculations were performed using the ORCA software and provided ultraviolet-visible spectra with a high level of accuracy.

Mehta, Kshitij [Oak Ridge National Laboratory (ORN↗
DOE OSTI2024

Scaling Ensembles of Data-Intensive Quantum Chemical Calculations for Millions of Molecules

Deep learning models are efficient computational tools that can accelerate the inverse design of molecules with desired functional properties by generating predictions at a fraction of the time required by traditional quantum chemical approaches. To ensure that a model maintains accuracy and transferability across broad regions of the chemical space explored during the inverse design, it must be trained on massively large volumes of simulation data. This requires running large-scale ensemble quantum chemical calculations on high-performance computing (HPC) systems for data collection. However, the efficient execution of such large ensemble calculations and the management of large volumes of output data require tools that can judiciously utilize computational resources and manage metadata overhead on the file system. Therefore, we present a high-performance, scalable, ensemble management framework for performing data-intensive quantum chemical electronic structure calculations for organic molecules. This framework provides abstractions to plug different ab initio, first principles, and first principles-based semi-empirical methods and executes them efficiently at large scale on HPC systems. It dynamically distributes tasks to resources and uses tiered storage for managing large collections of files. We employed this framework to process over ten million organic molecules and generate open-source datasets that provide UV-vis absorption spectra by running time-dependent density-functional tight-binding calculations. It is the largest database containing molecular optical spectra that were simulated with quantum chemical methods in a consistent manner.

Mehta, Kshitij↗
DOE OSTI2026

Memory-Aware External Facelist Calculation: A Data-Parallel Atomic Hash Counting Approach

Unstructured volumetric meshes serve as fundamental data representations in various scientific simulations and analyses. They play a crucial role in representing complex computational domains and are essential for important numerical techniques, such as finite element analysis. Whenever such a mesh is read from a file, streamed in-situ, or generated by algorithms, scientific visualization libraries rely on calculating the external surface of a geometry, named “external facelist”, to produce a polygonal mesh for rendering. Consequently, external facelist calculation has become one of the most widely used algorithms in the scientific visualization domain, necessitating optimal performance. In this paper, we explore relevant work on external facelist calculation algorithms in two common visualization libraries, VTK and Viskores, assess their performance and memory constraints, and introduce a novel memory-aware external facelist calculation algorithm employing an atomic hash counting approach. This algorithm fully leverages Viskores' data-parallel primitive operations, facilitating its execution across diverse many-core architectures. Our algorithm features the lowest memory footprint on the GPU and the second-lowest on the CPU among all evaluated methods, and it also delivers the fastest performance on both CPU and GPU. It has been made available under an open-source license in the VTK and Viskores visualization systems.

Tsalikis, Spiros [Kitware] (ORCID:0000000151137195↗
DOE OSTI2025

Ocpp 2.0.1. Interim Kpi Calculator

The project is split into four pieces. The first is a raw OCPP log parser. The second is a file splitter. The third is a message parser. The final piece is the Interim KPI calculator. The OCPP log parser was created from two different formats of raw OCPP 2.0.1 data. Its intended purpose is to extract device IDs and OCPP event messages from nontabular text logs. The parser looks for specific substrings in the logs to identify which of the two "standards" it should select from. The KPI generator does not perform any of its calculations in parallel. Instead, we opt for a naive batching approach. The splitter takes the file generated from the parser and creates many smaller files for each of the device IDs in the dataset. This allows the pandas queries in the log formatter to be iterate over a significantly smaller slice of data, increasing performance significantly. The message parser step takes messages from each of the files (containing distinct device IDs) and breaks the message out into pieces. The final result is a file with different columns specifying different attributes of the JSON message. The file is an aggregation of all different devices. This is the most complex portion of the code. The KPI calculator takes the parsed messages, as a single file, and calculates the KPI from that data. An excel file is produced with four sheets. These contain the metrics for Session Success, Charge Start Success, Charge End Success, and Charge Start Time. It includes the metrics for the different equations in the Interim KPI Implementation Guide as well as a weighted sum of the different equations for each KPI (excluding Charge End Success and Charge Start Time).

Quinn, Casey↗
DOE OSTI2024

Simple Battery Hazard Calculator

SAND2025-00455O Simple Battery Hazard Calculator is a software tool developed in Excel to analyze the electrical hazards of battery systems and determine the required personal protective equipment (PPE) for safe assembly and disassembly. It calculates potential arc flash incident energy, contact thermal hazard level, and recommends PPE based on parameters such as voltage, short circuit current, and over current protection trip time. The calculator uses a recognized method for calculating arc flash incident energy, which was previously inaccessible to most battery workers. 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.

Rosewater, David↗
DOE OSTI2020

A Filon-like integration strategy for calculating exact exchange in periodic boundary conditions: a plane-wave DFT implementation

An efficient and accurate approach for calculating exact exchange and other two-electron integrals has been developed for periodic electronic structure methods. Traditional approaches used for integrating over the Brillouin zone in band structure calculations, e.g. trapezoidal or Monkhorst-Pack, are not accurate enough for two-electron integrals. This is because their integrands contain multiple singularities over the double integration of the Brillouin zone, which with simple integration methods lead to very inaccurate results. A common approach to this problem has been to replace the Coulomb interaction with a screened Coulomb interaction that removes singularities from the integrands in the two-electron integrals, albeit at the inelegance of having to introduce a screening factor which must precomputed or guessed. Instead of introducing screened Coulomb interactions in an ad hoc way, the method developed in this work derives an effective screened potential using a Filon-like integration approach that is based only on the lattice parameters. This approach overcomes the limitations of traditionally defined screened Coulomb interactions for calculating two-electron integrals, and makes chemistry many-body calculations tractable in periodic boundary conditions. This method has been applied to several systems for which conventional DFT methods do not work well, including the reaction pathways for the addition of H 2 to phenol and Au$^{–}_{20}$ nanoparticle, and the electron transfer of a charge trapped state in the Fe(II) containing mica, annite.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Comparative Critical Mass Calculations for NNL and ENDF/B-VIII.0 Zirconium Hydride Thermal Neutron Scattering Laws

Zirconium hydride (ZrH x ) is a moderator material for TRIGA reactors and historical space reactor systems, such as SNAP-10A. Thermal neutron scattering laws (TSL) for two phases of this material, δ and ε, have been previously evaluated by Naval Nuclear Laboratory (NNL) and submitted to the National Nuclear Data Center (NNDC) for inclusion in the US national ENDF/B-VIII.1 nuclear data library. In contrast to the current ENDF/B-VIII.0 TSL evaluations, which consider only a single phase, the new evaluations are derived from separate ab initio calculations for both phases and include coherent elastic effects of the zirconium sublattice. To estimate the impact of these changes to the TSL evaluation of this material, comparative critical mass calculations were performed with MC21 for homogenous mixtures of high- enriched uranium (HEU) and ZrH x in bare and water reflected sphere configurations. These calculations yield an impact on the estimated critical mass as a function of 235 U loading density with maximum differences as large as 1% - 5% for over-moderated thermal spectrum systems. Consequently, the NNL TSL evaluations are anticipated to have a small impact on criticality calculations of thermal reactor systems regardless of the material phase. Nevertheless, characteristic differences exist in the predicted thermal spectra as function of energy for the two sets of TSL evaluations, which are attributed to difference in the underlying phonon density of states of hydrogen bound in ZrH x .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

PyDislocDyn: A Python code for calculating dislocation drag and other crystal properties

PyDislocDyn is a suite of python programs designed to perform various calculations for dislocation dynamics in the continuum limit. In particular, one of its main purposes is to calculate dislocation drag from phonon wind. Additional features include the averaging of elastic constants for polycrystals, the calculation of the dislocation field including its limiting velocities, and the calculation of dislocation self-energy and line tension.

36 MATERIALS SCIENCE↗
DOE OSTI2020

SPND Sensitivity Calculations Using MCNP and Experimental Data from ACRR

The use of the Monte Carlo N-Particle Transport Code (MCNP) to calculate detector sensitivity for Self-Powered Neutron Detectors (SPNDs) in the Annular Core Research Reactor (ACRR) could be a vital tool in the effort to optimize the design of next-generation SPNDs. Next-generation SPND designs, which consider specific materials and geometry, may provide experimenters with capabilities for advanced mixed field dosimetry. These detectors will need to be optimized for configuration, materials, and geometries and the ability to model and iterate must be available in order to decide on the ideal. SPND design. SPNDs were modeled in MCNP which closely resembled the dimensions and location of actual detectors used in the ACRR. Tallies were used to calculate detector sensitivity. Using metrics from a previous report, oscilloscope data from pulses were manipulated in a Matrix Laboratory computing environment (MATLAB) script to calculate experimental detector sensitivity. This report outlines the process in which experimental data from ACRR pulses verified results from tallies in an MCNP ACRR model. The sensitivity values from experiments and MCNP calculations agreed within one standard deviation. Parametric studies were also performed with MCNP to investigate the effects of materials and dimensions of different SPNDs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Acceleration of Thermochimica Calculations in Bison

The performance of the open-source thermochemical software library Thermochimica has been enhanced by developing a re-initialization algorithm to make use of data from previous calculations to reduce the number of convergence steps in the Gibbs energy minimization procedure. This algorithm has been tested in the context of stand-alone Thermochimica calculations, and speedups in the range of 2x-3x achieved for cases with chemistries resembling those of irradiated nuclear fuels. Routines to make use of this re-initialization procedure have been implemented in Bison, in which each node uses the results of the previous calculation at that node as the initial conditions for the following Thermochimica call. Examples based on calculating the diffusion of oxygen in UO$_2$ LWR fuel in 1D and 3D have been tested, with speedups up to 7.38x demonstrated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Connecting Activation Calculations with Gamma Spectroscopy

For many years Lawrence Livermore National Laboratory (LLNL) has used gamma spectroscopy to determine quantities of radioactive nuclides created by irradiation of stable species with neutrons and light charged particles. LLNL has also maintained activation codes and nuclear cross sections for the purpose of modeling activations under the experimental conditions, and the cross section sets used in those calculations typically include nuclear isomers. Connecting the measurements to activation calculations that include isomers is a non-trivial matter because an activation calculation must be followed by a procedure for bookkeeping of calculated isomer populations after the irradiation to simulate what radiochemists would infer about isotope populations from gamma spectroscopy performed days or weeks after the irradiation. Traditionally we have used the method of “isomer sums” for this bookkeeping. This report explains our formulation of isomer sums and application to specific cross section sets for different nuclides that were detected.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Validation of Spent Fuel Dose Rate Calculations Using NDA Measurement Campaign Data

The predicted radiation gamma dose rate is a key parameter to determine if a spent nuclear fuel assembly is self-protected from misuse as well as an indicator for worker/inspector safety. Recent studies showed that using state-of-the-art techniques for dose rate calculations, including the one adopted for this study, the numerical calculations of the gamma dose rates are three times lower than those from the most commonly-cited reference values for spent LWR fuel used by domestic and international agencies. With only a few publically available experimental data points worldwide for this type of benchmark study, the gamma flux measurements from the Spent Fuel Non Destructive Assay (NDA) measurement campaign can be used to help validate the physics approach used in calculating gamma dose rates. Specifically, this data would validate the gamma flux calculation from spent fuel, which is a key component of these multi-physics techniques for dose rate assessments; dose rate measurements were not part of the NDA campaign. These data are primarily the passive gamma measurements of the spent nuclear fuel assemblies from the Swedish Nuclear Fuel and Waste Management Company (SKB).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

A comprehensive final summary of the MCNP calculations performed concerning the UW 137 Cs release event

This document combines and supersedes earlier reports (e.g. LA-UR-19-32358 and LA-UR-20- 22216) that discussed MCNP-calculated conversion factors (count rate per unit activity per unit area) for collimated 2”x2” NaI detectors (Ludlum model 44-10 and Thermo model SPA3) in contact with various surfaces of interest to the University of Washington 137 Cs release event. These surfaces include: Exterior surface of round and rectangular ducts, Exterior surface of outdoor brick work, Exterior surface of cast iron pipes, Exterior surface of an O 2 filter assembly For each surface, the calculated conversion factors allowed inferences on 137 Cs activity on interior surfaces based on exterior contact measurements of the 662 keV photopeak count rate. The MCNP calculations included the effect of attached collimators (i.e. either a Pb/Sb-based Ludlum or a W-based “Rees” model) on detector response. As a consequence, the conversion factors only apply to the appropriate collimated 2”x2” NaI detector operated in pulse height analysis (PHA) mode. Other assumptions affecting the calculations and valid application of the conversion factors will be outlined below.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

System Model Calculations for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington

Two lined trenches (Trench 31 and 34) located within the 200 West Low-Level Burial Ground (LLBG) area of the U.S. Department of Energy’s (DOE) Hanford Site are designated for permanent disposal of low-level radioactive waste (LLW) and mixed low-level radioactive waste (MLLW). In accordance with DOE O 435.1, Radioactive Waste Management, radioactive waste shall be managed and disposed in a manner that is protective of worker, public health and safety, and the environment. DOE O 435.1 requires a site-specific radiological performance assessment (PA) that includes calculations of potential releases and subsequent doses to members of the public for a period of 1,000 years after closure of a low-level waste (LLW) disposal facility. The purpose of this environmental calculation file (ECF) is to document the various dose calculations performed to support DOE’s LLBG PA. The dose assessments are performed to evaluate the potential exposure of an all-pathway representative person to radionuclide contaminants of potential concern (COPCs) that may be released from the LLBG Trenches 31 and 34 to the point of calculation (POC) located at the outer edge of a 100 m buffer zone surrounding the trench boundary. This ECF uses inputs from other ECFs and their associated model package reports (MPRs) supporting the LLBG PA. The five major objectives of this ECF are to present the results of the following: (1) Deterministic dose assessments for both groundwater and atmospheric exposure scenarios; (2) Radon-222 (Rn-222) flux emanating at the surface from the disposed wastes; (3) Uncertainty analyses of the groundwater pathway annual dose due to uncertainty in input parameters; (4) Sensitivity analyses of the input parameters that could potentially impact the fate and transport of the contaminants and doses; (5) Doses arising from inadvertent intruder exposure scenarios.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗