Engineering PapersSearch

SEARCH · Engineering Papers

Results for “tensor times same vector”

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 19 records

Scalable Tensor Methods for Nonuniform Hypergraphs

While multilinear algebra appears natural for studying the multiway interactions modeled by hypergraphs, tensor methods for general hypergraphs have been stymied by theoretical and practical barriers. A recently proposed adjacency tensor is applicable to nonuniform hypergraphs, but is prohibitively costly to form and analyze in practice. We develop tensor times same vector (TTSV) algorithms for this tensor which improve complexity from $O(n^r)$ to a low-degree polynomial in $r$, where $n$ is the number of vertices and $r$ is the maximum hyperedge size. Our algorithms are implicit, avoiding formation of the order $r$ adjacency tensor. Here, we demonstrate the flexibility and utility of our approach in practice by developing tensor-based hypergraph centrality and clustering algorithms. We also show these tensor measures offer complementary information to analogous graph-reduction approaches on data, and are also able to detect higher-order structure that many existing matrix-based approaches provably cannot.

97 MATHEMATICS AND COMPUTING

Solar physics applications of computer graphics and image processing

Computer graphics devices coupled with computers and carefully developed software provide new opportunities to achieve insight into the geometry and time evolution of scalar, vector, and tensor fields and to extract more information quickly and cheaply from the same image data. Two or more different fields which overlay in space can be calculated from the data (and the physics), then displayed from any perspective, and compared visually. The maximum regions of one field can be compared with the gradients of another. Time changing fields can also be compared. Images can be added, subtracted, transformed, noise filtered, frequency filtered, contrast enhanced, color coded, enlarged, compressed, parameterized, and histogrammed, in whole or section by section. Today it is possible to process multiple digital images to reveal spatial and temporal correlations and cross correlations. Data from different observatories taken at different times can be processed, interpolated, and transformed to a common coordinate system.

Altschuler, M. D.

Search for charged-lepton flavour violation in top quark interactions with an up-type quark, a muon, and a $τ$ lepton in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A search for charged-lepton flavour violation (CLFV) in top quark (t) production and decay is presented. The search uses proton-proton collision data corresponding to 138 fb$^{-1}$ collected with the CMS experiment at $\sqrt{s}$ = 13 TeV. The signal consists of the production of a single top quark via a CLFV interaction or top quark pair production followed by a CLFV decay. The analysis selects events containing a hadronically decaying $τ$ lepton and a muon of opposite electric charge, as well as at least three jets, one of which is identified as originating from the fragmentation of a bottom quark. Machine learning classification techniques are used to distinguish signal from standard model background events. The results of this search are consistent with the standard model expectations. The upper limits at 95% confidence level on the branching fraction $\mathcal{B}$ for CLFV top quark decays to a muon, a $τ$ lepton, and an up or a charm quark are set at $\mathcal{B}$(t $\to$ $μτ$u) $\lt$ (0.04, 0.08, and 0.12) $\times$ 10$^{-6}$, and $\mathcal{B}$(t $\to$ $μτ$c) $\lt$ (0.81, 1.71, and 2.05) $\times$ 10$^{-6}$ for scalar, vector, and tensor-like operators, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Decoding the B → K ν ν excess at Belle II: Kinematics, operators, and masses

An excess in the branching fraction for B + → K + ν ν recently measured at Belle II may be a hint of new physics. We perform thorough likelihood analyses for different new physics scenarios such as B → K X with a new invisible particle X , or B → K χ χ through a scalar, vector, or tensor current with χ being a new invisible particle or a neutrino. We find that vector-current three-body decay with m X ≃ 0.6 GeV—which may be dark matter—is most favored, while two-body decay with m X ≃ 2 GeV is also competitive. The best-fit branching fractions for the scalar and tensor cases are a few times larger than for the two-body and vector cases. Past measurements provide further discrimination, although the best-fit parameters stay similar. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Determination of efficiencies, loss mechanisms, and performance degradation factors in chopper controlled dc vehical motors. Section 2: The time dependent finite element modeling of the electromagnetic field in electrical machines: Methods and applications

