Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensor product”

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

A Flexible Power Method for Solving Infinite Dimensional Tensor Eigenvalue Problems

We propose a flexible power method for computing the leftmost, i.e., algebraically smallest, eigenvalue of an infinite dimensional tensor eigenvalue problem, $H x = \lambda x$, where the infinite dimensional symmetric matrix $H$ exhibits a translational invariant structure. We assume the smallest eigenvalue of $H$ is simple and apply a power iteration of $e^{-H}$ with the eigenvector represented in a compact way as a translational invariant infinite Tensor Ring (iTR). Hence, the infinite dimensional eigenvector can be represented by a finite number of iTR cores of finite rank. In order to implement this power iteration, we use a small parameter $t$ so that the infinite matrix-vector operation $e^{-Ht}x$ can efficiently be approximated by the Lie product formula, also known as Suzuki--Trotter splitting, and we employ a low rank approximation through a truncated singular value decomposition on the iTR cores in order to keep the cost of subsequent power iterations bounded. We also use an efficient way for computing the iTR Rayleigh quotient and introduce a finite size iTR residual which is used to monitor the convergence of the Rayleigh quotient and to modify the timestep $t$. In this paper, we discuss 2 different implementations of the flexible power algorithm and illustrate the automatic timestep adaption approach for several numerical examples.

Beeumen, Roel Van↗
DOE OSTI2023

Conceptual design of DIII-D experiments to diagnose the lifetime of spin polarized fuel

Abstract In magnetic fusion experiments, the cross sections for the D-T and D- 3 He fusion reactions are increased by as much as 50% if the fuel remains spin polarized parallel to the magnetic field. The goal of this study is to assess the feasibility of lifetime measurements of spin polarization, in magnetic fusion relevant conditions, on the DIII-D tokamak using relative changes in charged fusion product (CFP) loss measurements that depend upon the differential fusion cross section d σ /dΩ. Relative measurements that capture changes in the escaping CFP pitch, poloidal, and energy distributions are studied in two realistic TRANSP calculated plasma scenarios: (a) vector-polarized 3 He and D pellets are injected into a hot hydrogen plasma to produce thermonuclear reactions, and (b) a tensor-polarized deuterium pellet is injected into an L-mode hydrogen background plasma that includes unpolarized 3 He neutral beam injection. Ideal CFP signals in both scenarios show substantial pitch sensitivity to polarization for 14.7 MeV proton detection at a poloidal angle of − 56 ∘ (on the outer wall in the ion ∇ B direction), and strong sensitivity to polarization for 3.6 MeV alpha flux detection by an array of poloidal detectors. Energy-resolved measurements of 14.7 MeV protons are also sensitive to the degree of polarization for the − 56 ∘ port in the beam-plasma scenario. A realistic assessment of CFP signals in the thermonuclear scenario show count rates in the range of 2 × 10 4 cps for pitch-resolved proton detection and 2 × 10 5 cps for alpha flux measurements. Reduced chi-squared χ r 2 calculations show polarization lifetime measurements are feasible for either proton or alpha detection of both enhanced or suppressed polarization for the thermonuclear scenario. Measurements of gamma rays produced in the weak D + 3 He → γ + 5 Li branch complement the CFP measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Search for charged-lepton flavor violation in top quark production and decay in pp collisions at $ \sqrt{s} $ = 13 TeV

Results are presented from a search for charged-lepton flavor violating (CLFV) interactions in top quark production and decay in pp collisions at a center-of-mass energy of 13 TeV. The events are required to contain one oppositely charged electron-muon pair in the final state, along with at least one jet identified as originating from a bottom quark. The data correspond to an integrated luminosity of 138 fb$^{-1}$, collected by the CMS experiment at the LHC. This analysis includes both the production (q → eμt) and decay (t → eμq) modes of the top quark through CLFV interactions, with q referring to a u or c quark. These interactions are parametrized using an effective field theory approach. With no significant excess over the standard model expectation, the results are interpreted in terms of vector-, scalar-, and tensor-like CLFV four-fermion effective interactions. Finally, observed exclusion limits are set at 95% confidence levels on the respective branching fractions of a top quark to an eμ pair and an up (charm) quark of 0.13 × 10$^{-6}$ (1.31 × 10$^{-6}$), 0.07 × 10$^{-6}$ (0.89 × 10$^{-6}$), and 0.25 × 10$^{-6}$ (2.59 × 10$^{-6}$) for vector, scalar, and tensor CLFV interactions, respectively.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Semi-inclusive deep-inelastic scattering on a polarized spin-1 target. I. Cross section and spin observables

