Engineering Papers⌕ Search

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Explicit simulation of the Brownian rotation of arbitrary shaped aerosol particles using quaternions

The shape of an aerosol particle strongly influences its mass and momentum transfer cross-sections, charging properties, and other physical properties. Here, we present an explicit time-stepping procedure to simulate the rotational Brownian motion of arbitrary shaped aerosol particles by solving Euler’s equation of rotation. A Langevin formulation of the rotation equations is used, wherein Brownian motion due to thermal collisions between a particle and background gas molecules is represented using a stochastic fluctuating torque and fluid resistance is included as a drag torque. To avoid singularities associated with describing the orientation of a shape with Euler angles, we employ a quaternion formulation that leads to first-order stochastic differential equations to describe the evolution of the angular position and angular velocity of a rigid body. We perform all the rotational dynamics calculations in the body-fixed frame of reference attached to the rotating shape whose basis vectors are the normalized eigenvectors of the inertia tensor of the particle. Numerical solutions to rotation under torque-free conditions, damped rotation without Brownian motion, and stochastic rotation for arbitrary shapes are presented and discussed. The presented method enables time-resolved simulation of Brownian rotation for direct comparison with experimentally measured trajectories or statistical measures. The second order accuracy of the used time-stepping procedure places a severe restriction on the timestep that can be used for obtaining accurate results. Animations of presented simulations are included for visualizing rotational motion at various gas pressures. To aid implementation, MATLAB ® codes are also provided. Extension to include translation Brownian motion is straightforward.

Roy, Mrittika↗
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.↗
NASA NTRS1993

Tensoral for post-processing users and simulation authors

The CTR post-processing effort aims to make turbulence simulations and data more readily and usefully available to the research and industrial communities. The Tensoral language, which provides the foundation for this effort, is introduced here in the form of a user's guide. The Tensoral user's guide is presented in two main sections. Section one acts as a general introduction and guides database users who wish to post-process simulation databases. Section two gives a brief description of how database authors and other advanced users can make simulation codes and/or the databases they generate available to the user community via Tensoral database back ends. The two-part structure of this document conforms to the two-level design structure of the Tensoral language. Tensoral has been designed to be a general computer language for performing tensor calculus and statistics on numerical data. Tensoral's generality allows it to be used for stand-alone native coding of high-level post-processing tasks (as described in section one of this guide). At the same time, Tensoral's specialization to a minute task (namely, to numerical tensor calculus and statistics) allows it to be easily embedded into applications written partly in Tensoral and partly in other computer languages (here, C and Vectoral). Embedded Tensoral, aimed at advanced users for more general coding (e.g. of efficient simulations, for interfacing with pre-existing software, for visualization, etc.), is described in section two of this guide.

Dresselhaus, Eliot↗
DOE OSTI2025

Impossibility of obtaining time-independent, three-dimensional, spherically symmetric densities of confined systems of relativistically moving constituents

The quantum-mechanical definition of probability, the uncertainty principle, and Poincaré invariance provide strong basic restrictions on the ability to define spatial densities associated with form factors describing the properties of confined systems of relativistically moving constituents. Despite this, many papers ignore one or more of these restrictions. Here I show how to obtain time-independent , two-dimensional densities that are consistent with the stated restrictions. This is done using the light-front, infinite momentum frame formalism. Two-dimensional density interpretations of the axial-vector form factor and all three gravitational form factors is obtained. The resulting mass radius is smaller than the charge radius. Additionally, an expression of a two-dimensional mass density related to the trace of the energy momentum tensor is obtained. I also show that all known methods for finding three-dimensional densities—using the Breit frame, Abel transformations, Wigner distributions, and spherically symmetric wave packets with vanishing spatial extent—violate the basic restrictions in different ways. Furthermore, the use of the latter leads to densities that vanish almost everywhere in space as time increases from an initial value.

form factors↗
NASA NTRS1986

Interactive color graphical postprocessing as a unifying influence in numerical analysis research

The general capabilities and the hardware environment of an interactive color graphical postprocessor for three-dimensional numerical analysis are described as well as the synergistic influence of this aspect of interactive graphics on the research environment. Areas of practical concern include finite element analysis, fracture propagation, surface and framed structures. The interactive color graphical postprocessor, POSTPRO3D, is capable of displaying the scalar, vector, and tensor results of finite difference, finite element, and boundary element analysis. In the case of fatigue cracking in steel crane-rail girders, it was shown that POSTPRO3D offers a high level of flexibility, interactivity and features which minimize response time.

Abel, J. F.↗
DOE OSTI2025

Thermal Radiation Transport with Tensor Trains

