Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A ceramic matrix composite thermal protection system for hypersonic vehicles

The next generation of hypersonic vehicles (NASP, SSTO) that require reusable thermal protection systems will experience acreage surface temperatures in excess of 1100 C. More important, they will experience a more severe physical environment than the Space Shuttle due to non-pristine launching and landing conditions. As a result, maintenance, inspection, and replacement factors must be more thoroughly incorporated into the design of the TPS. To meet these requirements, an advanced thermal protection system was conceived, designated 'TOPHAT'. This system consists of a toughened outer ceramic matrix composite (CMC) attached to a rigid reusable surface insulator (RSI) which is directly bonded to the surface. The objective of this effort was to evaluate this concept in an aeroconvective environment, to determine the effect of impacts to the CMC material, and to compare the results with existing thermal protection systems.

Riccitiello, Salvatore R.↗
NASA NTRS2001

Crystallization of bFGF-DNA Aptamer Complexes Using a Sparse Matrix Designed for Protein-Nucleic Acid Complexes

The Sparse Matrix approach for obtaining lead crystallization conditions has proven to be very fruitful for the crystallization of proteins and nucleic acids. Here we report a Sparse Matrix developed specifically for the crystallization of protein-DNA complexes. This method is rapid and economical, typically requiring 2.5 mg of complex to test 48 conditions. The method was originally developed to crystallize basic fibroblast growth factor (bFGF) complexed with DNA sequences identified through in vitro selection, or SELEX, methods. Two DNA aptamers that bind with approximately nanomolar affinity and inhibit the angiogenic properties of bFGF were selected for co-crystallization. The Sparse Matrix produced lead crystallization conditions for both bFGF-DNA complexes.

Cannone, Jaime J.↗
NASA NTRS1990

The mechanisms of dispersion strengthening and fracture in Al-based XD(tm) alloys, part 1

The influence of reinforcement size, volume fraction, and matrix deformation behavior on room and elevated temperature strength; the fracture toughness; and the fatigue crack growth rate of metal matrix composites of Al-4(pct)Cu-1.5(pct)Mg with TiB2 were examined. The influence of reinforcement volume fraction was also examined for pure aluminum with TiB2. Higher TiB2 volume fractions increased the tensile yield strength at both room and elevated temperatures, and reduced the elongation to fracture. Tensile tests also indicate that small particles provided a greater increase in strength for a given volume fraction than larger particles, whereas elongation to fracture appeared to be insensitive to reinforcement size. Interparticle spacing appears to be the factor that controls the strength of these alloys, with the exact nature of the dependence relying on the nature of dislocation slip in the matrix (planar vs. diffuse). The isothermal aging response of the precipitation strengthened Al-4(pct)Cu-1.5(pct)Mg alloys was not accelerated by the presence of TiB2. Cold work prior to artificial aging created additional geometrically necessary dislocations which serve as heterogeneous nucleation sites leading to accelerated aging, a finer precipitate size, and an increase in the strength of the alloy.

Aikin, R. M., Jr.↗
DOE OSTI2023

Irradiation Effects on Stability of δ-UZr2 phase in U-50 wt% Zr Alloy

U-50wt%Zr is a candidate metallic nuclear fuel with potential application in light water reactors due to its excellent thermal properties and high radiation tolerance. The Zr-rich UZr fuels possess greater swelling resistance and fission gas release characteristics compared with U-rich UZr fuels. In this current study, the δ-phase U-50wt%Zr is proton irradiated at room temperature to 1 displacement per atom (dpa) to provide insights on phase stability under irradiation conditions. High resolution characterization of Transmission Electron Microscopy (TEM) and Atom Probe Tomography (APT) characterization techniques are used to elucidate microstructural changes due to irradiation. TEM and APT results show highly oriented bcc β-Zr-rich platelet precipitates nucleating adjacent to α-U phases inside the UZr 2 matrix. Formation of this platelet morphology is characteristic of Widmänstatten structure which can be attributed to a variety of factors such as differences in thermal expansion coefficient between the phases, grain size, alloy composition, and cooling rate. The phases present are distinctly different than those observed through in situ annealing, but irradiation accelerates diffusion and phase separation kinetics. These microstructural changes in U-50wt%Zr are different from those achieved by pure thermodynamic or high temperature heavy ion irradiation experiments. Finally, our work, together with previous ones, highlight the necessity to study the U-Zr phase diagram under non-equilibrium thermodynamics conditions to support this material's deployment as a viable nuclear fuel form.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Estrogen depletion alters osteogenic differentiation and matrix production by osteoblasts in vitro

