Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “singular value decomposition”

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 109 records · Page 6
All sourcesNASA NTRSDOE 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 OSTI2025

Jacobian-based model diagnostics and application to equation oriented modeling of a carbon capture system

It can be difficult to identify the specific variables or equations responsible for convergence issues in large mathematical programming models. The Institute for the Design of Advanced Energy Systems Integrated Platform (IDAES-IP) contains a tool to identify poorly scaled constraints and variables by searching for rows and columns of the Jacobian matrix with small L2-norms. A singular value decomposition is then performed to identify degenerate sets of equations and remaining scaling issues. Here, this work presents a flowsheet developed for post-combustion carbon capture using a monoethanolamine (MEA) solvent system as a case study. This work takes the reader through the entire process of model diagnostics and reformulation, from a basic introduction to the mathematics behind these model diagnostics to the reformulations necessary to make the model numerically robust, including a significantly modified enhancement factor model.

IDAES↗
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

DESPERATE: A Python Library for Processing and Denoising NMR Spectra

NMR spectroscopy is an inherently insensitive technique with respect to the amount of observable signal. A common element in all NMR spectra is random thermal noise that is often characterized by a signal-to-noise ratio (SNR). SNR can be generically improved experimentally with repetitive signal averaging or during post-processing with apodization; the former of which often results in long experimental times and the latter results in the loss of spectral resolution. Denoising techniques can instead be used during post-processing to enhance SNR without compromising resolution. The most common approach relies on the singular-value decomposition (SVD) to discard noisy components of NMR data. SVD-based approaches work well, such as Cadzow and PCA, but are computationally expensive when used for large datasets that are often encountered in NMR (e.g., Carr-Purcell/Meiboom-Gill and nD datasets). Herein, we describe the implementation of a new wavelet transform (WT) routine for the fast and robust denoising of 1D and 2D NMR spectra. Several simulated and experimental datasets are denoised with both SVD-based Cadzow or PCA and WT’s, and the resulting SNR enhancements and spectral uniformity are compared. WT denoising offers similar and improved denoising compared with SVD and operates faster by several orders-of-magnitude in some cases. Further, all denoising and processing routines used in this work are included in a free and open-source Python library called DESPERATE.

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

Principal Component Analysis of azimuthal flow in intermediate-energy heavy-ion reactions

Principal Component Analysis (PCA) via Singular Value Decomposition (SVD) of large datasets is an adaptive exploratory method to uncover natural patterns underlying the data. Several recent applications of the PCA-SVD to event-by-event single-particle azimuthal angle distribution matrices in ultra-relativistic heavy-ion collisions at RHIC-LHC energies indicate that the sine and cosine functions chosen a priori in the traditional Fourier analysis are naturally the most optimal basis for azimuthal flow studies according to the data itself. We perform PCA-SVD analyses of mid-central Au+Au collisions at $E$ $beam$ / $A$ =1.23 GeV simulated using an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model to address the following two questions: (1) if the principal components of the covariance matrix of nucleon azimuthal angle distributions in heavy-ion reactions around 1 GeV/nucleon are naturally sine and/or cosine functions and (2) what if any advantages the PCA-SVD may have over the traditional flow analysis using the Fourier expansion for studying the EOS of dense nuclear matter. In conclusion, we find that (1) in none of our analyses the principal components come out naturally as sine and/or cosine functions, (2) while both the eigenvectors and eigenvalues of the covariance matrix are appreciably EOS dependent, the PCA-SVD has no apparent advantage over the traditional Fourier analysis for studying the EOS of dense nuclear matter using the azimuthal collective flow in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A learned closure method applied to phase mixing in a turbulent gradient-driven gyrokinetic system in simple geometry

