Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Irradiation Effect on Noble Metal Particles in Water Using in situ Liquid Cell STEM Observation

During geologic disposal of spent nuclear fuel (SNF) in an engineered nuclear waste repository, once all other barriers have degraded, oxidizing may occur at the solid-water interface owing to a self-generated radiolytic field. The repository design includes large quantities of iron (Fe), that is anticipated to corrode under an anoxic environment, and generate hydrogen (H 2 ) gas. This H 2 gas is thought to be able to suppress the dissolution of SNF through a catalytic reaction with noble metal particles (NMP) that are pre-existing in the SNF. This interaction leads to the decomposition of the major oxidant, hydrogen peroxide (H 2 O 2 ). In conclusion, these processes are described in the Fuel Matrix Degradation (FMD) model that is being used to predict SNF degradation rates. The NMP, therefore, plays an important role within the FMD model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

OR22-Neuromorphic Rad Detector-PD3Ra (Final Report)

In unattended monitoring scenarios, automated radiation detection algorithms must be able to detect low signal-to-noise ratio (SNR) anomalies in a potentially dynamic and noisy background and report these anomalies in a timely fashion. Dynamic and noisy backgrounds complicate the use of simple gross-counting algorithms because they can lead to either high false positive rates or low sensitivity. Algorithms that use the entire spectrum have been the most successful in this area; notable examples are the NSCRAD algorithm developed at Pacific Northwest National Laboratory and recently the nonnegative matrix factorization approach developed at Lawrence Berkeley National Laboratory (LBNL). These approaches use either spectral regions of interest or spectral decomposition to detect threat isotopes in the background.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Random projection using random quantum circuits

The random sampling task performed by Google's Sycamore processor gave us a glimpse of the “quantum supremacy era.” This has definitely shed some light on the power of random quantum circuits in this abstract task of sampling outputs from the (pseudo)random circuits. In this paper, we explore a practical near-term use of local random quantum circuits in dimensional reduction of large low-rank data sets. We make use of the well-studied dimensionality reduction technique called the random projection method. This method has been extensively used in various applications such as image processing, logistic regression, entropy computation of low-rank matrices, etc. We prove that the matrix representations of local random quantum circuits with sufficiently shorter depths [ ∼ O ( n ) ] serve as good candidates for random projection. We demonstrate numerically that their projection abilities are not far off from the computationally expensive classical principal components analysis on MNIST and CIFAR-100 image datasets. We also benchmark the performance of quantum random projection against the commonly used classical random projection in the tasks of dimensionality reduction of image data sets and computing von Neumann entropies of large low-rank density matrices. And finally, using variational quantum singular value decomposition, we demonstrate a near-term implementation of extracting the singular vectors with dominant singular values after quantum random projecting a large low-rank matrix to lower dimensions. All such numerical experiments unequivocally demonstrate the ability of local random circuits to randomize a large Hilbert space at sufficiently shorter depths with robust retention of properties of large data sets in reduced dimensions. Published by the American Physical Society 2024

Kumaran, Keerthi (ORCID:0009000949125721)↗
DOE OSTI2020

Localized Exponential Time Differencing Method for Shallow Water Equations: Algorithms and Numerical Study

Here, we explore the performance of the exponential time differencing (ETD) method applied to the rotating shallow water equations. Comparing with explicit time stepping of the same order accuracy in time, the ETD algorithms could reduce the computational time in many cases by allowing the use of large time step sizes while still maintaining numerical stability. To accelerate the ETD simulations, we propose a localized approach that synthesizes the ETD method and overlapping domain decomposition. By dividing the original problem into many subdomain problems of smaller sizes and solving them locally, the proposed approach could speed up the calculation of matrix exponential vector products. Several standard test cases for shallow water equations of one or multiple layers are considered. The results show great potential of the localized ETD method for high-performance computing because each subdomain problem can be naturally solved in parallel at every time step.

58 GEOSCIENCES↗
DOE OSTI2023

Development of Hydrogen Permeation Barrier and Matrix Interaction Protection Coatings for TZM Enclosure of Hydride Moderators (Microreactor Program M3 Report)

