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 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2022

SparseLU, A Novel Algorithm and Math Library for Sparse LU Factorization

Decomposing sparse matrices into lower and upper triangular matrices (sparse LU factorization) is a key operation in many computational scientific applications. We developed SparseLU, a sparse linear algebra library that implements a new algorithm for LU factorization on general sparse matrices. The new algorithm divides the input matrix into tiles to which OpenMP tasks are created for factorization computation, where only tiles that contain nonzero elements are computed. For comparative performance analysis, we used the reference library SuperLU. Testing was performed on synthetically generated matrices which replicate the conditions of the real-world matrices. SparseLU is able to reach a mean speedup of ~29× compared to SuperLU.

Valero Lara, Pedro↗

NASA NTRS1975

Diagrammatic technique for calculating matrix elements of collective operators in superradiance

Adopting the so-called genealogical construction, one can express the eigenstates of collective operators corresponding to a specified mode for an N-atom system in terms of those for an (N-1) atom system. Using these Dicke states as bases and using the Wigner-Eckart theorem, a matrix element of a collective operator of an arbitrary mode can be written as the product of an m-dependent factor and an m-independent reduced matrix element (RME). A set of recursion formulas for the RME is obtained. A graphical representation of the RME on the branching diagram for binary irreducible representations of permutation groups is then introduced. This gives a simple and systematic way of calculating the RME. This method is especially useful when the cooperation number r is close to N/2, where almost exact asymptotic expressions can be obtained easily. The result shows explicity the geometry dependence of superradiance and the relative importance of r-conserving and r-nonconserving processes.

Lee, C. T.↗

NASA NTRS1980

Effect of molecular anisotropy on the intensity and degree of polarization of light scattered from model atmospheres

Computations of the properties of sunlight scattered from models of the earth-atmosphere system are presented to show the effect of molecular anisotropy on the intensity, flux, and degree of polarization of the scattered light. The values of these parameters change significantly when the anisotropy factor is neglected in the molecular optical thickness and scattering phase matrix. However, if the Rayleigh scattering optical thickness is kept constant and the molecular anisotropy factor is included only in the Rayleigh phase matrix, the flux does not change, the intensity changes by a small amount, but the changes in the degree of polarization are still significant.

Bahethi, O. P.↗

NASA NTRS1975

Superradiant effects on pulse propagation in resonant media

Adopting the so-called genealogical construction, the eigenstates of collective operators can be expressed corresponding to a specified mode for an N-atom system in terms of those for an (N-1)-atom system. Matrix element of a collective operator of an arbitrary mode is presented which can be written as the product of an m-dependent factor and an m-independent reduced matrix element (RME). A set of recursion formulas for the RME was obtained. A graphical representation of the RME on the branching diagram for binary irreducible representations of permutation groups was then introduced. This gave a simple and systematic way of calculating the RME. Results show explicitly the geometry dependence of superradiance and the relative importance of r-conserving and r-nonconserving processes and clears up the chief difficulty encounted in the problem of N two-level atoms, spread over large regions, interacting with a multimode radiation field.

Lee, C.↗

DOE OSTI2021

Changes in rock matrix compressibility during deep CO 2 storage

Geologic carbon storage projects aim to permanently trap large volumes of CO 2 in reservoir rock sealed with low permeability layers. As high-pressure supercritical or liquid CO 2 is injected, hydromechanical and chemical processes caused by the reaction between the rock and acidic mixture of brine and CO 2 are initiated. The compressibility of reservoir rock needs to be properly characterized in order to assess the deformation and stability of the host formations, and there are a number of factors to be considered, including the matrix structure, solid, pores, and fluid. Here this study triggers from a fundamental question whether CO 2 treatment affects the compressibility of the rock matrix and its dominant composing solid minerals. Three different reservoir representatives are selected: Berea sandstone for silica-rich rock, and Apulian limestone and Indiana limestone for calcite-rich rock. Quartz and calcite are the main composing minerals of the reservoir rock, and their crystals are separately examined. Experimental methods are introduced for high-pressure CO 2 treatment of water-saturated materials, and measurements of the unjacketed and solid compressibilities are conducted. No change in the solid compressibility of the sandstone and quartz after CO 2 treatment is observed, whereas it increases by 18–21% for the limestones and by 15% for calcite. The latter observation is inconsistent with the ultrasonic wave velocities measurements and is believed to be due to the local dissolution of the calcite crystal surface. Our results show that only the solid matrix of the limestones becomes more compressible after CO 2 treatment. Consequent microimaging and mercury intrusion porosimetry analyses allowed observations of dissolution and precipitation of calcite, and creation of new connected and non-connected pores. Finally, the changes in limestone solid compressibilities and pore structure could significantly affect the rock properties and behavior during and after CO 2 injection and should be accounted for in the reservoir models.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