We present a novel tensor network algorithm to solve the time-dependent, gray thermal radiation transport equation. The method invokes a tensor train (TT) decomposition for the specific intensity. The efficiency of this approach is dictated by the rank of the decomposition. When the solution is “low rank,” the memory footprint of the specific intensity solution vector may be significantly compressed. The algorithm, following a step-then-truncate approach of a traditional discrete ordinates method, operates directly on the compressed state vector, thereby enabling large speedups for low-rank solutions. To achieve these speedups, we rely on a recently developed rounding approach based on the Gram-SVD. We detail how familiar S N algorithms for (gray) thermal transport can be mapped to this TT framework and present several numerical examples testing both the optically thick and thin regimes. The TT framework finds low-rank structure and supplies up to ≃60× speedups and ≃1000× compressions for problems demanding large angle counts, thereby enabling previously intractable SN calculations and supplying a promising avenue to mitigate ray effects.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1979

Finite-strain large-deflection elastic-viscoplastic finite-element transient response analysis of structures

A method of analysis for thin structures that incorporates finite strain, elastic-plastic, strain hardening, time dependent material behavior implemented with respect to a fixed configuration and is consistently valid for finite strains and finite rotations is developed. The theory is formulated systematically in a body fixed system of convected coordinates with materially embedded vectors that deform in common with continuum. Tensors are considered as linear vector functions and use is made of the dyadic representation. The kinematics of a deformable continuum is treated in detail, carefully defining precisely all quantities necessary for the analysis. The finite strain theory developed gives much better predictions and agreement with experiment than does the traditional small strain theory, and at practically no additional cost. This represents a very significant advance in the capability for the reliable prediction of nonlinear transient structural responses, including the reliable prediction of strains large enough to produce ductile metal rupture.

Rodal, J. J. A.↗
DOE OSTI2026

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗
DOE OSTI2026

Adiabatic fast passage spin manipulation measurements in solid polarized targets

Adiabatic fast passage (AFP) is a rapid method for reversing nuclear polarization and manipulating spin populations in polarized solid targets, avoiding the long repolarization times associated with dynamic nuclear polarization (DNP). We report AFP measurements in a 5 T, 1 K polarized-target system for irradiated 15 NH 3 , irradiated 14 ND 3 , and butanol-based materials prepared either with TEMPO doping or by irradiation. We also present a joint manipulated-lineshape analysis for spin-1 targets and demonstrate that vector and tensor polarizations can be extracted from AFP-manipulated deuteron NMR spectra even when the populations are not described by a single Boltzmann spin temperature. In conclusion, we report a reproducible polarization- and direction-dependent AFP response in a large irradiated 15 NH 3 sample. These ammonia results are presented as empirical observations under the specific sample–coil conditions of the experiment, with possible circuit-mediated mechanisms such as radiation damping or superradiant behavior discussed but not assigned as a definitive cause.

Adiabatic fast passage↗
NASA NTRS2002

Definition of Contravariant Velocity Components

This is an old issue in computational fluid dynamics (CFD). What is the so-called contravariant velocity or contravariant velocity component? In the article, we review the basics of tensor analysis and give the contravariant velocity component a rigorous explanation. For a given coordinate system, there exist two uniquely determined sets of base vector systems - one is the covariant and another is the contravariant base vector system. The two base vector systems are reciprocal. The so-called contravariant velocity component is really the contravariant component of a velocity vector for a time-independent coordinate system, or the contravariant component of a relative velocity between fluid and coordinates, for a time-dependent coordinate system. The contravariant velocity components are not physical quantities of the velocity vector. Their magnitudes, dimensions, and associated directions are controlled by their corresponding covariant base vectors. Several 2-D (two-dimensional) linear examples and 2-D mass-conservation equation are used to illustrate the details of expressing a vector with respect to the covariant and contravariant base vector systems, respectively.

Hung, Ching-Mao↗
DOE OSTI2024

Scalable learning of potentials to predict time-dependent Hartree–Fock dynamics

We propose a framework to learn the time-dependent Hartree–Fock (TDHF) inter-electronic potential of a molecule from its electron density dynamics. Although the entire TDHF Hamiltonian, including the inter-electronic potential, can be computed from first principles, we use this problem as a testbed to develop strategies that can be applied to learn a priori unknown terms that arise in other methods/approaches to quantum dynamics, e.g., emerging problems such as learning exchange–correlation potentials for time-dependent density functional theory. We develop, train, and test three models of the TDHF inter-electronic potential, each parameterized by a four-index tensor of size up to 60 × 60 × 60 × 60. Two of the models preserve Hermitian symmetry, while one model preserves an eight-fold permutation symmetry that implies Hermitian symmetry. Across seven different molecular systems, we find that accounting for the deeper eight-fold symmetry leads to the best-performing model across three metrics: training efficiency, test set predictive power, and direct comparison of true and learned inter-electronic potentials. All three models, when trained on ensembles of field-free trajectories, generate accurate electron dynamics predictions even in a field-on regime that lies outside the training set. To enable our models to scale to large molecular systems, we derive expressions for Jacobian-vector products that enable iterative, matrix-free training.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

