SEARCH · Engineering Papers

Results for “Tensor factorizations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Transmural Remodeling of Cardiac Microstructure in Aged Spontaneously Hypertensive Rats by Diffusion Tensor MRI

The long-standing high blood pressure (also known as hypertension) overworks the heart. Microstructural remodeling is a key factor of hypertensive heart disease progression. Diffusion tensor magnetic resonance imaging (DT-MRI) is a powerful tool for the rapid noninvasive nondestructive delineation of the cardiomyocyte organization. The spontaneously hypertensive rat (SHR) is a well-established model of genetic hypertension. The goal of this study was to employ high-resolution DT-MRI and the SHR animal model to assess the transmural layer-specific remodeling of myocardial microstructure associated with hypertension. Ex vivo experiments were performed on excised formalin-fixed hearts of aged SHRs (n = 4) and age-matched controls (n = 4). The DT-MRI-derived fractional anisotropy (FA), longitudinal diffusivity (λ L ), transversal diffusivity (λ T ), and mean diffusivity (MD) served as the readout parameters investigated at three transmural zones (i.e., endocardium, mesocardium, and epicardium). The helix angles (HAs) of the aggregated cardiomyocytes and the orientation of laminar sheetlets were also studied. Compared with controls, the SHRs exhibited decreased epicardial FA, while FA changes in the other two transmural regions were insignificant. No substantial differences were observed in the diffusivity parameters and the transmural course of HAs between the two groups. A consistent distribution pattern of laminar sheetlet orientation was not identified for either group. Our findings are in line with the known cellular microstructure from early painstaking histological studies. Biophysical explanations of the study outcomes are provided. In conclusion, our experimental findings indicate that the epicardial microstructure is more vulnerable to high blood pressure leading to more pronounced changes in this region during remodeling. DT-MRI is well-suited for elucidating these alterations. The revealed transmural nonuniformity of myocardial reorganization may shed light on the mechanisms of the microstructure-function relationship in hypertension progression. Our results provide insights into the management of patients with systemic arterial hypertension, thus prevent the progression toward heart failure. The findings of this study should be acknowledged by electromechanical models of the heart that simulate the specific cardiac pathology.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2017

Phenomenology of $ {\Lambda}_b\to {\Lambda}_c\tau {\overline{\nu}}_{\tau } $ using lattice QCD calculations

In a recent paper we studied the effect of new-physics operators with different Lorentz structures on the semileptonic Λ b → Λ c τ ν ¯ τ decay. This decay is of interest in light of the R(D (*)) puzzle in the semileptonic B ¯ → D ( ∗ ) τ ν ¯ τ decays. In this work we add tensor operators to extend our previous results and consider both model-independent new physics (NP) and specific classes of models proposed to address the R(D (*)) puzzle. We show that a measurement of R ( Λ c ) = B [ Λ b → Λ c τ ν ¯ τ ] / B [ Λ b → Λ c ℓ ν ¯ ℓ ] can strongly constrain the NP parameters of models discussed for the R(D (*)) puzzle. We use form factors from lattice QCD to calculate all Λ b → Λ c τ ν ¯ τ observables. The Λ b → Λ c tensor form factors had not previously been determined in lattice QCD, and we present new lattice results for these form factors here.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2014

Communication: Acceleration of coupled cluster singles and doubles via orbital-weighted least-squares tensor hypercontraction

We apply orbital-weighted least-squares tensor hypercontraction decomposition of the electron repulsion integrals to accelerate the coupled cluster singles and doubles (CCSD) method. Using accurate and flexible low-rank factorizations of the electron repulsion integral tensor, we are able to reduce the scaling of the most vexing particle-particle ladder term in CCSD from O(N{sup 6}) to O(N{sup 5}), with remarkably low error. Combined with a T{sub 1}-transformed Hamiltonian, this leads to substantial practical accelerations against an optimized density-fitted CCSD implementation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2012

EQUATION OF STATE FOR MASSIVE NEUTRON STARS

Using the relativistic Hartree-Fock approximation, we investigate the properties of neutron-star matter in detail. In the present calculation, we consider not only the tensor coupling of vector mesons to octet baryons and the form factors at interaction vertices but also the change in the internal (quark) structure of baryons in dense matter. The relativistic Hartree-Fock calculations are performed in two ways: one with coupling constants determined by SU(6) (quark model) symmetry and the other with coupling constants based on SU(3) (flavor) symmetry. For the latter case, we use the latest Nijmegen (ESC08) model. Then, it is very remarkable that the particle composition of the core matter in SU(3) symmetry is completely different from that in SU(6) symmetry. In SU(6) symmetry, all octet baryons appear in the density region below {approx}1.2 fm{sup -3}, while in the ESC08 model only the {Xi}{sup -} hyperon is produced. Furthermore, the medium modification of the internal baryon structure hardens the equation of state for the core matter. Taking all these effects into account, we can obtain the maximum neutron-star mass which is consistent with the recently observed mass, 1.97 {+-} 0.04 M{sub Sun} (PSR J1614-2230). We therefore conclude that the extension from SU(6) symmetry to SU(3) symmetry in meson-baryon couplings and the internal baryon structure variation in matter considerably enhance the mass of neutron stars. Furthermore, the effects of the form factor at the vertex and the Fock contribution, including the tensor coupling due to vector mesons, are indispensable for describing the core matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Exploration of Electrochemical Reactions at Organic–Inorganic Halide Perovskite Interfaces via Machine Learning in In Situ Time-of-Flight Secondary Ion Mass Spectrometry

