Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALCULATION METHODS”

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

The Koopman Operator: Capabilities and Recent Advances

In this paper, we provide an introduction to the Koopman Operator (KO) designed to be accessible to those not already familiar with the field. Our aim is to expose domain experts and controls practitioners to the concept and the capabilities the KO provides in the interest of promoting wider KO use; as such, we deliberately focus more on the uses and computational aspects of the KO than on the wealth of analytical results in the literature. Those researchers may then be able to pose new applications-based research questions for the KO community. We begin by defining the KO and describing its key components. Secondly, we discuss computational methods used to calculate the KO in practice. Thirdly, we provide an overview of different uses for the KO and its associated capabilities. Finally, we discuss some open areas of research for future KO work.

koopman operator, dynamical systems, optimal contr↗
DOE OSTI2025

FuelLib (Fuel Library) [SWR-25-26]

FuelLib is a library that utilizes the group contribution method (GCM) for calculating thermodynamic properties of hydro-carbon jet fuels. FuelLib utilizes the tables and functions of the GCM as proposed by Constantinou and Gani (1994) and Constantinou, Gani and O'Connel (1995), with additional physical properties discussed in Govindaraju & Ihme (2016). The code is based on Pavan B. Govindaraju's Matlab implementation of the GCM, and has been expanded to include additional thermodynamic properties and mixture properties. The fuel library contains gas chromatography (GC x GC) data for a variety of fuels ranging from simple single component fuels to complex jet fuels. The GC x GC data for POSF jet fuels comes from Edwards (2020).