We present a new method for formulating closures that learn from kinetic simulation data. We apply this method to phase mixing in a simple gyrokinetic turbulent system – temperature-gradient-driven turbulence in an unsheared slab. The closure, called the learned multi-mode (LMM) closure, is constructed by, first, extracting an optimal basis from a nonlinear kinetic simulation using singular value decomposition. Subsequent nonlinear fluid simulations are projected onto this basis and the results are used to formulate the closure. We compare the closure with other closures schemes over a broad range of the relevant two-dimensional parameter space (collisionality and gradient drive). We find that the turbulent kinetic system produces phase-mixing rates much lower than the linear expectations, which the LMM closure is capable of capturing. We also compare radial heat fluxes. A Hammett–Perkins closure, generalized to include collisional effects, is quite successful throughout the parameter space, producing ${\sim }14\,\%$ root-mean-square (r.m.s.) error. The LMM closure is also very effective: when trained at three (two) points (in a 35 point parameter grid), the LMM closure produces $8\,\%$ ( $12\,\%$ ) r.m.s. errors. The LMM procedure can be readily generalized to other closure problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Cluster Analysis of Combined EDS and EBSD Data to Solve Ambiguous Phase Identifications

A common problem in analytical scanning electron microscopy (SEM) using electron backscatter diffraction (EBSD) is the differentiation of phases with distinct chemistry but the same or very similar crystal structure. X-ray energy dispersive spectroscopy (EDS) is useful to help differentiate these phases of similar crystal structures but different elemental makeups. However, open, automated, and unbiased methods of differentiating phases of similar EBSD responses based on their EDS response are lacking. This paper describes a simple data analytics-based method, using a combination of singular value decomposition and cluster analysis, to merge simultaneously acquired EDS + EBSD information and automatically determine phases from both their crystal and elemental data. I use hexagonal TiB 2 ceramic contaminated with multiple crystallographically ambiguous but chemically distinct cubic phases to illustrate the method. Code, in the form of a Python 3 Jupyter Notebook, and the necessary data to replicate the analysis are provided as Supplementary material.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Net Present Value Optimization of a Natural Gas Combined Cycle Plant with CO 2 Capture using a Water-Lean Solvent Considering Transient Electricity Price for Multiple Regions

Global CO 2 emissions are increasing at about a 1.5% rate per year. Fossil fuel-based plants are one of the main contributors to this rise. In the power generation industry, fossil fuel plants are dominant, and many plants are under development. In this study, a natural gas combined cycle (NGCC) power plant with postcombustion capture using a leading water-lean solvent is considered. For optimal design and operating schedule, large-scale dynamic optimization is undertaken for net present value (NPV) optimization. The first principle dynamic model of NGCC is developed, including a model of the highly efficient H-class gas turbines. For computational tractability of the dynamic optimization problem, a reduced-order model is developed by using the Hankel singular value decomposition. A waterlean solvent, N-(2-ethoxyethyl)-3-morpholinopropan-1-amine, is used for carbon capture. A model of the capture system is developed in Aspen Plus, which is used to develop a reduced-order model by using ALAMO, a machine learning software. In addition, a reduced model of the CO 2 compression system with a dehydration unit is also considered. The integrated system is used for NPV optimization by using the Python-based PYOMO platform. The PCC process is analyzed for three configurations-conventional packed bed, rotating packed bed (RPB), and a combination of RPB and direct contact cooler. The NPV optimization is performed for 14 regional markets by considering year-long clustered and continuous locational marginal price data with a 1 h interval. Optimization results show that the PCC can achieve 90% CO 2 capture with a positive NPV for six regions. Sensitivity studies conducted by using the PCC configurations indicate that the process is economically feasible for 9 regions out of 14 regional electricity markets with NPV values in the range of 33−540 $MM.

cabon capture↗
DOE OSTI2026

Learning Latent Representations to Bridge Coarse-Grained and Atomistic Resolutions in Polymer Simulations

We present a machine-learning-based framework for learning reduced-order representations of polymer chain conformations across coarse-grained (CG) and united-atom (UA) fidelities. By employing linear singular value decomposition and nonlinear autoencoders, we compress high-dimensional polymer configurations into latent spaces with minimal loss of structural accuracy. Crucially, we demonstrate a near-perfect linear mapping between CG and UA latent spaces, enabling an efficient super-resolution back-mapping procedure that reconstructs high-fidelity UA configurations from CG simulations. While minor structural inaccuracies occur, they are effectively corrected through a brief molecular dynamics relaxation, forming a practical hybrid machine learning−physics scheme. This approach establishes the key structural prerequisites for accelerated polymer dynamics simulations: a compact and accurate latent encoding of polymer chain conformations and a validated multi-fidelity mapping that permits reconstruction of UA structures from CG configurations. The extension of this framework to explicit time evolution within the latent space, enabling dynamics to be propagated at CG fidelity and decoded to UA resolution only when required, represents a natural and well-motivated direction for future work.

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