It is important to develop a high-performance enclosure solution for hydride-based microreactor moderation components to effectively suppress hydride thermal decomposition and prevent hydrogen loss in harsh in-reactor environments. Currently, the DOE-NE Microreactor Program (MR program) is investigating the titanium-zirconium-molybdenum (TZM) canned yttrium hydride (YH 2-x ) solution, which provides acceptable hydrogen retainment performance up to ~700 °C. However, improvements in the enclosure are necessary for increasing long-term operation temperature beyond 700 °C and enhancing short-term survivability at even higher temperatures during power transients. One potential solution is the addition of a dedicated hydrogen permeation barrier, which can potentially improve the enclosure's performance. Furthermore, the hydride moderator components will be inserted into the microreactor matrix, which may consist of materials that are not chemically inert against TZM. Therefore, the TZM might need to be protected from the external environment through the application of another surface protection coating. To address these critical issues, two types of barrier coatings optimized for TZM-based moderator enclosure has been developed. A multilayer coating approach has been used to identify TZM-compatible barrier materials. Behavior of these coating designs against H 2 permeation and long-term high-temperature graphite matrix interactions has been examined in detail. Additionally, the developed coatings' behavior when exposed to both high-temperature thermal cycling and heavy-ion irradiation has also been investigated. Overall, the development of these barrier coatings is essential for achieving maturity for the metal hydride moderator technology. The results in this report will provide crucial insights into the behavior of the developed coatings and their potential for improving the long-term performance of TZM-based moderator enclosures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Foundations of variational discrete action theory

Variational wave functions and Green's functions are two important paradigms for solving quantum Hamiltonians, each having their own advantages. Here we detail the variational discrete action theory (VDAT), which exploits the advantages of both paradigms in order to approximately solve the ground state of quantum Hamiltonians. VDAT consists of two central components: the sequential product density matrix (SPD) ansatz and a discrete action associated with the SPD. The SPD is a variational ansatz inspired by the Trotter decomposition and characterized by an integer $\mathscr{N}$, recovering many well-known variational wave functions, in addition to the exact solution for $\mathscr{N}$ = ∞. The discrete action describes all dynamical information of an effective integer time evolution with respect to the SPD. We generalize the path integral to our integer time formalism, which converts a dynamic correlation function in integer time to a static correlation function in a compound space. We also generalize the usual many-body Green's function formalism to integer time, which results in analogous but distinct mathematical structures, yielding integer time versions of the generating functional, Dyson equation, and Bethe-Salpeter equation. We prove that the SPD can be exactly evaluated in the multiband Anderson impurity model (AIM) by summing a finite number of diagrams. For the multiband Hubbard model, we prove that the self-consistent canonical discrete action approximation (SCDA), which is the integer time analog of the dynamical mean-field theory, exactly evaluates the SPD for d = ∞. VDAT within the SCDA provides an efficient yet reliable method for capturing the local physics of quantum lattice models, which will have broad applications for strongly correlated electron materials. More generally, VDAT should find applications in various many-body problems in physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2009

Origins of Magnetite Nanocrystals in Martian Meteorite ALH84001

The Martian meteorite ALH84001 preserves evidence of interaction with aqueous fluids while on Mars in the form of microscopic carbonate disks. These carbonate disks are believed to have precipitated 3.9 Ga ago at beginning of the Noachian epoch on Mars during which both the oldest extant Martian surfaces were formed, and perhaps the earliest global oceans. Intimately associated within and throughout these carbonate disks are nanocrystal magnetites (Fe3O4) with unusual chemical and physical properties, whose origins have become the source of considerable debate. One group of hypotheses argues that these magnetites are the product of partial thermal decomposition of the host carbonate. Alternatively, the origins of mag- netite and carbonate may be unrelated; that is, from the perspective of the carbonate the magnetite is allochthonous. For example, the magnetites might have already been present in the aqueous fluids from which the carbonates were believed to have been deposited. We have sought to resolve between these hypotheses through the detailed characterized of the compo- sitional and structural relationships of the carbonate disks and associated magnetites with the orthopyroxene matrix in which they are embedded. Extensive use of focused ion beam milling techniques has been utilized for sample preparation. We then compared our observations with those from experimental thermal decomposition studies of sideritic carbonates under a range of plausible geological heating scenarios. We conclude that the vast majority of the nanocrystal magnetites present in the car- bonate disks could not have formed by any of the currently proposed thermal decomposition scenarios. Instead, we find there is considerable evidence in support of an alternative allochthonous origin for the magnetite unrelated to any shock or thermal processing of the carbonates.

Thomas-Keprta, Kathie L.↗
NASA NTRS1993

Application of the target decomposition theorem to a polarimetric random media model

