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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
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↗
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↗
DOE OSTI2023

Transverse momentum distributions of heavy hadrons and polarized heavy quarks

We initiate the study of transverse momentum-dependent (TMD) fragmentation functions for heavy quarks, demonstrate their factorization in terms of novel nonperturbative matrix elements in heavy-quark effective theory (HQET), and prove new TMD sum rules that arise from heavy-quark spin symmetry. We discuss the phenomenology of heavy-quark TMD FFs at B factories and find that the Collins effect, in contrast to claims in the literature, is not parametrically suppressed by the heavy-quark mass. We further calculate all TMD parton distribution functions for the production of heavy quarks from polarized gluons within the nucleon and use our results to demonstrate the potential of the future EIC to resolve TMD heavy-quark fragmentation in semi-inclusive DIS, complementing the planned EIC program to use heavy quarks as probes of gluon distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Computing rank‐revealing factorizations of matrices stored out‐of‐core

This paper describes efficient algorithms for computing rank-revealing factorizations of matrices that are too large to fit in main memory (RAM), and must instead be stored on slow external memory devices such as disks (out-of-core or out-of-memory). Traditional algorithms for computing rank-revealing factorizations (such as the column pivoted QR factorization and the singular value decomposition) are very communication intensive as they require many vector-vector and matrix-vector operations, which become prohibitively expensive when data is not in RAM. Randomization allows to reformulate new methods so that large contiguous blocks of the matrix are processed in bulk. The paper describes two distinct methods. The first is a blocked version of column pivoted Householder QR, organized as a “left-looking” method to minimize the number of the expensive write operations. The second method results employs a UTV factorization. It is organized as an algorithm-by-blocks to overlap computations and I/O operations. As it incorporates power iterations, it is much better at revealing the numerical rank. Numerical experiments on several computers demonstrate that the new algorithms are almost as fast when processing data stored on slow memory devices as traditional algorithms are for data stored in RAM.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Fabrication, thermal analysis, and heavy ion irradiation resistance of epoxy matrix nanocomposites loaded with silane-functionalized ceria nanoparticles

This paper describes a detailed understanding of how nanofillers function as radiation barriers within the polymer matrix, and how their effectiveness is impacted by factors such as composition, size, loading, surface chemistry, and dispersion. Here, we designed a comprehensive investigation of heavy ion irradiation resistance in epoxy matrix composites loaded with surface-modified ceria nanofillers, utilizing tandem computational and experimental methods to elucidate radiolytic damage processes and relate them to chemical and structural changes observed through thermal analysis, vibrational spectroscopy, and electron microscopy. A detailed mechanistic examination supported by FTIR spectroscopy data identified the bisphenol A moiety as a primary target for degradation reactions. Results of computational modeling by the Stopping Range of Ions in Matter (SRIM) Monte Carlo simulation were in good agreement with damage analysis from surface and cross-sectional SEM imaging. All metrics indicated that ceria nanofillers reduce the damage area in polymer nanocomposites, and that nanofiller loading and homogeneity of dispersion are key to effective damage prevention. The results of this study represent a significant pathway for engineered irradiation tolerance in a diverse array of polymer nanocomposite materials. Numerous areas of materials science can benefit from utilizing this facile and effective method to extend the reliability of polymer materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

A convergent genus expansion for the plateau

We conjecture a formula for the spectral form factor of a double-scaled matrix integral in the limit of large time, large density of states, and fixed temperature. The formula has a genus expansion with a nonzero radius of convergence. To understand the origin of this series, we compare to the semiclassical theory of “encounters” in periodic orbits. In Jackiw-Teitelboim (JT) gravity, encounters correspond to portions of the moduli space integral that mutually cancel (in the orientable case) but individually grow at low energies. At genus one we show how the full moduli space integral resolves the low energy region and gives a finite nonzero answer.

2D Gravity↗
DOE OSTI2022

Stability analysis of Monte Carlo burnup calculation based on PWR fuel rods

