SEARCH · Engineering Papers

Results for “factorization 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,459 records · Page 81

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Characterizing the magnetic fields of nearby molecular clouds using submillimeter polarization observations

Of all the factors that influence star formation, magnetic fields are perhaps the least well understood. The goal of this paper is to characterize the 3D magnetic field properties of nearby molecular clouds through various methods of statistically analysing maps of polarized dust emission. Our study focuses on nine clouds, with data taken from the Planck Sky Survey as well as data from the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry observations of Vela C. We compare the distributions of polarization fraction (p), dispersion in polarization angles ($\mathcal {S}$), and hydrogen column density (N H ) for each of our targeted clouds. To broaden the scope of our analysis, we compare the distributions of our clouds’ polarization observables with measurements from synthetic polarization maps generated from numerical simulations. We also use the distribution of polarization fraction measurements to estimate the inclination angle of each cloud’s cloud-scale magnetic field. We obtain a range of inclination angles associated with our clouds, varying from 16° to 69°. We establish inverse correlations between p and both $\mathcal {S}$ and NH in almost every cloud, but we are unable to establish a statistically robust $\mathcal {S}$ versus NH trend. By comparing the results of these different statistical analysis techniques, we are able to propose a more comprehensive view of each cloud’s 3D magnetic field properties. These detailed cloud analyses will be useful in the continued studies of cloud-scale magnetic fields and the ways in which they affect star formation within these molecular clouds.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics

Accurately modeling dense plasmas over wide-ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/$\mathrm{cm^3}$ to 16 g/$\mathrm{cm^3}$ and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature method. Here, we find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the sesame database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT-MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Data-Driven Performance Optimization of Gamma Spectrometers With Many Channels

In gamma spectrometers with variable spectroscopic performance across many channels (e.g., many pixels or voxels), a tradeoff exists between including data from successively worse-performing readout channels and increasing efficiency. Brute-force calculation of the optimal set of included channels is exponentially infeasible as the number of channels grows, and approximate methods are required. In this work, we present a data-driven framework for attempting to find near-optimal sets of included detector channels. The framework leverages non-negative matrix factorization (NMF) to learn the behavior of gamma spectra across the detector and clusters similarly-performing detector channels together. Performance comparisons are then made between spectra with channel clusters removed, which is more feasible than brute force. The framework is general and can be applied to arbitrary, user-defined performance metrics depending on the application. We apply this framework to optimizing gamma spectra measured by H3D M400 CdZnTe (CZT) spectrometers, which exhibit variable performance across their crystal volumes. In particular, we show several examples optimizing various performance metrics for uranium and plutonium gamma spectra in non-destructive assay (NDA) for nuclear safeguards, and explore trends in performance versus parameters such as clustering algorithm type. We also compare the NMF + clustering pipeline to several non-machine-learning (ML) algorithms, including several greedy algorithms. Although, we find that the NMF + clustering pipeline tends to find the best-performing set of detector voxels, significantly improving over the unoptimized spectra, but that a greedy accumulation of spectra segmented by detector depth can, in some cases, give similar performance improvements in much less computation time.

Energy resolution↗

DOE OSTI2021

Data-Driven Forgetting and Discount Factors for Vehicle Speed Forecasting in Ecological Adaptive Cruise Control

Abstract This paper investigates temporal correlations in human driving behavior using real-world driving to improve speed forecasting accuracy. These correlations can point to a measurement weighting function with two parameters: a forgetting factor for past speed measurements that the vehicle itself drove with, and a discount factor for the speeds of vehicles ahead based on information from vehicle-to-vehicle communication. The developed weighting approach is applied to a vehicle speed predictor using polynomial regression, a prediction method well-known in the literature. The performance of the developed approach is then assessed in both real-world and simulated traffic scenarios for accuracy and robustness. The new weighting method is applied to an ecological adaptive cruise control system, and its influence is analyzed on the prediction accuracy and the performance of the ecological adaptive cruise control in an electric vehicle powertrain model. The results show that the new prediction method improves energy saving from the eco-driving by up to 4.7% compared to a baseline least-square-based polynomial regression. This is a 10% improvement over the constant speed/acceleration model, a conventional speed predictor.

Automation & Control Systems↗

