SEARCH · Search NASA

Results for “SVD”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Use of Fisher's Ratio assisted multivariate curve resolution- alternating least squares for discovery-based analysis using ultrahigh pressure liquid chromatography-high resolution mass spectrometry

Non-targeted analysis of complex chemical mixtures can be difficult considering the convoluted nature of the matrix and the potential unknown chemical differences between samples or classes of samples. Ultrahigh pressure liquid chromatography coupled to quadrupole time-of-flight mass spectrometry (UHPLC-QTOF) is an ideal technique to probe chemical differences for a wide variety of samples. While UHPLC-QTOF can discover minute chemical differences down to low part per billion (ppb) concentrations with a high degree of confidence, the application of high-resolution mass spectrometry can yield massive amounts of information (∼ 10 gb per sample) that cannot be analyzed manually. Therefore, the application of chemometric techniques is mandatory for the interrogation of complex samples. Fisher's ratio (FR) assisted multivariate curve resolution-alternating least squares (MCR-ALS) was used to the discover and identify the chemical differences between two classes of materials: 1) a pond water matrix and 2) the matrix spiked with a pharmaceutical standard mix containing 17 compounds. Thirteen of the seventeen spiked compounds were discovered using FR analysis, and then five were successfully deconvoluted using MCR-ALS wherein the number of curves chosen were automatically determined using singular value decomposition (SVD). In conclusion, the use of an automated FR assisted MCR-ALS will aid in discovering trace levels of chemical components without the need for the researcher to provide potentially biased input which will aid in non-targeted workflow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Permutation-adapted complete and independent basis for atomic cluster expansion descriptors

Atomic cluster expansion (ACE) methods provide a systematic way to describe particle local environments of arbitrary body order. For practical applications it is often required that the basis of cluster functions be symmetrized with respect to rotations and permutations. Existing methodologies yield sets of symmetrized functions that are over-complete. These methodologies thus require an additional numerical procedure, such as singular value decomposition (SVD), to eliminate redundant functions. In this work, it is shown that analytical linear relationships for subsets of cluster functions may be derived using recursion and permutation properties of generalized Wigner symbols. From these relationships, subsets (blocks) of cluster functions can be selected such that, within each block, functions are guaranteed to be linearly independent. It is conjectured that this block-wise independent set of permutation-adapted rotation and permutation invariant (PA-RPI) functions forms a complete, independent basis for ACE. Along with the first analytical proofs of block-wise linear dependence of ACE cluster functions and other theoretical arguments, numerical results are offered to demonstrate this. The utility of the method is demonstrated in the development of an ACE interatomic potential for tantalum. Using the new basis functions in combination with Bayesian compressive sensing sparse regression, some high degree descriptors are observed to persist and help achieve high-accuracy models.

Angular momentum↗

DOE OSTI2022

Harmonizing solar induced fluorescence across spatial scales, instruments, and extraction methods using proximal and airborne remote sensing: A multi-scale study in a soybean field

Solar-induced chlorophyll fluorescence (SIF) has been widely used to track vegetation photosynthesis at different scales ranging from in-situ measurements to satellite products. Airborne platforms sample SIF data at a spatial scale intermediate between in-situ and satellite, matching that of ground measurement (e.g. flux tower footprints and other field sampling), enabling us to explore causes of SIF variation and validate satellite-based SIF products. However, harmonizing SIF across sensors and platforms (correcting for systematic errors to yield a consistent, comparable SIF product) is challenging because SIF can be retrieved in different absorption windows, with different instruments and methods complicating the comparison between different observational levels (i.e., ground, airborne, satellites) and between sites equipped with different instruments with varying optical properties (spectral resolution and sampling intervals, spatial resolution). Additionally, the spatial and temporal variability of atmospheric properties can influence the retrieval of the weak SIF signal. Because of these complications, direct comparisons of airborne and ground SIF across scales are rarely attempted. Here, in this study, we combined airborne SIF data with simultaneous ‘ground truth’ data collected by stationary and mobile platforms in a soybean field in Nebraska, USA. In this effort, we tested several SIF extraction methods, including Fraunhofer Line Discrimination (FLD), improved Fraunhofer Line Discrimination (iFLD), Spectral Fitting Method (SFM), SpecFit, and a Singular Vector Decomposition (SVD) method. The SpecFit method was sensitive to the 715–740 nm water bands and removing the water bands in the fitting process yielded better agreement between the airborne and ground SIF spectra. Accurate estimation of the ground level downwelling irradiance obtained by ground measurements over a calibration target improved agreement between airborne and ground SIF retrievals at the O 2 A band, and allowed us to derive a SIF dataset with improved agreement across platforms and sampling scales. This experimental approach provided a method for generating comparable SIF signals across instruments, methods and platforms, which is critical to understanding the SIF-GPP relationship at different scales and to cross-validate the diversity of platforms used for satellite products calibration and validation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Optimizing the differential connection schemes for detecting 3D magnetic perturbations in DIII-D