We develop the theoretical framework for semi-inclusive deep-inelastic scattering on a polarized spin-1 target and apply it to scattering on the polarized deuteron with spectator nucleon tagging. In Part I (this article) we present the general form of the semi-inclusive cross section and polarization observables for the spin-1 target. A relativistically covariant formulation in terms of four-vectors and invariant polarization parameters is employed. The target polarization is described by a spin density matrix with vector and tensor polarization. The spin and azimuthal angle dependence of the semi-inclusive cross section is derived and parametrized in terms of invariant structure functions. To validate the result, the structure functions are expressed as photon-target helicity amplitudes with known symmetry properties. The expressions presented here are kinematic (no assumptions about particle production dynamics) and valid in all regions of the deep-inelastic final state (current and target fragmentation regions). In Part II (following article), we consider deep-inelastic scattering on the polarized deuteron with spectator nucleon tagging as a special case of target fragmentation. The semi-inclusive structure functions are computed by separating nuclear and hadronic structure, and the polarization observables are explored as functions of the tagged nucleon momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Thermalization from quantum entanglement: Jet simulations in the massive Schwinger model

We investigate the emergence of thermalization in a quantum-field-theoretic model mimicking the production of jets in QCD: the massive lattice Schwinger model coupled to external sources. Specifically, we compute the expectation values of local operators as functions of time and compare them to their thermal counterparts, quantify the overlap between the evolving density matrix and the thermal one, and compare the dynamics of the energy-momentum tensor to predictions from relativistic hydrodynamics. Through these studies, we find that the system approaches thermalization at late times and elucidate the mechanisms by which quantum entanglement drives thermalization in closed field-theoretic systems. Our results show how thermodynamic behavior emerges in real time from unitary quantum dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Real-time High-resolution X-Ray Computed Tomography

Computed Tomography (CT) serves as a key imaging technology that relies on computationally intensive filtering and back-projection algorithms for 3D image reconstruction. While conventional high-resolution image reconstruction (> 2K3) solutions provide quick results, they typically treat reconstruction as an offline workload to be performed remotely on large-scale HPC systems. The growing demand for post-construction AI-driven analytics and the need for real-time adjustments call for high-resolution reconstruction solutions that are feasible on local computing resources, i.e. a multi-GPU server at most. In this paper, we propose a novel approach that utilizes Tensor Cores to optimize image reconstruction without sacrificing precision. We also introduce a framework designed to enable real-time execution of end-to-end distributed image reconstruction in a multi-GPU environment. Evaluations conducted on a single Nvidia A100 and H100 GPU show performance improvements of 1.91 × and 2.15 × compared to highly optimized production libraries. Furthermore, our framework, when deployed on 8-card Nvidia A100 GPU system, demonstrates the ability to reconstruct real-world datasets into 20483 volumes (32 GB) in slightly more than one minute and 40963 volumes (256 GB) in 7 minutes.

Wu, Du↗
DOE OSTI2022

Entropy variations and light ray operators from replica defects

We study the defect operator product expansion (OPE) of displacement operators in free and interacting conformal field theories using replica methods. We show that as n approaches 1 a contact term can emerge when the OPE contains defect operators of twist d - 2. For interacting theories and general states we give evidence that the only possibility is from the defect operator that becomes the stress tensor in the n → 1 limit. This implies that the quantum null energy condition (QNEC) is always saturated for CFTs with a twist gap. As a check, we show independently that in a large class of near vacuum states, the second variation of the entanglement entropy is given by a simple correlation function of averaged null energy operators as studied by Hofman and Maldacena. This suggests that sub-leading terms in the defect OPE are controlled by a defect version of the spin-3 non-local light ray operator and we speculate about the possible origin of such a defect operator. For free theories this contribution condenses to a contact term that leads to violations of QNEC saturation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Observation of an Axial-Vector State in the Study of the Decay ψ ( 3686 ) → ϕ η η ′

