Search NASA⌕ Search

SEARCH · Search NASA

Results for “calculation method”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Equilipy

There is high demand for high-throughput calculations of phase equilibria based on the CALPHAD method (CALculation of PHAse Diagram). High-throughput calculations are possible through HPC, currently by a commercial program called Thermo-Calc. However, the number of nodes/processors are limited to the number of purchased license (16 processors per license). This program provides a toolkit for high-throughput calculations of phase equilibria by the CALPHAD method. The program is prepared for Python environment, so that it can be easily installed and used together with other open-source programs. The program can be run in supercomputer.

Kwon, Sunyong [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2021

A Data-Driven Global Sensitivity Analysis Framework for Three-Phase Distribution System with PVs

Global sensitivity analysis (GSA) of distribution systems with respect to stochastic PV and load variations plays an important role in designing optimal voltage control schemes. This paper proposes a data-driven framework for GSA of distribution systems. In particular, two representative surrogate modeling-based approaches are developed, including the traditional Gaussian process-based and the analysis of variance (ANOVA) kernel ones. The key idea is to develop a surrogate model that captures the hidden global relationship between voltage and real and reactive power injections from the historical data. With the surrogate model, the Sobol indices can be conveniently calculated through either the sampling-based method or the analytical method to assess the global sensitivity of voltage to variations of PV and load power injections. The sampling-based method approximates the Sobol indices using Monte Carlo simulations while the analytical method calculates them by resorting to the ANOVA expansion framework. Comparison results with other model-based GSA methods on the unbalanced three-phase IEEE 37-bus and 123-bus distribution systems show that the proposed framework can achieve much higher computational efficiency with negligible loss of accuracy. The results on a real 240-node distribution system using actual smart meter data further validate the feasibility and scalability of the proposed framework.

14 SOLAR ENERGY↗
DOE OSTI2021

BATPAC--VERSION 5.0

Argonne National Laboratory has worked on electrochemical energy storage for several decades. The focus on lithium-ion chemistries started in the early 1990s, developing new materials, synthesis methods, and performance characterizations. Sponsored by the U.S. Department of Energy, Energy Efficiency and Renewable Energy, Vehicle Technologies Office (DOE-EERE-VTO), Argonne has led with many advances. The experimental activities were complemented with multi-scale modeling that ranged from the atomic to the system (manufacturing processes and automobiles) level. The Battery Performance and Cost (BatPaC) model is a calculation method based on Microsoft¿ Office Excel spreadsheets that have been developed at Argonne for estimating the performance and manufacturing cost of lithium-ion batteries for electric-drive vehicles including hybrid-electrics (HEV), plug-in hybrids (PHEV) and pure electrics. The effort is being funded by the Vehicle Technology Office (VTO), which is part of the Energy Efficiency and Renewable Energy (EERE) office of the U.S. Department of Energy (USDOE). BatPaC was first developed in 2007, was subsequently peer-reviewed, and it has served Argonne researchers and the greater battery community in studying the impact of material properties on performance at the pack level. With further developments, the model now allows the design of cells and battery packs for automotive applications, to meet performance requirements (power, energy, recharge time), and estimates the cost of manufacturing the designed batteries. Since the cost depends on the materials he design, and the manufacturing process, this bottom-up model/tool enables the user to study their effects. Designed or the lithium-ion cell and battery researcher, BatPaC helps answer many questions by being 1.Transparent in the assumptions made and the method of calculation 2.Capable of designing a battery specifically for the requirements of an application 3.Constrained by the physical limitations that govern battery performance 4. A bottom-up calculation approach to account for every cost factors. BatPaC predicts the impact of promising materials (and their properties) on the performance metrics relevant for the different applications. Researchers can use the specific capacities and the half-cell voltages of a particular set of electrode materials to calculate the mass and volume of a cell to develop a model incorporating the properties of all the other materials in the cell and the design of the cell enclosure. These calculations not only reveal the impact of an improved material but also enable researchers to calculate the material properties that would be needed to meet the performance criteria of a full battery pack. With this information, researchers can provide he battery industry with realistic expectations that will help it more successfully advance novel battery technologies an Applications.

AHMED, SHABBIR↗
DOE OSTI2022

Quantum Orbital Minimization Method for Excited States Calculation on a Quantum Computer

Herein we propose a quantum-classical hybrid variational algorithm, the quantum orbital minimization method (qOMM), for obtaining the ground state and low-lying excited states of a Hermitian operator. Given parametrized ansatz circuits representing eigenstates, qOMM implements quantum circuits to represent the objective function in the orbital minimization method and adopts a classical optimizer to minimize the objective function with respect to the parameters in ansatz circuits. The objective function has an orthogonality constraint implicitly embedded, which allows qOMM to apply a different ansatz circuit to each input reference state. We carry out numerical simulations that seek to find excited states of H 2 , LiH, and a toy model consisting of four hydrogen atoms arranged in a square lattice in the STO-3G basis with UCCSD ansatz circuits. Comparing the numerical results with existing excited states methods, qOMM is less prone to getting stuck in local minima and can achieve convergence with more shallow ansatz circuits.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Evaluating fast methods for static polarizabilities on extended conjugated oligomers

Given the importance of accurate polarizability calculations to many chemical applications, coupled with the need for efficiency when calculating the properties of sets of molecules or large oligomers, we present a benchmark study examining possible calculation methods for polarizable materials. We first investigate the accuracy of the additive model used in GFN2, a highly-efficient semi-empirical tight-binding method, and the D4 dispersion model, comparing its predicted additive polarizabilities to ωB97XD results for a subset of PubChemQC and a compiled benchmark set of molecules spanning polarizabilities from approximately 3 Å 3 to 600 Å 3 , with some compounds in the range of approximately 1200–1400 Å 3 . Although we find additive GFN2 polarizabilities, and thus D4, to have large errors with polarizability calculations on large conjugated oligomers, it would appear an empirical quadratic correction can largely remedy this. We also compare the accuracy of DFT polarizability calculations run using basis sets of varying size and level of augmentation, determining that a non-augmented basis set may be used for large, highly polarizable species in conjunction with a linear correction factor to achieve accuracy extremely close to that of aug-cc-pVTZ.

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

Neutron diffusion calculation in heterogeneous geometry based on local/global iteration using proper orthogonal decomposition

This study newly proposes a heterogeneous core calculation method based on local/global iteration using proper orthogonal decomposition (POD). By using the singular value decomposition (SVD) and the low-rank approximation, appropriate POD bases for expanding the neutron flux can be obtained from snapshot data of the neutron flux obtained by fine mesh calculations. By projection using the POD bases, the dimension of the target equation (e.g., discretized neutron diffusion equation) can be dramatically reduced. In the proposed method, POD is effectively applied to each single assembly calculation (local calculation). Furthermore, using the local/global iteration, the effective neutron multiplication factor and the neutron flux distribution in the whole core geometry can be obtained by combining the numerical results of the local calculation for each fuel assembly and the global calculation for the whole core. As a feasibility study, the proposed method is applied to a one-dimensional heterogeneous core analysis, and the accuracy is investigated by changing the total number of POD bases. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Challenges for density functional theory in simulating metal–metal singlet bonding: A case study of dimerized VO2

VO2 is renowned for its electric transition from an insulating monoclinic (M1) phase, characterized by V–V dimerized structures, to a metallic rutile (R) phase above 340 K. This transition is accompanied by a magnetic change: the M1 phase exhibits a non-magnetic spin-singlet state, while the R phase exhibits a state with local magnetic moments. Simultaneous simulation of the structural, electric, and magnetic properties of this compound is of fundamental importance, but the M1 phase alone has posed a significant challenge to the density functional theory (DFT). In this study, we show none of the commonly used DFT functionals, including those combined with on-site Hubbard U to treat 3d electrons better, can accurately predict the V–V dimer length. The spin-restricted method tends to overestimate the strength of the V–V bonds, resulting in a small V–V bond length. Conversely, the spin-symmetry-breaking method exhibits the opposite trends. Each of these two bond-calculation methods underscores one of the two contentious mechanisms, i.e., Peierls lattice distortion or Mott localization due to electron–electron repulsion, involved in the metal–insulator transition in VO2. To elucidate the challenges encountered in DFT, we also employ an effective Hamiltonian that integrates one-dimensional magnetic sites, thereby revealing the inherent difficulties linked with the DFT computations.

Chemistry↗
DOE OSTI2022

Analysis of SCALE Criticality and Sensitivity Calculations for Reflected HEU Cylinders [Abstract]

The SCALE code package offers several nuclear data libraries to support Monte Carlo (MC) transport, as well as MC-based derivation of $\kappa$ eff sensitivity and uncertainty (S/U) data. The CSAS sequence using the KENO MC code can utilize continuous-energy (CE) cross sections, or pre-generated multigroup (MG) cross section libraries. The use of MG libraries introduces bias into calculations in exchange for faster transport solutions. The TSUNAMI-3D sequence also utilizes KENO MC calculations. TSUNAMI-3D has two CE calculational methods: the Iterated Fission Probability (IFP) method, and the Contribution-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method. Previous work has shown poor agreement between CLUTCH and confirmatory direct perturbation calculations in specific applications, e.g., fissionable and polyethylene reflectors

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Deep Learning Methods for Symbolic Calculations in HEP

This project develops machine learning methods to accelerate symbolic calculations in high-energy physics. Using sequence-to-sequence transformer models, we construct frameworks to predict squared amplitudes and related quantities for Standard Model processes, including quantum electrodynamics, quantum chromodynamics, and electroweak interactions. The results demonstrate that deep learning can successfully learn complex symbolic relationships and provide a scalable approach to symbolic computation with potential applications in precision calculations and collider phenomenology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Hydrogen Risk Assessment Models (HyRAM) (Version 3.0 Technical Reference Manual)

The HyRAM software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen infrastructure and transportation systems. HyRAM is designed to facilitate the use of state-of-the-art science and engineering models to conduct robust, repeatable assessments of hydrogen safety, hazards, and risk. HyRAM includes generic probabilities for hydrogen equipment failures, probabilistic models for the impact of heat flux on humans and structures, and computationally and experimentally validated first-order models of hydrogen release and flame physics. HyRAM integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing hydrogen hazards (thermal effects from jet fires, overpressure effects from deflagrations), and assessing impact on people and structures. HyRAM is developed at Sandia National Laboratories for the U.S. Department of Energy to increase access to technical data about hydrogen safety and to enable the use of that data to support development and revision of national and international codes and standards. HyRAM is a research software in active development and thus the models and data may change. This report will be updated at appropriate developmental intervals. This document provides a description of the methodology and models contained in the HyRAM version 3.0. HyRAM 3.0 includes the new ability to model cryogenic hydrogen releases from liquid hydrogen systems, using a different property calculation method and different equations of state. Other changes include modifications to the ignition probability calculations, component leak frequency calculations, and addition of default impulse data.

08 HYDROGEN↗
DOE OSTI2026

Effects of screening and pressure ionization on the electron broadening of spectral lines in dense plasmas

Collisions between electrons and radiating atoms broaden spectral absorption and emission lines in dense plasmas. High densities also introduce screening and pressure ionization effects that distort the wave functions of both bound and free electrons. In order to study how dense plasma effects influence the electron broadening of spectral lines, this paper incorporates electron wave functions from an average-atom (AA) model to calculate the linewidth of the B III 2⁢𝑝−2⁢𝑠 transition at 𝑇 = 10 eV for mass densities ranging from 𝜌 = 10 −4 to 0.4⁢ g⁡/cc. The calculation method uses the impact approximation, allowing the linewidth to be written in terms of electron-collision cross sections and an interference term. Compared to an otherwise identical calculation that uses Coulomb free wave functions, the AA method is found to modify both the cross sections and the resulting linewidth at sufficiently high density by introducing screening and pressure-ionized bound states. Screening lowers the cross sections at low energies and near electron excitation thresholds, while pressure-ionized bound states introduce resonances into the continuum. Thus, as the density increases, the relative linewidth between the AA and Coulomb calculations follows a general decrease because of screening, with sharp increases at various intervals due to pressure ionization. Finally, the AA results are also compared with a common approach to introduce screening through the interaction potential and reduced models that use the Bethe formula for the inelastic electron-collision cross sections.

electronic excitation & ionization↗
DOE OSTI2022

twoaxistracking – a python package for simulating self-shading of two-axis tracking solar collectors

Self-shading in fields of two-axis tracking collectors typically ranges from 1% to 6% of the annual incident irradiation. It is thus essential to account for shading in order to obtain accurate yield estimates and financing for such solar projects. The present study presents the free and open-source Python package twoaxistracking for simulating self-shading in fields of two-axis tracking collectors. The package is freely available at: https://github.com/pvlib/twoaxistracking. The main steps of the method and mathematical formulation are described. Additionally, a demonstration of how to use the package is presented. The shading calculation method excels over previous methods found in the literature in that it can: handle arbitrary aperture geometries and distinguish between the total and active areas; account for sloped ground and collectors with different heights within the same field; reduce computation time by skipping calculations at high solar elevation angles.

14 SOLAR ENERGY↗
DOE OSTI2021

The dependence of the gradients of oxygen and nitrogen-to-oxygen on stellar age in MaNGA galaxies

For this work, we derived the oxygen abundance (O/H), the nitrogen-to-oxygen (N/O) abundance ratio, and their corresponding radial gradients for a sample of 1431 galaxies from the MaNGA DR15 survey using two different realisations of the strong line method: empirical R calibration and the Bayesian model-based HII-CHI-MISTRY (HCM) code. We find that both abundance calculation methods reveal a correlation between the O/H gradient and the stellar mass of a galaxy. This relation is non-linear, with the steepest average gradients in the intermediate mass range and flatter average gradients for high- and low-mass galaxies. The relation between the N/O gradient and the stellar mass is, on average, non-linear with the steepest gradients in the intermediate mass range (log(M/M ⊙ )~10), flatter gradients for high-mass galaxies, and the flattest gradients for low-mass galaxies. However, the general trend of steepening N/O gradients for higher masses, as reported in previous studies, remains evident. We find a dependence between the O/H and N/O gradients and the galaxy mean stellar age traced by the D(4000) index. For galaxies of lower masses, both gradients are, generally, steeper for intermediate values of D(4000) and flatter for low and high values of D(4000). Only the most massive galaxies do not show this correlation. We interpret this behaviour as an evolution of the metallicity gradients with the age of stellar population. Though the galaxies with a positive slope of the D(4000) radial gradient tend to have flatter O/H and N/O gradients, as compared to those with a negative D(4000) gradient.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Data-Driven Method for Adaptive Reserve Requirement Estimation via Probabilistic Net Load Forecasting

The University of Colorado Boulder Boulder, CO, 80309, USA With the increasing penetration of renewable energy, power systems are subject to more uncertainty. This makes power system reserve scheduling more challenging. Most of the current reserve requirement determination methods calculate reserve requirements based on historical data, which does not consider the real-time or future system uncertainty. In this paper, a data-driven method is developed to determine the non-spinning reserve requirement (NSRR) in the Electric Reliability Council of Texas (ERCOT) system. The method follows the procedure of the current ERCOT method while adaptively determining the NSRR based on probabilistic net load forecasts. Case studies with two years of ERCOT data show that the developed method significantly reduces the NSRR by introducing an adaptive temporal resolution and update rate. Sensitivity analysis with different forecasting and percentile thresholds indicates the flexibility of the developed method.

non-spinning reserve↗
DOE OSTI2020

Flexibility Estimation and Control of Thermostatically Controlled Loads with Lock Time for Regulation Service

A virtual battery model is a simple and general method to quantify aggregate flexibility from thermostatically controlled loads (TCLs), enabling grid operators to effectively coordinate a large number of flexible building loads with supplyside resources in power systems. Lockout controls are designed to avoid wear and tear resulting from short-cycling of hardware. The lock on/off time could significantly affect aggregate flexibility from TCLs to provide ancillary services and may even fail control algorithms designed without considering the lock time constraints. This paper focuses on flexibility estimation and control design for TCLs with lock time constraints to provide frequency regulation service. We first investigate the potential impacts of lock time on TCLs’ aggregate flexibility and control performance. Both control-dependent and control-independent power bounds are derived, based on either previous TCL switching operations or regulation signals. While the control-dependent method provides aggregate flexibility for a given control method, the control-independent method calculates the theoretical maximum of power bounds. Two control algorithms are proposed to better distribute flexibility over time and thereby improve signal tracking performance. The proposed methods are illustrated and validated through simulations.

Wang, Peng↗
DOE OSTI2020

A Method for Extrapolating CMM Data

This report provides a method for evaluating coordinate measuring machine (CMM) data and extrapolating the data to estimate the dimensions of a cylindrical component beyond the sampled measurements. The theoretical basis is provided, and a calculational method is presented. The report also provides detailed, formatted, printed, and plotted representations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Intensity modulated Ir-192 brachytherapy using high-Z 3D printed applicators

Gynecologic cancers are often asymmetric, yet current Ir-192 brachytherapy techniques provide only limited radial modulation of the dose. The shielded solutions investigated here solve this by providing the ability to modulate between highly asymmetric and radially symmetric dose distributions at a given location. To find applicator designs that can modulate between full dose and less than 50% dose, at the dimensions of the urethra, a 2D calculation algorithm was developed to narrow down the search space. Two shielding design types were then further investigated using Monte Carlo and Boltzmann-solver dose calculation algorithms. 3D printing techniques using ISO 10993 certified biocompatible plastics and 3D printable tungsten-loaded plastics were tested. It was also found that shadowing effects set by the shape of the shielding cannot be easily modulated out, hence careful design is required. The shielded applicator designs investigated here, allow for reduction of the dose by over 50% at 5 mm from the applicator surface in desired regions, while also allowing radially symmetric dose with isodose line deviations less than 0.5 mm from circular. The shielding designs were also chosen with treatment delivery time in mind. Treatment times for these shielded designs were found to be less than 1.4 times longer than a 6-channel unshielded cylinder for the equivalent fully symmetric dose distribution. The 2D calculation methods developed here provide a simple way to rapidly evaluate shielding designs, while the 3D printing techniques also allow for devices with novel shapes to be rapidly prototyped. Both TOPAS Monte Carlo and Acuros BV calculations show that significant dose shaping and organ at risk sparing can be achieved without significantly compromising the plan in regions that require the full dose.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Quantifying the Temperature Dependence of the Multi-Species, Multi-Reaction Model: Part II. Estimation of Entropy Coefficient for Meso-Carbon Micro-Bead Graphite

In Part 1 of this paper, the temperature impact on the Multi-Species, Multi Reaction (MSMR) model was studied. This was accomplished by acquiring data from slow rate lithiation and delithiation of a meso-carbon micro-bead (MCMB) graphite. Through this analysis, the temperature impact on the total fraction of available host sites in a particular MSMR gallery (X j ), the impact on the reference potential (U$^o_j$), and the impact on the parameter detailing the deviation from Nernstian behavior (ω j ) was determined. Here, the intercalation material is discussed, compared to traditional methods of acquiring the entropy coefficient, and comparison is made to previous mathematical estimates. Some of the challenges in using temperature dependent constant rate charge and discharge data as compared to the potentiodynamic entropy coefficient calculation method are also discussed, and a recommendation for future applications is proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