Highlights: • Estrogen depletion accelerates osteoblast differentiation in static culture. • Mechanical stimulation exacerbated the effects of estrogen depletion on osteoblast differentiation. • Estrogen depletion alters osteoblasts capacity to produce a good quality matrix. Recent studies have revealed that the effects of estrogen deficiency are not restricted to osteoclasts and bone resorption, but that bone matrix composition is altered and osteoblasts exhibit an impaired response to mechanical stimulation. In this study, we test the hypothesis that estrogen depletion alters osteogenic differentiation and matrix production by mechanically stimulated osteoblasts in vitro. MC3T3-E1 cells were pre-treated with estrogen for 14 days, after which estrogen was withdrawn or inhibited with Fulvestrant up to 14 days. Fluid shear stress (FSS) was applied using an orbital shaker. Under estrogen depletion in static culture, osteogenic marker (ALP) and gene expression (Runx2) were decreased at 2 and after 7 days of estrogen depletion, respectively. In addition, up to 7 day the inhibition of the estrogen receptor significantly decreased fibronectin expression (FN1) under static conditions. Under estrogen depletion and daily mechanical stimulation, changes in expression of Runx2 occurred earlier (4 days) and by 14 days, changes in matrix production (Col1a1) were reported. We propose that changes in osteoblast differentiation and impaired matrix production during estrogen depletion may contribute to the altered quality of the bone and act as a contributing factor to increased bone fragility in postmenopausal osteoporosis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Photoproduction of fully charmed tetraquark at electron-ion colliders

In this work, we investigate the inclusive photoproduction of the C -odd, S -wave fully charmed tetraquark at electron-ion colliders within the nonrelativistic QCD (NRQCD) factorization framework, at the lowest order in velocity and α s . The value of the NRQCD long-distance matrix element is estimated from two phenomenological potential models. Our studies reveal that the photoproduction of the 1 + − fully charmed tetraquark may be difficult to observe at HERA and the EicC; nevertheless, its observation prospects at the EIC appear to be bright. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Improved Statistical Analysis for the Neutrinoless Double-Beta Decay Matrix Element of 136Xe

Neutrinoless double beta decay nuclear matrix element (M0ν) for 136Xe was recently analyzed using a statistical approach (Phys. Rev. C 107, 045501 (2023)). In the analysis, three initial shell model effective Hamiltonians were randomly altered, and their results for 23 measured observables were used to infer credibility for the M0ν nuclear matrix element (NME) based on a Bayesian Model Averaging approach. In that analysis, a reasonable Gamow-Teller quenching factor of 0.7 was assumed for each starting effective Hamiltonian. Given that the result of the statistical analysis was sensible to this choice, we are here improving that analysis by assuming that the Gamow-Teller quenching factor is also randomly chosen within reasonabe limits for all three starting Hamiltonians. The outcomes are slightly higher expectation values and uncertainties for the M0ν NME.

Astronomy & Astrophysics↗
NASA NTRS1979

Normal and radial impact of composites with embedded penny-shaped cracks

A method is developed for the dynamic stress analysis of a layered composite containing an embedded penny-shaped crack and subjected to normal and radial impact. The material properties of the layers are chosen such that the crack lies in a layer of matrix material while the surrounding material possesses the average elastic properties of a two-phase medium consisting of a large number of fibers embedded in the matrix. Quantitatively, the time-dependent stresses near the crack border can be described by the dynamic stress intensity factors. Their magnitude depends on time, on the material properties of the composite and on the relative size of the crack compared to the composite local geometry. Results obtained show that, for the same material properties and geometry of the composite, the dynamic stress intensity factors for an embedded (penny-shaped) crack reach their peak values within a shorter period of time and with a lower magnitude than the corresponding dynamic stress intensity factors for a through-crack.

Sih, G. C.↗
NASA NTRS1990

Interface effects in interlaminar fracture of thermoplastic composites

Mode I interlaminar fracture toughness was determined on a series of unidirectional poly(phenylene oxide)/carbon fiber composites using the double cantilever beam test. Initial toughness for growth of a delamination from an insert depends on fiber type, and ranges from 190 to 440 J/sq m, with intermediate modulus fibers tending to give lower values than high-strain fibers. Low toughness values are attributed to poor fiber-matrix adhesion. As a delamination progress down the beam, fiber bridging increases the apparent toughness by up to a factor of 7, depending on the fiber.

Hinkley, J. A.↗
NASA NTRS1992

Role of interfacial thermal barrier in the transverse thermal conductivity of uniaxial SiC fiber-reinforced reaction bonded silicon nitride

