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 73 records · Page 4

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↗

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↗

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↗

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↗

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↗

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 ↗

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↗

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↗

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↗

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↗

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

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↗

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↗

Deciphering the small-angle scattering of polydisperse hard spheres using deep learning

We introduce a deep learning approach for analyzing the scattering function of the polydisperse hard sphere system. We use a variational autoencoder-based neural network to learn the bidirectional mapping between the scattering function and the system parameters, including the volume fraction and polydispersity. Such that the trained model serves both as a generator that produces a scattering function from the system parameters and an inferrer that extracts system parameters from the scattering function. We first generate a scattering dataset by carrying out molecular dynamics simulations of the polydisperse hard spheres modeled by the truncated-shifted Lennard-Jones model, then analyze the scattering function dataset using singular value decomposition to confirm the feasibility of dimensional compression. Then, we split the dataset into training and testing sets and train our neural network on the training set only. Our generator model produces a scattering function with significantly higher accuracy compared to the traditional Percus–Yevick approximation and β correction, and the inferrer model can extract the volume fraction and polydispersity with much higher accuracy than traditional model functions.

Ding, Lijie [ORNL] (ORCID:0000000227454606)↗

Chasing the multi-modal plasma response in MAST-U

Abstract Achieving edge localized modes (ELMs) suppression in spherical tokamaks by applying resonant magnetic perturbations (RMPs) has proven challenging. The poloidal spectrum of the applied RMP is a key parameter that has an impact on the capability to mitigate and eventually suppress ELMs. In this work the resistive magnetohydrodynamic code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) is used to evaluate the possibility of directly measuring the plasma response in MAST-U, and particularly its variation as function of the applied poloidal spectrum, in order to guide the experimental validation of the predicted best RMP configuration for ELM suppression. Toroidal mode number n = 2 RMP is considered to minimize the presence of sidebands, and to avoid the deleterious core coupling of n = 1. Singular Value Decomposition is used to highlight linearly independent structures in the simulated magnetic 3D fields and how those structures can be measured at the wall where the magnetic sensors are located. Alternative ways to measure the multimodal plasma response and how they can be used to infer the best RMP configuration to achieve ELM suppression are also presented, including the plasma displacement and the 3D footprints at the divertor plates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

21 cm Power Spectrum Analysis of North Celestial Pole Observations with the Tianlai Dish Pathfinder Array

The Tianlai Dish Pathfinder Array (TDPA) is a radio interferometer designed to test techniques for 21 cm intensity mapping in the post-reionization Universe as a means of measuring large-scale cosmic structure. Using nine nights of observations targeting the North Celestial Pole field, totaling approximately 107 hr of integration time, we analyze data in the frequency range 700–800 MHz (corresponding to redshift z ∼ 0.9). We do the data format conversion, radio frequency interference flagging, calibration, imaging and point source subtraction, and foreground removal via Singular Value Decomposition. The spherically averaged power spectrum Δ 2 (k) is obtained. Furthermore, this work successfully establishes and validates a comprehensive data analysis framework for the TDPA. We identify key improvements including sky model refinement, increased integration time, and pipeline optimization that will enable future detection of the 21 cm signal through auto-correlation and cross-correlation with optical galaxy surveys.

cosmology: large-scale structure of universe↗

Surrogate models for plasma displacement and current in 3D perturbed magnetohydrodynamic equilibria in tokamaks

Abstract A numerical database of over one thousand perturbed three-dimensional (3D) equilibria has been generated, constructed based on the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) computed plasma response to the externally applied 3D field sources in multiple tokamak devices. Perturbed 3D equilibria with the n = 1–4 ( n is the toroidal mode number) toroidal periodicity are computed. Surrogate models are created for the computed perturbed 3D equilibrium utilizing model order reduction (MOR) techniques. In particular, retaining the first few eigenstates from the singular value decomposition (SVD) of the data is found to produce reasonably accurate MOR-representations for the key perturbed quantities, such as the perturbed parallel plasma current density and the plasma radial displacement. SVD also helps to reveal the core versus edge plasma response to the applied 3D field. For the database covering the conventional aspect ratio devices, about 95% of data can be represented by the truncated SVD-series with inclusion of only the first five eigenstates, achieving a relative error (RE) below 20%. The MOR-data is further utilized to train neural networks (NNs) to enable fast reconstruction of perturbed 3D equilibria, based on the two-dimensional equilibrium input and the 3D source field. The best NN-training is achieved for the MOR-data obtained with a global SVD approach, where the full set of samples used for NN training and testing are stretched and form a large matrix which is then subject to SVD. The fully connected multi-layer perceptron, with one or two hidden layers, can be trained to predict the MOR-data with less than 10% RE. As a key insight, a better strategy is to train separate NNs for the plasma response fields with different toroidal mode numbers. It is also better to apply MOR and to subsequently train NNs separately for conventional and low aspect ratio devices, due to enhanced toroidal coupling of Fourier spectra in the plasma response in the latter case.

3D equilibrium↗

Investigating the role of edge neutrals in exciting tearing mode activity and achieving flat temperature profiles in LTX-β

Abstract We present observations, numerical simulations, and analysis from experiments in the Lithium Tokamak Experiment-Beta (LTX- β ) in which the electron temperature profile ( T e ( r )) shifts from flat to peaked and a tearing mode is also destabilized when the average density ( n e ave ) exceeds ∼10 19 m −3 . Flat T e ( r ) is obtained routinely in LTX- β , with a lithium coated, low-recycling first wall, once the external fueling is stopped and density decays [Boyle et al 2023 Nucl. Fusion 63 056020]. In the present experiment, flat T e profiles can be sustained while maintaining constant n e ave below a line averaged density threshold ( n e ave th ) of ∼10 19 m −3 . Above n e ave th , T e ( r ) shifts from flat to peaked and a tearing mode is destabilized. Due to low recycling, the achieved n e ave can be controlled precisely by external fueling and hence, a certain threshold of the edge neutral inventory from the external fueling is experimentally manifested through n e ave th . The goal of the present work is to investigate the role of edge neutrals in determining T e ( r ) and MHD stability in the unique low-recycling regime of LTX- β . Our hypothesis is that the peaking of T e ( r ) beyond n e ave th is due ultimately to the edge cooling by the cold neutrals beyond a critical fueling flux. At lower fueling flux, flat T e ( r ) results in broader pressure profile and lower resistivity, which in turn stabilizes the tearing mode. This hypothesis is supported by edge neutral density estimation by DEGAS 2 code. Mode analysis by singular value decomposition confirms the tearing mode structure to be m / n = 2/1 ( m and n being the poloidal and toroidal mode numbers). Linear tearing stability analysis with M3D-C1 predicts that plasmas with n e ave > 10 19 are highly susceptible to a n = 1 tearing mode. ORBIT simulations, however, confirmed that the tearing modes do not contribute to the loss of fast ions from neutral beam injection. This study shows for the first time that the neutral inventory at the edge could be one of the deciding factors for the achievability of the unique operation regime of flat T e ( r ) and the excitation of tearing activity that could be disruptive for the plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