Using ( 2712.4 ± 14.3 ) × 10 6 ψ ( 3686 ) events collected with the BESIII detector at BEPCII, a partial wave analysis of the decay ψ ( 3686 ) → ϕ η η ′ is performed with the covariant tensor approach. In addition to the established states h 1 ( 1900 ) and ϕ ( 2170 ) , an axial-vector state with a mass near 2.3 GeV / c 2 is observed for the first time. Its mass and width are measured to be 2316 ± 9 stat ± 3 0 syst MeV / c 2 and 89 ± 1 5 stat ± 2 6 syst MeV , respectively. The product branching fractions of B [ ψ ( 3686 ) → X ( 2300 ) η ′ ] B [ X ( 2300 ) → ϕ η ] and B [ ψ ( 3686 ) → X ( 2300 ) η ] B [ X ( 2300 ) → ϕ η ′ ] are determined to be ( 4.8 ± 1.3 stat ± 0.7 syst ) × 10 − 6 and ( 2.2 ± 0.7 stat ± 0.7 syst ) × 10 − 6 , respectively. The branching fraction B [ ψ ( 3686 ) → ϕ η η ′ ] is measured for the first time to be ( 3.14 ± 0.1 7 stat ± 0.2 4 syst ) × 10 − 5 . The first uncertainties are statistical and the second are systematic. Published by the American Physical Society 2025

Ablikim, M.↗
DOE OSTI2023

Even Faster Exact Exchange for Solids via Tensor Hypercontraction

Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals using the Gaussian plane wave approach. This is done to lower the computational scaling with respect to the number of basis functions (N) and k-points (N k ) at a fixed system size. Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e., a set of points, N ISDF ) to further reduce computation time and memory usage. This algorithm has linear scaling cost with k-points, no explicit dependence of N ISDF on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized k-point meshes, with additional gains for large k-point meshes. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. In conclusion, we perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm.

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

Static versioning in the polyhedral model

An approach is presented to enhancing the optimization process in a polyhedral compiler by introducing compile-time versioning, i.e., the production of several versions of optimized code under varying assumptions on its run-time parameters. We illustrate this process by enabling versioning in the polyhedral processor placement pass. We propose an efficient code generation method and validate that versioning can be useful in a polyhedral compiler by performing benchmarking on a small set of deep learning layers defined for dynamically-sized tensors.

Meister, Benoit J.↗
DOE OSTI2021

Complex Fault Geometry of the 2020 M ww 6.5 Monte Cristo Range, Nevada, Earthquake Sequence

On 15 May 2020 an M ww 6.5 earthquake occurred beneath the Monte Cristo Range in the Mina Deflection region of western Nevada. Rapid deployment of eight temporary seismic stations enabled detailed analysis of its productive and slowly decaying aftershock sequence (p=0.8), which included ~18,000 autodetected events in 3.5 months. Double-difference, waveform-based relative relocation of 16,714 earthquakes reveals a complex network of faults, many of which cross the inferred 35-km-long east–northeast-striking, left-lateral mainshock rupture. Seismicity aligns with left-lateral, right-lateral, and normal mechanism moment tensors of 128 of the largest earthquakes. The mainshock occurred near the middle of the aftershock zone at the intersection of two distinct zones of seismicity. In the western section, numerous subparallel, shallow, north-northeast-striking faults form a broad flower-structure-like fault mesh that coalesces at depth into a near-vertical, left-lateral fault. We infer the near-vertical fault to be a region of significant slip in the mainshock and an eastward extension of the left-lateral Candelaria fault. Near the mainshock hypocenter, seismicity occurs on a northeast-striking, west-dipping structure that extends north from the eastern Columbus Salt Marsh normal fault. Together, these two intersecting structures bound the Columbus Salt Marsh tectonic basin. East of this intersection and the mainshock hypocenter, seismicity occurs in a narrow, near-vertical, east-northeast-striking fault zone through to its eastern terminus. At the eastern end, the aftershock zone broadens and extends northwest toward the southern extension of the northwest-striking, right-lateral Petrified Springs fault system. The eastern section hosts significantly fewer aftershocks than the western section, but has more moment release. We infer that shallow aftershocks throughout the system highlight fault-fracture meshes that connect mapped fault systems at depth. Here, comparing earthquake data with surface ruptures and a simple geodetic fault model sheds light on the complexity of this recent M 6.5 Walker Lane earthquake.