With advances in polarimetric radar measurements of land surfaces, the need for understanding the underlying scattering mechanisms and dominant target features has become the focus of many studies. In particular, the maximum use of the polarimetric information to identify and/or separate parameters related to the surface features such as vegetation thickness, structure, water content, and soil surface characteristics will enhance the possibility of using polarimetric radars for monitoring the earth's surface from space. In this paper, Cloude's decomposition theorem is applied to a polarimetric random media model to simulate the radar measurements of vegetated canopies. The vegetated canopies are modeled as a three layer discrete random medium with leaves and branches in the first layer, tree trunks in the second layer, and a half space of homogeneous ground with rough interface as the bottom layer. The distorted born approximation (DBA) has been used to compute full Mueller matrix of the canopy, using canonical dielectric objects such as thin discs and cylinders as leaves, branches, and trunks, respectively. The Mueller matrix and the derived covariance matrix contain information on the second order statistics of radar signals at various polarizations from the canopy. To decompose the covariance matrix to its constituent targets, the eigenvalues and eigenvectors of the covariance matrix are computed in terms of the physical parameters of the canopy. In addition, each eigenvector explicitly shows the scattering mechanisms such as odd and even reflections in the canopy. Cloude's decomposition theorem is applied using the Pauli spin matrices as a basis and an expression for the degree of disorder or the entropy for the vegetated surface is found. Then, the physical parameters estimated from in situ measurements are used in the random media to obtain realistic covariance matrices. As a result, the sensitivity of the eigenvalue spectrums and the coefficients resulting from the target decomposition theorem to the physical parameters of the canopy are examined and the possible use of Cloude's theorem to estimate vegetation and soil parameters is discussed.

Saatchi, Sasan S.↗
NASA NTRS1983

Asymptotic unbounded root loci - Formulas and computation

A new geometric way of computing the asymptotic behavior of unbounded root loci of a strictly proper linear time-invariant control system as loop gain goes to infinity is presented. Properties of certain restricted linear maps and nested restrictions of linear maps are developed, and formulas are obtained for the leading coefficient of the asymptotic values of the unbounded multivariable root loci are obtained in terms of eigenvalues of those maps. Published results and a certain simple null structure assumption are used to relate these asymptotic values to the structure at infinity of the Smith-McMillan form of the open loop transfer function. Explicit matrix formulas for the more abstract derived formulas are given and additional geometric insights are developed with orthogonal projections and singular value decomposition. Formulas for the pivots of the unbounded root loci are calculated and shown to have the same form as the coefficients of the unbounded asymptotic root loci.

Sastry, S. S.↗
NASA NTRS1986

Use of laboratory simulated pyrolysis in tracing the history of sedimentary organic matter

Results from laboratory simulated pyrolyses experiments show that in addition to depth of burial, preservation of kerogen, and hence any morphologic structure in it, is also dependent on the mineral matrix with which it is associated. In the presence of clay minerals, and especially under dry conditions, extractable lipids released during kerogen decomposition are more rapidly destroyed than in the presence of calcite or chert matrices. The result is production of gas, polar bitumen and pyrobitumen and destruction of biomarkers. During such an early reorganization of the kerogen, the biomarker constituents can be destroyed, or unrecognizably altered. The above process of organic residues maturation appears to be inhibited in the presence of water and is significantly reduced where kerogen is hosted in limestones, dolomites or cherts. These minerals have been characteristically found to be the most reliable in yielding morphological fossils and small quantities of extractable bitumen in Archean and Proterozoic rocks. To understand the validity of chemical and morphological fossils, in the early geologic record, it will be necessary to understand the process of kerogen in sedimentary rocks. To test the role of various minerals on the preservation process, kerogen extracted from a variety of rocks has been heated together with montmorillonite, illite and calcite. The kinetics of the process has been monitored and the products quantitatively identified.

Kaplan, I. R.↗
NASA NTRS1987

Robust eigensystem assignment for flexible structures

An improved method is developed for eigenvalues and eigenvectors placement of a closed-loop control system using either state or output feedback. The method basically consists of three steps. First, the singular value of QR decomposition is used to generate an orthonormal basis that spans admissible eigenvector space corresponding to each assigned eigenvalue. Secondly, given a unitary matrix, the eigenvector set which best approximates the given matrix in the least-square sense and still satisfy eigenvalue cosntraints is determined. Thirdly, a unitary matrix is sought to minimize the error between the unitary matrix and the assignable eigenvector matrix. For use as the desired eigenvector set, two matrices, namely, the open-loop eigenvector matrix and its closest unitary matrix are proposed. The latter matrix generally encourages both minimum conditioning and control gains. In addition, the algorithm is formulated in real arithmetic for efficient implementation. To illustrate the basic concepts, numerical examples are included.

