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 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2022

How Bayesian methods can improve R -matrix analyses of data: The example of the d t reaction

The 3 H(d, n) 4 He reaction is of significant interest in nuclear astrophysics and nuclear applications. It is an important, early step in big-bang nucleosynthesis and a key process in nuclear fusion reactors. We use one- and two-level R-matrix approximations to analyze data on the cross section for this reaction at center-of-mass energies below 215 keV. We critically examine the data sets using a Bayesian statistical model that allows for both common-mode and additional point-to-point un- certainties. We use Markov Chain Monte Carlo sampling to evaluate this R-matrix-plus-statistical model and find two-level R-matrix results that are stable with respect to variations in the channel radii. The S factor at 40 keV evaluates to 25.36(19) MeV b (68% credibility interval). We discuss our Bayesian analysis in detail and provide guidance for future applications of Bayesian methods to R-matrix analyses. We also discuss possible paths to further reduction of the S-factor uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1999

Factors Controlling Stress Rupture of Fiber-Reinforced Ceramic Composites

The successful application of fiber-reinforced ceramic matrix composites (CMC) depends strongly on maximizing material rupture life over a wide range of temperatures and applied stresses. The objective of this paper is to examine the various intrinsic and extrinsic factors that control the high-temperature stress rupture of CMC for stresses below and above those required for cracking of the 0 C plies (Regions I and II, respectively). Using creep-rupture results for a variety of ceramic fibers and rupture data for CMC reinforced by these fibers, it is shown that in those cases where the matrix carries little structural load, CMC rupture conditions can be predicted very well from the fiber behavior measured under the appropriate test environment. As such, one can then examine the intrinsic characteristics of the fibers in order to develop design guidelines for selecting fibers and fiber microstructures in order to maximize CMC rupture life. For those cases where the fiber interfacial coatings are unstable in the test environment, CMC lives are generally worse than those predicted by fiber behavior alone. For those cases where the matrix can support structural load, CMC life can even be greater provided matrix creep behavior is properly controlled. Thus the achievement of long CMC rupture life requires understanding and optimizing the behavior of all constituents in the proper manner.

DiCarlo, J. A.↗

NASA NTRS1993

Acousto-Ultrasonic analysis of failure in ceramic matrix composite tensile specimens

Three types of acousto-ultrasonic (AU) measurements, stress-wave factor (SWF), lowest antisymmetric plate mode group velocity (VS), and lowest symmetric plate mode group velocity (VL), were performed on specimens before and after tensile failure. Three different Nicalon fiber architectures with ceramic matrices were tested. These composites were categorized as 1D (unidirectional fiber orientation) SiC/CAS glass ceramic, and 2D and 3D woven SiC/SiC ceramic matrix materials. SWF was found to be degraded after tensile failure in all three material categories. VS was found to be degraded only in the 1D SiC/CAS. VL was difficult to determine on the irregular specimen surfaces but appeared unchanged on all failed specimens. 3D woven specimens with heat-treatment at high temperature exhibited degradation only in SWF.

Kautz, Harold E.↗

DOE OSTI2020

Measurement of the Li 7 ( γ , t ) He 4 ground-state cross section between E γ = 4.4 and 10 MeV

Here, the Li 7 ( γ , t ) He 4 ground state cross section was measured for the first time using monoenergetic γ rays with energies between 4.4 and 10 MeV at the High Intensity Gamma-ray Source. The reaction is important for the primordial Li problem and for testing our understanding of the mirror α -capture reactions H 3 ( α , γ ) Li 7 and He 3 ( α , γ ) Be 7 . Although over the last 30 years most measurements of the H 3 ( α , γ ) Li 7 reaction have concentrated in an energy range below E γ = 3.65 MeV, measurements at higher energies could potentially restrict the extrapolation to astrophysically important energies. The experimental arrangement for measuring the Li 7 ( γ , t ) He 4 reaction included a large-area silicon detector array and several beam characterization instruments. The experimental astrophysical S factor of H 3 ( α , γ ) calculated from the present data was fitted using the R -matrix formalism. The results are in disagreement with previous experimental measurements in the same energy range but the extrapolated S factor agrees with the potential model calculation and lower energy experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1984

An efficient approximate factorization implicit scheme for the equations of gasdynamics

An efficient implicit finite-difference algorithm for the gas dynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved while maintaining the same favorable stability characteristics and generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows about a NACA 0012 airfoil. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations. It is shown that the method offers advantages over the conventional Beam and Warming scheme and can retrofit existing Beam and Warming codes with minimal effort.

