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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Getting to the Core of PARAFAC2, A Nonnegative Approach

In this paper, the authors present a novel method of performing PARAFAC2 factorization of three-way data using a compact representation of that data. In the standard PARAFAC2 algorithm, two modes of the data are recovered directly during the decomposition while the third mode is returned as a transformation matrix, which is then used to rotate sets of orthogonal third-mode basis factors into interpretable factors. In our new method, the data are first decomposed into a core matrix and orthogonal factor loading matrices in the first two modes as well as sets of orthogonal factors in the third mode (as in standard PARAFAC2). The core matrix is then decomposed using a the standard PARAFAC2 strategy to produce transformation matrices in all three modes. The algorithm is particularly useful for very large data sets and essentially permits imposition of nonnegativity in all three modes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Spectroscopy of 13 Be through isobaric analog states in 13 B

Spectroscopy of the exotic, neutron-unbound beryllium isotope, 13 Be, is still a puzzle despite significant experimental efforts. To observe T = 5/2 states in 13 B, establish spin parities and spectroscopic factors, and use isospin symmetry to inform spectroscopy of 13 Be, excitation functions for resonance elastic scattering of 12 Be on protons were measured in the c.m. energy range from 1 MeV to 5 MeV. Two T = 5/2 states in 13 B at excitation energies 18.25 MeV and 19.95 MeV were observed. Unambiguous 1/2 + and 5/2 + spin-parity assignments were made, and spectroscopic factors were established using R-matrix analysis of the measured 12 Be + p excitation functions. Here, we provide the first direct and unambiguous spin-parity assignments and spectroscopic factor measurements for the resonances in the A = 13, T = 5/2 isobaric multiplet. This puts us on a solid footing to identify the 2.3 MeV state in 13 Be as 5/2 + and infer the existence of an s-wave (1/2 + ) resonance at around 0.6 MeV with a relatively small spectroscopic factor.

6 ≤ A ≤ 19↗
DOE OSTI2020

Efficient Low-Order Refined Preconditioners for High-Order Matrix-Free Continuous and Discontinuous Galerkin Methods

In this paper, we design preconditioners for the matrix-free solution of high-order continuous and discontinuous Galerkin discretizations of elliptic problems based on finite element method--spectral element method (FEM-SEM) equivalence and additive Schwarz methods. The high-order operators are applied without forming the system matrix, making use of sum factorization for efficient evaluation. The system is preconditioned using a spectrally equivalent low-order (p = 1) finite element operator discretization on a refined mesh. The low-order refined mesh is anisotropic and not shape regular in the polynomial degree of the high-order operator, requiring specialized solvers to treat the anisotropy. We make use of an element-structured, geometric multigrid V-cycle with ordered ILU(0) smoothing. The preconditioner is parallelized through an overlapping additive Schwarz method that is robust in h and p. The method is extended to interior penalty and Bassi and Rebay (BR2) discontinuous Galerkin discretizations, for which it is also robust in the size of the penalty parameter. Finally, numerical results are presented on a variety of examples, verifying the uniformity of the preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Composition Quantification of SiGeSn Alloys Through Time-of-Flight Secondary Ion Mass Spectrometry: Calibration Methodologies and Validation With Atom Probe Tomography

Time-of-Flight Secondary Ion Mass Spectrometry (ToF-SIMS) is a powerful technique for elemental compositional analysis and depth profiling of materials. However, it encounters the problem of matrix effects that hinder its application. In this work, we introduce a pioneering ToF-SIMS calibration method tailored for SixGeySnz ternary alloys. SixGe1-x and Ge1-zSnz binary alloys with known compositions are used as calibration reference samples. Through a systematic SIMS quantification study of SiGe and GeSn binary alloys, we unveil a linear correlation between secondary ion intensity ratio and composition ratio for both SiGe and GeSn binary alloys, effectively mitigating the matrix effects. Extracted relative sensitivity factor (RSF) value from SixGe1-x (0.07 < x < 0.83) and Ge1-zSnz (0.066 < z < 0.183) binary alloys are subsequently applied to those of SixGeySnz (0.011 < x < 0.113, 0.863 < y < 0.935 and 0.023 < z < 0.103) ternary alloys for elemental compositions quantification. These values are cross-checked by Atom Probe Tomography (APT) analysis, an indication of the great accuracy and reliability of as-developed ToF-SIMS calibration process. Furthermore, the proposed method and its reference sample selection strategy in this work provide a low-cost as well as simple-to-follow calibration route for SiGeSn composition analysis, thus driving the development of next-generation multifunctional SiGeSn-related semiconductor devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Neutral excitations of quantum Hall states: A density matrix renormalization group study

