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 703 records · Page 39

A test for interfacial effects and stress transfer in ceramic matrix composites

The efforts to fabricate single embedded filament specimens of carbon and SiC fibers were unsuccessful largely due to the thermal stresses resulting from differences in thermal coefficient of expansion. Other factors appear to have been involved including embrittlement of the metal substrate by the H2 gas in the chemical vapor deposition flow stream and reaction layers formed between the silicon carbide and the metal substrate. The carbon fiber may have been attacked by the CVD reactant. It is concluded that these differential stresses are so large as to make the embedded fiber test impractical for the study of interphase effects and stress transfer in fiber ceramic matrix systems.

Bascom, Willard D.↗

Interfaces in polymer, ceramic, and metal matrix composites; Proceedings of the Second International Conference on Composite Interfaces (ICCI-II), Cleveland, OH, June 13-17, 1988

The present conference on interfacial factors in advanced composite materials discusses such topics in their development status and properties as silane coupling, the electrochemical and plasma-surface treatment of carbon fibers, fiber surface-analytical techniques, polymer molecular scale characterization by atom-probe field-ion microscopy, the study of carbon fiber-epoxy resin interface interactions by means of labeling techniques, and the surface tailoring of SiC by ion implantation. Also discussed are developments in laser light-scattered photoelasticity, microstructural evolutions under heat treatment and radiation damage, TEM for composites, the interfaces of carbon fiber-reinforced Al composites, the influence of the interface on macroscopic composite properties, and the status of the theoretical understanding of composite interfaces.

Ishida, Hatsuo↗

Polarimetric sensitivity of light-absorbing carbonaceous aerosols over ocean: A theoretical assessment

Visible-light-absorbing carbonaceous aerosols within the boundary layer affect the radiance and polarization states of the radiation at the top of the atmosphere. Remote sensing from suborbital and satellite-based platforms utilizes these radiance and polarization signals to retrieve the key properties of these aerosols. Recent retrieval algorithms have shown a progressive trend toward including multi-angular and multi-spectral polarimetric measurements to produce better retrieval accuracy in comparison to those using measurements based on a single viewing angle. Here, we perform a theoretical investigation of the top of atmosphere (TOA) radiance-related reflectance factor (bidirectional reflectance factor (BRF)) and the two types of polarimetry-related factors (polarized bidirectional reflectance factor (pBRF) and the degree of linear polarization (DoLP)) for different types of atmospheric light-absorbing carbonaceous aerosols as a function of particle size distribution. We selected three polarimetric bands corresponding to those utilized by NASA's Airborne Multiangle SpectroPolarimetric Imager (AirMSPI)—near-UV (470 nm), visible (660 nm), and near-infrared (865 nm)—for our simulations which were performed over ocean surface using the successive order of scattering (SOS) algorithm coupled to a Lorenz-Mie aerosol optics model. The analysis of particle phase matrix elements indicates a close relationship between the angular dependencies of DoLP and associated phase matrix components at the shortest polarimetric band (470 nm). Using Jacobian analysis, we find that the radiance- and polarimetry-related reflectance factors of weakly light-absorbing aerosols, such as brown carbon, are more sensitive to changes in particle size and imaginary refractive index in comparison with those of black carbon, which is strongly light-absorbing. Our results suggest that the DoLP data could be used by future retrieval algorithms for reliably estimating microphysical properties of absorbing carbonaceous aerosols with imaginary refractive index less than 0.4.

refractive index↗

Copper-Multiwall Carbon Nanotubes and Copper-Diamond Composites for Advanced Rocket Engines