Barth, T. J.↗

DOE OSTI2020

Electron transfer calculations between edge sharing octahedra in hematite, goethite, and annite

A key reaction underlying the charge transport in iron containing oxides, clays, micas is the Fe$^{2+}$-Fe$^{3+}$ exchange reaction between edge-sharing iron octahedra. These reactions facilitate conduction in these minerals by the thermally-activated hopping of small polarons across the lattice. Depending on the mineral and local charge state the small polaron can either encase an electron or hole. The probability for conduction of small polarons depends strongly on the height and adiabicity of the reaction barrier, with larger and more diabatic barriers yielding slow conduction associated with either weak coupling or a large prerequisite rearrangement of the lattice during charge transport. To model these reactions, a first principle electron transfer (ET) method was developed to model the small polaron hopping between the edge-sharing octahedra sites in hematite ($e^{-}$ polaron), goethite ($e^{-}$ polaron), and annite ($h^{+}$ polaron) bulk structures. The ET method is based on electronic structure methods (i.e., plane-wave Density Functional Theory) capable of performing calculations with periodic cells and large size systems efficiently while at the same time being accurate enough to be used in the estimation of the electron-transfer coupling matrix element, $V_{AB}$, and the electron transfer transmission factor, $\kappa_{el}$. Additionally, the calculations confirmed the existence of small polarons in all three minerals, and the reactions were predicted to be strongly adiabatic. It was found that transfer of a hole in the octahedral layer of annite had an adiabatic barrier of $0.311$ eV, and the transfer of an extra electron in hematite and goethite had adiabatic barriers of $0.242$ eV and $0.232$ eV respectively. The electronic coupling parameters,$V_{AB}$, were found to be $0.188$ eV, $0.196$ eV, and $0.102$ eV respectively for hematite, goethite, and annite. While similar bonding topologies pertain, the underlying basis for the differences is the subtle differences in local structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

metabCombiner 2.0: Disparate Multi-Dataset Feature Alignment for LC-MS Metabolomics

Liquid chromatography–high-resolution mass spectrometry (LC-HRMS), as applied to untargeted metabolomics, enables the simultaneous detection of thousands of small molecules, generating complex datasets. Alignment is a crucial step in data processing pipelines, whereby LC-MS features derived from common ions are assembled into a unified matrix amenable to further analysis. Variability in the analytical factors that influence liquid chromatography separations complicates data alignment. This is prominent when aligning data acquired in different laboratories, generated using non-identical instruments, or between batches from large-scale studies. Previously, we developed metabCombiner for aligning disparately acquired LC-MS metabolomics datasets. Here, we report significant upgrades to metabCombiner that enable the stepwise alignment of multiple untargeted LC-MS metabolomics datasets, facilitating inter-laboratory reproducibility studies. To accomplish this, a “primary” feature list is used as a template for matching compounds in “target” feature lists. We demonstrate this workflow by aligning four lipidomics datasets from core laboratories generated using each institution’s in-house LC-MS instrumentation and methods. We also introduce batchCombine, an application of the metabCombiner framework for aligning experiments composed of multiple batches. metabCombiner is available as an R package on Github and Bioconductor, along with a new online version implemented as an R Shiny App.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1978

Conservative implicit schemes for the full potential equation applied to transonic flows

Implicit approximate factorization techniques (AF) were investigated for the solution of matrix equations resulting from finite difference approximations to the full potential equation in conservation form. For transonic flows, an artificial viscosity, required to maintain stability in supersonic regions, was introduced by an upwind bias of the density. Two implicit AF procedures are presented and their convergence performance is compared with that of the standard transonic solution procedure, successive line overrelaxation (SLOR). Subcritical and supercritical test cases are considered. The results indicate that the AF schemes are substantially faster than SLOR.

Holst, T. L.↗

NASA NTRS1979

Fast, conservative schemes for the full potential equation applied to transonic flows

Implicit approximate factorization techniques (AF) are investigated for the solution of matrix equations resulting from finite-difference approximations to the full potential equation in conservation form. For transonic flows, an artificial viscosity, required to maintain stability in supersonic regions, is introduced by an upwind bias of the density. Two implicit AF procedures are presented, and their convergence performance is compared with that of the standard transonic solution procedure: successive line overrelaxation (SLOR). Subcritical and supercritical test cases are considered. Results indicate a substantial improvement in convergence rate for AF schemes relative to SLOR.