The transverse thermal conductivity of reaction-bonded Si3N4 is significantly affected by an interfacial barrier at the interface formed with SiC reinforcing fibers. A comparative study of composites with and without reinforcing-fiber carbon coatings found the coating to reduce effective thermal conductivity by a factor of about 2; this, however, is partially due to a thermal expansion-mismatch gap between fiber and matrix. HIPing of composites with coated fibers led to an enhancement of thermal conductivity via improved interfacial thermal contact and greater grain size and crystallinity of the fibers.

Bhatt, H.↗
NASA NTRS1999

Quantifying Uncertainties in the Thermo-Mechanical Properties of Particulate Reinforced Composites

The present paper reports results from a computational simulation of probabilistic particulate reinforced composite behavior. The approach consists use of simplified micromechanics of particulate reinforced composites together with a Fast Probability Integration (FPI) technique. Sample results are presented for a Al/SiC(sub p)(silicon carbide particles in aluminum matrix) composite. The probability density functions for composite moduli, thermal expansion coefficient and thermal conductivities along with their sensitivity factors are computed. The effect of different assumed distributions and the effect of reducing scatter in constituent properties on the thermal expansion coefficient are also evaluated. The variations in the constituent properties that directly effect these composite properties are accounted for by assumed probabilistic distributions. The results show that the present technique provides valuable information about the scatter in composite properties and sensitivity factors, which are useful to test or design engineers.

Mital, Subodh K.↗
DOE OSTI2023

Sparse Approximate Multifrontal Factorization with Composite Compression Methods

This article presents a fast and approximate multifrontal solver for large sparse linear systems. In a recent work by Liu et al., we showed the efficiency of a multifrontal solver leveraging the butterfly algorithm and its hierarchical matrix extension, HODBF (hierarchical off-diagonal butterfly) compression to compress large frontal matrices. The resulting multifrontal solver can attain quasi-linear computation and memory complexity when applied to sparse linear systems arising from spatial discretization of high-frequency wave equations. To further reduce the overall number of operations and especially the factorization memory usage to scale to larger problem sizes, in this article we develop a composite multifrontal solver that employs the HODBF format for large-sized fronts, a reduced-memory version of the nonhierarchical block low-rank format for medium-sized fronts, and a lossy compression format for small-sized fronts. This allows us to solve sparse linear systems of dimension up to 2.7 × larger than before and leads to a memory consumption that is reduced by 70% while ensuring the same execution time. The code is made publicly available in GitHub.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2011

Linearized T-Matrix and Mie Scattering Computations

We present a new linearization of T-Matrix and Mie computations for light scattering by non-spherical and spherical particles, respectively. In addition to the usual extinction and scattering cross-sections and the scattering matrix outputs, the linearized models will generate analytical derivatives of these optical properties with respect to the real and imaginary parts of the particle refractive index, and (for non-spherical scatterers) with respect to the ''shape'' parameter (the spheroid aspect ratio, cylinder diameter/height ratio, Chebyshev particle deformation factor). These derivatives are based on the essential linearity of Maxwell's theory. Analytical derivatives are also available for polydisperse particle size distribution parameters such as the mode radius. The T-matrix formulation is based on the NASA Goddard Institute for Space Studies FORTRAN 77 code developed in the 1990s. The linearized scattering codes presented here are in FORTRAN 90 and will be made publicly available.

Analytic linearization↗
DOE OSTI2021

Nodal Admittance Matrix Based Area Partition Method for Small-Signal Stability Analysis of Large-Scale Power Electronics Based Power Systems

In power electronics-based power systems (PEPSs), small-signal stability is an important factor for system design and operation, where the impedance-based approach is often used. However, unlike small-scale PEPSs with simple and straightforward impedance models, it would take much more efforts to derive the large-scale PEPSs impedance model, which is very complicated and sometimes may get wrong results due to the elimination of right-half plane (RHP) poles during the impedance aggregation process. To simplify the derivation procedure and analyze the small-signal stability of large-scale PEPSs, this paper proposes a nodal admittance matrix (NAM) based area partition method. In this method, the large-scale PEPS is divided into several sub-areas, and the stability is analyzed within the sub-area first, and then the interconnection stability among these sub-areas is analyzed. The proposed method is scalable and can help to locate the weakest areas/converters that may cause instability in the whole system. In this paper, the concept of the proposed method and its application to an example system are introduced. Experimental results are also given to validate the effectiveness of the proposed method.

Qiao, Liang↗
DOE OSTI2021

Revisiting matrix-based inversion of scanning mobility particle sizer (SMPS) and humidified tandem differential mobility analyzer (HTDMA) data

