Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensor 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Explosive Yield Estimation Using Regional Seismic Moment Tensors

Here, we use the Pasyanos and Chiang (2022) data set to calculate the seismic moment M 0 for each explosion and use the measured explosive yield W to validate the W~M 0 relationship in Denny and Johnson (1991; hereafter, DJ91). The M 0 is corrected by transforming to a potency tensor and applying more appropriate near-source geophysical parameter values in the moment estimate. The mean residual between observed and predicted yield is near zero; however, the standard deviation of the residuals results in an F-value (a 95% confidence factor) of about 5. We re-estimate the coefficients in the DJ91 model and find similar values and only a slight improvement in the F-value. Next, we embark on a similar model selection process as DJ91, allowing for non-cube-root yield scaling and other plausible near-source elastic moduli. As was found by DJ91, the yield dependence is not significantly different from unity, and a cube root assumption is valid. Therefore, we yield scale the seismic moment and test the significance of all plausible explanatory variables. Isotropic moment performs better in the response variable than total moment. The preference for isotropic moment could be due to its relationship to volume change, which would be more directly affected by explosive yield. Surprisingly, we find that the overburden pressure, which is a function of depth, is not a significant parameter in the model. We hypothesize that this is due to the competing depth effects on source asymmetry and the incorporation of depth in the Green’s functions used to calculate the seismic moment tensors. Importantly, this emphasizes that only seismic moment tensor-derived moments should be used in these models. After removing insignificant model parameters, we are left with a simple model to predict explosive yield $\widehat{W}$ in kt from isotropic moment M I in N·m, $\widehat{W}$=κ –1.4132 10 0.035626GP M I , in which κ and GP are the near-source bulk modulus and gas porosity in Pa and %, respectively. The F-value for this model is approximately 3.

58 GEOSCIENCES↗
DOE OSTI2025

Multimodality in the Search for New Physics in Pulsar Timing Data and the Case of Kination-amplified Gravitational-wave Background from Inflation

We investigate the kination-amplified inflationary gravitational-wave background (GWB) interpretation of the signal recently reported by various pulsar timing array (PTA) experiments. Kination is a post-inflationary phase in the expansion history dominated by the kinetic energy of some scalar field, characterized by a stiff equation of state w = 1. Within the inflationary GWB model, we identify two modes that can fit the current data sets (NANOGrav and EPTA) with equal likelihood: the kination-amplification (KA) mode and the ordinary, no-kination-amplification (no-KA) mode. The multimodality of the likelihood motivates a Bayesian analysis with nested sampling. We analyze the free spectra of current PTA data and mock free spectra constructed with higher signal-to-noise ratios using nested sampling. The analysis of the mock spectrum designed to be consistent with the best fit to the NANOGrav 15 yr (NG15) data successfully reveals the expected bimodal posterior for the first time while excluding the reheating mode that appears in the fit to the current NG15 data, making a case for our correct and comprehensive treatment of potential multimodal posteriors arising from future PTA data sets. The resultant Bayes factor is $\mathcal{B}$ $\equiv$ Z no–KA /Z KA = 2.9 ± 1.9, indicating comparable statistical significance between the two modes. Given the theoretical model-building challenges of producing highly blue-tilted primordial tensor spectra, the KA mode has the advantage of requiring less blue primordial spectra, compared with the no-KA mode. The synergy between future cosmic microwave background polarization, pulsar timing, and laser interferometer measurements of gravitational waves will help resolve the ambiguity implied by the multimodal posterior in PTA-only searches.

Cosmology↗
DOE OSTI2020

Measurements of B-mode polarization of the cosmic microwave background from 500 square degrees of SPTpol data