DOE OSTI2020

Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials (Second Annual Progress Report)

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of nuclear reactor parts. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process, which can introduce internal material flaws, such as pores and anisotropy. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures due to exposure to high temperature, radiation and corrosive environment in a nuclear reactor. Thermal tomography (TT) provides a capability for non-destructive evaluation of sub-surface defects in arbitrary size structures. Thermal tomography is a computational method for heat diffusion-based imaging of solids, which provides 3D visualization of data from flash thermography measurements. We investigate thermal tomography imaging and nondestructive evaluation of stainless steel and nickel super alloy metallic structures produced with laser powder bed fusion (LPBF) additive manufacturing (AM) process. Metallic structures produced with LPBF contain defects, and there are limited capabilities to evaluate these structures non-destructively. Thermal tomography reconstruction of 3D apparent spatial effusivity provides information about AM structure geometry and internal material flaws. We study performance of thermal tomography in imaging of metallic structures through COMSOL computer simulations of transient heat transfer, and through reconstruction of data obtained from experimental measurements. Reconstruction of internal defects is investigated using a stainless steel 316L specimen with flat bottom hole (FBH) indentations, and Inconel 718 plate produced with laser powder bed fusion (LPBF) method, which contains imprinted hemispherical shape low density regions containing non-sintered metallic powder. The FBH’s have the same sizes as the imprinted defects in the LPBF specimens, but offer better imaging contrast. Thermal tomography reconstructions provide visualizations of internal defects, and allow for estimation of their sizes and locations. Detection sensitivity of TT is limited by noises. We investigate separation of signal from noise in thermography images using several machine learning (ML) methods, including new spatio-temporal blind source separation (STBSS) and spatio-temporal sparse dictionary learning (STSDL) methods. Performance of the ML methods is benchmarked using thermography data obtained from imaging stainless steel 316L and Inconel 718 specimens produced LPBF method with imprinted calibrated porosity defects. The ML methods are ranked by F-score and execution runtime. Finally, we investigate TT of AM stainless steel 316L specimen with imprinted internal porosity defects using relatively low-cost, small form factor infrared (IR) camera based on uncooled micro bolometer detector. Sparse coding related K-means singular value decomposition (SVD) machine learning, image processing algorithms are developed to improve quality of TT images through removal of Additive white Gaussian noise without blurring the images. Following initial qualification of an AM component for deployment in a nuclear reactor, a compact TT can also be used for in-service nondestructive evaluation (NDE). With capability to perform in-service NDE of the AM component lifecycle, TT data can be used for development of a component digital twin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Reducing Memory Consumption in Calico with Shared Memory

This document details the work to reduce memory consumption in Calico. Calico is SimTools’ Constructive Solid Geometry (CSG) and geometry painting library. It is primarily used to paint material volume fractions in the Eulerian meshes of the physics codes. Calico provides point-in-body checks for the geometry supplied by an Oso model, which are then aggregated by the host codes. In addition, Calico can be used to build Oso models and is used by Ingen for that purpose. Oso models, and thus Calico, provide support for various CSG primitives such as spheres, cylinders, surfaces generated by rotating tabular curve data, and STL files as well as binary combinations of those primitives. Prior to refactoring Calico will run out of memory on CTS-1 machines when 36 MPI ranks are used per node when reading STL models on the order of 1.5 GB. This limitation is a bottleneck in designer workflow. This problem has been alleviated through the use of data structures to both reduce memory consumption and to leverage MPI-3 shared memory. This report details the data structures targeted for refactoring in Calico, the methods and implementation details for reducing memory consumption and leveraging shared memory, and results for one test problem. Results show a memory reduction when loading a 1 GB STL file by a factor of 27.5, from 93.4 to 3.4 GB.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Light-Front Field Theory on Current Quantum Computers

We present a quantum algorithm for simulation of quantum field theory in the light-front formulation and demonstrate how existing quantum devices can be used to study the structure of bound states in relativistic nuclear physics. Specifically, we apply the Variational Quantum Eigensolver algorithm to find the ground state of the light-front Hamiltonian obtained within the Basis Light-Front Quantization (BLFQ) framework. The BLFQ formulation of quantum field theory allows one to readily import techniques developed for digital quantum simulation of quantum chemistry. This provides a method that can be scaled up to simulation of full, relativistic quantum field theories in the quantum advantage regime. As an illustration, we calculate the mass, mass radius, decay constant, electromagnetic form factor, and charge radius of the pion on the IBM Vigo chip. This is the first time that the light-front approach to quantum field theory has been used to enable simulation of a real physical system on a quantum computer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Optimal Renewable Resource Allocation and Load Scheduling of Resilient Communities