The existing Monte Carlo burnup algorithms may give some unstable results when using multiple burnable cells during burnup calculations. A very small time step is required to maintain the stability of the system, which pose some limits for the realistic applications. In this paper, we firstly repeat the instability issue by using a Monte Carlo code based on a typical PWR fuel rod model. Then, a theoretical stability analysis was carried out by using the perturbation theory in order to get a better understanding about the property of system stability. We determine a perturbation matrix which relates some main influence factors, such as the time step and the height of fuel rod, it can help people estimate the stability of burnup calculation system and the results are quantitatively consistent with the Monte Carlo code. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

The ITensor Software Library for Tensor Network Calculations

ITensor is a system for programming tensor network calculations with an interface modeled on tensor diagrams, allowing users to focus on the connectivity of a tensor network without manually bookkeeping tensor indices. The ITensor interface rules out common programming errors and enables rapid prototyping of algorithms. After discussing the philosophy behind the ITensor approach, we show examples of each part of the interface including Index objects, the ITensor product operator, tensor factorizations, tensor storage types, algorithms for matrix product state (MPS) and matrix product operator (MPO) tensor networks, quantum number conserving block sparse tensors, and the NDTensors library. We also review publications that have used ITensor for quantum many-body physics and for other areas where tensor networks are increasingly applied. To conclude we discuss promising features and optimizations to be added in the future.

Fishman, Matthew↗
DOE OSTI2022

Codebase release 0.3 for ITensor

ITensor is a system for programming tensor network calculations with an interface modeled on tensor diagrams, allowing users to focus on the connectivity of a tensor network without manually bookkeeping tensor indices. The ITensor interface rules out common programming errors and enables rapid prototyping of algorithms. After discussing the philosophy behind the ITensor approach, we show examples of each part of the interface including Index objects, the ITensor product operator, tensor factorizations, tensor storage types, algorithms for matrix product state (MPS) and matrix product operator (MPO) tensor networks, quantum number conserving block sparse tensors, and the NDTensors library. We also review publications that have used ITensor for quantum many-body physics and for other areas where tensor networks are increasingly applied. To conclude we discuss promising features and optimizations to be added in the future.

Fishman, Matthew↗
DOE OSTI2024

Assessment of a Grade 91 Steel Forging and Seamless Pipe Section After 141,000-Hours of Operation in a Superheat Outlet Header

In this work, two unique heats of 9Cr creep strength enhanced ferritic (CSEF) steels extracted from a retired superheat outlet header after 141,000 hours of service were evaluated. These two CSEF steels were a forging manufactured to SA-182 F91 (F91) reducer and a seamless pipe produced to SA-335 P91 (P91) pipe. Their creep deformation and fracture behavior were assessed using a lever arm creep frame integrated with in-situ high-temperature digital image correlation (DIC) system. Critical metallurgical and microstructure factors, including composition, service damage, grain matrix degradation, precipitates, and inclusions were quantitatively characterized to link the performance of the two service aged F91 and P91 CSEF steels. The creep test results show the F91 and P91 steels exhibit a large variation in creep strength and creep ductility. The F91 steel fractured at 572 hours while P91 steel fractured at 1,901 hours when subjected to a test condition of 650 °C and 100 MPa. The nominal creep strains at fracture were 12.5% (F91) and 14.5% (P91), respectively. The high-resolution DIC strain measurements reveal the local creep strain in F91 was about 50% while the local creep strain in P91 was >80%. The characterization results show that the F91 steel possessed pre-existing creep damage from its time in service, a higher fraction of inclusions, and a faster matrix grain coarsening rate. These features contribute to the observed reduction in performance for the F91 steel. The context for these findings, and the importance of metallurgical risk in an integrated life management approach will be emphasized.

Wang, Yiyu↗
DOE OSTI2021

Origin and resummation of threshold logarithms in the lattice QCD calculations of PDFs

