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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Inclusive prompt photon-jet correlations as a probe of gluon saturation in electron-nucleus scattering at small x

We compute the differential cross-section for inclusive prompt photon+quark production in deeply inelastic scattering of electrons off nuclei at small x (e + A DIS) in the framework of the Color Glass Condensate effective field theory. The result is expressed as a convolution of the leading order (in the strong coupling α s ) impact factor for the process and universal dipole matrix elements, in the limit of hard photon transverse momentum relative to the nuclear saturation scale Q s,A (x). We perform a numerical study of this process for the kinematics of the Electron-Ion Collider (EIC), exploring in particular the azimuthal angle correlations between the final state photon and quark. We observe a systematic suppression and broadening pattern of the back-to-back peak in the relative azimuthal angle distribution, as the saturation scale is increased by replacing proton targets with gold nuclei. Our results suggest that photon+jet final states in inclusive e + A DIS at high energies are in general a promising channel for exploring gluon saturation that is complementary to inclusive and diffractive dijet production. They also provide a sensitive empirical test of the universality of dipole matrix elements when compared to identical measurements in proton-nucleus collisions. However because photon+jet correlations at small x in EIC kinematics require jet reconstruction at small k ⊥ , it will be important to study their feasibility relative to photon-hadron correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Sterile neutrinos and neutrinoless double beta decay in effective field theory

We investigate neutrinoless double beta decay (0νββ) in the presence of sterile neutrinos with Majorana mass terms. These gauge-singlet fields are allowed to interact with Standard-Model (SM) fields via renormalizable Yukawa couplings as well as higher-dimensional gauge-invariant operators up to dimension seven in the Standard Model Effective Field Theory extended with sterile neutrinos. At the GeV scale, we use Chiral effective field theory involving sterile neutrinos to connect the operators at the level of quarks and gluons to hadronic interactions involving pions and nucleons. This allows us to derive an expression for 0νββ rates for various isotopes in terms of phase-space factors, hadronic low-energy constants, nuclear matrix elements, the neutrino masses, and the Wilson coefficients of higher-dimensional operators. The required hadronic low-energy constants and nuclear matrix elements depend on the neutrino masses, for which we obtain interpolation formulae grounded in QCD and chiral perturbation theory that improve existing formulae that are only valid in a small regime of neutrino masses. The resulting framework can be used directly to assess the impact of 0νββ experiments on scenarios with light sterile neutrinos and should prove useful in global analyses of sterile-neutrino searches. We per- form several phenomenological studies of 0νββ in the presence of sterile neutrinos with and without higher-dimensional operators. We find that non-standard interactions involving sterile neutrinos have a dramatic impact on 0νββ phenomenology, and next-generation experiments can probe such interactions up to scales of $\mathcal O$ (100) TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Additive manufacturing of cobalt-based alloy on tool steel by directed energy deposition

Cladding hard-surfacing alloys on tool steel is an effective approach to enhance the surface properties of tool steel. In this study, a Co-based alloy was deposited on tool steel by Directed Energy Deposition (DED) following a three-factor three-level design of experiment matrix with varied laser power, scan speed, and powder flow rate. The microstructure of the deposits was characterized using scanning electron microscopy (SEM). The residual stress on the surface of the samples was measured by the X-ray diffraction (XRD) sin2? technique. The parameters that produced promising deposits were used to fabricate samples for tensile test, four-point bending test, Charpy impact test, and hardness measurement. Our result reveals that the processing parameters have a significant role in the residual stress of the coatings. Residual stress reduces with the increase of laser energy density. Cracks were found at samples with energy density below a threshold. Tensile testing of the coating/substrate combined structure reveals fracture at the coatings with an ultimate tensile strength of 633.9 ± 54.7 MPa. The bi-material interface survived the tensile test, indicating a strong interfacial bond. The four-point bending test of coating/substrate laminates shows an ultimate flexure strength of 860.6 ± 36.9 MPa. Cracks initiated from the coatings ignored the interface and penetrated the substrate, suggesting a solid bi-material bond. Charpy impact test shows the absorbed energy of coating/substrate laminates is more than doubled that of the substrate.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ultrafine Interwoven Dendritic Cu 2 Se/CuFeSe 2 Composites with Enhanced Thermoelectric Performance