This paper presents a methodology for enhancing community resilience through optimal renewable resource allocation and load scheduling in order to minimize unserved load and thermal discomfort. The proposed control architecture distributes the computational effort and is easier to be scaled up than traditional centralized control. The decentralized control architecture consists of two layers: The community operator layer (COL) allocates the limited amount of renewable energy resource according to the power flexibility of each building. The building agent layer (BAL) addresses the optimal load scheduling problem for each building with the allowable load determined by the COL. Both layers are formulated as a model predictive control (MPC) based optimization. Simulation scenarios are designed to compare different combinations of building weighting methods and objective functions to provide guidance for real-world deployment by community and microgrid operators. The results indicate that the impact of power flexibility is more prominent than the weighting factor to the resource allocation process. Allocation based purely on occupancy status could lead to an increase of PV curtailment. Further, it is necessary for the building agent to have multi-objective optimization to minimize unserved load ratio and maximize comfort simultaneously.

model predictive control↗

DOE OSTI2022

The Global Structure of the Milky Way’s Stellar Halo Based on the Orbits of Local Metal-poor Stars

We analyze the global structure of the Milky Way (MW)'s stellar halo, including its dominant subcomponent, Gaia-Sausage-Enceladus (GSE). The method for reconstructing the global distribution of this old stellar component is to employ the superposition of the orbits covering the large MW's space, where each of the orbit-weighting factors is assigned following the probability that the star is located at its currently observed position. The selected local, metal-poor sample with [Fe/H] <–1, using Gaia Early Data Release 3 and Sloan Digital Sky Survey Data Release 16, shows that the global shape of the halo is systematically rounder at all radii in more metal-poor ranges, such that an axial ratio, q, is nearly 1 for [Fe/H] <–2.2 and ~0.7 for –1.4 < [Fe/H] < –1.0. It is also found that a halo in the relatively metal-rich range of [Fe/H] >–1.8 actually shows a boxy/peanut-like shape, suggesting a major merger event. The distribution of azimuthal velocities shows a disk-like flattened structure at –1.4 < [Fe/H] < –1.0, which is thought to be the metal-weak thick disk. For the subsample of stars showing GSE-like kinematics, at [Fe/H] >–1.8, its global density distribution has an axis ratio of 0.9, which is more spherical than the general halo sample, and an outer ridge at r ~ 20 kpc. This spherical shape is consistent with the features of accreted halo components, and the ridge suggests that the orbit of GSE's progenitor had an apocenter of ~20 kpc. Implications for the formation of the stellar halo are also presented.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Results of ISO/TS 6336-22 Evaluating Full Contact Zone

ISO/TS 6336-22 (Calculation of load capacity of spur and helical gears - Part 22: Calculation of micropitting load capacity) is the ISO technical specification containing a proposal for calculations of the risk of micropitting in gear sets. Micropitting is a Hertzian fatigue phenomenon that appears as ultra-fine cracking and pitting on the flanks of gear teeth. Since progressive micropitting can lead to macropitting and flank damage, critical applications such as wind turbines, marine drives, and high-speed gear drives seek to accurately predict whether their designs are susceptible to this damage. ISO/TS 6336-22 assesses micropitting risk through a safety factor which is calculated as the minimum specific film thickness in the contact zone divided by a permissible specific film thickness. In the previous paper, the calculations were performed using the simplified method (Method B) that evaluates points on the path of contact. This was done for three gear sets that experienced micropitting in operation. The minimum specific film thickness for the two field cases was very high, which indicates that the gears were operating in the full elastohydrodynamic lubrication (EHL) regime. A more accurate calculation for these cases (Method A) calculates the specific film thickness across the entire contact zone. This paper applies this Method to the case study from the previous paper. The results are compared to micropitting observed in operation. Results are also compared to the results of the previous paper. Conclusions are made regarding the accuracy of both Methods compared to the field cases and relative to each other.

ENGINEERING,WIND ENERGY↗

DOE OSTI2022

Results of ISO/TS 6336-22 Evaluating Full Contact Zone: Preprint

