SEARCH · Engineering Papers

Results for “Matrix factorization”

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 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Tuning the Properties of Nanocomposites by Trapping Them in Deep Metastable States

We report the interfacial region formed in a polymer matrix around nanoparticles (NPs) controls many macroscopic properties of polymer nanocomposites (PNCs). However, understanding the factors controlling structure and properties of interfacial region remains a challenge. We demonstrated that the initial trapping of polymer chains at surfaces of NPs in solution strongly affects the macroscopic properties of PNCs. The most surprising result is that the differences in properties of PNCs persist even after an extremely long thermal annealing time. We ascribe the observed changes to the formation of interfacial layer that is trapped in a deep metastable state already in solution. Furthermore, presented analysis suggests that PNC equilibration time is defined by the chain desorption time that can be extremely long and, in some cases, even not accessible on a reasonable experimental time scale. These results highlight the importance of polymer solution concentrations on the formation of interfacial layer and the macroscopic properties of PNCs.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Discontinuous Galerkin discretization for quantum simulation of chemistry

All-electron electronic structure methods based on the linear combination of atomic orbitals method with Gaussian basis set discretization offer a well established, compact representation that forms much of the foundation of modern correlated quantum chemistry calculations—on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block-diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from O(N 4.5 ) in molecular orbital representations to O(N 2.6 ) in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1–2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Double copy root of Hawking thermality

The Hawking radiation spectrum from a collapsing null shell can be derived via the double copy of a simpler gauge theory calculation. Analyzing the non-abelian Yang-Mills root of this process, we demonstrate that the radiation spectrum is thermal in the color charge eigenvalue 𝜆, not energy. Considering the 𝑆⁢𝑈⁡(𝑁 𝑐 ) gauge theory in the large 𝑁 𝑐 limit, we find the differential spectrum ⅆ⁢𝑁/ⅆ⁢𝜆 is a product of the gravitationally familiar Planck-like factor and the color phase space density, modeled here as the Wigner semicircle from random matrix theory. This reveals that apparent energy thermality in gravity is the direct dual of charge thermality in its underlying non-abelian gauge theory.