Optimizing the differential pair connection scheme (i.e. the set of pairs) of a toroidal array of magnetic sensors dedicated to measuring slowly rotating asymmetric fields can enhance the mode number detection capability and failure-resilience. In this work the condition number obtained from singular value decomposition (SVD) of the design matrix is used as a metric to evaluate the quality of a connection scheme. A large number of possible pair connections are usually available, so evaluating all of them may require extensive use of computational resources and can be very time consuming. Alternative methods to reduce the number of pairs evaluated without losing the capabilities of toroidal mode detection are presented in this paper. Here, three examples of the applications of such analysis for the 3D magnetic diagnostic system of DIII-D are also presented: the addition of two new toroidal arrays with n > 3 detection capabilities, the modification of an existing toroidal array in the low field side of the machine to accommodate the addition of a helicon antenna, and the design of changes in several toroidal arrays in the high field side to accommodate the addition of a lower hybrid current drive antenna on the center post.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Antineutrino energy spectrum unfolding based on the Daya Bay measurement and its applications

The prediction of reactor antineutrino spectra will play a crucial role as reactor experiments enter the precision era. The positron energy spectrum of 3.5 million antineutrino inverse beta decay reactions observed by the Daya Bay experiment, in combination with the fission rates of fissile isotopes in the reactor, is used to extract the positron energy spectra resulting from the fission of specific isotopes. This information can be used to produce a precise, data-based prediction of the antineutrino energy spectrum in other reactor antineutrino experiments with different fission fractions than Daya Bay. The positron energy spectra are unfolded to obtain the antineutrino energy spectra by removing the contribution from detector response with the Wiener-SVD unfolding method. Consistent results are obtained with other unfolding methods. A technique to construct a data-based prediction of the reactor antineutrino energy spectrum is proposed and investigated. Given the reactor fission fractions, the technique can predict the energy spectrum to a 2% precision. In addition, we illustrate how to perform a rigorous comparison between the unfolded antineutrino spectrum and a theoretical model prediction that avoids the input model bias of the unfolding method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Joint Determination of Reactor Antineutrino Spectra from U 235 and Pu 239 Fission by Daya Bay and PROSPECT

A joint determination of the reactor antineutrino spectra resulting from the fission of 235U and 239Pu has been carried out by the Daya Bay and PROSPECT collaborations. This Letter reports the level of consistency of 235U spectrum measurements from the two experiments and presents new results from a joint analysis of both data sets. The measurements are found to be consistent. The combined analysis reduces the degeneracy between the dominant 235U and 239Pu isotopes and improves the uncertainty of the 235U spectral shape to about 3%. The 235U and 239Pu antineutrino energy spectra are unfolded from the jointly deconvolved reactor spectra using the Wiener-SVD unfolding method, providing a data-based reference for other reactor antineutrino experiments and other applications. This is the first measurement of the 235U and 239Pu spectra based on the combination of experiments at low- and highly enriched uranium reactors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

A Predictor-Corrector Strategy for Adaptivity in Dynamical Low-Rank Approximations

Here, in this paper, we present a predictor-corrector strategy for constructing rank-adaptive, dynamical low-rank approximations (DLRAs) of matrix-valued ODE systems. The strategy is a compromise between (i) low-rank step-truncation approaches that alternately evolve and compress solutions and (ii) strict DLRA approaches that augment the low-rank manifold using subspaces generated locally in time by the DLRA integrator. The strategy is based on an analysis of the error between a forward temporal update into the ambient full-rank space, which is typically computed in a step-truncation approach before recompressing, and the standard DLRA update, which is forced to live in a low-rank manifold. We use this error, without requiring its full-rank representation, to correct the DLRA solution. A key ingredient for maintaining a low-rank representation of the error is a randomized SVD, which introduces some degree of stochastic variability into the implementation. The strategy is formulated and implemented in the context of discontinuous Galerkin spatial discretizations of PDEs and applied to several versions of DLRA methods found in the literature as well as a new variant. Numerical experiments comparing the predictor-corrector strategy to other methods demonstrate robustness to overcome shortcomings of step truncation or strict DLRA approaches: The former may require more memory than is strictly needed, while the latter may miss transients solution features that cannot be recovered. The effect of randomization, tolerances, and other implementation parameters is also explored.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Parallel Randomized Tucker Decomposition Algorithms