ISO/TS 6336-22 (Calculation of load capacity of spur and helical gears — Part 22: Calculation of micropitting load capacity) is the ISO technical specification containing a proposal for calculations of the risk of micropitting in gear sets. Micropitting is a Hertzian fatigue phenomenon that appears as ultra-fine cracking and pitting on the flanks of gear teeth. Since progressive micropitting can lead to macropitting and flank damage, critical applications such as wind turbines, marine drives, and high-speed gear drives seek to accurately predict whether their designs are susceptible to this damage. ISO/TS 6336-22 assesses micropitting risk through a safety factor which is calculated as the minimum specific film thickness in the contact zone divided by a permissible specific film thickness. In the previous paper, the calculations were performed using the simplified method (Method B) that evaluates points on the path of contact. This was done for three gear sets that experienced micropitting in operation. The minimum specific film thickness for the two field cases was very high, which indicates that the gears were operating in the full elastohydrodynamic lubrication (EHL) regime. A more accurate calculation for these cases (Method A) calculates the specific film thickness across the entire contact zone. This paper applies this Method to the case study from the previous paper. The results are compared to micropitting observed in operation. Results are also compared to the results of the previous paper. Conclusions are made regarding the accuracy of both Methods compared to the field cases and relative to each other.

gearbox↗

DOE OSTI2022

BETO 2021 Peer Review - 1.3.2.001 - Algae Biomass Composition

Addressing critical improvements in biomass productivity and associated biochemical composition is a priority for the economic and sustainable commercial development of biofuels and bioproducts from algae. Capitalizing on pathways that integrate engineering approaches with fundamental biochemistry of photosynthetic organisms will lead to a better understanding of the complex nexus of algae growth rates, productivity and composition. This project focuses on identifying the critical factors for economic development of fuel and bioproduct technologies. Algal compositional characteristics form the foundation of robust economic and business models. This project supports that foundation by developing and validating accurate compositional methods and disseminating them to the greater community. Simultaneously, we build a deep understanding of the dynamic biochemical composition and carbon allocation for biomass value and conversion yields. A co-product portfolio developed under this project demonstrated a 30% increase in intrinsic value. Additionally, an integrated pipeline of molecular diversity mapping for product discovery with quantitative demonstration across species was deployed over the BETO algae program. The advances made here are highly relevant to BETO's multi-year program targets of reducing costs and integrating dynamic biomass composition with conversion processes to provide options for bioproducts, all leveraging the molecular diversity of algae.

algal compositional characteristics↗

DOE OSTI2022

Results of ISO/TS 6336-22 Evaluating Full Contact Zone

gearbox↗

DOE OSTI2024

Systems and methods for assessing cable insulation

Methods for assessing cable insulation are provided that can include engaging a first tine gap interdigital capacitor assembly with the exterior of the cable jacket to determine a first capacitance and dissipation factor value of the insulation using the first tine gap interdigital capacitor assembly; engaging a second tine gap interdigital capacitor assembly with the exterior of the cable jacket to determine a second capacitance and dissipation factor value of the insulation using the second tine gap interdigital capacitor assembly.Systems for assessing cable insulation are also provided. The systems can include a cable jacket encompassing cable insulation and conductive material; at least two tine gap interdigital capacitor assemblies engaged with the cable jacket; and processing circuitry operatively engaged with the capacitor assemblies.

Glass, III, Samuel W.↗

DOE OSTI2025

EVALUATION OF HRA METHODOLOGIES FOR APPLICATION IN SDP WORK

This study critically evaluates human reliability analysis (HRA) methodologies applicable to regulatory probabilistic safety assessment (PSA) model, with a particular focus on their role in supporting the significance determination process (SDP) in nuclear safety assessment. Firstly, three widely utilized HRA methods – IDHEAS-ECA, SPAR-H, and ASEP/THERP – were qualitatively and quantitatively assessed. Qualitative assessments were conducted using attributes from the NEA/CSNI/R(2015)1 report, while quantitative evaluations employed regression and correlation analyses to compare predicted human error probabilities (HEPs) against empirical data. Results reveal distinct strengths, for example, IDHEAS-ECA’s robust predictive accuracy and K-HRA’s alignment with operational practices. In addition, dependency analysis and recovery analysis were critically evaluated. For dependency analysis, the methods’ handling of inter-task dependencies and their impact on HEPs were examined, while recovery analysis highlighted strategies for mitigating failure events. Furthermore, strategies were proposed to evaluate performance-shaping factors under conditions of reduced human performance, such as stress, fatigue, or cognitive overload, addressing specific challenges faced in SDP evaluations. Human errors from KINS’s operational performance information system event reports were evaluated as a case study. This study identifies gaps and provides actionable insights to ensure their validity and applicability in SDP HRA applications. This paper is a part of research conducted by KINS, and it should be noted that this result does not represent the regulatory position of KINS.