Carrasco, John Joseph M. [Northwestern University,↗

DOE OSTI2025

Gauge-fixing quantum density operators at scale

We provide a theory, algorithms, and simulations of nonequilibrium quantum systems using a one-dimensional (1D) completely positive (CP), matrix-product (MP) density-operator (𝜌) representation. By generalizing the matrix product state's orthogonality center, to additionally store positive classical mixture correlations, the MP⁢𝜌 factorization naturally emerges. In this setting, we analytically and numerically examine the virtual gauge freedoms associated with the representation of quantum density operators. Based on this perspective, we simplify algorithms in certain limits to speed up the integration of the canonical-form master-equation dynamics. This enables us to quickly evolve under the dynamics of two-body quantum channels without resorting to optimization-based methods. In addition to this technical advance, we also scale up numerical examples and discuss implications for accurately modeling hardware architectures and predicting their performance in the near term. This includes an example of the quantum to classical transition of informationally leaky, i.e., decohering, qubits. In this setting, because of loss from environmental interactions, nonlocal complex coherence correlations are converted into global incoherent classical statistical mixture correlations. Lastly, the representation of both global and local correlations is discussed. We expect this work to have applications in additional nonequilibrium settings, beyond qubit engineering.

Gangapuram, Amit Jamadagni [Oak Ridge National Lab↗

DOE OSTI2021

Neutron-neutron scattering length from the 6 He(p,pα)nn reaction

We propose a novel method to measure the neutron-neutron scattering length using the 6 He(p,pα)nn reaction in inverse kinematics at high energies. The method is based on the final-state interaction (FSI) between the neutrons after the sudden knockout of the α particle. We show that the details of the neutron-neutron relative-energy distribution allow for a precise extraction of the s-wave scattering length. Here, we present the state of the art in regard to the theory of this distribution. The distribution is calculated in two steps. First, we calculate the ground-state wave function of 6 He as a αnn three-body system. For this purpose we use Halo effective field theory, which also provides uncertainty estimates for the results. We compare our results at this stage to model calculations done with the computer code face. In a second step we determine the effects of the nn FSI using the nn t-matrix. We compare these FSI results to approximate FSI approaches based on standard FSI enhancement factors. While the final distribution is sensitive to the nn scattering length, it depends only weakly on the effective range. Throughout we emphasize the impact of theoretical uncertainties on the neutron-neutron relative-energy distribution, and discuss the extent to which those uncertainties limit the extraction of the neutron-neutron scattering length from the reaction 6 He(p,pα)nn.

6 ≤ A ≤ 19↗

DOE OSTI2023

FWP FEAA149: “Next Generation Environmental Barrier Coatings”

Environmental barrier coatings (EBCs) are required coatings for utilization of SiC/SiC ceramic matrix composite (CMC) components in gas turbines, where the EBC represents the life-limiting factor for such components. EBC/CMC systems have shown success in aero-engine applications with increased turbine inlet temperatures and improved efficiencies, which are achieved through higher temperature stability, lower density, and decreased reliance on cooling air compared to traditional superalloys. While industrial gas turbines (IGTs) do not currently utilize SiC/SiC CMCs, the current shift towards low-carbon or carbon-free fuel sources for power generation could result in a need for EBC/CMC components with higher temperature capabilities. In this work, three tasks were outlined to improve understanding of EBC lifetimes to encourage use in IGTs with carbon-free fuel such as hydrogen: 1. Define the bond coating oxidation kinetics and EBC failure criteria, 2. Measure thermal expansion coefficients of each layered material, and 3. Perform advanced characterization and modeling to assess EBC lifetimes. Cyclic steam oxidation tests were conducted on various EBC/Si/SiC chemistries and EBC/SiC architectures to define substrate oxidation kinetics and EBC failure modes. An open-source code was developed to quantify the undulating thermally grown oxide thickness with thousands of measurements from specimen cross-section images. Bond coating oxidation kinetics were determined and used to develop a kinetic and thermodynamic model for predicting EBC lifetimes. High-temperature Raman spectroscopy was utilized for determining the SiO 2 thermally grown oxide phase transformation as the life-limiting feature for EBCs. Model efforts supported the claim that the SiO 2 phase transformation causes elevated stress during thermal cycling with associated cracking that decreases the adhesion strength of the EBC, eventually resulting in coating spallation. The finite element model subroutine will be made publicly available upon internal review. Further development of an EBC lifetime model for IGTs involves definition of a critical SiO 2 thickness for EBC spallation and must also consider both environmental (gas velocity, pressure, etc.) and specimen (EBC dopants, layer architectures, etc.) effects into predicted bond coating oxidation kinetics for long-lifetime components.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A Challenge for Discrimination of Color-Singlet versus Color-Octet Quarkonium Production

The precise mechanism for production of quarkonium at hadron colliders is still an open question. Within non-relativistic quantum chromodynamics, quarkonium production cross sections can be factorized into short-distance, perturbative contributions and universal, non-perturbative, long-distance matrix elements, and then summed over quantum numbers of the heavy quark pair. In principle, at short-distances, the heavy quark pair can be either in a color-singlet or color-octet state, and it is desirable to establish the relative contributions to compare with data and to make predictions in different experimental environments. From the explicit form of the lowest-order perturbative matrix elements for color-singlet and color-octet production, we show that the structure of the optimal observable for discrimination on phase space, the likelihood ratio, has strong dependence on the angular momentum state of the heavy quark pair. This presents an obstruction for construction of a single, robust discrimination observable that can be applied to production of an arbitrary quarkonium state.

Larkoski, Andrew↗

NASA NTRS1973

A penny shaped crack in a filament-reinforced matrix. 2: The crack problem

The elastostatic interaction problem between a penny-shaped crack and a slender inclusion or filament in an elastic matrix was formulated. For a single filament as well as multiple identical filaments located symmetrically around the crack the problem is shown to reduce to a singular integral equation. The solution of the problem is obtained for various geometries and filament-to-matrix stiffness ratios, and the results relating to the angular variation of the stress intensity factor and the maximum filament stress are presented.

Pacella, A. H.↗

NASA NTRS1982

Sampling errors in the estimation of empirical orthogonal functions

Empirical Orthogonal Functions (EOF's), eigenvectors of the spatial cross-covariance matrix of a meteorological field, are reviewed with special attention given to the necessary weighting factors for gridded data and the sampling errors incurred when too small a sample is available. The geographical shape of an EOF shows large intersample variability when its associated eigenvalue is 'close' to a neighboring one. A rule of thumb indicating when an EOF is likely to be subject to large sampling fluctuations is presented. An explicit example, based on the statistics of the 500 mb geopotential height field, displays large intersample variability in the EOF's for sample sizes of a few hundred independent realizations, a size seldom exceeded by meteorological data sets.

North, G. R.↗

NASA NTRS1985

Interaction between a crack and a soft inclusion

With the application to weld defects in mind, the interaction problem between a planar-crack and a flat inclusion in an elastic solid is considered. The elastic inclusion is assumed to be sufficiently thin so that the thickness distribution of the stresses in the inclusion may be neglected. The problem is reduced to a system of four integral equations having Cauchy-type dominant kernels. The stress intensity factors are calculated and tabulated for various crack-inclusion geometries and the inclusion to matrix modulus ratios, and for general homogeneous loadiong conditions away from the crack-inclusion region.

Xue-Hui, L.↗

NASA NTRS1987

Chemistry-split techniques for viscous reactive blunt body flow computations

The weak-coupling structure between the fluid and species equations has been exploited and resulted in three, closely related, time-iterative implicit techniques. While the primitive variables are solved in two separated groups and each by an Alternating Direction Implicit (ADI) factorization scheme, the rate-species Jacobian can be treated in either full or diagonal matrix form, or simply ignored. The latter two versions render the split technique to solving for species as scalar rather than vector variables. The solution is completed at the end of each iteration after determining temperature and pressure from the flow density, energy and species concentrations. Numerical experimentation has shown that the split scalar technique, using partial rate Jacobian, yields the best overall stability and consistency. Satisfactory viscous solutions were obtained for an ellipsoidal body of axis ratio 3:1 at Mach 35 and an angle of attack of 20 degrees.

Li, C. P.↗

NASA NTRS1989

Global properties of pseudospectral methods

The present application of polynomial interpolation methods to function-approximation and numerical solutions for hyperbolic and elliptic PDEs allows the explicit construction of the derivative matrix for a general sequence of collocation points. An evaluation of the effect of several factors on the performance of these methods indicates an inability to interpret global methods in terms of local ones; the accuracy of the approximation will differ when the function's large gradients occur near the center of the region or near the boundary, irrespective of the boundary vicinity's collocation-point density.

Solomonoff, A.↗

NASA NTRS1993

Optimum Interfacial Layers For Fiber/Metal Composites

Report presents theoretical stress-analysis study of some of factors affecting choice of thin layer of material placed at interface between each fiber and matrix of ceramic-fiber/metal-matrix composite. Effects of thickness, modulus of elasticity, and coefficient of thermal expansion considered.

Lerch, Bradley A.↗

NASA NTRS1998

Constitutive Modeling and Testing of Polymer Matrix Composites Incorporating Physical Aging at Elevated Temperatures

Advanced polymer matrix composites (PMC's) are desirable for structural materials in diverse applications such as aircraft, civil infrastructure and biomedical implants because of their improved strength-to-weight and stiffness-to-weight ratios. For example, the next generation military and commercial aircraft requires applications for high strength, low weight structural components subjected to elevated temperatures. A possible disadvantage of polymer-based composites is that the physical and mechanical properties of the matrix often change significantly over time due to the exposure of elevated temperatures and environmental factors. For design, long term exposure (i.e. aging) of PMC's must be accounted for through constitutive models in order to accurately assess the effects of aging on performance, crack initiation and remaining life. One particular aspect of this aging process, physical aging, is considered in this research.

Veazie, David R.↗

NASA NTRS2003

X-ray Polarimetry with an Active-Matrix Pixel Proportional Counter

We report the first results from an X-ray polarimeter with a micropattern gas proportional counter using an amorphous silicon active matrix readout. With 100% polarized X-rays at 4.5 keV, we obtain a modulation factor of 0.33+/- 0.03, confirming previous reports of the high polarization sensitivity of a finely segmented pixel proportional counter. The detector described here has a geometry suitable for the focal plane of an astronomical X-ray telescope. Amorphous silicon readout technology will enable additional extensions and improvements.

Black, J. K.↗

DOE OSTI2022

Introduction of the Adding and Doubling Method for Solving Bateman Equations for Nuclear Fuel Depletion

This paper introduces and evaluates the Adding and Doubling Method (ADM) for solving the Bateman equations for depletion systems with varying numbers of nuclides and compares it to the Chebyshev Rational Approximation Method (CRAM), both implemented in the reactor physics analysis application Griffin. ADM, when applied to the Crank-Nicolson Finite Difference method, can produce results comparable in accuracy and precision to CRAM with comparable run times for systems with 35 or 297 nuclides. For systems with more than 300 nuclides, the matrix-matrix operations required by ADM are significantly more costly than the matrix-vector operations required by CRAM, making CRAM the more efficient method for systems with large numbers of nuclides. ADM is an accurate method that maintains other advantages over CRAM in that it does not depend on pre-generated coefficients or require complex number operations. ADM also manages to outperform CRAM by a factor of more than 250 in terms of run time for depletion systems that require multiple Bateman solves while the depletion matrix and time step size remain constant over all depletion intervals.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1985

Adaptive Inverse Control for Rotorcraft Vibration Reduction

This thesis extends the Least Mean Square (LMS) algorithm to solve the mult!ple-input, multiple-output problem of alleviating N/Rev (revolutions per minute by number of blades) helicopter fuselage vibration by means of adaptive inverse control. A frequency domain locally linear model is used to represent the transfer matrix relating the higher harmonic pitch control inputs to the harmonic vibration outputs to be controlled. By using the inverse matrix as the controller gain matrix, an adaptive inverse regulator is formed to alleviate the N/Rev vibration. The stability and rate of convergence properties of the extended LMS algorithm are discussed. It is shown that the stability ranges for the elements of the stability gain matrix are directly related to the eigenvalues of the vibration signal information matrix for the learning phase, but not for the control phase. The overall conclusion is that the LMS adaptive inverse control method can form a robust vibration control system, but will require some tuning of the input sensor gains, the stability gain matrix, and the amount of control relaxation to be used. The learning curve of the controller during the learning phase is shown to be quantitatively close to that predicted by averaging the learning curves of the normal modes. For higher order transfer matrices, a rough estimate of the inverse is needed to start the algorithm efficiently. The simulation results indicate that the factor which most influences LMS adaptive inverse control is the product of the control relaxation and the the stability gain matrix. A small stability gain matrix makes the controller less sensitive to relaxation selection, and permits faster and more stable vibration reduction, than by choosing the stability gain matrix large and the control relaxation term small. It is shown that the best selections of the stability gain matrix elements and the amount of control relaxation is basically a compromise between slow, stable convergence and fast convergence with increased possibility of unstable identification. In the simulation studies, the LMS adaptive inverse control algorithm is shown to be capable of adapting the inverse (controller) matrix to track changes in the flight conditions. The algorithm converges quickly for moderate disturbances, while taking longer for larger disturbances. Perfect knowledge of the inverse matrix is not required for good control of the N/Rev vibration. However it is shown that measurement noise will prevent the LMS adaptive inverse control technique from controlling the vibration, unless the signal averaging method presented is incorporated into the algorithm.

Jacklin, Stephen A.↗

NASA NTRS1999

Gravity and Skeletal Growth

Two simultaneous experiments were performed using 5-week-old male Sprague Dawley rats; in one study, the rats were flown in low earth orbit; in the other study, the hindlimbs of the growing rats were elevated to prevent weight bearing. Following 9 d of unloading, weight bearing was restored for 4, 28, and 76 hrs. Afterwards, additional hindlimb unloading experiments were performed to evaluate the skeletal response to 0, 2, 4, 6, 8, 10, 12, 16, and 24 hrs of restored weight bearing following 7 d of unloading. Cancellous and cortical bone histomorphometry were evaluated in the left tibia at the proximal metaphysis and in the left femur at mid-diaphysis, respectively. Steady-state mRNA levels for bone matrix proteins and skeletal signaling peptides were determined in total cellular RNA extracted from trabeculae from the right proximal tibiametaphysis and periosteum from the right femur. Spaceflight and hindlimb unloading each resulted in cancellous osteopenia, as well as a tendency towards decreased periosteal bone formation. Both models for skeletal unloading resulted in site specific reductions in mRNA levels for transforming growth factor-beta (sub 1) (TGF-beta) osteocalcin (OC), and prepro-alpha (I) subunit of type 1 collagen (collagen) and little or no changes in mRNA levels for glyceraldehyde-3-phosphate dehydrogenase (GAP) and insulin-like growth factor I (IGF-I). Restoration of normal weight bearing resulted in transient increases in mRNA levels for the bone matrix proteins and TGF-beta in the proximal metaphysis and periosteum and no changes in either GAP or IGF-I mRNA levels. The timecourse for the response differed between the two skeletal compartments; the tibial metaphysis responded much more quickly to reloading. These results suggest that the skeletal adaptation to acute physiological changes in mechanical usage are mediated, in part, by changes in mRNA levels for bone matrix proteins and TGF-beta.

Morey-Holton, Emily↗