Holst, T. L.↗

NASA NTRS1984

Water vapor profile retrieval simulation studies for the HIRS/MSU and AMTS/MSU sounders

A combined physical statistical method for the retrieval of atmospheric water vapor column densities has been developed at the Goddard Laboratory for Atmospheric Sciences. In this method differences between observed and calculated clear column brightness temperatures are used in a regression formalism to obtain estimates of integrated column densities in given atmospheric layers. The physics of the radiative transfer problem is contained in the calculated clear column brightness temperatures and in accounting for the effects of clouds on the observed radiances. The regression matrix is equivalent to a set of empirical sensitivity factors.

Reuter, D.↗

NASA NTRS1985

Microstructure and fracture in SiC whisker reinforced 2124 aluminum composite

The microstructures of extruded and hot-rolled 2124 Al-15 percent (by weight) SiC whisker composites have been investigated, experimentally. Among the specific factors studied were: the strength of the whisker-matrix interfaces; (2) the presence of oxides; (3) the presence of defective whiskers; (4) and the presence of distribution of intermetallic compounds, impurities in the SiC(w) powder, and microstructural inhomogeneities. Modifications in the microstructure of the SiC/AL composites due to hot rolling and extrusion are illustrated in a series of microphotographs. It was found that hot rolling along the axis of extrusion was associated with some types of whisker damage, while the whiskers still retain their original orientation. Hot-rolling perpendicular to the axis of extrusion, however, tended to rotate the whiskers and produced a nearly isotropic material. Whisker free zones were virtually eliminated or reduced in size by hot rolling. In situ Auger fractography of the composite showed that the interfacial bonding between the SiC and the Al matrix was good and that Al2O2 had no significant influence on the fracture mechanics of the composite.

Nieh, T. G.↗

NASA NTRS1988

Data traffic reduction schemes for sparse Cholesky factorizations

Load distribution schemes are presented which minimize the total data traffic in the Cholesky factorization of dense and sparse, symmetric, positive definite matrices on multiprocessor systems with local and shared memory. The total data traffic in factoring an n x n sparse, symmetric, positive definite matrix representing an n-vertex regular 2-D grid graph using n (sup alpha), alpha is equal to or less than 1, processors are shown to be O(n(sup 1 + alpha/2)). It is O(n(sup 3/2)), when n (sup alpha), alpha is equal to or greater than 1, processors are used. Under the conditions of uniform load distribution, these results are shown to be asymptotically optimal. The schemes allow efficient use of up to O(n) processors before the total data traffic reaches the maximum value of O(n(sup 3/2)). The partitioning employed within the scheme, allows a better utilization of the data accessed from shared memory than those of previously published methods.

Naik, Vijay K.↗

NASA NTRS1989

Communication requirements of sparse Cholesky factorization with nested dissection ordering

Load distribution schemes for minimizing the communication requirements of the Cholesky factorization of dense and sparse, symmetric, positive definite matrices on multiprocessor systems are presented. The total data traffic in factoring an n x n sparse symmetric positive definite matrix representing an n-vertex regular two-dimensional grid graph using n exp alpha, alpha not greater than 1, processors are shown to be O(n exp 1 + alpha/2). It is O(n), when n exp alpha, alpha not smaller than 1, processors are used. Under the conditions of uniform load distribution, these results are shown to be asymptotically optimal.

Naik, Vijay K.↗

NASA NTRS1995

A far-wing line shape theory which satisfies the detailed balance principle