99 - GENERAL AND MISCELLANEOUS↗

DOE OSTI2022

Artificial intelligence inferred microstructural properties from voltage–capacity curves

Abstract The quantification of microstructural properties to optimize battery design and performance, to maintain product quality, or to track the degradation of LIBs remains expensive and slow when performed through currently used characterization approaches. In this paper, a convolution neural network-based deep learning approach (CNN) is reported to infer electrode microstructural properties from the inexpensive, easy to measure cell voltage versus capacity data. The developed framework combines two CNN models to balance the bias and variance of the overall predictions. As an example application, the method was demonstrated against porous electrode theory-generated voltage versus capacity plots. For the graphite|LiMn $$_2$$ 2 O $$_4$$ 4 chemistry, each voltage curve was parameterized as a function of the cathode microstructure tortuosity and area density, delivering CNN predictions of Bruggeman’s exponent and shape factor with 0.97 $$R^2$$ R 2 score within 2 s each, enabling to distinguish between different types of particle morphologies, anisotropies, and particle alignments. The developed neural network model can readily accelerate the processing-properties-performance and degradation characteristics of the existing and emerging LIB chemistries.

25 ENERGY STORAGE↗

DOE OSTI2025

Dehydroxylation kinetics of kaolinite and montmorillonite examined using isoconversional methods

The use of calcined clays as supplementary cementitious materials (SCMs) in concrete is a promising strategy towards decarbonizing the cement and concrete industry. This is especially relevant considering the ever-increasing demand for concrete. Comprehensive understanding of the kinetics of calcination is essential towards maximizing the potential reactivity of clay minerals while ensuring energy efficiency. In this study, the kinetics of the dehydroxylation of kaolinite and montmorillonite are investigated under non-isothermal conditions at constant heating rate. Activation energies ( E a ) are determined via Friedman differential and advanced Vyazovkin incremental methods over the isoconversional range; these are devoid of computational approximations, thus allowing kinetic analysis without assuming a specific reaction model. Kinetic equations—in the differential form as well as a combination of differential and integral forms are compared against the experimentally determined reaction models to identify the most probable dehydroxylation mechanism for kaolinite and montmorillonite. A reaction order mechanism is established for dehydroxylation of kaolinite, while montmorillonite is noted to undergo dehydroxylation via a single-step reversible diffusion-controlled process. Kinetic triplet—comprising activation energy, reaction model and pre-exponential factor—is used to predict isothermal calcination conditions, which is further verified using analytical techniques. Heat release rates of clay-portlandite blends from isothermal calorimetry are used within a thermodynamic framework to quantify reactivity of the calcined clays. Here, the study demonstrates a general approach based on isoconversional methods to predict calcination conditions for different clays that can be used in efficient and optimized production of blended cements or SCMs.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Accelerating multigrid with streaming chiral SVD for Wilson fermions in lattice QCD

A modification to the setup algorithm for the multigrid preconditioner of Wilson fermions in lattice QCD is presented. A larger basis of test vectors than that used in regular multigrid is calculated by the smoother and truncated by singular value decomposition on the chiral components of the test vectors. The truncated basis is used to form the prolongation and restriction matrices of the multigrid hierarchy. This modification of the setup method is demonstrated to increase the convergence of linear solvers on an anisotropic lattice with m π ≈ 239 MeV from the Hadron Spectrum Collaboration and an isotropic lattice with m π ≈ 220 MeV from the MILC Collaboration. The lattice volume dependence of the method is also examined. Increasing the number of test vectors improves speedup up to a point, but storing these vectors becomes impossible in limited memory resources such as GPUs. To address storage cost, we implement a streaming singular value decomposition of the basis of test vectors on the chiral components and demonstrate a decrease in the number of fine level iterations by a factor of 1.7 for m q ≈ m crit

Iterative methods↗