58 GEOSCIENCES↗
DOE OSTI2023

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

We present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. We inverted waveforms from 192 Global Centroid Moment Tensor earthquakes (MW 5.5-7.0) recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region. Spatial coverage of the available data is highly uneven due to earthquakes clustered along plate boundaries and sparse coverage of open seismic networks in the region. We considered three possible starting models: the SPiRaL global model (Simmons et al., 2021); MEC-1 (Kaviani et al., 2020); and CSEM2.0 (Noe et al., 2023). Because the SPiRaL model provides good fits to the observed waveforms measured by the time-bandwidth product of selected windows in several period bands, provides all the necessary parameters and covers the entire domain we used it for the starting model with the period band 50-100 seconds. Inversion iterations proceeded using time-frequency phase misfits in six stages and 54 total iterations reducing the minimum period to 30 seconds. Our final model, MESWA, provides improved waveform fits compared to the starting model for both the data used in the inversion and an independent validation data set of 66 events. Two metrics of waveform fit (the time-frequency phase misfit used in the optimization and normalized L2 misfit) were both reduced by nearly 60% for both data sets and MESWA provides significantly larger misfit reductions relative to the SPiRaL model than the MEC-1 or CSEM models. We also find that MESWA provides a larger time-bandwidth product of selected windows indicating that more information content of the observed waveforms is explained by MESWA than the other models. Our new model reveals tectonic features imaged by other studies and methods but in a new holistic model of shear and compressional wavespeeds (v S and v P , respectively) with anisotropy covering the crust and uppermost mantle of a larger domain. MESWA has smaller scale-length features and tends to sharpen some features relative to the SPiRaL starting model. Examples include: low crustal v S in the TurkishIranian Plateau, Zagros Mountains, Afghan Central Blocks and Sulaiman Fold Belt; low mantle vSfollowing divergent (Gulf of Aden, Red Sea) and transform (Dead Sea Fault) margins of the Arabian Plate; low and high v S in the mantle beneath the Arabian Shield and Platform, respectively. Low vS is imaged below Cenozoic volcanic centers of the Arabian Peninsula, the so-called Mecca-Madina-Nafud (MMN) Line. Positive anisotropy (v SH > v SV ) is inferred for asthenospheric depths across the region except where up/downwelling may influence fabric alignment (e.g. Afar, Red Sea, Arabian Shield). Elevated vS tracks Makran subduction under southeast Iran. MESWA resembles the SPiRaL model in its long-wavelength structure, but enhances shorter wavelengths features on the order of 200 km and smaller. The resulting model could be used for as a starting model for further improvements, say using waveforms from in-country seismic networks that are not openly available or smaller-scale studies targeting shorter period waveforms. The model also could be used for source characterization and moment tensor inversion to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗
DOE OSTI2023

Simulating lossy Gaussian boson sampling with matrix-product operators

Gaussian boson sampling, a computational model that is widely believed to admit quantum supremacy, has already been experimentally demonstrated and is claimed to surpass the classical simulation capabilities of even the most powerful supercomputers today. However, whether the current approach limited by photon loss and noise in such experiments prescribes a scalable path to quantum advantage is an open question. Here, to understand the effect of photon loss on the scalability of Gaussian boson sampling, we analytically derive the asymptotic operator entanglement entropy scaling, which relates to the simulation complexity. As a result, we observe that efficient tensor network simulations are likely possible under the N out ∝ √N scaling of the number of surviving photons N out in the number of input photons N. We numerically verify this result using a tensor network algorithm with U⁡(1) symmetry, and we overcome previous challenges due to the large local Hilbert-space dimensions in Gaussian boson sampling with hardware acceleration. Additionally, we observe that increasing the photon number through larger squeezing does not increase the entanglement entropy significantly. Finally, we numerically find the bond dimension necessary for fixed accuracy simulations, providing more direct evidence for the complexity of tensor networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Towards a Quantum Algorithm for the Incompressible Nonlinear Navier-Stokes Equations