As an effort to develop high-performance and stable copper chalcogenide thermoelectrics, several (1-x)Cu 2 Se/(x)CuFeSe 2 composites were successfully synthesized via a one-step solid-state transformation of CuSe 2 template. Uniform ultrafine dendritic structures consisting of interweaving Cu 2 Se and CuFeSe 2 nanofibers were observed in samples with low Fe content (0.05 ≤ x ≤ 0.1). Increasing the Fe content to x = 0.5 led to phase segregation into Cu 2 Se-rich region with embedded CuFeSe 2 fine structures interwoven with CuFeSe 2 -rich region containing Cu 2 Se nanofibers. The formation of such an entwined dendritic structure is believed to arise from the temperature-dependent solubility of CuFeSe 2 in the Cu 2 Se matrix. The dynamic dissolution of CuFeSe 2 into the Cu 2 Se at high temperatures leads to temperature-dependent doping of the Cu 2 Se matrix, enabling drastic enhancements of thermoelectric power factor at high temperatures. As a result, such a strategy is expected to be a powerful tool for properties modulation in thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Signatures of the attractive interaction in spin spectra of one-dimensional cuprate chains

Identifying the minimal model for cuprates is crucial for explaining the high- T c pairing mechanism. Recent photoemission experiments have suggested a significant near-neighbor attractive interaction V in cuprate chains, favoring pairing instability. To determine its strength, we systematically investigate the dynamical spin structure factors S ( q , ω ) using the density matrix renormalization group. Our analysis quantitatively reveals a notable softening in the two-spinon continuum, particularly evident in the intense spectrum at large momentum. This softening is primarily driven by the renormalization of the superexchange interaction, as determined by a comparison with the slave-boson theory. We also demonstrate the feasibility of detecting this spectral shift in thin-film samples using resonant inelastic x-ray scattering. Therefore, this provides a distinctive fingerprint for the attractive interaction, motivating future experiments to unveil essential ingredients in cuprates. Published by the American Physical Society 2024

Shen, Zecheng↗
DOE OSTI2022

The TPR domain of PgaA is a multifunctional scaffold that binds PNAG and modulates PgaB-dependent polymer processing

The synthesis of exopolysaccharides as biofilm matrix components by pathogens is a crucial factor for chronic infections and antibiotic resistance. Many periplasmic proteins involved in polymer processing and secretion in Gram-negative synthase dependent exopolysaccharide biosynthetic systems have been individually characterized. The operons responsible for the production of PNAG, alginate, cellulose and the Pel polysaccharide each contain a gene that encodes an outer membrane associated tetratricopeptide repeat (TPR) domain containing protein. While the TPR domain has been shown to bind other periplasmic proteins, the functional consequences of these interactions for the polymer remain poorly understood. Herein, we show that the C-terminal TPR region of PgaA interacts with the de-N-acetylase domain of PgaB, and increases its deacetylase activity. Additionally, we found that when the two proteins form a complex, the glycoside hydrolase activity of PgaB is also increased. To better understand structure-function relationships we determined the crystal structure of a stable TPR module, which has a conserved groove formed by three repeat motifs. Tryptophan quenching, mass spectrometry analysis and molecular dynamics simulation studies suggest that the crystallized TPR module can bind PNAG/dPNAG via its electronegative groove on the concave surface, and potentially guide the polymer through the periplasm towards the porin for export. Our results suggest a scaffolding role for the TPR domain that combines PNAG/dPNAG translocation with the modulation of its chemical structure by PgaB.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Analysis of Covariance Intersection For Triangulation