Abstract. Tikhonov regularization is a tool for reducing noise amplification during data inversion. This work introduces RegularizationTools.jl, a general-purpose software package for applying Tikhonov regularization to data. The package implements well-established numerical algorithms and is suitable for systems of up to ∼ 1000 equations. Included is an abstraction to systematically categorize specific inversion configurations and their associated hyperparameters. A generic interface translates arbitrary linear forward models defined by a computer function into the corresponding design matrix. This obviates the need to explicitly write out and discretize the Fredholm integral equation, thus facilitating fast prototyping of new regularization schemes associated with measurement techniques. Example applications include the inversion involving data from scanning mobility particle sizers (SMPSs) and humidified tandem differential mobility analyzers (HTDMAs). Inversion of SMPS size distributions reported in this work builds upon the freely available software DifferentialMobilityAnalyzers.jl. The speed of inversion is improved by a factor of ∼ 200, now requiring between 2 and 5 ms per SMPS scan when using 120 size bins. Previously reported occasional failure to converge to a valid solution is reduced by switching from the L-curve method to generalized cross-validation as the metric to search for the optimal regularization parameter. Higher-order inversions resulting in smooth, denoised reconstructions of size distributions are now included in DifferentialMobilityAnalyzers.jl. This work also demonstrates that an SMPS-style matrix-based inversion can be applied to find the growth factor frequency distribution from raw HTDMA data while also accounting for multiply charged particles. The outcome of the aerosol-related inversion methods is showcased by inverting multi-week SMPS and HTDMA datasets from ground-based observations, including SMPS data obtained at Bodega Marine Laboratory during the CalWater 2/ACAPEX campaign and co-located SMPS and HTDMA data collected at the US Department of Energy observatory located at the Southern Great Plains site in Oklahoma, USA. Results show that the proposed approaches are suitable for unsupervised, nonparametric inversion of large-scale datasets as well as inversion in real time during data acquisition on low-cost reduced-instruction-set architectures used in single-board computers. The included software implementation of Tikhonov regularization is freely available, general, and domain-independent and thus can be applied to many other inverse problems arising in atmospheric measurement techniques and beyond.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

The Efficiency of Various Computers and Optimizations in Performing Finite Element Computations

With the advent of computers with many processors, it becomes unclear how to best exploit this advantage. For example, matrices can be inverted by applying several processors to each vector operation, or one processor can be applied to each matrix. The former approach has diminishing returns beyond a handful of processors, but how many processors depends on the computer architecture. Applying one processor to each matrix is feasible with enough ram memory and scratch disk space, but the speed at which this is done is found to vary by a factor of three depending on how it is done. The cost of the computer must also be taken into account. A computer with many processors and fast interprocessor communication is much more expensive than the same computer and processors with slow interprocessor communication. Consequently, for problems that require several matrices to be inverted, the best speed per dollar for computers is found to be several small workstations that are networked together, such as in a Beowulf cluster. Since these machines typically have two processors per node, each matrix is most efficiently inverted with no more than two processors assigned to it.

Marcus, Martin H.↗
DOE OSTI2024

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2014

Injury Surveillance Among NASA Astronauts Using the Barell Injury Diagnosis Matrix

Astronauts perform physically demanding tasks and risk incurring musculoskeletal injuries during both groundbased training and missions. Increased injury rates throughout the history of the U.S. space program have been attributed to numerous factors, including an aging astronaut corps, increased Weightless Environment Training Facility (WETF) and Neutral Buoyancy Laboratory (NBL) training to construct the International Space Station, and improved clinical operations that promote injury prevention and reporting. With NASA program changes through the years (including retirement of the Shuttle program) and an improved training environment (including a new astronaut gym), there is no surveillance program to systematically track injury rates. A limited number of research projects have been conducted over the past 20 years to evaluate musculoskeletal injuries: (1) to evaluate orthopedic injuries from 1987 to 1995, (2) to describe upper extremity injuries, (3) to evaluate EVA spacesuit training related injuries, and (4) to evaluate in-flight musculoskeletal injuries. Nevertheless, there has been no consistently performed comprehensive assessment of musculoskeletal injuries among astronauts. The Barell Injury Diagnosis Matrix was introduced at the 2001 meeting of the International Collaborative Effort (ICE) on Injury Statistics. The Matrix proposes a standardized method of classifying body region by nature of injury. Diagnoses are coded using the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) coding system. The purpose of this study is to assess the usefulness and complexity of the Barell Injury Diagnosis Matrix to classify and track musculoskeletal injuries among NASA astronauts.

Murray, J. D.↗