Juang, Jer-Nan↗
NASA NTRS1989

An interpretation and solution of ill-conditioned linear equations

Data insufficiency, poorly conditioned matrices and singularities in equations occur regularly in complex optimization, correlation, and interdisciplinary model studies. This work concerns itself with two methods of obtaining certain physically realistic solutions to ill-conditioned or singular algebraic systems of linear equations arising from such studies. Two efficient computational solution procedures that generally lead to locally unique solutions are presented when there is insufficient data to completely define the model, or a least-squares error formulation of this system results in an ill-conditioned system of equations. If it is assumed that a reasonable estimate of the uncertain data is available in both cases cited above, then we shall show how to obtain realistic solutions efficiently, in spite of the insufficiency of independent data. The proposed methods of solution are more efficient than singular-value decomposition for dealing with such systems, since they do not require solutions for all the non-zero eigenvalues of the coefficient matrix.

Ojalvo, I. U.↗
NASA NTRS2012

The Role of (Delta)C-13 in the Search for Reduced Organics on the Surface of Mars

The capabilities of the Sample Analysis at Mars (SAM) instrument suite on the Mars Science Laboratory (MSL) to detect trace amounts of organic carbon compounds are unprecedented, and MSL may be the first mission to reveal the presence of organic carbon on Mars. The search for reduced organic carbon on Mars is inextricably tied to: a) the preservation potential of the environment from which we take a solid sample, and b) the evolved gas analysis (EGA) techniques used by SAM to release volatiles from this solid sample. Several prospective targets have been identified for sample analysis at Gale Crater. Stratigraphic sequences of phyllosilicates and sulfates at Gale are thought to represent a period of global climate transition from a moderate pH lacustrine environment to an evaporitic environment, both of which could sequester organic carbon (Thomson et al. 2011). The sediment mound in Gale Crater contains a range of lithologies suggesting changes in redox conditions, and evidence of both lacustrine and fluvial depositional processes, which may have transported organic carbon from the layer in which it formed and resulted in its preservation elsewhere within the crater (Anderson and Bell, 2010). Inverted channel fills suggest erosion resistant material that could serve to preserve organics originally deposited in a low energy aqueous environment. The lithology sampled will affect not only the preservation of organics, but also our ability to detect organics during our evolved gas analysis, based on the sample matrix. For example, reduced organics may be trapped in the mineral structure, and thermal evolution of these organics will occur during thermal decomposition of the host mineral. If organics are occluded in minerals that have very high thermal decomposition temperatures, they may be, in effect, "too well preserved," and difficult to detect during EGA. Alternatively, the possible presence of perchlorate, or other strong oxidants in surface regolith, may result in destruction of structural information identifying organic molecules before reaching the QMS on SAM via oxidation to C02 during heating. If this is the case, the stable carbon isotopic composition (delta 13C) of the C02 evolved and measured by the Tunable Laser Spectrometer (TLS) on SAM may help identify the presence of organics. On Earth, biological activity can cause large fractionations of 13C/12C, which can preserved in sedimentary deposits and distinguish the organic products of biotic processes from inorganic atmospheric and geological reservoirs. It is plausible that similar fractionations could occur on Mars and be preserved in reduced organic matter in sediments. Bulk delta 13C measurements alone may not reveal a signature of trace organic carbon that may be present along with inorganic carbon. If both organic and inorganic carbon compounds are present, it may be possible to detect the organic carbon by comparing the 013C of pyrolysis and combustion experiments. The TLS on SAM is capable of obtaining high precision measurements of delta13C from C02 evolved during pyrolysis and combustion of solid regolith samples. Because carbonates are expected to be present at abundances of 0.1-1 % in Martian soil, and organics in the ppb range (Webster and Mahaffy, 2011), analog samples must represent this mix of reduced organic carbon and carbonate. The work presented here will examine the use of delta13C of C02 produced during combustion of bulk Mars analog samples as a proxy for detection of reduced organic carbon.

Stern, J. C.↗
DOE OSTI2026

Definitive Assessment of the Accuracy, Variationality, and Convergence of Relativistic Coupled Cluster and Density Matrix Renormalization Group in 100-Orbital Space