The Tucker tensor decomposition is a natural extension of the singular value decomposition (SVD) to multiway data. Here, we propose to accelerate Tucker tensor decomposition algorithms by using randomization and parallelization. We present two algorithms that scale to large data and many processors, significantly reduce both computation and communication cost compared to previous deterministic and randomized approaches, and obtain nearly the same approximation errors. The key idea in our algorithms is to perform randomized sketches with Kronecker-structured random matrices, which reduces computation compared to unstructured matrices and can be implemented using a fundamental tensor computational kernel. We provide probabilistic error analysis of our algorithms and implement a new parallel algorithm for the structured randomized sketch. Our experimental results demonstrate that our combination of randomization and parallelization achieves accurate Tucker decompositions much faster than alternative approaches. We observe up to a 16X speedup over the fastest deterministic parallel implementation on 3D simulation data.

Tucker decompositions↗

DOE OSTI2020

Ultrafast structural response to charge redistribution within a photosinthetic reaction centre

This data set contains stream files, difference maps and SVD maps obtained from data processing of time-resolved serial femtosecond crystallography measurements performed at LCLS XFEL, used to construct the images and perform the analysis published in the paper "Ultrafast structural response to charge redistribution within a photosinthetic reaction centre". Two (2) time points plus a dark state from 2015 ('a' label) as well as five (5) time points plus a dark state from 2016 ('b' label) are provided.

CXI↗

DOE OSTI2023

DT-HYDRO

The software solves the time dependent, one dimensional (1D) coupled mass and momentum balance equations governing the elastic flow of water through the penstock, turbine and draft tube in a hydroelectric facility using a high order finite volume based method. The numerical method is based on the Kurganov-Tadmor central method paired with the Monotonic Upstream-centered Scheme for Conservation Laws (MUSCL). This solution method accurately resolves the fast transient behavior of the flow, including water hammer. Additionally, the software estimates the full 3D flow field within the turbine chamber in real time, a feat that is made possible by leveraging pre-computed CFD results by utilizing a reduced order modeling method based on an efficient randomized singular value decomposition (SVD) driven proper orthogonal decomposition (POD) with POD-mode weight regression. The reduced order model of the 3D flow is directly coupled to the 1D elastic flow model so the entire flow field through the penstock and turbine system is resolved quickly and with high fidelity.