An Improved Fracture Mechanics-Informed Multiscale Thermomechanical Damage Model for Ceramic Matrix Composites

This paper extends recent work done by the authors in modeling length scale-dependent damage behavior of ceramic matrix composites (CMCs) to include effects of local anisotropy introduced by matrix cracking. This model captures scale-dependent damage initiation and propagation behavior of the brittle matrix by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress/strain constitutive relationships at each length scale. The damage ISV captures the effects of matrix cracking and growth by using fracture mechanics and the self-consistent scheme to determine the reduced stiffness of the cracked matrix. Matrix cracks, which activate when stress intensity factors near manufacturing induced cavities exceed the fracture toughness of the material, are assumed to be transversely isotropic in the plane of the crack, and matrix anisotropy occurs when the damaged stiffness tensor is rotated from the crack plane to the global axes. The crack progression and temporal evolution of the damage ISV are governed by fracture mechanics and crack growth kinetics. The model effectively captures first matrix cracking, which is the first significant deviation from linear elasticity. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide (SiC) and a 2D woven five-harness satin (5HS) carbon fiber SiC matrix (C/SiC) CMC.

Skinner, Travis↗

NASA NTRS1991

Darkening in gas-rich ordinary chondrites: Spectral modelling and implications for the regoliths of ordinary chondrite parent bodies

The dark fine grained matrix of gas-rich ordinary chondrites replicates many of the physical, morphological, and spectral characteristics of the highly shocked and optically altered black chondrites. Spectral mixture modeling shows that the darkening and spectral attenuation seen in the dark matrix can be simulated with realistic mass fractions of light host material and black chondritic material. All these factors point to the conclusion that the dark matrix of gas-rich ordinary chondrites is dark due to the same processes that darkens black chondrites, shock distributed small particle size FeNi metal and troilite. Because the darkening is not seen in any of the non-gas-rich light portions and is only seen in the gas rich grains of the meteorite, the shock darkening would have to occur as part of the matrix's exposure to regolith processes. Since all gas-rich grains are darkened, it follows that darkening is not only common, but pervasive in asteroidal regoliths. These results imply that the upper, optically active layer of an ordinary chondrite parent body should have the spectral characteristics of a black chondrite, which are a dark, relatively featureless spectrum with modest red slope in the infrared. These are the characteristics of special type C asteroids.

Britt, Daniel T.↗

DOE OSTI2025

Pion and Kaon decays from lattice QCD and |Vud|, |Vus|

We review lattice calculations of the form factors for Kaon semileptonic decays and the leptonic decay constants of Kaons and Pions. We discuss the determinations of the CKM matrix elements $|V_{ud}|$ and $|V_{us}|$ using the form factors and experimental results. We use the CKM matrix elements to test first-row unitarity of the CKM matrix.

Simone, J. N. [Fermilab] (ORCID:0000000185153337)↗

DOE OSTI2024

Reinforcement induced microcracking during the conversion of polymer-derived ceramics

This study investigates the role of discontinuous reinforcement on the microcracking of a preceramic polymer matrix during polymer to ceramic conversion. The material system comprised a carbosilane-based, preceramic polymer reinforced with mullite particles. The preceramic resin slurry was formulated for photopolymerization on a digital light projection 3D printer. Micro X-ray computed tomography, using a synchrotron light source, monitored a printed composite part during pyrolysis. The conversion of the carbosilane matrix into Si(O)C generated microcracks in the matrix that radially extended from particles. The likelihood of microcrack formation positively correlated with particle size. The largest particles (>10 4 µm 3 volume or >60 µm side length) always abutted microcracks, while the matrix surrounding smaller particles (<5 µm) was free of microcracks. Numerical calculations showed the particles resist the shrinking matrix during its conversion. This created tensile stresses within the matrix. The driving force for microcrack growth was also analyzed with regard to the influence of matrix material parameters, particle form-factor and crack length. These results reveal the need to consider the form factor of discontinuous reinforcements, with regard to the material properties of the polymer-derived ceramic, in order to reduce or eliminate defects in the final part.