In this document I will discuss different implementations of Covariance Intersection (CI) for object triangulation as well as the robustness of CI. For CI methods we will compare the performance of different methods and discuss edge cases which must be considered. For robustness we will focus on the impact of removing sensors on the final fused estimate. Here we look at factors which influence the final fused covariance matrix.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Binning is Sinning: Redemption for Hubble Diagram Using Photometrically Classified Type Ia Supernovae

Bayesian Estimation Applied to Multiple Species (BEAMS) is implemented in the BEAMS with Bias Corrections (BBC) framework to produce a redshift-binned Hubble diagram (HD) for Type Ia supernovae (SNe Ia). BBC corrects for selection effects and non–SN Ia contamination, and systematic uncertainties are described by a covariance matrix with dimension matching the number of BBC redshift bins. For spectroscopically confirmed SN Ia samples, a recent "Binning is Sinning" article showed that an unbinned HD and covariance matrix reduces the systematic uncertainty by a factor of ~1.5 compared to the binned approach. Here we extend their analysis to obtain an unbinned HD for a photometrically identified sample processed with BBC. To test this new method, we simulate and analyze 50 samples corresponding to the Dark Energy Survey (DES) with a low-redshift anchor; the simulation includes SNe Ia, and contaminants from core-collapse SNe and peculiar SNe Ia. The analysis includes systematic uncertainties for calibration and measures the dark energy equation of state parameter (w). Compared to a redshift-binned HD, the unbinned HD with nearly 2000 events results in a smaller systematic uncertainty, in qualitative agreement with BHS21, and averaging results among the 50 samples we find no evidence for a w-bias. To reduce computation time for fitting an unbinned HD with large samples, we propose an HD-rebinning method that defines the HD in bins of redshift, color, and stretch; the rebinned HD results in similar uncertainty as the unbinned case, and shows no evidence for a w-bias.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Scaled ILU Smoothers for Navier-Stokes Pressure Projection

Incomplete LU (ILU) smoothers are effective in the algebraic multigrid (AMG) V-cycle for reducing high-frequency components of the error. However, the requisite direct triangular solves are comparatively slow on GPUs. Previous work has demonstrated the advantages of Jacobi iteration as an alternative to direct solution of these systems. Depending on the threshold and fill-level parameters chosen, the factors can be highly nonnormal and Jacobi is unlikely to converge in a low number of iterations. We demonstrate that row scaling can reduce the departure from normality, allowing us to replace the inherently sequential solve with a rapidly converging Richardson iteration. There are several advantages beyond the lower compute time. Scaling is performed locally for a diagonal block of the global matrix because it is applied directly to the factor. Further, an ILUT Schur complement smoother maintains a constant GMRES iteration count as the number of MPI ranks increases, and thus parallel strong-scaling is improved. Our algorithms have been incorporated into hypre, and we demonstrate improved time to solution for linear systems arising in the Nalu-Wind and PeleLM pressure solvers. For large problem sizes, GMRES+AMG executes at least five times faster when using iterative triangular solves compared with direct solves on massively parallel GPUs.

algebraic multigrid↗
DOE OSTI2025

Towards excitations and dynamical quantities in correlated lattices with density matrix embedding theory

Density matrix embedding theory (DMET) provides a framework to describe ground-state expectation values in strongly correlated systems, but its extension to dynamical quantities is still an open problem. We show one route to obtaining excitations and dynamical spectral functions by using the techniques of DMET to approximate the matrix elements that arise in a single-mode inspired excitation ansatz. We demonstrate this approach in the one-dimensional Hubbard model, comparing the neutral excitations, single-particle density of states, charge, and spin dynamical structure factors to benchmarks from the Bethe ansatz and density matrix renormalization group. Finally, our work highlights the potential of these ideas in building computationally efficient approaches for dynamical quantities.

1-dimensional systems↗
DOE OSTI2021

Analysis Report for Determination of pH Correction Factors for Brines under AP-157 Rev.1