An Introduction to Tensors for Students of Physics and Engineering

Tensor analysis is the type of subject that can make even the best of students shudder. My own post-graduate instructor in the subject took away much of the fear by speaking of an implicit rhythm in the peculiar notation traditionally used, and helped us to see how this rhythm plays its way throughout the various formalisms. Prior to taking that class, I had spent many years "playing" on my own with tensors. I found the going to be tremendously difficult but was able, over time, to back out some physical and geometrical considerations that helped to make the subject a little more transparent. Today, it is sometimes hard not to think in terms of tensors and their associated concepts. This article, prompted and greatly enhanced by Marlos Jacob, whom I've met only by e-mail, is an attempt to record those early notions concerning tensors. It is intended to serve as a bridge from the point where most undergraduate students "leave off" in their studies of mathematics to the place where most texts on tensor analysis begin. A basic knowledge of vectors, matrices, and physics is assumed. A semi-intuitive approach to those notions underlying tensor analysis is given via scalars, vectors, dyads, triads, and higher vector products. The reader must be prepared to do some mathematics and to think. For those students who wish to go beyond this humble start, I can only recommend my professor's wisdom: find the rhythm in the mathematics and you will fare pretty well.

Kolecki, Joseph C.↗
NASA NTRS2010

Introduction to Vector Field Visualization

Vector field visualization techniques are essential to help us understand the complex dynamics of flow fields. These can be found in a wide range of applications such as study of flows around an aircraft, the blood flow in our heart chambers, ocean circulation models, and severe weather predictions. The vector fields from these various applications can be visually depicted using a number of techniques such as particle traces and advecting textures. In this tutorial, we present several fundamental algorithms in flow visualization including particle integration, particle tracking in time-dependent flows, and seeding strategies. For flows near surfaces, a wide variety of synthetic texture-based algorithms have been developed to depict near-body flow features. The most common approach is based on the Line Integral Convolution (LIC) algorithm. There also exist extensions of LIC to support more flexible texture generations for 3D flow data. This tutorial reviews these algorithms. Tensor fields are found in several real-world applications and also require the aid of visualization to help users understand their data sets. Examples where one can find tensor fields include mechanics to see how material respond to external forces, civil engineering and geomechanics of roads and bridges, and the study of neural pathway via diffusion tensor imaging. This tutorial will provide an overview of the different tensor field visualization techniques, discuss basic tensor decompositions, and go into detail on glyph based methods, deformation based methods, and streamline based methods. Practical examples will be used when presenting the methods; and applications from some case studies will be used as part of the motivation.

Kao, David↗
NASA NTRS1995

Global Dynamic Modeling of Space-Geodetic Data

The proposal had outlined a year for program conversion, a year for testing and debugging, and two years for numerical experiments. We kept to that schedule. In first (partial) year, author designed a finite element for isostatic thin-shell deformation on a sphere, derived all of its algebraic and stiffness properties, and embedded it in a new finite element code which derives its basic solution strategy (and some critical subroutines) from earlier flat-Earth codes. Also designed and programmed a new fault element to represent faults along plate boundaries. Wrote a preliminary version of a spherical graphics program for the display of output. Tested this new code for accuracy on individual model plates. Made estimates of the computer-time/cost efficiency of the code for whole-earth grids, which were reasonable. Finally, converted an interactive graphical grid-designer program from Cartesian to spherical geometry to permit the beginning of serious modeling. For reasons of cost efficiency, models are isostatic, and do not consider the local effects of unsupported loads or bending stresses. The requirements are: (1) ability to represent rigid rotation on a sphere; (2) ability to represent a spatially uniform strain-rate tensor in the limit of small elements; and (3) continuity of velocity across all element boundaries. Author designed a 3-node triangle shell element which has two different sets of basis functions to represent (vector) velocity and all other (scalar) variables. Such elements can be shown to converge to the formulas for plane triangles in the limit of small size, but can also applied to cover any area smaller than a hemisphere. The difficult volume integrals involved in computing the stiffness of such elements are performed numerically using 7 Gauss integration points on the surface of the sphere, beneath each of which a vertical integral is performed using about 100 points.

Bird, Peter↗
DOE OSTI2025

Beyond-classical computation in quantum simulation