This paper reports on the research effort to improve the thermal conductivity of the copper-based alloy NARloy-Z (Cu-3 wt.%Ag-0.5 wt.% Zr), the state-of-the-art alloy used to make combustion chamber liners in regeneratively-cooled liquid rocket engines, using nanotechnology. The approach was to embed high thermal conductivity multiwall carbon nanotubes (MWCNTs) and diamond (D) particles in the NARloy-Z matrix using powder metallurgy techniques. The thermal conductivity of MWCNTs and D have been reported to be 5 to 10 times that of NARloy-Z. Hence, 10 to 20 vol. % MWCNT finely dispersed in NARloy-Z matrix could nearly double the thermal conductivity, provided there is a good thermal bond between MWCNTs and copper matrix. Quantum mechanics-based modeling showed that zirconium (Zr) in NARloy-Z should form ZrC at the MWCNT-Cu interface and provide a good thermal bond. In this study, NARloy-Z powder was blended with MWCNTs in a ball mill, and the resulting mixture was consolidated under high pressure and temperature using Field Assisted Sintering Technology (FAST). Microstructural analysis showed that the MWCNTs, which were provided as tangles of MWCNTs by the manufacturer, did not detangle well during blending and formed clumps at the prior particle boundaries. The composites made form these powders showed lower thermal conductivity than the base NARloy-Z. To eliminate the observed physical agglomeration, tangled multiwall MWCNTs were separated by acid treatment and electroless plated with a thin layer of chromium to keep them separated during further processing. Separately, the thermal conductivities of MWCNTs used in this work were measured, and the results showed very low values, a major factor in the low thermal conductivity of the composite. On the other hand, D particles embedded in NARloy-Z matrix showed much improved thermal conductivity. Elemental analysis showed migration of Zr to the NARloy-Z-D interface to form ZrC, which appeared to provide a low contact thermal resistance. These results are consistent with the quantum mechanics-based model predictions. NARloy-Z-D composites have relatively high thermal conductivities and are promising for further development.

Bhat, Biliyar N.↗

Copper Multiwall Carbon Nanotubes and Copper-Diamond Composites for Advanced Rocket Engines

This paper reports on the research effort to improve the thermal conductivity of the copper-based alloy NARloy-Z (Cu-3 wt.%Ag-0.5 wt.% Zr), the state-of-the-art alloy used to make combustion chamber liners in regeneratively-cooled liquid rocket engines, using nanotechnology. The approach was to embed high thermal conductivity multiwall carbon nanotubes (MWCNTs) and diamond (D) particles in the NARloy-Z matrix using powder metallurgy techniques. The thermal conductivity of MWCNTs and D have been reported to be 5 to 10 times that of NARloy-Z. Hence, 10 to 20 vol. % MWCNT finely dispersed in NARloy-Z matrix could nearly double the thermal conductivity, provided there is a good thermal bond between MWCNTs and copper matrix. Quantum mechanics-based modeling showed that zirconium (Zr) in NARloy-Z should form ZrC at the MWCNT-Cu interface and provide a good thermal bond. In this study, NARloy-Z powder was blended with MWCNTs in a ball mill, and the resulting mixture was consolidated under high pressure and temperature using Field Assisted Sintering Technology (FAST). Microstructural analysis showed that the MWCNTs, which were provided as tangles of MWCNTs by the manufacturer, did not detangle well during blending and formed clumps at the prior particle boundaries. The composites made form these powders showed lower thermal conductivity than the base NARloy-Z. To eliminate the observed physical agglomeration, tangled multiwall MWCNTs were separated by acid treatment and electroless plated with a thin layer of chromium to keep them separated during further processing. Separately, the thermal conductivities of MWCNTs used in this work were measured, and the results showed very low values, a major factor in the low thermal conductivity of the composite. On the other hand, D particles embedded in NARloy-Z matrix showed much improved thermal conductivity. Elemental analysis showed migration of Zr to the NARloy-Z-D interface to form ZrC, which appeared to provide a low contact thermal resistance. These results are consistent with the quantum mechanics-based model predictions. NARloy-Z-D composites have relatively high thermal conductivities and are promising for further development.

Bhat, Biliyar N.↗

A fast efficient implicit scheme for the gasdynamic equations using a matrix reduction technique

An efficient implicit finite-difference algorithm for the gasdynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved without loss of the stability characteristics 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 and used to show that the method offers computational advantages over the conventional Beam and Warming scheme. Existing Beam and Warming codes can be retrofit with minimal effort. 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.

Barth, T. J.↗

Structure and Sulfur: Tuning the Viscoelastic and Surface Properties of Natural Keratin Fibers