This Analysis Report (AR) documents the determination of pH correction factors for the observed pH readings. The correction factor converts the observed pH reading recorded from the brines used in geochemical studies in support of the Waste Isolation Pilot Plant (WIPP) to a corrected pH value. The data analysis in this AR falls under AP-157 Rev. 1 Analysis Plan for Determination of pH Correction Factors in Brines (Kirkes et al, 2021). Measurement of pH in some solutions can be challenging due to numerous factors such as high ionic strength, elevated or lowered temperature, complex matrix composition, etc. (Knauss et al., 1990; and Rai et al., 1995). The measured pH can be corrected by applying the correction factor, empirically obtained from a specific test solution.

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

A new upper bound for the growth factor in Gaussian elimination with complete pivoting

Abstract The growth factor in Gaussian elimination measures how large the entries of an LU factorization can be relative to the entries of the original matrix. It is a key parameter in error estimates, and one of the most fundamental topics in numerical analysis. We produce an upper bound of for the growth factor in Gaussian elimination with complete pivoting — the first improvement upon Wilkinson's original 1961 bound of .

Bisain, Ankit [Department of Mathematics Massachus↗
DOE OSTI2026

Extracting $B_s\to D_s^*\ellν_\ell$ form factors

Semileptonic $B_{(s)}$ decays are of great phenomenological interest because they allow to determine e.g. CKM matrix elements or test lepton flavor universality. Taking advantage of already existing lattice data, we demonstrate the analysis steps to extract the four form factors describing exclusive semileptonic $B_s\to D_s^*\ellν_\ell$ decays using the narrow width approximation. Our data are based on RBC/UKQCD's set of 2+1 flavor gauge field ensembles with Shamir domain-wall fermion and Iwasaki gauge field action featuring inverse lattice spacings of $a^{-1}=$1.785, 2.383, and 2.785 GeV as well as pion masses between 268 and 433 MeV. Light, strange and charm quarks are simulated using domain-wall fermions, whereas bottom quarks are generated with the relativistic heavy quark (RHQ) action.

Boushmelev, Anastasia [Siegen U.] (ORCID:000000033↗
DOE OSTI2023

Tangent functional connectomes uncover more unique phenotypic traits

Functional connectomes (FCs) containing pairwise estimations of functional couplings between pairs of brain regions are commonly represented by correlation matrices. As symmetric positive definite matrices, FCs can be transformed via tangent space projections, resulting into tangent-FCs. Tangent-FCs have led to more accurate models predicting brain conditions or aging. Motivated by the fact that tangent-FCs seem to be better biomarkers than FCs, we hypothesized that tangent-FCs have also a higher fingerprint. We explored the effects of six factors: fMRI condition, scan length, parcellation granularity, reference matrix, main-diagonal regularization, and distance metric. Our results showed that identification rates are systematically higher when using tangent-FCs across the “fingerprint gradient” (here including test-retest, monozygotic and dizygotic twins). Highest identification rates were achieved when minimally (0.01) regularizing FCs while performing tangent space projection using Riemann reference matrix and using correlation distance to compare the resulting tangent-FCs. Such configuration was validated in a second dataset (resting-state).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Active operator learning with predictive uncertainty quantification for partial differential equations

With the increased prevalence of neural operators being used to provide rapid solutions to partial differential equations (PDEs), understanding the accuracy of model predictions and the associated error levels is necessary for deploying reliable surrogate models in scientific applications. Existing uncertainty quantification (UQ) frameworks employ ensembles or Bayesian methods, which can incur substantial computational costs during both training and inference. Here, we propose a lightweight predictive UQ method tailored for Deep operator networks (DeepONets) that also generalizes to other operator networks. Numerical experiments on linear and nonlinear PDEs demonstrate that the framework’s uncertainty estimates are unbiased and provide accurate out-of-distribution uncertainty predictions with a sufficiently large training dataset. Our framework provides fast inference and uncertainty estimates that can efficiently drive outer-loop analyses that would be prohibitively expensive with conventional solvers. We demonstrate how predictive uncertainties can be used in the context of Bayesian optimization and active learning problems to yield improvements in accuracy and data-efficiency for outer-loop optimization procedures. In the active learning setup, we extend the framework to Fourier Neural Operators (FNO) and describe a generalized method for other operator networks. To enable real-time deployment, we introduce an inference strategy based on precomputed trunk outputs and a sparse placement matrix, reducing evaluation time by more than a factor of five. Our method provides a practical route to uncertainty-aware operator learning in time-sensitive settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A matrix-free approach for finite-strain hyperelastic problems using geometric multigrid

This paper investigates matrix-free algorithms for problems in quasi-static finite-strain hyperelasticity. Iterative solvers with matrix-free operator evaluation have emerged as an attractive alternative to sparse matrices in the fluid dynamics and wave propagation communities because they significantly reduce the memory traffic, the limiting factor in classical finite element solvers. Specifically, we study different matrix-free realizations of the finite element tangent operator and determine whether generalized methods of incorporating complex constitutive behavior might be feasible. Furthermore, in order to improve the convergence behavior of iterative solvers, we also propose a method by which to construct level tangent operators and employ them to define a geometric multigrid preconditioner. Additionally, the performance of the matrix-free operator and the geometric multigrid preconditioner is compared to the matrix-based implementation with an algebraic multigrid (AMG) preconditioner on a single node for a representative numerical example of a heterogeneous hyperelastic material in two and three dimensions. We find that matrix-free methods for finite-strain solid mechanics are very promising, outperforming linear matrix-based schemes by two to five times, and that it is possible to develop numerically efficient implementations that are independent of the hyperelastic constitutive law.

42 ENGINEERING↗
DOE OSTI2025

Form factors for semi-leptonic $B_{(s)} \to D_{(s)}^∗\ell \nu_\ell$ decays

Semileptonic $B_{(s)}$ decays are of great phenomenological interestbecause they allow to extract CKM matrix elements or test lepton flavouruniversality. Taking advantage of existing data, we explore extractingform factors for vector final states using the narrow widthapproximation. Based on RBC/UKQCD's set of 2+1 flavour gauge fieldensembles with Shamir domain-wall fermion and Iwasaki gauge fieldaction, we study semileptonic $B_{(s)}$ decays using domain-wallfermions for light, strange and charm quarks, whereas bottom quarks aresimulated with the relativistic heavy quark (RHQ) action. Exploratoryresults for $B_s \to D_s^* \ell \nu_\ell$ are presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Pass-efficient methods for compression of high-dimensional turbulent flow data

The future of high-performance computing, specifically on future Exascale computers, will presumably see memory capacity and bandwidth fail to keep pace with data generated, for instance, from massively parallel partial differential equation (PDE) systems. Current strategies proposed to address this bottleneck entail the omission of large fractions of data, as well as the incorporation of in situ compression algorithms to avoid overuse of memory. To ensure that post-processing operations are successful, this must be done in a way that a sufficiently accurate representation of the solution is stored. Moreover, in situations where the input/output system becomes a bottleneck in analysis, visualization, etc., or the execution of the PDE solver is expensive, the number of passes made over the data must be minimized. In the interest of addressing this problem, this work focuses on the utility of pass-efficient, parallelizable, low-rank, matrix decomposition methods in compressing high-dimensional simulation data from turbulent flows. Additionally, a particular emphasis is placed on using coarse representation of the data – compatible with the PDE discretization grid – to accelerate the construction of the low-rank factorization. This includes the presentation of a novel single-pass matrix decomposition algorithm for computing the so-called interpolative decomposition. The methods are described extensively and numerical experiments on two turbulent channel flow data are performed. In the first (unladen) channel flow case, compression factors exceeding 400 are achieved while maintaining accuracy with respect to first- and second-order flow statistics. In the particle-laden case, compression factors of 100 are achieved and the compressed data is used to recover particle velocities. These results show that these compression methods can enable efficient computation of various quantities of interest in both the carrier and disperse phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