36 MATERIALS SCIENCE↗

NASA NTRS1988

A spatial operator algebra for manipulator modeling and control

A powerful new spatial operator algebra for modeling, control, and trajectory design of manipulators is discussed along with its implementation in the Ada programming language. Applications of this algebra to robotics include an operator representation of the manipulator Jacobian matrix; the robot dynamical equations formulated in terms of the spatial algebra, showing the complete equivalence between the recursive Newton-Euler formulations to robot dynamics; the operator factorization and inversion of the manipulator mass matrix which immediately results in O(N) recursive forward dynamics algorithms; the joint accelerations of a manipulator due to a tip contact force; the recursive computation of the equivalent mass matrix as seen at the tip of a manipulator; and recursive forward dynamics of a closed chain system. Finally, additional applications and current research involving the use of the spatial operator algebra are discussed in general terms.

Rodriguez, G.↗

DOE OSTI2020

Many-Body Level Statistics of Single-Particle Quantum Chaos

We consider a noninteracting many-fermion system populating levels of a unitary random matrix ensemble (equivalent to the q = 2 complex Sachdev-Ye-Kitaev model)—a generic model of single-particle quantum chaos. We study the corresponding many-particle level statistics by calculating the spectral form factor analytically using algebraic methods of random matrix theory, and match it with an exact numerical simulation. Despite the integrability of the theory, the many-body spectral rigidity is found to have a surprisingly rich landscape. In particular, we find a residual repulsion of distant many-body levels stemming from single-particle chaos, together with islands of level attraction. These results are encoded in an exponential ramp in the spectral form factor, which we show to be a universal feature of nonergodic many-fermion systems embedded in a chaotic medium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1985

Fracture surface analysis in composite and titanium bonding

To understand the mechanical properties of fiber-reinforced composite materials, it is necessary to understand the mechanical properties of the matrix materials and of the reinforcing fibers. Another factor that can affect the mechanical properties of a composite material is the interaction between the fiber and the matrix. In general, composites with strong fiber matrix bonding will give higher modulus, lower toughness composites. Composites with weak bonding will have a lower modulus and more ductility. The situation becomes a bit more complex when all possibilities are examined. To be considered are the following: the properties of the surface layer on the fiber, the interactive forces between polymer and matrix, the surface roughness and porosity of the fiber, and the morphology of the matrix polymer at the fiber surface. In practice, the surface of the fibers is treated to enhance the mechanical properties of a composite. These treatments include anodization, acid etching, high temperature oxidation, and plasma oxidation, to name a few. The goal is to be able to predict the surface properties of carbon fibers treated in various ways, and then to relate surface properties to fiber matrix bonding.

Devilbiss, T. A.↗

NASA NTRS1992

An approximately factored incremental strategy for calculating consistent discrete aerodynamic sensitivity derivatives

An incremental strategy is presented for iteratively solving very large systems of linear equations, which are associated with aerodynamic sensitivity derivatives for advanced CFD codes. It is shown that the left-hand side matrix operator and the well-known factorization algorithm used to solve the nonlinear flow equations can also be used to efficiently solve the linear sensitivity equations. Two airfoil problems are considered as an example: subsonic low Reynolds number laminar flow and transonic high Reynolds number turbulent flow.

Korivi, V. M.↗

NASA NTRS1993

Modeling fatigue crack growth in cross ply titanium matrix composites

In this study, the fatigue crack growth behavior of fiber bridging matrix cracks in cross-ply SCS-6/Ti-15-3 and SCS-6/Timetal-21S laminates containing center holes was investigated. Experimental observations revealed that matrix cracking was far more extensive and wide spread in the SCS-6/Ti-15-3 laminates compared to that in the SCS-6/Timetal-21S laminates. In addition, the fatigue life of the SCS-6/Ti-15-3 laminates was significantly longer than that of the SCS-6/Timetal-21S laminates. The matrix cracking observed in both material systems was analyzed using a fiber bridging (FB) model which was formulated using the boundary correction factors and weight functions for center hole specimen configurations. A frictional shear stress is assumed in the FB model and was used as a curve fitting parameter to model matrix crack growth data. The higher frictional shear stresses calculated in the SCS-6/Timetal-21S laminates resulted in lower stress intensity factors in the matrix and higher axial stresses in the fibers compared to those in the SCS-6/Ti-15-3 laminates at the same applied stress levels.