The instability of hybrid organic–inorganic perovskite (HOIP) devices is one of the significant challenges preventing commercialization. Exploring these phenomena is severely limited by the complexity of the intrinsic electrochemistry of HOIPs, the presence of multiple volatile and mobile ionic species, and the possible role of environmentally induced reactions at surfaces and triple-phase junctions. Here, in situ studies of the electrochemistry of methylammonium lead bromide perovskite with the Au electrode interface are reported via light- and voltage-dependent time-of-flight secondary ion mass spectrometry (ToF-SIMS) imaging of lateral perovskite heterostructures. While ToF-SIMS allows for the visualization of the chemical composition along the surface and its evolution with light and electrical bias, the interpretation of the multidimensional data obtained is often limited due to strong correlations between chemical signatures and the need to track multiple peaks at once. Here, a machine learning workflow combining the Hough transform and non-negative matrix factorization and no-negative tensor decomposition is developed to avoid this limitation and extract salient features of associated chemical changes and to separate the light- and voltage-dependent dynamics. Finally, combining these in situ characterizations and the machine learning workflow provides comprehensive information on the chemical nature of moving species, ion accumulation, and interfacial electrochemical reactions in HOIP devices.

36 MATERIALS SCIENCE↗

NASA NTRS2003

Macroscopic Models of Radiative Transfer as Applied to Computation of the Radiation Field in the Solar Atmosphere