The time dependent solution of the magnetic field is introduced as a method for accounting for the variation, in time, of the machine parameters in predicting and analyzing the performance of the electrical machines. The method of time dependent finite element was used in combination with an also time dependent construction of a grid for the air gap region. The Maxwell stress tensor was used to calculate the airgap torque from the magnetic vector potential distribution. Incremental inductances were defined and calculated as functions of time, depending on eddy currents and saturation. The currents in all the machine circuits were calculated in the time domain based on these inductances, which were continuously updated. The method was applied to a chopper controlled DC series motor used for electric vehicle drive, and to a salient pole sychronous motor with damper bars. Simulation results were compared to experimentally obtained ones.

Hamilton, H. B.

Riemann Curvature tensor in nonholonomic coordinates and non-Riemannian space-times

The calculation of the Riemann Curvature from the deviation of a vector undergoing parallel transport around a closed loop takes a very simple form when expressed in generalized geometrical notation. The parallel transport of a vector and the type of closed loop used in the calculation are discussed. The method generalizes similar work of Morganstern (1977) to nonholonomic coordinates and non-Riemannian space-times.

Smalley, L. L.

NEML2: A High Performance Library for Constitutive Modeling

NEML2, the New Engineering Material model Library, version 2, is an offshoot of NEML, an earlier material modeling code developed at Argonne National Laboratory. NEML2 extends the key philosophy of its predecessor, i.e., material models are flexible, modular, and can be built from smaller blocks. It also provides modern features that do not exist in the framework of its predecessor such as material model vectorization, automatic differentiation, device-portable just-in-time compilation, operator fusion, lazy tensor evaluation, etc. Moreover, NEML2 can seamlessly integrate with the popular machine learning package PyTorch to take advantage of modern and fast-growing machine learning techniques. In this fiscal year, the development of core library features and capabilities are complete. The purpose of this report is not to serve as a verbatim copy of the software API reference (which is available online at https://reverendbedford.github.io/neml2/). Instead, this report documents the motivation, implementation, design choices, and usage of each core capability as well as their applications in solving practical engineering problems. This report is compiled based on the NEML2 major release 2.0.0.

36 MATERIALS SCIENCE

Numerical solutions of the complete Navier-Stokes equations

Using ideas from the kinetic theory, the Navier-Stokes equations are modified in such a way that they can be cast as a set of first order hyperbolic equations. This is achieved by incorporating time dependent terms into the definition of the stress tensor and the heat flux vectors. The boundary conditions are then determined from the theory of characteristics. Because the resulting equations reduce to the traditional Navier-Stokes equations when the steady state is reached, the present approach provides a straightforward scheme for the determination of inflow and outflow boundary conditions. The method is validated by comparing its predictions with known exact solutions of the steady Navier-Stokes equations.

Hassan, H. A.

Boundary conditions for the Navier-Stokes equations

Using ideas from kinetic theory, the Navier-Stokes equations are modified in such a way that they can be cast as a set of first order hyperbolic equations. This is achieved by incorporating time dependent terms into the definition of the stress tensor and the heat flux vectors. The boundary conditions are then determined from the theory of characteristics. Because the resulting equations reduce to the traditional Navier-Stokes equations when the steady state is reached, the present approach provides a straightforward scheme for the determination of inflow and outflow boundary conditions. The method is validated by comparing its predictions with exact and numerical solutions of the steady Navier-Stokes equations.

Rudisill, Edgar N., Jr.

Turbulence Modeling Effects on the Prediction of Equilibrium States of Buoyant Shear Flows

The effects of turbulence modeling on the prediction of equilibrium states of turbulent buoyant shear flows were investigated. The velocity field models used include a two-equation closure, a Reynolds-stress closure assuming two different pressure-strain models and three different dissipation rate tensor models. As for the thermal field closure models, two different pressure-scrambling models and nine different temperature variance dissipation rate, Epsilon(0) equations were considered. The emphasis of this paper is focused on the effects of the Epsilon(0)-equation, of the dissipation rate models, of the pressure-strain models and of the pressure-scrambling models on the prediction of the approach to equilibrium turbulence. Equilibrium turbulence is defined by the time rate (if change of the scaled Reynolds stress anisotropic tensor and heat flux vector becoming zero. These conditions lead to the equilibrium state parameters. Calculations show that the Epsilon(0)-equation has a significant effect on the prediction of the approach to equilibrium turbulence. For a particular Epsilon(0)-equation, all velocity closure models considered give an equilibrium state if anisotropic dissipation is accounted for in one form or another in the dissipation rate tensor or in the Epsilon(0)-equation. It is further found that the models considered for the pressure-strain tensor and the pressure-scrambling vector have little or no effect on the prediction of the approach to equilibrium turbulence.

Zhao, C. Y.

Far-from-equilibrium slow modes and momentum anisotropy in an expanding plasma

The momentum distribution of particle production in heavy-ion collisions encodes information about thermalization processes in the early-stage quark-gluon plasma. We use kinetic theory to study the far-from-equilibrium evolution of an expanding plasma with an anisotropic momentum-space distribution. We identify slow and fast degrees of freedom in the far-from-equilibrium plasma from the evolution of moments of this distribution. At late times, the slow modes correspond to hydrodynamic degrees of freedom and are naturally gapped from the fast modes by the inverse of the relaxation time, τ R − 1 . At early times, however, there are an infinite number of slow modes with a gap inversely proportional to time, τ − 1 . From the evolution of the slow modes we generalize the paradigm of the far-from-equilibrium attractor to vector and tensor components of the energy-momentum tensor, and even to higher moments of the distribution function that are not part of the hydrodynamic evolution. We predict that initial-state momentum anisotropy decays slowly in the far-from-equilibrium phase and may persist until the relaxation time. Published by the American Physical Society 2024

Brewer, Jasmine (ORCID:0000000230840663)

Polar motions excited by a convecting viscous mantle

The role played by time-dependent mantle convection on exciting long-term polar motions is examined by means of a viscous model. For sufficiently low effective viscosity, polar-wander speeds of 0(1 deg/Myr) would require contributions from the off-diagonal elements of the inertia tensor, with magnitudes around 100,000 times smaller than those found in viscoelastic models used for postglacial rebound. Contributions from the large-scale mantle flow to the angular momentum vector can be comparable to those due to changes in moment of inertia tensor for nonlinear rheology or for young planet. The relative roles of the two contributions depend on the nonlinear rheology of the mantle and its convective vigor.

Moser, Jiri

Development of the Tensoral Computer Language

The research scientist or engineer wishing to perform large scale simulations or to extract useful information from existing databases is required to have expertise in the details of the particular database, the numerical methods and the computer architecture to be used. This poses a significant practical barrier to the use of simulation data. The goal of this research was to develop a high-level computer language called Tensoral, designed to remove this barrier. The Tensoral language provides a framework in which efficient generic data manipulations can be easily coded and implemented. First of all, Tensoral is general. The fundamental objects in Tensoral represent tensor fields and the operators that act on them. The numerical implementation of these tensors and operators is completely and flexibly programmable. New mathematical constructs and operators can be easily added to the Tensoral system. Tensoral is compatible with existing languages. Tensoral tensor operations co-exist in a natural way with a host language, which may be any sufficiently powerful computer language such as Fortran, C, or Vectoral. Tensoral is very-high-level. Tensor operations in Tensoral typically act on entire databases (i.e., arrays) at one time and may, therefore, correspond to many lines of code in a conventional language. Tensoral is efficient. Tensoral is a compiled language. Database manipulations are simplified optimized and scheduled by the compiler eventually resulting in efficient machine code to implement them.

Ferziger, Joel

Polar Hydra Data Analysis

The science activities are: 1) Hydra is still operating successfully on orbit. 2) A large amount of analysis and discovery has occurred with the Hydra ground data processing this past year. 3) Full interdetector calibration has been implemented and documented. This intercalibration was necessitated by the incorrect installation of bias resistors in the pre-acceleration stage to the electron channeltrons. This had the effect of making the counting efficiency for electrons energy dependent as well as channeltron specific. The nature of the error had no impact on the ion detection efficiency since they have a different bias arrangement. This intercalibration is so effective, that the electron and ion moment densities are routinely produced with a level of agreement better than 20%. 4) The data processing routinely removes glint in the sensors and produces public energy time spectrograms on the web overnight. 6) Routine, but more intensive computer processing codes are operational that determine for electrons and ions, the density, the flow vector, the pressure tensor and the heat flux by numerical integration. These codes use the magnetic field to sustain the quality of their output. To gain access to this high quality magnetic field within our data stream we have monitored Russell's web page for zero levels and timing files (since his data acquisition is not telemetry synchronous) and have a local reconstruction of B for our use. We have also detected a routine anomaly in the magnetometer data stream that we have documented to Chris Russell and developed an editing algorithm to intercept these "hits" and remove them from the geophysical analysis.