Montgomery, David [National Renewable Energy Labor↗
DOE OSTI2026

HygroThermFEM v1.0

HygroThermFEM is a Finite Element Method-based numerical calculation engine for solving 2-D heat and moisture transfer problems. This numerical engine is used in the THERM software tool, and its primary purpose is for the analysis of building envelopes (e.g., windows, walls, roofs, foundations, etc.). However, the engine can also be used for any heat and moisture transfer problems that require solving fundamental 2-D energy and mass transfer equations. Fluid flow solutions (Navier-Stokes momentum equations) are not included, but the correlations for various convection heat transfer situations are provided, including the translation of complex cavity geometries into those for which correlations are applicable. The calculation engine is written in C++ and includes an API for connecting to third-party tools.

Vidanovic, Dragan [Lawrence Berkeley National Labo↗
DOE OSTI2022

Neutronic Safety Analysis of Pu-238 Production at Idaho National Laboratory

This analysis was completed to support the irradiation of plutonium fuel services (PFS) targets in the NEFT in the Advanced Test Reactor (ATR) as a part of the campaign to restart domestic production of plutonium-238 used in radioisotope power systems (RPS) by the National Aeronautical and Space Administration (NASA) and Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Nuclear Infrastructure Program (NE-3). Referred to as the PFS-ATR-GEN1-NEFT experiment, the assembly was designed to hold 46 PFS targets in the NEFT. The scope of this paper is to outline the MOPY method used to calculate the heat generation rates (HGRs), flux, fission density, and quantify the viability of the target design for Pu-238 production in ATR.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

Hanford Site Composite Analysis Data Package: Exposure Scenarios and Radionuclide Specific Dose Conversion Factors

This data package summarizes the exposure assumptions, equations, and methods used to calculate radionuclide-specific unit dose factors and the radiological doses for both groundwater and atmospheric pathways as a part of the updated Hanford Site Composite Analysis (CA). An All Pathways Representative Person exposure scenario is considered to evaluate exposure via both groundwater and atmospheric transport pathways. The radiological dose assessments for both groundwater and atmospheric pathways are included in the performance assessments (PAs) for various Waste Management Areas at the Hanford Site. This data package calculates exposure route-specific and total unit dose factors for CA-specific radionuclides of concern based on the exposure assumptions used in the CA and PAs. This data package also presents the results and comparison of the radionuclide-specific unit dose factors based on the exposure assumptions used in the revised Hanford Site Composite Analysis and Hanford Site performance assessments for low-level waste disposal facilities.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Uncertainty Analysis and Software Verification

Uncertainty analyses are an important part of calibrations and testing. They allow researchers insight on how to reduce and mitigate errors in testing. Each device used in testing introduces error in a system, as well as other sources such as environmental conditions, electronics, analog to digital, and random errors. Each source is carefully examined to identify how much error it introduces to a system. These sources are then combined using various methods of uncertainty calculations. To verify and validate software, a manual calculation is required to ensure the software is performing as intended. Using Excel to verify the calculations, we can identify discrepancies within the software. The root of the sum of the squares uncertainty (RSS) is used to find the combined uncertainty of a device at one and two standard deviations of the mean. Calibrations on accelerometers are performed using a vibration system along with a back-toback reference accelerometer. The vibration system takes a reference point at 100Hz frequency at 10g amplitude. The sensitivities are collected at each dialed in frequency. The sensitivity of the device represents the electrical output of the UUT (mV, pC, etc.) per unit of acceleration (g). The full history of sensitivities of selected accelerometers are used to find the averaged, standard deviation, and uncertainty of the device at each frequency tested. The uncertainty calculations from excel and the software are then compared.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Engineering Assessment of UO 2 and Cladding Behavior under High Burnup LOCA Conditions

To maximize the data extracted from a limited number of high-burnup fuel rod samples, several modeling efforts were performed to elucidate the fuel and the cladding responses of these fuels under transient conditions. These objectives were (1) to determine the role of the fuel stress state in fuel pulverization, (2) to ascertain the differences between conditions of cladding burst during experiments and those expected during a commercial reactor transient, and (3) to develop a method to conservatively calculate the geometry of the cladding rupture’s opening to inform fuel dispersal susceptibility.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Hanford Site Composite Analysis Data Package: Exposure Scenarios and Radionuclide Specific Dose Conversion Factors.

This data package summarizes the exposure assumptions, equations, and methods used to calculate radionuclide-specific unit dose factors and the radiological doses for both groundwater and atmospheric pathways as a part of the revised Hanford Site Composite Analysis. An All-Pathways Representative Person exposure scenario is considered to evaluate exposure via both groundwater and atmospheric transport pathways. The radiological dose assessments for both groundwater and atmospheric pathways are included in the performance assessments for various Waste Management Areas at the Hanford Site. This data package calculates exposure route-specific and total unit dose factors for composite-analysis-specific radionuclides of concern based on the exposure assumptions used in the composite analysis and performance assessments. This data package presents the results and comparison of the radionuclide-specific unit dose factors based on the exposure assumptions used in the revised composite analysis and various performance assessments.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

The Analysis and Support for LANL Technology Testbed: High Altitude Balloon Gondola and Orion Eagle Gamma Ray Detection Instrument

The Analysis and Support for LANL Technology Testbed: High Altitude Balloon Gondola and Orion Eagle Gamma-Ray Detection Instrument is a project where Los Alamos National Laboratory (LANL), in conjunction with the University of Michigan, are looking into gamma-ray detection in a space-like environment. The purpose of the detection is to work on instruments that can assist with planetary exploration and nuclear deterrents using a High-Altitude Balloon (HAB). For this paper, the main focus is on the mechanical aspects related to the structural analysis of the gondola. It is vital to understand the structural analysis because it helps determine if the gondola can survive the changes in the environment and the different forces or shocks that it will experience when the parachute is deployed and at landing. Additionally, NASA's Balloon Program Office (BPO) has specific requirements to meet to achieve flight certifications. Part of those certifications requires the gondola to simulate different scenarios at specific conditions with various forces applied and orientations. Therefore, it is necessary to adjust the calculations of each method and the force applied to get the right results. Once that is done, the information is inputted into the structural analysis modeling system. The system will aid the mechanical team in assessing the conditions of the gondola at each scenario and determine its margin of safety. After the gondola and instrument are certified, the payload will be able to fly. The HAB and payload are launched from NASA's Columbia Scientific Balloon Facility (CSBF) at Fort Sumner, NM. The official flight was on September 26th, 2021.

42 ENGINEERING↗
DOE OSTI2022

Geometry, Disorder and Phase Transitions in Topological States of Matter

The quantum Hall effect is the birthplace of topological states of matter, a major theme at the forefront of condensed matter physics in the past two decades. The fractional quantum Hall (FQH) effect revolutionized our understanding of phases of electronic matter. FQH states support exotic fractionally charged excitations that obey Abelian or non-Abelian fractional statistics, which are topological excitations that result from the underlying topological order. During this project, our group discovered a previously unrecognized geometric degree of freedom of incompressible FQH states and studied that for a variety of gapped FQH states. We brought this new concept into direct contact with experiments for the first time by generalizing it to Fermi-liquid states of composite fermions. Using the newly formulated powerful infinite Density Matrix Renormalization Group method, our numerical calculations yielded a parameter free prediction that was found to be in excellent agreement with experimental findings on electron systems in semiconductor heterostructures. In parallel, we performed extensive numerical studies on different, competing phases at various Landau level filling factors, and quantum phase transitions that result from such a competition, e.g. Abelian-non-Abelian phase transitions in bilayer systems. We studied geometrical excitations dubbed “gravitons” (because of their analogy with excitations in the theory of gravitation) and ways to excite and detect them, and explored how they couple with topological excitations. In graphene-based chiral materials, we realized the ability to tune through different incompressible and compressible states in a single Landau level, and found appropriate experimental parameters for the exploration of universal Luttinger liquid behavior not obtained in semiconductor-based electron systems. We showed that topological systems had a very different response from nontopological systems to strong disorder (many-body localization) as well as periodic drives.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Progress Report of Alloy 709 Performance in Liquid Sodium

With improved high-temperature properties, Alloy 709 shows great potential for advanced reactor applications. To understand the effects of sodium on A709, small tensile specimens were machined from two commercial heats and subject to sodium exposure in two forced convection loops operated by Argonne. Thermal aging experiments were also conducted in parallel to isolate any sodium-specific effects. An increasing trend in yield strength and a decrease trend in ultimate tensile strength can be seen with the increasing exposure time, suggesting microstructural changes at these temperatures. Nonetheless, the thermally aged and sodium exposed samples behaved similarly, implying no sodium-specific effect at these temperatures. In contrast, 316H SS exposed to sodium at 650°C showed a strong sodium effect. The different responses of A709 and 316H SS to sodium exposure is of interest and requires additional effort to understand. The cold traps of the sodium loop were replaced early in FY23. During a maintenance operation, an unexpected sodium leak incident occurred at one of the sodium loops, triggering an emergency response. An investigation into the incident revealed that the leak happened at a VCR fitting on the recently replaced cold trap. The metal gasket of the leaked VCR fitting was not compressed sufficiently. A dimensional mismatch between the new component and old piping configuration was responsible for a significant misalignment. A new cold trap with improved design is in production. The new design provides additional flexibility and an access point for post-installation helium leak check and for venting of trapped gas. These design features will make the system more error-tolerant and the installation quality verifiable. A thermodynamic analysis was performed to understand the carburization-decarburization behavior of A709 in sodium environment. The equilibrium carbon activity in the alloy was evaluated with M 23 C 6 , TiC, and NbC. The carburization-decarburization boundary was calculated with different methods developed previously. The results suggests that A709 will be carburized in SMT-1 and SMT-2, and can be decarburized in low-carbon-activity environments such as EBR-II at temperatures above 620-650°C. Nonetheless, A709 should be less likely to be decarburized than 316H SS when exposed to the same low-carbon-activity environment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid– solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Toward a microscopic picture of hadronization and multi-parton processes

This project advanced the understanding of how quarks and gluons produced in high-energy collisions transform into the hadrons observed in particle detectors, a fundamental process known as quantum chromodynamics (QCD) hadronization. By combining theoretical calculations, quantum simulation methods, and modern AI techniques, the research developed new tools to study multi-parton dynamics and nonperturbative effects that are essential for interpreting data from current and future nuclear physics experiments. Key outcomes include new theoretical frameworks for jet and hadron measurements, pioneering quantum simulation algorithms for real-time dynamics in field theories, and the development of advanced machine-learning models, such as diffusion models and explainable classifiers, to simulate and analyze collider events. These results are directly relevant to experiments at Jefferson Lab, Brookhaven National Laboratory, and the future Electron-Ion Collider, and they also have a broader impact in areas such as quantum information science and data-driven modeling of complex systems. The project supported the training of graduate students and postdoctoral fellows and contributed to the broader scientific community through publications, workshops, and collaborative activities. Overall, this work provides new insights into the microscopic mechanisms of hadron formation and establishes a foundation for future studies at the intersection of nuclear physics, artificial intelligence, and quantum computing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Efficient Ensemble-Based Stochastic Gradient Methods for Optimization Under Geological Uncertainty

Ensemble-based stochastic gradient methods, such as the ensemble optimization (EnOpt) method, the simplex gradient (SG) method, and the stochastic simplex approximate gradient (StoSAG) method, approximate the gradient of an objective function using an ensemble of perturbed control vectors. These methods are increasingly used in solving reservoir optimization problems because they are not only easy to parallelize and couple with any simulator but also computationally more efficient than the conventional finite-difference method for gradient calculations. In this work, we show that EnOpt may fail to achieve sufficient improvement of the objective function when the differences between the objective function values of perturbed control variables and their ensemble mean are large. On the basis of the comparison of EnOpt and SG, we propose a hybrid gradient of EnOpt and SG to save on the computational cost of SG. We also suggest practical ways to reduce the computational cost of EnOpt and StoSAG by approximating the objective function values of unperturbed control variables using the values of perturbed ones. We first demonstrate the performance of our improved ensemble schemes using a benchmark problem. Results show that the proposed gradients saved about 30–50% of the computational cost of the same optimization by using EnOpt, SG, and StoSAG. As a real application, we consider pressure management in carbon storage reservoirs, for which brine extraction wells need to be optimally placed to reduce reservoir pressure buildup while maximizing the net present value. Results show that our improved schemes reduce the computational cost significantly.

58 GEOSCIENCES↗
DOE OSTI2022

Solar data uncertainty impacts on MCMC methods for r-process nucleosynthesis

In recent work, we developed a Markov Chain Monte Carlo (MCMC) procedure to predict the ground state masses capable of forming the observed Solar r-process rare-earth abundance peak. By applying this method to nucleosynthesis calculations which make use of distinct astrophysical conditions and comparing our results to the latest precision mass measurements, we are able to shed light on the conditions/masses capable of producing a rare-earth peak which matches Solar data. Here we examine how our mass predictions change when using a few different sets of r-process Solar abundance residuals that have been reported in the literature. We explore how the differing error estimates of these Solar evaluations propagate through the Markov Chain Monte Carlo to our mass predictions. We find that Solar data which reports the rare-earth peak to have its highest abundance at mass number A = 162 can require distinctly different mass predictions from data with the peak centered at A = 164. Nevertheless, we find that two important general conclusions from past work, regarding the inconsistency of ‘cold’ astrophysical outflows with current mass measurements and the need for local stability at N = 104 in ‘hot’ scenarios, remain robust in the face of differing Solar data evaluations. Additionally, we show that the masses our procedure finds capable of producing a peak at A < 164 are not in line with the latest precision mass measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Ab Initio Phase Diagram of Chromium to 2.5 TPa

Chromium possesses remarkable physical properties such as hardness and corrosion resistance. Chromium is also a very important geophysical material as it is assumed that lighter Cr isotopes were dissolved in the Earth’s molten core during the planet’s formation, which makes Cr one of the main constituents of the Earth’s core. Unfortunately, Cr has remained one of the least studied 3d transition metals. In a very recent combined experimental and theoretical study (Anzellini et al., Scientific Reports, 2022), the equation of state and melting curve of chromium were studied to 150 GPa, and it was determined that the ambient body-centered cubic (bcc) phase of crystalline Cr remains stable in the whole pressure range considered. However, the importance of the knowledge of the physical properties of Cr, specifically its phase diagram, necessitates further study of Cr to higher pressure. In this work, using a suite of ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries, we obtain the theoretical phase diagram of Cr to 2.5 TPa. We calculate the melting curves of the two solid phases that are present on its phase diagram, namely, the lower-pressure bcc and the higher-pressure hexagonal close-packed (hcp) ones, and obtain the equation for the bcc-hcp solid–solid phase transition boundary. We also obtain the thermal equations of state of both bcc-Cr and hcp-Cr, which are in excellent agreement with both experimental data and QMD simulations. We argue that 2180 K as the value of the ambient melting point of Cr which is offered by several public web resources (“Wikipedia,” “WebElements,” “It’s Elemental,” etc.) is most likely incorrect and should be replaced with 2135 K, found in most experimental studies as well as in the present theoretical work.

equation of state↗
DOE OSTI2020

System and method for determining linear density of carbon fiber

A method for determining a linear density of a carbon fiber tow. The method includes providing a pulse of air directed toward the carbon fiber tow, determining, via a first sensor, an air pressure of the pulse of air, and measuring, via a second sensor, a displacement of the carbon fiber tow in response to the pulse of air directed toward the carbon fiber tow. The method further includes calculating, via a controller, a linear density of the carbon fiber tow based on the air pressure and the displacement, and outputting, via the controller, the linear density.

Koester, David Joseph↗
DOE OSTI2022

Customer outcomes in Pay-As-You-Save programs

We review the energy and financial outcomes of households participating in several programs based on successive versions of the Pay As You Save¯ (PAYS¯) system. PAYS¯ programs offer non-debt financing for energy efficiency (and sometimes other technologies) in residential buildings through a tariff attached to the home’s utility meter, designed to be offset by project savings. We find that the five programs we study generally serve customers living in zip codes with levels of income and education below the national average and unemployment rates above the national average, demonstrating their potential to improve equity in energy efficiency adoption. Using weather-normalized analysis of energy consumption data, we show that most customers of Midwest Energy’s program reduce annual electricity and gas consumption, averaging 15% and 26% reductions respectively. Changes in energy consumption calculated using this method represent a combination of project effects and changes in occupant behavior. These results are similar to existing analyses of PAYS¯ programs in North Carolina, Arkansas, and Tennessee. About half of participating Midwest households generate sufficient energy cost savings to cover their monthly tariff. Various factors, including changes in occupant behavior, program error, causes independent of the customer or program, or some combination thereof may explain lower-than-expected cost reductions in some projects. Given the inherent variability in annual household electricity consumption, we feel these programs are enabling energy efficiency improvements and their attendant co-benefits, including occupant health and comfort and reduced carbon emissions, while reasonably balancing energy savings and tariff costs. Pairing PAYS¯ with additional financial assistance, as well as promoting cost-effective measures such as air and duct sealing, could further broaden program participation by enabling additional projects to meet PAYS¯ program eligibility rules.

Deason, Jeff↗