In this paper a new three dimensional half-moment model for radiative transfer is presented for a gray medium. It describes the evolution of the zeroth and first directional half moments of the radiative intensity. The closure is provided, similarly to Dubroca and Klar, by the maximum entropy concept. This work generalizes that model to three dimensions. The model presented here (the derivation being done in Ripoll and Wray, called the M(sup 1/2)(sub 1) model, is a hyperbolic system consisting of a total of eight equations in three dimensions, four equations for each direction. Each half model has the classical form of a macroscopic moment model in which the pressure tensor is constructed from the well-known Eddington tensor with a particular Eddington factor. Moreover, different source and border terms occur. The latter introduce couplings between the macroscopic and microscopic quantities and between the + and - streams, through the intensity in the plane perpendicular to the flux. The main theoretical application of the half moment model, treated in this paper, is its reduction to a full moment model, called M(sup +)(sub 1), for the particular but important case of a hot, opaque source radiating in a cold transparent (or semi-transparent) medium for very specific applications, such as stellar interiors or atmospheres, or combustion problems. The structure of the paper is as follows. In section 2, the model M(sup 1/2)(sub 1) is presented. In section 3, for the particular case of a hot, opaque source radiating into a cold medium, the half moment model is reduced to the M(sup +)(sub 1) model. In section 4, we first solve a simple and academic problem to validate the models, followed by a simplified solar atmosphere.

Ripoll, J.-F.↗

DOE OSTI2023

The gravitational form factors of the electron in quantum electrodynamics

Not Available

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

spdlayers

Symmetric Positive Definite (SPD) enforcement layers for PyTorch. Regardless of the input, the output of these layers will always be a SPD tensor! The `Cholesky` layer uses a cholesky factorization to enforce SPD, and the `Eigen` layer uses an eigendecomposition to enforce SPD. Both layers take in some tensor of shape `[batch_size, input_shape]` and output a SPD tensor of shape`[batch_size, output_shape, output_shape]`. The relationship between input and output is defined by the following. ```python input_shape = sum([i for i in range(output_shape + 1)]) ``` The layers have no learnable parameters, and merely serve to transform a vector space to a SPD matrix space.

Jekel, CharlesF.↗

DOE OSTI2018

Parallel Nonnegative CP Decomposition of Dense Tensors

The CP tensor decomposition is a low-rank approximation of a tensor. We present a distributed-memory parallel algorithm and implementation of an alternating optimization method for computing a CP decomposition of dense tensors that can enforce nonnegativity of the computed low-rank factors. The principal task is to parallelize the Matricized-Tensor Times Khatri-Rao Product (MTTKRP) bottleneck subcomputation. The algorithm is computation efficient, using dimension trees to avoid redundant computation across MTTKRPs within the alternating method. Our approach is also communication efficient, using a data distribution and parallel algorithm across a multidimensional processor grid that can be tuned to minimize communication. We benchmark our software on synthetic as well as hyperspectral image and neuroscience dynamic functional connectivity data, demonstrating that our algorithm scales well to 100s of nodes (up to 4096 cores) and is faster and more general than the currently available parallel software.

Ballard, Grey↗

DOE OSTI2021

Multi-Area Model-Free State Estimation via Distributed Tensor Decomposition

This paper proposes a model-free method for distribution system state estimation based on tensor completion using canonical polyadic decomposition. In particular, we consider a setting where the network is divided into multiple areas. The measured physical quantities at buses located in the same area are processed by an area controller. A third-order tensor is constructed to collect these measured quantities. The measurements are analyzed locally to recover the full state information of the network. A closed-form iterative algorithm based on the alternating direction method of multipliers is developed to obtain the low-rank factors of the whole network state tensor where information exchange happens only between neighboring areas. To demonstrate the efficacy of the developed algorithm, numerical simulations are carried out using an IEEE test system.

alternating direction method of multipliers↗

NASA NTRS2009

Stress Concentration and Fracture at Inter-variant Boundaries in an Al-Li Alloy

Delamination fracture has limited the use of lightweight Al-Li alloys. Studies of secondary, delamination cracks in alloy 2090, L-T fracture toughness samples showed grain boundary failure between variants of the brass texture component. Although the adjacent texture variants, designated B(sub s1) and B(sub s2), behave similarly during rolling, their plastic responses to mechanical tests can be quite different. EBSD data from through-thickness scans were used to generate Taylor factor maps. When a combined boundary normal and shear tensor was used in the calculation, the delaminating grains showed the greatest Taylor Factor differences of any grain pairs. Kernel Average Misorientation (KAM) maps also showed damage accumulation on one side of the interface. Both of these are consistent with poor slip accommodation from a crystallographically softer grain to a harder one. Transmission electron microscopy was used to confirm the EBSD observations and to show the role of slip bands in the development of large, interfacial stress concentrations. A viewgraph presentation accompanies the provided abstract.

Crooks, Roy↗

DOE OSTI2025

A Linear-Complexity Tensor Butterfly Algorithm for Compressing High-Dimensional Oscillatory Integral Operators

This paper presents a multilevel tensor compression algorithm called tensor butterfly algorithm for efficiently representing large-scale and high-dimensional oscillatory integral operators, including Green's functions for wave equations and integral transforms such as Radon transforms and Fourier transforms. The proposed algorithm leverages a tensor extension of the so-called complementary low-rank property of existing matrix butterfly algorithms. The algorithm partitions the discretized integral operator tensor into subtensors of multiple levels and factorizes each subtensor at the middle level as a Tucker-type interpolative decomposition, whose factor matrices are formed in a multilevel fashion. For a d-dimensional (d > 1) integral operator discretized into a 2d-mode tensor with n2d entries, the overall CPU time and memory requirement scale as O(nd), in stark contrast to the O(nd log n) complexity of existing matrix algorithms such as matrix butterfly algorithms and fast Fourier transforms (FFTs), where n is the number of points per direction. When comparing with other tensor algorithms such as quantized tensor train (QTT), the proposed algorithm also shows superior CPU and memory performance for tensor contraction. Remarkably, the tensor butterfly algorithm can efficiently model high-frequency Green's function interactions between two unit cubes, each spanning 512 wavelengths per direction, which represents problems of scale over 512× larger than that existing butterfly algorithms can handle, with the same amount of computation resources. On the other hand, for a problem representing 64 wavelengths per direction, which is the largest size existing algebraic matrix algorithms can handle, our tensor butterfly algorithm exhibits 200x speedups and 30× memory reduction compared with existing ones. Moreover, the tensor butterfly algorithm also permits O(nd)-complexity FFTs and Radon transforms up to d = 6 dimensions.

Kielstra, P Michael↗

DOE OSTI2011

Visser's massive graviton bimetric theory revisited

A massive gravity theory was proposed by Visser in the late 1990s. This theory, based on a background metric b{sub {alpha}{beta}} and on an usual dynamical metric g{sub {alpha}{beta}} has the advantage of being free of ghosts as well as discontinuities present in other massive theories proposed in the past. In the present investigation, the equations of Visser's theory are revisited with particular care on the related conservation laws. It will be shown that a multiplicative factor is missing in the graviton tensor originally derived by Visser, which has no incidence on the weak field approach but becomes important in the strong field regime when, for instance, cosmological applications are considered. In this case, contrary to some previous claims found in the literature, we conclude that a nonstatic background metric is required in order to obtain a solution able to mimic the {Lambda}CDM cosmology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

An experimental program with high duty-cycle polarized and unpolarized positron beams at Jefferson Lab

Positron beams, both polarized and unpolarized, are identified as important ingredients for the experimental programs at the next generation of lepton accelerators. In the context of the hadronic physics program at Jefferson Lab (JLab), positron beams are complementary, even essential, tools for a precise understanding of the electromagnetic structure of nucleons and nuclei, in both the elastic and deep-inelastic regimes. For instance, elastic scattering of polarized and unpolarized electrons and positrons from the nucleon enables a model independent determination of its electromagnetic form factors. Also, the deeply-virtual scattering of polarized and unpolarized electrons and positrons allows unambiguous separation of the different contributions to the cross section of the lepto-production of photons and of lepton-pairs, enabling an accurate determination of the nucleons and nuclei generalized parton distributions, and providing an access to the gravitational form factors of the energy-momentum tensor. Furthermore, positron beams offer the possibility of alternative tests of the Standard Model of particle physics through the search of a dark photon, the precise measurement of electroweak couplings, and the investigation of charged lepton flavor violation. This document discusses the perspectives of an experimental program with high duty-cycle positron beams at JLab.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2013

Physics-based seismic hazard analysis on petascale heterogeneous supercomputers, In: SC '13 Proceedings of the International Conference on High Performance Computing, Networking, Storage and Analysis

We have developed a highly scalable and efficient GPU-based finite-difference code (AWP) for earthquake simulation that implements high throughput, memory locality, communication reduction and communication/computation overlap and achieves linear scalability on Cray XK7 Titan at ORNL and NCSA's Blue Waters system. We simulate realistic 0-10 Hz earthquake ground motions relevant to building engineering design using high-performance AWP. Moreover, we show that AWP provides a speedup by a factor of 110 in key strain tensor calculations critical to probabilistic seismic hazard analysis (PSHA). These performance improvements to critical scientific application software, coupled with improved co-scheduling capabilities of our workflow-managed systems, make a statewide hazard model a goal reachable with existing supercomputers. The performance improvements of GPU-based AWP are expected to save millions of core-hours over the next few years as physics-based seismic hazard analysis is developed using heterogeneous petascale supercomputers.

Computer Science↗

DOE OSTI2013

Chiral Dynamics and Transverse Peripheral Densities

We study the transverse charge and current densities in the periphery of the nucleon where they are governed by universal chiral dynamics. The large?distance asymptotics of the densities are derived from the spectral functions of the Dirac and Pauli form factors near the threshold at t = 4M?2 obtained from invariant chiral perturbation theory. In the equivalent light?front formulation, the peripheral densities arise from large?size ?N configurations in the nucleon?s light?front wave function obtained from the chiral Lagrangian. We compute the chiral component of the charge and current densities and explain their order?of?magnitude in the chiral expansion. We also discuss the role of ?? intermediate states in the large Nc limit and the heavy?baryon expansion of peripheral densities, and comment on progress extending an analogous analysis concerning transverse densities from form factors of the energy-momentum tensor and its connections with studies on GPDs.

Granados, Carlos↗

DOE OSTI2025

A two-loop four-point form factor at function level

Recently, the maximally-helicity-violating four-point form factor for the chiral stress-energy tensor in planar $\mathcal{N}$ = 4 super Yang-Mills was computed to three loops at the level of the symbol associated with multiple polylogarithms. It exhibits antipodal self-duality, or invariance under the combined action of a kinematic map and reversing the ordering of letters in the symbol. Here we lift the two-loop form factor from symbol level to function level. We provide an iterated representation of the function’s derivatives (coproducts). In order to do so, we find a three-parameter limit of the five-parameter phase space where the symbol’s letters are all rational. We also use function-level information about dihedral symmetries and the soft, collinear, and factorization limits, as well as limits governed by the form-factor operator product expansion (FFOPE). We provide plots of the remainder function on several kinematic slices, and show that the result is compatible with the FFOPE data. We further verify that antipodal self-duality is valid at two loops beyond the level of the symbol.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

A note on NMHV form factors from the Graßmannian and the twistor string

In this note we investigate Graßmannian formulas for form factors of the chiral part of the stress-tensor multiplet in N = 4 superconformal Yang-Mills theory. We present an all-n contour for the G(3, n + 2) Graßmannian integral of NMHV form factors derived from on-shell diagrams and the BCFW recursion relation. In addition, we study other G(3, n + 2) formulas obtained from the connected prescription introduced recently. We find a recursive expression for all n and study its properties. For n ≥ 6, our formula has the same recursive structure as its amplitude counterpart, making its soft behaviour manifest. Finally, we explore the connection between the two Graßmannian formulations, using the global residue theorem, and find that it is much more intricate compared to scattering amplitudes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