Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated data”

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 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Analysis of Interpretable Data Representations for 4D-STEM Using Unsupervised Learning

Abstract Understanding the structure of materials is crucial for engineering devices and materials with enhanced performance. Four-dimensional scanning transmission electron microscopy (4D-STEM) is capable of mapping nanometer-scale local crystallographic structure over micron-scale field of views. However, 4D-STEM datasets can contain tens of thousands of images from a wide variety of material structures, making it difficult to automate detection and classification of structures. Traditional automated analysis pipelines for 4D-STEM focus on supervised approaches, which require prior knowledge of the material structure and cannot describe anomalous or deviant structures. In this article, a pipeline for engineering 4D-STEM feature representations for unsupervised clustering using non-negative matrix factorization (NMF) is introduced. Each feature is evaluated using NMF and results are presented for both simulated and experimental data. It is shown that some data representations more reliably identify overlapping grains. Additionally, real space refinement is applied to identify spatially distinct sample regions, allowing for size and shape analysis to be performed. This work lays the foundation for improved analysis of nanoscale structural features in materials that deviate from expected crystallographic arrangement using 4D-STEM.

Bruefach, Alexandra (ORCID:0000000209323477)↗
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 OSTI2026

ACT University Collaboration Proposal: Role of manufacturing defects on material failure under dynamic loading for developing enhanced failure models and theories (Final Report)

Research on the tensile behavior of additively manufactured 316L stainless steel coupon specimens at increasing strain rates was conducted over the past year at Penn State. Dynamic loading rates into the 1000 1/s loading rates were performed on specimens with nearly 1.32” gage lengths. Stress-strain plots show ductile and plastic behavior well beyond 12% strain with necked specimens, having smaller effective gage lengths showing up to 56% ultimate strain. Additional tests performed on compact tension specimens helped with simulation work to understand the deformation behavior. Using finite elements, it was possible to determine the feasibility of a comprehensive experimental validation study towards an improved model for AM material failure – specifically, the Bai-Wierzbicki approach which accounts for lode angle and triaxiality. A non-significant number of tests are projected for eight specific failure nodes, with additional replicates to provide strain-rate capabilities to the existing formulation. One such correction factor is explored for quasi-static, notched specimens. Finally, stress intensity factor was explored using a set of compact-tension specimens which would also provide useful validation data for any simulation work.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spectral Effects in Albedo and Rearside Irradiance Measurement for Bifacial Performance Estimation

We investigate the impact of spectral dependence of ground surface reflectivity on albedo and rearside irradiance measurements necessary for bifacial photovoltaic (PV) module performance estimation and monitoring. Because PV modules are spectrally selective, albedo and irradiance measurements performed with common irradiance sensors may require spectral mismatch corrections when used for performance prediction. We investigate via simulation the differences in spectrally responsive albedo measured with thermopile pyranometers and crystalline silicon PV reference cells in comparison to a typical crystalline-silicon bifacial PV module. Simulations are performed for nine different representative ground surface materials using simulated solar spectra together with spectral reflectivity data distributed with the SMARTS simulation software. For the materials considered, the results show that albedo spectral mismatch relative to the bifacial module is distributed over a range of ±9.2% for thermopile pyranometers versus only ±3.7% for a typical PV reference cell. We consider the impact of this spectrally-responsive albedo mismatch on bifacial PV module rearside irradiance measurements. Using synthesized rearside spectral irradiance distributions, we find that for the nine different ground surface materials the predicted rearside irradiance measurement deviates from the effective irradiance observed by the PV module by on the order of 16.5 W/m 2 for the pyranometer and 3.6 W/m 2 for the PV reference cell. We discuss the implications for bifacial albedo and irradiance measurement.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2026

DESI DR2 Reference Mocks: Clustering results from UCHUU ELGs and QSOs