Natural keratin fibers, such as wool, possess a complex hierarchical structure that governs their mechanical properties and surface energy. However, the extent to which these characteristics are influenced by combined contributions of structural variations (e.g., fiber diameter, intermediate filament (IF) packing) and chemical composition (e.g., disulfide bond density) remains poorly understood. In this study, we investigate wool fibers from five sheep breeds (Merino, Polwarth, Cheviot, Eider, and Devon) to elucidate how these factors influence viscoelasticity and surface interactions. Using a multimodal approach integrating interfacial and bulk characterization methods, including inverse gas chromatography (IGC), atomic force microscopy-infrared spectroscopy (AFM-IR), X-ray photoelectron spectroscopy (XPS), uniaxial tensile testing, and synchrotron small-angle X-ray scattering (SAXS), we show that the nanometer-thick 18-methyleicosanoic acid (18-MEA) layer is consistently present across all wool types and plays a key role in governing hydrophobicity and surface heterogeneity. A controlled isothermal treatment at 200 °C, designed to cleave disulfide bonds, results in a nearly 40% reduction in specific surface area across all fiber types, accompanied by a significant decrease in tensile strength and 80% reduction in elongation at break for Merino and Devon wool, but limited influence on the mechanical properties of Eider fibers. Furthermore, rate-dependent tensile testing within the elastic regime reveals distinct viscoelastic responses among the fiber types, suggesting that the sulfur-rich protein matrix surrounding IFs and its structure contribute actively to stress partitioning. Altogether, when combined with conclusions from SAXS measurements of IF spacing, our work offers compelling insights into the role of the keratin-associated protein (KAP) matrix in shaping wool fiber mechanics. Differences in mechanical behavior among wool types, despite similar IF spacing or sulfur content, highlight the importance of matrix composition and cross-linking density, suggesting that the molecular architecture of the KAP network may be a dominant factor in determining fiber performance.

X-ray scattering↗

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Micromechanical predictions of crack initiation, propagation and crack growth resistance in boron/aluminum composites

An elastoplastic, axisymmetric finite element model has been used to predict the initiation and propagation of a crack in a composite model consisting of a single broken boron fiber embedded in an annular sheath of aluminum matrix. The accuracy of the axisymmetric finite element model for crack problems has been established by solving the classical problem of a penny-shaped crack in a thick cylindrical rod under axial tension. Also, the stress intensity factors predicted by the present numerical model are compared with continuum results. A constant displacement boundary condition applied during an increment of crack growth permits a substantial amount of stable crack growth in the matrix material. The concept of Crack Growth Resistance Curves (KR-curves) has been used to determine the point of crack instability

Mahishi, J. M.↗

Low-energy physics of isotropic spin-1 chains in the critical and Haldane phases

Using a matrix product state algorithm with infinite boundary conditions, we compute high-resolution dynamic spin and quadrupolar structure factors in the thermodynamic limit to explore the low-energy excitations of isotropic bilinear-biquadratic spin-1 chains. Haldane mapped the spin-1 Heisenberg antiferromagnet to a continuum field theory, the nonlinear sigma model (NLσM). In this work, we find that the NLσM fails to capture the influence of the biquadratic term and provides only an unsatisfactory description of the Haldane phase physics. But several features in the Haldane phase can be explained by noninteracting multimagnon states. The physics at the Uimin-Lai-Sutherland point is characterized by multisoliton continua. Moving into the extended critical phase, we find that these excitation continua contract, which we explain using a field-theoretic description. New excitations emerge at higher energies and, in the vicinity of the purely biquadratic point, they show simple cosine dispersions. Using block fidelities, we identify them as elementary one-particle excitations and relate them to the integrable Temperley-Lieb chain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spectroscopic signatures of next-nearest-neighbor hopping in the charge and spin dynamics of doped one-dimensional antiferromagnets

In this work, we study the impact of next-nearest-neighbor (nnn) hopping on the low-energy collective excitations of strongly correlated doped antiferromagnetic cuprate spin chains. Specifically, we use exact diagonalization and the density matrix renormalization group method to study the single-particle spectral function, the dynamical spin and charge structure factors, and the Cu L-edge resonant inelastic x-ray scattering (RIXS) intensity of the doped t-t'-J model for a set of t' values. We find evidence for the breakdown of spin-charge separation as |t'| increases and identify its fingerprints in the dynamical response functions. The inclusion of nnn hopping couples the spinon and holon excitations, resulting in the formation of a spin-polaron, where a ferromagnetic spin polarization cloud dresses the doped carrier. The spin-polaron manifests itself as additional spectral weight in the dynamical correlation functions, which appear simultaneously in the spin- and charge-sensitive channels. We also demonstrate that RIXS can provide a unique view of the spin-polaron, due to its sensitivity to both the spin and charge degrees of freedom.