Gurecky, William [Oak Ridge National Laboratory (O↗

DOE OSTI2025

AutoDiffSVDCompression [SWR-25-162]

This software contains the code and results for the paper Reducing Memory Usage of Reverse Mode AD by Compression. In particular, the code includes implementation of the necessary structures and functions for using the Singular Value Decomposition (SVD) to compress forward pass values needed for reverse mode AD as well as the necessary numerical solvers for Burgers’ Equation and numerical solution of a control problem. In addition to the code, the repository also contains the scripts, results, and logs used in the aforementioned publication. This repository is meant to enable the reproducibility of the results in the publication.

Maack, Jonathan [National Renewable Energy Laborat↗

DOE OSTI2025

Geometric Interpretation of the Cluster Location Problem Part II: Application to the Pahala, Hawaii, Earthquake Sequence

In the companion “Theory” article, we presented a new framing of the seismic location problem in terms of differential geometry (Harris et al., 2025). From that viewpoint, we developed a “project and correct” approach for estimating the relative locations of earthquakes. Here, in this study, we use project and correct to estimate high-precision relative locations of events from an earthquake sequence beneath the town of Pahala, Hawaii, using high-precision correlation-derived picks. The sequence was active from 2020 through 2022 and produced many highly correlated signals at Hawaii Volcano Observatory (HVO) stations on the island of Hawaii. The data we inverted consisted of 2882 events with observations at 5 HVO stations. For comparison with the travel-time image, we also produced conventional hypocenter solutions using both the Bayesloc program (Myers et al., 2007, 2009) and a purpose-built double-difference code. There were obvious structural elements in the resulting image, the resolution of which we used to test the performance of the project and the correct algorithm. For the projection step, we first produced a 3D local basis using an singular value decomposition (SVD) of the 2882 groups of times. Projection of the travel-time vectors into this basis resulted in an image with structures similar to those produced by our conventional locators, but with distortion as predicted by theory. Removing the distortion requires an inverse operator generated from the metric tensor at the geometric centroid of the events. We compared two approaches to obtaining such an inverse operator. The first uses an estimate of the geographic centroid of the event cloud from the centroid of the travel-time data. The second approach uses the centroid of the conventionally produced locations. The first approach produces a corrected image very similar to the conventional results, but with a rotation. The corrected image produced using the conventionally derived centroid is a near-exact match to the conventional locations.

Dodge, Douglas A. [Lawrence Livermore National Lab↗

DOE OSTI2026

ORBIT CORRECTION STUDIES ON THE ELECTRON TRANSPORT LINE FROM RCS TO ESR

A dedicated electron transfer line from Rapid Cycling Synchrotron (RCS) to Electron Storage Ring (ESR), referred to as the RTE line has been designed for the Electron-Ion Collider (EIC). The beamline follows a straight-line geometry, with a length of 133 m, and is consists with two matching sections and a FODO section for beam diagnostics. Imperfections with magnet alignments introduce orbit distortions, making orbit correction scheme a critical component in the design. To facilitate orbit correction, each quadrupole magnet is equipped with a pair of beam position monitors (BPMs) and kickers. The Singular Value Decomposition (SVD) algorithm is used for orbit correction and tolerance studies. This paper presents the ongoing progress in the optics design and error correction scheme of the RTE line.

Neththikumara, I. [Thomas Jefferson National Accel↗

DOE OSTI2020

Development of Thermal Tomography Imaging System for In-Service Nondestructive Evaluation of Additively Manufactured Components

36 MATERIALS SCIENCE↗

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 OSTI2021

(U) A Linear Response Model Predicts Reactivity From a Density Profile

We tested the ability to predict the system reactivity, described by alpha, given a density profile using a simple linear system response. We generated a suite of 1-dimensional density profiles that consisted of nominal density, a discontinuity, and a decay. These profiles were prescribed a functional form and the mass was conserved in all cases. From these density profiles, we calculated the alpha value of the 3-dimensional system.We calculated a linear response function given a training set of the 1-dimensional density profiles, and the system reactivity described by alpha. We tested the robustness of the response function using the remaining test data. Our results showed very good agreement between the predicted and calculated test values, where the distribution of alpha differences was centered about zero and had a standard deviation of 0.005 gens/shake. The predicted and calculated alpha values did not significantly differ (t=-0.0009 p<0.99). We used Singular Value Decomposition (SVD) to reduce the matrix rank by retaining95% of the cumulative singular value contributions. This reduced the matrix rank by 91.7%. We generated the linear response matrix and calculated the difference between the predicted and calculated alpha values. Using the reduced order matrix, we showed good agreement between the predicted and calculated alpha values where the distribution of differences was centered near zero, the standard deviation was 0.006 gens/shake, and the statistical t-test showed good agreement (t=0.02, p<0.98). These results show a linear relationship between a series of 1-dimensional density profiles,where the mass was conserved, and the system reactivity. The next steps of this work will be to investigate the linear response using 2-dimensional density profiles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Muon-Neutrino Charged-Current Cross Sections from MicroBooNE: First Simultaneous Measurements of Final States with and without Protons for Muon-Neutrino Scattering on Argon

A detailed understanding of muon neutrino charged-current interactions on argon is crucial to the study of neutrino oscillations in current and future experiments using liquid argon time projection chambers. To help fill this need, MicroBooNE has produced a comprehensive set of cross section measurements which simultaneously probe the leptonic and hadronic systems by dividing the inclusive channel into final states with and without protons. Data-driven model validation utilizing the conditional constraint formalism is employed to detect mismodeling that may bias the nominal flux averaged cross section results, which are extracted with the Wiener-SVD unfolding method. The results are compared to widely used event generator predictions revealing significant mismodeling of final states without protons, possibly due to insufficient treatment of final state interactions. These are first differential muon neutrino-argon cross section measurements made simultaneously for final states with and without protons and provide novel information that will help stimulate the improvement of event generator modeling.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Measurement of Neutral-Current $\pi^0$ Cross Section in MicroBooNE using Wire-Cell Reconstruction

This note presents the extraction of flux-averaged total and differential neutral-current $π^0$ cross section measurements, based on the Wiener-SVD unfolding method, using the MicroBooNE Liquid Argon Time Projection Chamber (LArTPC) and the Wire-Cell reconstruction paradigm. This measurement is critical for single-photon searches as neutral-current pion production represents one of the largest sources of background. We report the results for an exposure corresponding to 5.327 × 10 19 protons-on-target (POT).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