Heterogeneity in M. tuberculosis β-lactamase inhibition by Sulbactam

Abstract For decades, researchers have elucidated essential enzymatic functions on the atomic length scale by tracing atomic positions in real-time. Our work builds on possibilities unleashed by mix-and-inject serial crystallography (MISC) at X-ray free electron laser facilities. In this approach, enzymatic reactions are triggered by mixing substrate or ligand solutions with enzyme microcrystals. Here, we report in atomic detail (between 2.2 and 2.7 Å resolution) by room-temperature, time-resolved crystallography with millisecond time-resolution (with timepoints between 3 ms and 700 ms) how the Mycobacterium tuberculosis enzyme BlaC is inhibited by sulbactam (SUB). Our results reveal ligand binding heterogeneity, ligand gating, cooperativity, induced fit, and conformational selection all from the same set of MISC data, detailing how SUB approaches the catalytic clefts and binds to the enzyme noncovalently before reacting to a trans- enamine. This was made possible in part by the application of singular value decomposition to the MISC data using a program that remains functional even if unit cell parameters change up to 3 Å during the reaction.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Classification of gaseous UF 6 assay by femtosecond LIBS in the 424.4 nm spectral region using numerical HOGSVD-DTW features

This technical note presents experimental results using numerical features of fs-LIBS data to classify the assay value of a gaseous UF 6 material. Here, the data-driven feature vectors are computed by Higher Order Generalized Singular Value Decomposition (HOGSVD) and Dynamic Time Warp (DTW). The method achieves 96.97% accuracy in spectral classification testing with fs-LIBS samples obtained from a UF 6 material with five known assay values ranging from 0.287% to 61.740%, with 100% accuracy for the four largest assay values ranging from 4.615% to 61.740%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

The overlapping fragment approach for non-orthogonal configuration interaction with fragments

The non-orthogonal configuration interaction with fragments (NOCI-F) approach is extended opening the possibility to study intramolecular processes and materials with covalent or ionic lattices. So far, NOCI-F has been applied to study intermolecular energy and electron transfer employing ensembles of fragments that do not have atoms or bonds in common. The here presented approach divides the target system into two overlapping fragments that share one or more atoms and/or one or more bonds. After the construction of a collection of (multiconfigurational) fragment wave functions in a state specific optimization procedure, the fragment wave functions are combined to form many-electron basis functions for the non-orthogonal configuration interaction of the whole system. The orbitals in the overlapping fragment are defined by a corresponding orbital transformation of the fragment orbitals through a singular value decomposition. The overlapping fragments approach is first illustrated for a model system and then used to highlight some possible applications of NOCI with overlapping fragments. In conclusion, the results of excited state diffusion in transition metal oxide, intramolecular singlet fission and magnetic interactions in organic biradicals and ionic compounds are discussed.

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

Deconvolution of dynamic heterogeneity in protein structure

Heterogeneity is intrinsic to the dynamic process of a chemical reaction. As reactants are converted to products via intermediates, the nature and extent of heterogeneity vary temporally throughout the duration of the reaction and spatially across the molecular ensemble. The goal of many biophysical techniques, including crystallography and spectroscopy, is to establish a reaction trajectory that follows an experimentally provoked dynamic process. It is essential to properly analyze and resolve heterogeneity inevitably embedded in experimental datasets. We have developed a deconvolution technique based on singular value decomposition (SVD), which we have rigorously practiced in diverse research projects. In this review, we recapitulate the motivation and challenges in addressing the heterogeneity problem and lay out the mathematical foundation of our methodology that enables isolation of chemically sensible structural signals. We also present a few case studies to demonstrate the concept and outcome of the SVD-based deconvolution. Finally, we highlight a few recent studies with mechanistic insights made possible by heterogeneity deconvolution.

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

Exploring Hilbert space on a budget: Novel benchmark set and performance metric for testing electronic structure methods in the regime of strong correlation