1D spin chains↗

Brief Announcement: Communication Optimal Sparse LU Factorization for Planar Matrices

We introduce a new parallel algorithm for solving sparse LU factorization of planar matrices, which commonly arise in the finite element method for 2D PDEs. Existing scalable methods, such as the multifrontal approach with subtree-to-subcube mapping by Gupta et al. [1] and right-looking with 3D mapping by Sao et al. [2] fail to achieve optimal communication costs for these matrices. Our new algorithm combines 3D mapping and subtree-to-subcube mapping to minimize communication costs while allowing trade-offs between extra memory and reduced communication. We demonstrate that our proposed algorithm attains the communication lower bound up to a factor of O(log log n) in the memory-optimal case and up to a factor of O(log P) in the memory-independent case for an n-dimensional planar sparse matrix on P processors.

Sao, Piyush↗

Advanced Ab Initio Methods for Nuclear Structure (Final Report)

Over the past decade, there has been enormous progress in the description of nuclear structure from first principles, using interactions from Chiral Effective Field Theory that are rooted in Quantum Chromodynamics, the fundamental theory of strong interactions, and many-body methods that solve the Schrödinger equation with systematically improvable approximations. Amongst those are the family of In-Medium Similarity Renormalization Group (IMSRG) framework developed by the PI and his co-workers. While these methods scale polynomially in the size N of the single-particle basis, computational efforts still grows dramatically as we increase the degrees of freedom for the nucleons by relaxing symmetries or introducing continuum couplings (see below), or as we push to improved truncations to provide precise inputs for experimental efforts, in particular in fundamental symmetry searches. In order to address this growing computational cost, one focus area of this award was the exploration of compression and factorization methods. The key to success or failure is the presence of low-rank structures within the matrix elements of NN and 3N interactions or the IMSRG evolution operator, and a means to reformulate the method that will let us exploit them.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A high order accurate finite element algorithm for high Reynolds number flow prediction

A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.

Baker, A. J.↗

The viscoelastic behavior of the principal compliance matrix of a unidirectional graphite/epoxy composite

The time-temperature response of the principal compliances of a unidirectional graphite/epoxy composite was determined. It is shown that two components of the compliance matrix are time and temperature independent and that the compliance matrix is symmetric for the viscoelastic composite. The time-temperature superposition principle is used to determine shift factors which are independent of fiber orientation, for fiber angles that vary from 10 D to 90 D with respect to the load direction.

Morris, D. H.↗

The viscoelastic behavior of the principal compliance matrix of a unidirectional graphite/epoxy composite

A testing program was conducted to determine the time-temperature response of the principal compliances of a unidirectional graphite/epoxy composite. It was shown that the two components of the compliance matrix are time and temperature independent. In addition, the compliance matrix is symmetric for the viscoelastic composite. The time-temperature superposition principle is used to determine shift factors. It was shown that the shift factors are independent of fiber orientation for fiber angles that vary from 10 to 90 deg with respect to the load direction.

Morris, D. H.↗

Impact damage and residual tension strength of a thick graphite/epoxy rocket motor case

Impacters of various masses were dropped from various heights onto thick graphite/epoxy filament-wound cylinders. The cylinders represented filament-wound cases made for the booster motors of the Space Shuttle. Tups of various shapes were affixed to the impacters. Some of the cylinders were filled with inert propellant, and some were empty. The cylinders were impacted numerous times around the circumference and then cut into tension coupons, each containing an impact site. The size of the damage and the residual tension strength were measured. For hemispherical tups, strength was reduced as much as 30 percent by nonvisible damage. The damage consisted of matrix cracking and broken fibers. Analytical methods were used to predict the damage and residual tension strength. A factor of safety to account for nonvisible damage was determined. For corner and rod shaped tups, any damage that resulted in strength loss was readily visible.

Poe, C. C., Jr.↗