Bakuckas, J. G., Jr.↗

DOE OSTI2024

Timelike meson form factors beyond the elastic region from lattice QCD

We present a calculation of the vector-isovector timelike form factors of the pion and the kaon using lattice quantum chromodynamics. We calculate two-point correlation functions with m π ∼ 280 MeV , extracting both the finite-volume spectrum and matrix elements for these states created from the vacuum by a vector current. After determining the coupled-channel π π , K K ¯ scattering amplitudes, we perform the necessary correction for the significant finite-volume effects present in the current matrix elements, leading to the timelike form factors. We find these to be dominated by the presence of the ρ resonance, and we extract its decay constant by an analytic continuation of the amplitudes to the resonance pole. In addition, the spacelike pion form factor is determined on the same lattice configurations, and a dispersive parametrization is used to simultaneously describe the spacelike and elastic timelike regions. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2005

Toughened Matrix SiC Fiber Reinforced Composites

First matrix cracking stress is a critical parameter for application of Sic fiber reinforced composites in highly stressed, environmentally demanding applications such as turbine blades. High matrix fracture toughness is a key property that contributes to high composite fracture stress. Silicon nitride offers reduced matrix elastic modulus, lower coefficient of thermal expansion, and potentially high fracture toughness compared to Sic matrices. All of these factors can be used to advantage to increase matrix fracture stress. As a first model system we are pursuing toughened silicon nitride matrix composites reinforced with SCS-9 fibers. Fabrication is by tape casting the matrix plies and tape lay-up with fiber plies followed by hot pressing at 1800 C. Progress toward this end will be reported.

Levine, Stanley R.↗

DOE OSTI2025

Conformal collider physics meets LHC data

The remarkably high energies of the Large Hadron Collider (LHC) have allowed for the first measurements of the shapes and scalings of multipoint correlators of energy flow operators, ⟨ Ψ | E ( n → 1 ) E ( n → 2 ) ⋯ E ( n → k ) | Ψ ⟩ , providing new insights into the Lorentzian dynamics of quantum chromodynamics (QCD). In this letter, we use recent advances in effective field theory to derive a rigorous factorization theorem for the light-ray density matrix, ρ = | Ψ ⟩ ⟨ Ψ | , inside high transverse momentum jets at the LHC. Using the light-ray operator product expansion, the scaling behavior of multipoint correlators can be computed from the expectation value of the twist-2 spin- J light-ray operators, O [ J ] , in this state, Tr [ ρ O [ J ] ] . We compute the light-ray density matrix at next-to-leading order, and combine this with results for the next-to-leading logarithmic scaling behavior of the correlators up to six-points, comparing with CMS open data. This theoretical accuracy allows us to resolve the quantum scaling dimensions of QCD light-ray operators inside jets at the LHC. Our factorization theorem for the light-ray density matrix at the LHC completes the link between recent developments in the study of energy correlators and LHC phenomenology, opening the door to a wide variety of precision jet substructure studies. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Ultralow Thermal Conductivity in Nanoporous Crystalline Fe 3 O 4

While there is no known fundamental lower limit to the thermal conductivity of a material, the lowest thermal conductivities are typically found in amorphous and strongly disordered materials, not highly crystalline materials. Here, we demonstrate a surprising nanostructuring route to ultralow thermal conductivity in a large-unit-cell oxide crystal (Fe 3 O 4 ) containing close-packed nanoscale pores. The electrical conductivity of this material reduces by a factor of 5 relative to dense v, independent of pore size. In contrast, thermal conductivity has a strong dependence on pore size with a factor of 40 of suppression relative to dense Fe 3 O 4 for 40 nm pores vs a factor of 5 for 500 nm pores. The matrix thermal conductivity of Fe 3 O 4 containing 40 nm pores falls below the predicted minimum thermal conductivity by a factor of 3. Finally, we attribute this to strong acoustic phonon scattering and intrinsically limited contributions to thermal conductivity from optical phonons with small dispersion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