High-redshift galaxy clustering provides a powerful probe of the growth of structure, testing models of dark matter, dark energy, and galaxy formation during the epoch when the Universe was rapidly evolving. Emission line galaxies (ELGs) and quasars (QSOs) are used as tracers of dark matter by the Dark Energy Spectroscopic Instrument (DESI) to probe this redshift regime. We present results from ELG and QSO mock catalogs created from the Uchuu N-body simulation and tuned to DESI Data Release 2 (DR2) clustering. Employing a modified subhalo abundance matching (SHAM) technique, we populate Uchuu halos and subhalos with QSOs between 0.8 < z < 2.1. For ELGs, we modify this method to select satellite galaxies with low velocities relative to their associated central halos, and populate a separate set of Uchuu halos and subhalos with ELGs between 0.8 < z < 1.6. In this paper, we reproduce the redshift evolution of number density and clustering statistics across the fitted range of scales. We also measure the large-scale clustering bias of both the data and mock samples. These results improve simulated lightcone construction from cosmological models and enhance our understanding of the galaxy-halo connection.

Vaisakh, R. [Southern Methodist U.] (ORCID:0009000↗
DOE OSTI2020

Depletion capabilities in the OpenMC Monte Carlo particle transport code

A depletion solver has been implemented in OpenMC and is described herein. The depletion solver is implemented in Python and interfaces with OpenMC’s transport solver through a C++ application programming interface, which enables an in-memory transport-depletion coupling. Multiple integration methods for advancing in time have been implemented and exhibit tradeoffs in cost, accuracy, and memory use. For all time integration methods, evaluation of the matrix exponential is performed by using the incomplete partial fraction form of the Chebyshev rational approximation method. Simulations of a pressurized water reactor (PWR) pincell and a sodium-cooled fast reactor (SFR) assembly were carried out with OpenMC and Serpent. For both problems, the use of a high-fidelity depletion chain results in predictions of k eff that agree within 20–30 pcm between OpenMC and Serpent. Predicted actinide concentrations were found to agree to a fraction of a percent, and most fission product concentrations were found to agree within 1%. Here, the few cases where larger differences were observed can be attributed either to differences in how the energy dependence of fission product yields is handled or deficiencies in the nuclear data used. OpenMC simulations of the PWR and SFR problems using a simplified 228-nuclide depletion chain demonstrate that it achieves accuracy close to that of the full, high-fidelity depletion chain with respect to the studied responses.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Evaluation of a practical approach for field scale moisture flow modeling in heterogeneous media at a semiarid site

Abstract A practical approach for modeling field‐scale moisture flow in a highly heterogeneous unsaturated medium is described in this study. The validity of this approach is demonstrated through comparison of the numerical simulations with field observations at a semiarid site located in southcentral Washington State. The methodology is based on upscaling the core scale hydraulic properties and combining power‐law and tensorial connectivity‐tortuosity (PA‐TCT) approaches to derive macroscopic anisotropy parameters for each hydrostratigraphic unit (HSU) identified in the field. Each heterogeneous HSU is approximated by an equivalent homogeneous medium (EHM) model for which PA‐TCT parameters are used in the flow simulations. The available field data on moisture content and matric potential are compared with steady‐state flow simulations based on the mean form of Richards' equation. While the homogenization or averaging of heterogeneities, embedded in the EHM modeling approximation, cannot capture all of the field‐scale variability, the simulated steady‐state moisture and matric potential profiles capture well the central tendency of the field data. This approach is deemed practical for assessing the fate and transport of contaminants in highly heterogeneous unsaturated media at the transport scale of hundreds of meters.

Khaleel, Raziuddin↗
DOE OSTI2024

ZFP: A compressed array representation for numerical computations

HPC trends favor algorithms and implementations that reduce data motion relative to FLOPS. We investigate the use of lossy compressed data arrays in place of traditional IEEE floating point arrays to store the primary data of calculations. Simulation is fundamentally an exercise in controlled approximation, and error introduced by finite-precision arithmetic (or lossy compression) is just one of several sources of error that need to be managed to ensure sufficient accuracy in a computed result. We describe ZFP, a compressed numerical format designed for in-memory storage of multidimensional arrays, and summarize theoretical results that demonstrate that the error of repeated lossy compression can be bounded and controlled. Furthermore, we establish a relationship between grid resolution and compression-induced errors and show that, contrary to conventional floating point, ZFP reduces finite-difference errors with finer grids. We present example calculations that demonstrate data reduction by 4x or more with negligible impact on solution accuracy. Our results further demonstrate several orders-of-magnitude increase in accuracy using ZFP over IEEE floating point and Posits for the same storage budget.

Lindstrom, Peter↗
DOE OSTI2020

Effects of temperature and density evolution in MHD simulations of HIT-SI

The helicity injected torus-steady inductive (HIT-SI) experiment uses steady inductive helicity injection to form a spheromak equilibrium and sustain the structure against resistive decay. Helicity injection is performed using two half-tori “injectors” connected to the main plasma volume, whose fields are oscillated in an AC manner. The properties of the sustained spheromak equilibrium have been experimentally observed to vary with the frequency of the injector oscillation, producing higher current gains and more-symmetric and outwardly shifted current centroids with higher frequency. A computational scan of injector frequency using the 3D MHD code PSI-Tet, which models the entire HIT-SI plasma volume including the injectors, has been performed, including a comparison of the results using the full Hall MHD model to those obtained using a simplified “zero-beta” (constant temperature and density) model. The results of both PSI-Tet models are also compared with experimental data and with simulations using the NIMROD code, which does not model the injector regions. The results of the PSI-Tet simulations show that the average temperature and current gain increase with injector frequency, in agreement with experimental trends. The simulations also show qualitative changes in the dynamics of several quantities with increasing injector frequency, such as density oscillations and current evolution. However, the outward shift and symmetrizing of the current centroid, observed experimentally, are not observed in these MHD simulations, indicating that unresolved or excluded dynamics may be important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage (CRADA Report)

NREL and Element 16 collaborated on sulfur thermal energy storage modeling using NREL’s high performance computing (HPC) resources to assist its application in industrial processes. Industrial process heat (IPH) accounts for ~70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately, 1500 TWht (~60% Terawatt hour thermal) of IPH demand is in the temperature range of 100-300°C. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. This HPC4EI project facilitated Element 16’s development of low-cost and high-impact molten sulfur TES for dispatchable IPH. The development of a high-fidelity model validated by experimental data and HPC simulations enabled the successful resolution of the complex interplay between fluid dynamics and heat transfer processes during transient operation of sulfur TES, overcoming the numerical challenges posed by the non-linear temperature-dependent physical properties of sulfur. The project helped accelerate Element 16’s molten sulfur TES product design and support its broad applications.

25 ENERGY STORAGE↗
DOE OSTI2024

Review and meta-analysis of electron temperatures from high-intensity laser–solid interactions

The accelerated electron spectrum from high-intensity laser–solid interaction is often conveniently described using a Boltzmann distribution, whose temperature is known within the field as the hot-electron temperature. The importance of the electron temperature is highlighted by the sheer number of experimental and simulation studies on the subject over the past three decades. Recently, multi-kJ, multi-ps pulses have yielded electron spectra with temperatures far beyond the expected ponderomotive result. Expressions that predict the electron temperature considering laser parameters beyond intensity and wavelength have been developed, albeit using small datasets. In this review, we present what is, to the best of our knowledge, the largest dataset of electron temperatures gathered from experimental measurements and particle-in-cell simulations. This dataset allows us to compare existing analytical and empirical hot-electron temperature scaling models over a wide parameter range. We also develop new scaling models that incorporate the laser pulse duration of the laser and the plasma scale length. Three models that include pulse-duration and scale length dependence are especially successful at predicting both simulated and experimental data. The dataset will soon be made publicly available to encourage further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Uncertainty Quantification of Global Net Methane Emissions From Terrestrial Ecosystems Using a Mechanistically Based Biogeochemistry Model

Quantification of methane (CH4) emissions from wetlands and its sinks from uplands is still fraught with large uncertainties. Here, a methane biogeochemistry model was revised, parameterized, and verified for various wetland ecosystems across the globe. The model was then extrapolated to the global scale to quantify the uncertainty induced from four different types of uncertainty sources including parameterization, wetland type distribution, wetland area distribution, and meteorological input. We found that global wetland emissions are 212 ± 62 and 212 ± 32 Tg CH4 year -1 (1Tg = 10 12 g) due to uncertain parameters and wetland type distribution, respectively, during 2000–2012. Using two wetland distribution data sets and three sets of climate data, the model simulations indicated that the global wetland emissions range from 186 to 212 CH 4 year -1 for the same period. The parameters were the most significant uncertainty source. After combining the global methane consumption in the range of -34 to -46 Tg CH 4 year -1 , we estimated that the global net land methane emissions are 149–176 Tg CH 4 year -1 due to uncertain wetland distribution and meteorological input. Spatially, the northeast United States and Amazon were two hotspots of methane emission, while consumption hotspots were in the Eastern United States and eastern China. During 1950–2016, both wetland emissions and upland consumption increased during El Niño events and decreased during La Niña events. This study highlights the need for more in situ methane flux data, more accurate wetland type, and area distribution information to better constrain the model uncertainty.

58 GEOSCIENCES↗
DOE OSTI2021

Downscaling SSP-consistent global spatial urban land projections from 1/8-degree to 1-km resolution 2000–2100

Abstract Long-term, spatial urban land projections that simultaneously offer global coverage and local-scale empirical accuracy are rare. Recently a set of such projections was produced using data-science-based simulations and the Shared Socioeconomic Pathways (SSPs). These projections update at decadal time intervals from 2000 to 2100 with a spatial resolution of 1/8 degree, while many socio-environmental studies customarily run their analysis and modelling at finer spatial resolutions, e.g. 1-km. Here we develop and validate an algorithm to downscale the 1/8-degree spatial urban land projections to the 1-km resolution. The algorithm uses an iterative process to allocate the decadal amount of urban land expansion originally projected for each 1/8-degree grid to its constituent 1-km grids. The results are a set of global maps showing urban land fractions at the 1-km resolution, updated at decadal intervals from 2000 to 2100, under five different urban land expansion scenarios consistent with the SSPs. The data can support studies of potential interactions between future urbanization and environmental changes across spatial and temporal scales.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Validation and Independent Uncertainty Analysis of the MIX-SOL-THERM-003 ICSBEP Benchmark

The International Criticality Safety Benchmark Evaluation Project (ICSBEP) was started in 1992 by the United States Department of Energy and later in 1995 became an international project with contributions from 22 countries. The project is now organized by the OECD (Organisation for Economic Co-operation and Development) Nuclear Energy Agency (NEA). In its most recent iteration, the ICSBEP handbook contains over five thousand evaluations of critical, near-critical, and subcritical experiments conducted in facilities all around the world. These benchmarks serve as valuable information for criticality safety engineers who can use them to validate calculation techniques and establish minimum subcritical margins for operations with fissionable materials. The benchmarks in the handbook are categorized by their fissile material composition, material form (oxide, solution, or metal), and fission energy spectra. This is especially useful for those looking for benchmarks similar to a system they are working on to compare methods and identify trends. The ICSBEP Handbook Uncertainty Guide is document outlining recommended practices and methods for determining uncertainties in these benchmarks. Quantifying these uncertainties thoroughly is crucial as it allows a higher degree of confidence that data used from them is valid and relevant. The guide stresses the importance of a thorough and well documented uncertainty analysis when evaluating an experiment. All measured values of a system, whether they be dimensions or material compositions, have a certain amount of uncertainty associated with them and can be analyzed one by one to determine their effects on the system. Many evaluated benchmarks in the handbook present this in detail, however some do not, mostly earlier evaluations performed in the 1990’s and early 2000’s. Recently at Los Alamos National Laboratory (LANL), the Nuclear Criticality Safety Division (NCSD) of LANL has been validating MCNP6.2 ® input files of criticality benchmarks for use by Whisper, a criticality safety code developed at LANL. This effort is also part of the OECD NEA Working Party on International Nuclear Data Evaluation Co-operation (WPEC) Subgroup 45, also known as Validation of Nuclear Data Libraries (VaNDaL). The goal of VaNDaL is to compile a set of validated simulation inputs for use in validating nuclear data and simulation codes. As part of these efforts, one of the benchmarks reviewed was the MIX SOL-THERM-003 ICSBEP benchmark. This paper provides an independent uncertainty analysis of this benchmark experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Sparse-Data Deep Learning Strategies for Radiographic Non-Destructive Testing

Radiography is an imaging technique used in a variety of applications, such as medical diagnosis, airport security, and nondestructive testing. We present a deep learning system for extracting information from radiographic images. We perform various prediction tasks using our system, including material classification and regression on the dimensions of a given object that is being radiographed. Our system is designed to address the sparse-data issue for radiographic nondestructive testing applications. It uses a radiographic simulation tool for synthetic data augmentation, and it uses transfer learning with a pre-trained convolutional neural network model. Using this system, our preliminary results indicate that the object geometry regression task saw an improvement of 70% in the R-squared value when using a multi-regime model. In addition, we increase the performance of the object material classification tasks by utilizing data from different imaging systems. In particular, using neutron imaging improved the material classification accuracy by 20% when compared to x-ray imaging.

convolutional neural networks↗
DOE OSTI2022

Directed gas phase preparation of ethynylallene (H 2 CCCHCCH; X 1 A') via the crossed molecular beam reaction of the methylidyne radical (CH; X 2 Π) with vinylacetylene (H 2 CCHCCH; X 1 A')

The gas-phase bimolecular reaction of the methylidyne (CH; X 2 Π) radical with vinylacetylene (H 2 CCHCCH; X 1 A') was conducted at a collision energy of 20.3 kJ mol –1 under single collision conditions exploiting the crossed molecular beam experimental results merged with ab initio electronic structure calculations and ab initio molecular dynamics (AIMD) simulations. The laboratory data reveal that the bimolecular reaction proceeds barrierlessly via indirect scattering dynamics through long-lived C 5 H 5 reaction intermediate(s) ultimately dissociating to C 5 H 4 isomers along with atomic hydrogen with the latter predominantly originating from the vinylacetylene reactant as confirmed by the isotopic substitution experiments in the D1-methylidyne–vinylacetylene reaction. Combined with ab initio calculations of the potential energy surface (PES) and statistical Rice–Ramsperger–Kassel–Marcus (RRKM) calculations, the experimental determined reaction energy of –146 ± 26 kJ mol –1 along with the distribution minimum of T(θ) at 90° and isotopic substitution experiments suggest ethynylallene (p1; ΔrG = –230 ± 4 kJ mol –1 ) as the dominant product. The ethynylallene (p1) may be formed with extensive rovibrational excitation, which would result in a lower maximum translational energy. Further, AIMD simulations reveal that the reaction dynamics leads to p1 (ethynylallene, 75%) plus atomic hydrogen with the dominant initial complex being i1 formed by methylidyne radical addition to the double C[double bond, length as m-dash]C bond in vinylacetylene. Altogether, combining the crossed molecular beam experimental results with ab initio electronic structure calculations and ab initio molecular dynamics (AIMD) simulations, ethynylallene (p1) is expected to represent the dominant product in the reaction of the methylidyne (CH; X 2 Π) radical with vinylacetylene (H 2 CCHCCH; X 1 A').

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