In this work, we present novel concepts for quantum algorithms to solve transient, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the incompressible Navier-Stokes equations as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts solving nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. We propose a new framework based on matrix product states (MPSs) and matrix product operators (MPOs), in addition to the Krylov subspace methods. For example, the solution variables of the Navier-Stokes equations are represented by MPSs, and the linear and nonlinear terms are processed by MPOs. The time evolution of the operators is attained by a fast-forwarding algorithm using Krylov subspace methods. Furthermore, we discuss various techniques for efficient encoding of MPSs, measurement reduction for MPOs, and use of tensor operations to treat multi-variate, multi-physics characteristics of Navier-Stokes.

Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗
NASA NTRS1994

Experimental study of the flow field inside a whirling annular seal

The flow field inside a whirling annular seal has been measured using a 3-D Laser Doppler Anemometer (LDA) system. The seal investigated has a clearance of 1.27 mm, a length of 37.3 mm and is mounted on a drive shaft with a 50 percent eccentricity ratio. This results in the rotor whirling at the same speed as the shaft rotation (whirl ratio = 1.0). The seal is operated at a Reynolds number of 12000 and a Taylor number of 6300 (3600 rpm). The 3-D LDA system is equipped with a rotary encoding system which is used to produce phase averaged measurements of the entire mean velocity vector field and Reynolds stress tensor field from 0.13 mm upstream to 0.13 mm downstream of the seal. The mean velocity field reveals a highly three-dimensional flow field with large radial velocities near the inlet of the seal as well as a recirculation zone on the rotor surface. The location of maximum mean axial velocity migrates from the pressure side of the rotor at the inlet to the suction side at the exit. Turbulence production is a maximum near the seal inlet as indicated by the rapid increase of the turbulence kinetic energy (kappa). However, turbulence production and dissipation attain equilibrium fairly quickly with kappa remaining relatively constant over the last half of the seal.

Morrison, Gerald L.↗
NASA NTRS1992

Experimental study of the flow field inside a whirling annular seal

The flow field inside a whirling annular seal was measured using a 3-D Laser Doppler Anemometer (LDA) system. The seal investigated has a clearance of 1.27 mm, a length of 37.3 mm, and is mounted on a drive shaft with a 50 percent eccentricity ratio. This results in the rotor whirling at the same speed as the shaft rotation (whirl ratio = 1.0). The seal is operated at Reynolds number of 12,000 and a Taylor number of 6,300 (3,600 rpm). The 3-D LDA system is equipped with a rotary encoding system which is used to produce phase averaged measurements of the entire mean velocity vector field and Reynolds stress tensor field from 0.13 mm upstream to 0.13 mm downstream of the seal. The mean velocity field reveals a highly three dimensional flow field with large radial velocities near the inlet of the seal as well as a recirculation zone on the rotor surface. The location of maximum mean axial velocity migrates from the pressure side of the rotor at the inlet to the suction side at turbulence kinetic energy. However, turbulence production and dissipation attain equilibrium fairly quickly with remaining relatively constant over the last half of the seal.

Morrison, Gerald L.↗
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 OSTI2020

Quadrupolar interactions between acceptor pairs in $\mathcal{p}$-doped semiconductors

Here, we consider the interaction between acceptor pairs in doped semiconductors in the limit of large interacceptor separation relevant for low doping densities. Modeling individual acceptors via the spherical model of Baldereschi and Lipari, we calculate matrix elements of the quadrupole tensor between the four degenerate ground states and show that the acceptor has a nonzero quadrupole moment. As a result, the dominant contribution to the large-separation acceptor-acceptor interaction comes from direct (charge-density) terms rather than exchange terms. The quadrupole is the leading nonzero moment, so the electric quadrupole-quadrupole interaction dominates for large separation. We calculate the matrix elements of the quadrupole-quadrupole interaction Hamiltonian in a product-state basis and diagonalize, obtaining a closed-form expression for the energies and degeneracies of the sixteen-state energy spectrum. All dependence on material parameters enters via an overall prefactor, resulting in surprisingly simple and universal results. This simplicity is due, in part, to a mathematical happenstance, the nontrivial vanishing of a particular Wigner 6-$\mathcal{j}$ symbol, $^\{\frac{^2_3}{\small2} \frac{^2_3}{\small2} \frac{^2_3}{\small2}$$^\}$ = 0. Results are relevant to the control of two-qubit interactions in quantum computing implementations based on acceptor spins, as well as calculations of the thermodynamic properties of insulating $\mathcal{p}$-type semiconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