Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “calculator”

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

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

VERA-Grizzly Ex-Core Calculations: Watts Bar Unit 1 Cycles 1-2

The critical structures that comprise light-water reactor (LWR) nuclear power plants are subjected to operating environments that can challenge their integrity. Structures in close proximity to the reactor core, such as the reactor pressure vessel (RPV) and the biological shield wall, are subjected to high levels of radiation emanating from the core, as well as elevated temperatures. As the US fleet of operating LWRs ages, the effects of these operating environments on the integrity of these structures must be considered to ensure their continued safe operation. Extending the lifetime of commercial reactors and maintaining the aging reactor fleet require accurate prediction of the exposure of ex-core components to neutron and photon radiation. In particular, concrete degradation studies must be performed to evaluate the safety and long-term operation of reactors with lifetime extensions. The concrete reactor bioshield is important for providing radiological protection during operation and must last for the entire lifetime of the reactor. Recent interest in lifetime extensions furthers the need to accurately simulate concrete material degradation in the reactor bioshield. As a result of this need, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has funded this study to couple its tools, Virtual Environment for Reactor Applications (VERA) and Grizzly. VERA allows users to set up models to calculate time-dependent and fully coupled solutions (with thermal feedback) for ex-core quantities of interest such as vessel and coupon fluence and detector responses for multiple statepoints and cycles. Grizzly is a finite-element application based on the Multiphysics Object Oriented Simulation Environment (MOOSE) framework that is used to enable aging materials calculations. This report highlights the work performed to calculate the fluence in the vessel and concrete for Watts Bar Nuclear Plant Unit 1 (WBN1) Cycles 1 and 2. The fluences obtained from VERA were successfully transferred to Grizzly using a Python script. Four simulations were run with Grizzly: (1) the Mazars model with the initial Young’s modulus being the instantaneous modulus, (2) the Mazars model with the initial Young’s modulus being the delayed modulus, (3) the Mazars model with the initial Young’s modulus being the delayed modulus with the addition of the effects of micro-damage caused by irradiation, and (4) the Mazars model with the initial Young’s modulus being the instantaneous modulus, and with the addition of micro-damage and creep. Details regarding the methods used to obtain the fluence and the statistical errors associated with the VERA Monte Carlo Shift calculations are discussed in greater detail in this report. The results obtained from the four Grizzly models are also presented in this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Sister Rod Destructive Examinations (FY2021) Appendix G: CRIFT Uncertainty Calculations

This document discusses the uncertainty of measurements associated with the CIRFT system and the resulting uncertainty in the acquired measurements and calculated information. The primary parameters measured during the CIRFT test include the deflection magnitude of the specimen, the load applied in bending, and the number of cycles to specimen gross failure. The associated primary calculated parameters are the bending radius/curvature and bending moment. Secondary calculated parameters include specimen stress and strain amplitudes and flexural rigidity. The uncertainty associated with the primary and secondary calculated parameters are investigated based on the known uncertainty in the measured parameters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Extending the Air and Moisture Leakage Calculator to add Residential Buildings and Additional Commercial Buildings

The DOE Windows and Building Envelope Research and Development Roadmap for Emerging Technologies shows that in 2010, infiltration was responsible for 4 quads of space conditioning primary energy use in the residential and commercial sectors. The relative contribution of air leakage in building heating and cooling load is increasing with improvement in the thermal resistance of building envelopes. Advanced air barrier technologies and construction practices have been developed to reduce air leakage in buildings. However, limited information on the impact of air barrier technologies on energy consumption and the durability of buildings has hindered their adoption. In the past Oak Ridge National Laboratory (ORNL), the National Institute of Standards and Technology (NIST), Air Barrier Association of America (ABBA), and U.S.-China Clean Energy Research Center for Building Energy Efficiency (CERC-BEE) collaborated to develop an online calculator that estimates the potential energy and cost savings in major U.S., Canadian and Chinese cities from improvement in air tightness in commercial buildings. In 2018–2019, the calculator was expanded to add moisture transfer calculations given that air leakage through the building envelope can have a significant impact on moisture transfer and associated impacts. In this study, the calculator is expanded further by adding data for two additional commercial buildings (strip mall and primary school) and a residential building. The team investigated the impact of airtightness on energy consumption and moisture transfer of the added buildings. The study includes the analysis of air tightness in 52 major cities in the U.S. and 5 cities in Canada.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