Quantum computers hold the promise of solving certain problems that lie beyond the reach of conventional computers. However, establishing this capability, especially for impactful and meaningful problems, remains a central challenge. Here, we show that superconducting quantum annealing processors can rapidly generate samples in close agreement with solutions of the Schrödinger equation. We demonstrate area-law scaling of entanglement in the model quench dynamics of two-, three-, and infinite-dimensional spin glasses, supporting the observed stretched-exponential scaling of effort for matrix-product-state approaches. We show that several leading approximate methods based on tensor networks and neural networks cannot achieve the same accuracy as the quantum annealer within a reasonable time frame. Thus, quantum annealers can answer questions of practical importance that may remain out of reach for classical computation.

King, Andrew D. [D-Wave Quantum Inc., Burnaby, BC ↗
NASA NTRS2014

The Primordial Inflation Polarization Explorer (PIPER)

The Primordial Inflation Polarization ExploreR (Piper) is a balloon-borne cosmic microwave background (CMB) polarimeter designed to search for evidence of inflation by measuring the large-angular scale CMB polarization signal. Bicep2 recently reported a detection of B-mode power corresponding to the tensor-to-scalar ratio r = 0.2 on approximately 2 degree scales. If the Bicep2 signal is caused by inflationary gravitational waves (IGWs), then there should be a corresponding increase in B-mode power on angular scales larger than 18 degrees. Piper is currently the only suborbital instrument capable of fully testing and extending the Bicep2 results by measuring the B-mode power spectrum on angular scales theta θ = approximately 0.6 deg to 90 deg, covering both the reionization bump and recombination peak, with sensitivity to measure the tensor-to-scalar ratio down to r = 0.007, and four frequency bands to distinguish foregrounds. Piper will accomplish this by mapping 85% of the sky in four frequency bands (200, 270, 350, 600 GHz) over a series of 8 conventional balloon flights from the northern and southern hemispheres. The instrument has background-limited sensitivity provided by fully cryogenic (1.5 K) optics focusing the sky signal onto four 32×40-pixel arrays of time-domain multiplexed Transition-Edge Sensor (TES) bolometers held at 140 milli-Kelvin. Polarization sensitivity and systematic control are provided by front-end Variabledelay Polarization Modulators (VPMs), which rapidly modulate only the polarized sky signal at 3 Hz and allow Piper to instantaneously measure the full Stokes vector (I,Q,U,0V) for each pointing. We describe the Piper instrument and progress towards its first flight.

polarimeter↗
DOE OSTI2025

Numerically exact configuration interaction at quadrillion-determinant scale

The combinatorial growth of configuration interaction (CI) has long limited this formally exact quantum chemistry method to only the smallest molecules. Here, we report a numerically exact CI calculation exceeding one quadrillion (10 15 ) determinants, made possible by a lossless categorical compression strategy within the small-tensor-product distributed active space (STP-DAS) framework. This approach overcomes the traditional memory bottlenecks of CI by a numerically exact compression of the wavefunction representation and reformulating the most computationally demanding matrix–vector operations. Using this method, we performed a fully relativistic CI calculation of the ground state of HBrTe with over 10 15 complex-valued determinants in just 34.5 h on 1000 computing nodes—the largest CI calculation ever reported. We further achieved fast computation for systems with hundreds of billions of determinants on only a few compute nodes. Extensive benchmarks confirm that the method retains full numerical exactness while cutting memory and computational cost by orders of magnitude. Compared to previous state-of-the-art CI calculations, this work achieves a 1000 times increase in CI space, a 10 6 -fold increase in floating-point operations performed, and a 10 6 -fold improvement in computational speed.

Computational chemistry↗
DOE OSTI2025

Accelerating GNNs on GPU Sparse Tensor Cores through N:M Sparsity-Oriented Graph Reordering

Recent advancements in GPU hardware support have introduced the capability to leverage N:M sparse patterns for substantial performance gains. Graphs in Graph Neural Networks (GNNs) are typically sparse, but the sparsity is often irregular, not conforming to such sparse patterns. In this paper, we propose a novel graph reordering algorithm, the first of its kind, to reshape irregular graph data into the N:M structured sparse pattern at the tile level, allowing linear-algebra-based graph operations in GNNs to benefit from the N:M sparse hardware. The optimization is lossless, maintaining the accuracy of GNN. It can remove 98-100\% violations of the N:M sparse patterns at the vector level, and increase the proportion of conforming graphs in SuiteSparse collection from 5-9\% to 88.7-93.5\%. On A100 GPUs, the optimization accelerates Sparse Matrix Matrix (SpMM) by up to 43X (2.3X -- 7.5X on average) and speeds up the key graph operations in GNNs on real graphs by as much as 8.6X (3.5X on average).

artificial intelligence, graph neural networks↗