We report a B-mode power spectrum measurement from the cosmic microwave background (CMB) polarization anisotropy observations made using the SPTpol instrument on the South Pole Telescope. This work uses 500 deg 2 of SPTpol data, a five-fold increase over the last SPTpol B-mode release. As a result, the bandpower uncertainties have been reduced by more than a factor of two, and the measurement extends to lower multipoles: 52 < ` < 2301. Data from both 95 and 150 GHz are used, allowing for three cross-spectra: 95 GHz x 95 GHz, 95 GHz x 150 GHz, and 150 GHz x 150 GHz. B-mode power is detected at very high significance; we find P(BB < 0) = 5.8 × 10 -71 , corresponding to a 18.1 σ detection of power. With a prior on the galactic dust from Planck, WMAP and BICEP2/Keck observations, the SPTpol B-mode data can be used to set an upper limit on the tensor-to-scalar ratio, r < 0.44 at 95% confidence (the expected 1 σ constraint on r given the measurement uncertainties is 0.22). We find the measured B-mode power is consistent with the Planck best-fit ΛCDM model predictions. Scaling the predicted lensing B-mode power in this model by a factor A lens , the data prefer A lens = 1.17 ± 0.13. These data are currently the most precise measurements of B-mode power at ℓ > 320.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Field validation of data-driven BSDF and peak extraction models for light-scattering fabric shades

Shading and daylighting systems affect cooling, heating, and lighting energy use by modulating solar radiation through the building façade. Characterizing shading systems holistically and accurately helps designers and engineers evaluate shading systems to achieve energy and non-energy performance goals. These complex fenestration systems can be modeled using Bidirectional Scattering Distribution Functions (BSDF), which map incident radiation to hemispherical distributions of outgoing radiation. Data-driven, tabulated BSDFs are derived from interpolated goniophotometer measured data, then sampled during the raytracing calculation. A peak extraction (PE) algorithm was developed to circumvent limits in BSDF angular resolution, where the specular peak is extracted during simulation by evaluating the BSDF in the through direction and surrounding region. The objective of this study was to validate this measurement and modeling workflow using field monitored data from a full scale testbed with eleven installed fabrics of different weaves, openness factors, and colors and assess the accuracy of the workflow under different adaptation and contrast conditions. Test conditions were limited to clear sky conditions with the sun in the field of view. Results showed that, for tensor tree datasets, vertical illuminance, solar luminance (2.5° apex), and daylight glare probability (DGP) were predicted to within a mean bias error (MBE) error of -456 lx (-12.3%), -3.46e5 (-38.4%), and -0.042 (-7.8%) when full PE occurred. With a binary classification of glare/ no glare, DGP was predicted accurately with a true positive rate of 0.98 and true negative rate of 1.0 using tensor tree data and less accurately with Klems BSDF data, particularly for cases of no glare. The workflow may be of insufficient accuracy to distinguish borderline performance between fabrics using the four-point glare scale, particularly under low adaptation, high contrast daylit conditions. Errors were due to reductions in peak shape and intensity across the BSDF interpolation and data reduction workflow. Future work is needed to better preserve measurement fidelity during interpolation and sampling, which in turn will improve PE performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Impact of anisotropy on TRISO fuel performance

Manufacturing of tristructural isotropic (TRISO) particles involves the deposition of pyrolytic carbon (PyC) and silicon carbide (SiC) layers using the fluidized bed chemical vapor deposition (CVD) process. The CVD process is known to generate polycrystalline layers with crystallographic textures, which imparts anisotropic thermophysical properties to the layers. Past studies have shown the risk for particle failure increases with an increase in anisotropy. The limit beyond which the anisotropy of PyC layers becomes unacceptable due to failure risk has been identified as a high-priority knowledge gap. This work presents a first systematic study on the effects of anisotropic thermal and mechanical properties on TRISO fuel performance. This computational study, performed using the fuel performance code BISON, investigates how the anisotropy in elasticity and thermal properties affect the stresses, temperature, and failure of a TRISO particle. The influence of other factors, such as operating temperature and particle geometry on the anisotropy effects, also has been analyzed. The studies utilize the recently published anisotropic elasticity and thermal behavior models for TRISO PyC and SiC layers implemented using tensors with full anisotropic capability. The spherical TRISO particles with anisotropic properties were found to have greater maximum tensile stress and significantly higher failure probability than the spherical particles with isotropic properties. In conclusion, the fuel performance predicted using these recently developed models was found to be comparable with the performance obtained using the historical models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Constraining scale dependent growth with redshift surveys

Ongoing and future redshift surveys have the capability to measure the growth rate of large scale structure at the percent level over a broad range of redshifts, tightly constraining cosmological parameters. Beyond general relativity, however, the growth rate in the linear density perturbation regime can be not only redshift dependent but scale dependent, revealing important clues to modified gravity. We demonstrate that a fully model independent approach of binning the gravitational strength G eff (k,z) matches scalar-tensor results for the growth rate fσ 8 (k,z) to 0.02%–0.27% rms accuracy. For data of the quality of the Dark Energy Spectroscopic Instrument (DESI) we find the bin values can be constrained to 1.4%–28%. We also explore the general scalar-tensor form, constraining the amplitude and past and future scalaron mass/shape parameters. Perhaps most interesting is the strong complementarity of low redshift peculiar velocity data with DESI-like redshift space distortion measurements, enabling improvements up to a factor 6–7 on 2D joint confidence contour areas. Finally, we quantify some issues with gravity parametrizations that do not include all the key physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Nuclear mass radius and pressure in the Skyrme model

We compute the mass radius, scalar radius, tensor radius, baryon number radius, and mechanical radius of nuclei with baryon number B = 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , 32, 108 in the Skyrme model. The relations between these radii and the nuclear gravitational form factors are investigated. We also compute the ‘pressure’ distribution and find that it is negative in the core region for all the nuclei with B > 1 . This suggests that the way mechanical stability is achieved in nuclei is qualitatively different than in the nucleon. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Drell-Yan angular lepton distributions at small x from TMD factorization.

The Drell-Yan process is studied in the framework of TMD factorization in the Sudakov region s » Q 2 » \( {q}_{\perp}^2 \) corresponding to recent LHC experiments with Q 2 of order of mass of Z-boson and transverse momentum of DY pair ~ few tens GeV. The DY hadronic tensors are expressed in terms of quark and quark-gluon TMDs with \( \frac{1}{Q^2} \) and \( \frac{1}{N_c^2} \) accuracy. It is demonstrated that in the leading order in N c the higher-twist quark-quark-gluon TMDs reduce to leading-twist TMDs due to QCD equation of motion. The resulting hadronic tensors depend on two leading-twist TMDs: f 1 responsible for total DY cross section, and Boer-Mulders function \( {h}_1^{\perp } \) . The corresponding qualitative and semi-quantitative predictions seem to agree with LHC data on five angular coefficients A 0 – A 4 of DY pair production. The remaining three coefficients A 5 – A 7 are determined by quark-quark-gluon TMDs multiplied by extra \( \frac{1}{N_c} \) so they appear to be relatively small in accordance with LHC results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Spatial Signatures of Electron Correlation in Least-Squares Tensor Hypercontraction

Least Squares Tensor Hypercontraction (LS-THC) has received some attention in recent years as an approach to reduce the significant computational costs of wavefunc- tion based methods in quantum chemistry. However, previous work has demonstrated that the LS-THC factorization performs disproportionately worse in the description of wavefunction components (e.g. cluster amplitudes T 2 ) than Hamiltonian compo- nents (e.g. electron repulsion integrals (pq|rs)). This work develops novel theoretical methods to study the source of these errors in the context of the real-space T 2 kernel, and reports, for the first time, the existence of a “correlation feature” in the errors of the LS-THC representation of the “exchange-like” correlation energy EX and T 2 that is remarkably consistent across ten molecular species, three correlated wavefunctions, and four basis sets. This correlation feature portends the existence of a “pair-point kernel” missing in the usual LS-THC representation of the wavefunction, which critically depends upon pairs of grid points situated close to atoms and with inter-pair distances between one and two Bohr radii. These findings point the way for future LS-THC developments to address these shortcomings.

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

Quantifying unitary flow efficiency and entanglement for many-body localization

Holography tells us that there should be a mapping between a boundary Hamiltonian and bulk degrees of freedom. We probe the bulk geometry of the Wegner Wilson Flow (WWF) in the context of many-body localization by addressing efficiency and bulk entanglement growth measures through approximating upper bounds on the boundary entanglement entropy. Here, we connect these upper bounds to the Fubini-Study metric and clarify how a central quantity, the information fluctuation complexity, distinguishes bulk unitary rotation from entanglement production. We also give a short new proof of the small incremental entangling theorem in the absence of ancillas, achieving a dimension independent, universal factor of c = 2.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Unleashing the power of EFT in neutrino-nucleus scattering

Neutrino physics is advancing into a precision era with the construction of new experiments, particularly in the few GeV energy range. Within this energy range, neutrinos exhibit diverse interactions with nucleons and nuclei. This study delves in particular into neutrino-nucleus quasi-elastic cross sections, taking into account both standard and, for the first time, non-standard interactions, all within the framework of effective field theory (EFT). The main uncertainties in these cross sections stem from uncertainties in the nucleon-level form factors, and from the approximations necessary to solve the nuclear many-body problem. We explore how these uncertainties influence the potential of neutrino experiments to probe new physics introduced by left-handed, right-handed, scalar, pseudoscalar, and tensor interactions. For some of these interactions the cross section is enhanced, making long-baseline experiments an excellent place to search for them. Our results, including tabulated cross sections for all interaction types and all neutrino flavors, can serve as the foundation for such searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The Beta-decay Paul Trap Mk IV: Design and commissioning

Here, the Beta-decay Paul Trap is an open-geometry, linear trap used to measure the decays of 8 Li and 8 B to search for a tensor contribution to the weak interaction. In the latest 8 Li measurement of Burkey et al. (2022), β scattering was the dominant experimental systematic uncertainty. The Beta-decay Paul Trap Mk IV reduces the prevalence of β scattering by a factor of 4 through a redesigned electrode geometry and the use of glassy carbon and graphite as electrode materials. The trap has been constructed and successfully commissioned with 8 Li in a new data campaign that collected 2.6 million triple coincidence events, an increase in statistics by 30% with 4 times less β scattering compared to the previous 8 Li data set.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Strong magnetoelastic coupling in Mn 3 X ( X = Ge , Sn)

Here we measure the full elastic tensors of Mn 3 Ge and Mn 3 Sn as a function of temperature through their respective antiferromagnetic phase transitions. Large discontinuities in the bulk moduli at the Néel transitions indicate strong magnetoelastic coupling in both compounds. Strikingly, the discontinuities are nearly a factor of 10 larger in Mn 3 Ge than in Mn 3 Sn. We use the magnitudes of the discontinuities to calculate the pressure derivatives of the Néel temperature, which are 39 K/GPa 14.3 K/GPa for Mn 3 Ge and Mn 3 Sn, respectively. We measured the in-plane shear modulus c 66 , which couples strongly to the magnetic order, in magnetic fields up to 18 T and found quantitatively similar behavior in both compounds. Recent measurements have demonstrated strong piezomagnetism in Mn 3 Sn : Our results suggest that Mn 3 Ge may be an even better candidate for this effect.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Design Choices in Anomaly Detection for Industrial Control Systems: Insights from Gas Pipeline Data

Industrial control systems (ICS) remain vulnerable to increasingly sophisticated cyberattacks, yet evaluating anomaly detection models in these environments is challenging due to temporal dependencies, missing-not-at-random patterns, and extremely imbalanced datasets. These factors make common practices—especially random data splits and naïve imputation—prone to severe temporal leakage, which can inflate reported performance and obscure real-world limitations. In this work, we systematically examine classical machine learning models, temporal deep learning architecture, and tensor-decomposition–based methods on a gas-pipeline dataset using a fully temporally separated evaluation pipeline designed to mimic realistic deployment conditions. Our findings show that proper temporal handling and MNAR-aware preprocessing significantly alter the relative performance of popular anomaly-detection methods, providing practical guidance for designing reliable, leakage-resistant ICS intrusion-detection systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A multisite decomposition of the tensor network path integrals

Tensor network decompositions of path integrals for simulating open quantum systems have recently been proven to be useful. However, these methods scale exponentially with the system size. This makes it challenging to simulate the non-equilibrium dynamics of extended quantum systems coupled with local dissipative environments. In this work, we extend the tensor network path integral (TNPI) framework to efficiently simulate such extended systems. The Feynman–Vernon influence functional is a popular approach used to account for the effect of environments on the dynamics of the system. In order to facilitate the incorporation of the influence functional into a multisite framework (MS-TNPI), we combine a matrix product state (MPS) decomposition of the reduced density tensor of the system along the sites with a corresponding tensor network representation of the time axis to construct an efficient 2D tensor network. The 2D MS-TNPI network, when contracted, yields the time-dependent reduced density tensor of the extended system as an MPS. The algorithm presented is independent of the system Hamiltonian. We outline an iteration scheme to take the simulation beyond the non-Markovian memory introduced by solvents. Applications to spin chains coupled to local harmonic baths are presented; we consider the Ising, XXZ, and Heisenberg models, demonstrating that the presence of local environments can often dissipate the entanglement between the sites. We discuss three factors causing the system to transition from a coherent oscillatory dynamics to a fully incoherent dynamics. The MS-TNPI method is useful for studying a variety of extended quantum systems coupled with solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