This work explores the ability of classical electronic structure methods to efficiently represent (compress) the information content of full configuration interaction (FCI) wave functions. We introduce a benchmark set of four hydrogen model systems of different dimensionalities and distinctive electronic structures: a 1D chain, a 1D ring, a 2D triangular lattice, and a 3D close-packed pyramid. To assess the ability of a computational method to produce accurate and compact wave functions, we introduce the accuracy volume, a metric that measures the number of variational parameters necessary to achieve a target energy error. Using this metric and the hydrogen models, we examine the performance of three classical deterministic methods: (i) selected configuration interaction (sCI) realized both via an a posteriori (ap-sCI) and variational selection of the most important determinants, (ii) an a posteriori singular value decomposition (SVD) of the FCI tensor (SVD-FCI), and (iii) the matrix product state representation obtained via the density matrix renormalization group (DMRG). We find that the DMRG generally gives the most efficient wave function representation for all systems, particularly in the 1D chain with a localized basis. For the 2D and 3D systems, all methods (except DMRG) perform best with a delocalized basis, and the efficiency of sCI and SVD-FCI is closer to that of DMRG. For larger analogs of the models, the DMRG consistently requires the fewest parameters but still scales exponentially in 2D and 3D systems, and the performance of SVD-FCI is essentially equivalent to that of ap-sCI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2022

Comparing the shock sensitivity of insensitive energetic materials

We present a tabletop method to study the shock sensitivity of plastic-bonded explosives that are considered shock insensitive using high dynamic range optical emission spectroscopy with laser-launched km/s flyer plate impacts (2, 3, and 4 km/s), which measures the spectral radiance (the emission spectrum vs a calibrated intensity standard) every 0.8 ns in the nanosecond and microsecond regimes. The explosives were TATB (1,3,5 trinitro, 2,4,6 triamino benzene), FOX-7 (1,1-diamino-2,2-dinitroethylene), LLM-105 (2,6-diamino-3,5-dinitropyrazine-1-oxide), and NTO (nitrotriazolone), all with 20% Sylgard 182 binder. Time resolving emission from shocked explosives produces a unique fingerprint at each shock pressure, and the emission temporal profile can be used to understand each material's shock response and the underlying mechanisms within the framework of the well-known shock ignition and growth models. Ignition was characterized by the emission intensity during nanosecond hot spot ignition. The growth process, which occurs after the shock has unloaded, was characterized by the time derivative of microsecond emission intensity. Singular-value decomposition was used to determine a characteristic spectral radiance over the entire time range, and this radiance was fit to a graybody model to determine a characteristic temperature Tc. Even though these insensitive explosives have similar shock sensitivities, our method can clearly distinguish their shock response. FOX-7 and NTO were the most sensitive, with FOX-7 having the greatest shock ignition sensitivity, growth rate, and Tc. LLM-105 was much less sensitive than FOX-7 and NTO, and TATB consistently showed the lowest sensitivity by any measure.

Akhtar, Meysam (ORCID:0000000243897837)↗
DOE OSTI2023

Atomic-scale electronic inhomogeneity in single-layer iron chalcogenide alloys revealed by machine learning of STM/S data

Chemical pressure from the isovalent substitution of Se by a larger Te atom in the epitaxial film of iron chalcogenide FeSe can effectively tune its superconducting, topological, and magnetic properties. However, such substitution during epitaxial growth inherently leads to defects and structural inhomogeneity, making the determination of alloy composition and atomic sites for the substitutional Te atoms challenging. Here, we utilize machine learning to distinguish between Se and Te atoms in scanning tunneling microscopy images of single-layer FeSe1−xTex on SrTiO3(001) substrates. Defect locations are first identified by analyzing spatial-dependent dI/dV tunneling spectra using the K-means clustering method. After excluding the defect regions, the remaining dI/dV spectra are further analyzed using the singular value decomposition method to determine the Se/Te ratio. Our findings demonstrate an effective and reliable approach for determining alloy composition and atomic-scale electronic inhomogeneity in superconducting single-layer iron chalcogenide films.

Materials Science↗
DOE OSTI2025

Sensitivity of the magnetic axis to variations in the magnetic field

Lagrangian variational principles are used to locate the magnetic axis and to determine its sensitivity to variations in the magnetic field. Only the change in the magnetic field on axis is required to determine, to lowest order, the location of the new magnetic axis. Using a simple example, namely, the magnetic field produced by a set of current-carrying filaments (a.k.a. coils), the variation in the coil geometry that induces the largest change in the magnetic axis is determined using singular value decomposition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