Many present lattice QCD approaches to calculate the parton distribution functions (PDFs) rely on a factorization formula or effective theory expansion of certain Euclidean matrix elements in boosted hadron states. In the quasi- and pseudo-PDF methods, the matching coefficient in the factorization or expansion formula includes large logarithms near the threshold, which arise from the subtle interplay of collinear and soft divergences of an underlying 3D momentum distribution. We use the standard prescription to resum such logarithms in the Mellin-moment space at next-to-leading logarithmic accuracy, which also accounts for the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution, and we show that it can suppress the PDF at large x . Unlike the deep inelastic scattering and Drell-Yan cross sections, the resummation formula is away from the Landau pole. We then apply our formulation to reanalyze the recent lattice results for the pion valence PDF, and find that within the current data sensitivity, the effect of threshold resummation is marginal for the accessible moments and the PDF at large x.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Geometry-aware training of factorized layers in tensor Tucker format

Reducing parameter redundancies in neural network architectures is crucial for achieving feasible computational and memory requirements during train and inference of large networks. Given its easy implementation and flexibility, one promising approach is layer factorization, which reshapes weight tensors into a matrix format and parameterizes it as the product of two rank-r matrices. However, this family of approaches often requires an initial full-model warm-up phase, prior knowledge of a feasible rank, and it is sensitive to parameter initialization.In this work, we introduce a novel approach to train the factors of a Tucker decomposition of the weight tensors. Our training proposal proves to be optimal in locally approximating the original unfactorized dynamics and stable for the initialization. Furthermore, the rank of each mode is dynamically updated during training.We provide a theoretical analysis of the algorithm, showing convergence, approximation and local descent guarantees. The method's performance is further illustrated through a variety of experiments, showing remarkable training compression rates and comparable or even better performance than the full baseline and alternative layer factorization strategies.