Accuracy, variationality, and convergence underpin the reliability of modern electronic structure methods, yet definitive benchmarks in the relativistic regime remain elusive due to the absence of numerically exact full configuration interaction (CI) references. Recent algorithmic advances in the CI framework, enabled by the small-tensor-product (STP) decomposition approach, have dramatically extended the tractable size of the configuration space, making numerically exact CI calculations feasible in large active spaces previously beyond reach. In this paper, we employ the recently developed STP-CI framework to perform large-scale numerically exact CI calculations and directly benchmark relativistic coupled cluster and density matrix renormalization group methods. Definitive benchmarking of approximate relativistic electronic structure methods is ensured through the application of the gap theorem, which provides rigorous error bounds on the CI reference and establishes a controlled standard for assessing accuracy, variationality, and convergence.

Chemical calculations↗
NASA NTRS2010

Probabilistic Round Trip Contamination Analysis of a Mars Sample Acquisition and Handling Process Using Markovian Decompositions

A method for evaluating the probability of a Viable Earth Microorganism (VEM) contaminating a sample during the sample acquisition and handling (SAH) process of a potential future Mars Sample Return mission is developed. A scenario where multiple core samples would be acquired using a rotary percussive coring tool, deployed from an arm on a MER class rover is analyzed. The analysis is conducted in a structured way by decomposing sample acquisition and handling process into a series of discrete time steps, and breaking the physical system into a set of relevant components. At each discrete time step, two key functions are defined: The probability of a VEM being released from each component, and the transport matrix, which represents the probability of VEM transport from one component to another. By defining the expected the number of VEMs on each component at the start of the sampling process, these decompositions allow the expected number of VEMs on each component at each sampling step to be represented as a Markov chain. This formalism provides a rigorous mathematical framework in which to analyze the probability of a VEM entering the sample chain, as well as making the analysis tractable by breaking the process down into small analyzable steps.

Hudson, Nicolas↗
DOE OSTI2026

Insights into Native Single-Atom Electrocatalyst Site Structures

Single-atom electrocatalysts consisting of metal atoms embedded in a carbon matrix are promising next-generation catalysts for green hydrogen production and utilization, CO2 reduction, low-temperature CO oxidation, ammonia production, plastic decomposition, and electrochemical energy storage. The origins of activity and stability for the single-atom sites are still debatable, however, because of constrained insights into their local structure resulting from idealized models and experiments derived from a large number of individual sites. Insights into structural variations around single atomic sites are therefore critical for the continued development of these next-generation catalysts. While electron microscopy commonly provides atomic-scale information about these materials, the beam sensitivity of individual sites makes structural determination by conventional low-voltage (60 keV) techniques challenging. Here, we introduce ultralow-voltage electron ptychography, performed at 30 keV, that enables determination of the lattice structure around individual metal sites in a well-defined single-atom electrocatalyst system while essentially eliminating knock-on structural modifications. Pairing these atomic-scale, site-specific measurements with computational methods will broaden our understanding of the activity and stability of these materials, which will accelerate the development of the next generation of catalysts.

Zachman, Michael [ORNL] (ORCID:0000000319101357)↗
NASA NTRS1988

Eigenstructure assignment for the control of highly augmented aircraft

Eigenstructure assignment is utilized to design flight control laws for aircraft with many control effectors. It is shown that a previous eigenstructure design for the flight propulsion control coupling (FPCC) aircraft lateral dynamics with three control surfaces exhibits a lack of stability robustness because the control distribution matrix is nearly rank-deficient. A method is mapped back, reducing the control space to two dimensions by using the singular value decomposition. After the design is complete, the controller is mapped back to the original three-dimensional control space. This design approach yields a controller with both smaller gains and improved multivariable stability margins at the aircraft inputs. An interesting characteristics of the control mapping, as applied to the given example, is that the most effective inputs have the larger gains while the less effective inputs have smaller gains.

Sobel, Kenneth M.↗
NASA NTRS1993

Preconditioned domain decomposition scheme for three-dimensional aerodynamic sensitivity analysis

A preconditioned domain decomposition scheme is introduced for the solution of the 3D aerodynamic sensitivity equation. This scheme uses the iterative GMRES procedure to solve the effective sensitivity equation of the boundary-interface cells in the sensitivity analysis domain-decomposition scheme. Excluding the dense matrices and the effect of cross terms between boundary-interfaces is found to produce an efficient preconditioning matrix.

Eleshaky, Mohammed E.↗