Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensor”

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

Probing the Electronic Structure of a Thorium Nitride Complex by Solid-State 15N NMR Spectroscopy

The solid-state 15N NMR powder spectra of the thorium nitride complex, [K(18-crown-6)(THF)2][(R2N)3Th(m-15N)(Th(NR2)3] ([K][1-15N], R = SiMe3), and the thorium amide complex, [Th(NR2)3(15NH2)] (2-15N) were recorded. The spectrum for [K][1-15N] represents the first reported solid-state 15N NMR data for an actinide nitride complex. The experimentally measured tensor spans are 847 ppm for [K][1-15N] and 237 ppm for 2-15N. Both shielding tensors exhibit a near-zero asymmetry parameter, which for [K][1-15N] is consistent with a local rotational symmetry of its 15N-labelled nitride ligand. For 2-15N, the lack of asymmetry can be rationalized by a quasi-free Th-NH2 bond rotation in the solid-state. DFT calculations overestimate the tensor span somewhat for [K][1-15N], but provide isotropic shifts in good agreement with both the solid-state and solution values for both complexes. Natural localized molecular orbital (NLMO) analyses of the nuclear shielding reveal that the larger tensor span in [K][1-15N] vs. 2-15N is primarily a consequence of more pronounced covalency of the N-Th bonds, and large spin-orbit coupling due to significant Th 5f orbital contribution to those bonds, impacting the principal components of the shielding tensor perpendicular to the Th-N-Th axis. Overall, our analysis confirms the involvement of the 5f orbitals in Th-N multiple bonds, and further demonstrates the value of solid-state NMR spectroscopy for interrogating actinide-ligand bonding.

actinides, covalent bonding, Nuclear Magnetic Reso↗
DOE OSTI2021

Helicity, spin, and infra-zilch of light: A Lorentz covariant formulation

Highlights: • Spin and orbital parts of angular momentum. • The helicity tensor contains the helicity, spin and infra-zilch currents. • The helicity tensor is a Noether current for a duality-symmetric action. • Symmetry generator for the helicity tensor. • The helicity tensor is conserved in Lorenz gauge. In this paper, a novel conserved Lorentz covariant tensor, termed the helicity tensor, is introduced in Maxwell theory. The conservation of the helicity tensor expresses the conservation laws contained in the helicity array, introduced by Cameron et al. (2012), including helicity, spin, as well as the spin-flux or infra-zilch. The Lorentz covariance of the helicity tensor is in contrast to previous formulations of the helicity hierarchy of conservation laws, which required the non-Lorentz covariant transverse gauge. The helicity tensor is shown to arise as a Noether current for a variational symmetry of a duality-symmetric Lagrangian for Maxwell theory. This symmetry transformation generalizes the duality symmetry and includes the symmetry underlying the conservation of the spin part of the angular momentum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2023

Robust In-Situ Strain Measurements to Monitor CO 2 Storage

The goal of this project was to develop and demonstrate robust instrumentation to monitor the in-situ strain tensor in order to improve the reliability and security of CO 2 storage in geologic formations. We met the original goals of the project and the major overarching accomplishment is the advancement of strain tensor monitoring from an intriguing concept to a commercially available technology with a solid foundation of novel instruments supported by theoretical analyses and validation experiments. The main accomplishments of the project are summarized below. We designed, built and evaluated nine new optical fiber strainmeters and tiltmeters using Michelson interferometers to measure deformation with ultra-high resolution at both shallow and deep point locations in the subsurface. These are the robust strainmeters that motivated the title of the project. We designed, built and evaluated a novel method of measuring distributed strain in optical fibers with nanostrain resolution, and cm-scale location, and sampling into the seismic band. The new method is called Coherence-length-gated Microwave Photonics Interfereometry (CMPI). CMPI technology has advantages over existing commercial DAS and DSS methods. We developed and demonstrated capabilities to deploy instruments in the field and used them to measure strain caused by ambient signals like barometric pressure and tides, as well as induced signals like surface loading and pore pressure changes from pumping tests. We deployed a working strainmeter at 1,700 ft depth, slightly above an active reservoir. This is to our knowledge the greatest depth a strainmeter has been deployed and the techniques we used can readily be extended to greater depths. Optical fiber borehole tensor strainmeter techology was advanced from a TRL 4 at the start, to a TRL of 7 at the conclusion of the project. The project included advances in simulations and theoretical analyses. We developed and demonstrated a computational workflow that uses machine learning to reduce the computational requirements and make it practical to use Bayesian inversion to solve large numerical poroelastic analyses needed to interpret strain tensor field data. We evaluated the strain tensor fields and time series that would be caused by leaks of CO 2 or other fluids from reservoirs. These simulations demonstrated that signals from leaks could be measured with instruments developed for the project, opening a potentially new method for ensuring storage security. We showed that strains in caprock can be used to estimate pressure in a reservoir. This avoids the need to drill monitoring wells into the reservoir, and it expands the capabilities of monitoring in the caprock. The project includes a derivation and application of a novel analytical solution to the strains in the vicinity of a pressurized poroelastic inclusion. This solution explains field data measured during injeciton tests at the North Avant Field, and it will simplify future interpretation of strain tensor data. The project included a broad range of experiments, and of the most significant is the characterization of the strain tensor at an array three strainmeters during six injection tests at the North Avant Field, Oklahoma. This demonstrated repeatability of the strain signal measured by the new instruments developed for the project, and it showed similarities between the strain signal at shallow depths and pressure in the underlying reservoir. We also demonstrated that useful strain data can be measured at reservoir depths. This confirms that strain tensor data can be measured throughout the caprock over a reservoir. The project demonstrated the feasibility of using the strain tensor and distributed strain measured in caprock during a variety of different well tests where the pumping rate was constant, sinusoidal and positive, or a periodic square wave with zero net rate. This further strengthens the validity of using strain data to characterize reservoirs and aquifers. We also demonstrated that strain caused be fluctuations of air pressure and water pressure in the vadose zone can be measured and interpreted, suggesting that high resolution distributed strain measurements hold promise for monitoring the vadose zone. The project partially supported nine graduate students in the Environmental Engineering, Hydrogeology, Electrical Engineering programs at Clemson University. The research was described in nine journal papers, 23 talks and conference abstracts. Additional journal papers are in preparation. A new company called Tensora was started to provide strainmeter technology for commercial applications.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Brown-York charges at null boundaries

The Brown-York stress tensor provides a means for defining quasilocal gravitational charges in subregions bounded by a timelike hypersurface. We consider the generalization of this stress tensor to null hypersurfaces. Such a stress tensor can be derived from the on-shell subregion action of general relativity associated with a Dirichlet variational principle, which fixes an induced Carroll structure on the null boundary. The formula for the mixed-index tensor T i j takes a remarkably simple form that is manifestly independent of the choice of auxiliary null vector at the null surface, and we compare this expression to previous proposals for null Brown-York stress tensors. The stress tensor we obtain satisfies a covariant conservation equation with respect to any connection induced from a rigging vector at the hypersurface, as a result of the null constraint equations. For transformations that act covariantly on the boundary structures, the Brown-York charges coincide with canonical charges constructed from a version of the Wald-Zoupas procedure. For anomalous transformations, the charges differ by an intrinsic functional of the boundary geometry, which we explicity verify for a set of symmetries associated with finite null hyper-surfaces. Applications of the null Brown-York stress tensor to symmetries of asymptotically flat spacetimes and celestial holography are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Model-Free State Estimation Using Low-Rank Canonical Polyadic Decomposition

As electric grids experience high penetration levels of renewable generation, fundamental changes are required to address real-time situational awareness. Here, we utilize unique traits of tensors to devise a model-free situational awareness and energy forecasting framework for distribution networks. This work formulates the state of the network at multiple time instants as a three-way tensor; hence, recovering full state information of the network is tantamount to estimating all the values of the tensor. Given measurements received from µphasor measurement units and/or smart meters, the recovery of unobserved quantities is carried out using the low-rank canonical polyadic decomposition of the state tensor—that is, the state estimation task is posed as a tensor imputation problem utilizing observed patterns in measured (sampled) quantities. Two structured sampling schemes are considered, namely, asynchronous slab and fiber sampling. For both schemes, we present sufficient conditions on the number of sampled slabs and fibers that guarantee identifiability of the factors of the state tensor. Numerical results demonstrate the ability of the proposed framework to achieve high estimation accuracy in multiple sampling scenarios.

42 ENGINEERING↗
DOE OSTI2020

Critique on “A Revisit of Navier-Stokes Equation”

In a recently published work (Sheng, 2020), Sheng identified inconsistencies in the concepts and principles leading to the Navier–Stokes equations and proposed a new asymmetric “friction tensor” to resolve these inconsistencies. We show that this friction tensor can be interpreted as the non-isotropic part of a vorticity-influenced viscous stress tensor following the analysis in Berdahl and Strang (1986). The asymmetric stress tensor however violates the principles of mechanical equilibrium in fluid flows that are free of external or body couples and is therefore non-physical in classical continuum hydrodynamics. We refute the arguments in favour of the asymmetric stress tensor in Sheng (2020) whilst also demonstrating that the friction tensor leads to the incorrect form of the Navier–Stokes equations in non-Cartesian co-ordinates and its ability to recover the correct governing equations in Cartesian co-ordinates is fortuitous.

asymmetric↗
DOE OSTI2025

SymProp: Scaling Sparse Symmetric Tucker Decomposition via Symmetry Propagation

Sparse symmetric tensors are an important class of tensors, and their decompositions serve as powerful tools for revealing low-rank structures. This paper introduces SymProp, a novel approach for scaling sparse symmetric Tucker decomposition by propagating symmetry through intermediate computations. SymProp optimizes two key computational kernels: Sparse Symmetric Tensor Times Same Matrix chain (S3 TTMc) for Higher-Order Orthogonal Iteration (HOOI) and Sparse Symmetric Tensor Times Same Matrix chain Times Core (S3 TTMcTC) for Higher-Order QR Iteration (HOQRI). Our method employs a metaprogramming-based index iteration approach to efficiently handle the upper triangular parts of intermediate dense symmetric tensors. SymProp achieves up to 50.9× speedup over SPLATT and up to 360.8× over Compressed Sparse Symmetric (CSS) format on the S3 TTMc operation. Moreover, our S3 TTMc and S3 TTMcTC implementations support tensor orders four levels higher than state-of-the-art methods. Our HOQRI demonstrates superior scalability and up to a 33.6× speedup over optimized HOOI. By enabling more scalable Tucker decompositions for higher orders, decomposition ranks, and dimension sizes, SymProp opens new possibilities for analyzing complex hypergraph structures in fields such as network science, data mining, and machine learning.

Li, Zecheng [North Carolina State University]↗
DOE OSTI2021

pnnl/HiParTI

A Hierarchical Parallel Tensor Infrastructure (HiParTI), is to support fast essential sparse tensor operations and tensor decompositions on multicore CPU and GPU architectures. It consists of sparse tensor decompositions, CANDECOMP/PARAFAC (CP) and Tucker decompositions, fundamental tensor operations, and tensor transformations

Li, Jiajia↗
DOE OSTI2021

PyMTensor 1.0.0

SAND2021-1654 O The Python Material Tensor code (PyMTensor) applies crystal symmetries to arbitrary material tensors. It uses symbolic algebra to exactly determine the zero tensor components and the relationships between the nonzero components for all crystallographic point groups. It is capable of computing components of higher-order material tensors needed when treating nonlinear effects. It has a simple and flexible interface for inputting tensor rank information and any known symmetries between tensor indices.Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Jensen, DanielS↗
DOE OSTI2026

BoBa

BoBa is a C++ software library for working with large matrices, tensors, and tensor decompositions. The library provides tools for dense matrix and tensor operations, tensor decompositions, and tensor decomposition methods that support modern CPU and GPU architectures. It includes portable abstractions for linear algebra, tensor algebra, and multidimensional computation. BoBa is intended for scientific computing applications that involve large multidimensional data sets or high dimensional mathematical models. Its capabilities support tasks such as data compression, linear algebra, efficient numerical computation, and the development of scalable algorithms for heterogeneous hardware. Tutorials, tests, and example applications are included to help users learn and apply the library.

Yao, Jin [Lawrence Livermore National Laboratory (↗
DOE OSTI2024

Collocation methods for nonlinear differential equations on low-rank manifolds

We introduce new methods for integrating nonlinear differential equations on low-rank manifolds. These methods rely on interpolatory projections onto the tangent space, enabling low-rank time integration of vector fields that can be evaluated entry-wise. A key advantage of our approach is that it does not require the vector field to exhibit low-rank structure, thereby overcoming significant limitations of traditional dynamical low-rank methods based on orthogonal projection. To construct the interpolatory projectors, we develop a sparse tensor sampling algorithm based on the discrete empirical interpolation method (DEIM) that parameterizes tensor train manifolds and their tangent spaces with cross interpolation. Using these projectors, we propose two time integration schemes on low-rank tensor train manifolds. The first scheme integrates the solution at selected interpolation indices and constructs the solution with cross interpolation. The second scheme generalizes the well-known orthogonal projector-splitting integrator to interpolatory projectors. We demonstrate the proposed methods with applications to several tensor differential equations arising from the discretization of partial differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A continuum framework for coupled solid deformation–fluid flow through anisotropic elastoplastic porous media

We present a continuum framework for coupled solid deformation–fluid flow in anisotropic elastoplastic porous media. A thermodynamic formulation of the coupled processes gives rise to an anisotropic Biot tensor that is a function of drained elastic tangent moduli tensor of the solid skeleton and the intrinsic bulk modulus of the solid constituent. Two effective stress measures emerge from the formulation, namely, σ', which is energy-conjugate to the elastic strain, and σ'', which is energy-conjugate to the plastic strain. For the special case of transverse isotropy that is commonly encountered in natural rocks, the Biot tensor can be expressed in terms of its normal and tangential components to the bedding plane, along with a microstructure tensor. Apart from its thermodynamic consistency, an advantage of this new formulation is that standard mixed finite element formulation can be employed to discretize the domain and solve initial boundary-value problems. We conduct plane strain simulations of coupled solid deformation–fluid flow in a transversely isotropic porous medium to demonstrate the impacts of material anisotropy, stress history, and the Biot tensor on the system response.

42 ENGINEERING↗
DOE OSTI2024

Cu(II) Stability and UV-Induced Electron Transfer in a Metal–Organic Hybrid: An EPR, DFT, and Crystallographic Characterization of Copper-Doped Zinc Creatininium Sulfate

Single-crystal X-ray diffraction and electron paramagnetic resonance (EPR) spectroscopic experiments, complemented by quantum chemical DFT calculations, were carried out on the copper-doped metal–organic hybrid and Tutton salt analogue zinc creatininium sulfate to determine its crystal structure, to characterize the electronic structure of the doped Cu(II) binding site, and to propose a pathway for an excited-state, proton-coupled electron transfer (PCET) process in UV-exposed crystals. The crystal structure is isomorphous to that of cadmium creatininium sulfate, which has the transition ion, not in direct coordination with the creatinine, but forming a hexahydrate complex, which is bridged to a creatininium through an intervening sulfate ion. The EPR g (2.446, 2.112, 2.082) and copper hyperfine (A Cu : -327, -59.6, 10.8 MHz) tensor parameters are consistent with doped copper replacing host zinc in the metal–hexahydrate complex. These parameters are similar to those observed for copper hexahydrate in doped Tutton salt systems at low temperature, where the unpaired electron occupies mainly the copper 3d x 2 –y 2 orbital. At room temperature in the Tutton systems, vibration couplings stemming from a dynamic Jahn–Teller effect cause tensor averaging which results in a reduction in their maximum g-tensor and hyperfine tensor values. However, like for the doped isomorphous Cd creatinine crystal, the Cu(II) EPR exhibits little, or no room temperature averaging compared to its low temperature pattern. Samples exposed to 254 nm UV light generate a carbon-centered free radical species, characterized by an isotropic g-tensor (g = 2.0029) and an alpha-proton hyperfine coupling (-24 -14 +4 G). These parameters identify it as a creatinine radical cation formed by the oxidative release of one of its C2 methylene hydrogens. DFT calculations confirm the unpaired electronic structures of both the Cu(II) site and free radical. The growth in radical concentration with an increase in the UV exposure time coincides with a decrease in the copper EPR signal, indicating a coupled light-induced oxidation reduction process. A comparison of the crystal structure with the EPR parameters and DFT results provides evidence for a UV-induced PCET.

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

Ab initio calculations of third-order elastic coefficients

Third-order elasticity (TOE) theory predicts strain-induced changes in second-order elastic coefficients (SOECs) and can model elastic wave propagation in stressed media. Although third-order elastic tensors have been determined based on first principles in previous studies, their current definition is based on an expansion of thermodynamic energy in terms of the Lagrangian strain near the natural, or zero pressure, reference state. This definition is inconvenient for predictions of SOECs under significant initial stresses. Therefore, when TOE theory is necessary to study the strain dependence of elasticity, the seismological community has resorted to an empirical version of the theory. This study reviews the thermodynamic definition of the third-order elastic tensor and proposes using an “effective” third-order elastic tensor. An explicit expression for the effective third-order elastic tensor is given and verified. Additionally, we extend the ab initio approach to calculate third-order elastic tensors under finite pressure and apply it to two cubic systems, namely, NaCl and MgO. As applications and validations, we evaluate (a) strain-induced changes in SOECs and (b) pressure derivatives of SOECs based on ab initio calculations. Good agreement between third-order elasticity-based predictions and numerically calculated values confirms the validity of our theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A PEPS Plugin for TNQVM

This work introduces an extension to the Tensor Network Quantum Virtual Machine (TNQVM) tool, enhancing the existing stack of ExaScale Tensor Network (ExaTN), ExaScale Accelerator (XACC), and TNQVM. It features a new plugin that enables efficient simulation of a Projected Entangled Pair State (PEPS), a 2D tensor network. To improve simulation efficiency for PEPS, we have implemented the snake boundary contraction algorithm. By integrating this capability into the existing stack, we enhance the overall functionality and versatility of the framework. We tested this new PEPS topology for a simple GHZ bell-pair generation quantum circuit and saw that its runtime is very close to that of the MPS topology. We estimate that the real potential of the PEPS topology becomes discernible when quantum circuits with multidimensional entanglement are simulated using tensor networks. In such cases, 1D tensor networks fail to represent or contract them efficiently.

Chundury, Srikar↗
DOE OSTI2022

Unraveling the Wrinkle in Time-Variable Sources with Lunes and Synthetic Seismic Data

In this report, we describe how to estimate the time-variable components of the seismic moment tensor and compare these estimates to the more conventional analysis that incorporates an assumption of the source time function (STF) across all components of the seismic moment tensor. The advantage of our method is that we are able to independently estimate the time-evolution of each component of the seismic moment tensor, which may help to resolve the complex source phenomena associated with buried explosions. By performing an eigen decomposition of the time-evolving seismic moment tensor components, we are able to plot the seismic mechanism as a trajectory on a lune diagram. This technique enables interpretation of the seismic mechanism as a function of time, as opposed to the more conventional analysis which assumes that the seismic mechanism is time invariant. Finally, we describe the differences between the seismic moment and the seismic moment rate STFs, how to implement each one in inversion schemes, and the relative strengths/weaknesses of each. Our key take-away is that we are able to distinguish nearly-overlapping sources with highly different mechanisms, such as an explosion immediately following an earthquake, by estimating moment rate from seismic data through a STF-invariant inversion for the full time-variable moment tensor.

58 GEOSCIENCES↗
DOE OSTI2024

Georgia Tech Accelerated, Compressed, and Regularized Compute of Kinetic-based PDEs (Final Report)

This report summarizes the collaborative effort between Lawrence Livermore National Laboratory and Georgia Tech to enhance the BoBa library for tensor train computation in PDE solvers, with a target on kinetic equations and their continuum limits. We aimed to reduce computational cost and memory usage by replacing traditional array-based computations with tensor trains. We examined the compressibility of time-evolving solutions to the Euler equations with discontinuities. We also explored using the first invsicid and linear regularization of the compressible flow equations via the information geometric regularization (IGR). We explored this in a tensor train formulation. To identify that inverse terms in the IGR equations pose problems for tensor train formulations and investigate efficient methods for batched inversion of tensor trains.

97 MATHEMATICS AND COMPUTING↗