Zangrando, Emanuele [Gran Sasso Science Institute ↗
DOE OSTI2021

Refinement of Salt Dissolution Inhibitor Requirements (Final Report)

At Savannah River Site (SRS), High-Level Waste is stored in below-grade carbon steel tanks. This waste in part consists of sludge, salt cake, and/or supernate. Preparation of this waste for future processing involves dissolution of the salt cake layer. The salt dissolution process can create conditions that leave the carbon steel tanks susceptible to localized corrosion. The salt to be dissolved contains high concentrations of nitrate, that once released, create an environment that may be conducive to pitting corrosion and/or stress corrosion cracking (SCC) of carbon steel. The salt dissolution process also liberates interstitial liquid trapped between the salt crystals. This liquid is initially high in nitrite and hydroxide concentration. High pH and greater ratios of nitrite to nitrate act as inhibitors to minimize corrosion of carbon steel in high nitrate environments. However, as dissolution proceeds, the concentration of nitrate will increase, while the hydroxide and nitrite concentration of the interstitial liquid will deplete and become insufficient to prevent the onset of corrosion attack. Tank blending and addition of inhibitors are used to ensure adequate concentrations of hydroxide and nitrite. However, this is not desirable during salt dissolution as it can reduce process efficiency and increase the amount of waste that needs processing. It has been proposed that these corrosion control limits be revisited to evaluate the corrosion susceptibility of carbon steel in environments that more closely resemble current operating conditions at SRS. An experimental matrix was designed to evaluate the use of the pitting factor for supernate chemistries characteristic to SRS, particularly during the salt dissolution process. Two electrochemical methods were identified to determine the susceptibility of A537 and A285 low-carbon steels to pitting corrosion with this chemistry envelope at temperatures up to 75 °C. The predominant electrochemical test method was Cyclic Potentiodynamic Polarization (CPP) studies. Through CPP, the pitting factor was used, based on Hanford Site corrosion studies, to accurately identify pitting susceptibility within the compositional range studied with some conservatism. Additionally, sulfate was determined to have no statistically significant influence, at concentrations up to 0.6 M, on pitting behavior in more concentrated solutions where other aggressive species govern pitting susceptibility. Where CPP was inconclusive, Modified ASTM G192 was successfully used to evaluate pitting susceptibility conditions and allowed for a pass/fail result to be determined. In all cases, the pitting factor was determined to be applicable to the simulants tested, with this metric accurately predicting incidences in which pitting occurred. Based upon the findings in this work, a pitting factor of 1.2 is being proposed to build in a safety factor and remain consistent with the Hanford Site approach. Additionally, a minimum pH limit of 12 is proposed to ensure carbon steel passivity and localized corrosion the primary degradation mechanism. Susceptibility to SCC was evaluated using a reduced matrix of tests at 75 °C. No failures due SCC were observed at open circuit. In addition, tests polarized anodically by 200 mV only resulted in failures for trials with pitting factors less than 0.86. However, a test with a passing condition based upon the pitting factor metric (pitting factor = 1.40) did exhibit a failure with an applied potential of +300 mV vs. OCP. This result is contrary to the prediction based upon the pitting factor, however, a polarization of 300 mV, or even 200 mV, from open circuit is substantial. The relationship between these testing parameters and service environment/conditions and the desired level of conservatism in the metric should be further evaluated in the determination of the significance of this result. While the pitting factor accurately predicted susceptibility to SCC at temperatures up to 75 °C and with positive overpotentials up to 200 mV, the relatively small sample matrix and failure of a passing pitting factor with a 300 mV polarization resulted in an inconclusive determination of whether the pitting factor may be used for predicting susceptibility to SCC at temperatures between 50 °C and 75 °C. As such, additional testing is recommended to evaluate the validity of the pitting factor for SCC susceptibility prediction at temperatures between 50 °C and 75 °C.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Measurement of Angular Coefficients of B ¯ → D * ℓ ν ¯ ℓ : Implications for | V c b | and Tests of Lepton Flavor Universality

We measure the complete set of angular coefficients J i for exclusive B ¯ → D * ℓ ν ¯ ℓ decays ( ℓ = e , μ ). Our analysis uses the full 711 fb − 1 Belle dataset with hadronic tag-side reconstruction. The results allow us to extract the form factors describing the B ¯ → D * transition and the Cabibbo-Kobayashi-Maskawa matrix element | V cb | . Using recent lattice QCD calculations for the hadronic form factors, we find | V cb | = ( 40.7 ± 0.7 ) × 10 − 3 using the Boyd-Grinstein-Lebed parametrization, compatible with determinations from inclusive semileptonic decays. We search for lepton flavor universality violation as a function of the hadronic recoil parameter w and investigate the differences of the electron and muon angular distributions. We find no deviation from standard model expectations. Published by the American Physical Society 2024

Physics↗
DOE OSTI2020

The soft $ \mathcal{S} $-matrix in gravity

The gravitational $ \mathcal{S}$-matrix defined with an infrared (IR) cutoff factorizes into hard and soft factors. The soft factor is universal and contains all the IR and collinear divergences. Here we show, in a momentum space basis, that the intricate expression for the soft factor is fully reproduced by two boundary currents, which live on the celestial sphere. The first of these is the supertranslation current, which generates spacetime supertranslations. The second is its symplectic partner, the Goldstone current for spontaneously broken supertranslations. The current algebra has an off-diagonal level structure involving the gravitational cusp anomalous dimension and the logarithm of the IR cutoff. It is further shown that the gravitational memory effect is contained as an IR safe observable within the soft $ \mathcal{S}$-matrix.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Accessing the gravitational form factors of the nucleon and nuclei through a massive graviton

In contrast to the electromagnetic form factors of the nucleon and nuclei that have been extensively studied in electron scattering, there is no known way to directly measure the gravitational form factors (GFFs), the off-forward hadronic matrix element of the QCD energy-momentum tensor. I suggest exploring the possibility to access the GFFs of the proton and nuclei in conjunction with massive graviton searches at future TeV-scale lepton-ion colliders. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Nuclear Matrix Elements for Neutrinoless Double-Beta Decay

Neutrinoless double-beta decay ($0\nu\beta\beta$) is a rare hypothesised process that, if discovered, would establish that the neutrino is Majorana, that is, it is its own antiparticle. Interpretation of experimental results relies on knowledge of nuclear matrix elements, whose large model uncertainty is the limiting factor in comparing measured (bounds on) half-lives to the neutrino mass. Nuclear effective field theory and lattice QCD have the potential to compute these matrix elements with better control over uncertainties, enhancing the discovery potential of next-generation $0\nu\beta\beta$ experiments. This work will survey various lattice QCD double-beta decay calculations and discuss their implications.

Grebe, Anthony V. [Fermilab] (ORCID:00000003103201↗