We use the dynamical structure factors of the quantum Hall states at ν = 1 / 3 and 1 / 2 in the lowest Landau level to study their excitation spectrum. Using the density matrix renormalization group in combination with the time-dependent variational principle on an infinite cylinder geometry, we extract the low-energy properties. At ν = 1 / 3 , a sharp magnetoroton mode and the two-roton continuum are present and the finite-size effects can be understood using the fractional charge of the quasiparticle. At ν = 1 / 2 , we find low-energy modes with linear dispersion and the static structure factor s ¯ ( q ) ~ ( q ℓ ) 3 in the limit q ℓ → 0 . The properties of these modes agree quantitatively with the predictions of the composite-fermion theory placed on the infinite cylinder.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Independence of Environmental and Mechanical Damages on Silicone Adhesive Stored in a Thermo-Oxidative Environment

Abstract Elastomeric polymer materials hold a special place in today’s industrial sectors with respect to structural application needs. Low quality materials are detrimental to industries like aerospace, and automotive engineering. Assessment of the durability of elastomers for structural applications has been of much interest among the literary circles for decades and new materials keep outperforming the existing ones. Polymeric adhesives are one of the most abundantly used materials in these industries. All polymeric materials get damaged when in contact with aggressive environments in presence of high temperature and oxygen. Commonly referred to as thermo-oxidation, this environment exposes the material to heat and oxidation reactions in presence of oxygen. Resultantly, during service life, the damage to the polymer matrix is primarily caused by two factors: mechanical damage and environmental aging. Environmental aging is an irreversible phenomenon caused by changes in the molecular structure while mechanical damage maintains the shape of the polymer matrix, and the deterioration is mostly due to polymer chain mobility. Environmental aging can be caused by a single environmental agent or by a synergized impact of several environmental elements. Increasing temperature is found to be proportional to decreasing tensile strength and toughness of material. The rate and extent of degradation can be accessed by scrutinizing the changes in constitutive behavior of material through mechanical and chemical properties. Accelerated thermal aging is among the most common modes of process related degradation, leading generally to chain scission, and crosslinking phenomena and reduced resistance to fracture stress, and strain. In this experimental study, our goal is to separate the environmental degradation from mechanical damage. A silicone-based adhesive was aged in thermo-oxidative at (0%RH) aging environment. The damage and decay mechanisms have been used to draw a distinction between environmental degradation and mechanical damage. Material characterization included uniaxial tensile test (failure and cyclic) and scanning electron microscopy (SEM) tests on as-received and aged samples. Aging was conducted at three different temperatures (60°C, 80°C and 95°C) and six different exposure durations (1, 3, 10, 30, 90 and 200 days). This work confirms that environmental damage is superposed on top of the mechanical damage, and thus, they are separable. The total mechanical and environmental damage is a synergized effect of all exposure conditions and parameters involved i.e., aging time, temperature, and oxygen. The chemistry and mechanics of the polymer degradation were found to be in good agreement with each other.

Alazhary, Sharif↗
DOE OSTI2025

On-shell recursion and holomorphic HQET for heavy quark hadronic resonances

We develop a new theoretical framework for the treatment of heavy quark (HQ) resonances within heavy quark effective theory (HQET). This framework uses on-shell recursion techniques to express the resonant amplitude as a product of on-shell subamplitudes, which allows one to employ a form-factor representation of the hadronic matrix elements and to obtain an HQ expansion, but at the price of introducing complex momenta. We construct a generalized “holomorphic HQET” onto which such complex-momentum matrix elements can be matched, and we show that PT symmetry ensures the Isgur-Wise functions (and the perturbative corrections) become holomorphic functions of the complex recoil parameter with real coefficients. They are thus an analytic continuation of the standard HQET description. This framework admits a HQ hadron (strong decay) width expansion. At second order, we show it is compatible with data for the B12∗$$ {B}_{1(2)}^{\left(\ast \right)} $$ and D12∗$$ {D}_{1(2)}^{\left(\ast \right)} $$ HQ doublets. Taking the B¯→D1∗1−→Dπlν$$ \overline{B}\to \left({D}_1^{\ast}\left({1}^{-}\right)\to D\pi \right) l u $$ system as an example, we compute the holomorphic HQET expansion to first order, as well as the complex-momentum on-shell subamplitudes. A toy numerical study of the resulting differential rates demonstrates that this framework generates HQ resonance lineshapes with large tails, resembling those seen in data.

Manzari, Claudio Andrea↗
DOE OSTI2021

Learning constitutive relations using symmetric positive definite neural networks

In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Partial inhibition of activin receptor-like kinase 4 alleviates bladder fibrosis caused by bladder outlet obstruction