Scudder, J. D.

Topology of three-dimensional, variable density flows

This paper is concerned with the interpretation of unsteady, variable-density flow fields. The topology of the flow is determined by finding critical points and identifying the character of local solution trajectories. The time evolution of the flow is studied by following the paths of the critical points in the three-dimensional space of invariants of the local deformations tensor. The methodology can be applied to any smooth vector field and its associated gradient tensor including the vorticity and pressure gradient fields. This approach provides a framework for describing the geometry of complex flow patterns. Concisely summarizing that geometry in the space of invariants of the local gradient tensor may be a useful way of gaining insight into time-dependent processes described by large computational data bases. Applications to the descriptions of a flickering diffusion flame and a compressible wake are discussed.

Cantwell, Brian

Geometric Interpretation of the Cluster Location Problem Part II: Application to the Pahala, Hawaii, Earthquake Sequence

In the companion “Theory” article, we presented a new framing of the seismic location problem in terms of differential geometry (Harris et al., 2025). From that viewpoint, we developed a “project and correct” approach for estimating the relative locations of earthquakes. Here, in this study, we use project and correct to estimate high-precision relative locations of events from an earthquake sequence beneath the town of Pahala, Hawaii, using high-precision correlation-derived picks. The sequence was active from 2020 through 2022 and produced many highly correlated signals at Hawaii Volcano Observatory (HVO) stations on the island of Hawaii. The data we inverted consisted of 2882 events with observations at 5 HVO stations. For comparison with the travel-time image, we also produced conventional hypocenter solutions using both the Bayesloc program (Myers et al., 2007, 2009) and a purpose-built double-difference code. There were obvious structural elements in the resulting image, the resolution of which we used to test the performance of the project and the correct algorithm. For the projection step, we first produced a 3D local basis using an singular value decomposition (SVD) of the 2882 groups of times. Projection of the travel-time vectors into this basis resulted in an image with structures similar to those produced by our conventional locators, but with distortion as predicted by theory. Removing the distortion requires an inverse operator generated from the metric tensor at the geometric centroid of the events. We compared two approaches to obtaining such an inverse operator. The first uses an estimate of the geographic centroid of the event cloud from the centroid of the travel-time data. The second approach uses the centroid of the conventionally produced locations. The first approach produces a corrected image very similar to the conventional results, but with a rotation. The corrected image produced using the conventionally derived centroid is a near-exact match to the conventional locations.

Dodge, Douglas A. [Lawrence Livermore National Lab

Optical Computation Of Matrices From Vectors

Proposed optical apparatus generates rectangular pattern of light and dark areas, brightnesses of which represent elements of matrix product of two vectors. Photorefractive effect gives rise to four-wave mixing, which generates output beam modulated spatially by matrix product xy(SupT). Optical multipliers used as real-time analog tensor generators, components of neural networks, and generators of patterns to steer electromagnetic beams by diffraction or to effect temporary interconnections in very-large-scale integrated circuits.

Liu, Hua-Kuang

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