A far-wing theory in which the validity of the detailed balance principle is maintained in each step of the derivation is presented. The role of the total density matrix including the initial correlations is analyzed rigorously. By factoring out the rapidly varying terms in the complex-time development operator in the interaction representation, better approximate expressions can be obtained. As a result, the spectral density can be expressed in terms of the line-coupling functions in which two coupled lines are arranged symmetrically and whose frequency detunings are omega - 1/2(omega(sub ji) + omega (sub j'i'). Using the approximate values omega - omega(sub ji) results in expressions that do not satisfy the detailed balance principle. However, this principle remains satisfied for the symmetrized spectral density in which not only the coupled lines are arranged symmetrically, but also the initial and final states belonging to the same lines are arranged symmetrically as well.

Ma, Q.↗

NASA NTRS2015

3D Progressive Damage Modeling for Laminated Composite Based on Crack Band Theory and Continuum Damage Mechanics

A simple continuum damage mechanics (CDM) based 3D progressive damage analysis (PDA) tool for laminated composites was developed and implemented as a user defined material subroutine to link with a commercially available explicit finite element code. This PDA tool uses linear lamina properties from standard tests, predicts damage initiation with an easy-to-implement Hashin-Rotem failure criteria, and in the damage evolution phase, evaluates the degradation of material properties based on the crack band theory and traction-separation cohesive laws. It follows Matzenmiller et al.'s formulation to incorporate the degrading material properties into the damaged stiffness matrix. Since nonlinear shear and matrix stress-strain relations are not implemented, correction factors are used for slowing the reduction of the damaged shear stiffness terms to reflect the effect of these nonlinearities on the laminate strength predictions. This CDM based PDA tool is implemented as a user defined material (VUMAT) to link with the Abaqus/Explicit code. Strength predictions obtained, using this VUMAT, are correlated with test data for a set of notched specimens under tension and compression loads.

Wang, John T.↗

NASA NTRS2018

Microscale Fracture Testing of HfO2-Si Environmental Barrier Coatings

Environmental barrier coatings (EBCs) play a vital role in protecting advanced turbine components constructed from ceramic matrix composites. EBC performance is evaluated based on many factors such as mechanical and thermal properties and chemical stability. This study focuses on the mechanical properties of individual phases within a HfO2-Si EBC through the combined use of chevron-notched microcantilevers and finite element modeling. The Si- and HfO2-rich phases were found to have fracture toughness K(sub IC) values of 1.36+/-0.07 MPa·m1/2 and 1.99+/-0.36 MPa·m1/2, respectively. Despite irregular notch geometry, calculated K(sub IC) values agreed reasonably with literature values of monolithic samples for the respective phases. The relatively large scatter in the HfO2-rich samples is attributed to heterogeneity in fracture feature size and distribution of the HfO2 and HfSiO4 phases. Comparisons with bulk single-edge V-notched beams suggest bulk behavior is dictated by the fracture behavior of the Sirich phase with no contribution from the HfO2-rich phase.

Robertson, Andrew L.↗

NASA NTRS2005

Tissue Engineering Using Transfected Growth-Factor Genes

A method of growing bioengineered tissues includes, as a major component, the use of mammalian cells that have been transfected with genes for secretion of regulator and growth-factor substances. In a typical application, one either seeds the cells onto an artificial matrix made of a synthetic or natural biocompatible material, or else one cultures the cells until they secrete a desired amount of an extracellular matrix. If such a bioengineered tissue construct is to be used for surgical replacement of injured tissue, then the cells should preferably be the patient s own cells or, if not, at least cells matched to the patient s cells according to a human-leucocyteantigen (HLA) test. The bioengineered tissue construct is typically implanted in the patient's injured natural tissue, wherein the growth-factor genes enhance metabolic functions that promote the in vitro development of functional tissue constructs and their integration with native tissues. If the matrix is biodegradable, then one of the results of metabolism could be absorption of the matrix and replacement of the matrix with tissue formed at least partly by the transfected cells. The method was developed for articular chondrocytes but can (at least in principle) be extended to a variety of cell types and biocompatible matrix materials, including ones that have been exploited in prior tissue-engineering methods. Examples of cell types include chondrocytes, hepatocytes, islet cells, nerve cells, muscle cells, other organ cells, bone- and cartilage-forming cells, epithelial and endothelial cells, connective- tissue stem cells, mesodermal stem cells, and cells of the liver and the pancreas. Cells can be obtained from cell-line cultures, biopsies, and tissue banks. Genes, molecules, or nucleic acids that secrete factors that influence the growth of cells, the production of extracellular matrix material, and other cell functions can be inserted in cells by any of a variety of standard transfection techniques.

Madry, Henning↗

NASA NTRS1990

Finite element analysis and optimization of composite structures

Linearly elastic fiber reinforced composite discs and laminates in plane stress with variable local orientation and concentration of one or two fiber fields embedded in the matrix material, are considered. The thicknesses and the domain of the discs or laminates are assumed to be given, together with prescribed boundary conditions and in-plane loading along the edge. The problem under study consists in determining throughout the structural domain the optimum orientations and concentrations of the fiber fields in such a way as to maximize the integral stiffness of the composite disc or laminate under the seven loading. Minimization of the integral stiffness can also be performed. The optimization is performed subject to a prescribed bound on the total cost or weight of the composite that for given unit cost factors or specific weights determines the amounts of fiber and matrix materials in the structure. Examples are presented.

Thomsen, Jan↗