The bladder undergoes profound structural alterations after bladder outlet obstruction (BOO), characterized by hypertrophy of the bladder wall and accumulation of extracellular matrix (ECM). Transforming growth factor-β (TGF-β) has been found to promote fibrosis of the bladder induced by partial bladder outlet obstruction (pBOO). Activin receptor-like kinase 4 (ALK4) is a downstream receptor of the TGF-β superfamily. However, the role of the ALK4–Smad2/3 pathway in the pathogenesis of bladder fibrosis caused by pBOO remains unknown. This study focused on learning the role of ALK4 in the process of bladder fibrosis caused by pBOO. The pBOO mice models showed that ALK4 expression was found to upregulate in the wild-type bladder 6 weeks after pBOO compared to control group. Then, mice with heterozygous knockout of the ALK4 gene (ALK4+/−) were generated. Histological analysis and Western blot (WB) results showed significant suppression of collagen expression in the bladders of ALK4+/− mice after pBOO compared with WT mice. WB also showed that ALK4+/− mice demonstrated significant suppression of phosphorylated Smad2/3 (p-Smad2/3) expression in the bladder 6 weeks after pBOO but not of phosphorylated extracellular signal-regulated kinase, c-Jun N-terminal kinase or protein 38 (p-ERK, p-JNK, p-P38) expression. This effect might have occurred through partial inactivation of the Smad2/3 signaling pathway. In vitro, ALK4 overexpression promoted collagen production in cultured BSMCs and activated the Smad2/3 signaling pathway. Taken together, our results demonstrated that ALK4 insufficiency alleviated bladder fibrosis in a mouse model of pBOO partly by suppressing Smad2/3 activity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Dimerization and spin decoupling in a two-leg Heisenberg ladder with frustrated trimer rungs

We study the antiferromagnetic spin-half Heisenberg ladder in the presence of an additional frustrating rung spin which is motivated and relevant also for the description of real two-dimensional materials such as the two-dimensional trimer magnet Ba 4 Ir 3 O 10 . We study the zero-temperature phase diagram, where we combine numerical and analytical methods into an overall consistent description. All numerical simulations are also accompanied by studies of the dynamical spin structure factor obtained via the density matrix renormalization group. Overall, we find in the regime of strong rung coupling a gapped dimerized phase related to competing symmetry sectors in Hilbert space that ultimately results in frustration-driven spin-Peierls transition. In the weak rung-coupling regime, the system is uniform, yet shows a gapped spinon continuum together with a sharp coherent low-energy branch which renders the system critical overall. In either case, the additional rung spin quickly get sidelined and nearly decouple in the regime when their bare coupling to the ladder is weakened

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

IRIS-DMEM: Efficient Memory Management for Heterogeneous Computing

This paper proposes an efficient data memory management approach for the Intelligent RuntIme System (IRIS) heterogeneous computing framework along with new data transfer policies. IRIS provides a task-based programming model for extreme heterogeneous computing (e.g., CPU, GPU, DSP, FPGA) with support for today's most important programming languages (e.g., OpenMP, OpenCL, CUDA, HIP, OpenACC). However, the IRIS framework either forces the programmer to introduce data transfer commands for each task or relies on suboptimal memory management for automatic and transparent data transfers. The work described here extends IRIS with novel heterogeneous memory handling and introduces novel data transfer policies by employing the Distributed data MEMory handler (DMEM) for efficient and optimal movement of data among the various computing resources. The proposed approach achieves performance gains of up to 7× for tiled LU factorization and tiled DGEMM (i.e., matrix multiplication) benchmarks. Moreover, this approach also reduces data transfers by up to 71% when compared to previous IRIS heterogeneous memory management handlers. This work compares the performance results of the IRIS framework's novel DMEM with the StarPU runtime and MAGMA math library for GPUs. Experiments show a performance gain of up to 1.95× over StarPU and 2.1× over MAGMA.

Miniskar, Narasinga Rao↗
DOE OSTI2020

Tucker-1 Boolean Tensor Factorization with Quantum Annealers

Quantum annealers are an emerging computational architecture that have the potential to address some challenging computational issues that will be left unresolved as we approach the end of the Moore's Law era of computing. D-Wave quantum annealers are designed to solve a challenging set of problems - quadratic unconstrained binary optimization problems. This makes them a natural fit for solving problems with binary or Boolean variables. Here, we explore the use of a quantum annealer to solve Boolean tensor factorization. The goal of Boolean tensor factorization is to represent a high-dimensional tensor filled with Boolean values as a product of Boolean matrices and a Boolean core tensor. We show that a particular Boolean tensor factorization problem (called Tucker-1 factorization) can be decomposed into a sequence of quadratic unconstrained binary optimization problems that can be solved with a D-Wave 2000Q quantum annealer. While quantum annealers specifically and quantum computers in general are at a fairly early stage in their development, they are currently capable of solving these Boolean tensor factorization problems. Importantly, our results show that for fairly small tensors, we are frequently able to obtain an accurate (sometimes exact) factorization using quantum annealing.

97 MATHEMATICS AND COMPUTING↗
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↗
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↗
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↗
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↗