Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Nucleon isovector axial form factors

We present results for the isovector axial vector form factors obtained using thirteen 2 + 1 + 1 -flavor highly improved staggered quark (HISQ) ensembles generated by the MILC collaboration. The calculation of nucleon two- and three-point correlation functions has been done using Wilson-clover fermions. In the analysis of these data, we quantify the sensitivity of the results to strategies used for removing excited state contamination and invoke the partially conserved axial current relation between the form factors to choose between them. Our data driven analysis includes removing contributions from multihadron N π states that make significant contributions. Our final results are g A = 1.292 ( 53 ) stat ( 24 ) sys for the axial charge; g S = 1.085 ( 50 ) stat ( 103 ) sys and g T = 0.991 ( 21 ) stat ( 10 ) sys for the scalar and tensor charges; ⟨ r A 2 ⟩ = 0.439 ( 56 ) stat ( 34 ) sys fm 2 for the mean squared axial charge radius, g P * = 9.03 ( 47 ) stat ( 42 ) sys for the induced pseudoscalar charge; and g π N N = 14.14 ( 81 ) stat ( 85 ) sys for the pion-nucleon coupling. We also provide a parametrization of the axial form factor G A ( Q 2 ) over the range 0 ≤ Q 2 ≤ 1 GeV 2 for use in phenomenology and a comparison with other lattice determinations. We find that the various lattice data agree within 10% but are significantly different from the extraction of G A ( Q 2 ) from the ν -deuterium scattering data. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Results from the Atacama B-mode Search (ABS) experiment

We report the Atacama B-mode Search is an experiment designed to measure the cosmic microwave background polarization at large angular scales (0ℓ>4). It observes at 145 GHz from a site at 5,190 m elevation in northern Chile. The noise equivalent polarization temperature, or NEQ, is 41 μK√s. One of the unique features of ABS is its use of a rapidly rotating ambient-temperature half-wave plate (HWP) {as the first optical element}. {The HWP spins} at 2.55 Hz to modulate the incident polarized signal at frequencies above where instrument white noise dominates over atmospheric fluctuations and other sources of low-frequency noise. We report here on the analysis of data from a 2,400 deg 2 region of sky. We perform a blind analysis to reduce potential bias. After unblinding, we find agreement with the Planck TE and EE measurements on the same region of sky, {with a derived calibration factor of 00.89 ± 0.1}. We marginally detect polarized dust emission {(at 3.2 σ for EE and 2.2 σ for BB)} and give an upper limit on the tensor-to-scalar ratio of r<2.3 (95% confidence level) with the equivalent of 100 on-sky days of observation. Finally, we also present a new measurement of the polarization of Tau A and introduce new methods for calibration and data analysis associated with HWP-based observations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

The effects of gravity darkening on the ultraviolet continuum polarization produced by circumstellar disks

We investigate the effects of gravity darkening on the UV continuum polarization produced by an axisymmetric disk that surrounds a rapidly rotating star. Although the model is a single scattering approximation, we do include the effects of attenuation (electron scattering plus hydrogen bound-free absorption) by the disk, using an approach similar to that of Sobolev (1963). Because of the gravity darkening of the star and the attenuation within the disk, the radiation field is not axially symmetric about the radius vector. This implies that the polarization source functions are no longer provided by the finite disk depolarization factors of Cassinelli, Nordsieck, & Murison (1987), which are functions of the intensity moments in a spherically symmetric atmosphere. We reformulate the polarization source functions using generalized intensity moment tensors (J, H(sub i), K(sub ij)) that are valid for an arbitrary radiation field and envelope geometry. We find that the polarization source functions are simplest when using intensity moments in the observer's reference frame. On the other hand, the intensity moments are most easily evaluated in the stellar reference frame. Using the rotation transformation properties of the generalized intensity moments, we relate the observer's moments to those evaluated in the stellar reference frame. Our procedure for determining the polarization source functions thus merely involves choosing a set of Euler angles for the coordinate rotations, and then evaluating the associated rotation matrix. The geometrical complications of polarization transfer are thus reduced to obtaining a coordinate rotation matrix.

Bjorkman, J. E.↗
DOE OSTI2021

Systematic Generation of Moment Invariant Bases for 2D and 3D Tensor Fields

Moment invariants have been successfully applied to pattern detection tasks in 2D and 3D scalar, vector, and matrix valued data. However so far no flexible basis of invariants exists, i.e., no set that is optimal in the sense that it is complete and independent for every input pattern. Here, we prove that a basis of moment invariants can be generated that consists of tensor contractions of not more than two different moment tensors each under the conjecture of the set of all possible tensor contractions to be complete. This result allows us to derive the first generator algorithm that produces flexible bases of moment invariants with respect to orthogonal transformations by selecting a single non-zero moment to pair with all others in these two-factor products. Since at least one non-zero moment can be found in every non-zero pattern, this approach always generates a complete set of descriptors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

The lapse function in Friedmann—Lemaître–Robertson–Walker cosmologies

Highlights: • The Hubble flow is non-inertial for all cosmic equations-of-state other than p=−ρ∕3. • The Local Flatness Theorem is therefore violated when the lapse function is set equal to a constant. • Thus, the symmetries in the FLRW metric appear to be valid only when ρ+3p=0, producing an expansion factor a(t) t. - Abstract: The Friedmann–Lemaître–Robertson–Walker (FLRW) metric, the backbone of modern cosmology, is founded on the cosmological principle, which assumes homogeneity and isotropy throughout the cosmos. One of its simplifications is the choice of lapse function g{sub tt}=1, regardless of which stress–energy tensor T{sup μν} is used in Einstein’s field equations. It is sometimes argued that this selection is justified by gauge freedom, given that g{sub tt} in FLRW may be a function solely of t, not of the spatial coordinates, permitting a redefinition of the time. We show in this paper, however, that the comoving frame in the Hubble expansion is non inertial for all but a few special cases of the expansion factor a(t). Changing the gauge changes the frame of reference and cannot alter this property of the expansion profile, since it would of necessity reformat the metric using the coordinates of a non-comoving observer. We therefore suggest that the pre-selection of g{sub tt}=1, independently of the equation of state in the cosmic fluid, incorrectly avoids the time dilation that ought to be present relative to the actual free-falling frame when a{sup ̈}≠0.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Quantitative linkage between the stress at dislocation channel – Grain boundary interaction sites and irradiation assisted stress corrosion crack initiation

Localized deformation has emerged as a key factor in the crack initiation process for irradiated steels, as cracks are observed to nucleate preferentially at these sites. Using high resolution electron backscatter diffraction (HREBSD), the local stress tensor surrounding the dislocation channel-grain boundary interaction sites was quantified and coupled with fully determined grain boundary plane orientation information to determine, for the first time, the relationship between grain boundary normal stress and intergranular crack initiation in irradiated austenitic stainless steel. A Fe-13Cr-15Ni alloy was strained in simulated boiling water reactor, normal water chemistry after quantifying the residual stress tensor at discontinuous dislocation channel – grain boundary interaction sites where grain boundaries were determined to be well oriented with respect to the loading axis. Local stresses at the grain boundary were observed to reach magnitudes greater than 1.5 GPa at a distance of 200 nm from the intersection between the dislocation channel and the grain boundary. A pseudo-threshold stress of 0.9 GPa was measured, below which no cracking was observed. As the stress acting normal to the grain boundary increased above this value, the susceptibility to cracking increased with the cracking fraction reaching 100% at the high end of the stress range. This study shows for the first time that not only does intersection between discontinuous dislocation channels and grain boundaries result in peak local stresses, while the magnitude of the local tensile stress drives the crack initiation process.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Precise Measurement of Deuteron Tensor Analyzing Powers with BLAST

We report a precision measurement of the deuteron tensor analyzing powers T{sub 20} and T{sub 21} at the MIT-Bates Linear Accelerator Center. Data were collected simultaneously over a momentum transfer range Q=2.15-4.50 fm{sup -1} with the Bates Large Acceptance Spectrometer Toroid using a highly polarized deuterium internal gas target. The data are in excellent agreement with calculations in a framework of effective field theory. The deuteron charge monopole and quadrupole form factors G{sub C} and G{sub Q} were separated with improved precision, and the location of the first node of G{sub C} was confirmed at Q=4.19{+-}0.05 fm{sup -1}. The new data provide a strong constraint on theoretical models in a momentum transfer range covering the minimum of T{sub 20} and the first node of G{sub C}.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Mixed-precision iterative refinement using tensor cores on GPUs to accelerate solution of linear systems

Double-precision floating-point arithmetic (FP64) has been the de facto standard for engineering and scientific simulations for several decades. Problem complexity and the sheer volume of data coming from various instruments and sensors motivate researchers to mix and match various approaches to optimize compute resources, including different levels of floating-point precision. In recent years, machine learning has motivated hardware support for half-precision floating-point arithmetic. A primary challenge in high-performance computing is to leverage reduced-precision and mixed-precision hardware. We show how the FP16/FP32 Tensor Cores on NVIDIA GPUs can be exploited to accelerate the solution of linear systems of equations Ax = b without sacrificing numerical stability. The techniques we employ include multiprecision LU factorization, the preconditioned generalized minimal residual algorithm (GMRES), and scaling and auto-adaptive rounding to avoid overflow. We also show how to efficiently handle systems with multiple right-hand sides. On the NVIDIA Quadro GV100 (Volta) GPU, we achieve a 4×-5× performance increase and 5× better energy efficiency versus the standard FP64 implementation while maintaining an FP64 level of numerical stability.

GMRES↗
DOE OSTI2025

𝐷-shell mixing in light baryons and its effect on the orbital motion

The standard description of the nucleon in the nonrelativistic quark model is an 1S,L = 0 state without orbital motion. Yet, there are several indications from phenomenology that an admixture of states with nonzero orbital motion maybe substantial. In this paper we focus on the “second shell” of the nucleon excitations (D-shell), for which we give a modern description of the wave functions. We follow it by investigating what we call a “maximal mixing” scenario, assuming a hypothetical long-range tensor force. We give the explicit wave functions for all states, before and after mixing, and reassess many predictions such as the magnetic moments, the standard and transitional form-factors from the nucleon to N*. Unexpectedly, in this scenario we can reproduce the long-puzzling features of the Roper resonance N*(1440). But even in this extreme case, the admixture of the 1D,L = 2 state to a nucleon remains significantly smaller than expected from phenomenology.

Baryons↗
DOE OSTI2011

Quantum modifications to gravity waves in de Sitter spacetime

We treat a model in which tensor perturbations of de Sitter spacetime, represented as a spatially flat model, are modified by the effects of the vacuum fluctuations of a massless conformally invariant field, such as the electromagnetic field. We use the semiclassical theory of gravity with the expectation value of the conformal field stress tensor as a source. We first study the stability of de Sitter spacetime by searching for growing, spatially homogeneous modes, and conclude that it is stable within the limits of validity of the semiclassical theory. We next examine the modification of linearized plane gravity waves by the effects of the quantum stress tensor. We find a correction term which is of the same form as the original wave, but displaced in phase by {pi}/2, and with an amplitude which depends upon an initial time. The magnitude of this effect is proportional to the change in scale factor after this time. We discuss alternative interpretations of this time, but pay particular attention to the view that this is the beginning of inflation. So long as the energy scale of inflation and the proper frequency of the mode at the beginning of inflation are well below the Planck scale, the fractional correction is small. However, modes which are trans-Planckian at the onset of inflation can undergo a significant correction. The increase in amplitude can potentially have observable consequences through a modification of the power spectrum of tensor perturbations in inflationary cosmology. This enhancement of the power spectrum depends upon the initial time, and is greater for shorter wavelengths.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Utah FORGE 3-2514: A Strain Sensing Array to Characterize Deformation at the FORGE Site - Workshop Presentation

This is a presentation on the Strain Sensing Array to Characterize Deformation at the FORGE Site project by Clemson University, presented by Lawrence Murdoch. The project's objective was to evaluate the feasibility of measuring and interpreting tensor strain data to improve the performance of EGS. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Quantum Fourier transform revisited

Summary The fast Fourier transform (FFT) is one of the most successful numerical algorithms of the 20th century and has found numerous applications in many branches of computational science and engineering. The FFT algorithm can be derived from a particular matrix decomposition of the discrete Fourier transform (DFT) matrix. In this paper, we show that the quantum Fourier transform (QFT) can be derived by further decomposing the diagonal factors of the FFT matrix decomposition into products of matrices with Kronecker product structure. We analyze the implication of this Kronecker product structure on the discrete Fourier transform of rank‐1 tensors on a classical computer. We also explain why such a structure can take advantage of an important quantum computer feature that enables the QFT algorithm to attain an exponential speedup on a quantum computer over the FFT algorithm on a classical computer. Further, the connection between the matrix decomposition of the DFT matrix and a quantum circuit is made. We also discuss a natural extension of a radix‐2 QFT decomposition to a radix‐ d QFT decomposition. No prior knowledge of quantum computing is required to understand what is presented in this paper. Yet, we believe this paper may help readers to gain some rudimentary understanding of the nature of quantum computing from a matrix computation point of view.

Camps, Daan↗
DOE OSTI2023

Distributed out-of-memory NMF on CPU/GPU architectures

We propose an efficient distributed out-of-memory implementation of the non-negative matrix factorization (NMF) algorithm for heterogeneous high-performance-computing systems. The proposed implementation is based on prior work on NMFk, which can perform automatic model selection and extract latent variables and patterns from data. In this work, we extend NMFk by adding support for dense and sparse matrix operation on multi-node, multi-GPU systems. The resulting algorithm is optimized for out-of-memory problems where the memory required to factorize a given matrix is greater than the available GPU memory. Memory complexity is reduced by batching/tiling strategies, and sparse and dense matrix operations are significantly accelerated with GPU cores (or tensor cores when available). Input/output latency associated with batch copies between host and device is hidden using CUDA streams to overlap data transfers and compute asynchronously, and latency associated with collective communications (both intra-node and inter-node) is reduced using optimized NVIDIA Collective Communication Library (NCCL) based communicators. Benchmark results show significant improvement, from 32X to 76x speedup, with the new implementation using GPUs over the CPU-based NMFk. Good weak scaling was demonstrated on up to 4096 multi-GPU cluster nodes with approximately 25,000 GPUs when decomposing a dense 340 Terabyte-size matrix and an 11 Exabyte-size sparse matrix of density 10 -6 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

The inversion layer of electric fields and electron phase-space-hole structure during two-dimensional collisionless magnetic reconnection

Based on two-dimensional fully kinetic simulations that resolve the electron diffusion layer in undriven collisionless magnetic reconnection with zero guide field, this paper reports the existence and evolution of an inversion layer of bipolar electric fields, its corresponding phase-space structure (an electron-hole layer), and the implication to collisionless dissipation. The inversion electric field layer is embedded in the layer of bipolar Hall electric field and extends throughout the entire length of the electron diffusion layer. The electron phase-space hole structure spontaneously arises during the explosive growth phase when there exist significant inflows into the reconnection layer, and electrons perform meandering orbits across the layer while being cyclotron-turned toward the outflow directions. The cyclotron turning of meandering electrons by the magnetic field normal to the reconnection layer is shown to be a primary factor limiting the current density in the region where the reconnection electric field is balanced by the gradient (along the current sheet normal) of the off-diagonal electron pressure-tensor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

An orthogonal coordinate grid following the three-dimensional viscous flow over a concave surface

Swept wings designed for laminar flow control exhibit both centrifugal and crossflow instabilities which produce streamwise vortices that can lead to early transition from laminar to turbulent flow in the presence of Tollmien-Schlichting waves. This paper outlines an iterative algorithm for generation of an orthogonal, curvilinear, coordinate grid following the streamlines of the three-dimensional viscous flow over a swept, concave surface. The governing equations for the metric tensor are derived from the Riemann-Christoffel tensor for an Euclidian geometry. Unit vectors along streamline, normal and binormal directions are determined. The governing equations are not solved directly, but are employed only as compatibility equations. The scale factor for the streamline coordinate is obtained by an iterative integration scheme on a 200 x 100 x 5 grid, while the other two scale factors are determined from definitions. Sample results are obtained which indicate that the compatibility equation error decreases linearly with grid step size. Grids smaller than 200 x 100 x 5 are found to be inadequate to resolve the grid curvature.

Dagenhart, J. R↗
NASA NTRS1992

LES versus DNS: A comparative study

We have performed Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) of forced isotropic turbulence at moderate Reynolds numbers. The subgrid scale model used in the LES is based on an eddy viscosity which adjusts instantaneously the energy spectrum of the LES to that of the DNS. The statistics of the large scales of the DNS (filtered DNS field or fDNS) are compared to that of the LES. We present results for the transfer spectra, the skewness and flatness factors of the velocity components, the PDF's of the angle between the vorticity and the eigenvectors of the rate of strain, and that between the vorticity and the vorticity stretching tensor. The above LES statistics are found to be in good agreement with those measured in the fDNS field. We further observe that in all the numerical measurements, the trend was for the LES field to be more gaussian than the fDNS field. Future research on this point is planned.

Shtilman, L.↗
DOE OSTI2010

Kaluza-Klein two-brane-worlds cosmology at low energy

We study two (4+n)-dimensional branes embedded in (5+n)-dimensional spacetime. Using the gradient expansion approximation, we find that the effective theory is described by (4+n)-dimensional scalar-tensor gravity with a specific coupling function. Based on this theory we investigate the Kaluza-Klein two-brane-worlds cosmology at low energy, in both the static and the nonstatic internal dimensions. In the case of the static internal dimensions, the effective gravitational constant in the induced Friedmann equation depends on the equations of state of the brane matter, and the dark radiation term naturally appears. In the nonstatic case we take a relation between the external and internal scale factors of the form b(t)=a{sup {gamma}(t)} in which the brane world evolves with two scale factors. In this case, the induced Friedmann equation on the brane is modified in the effective gravitational constant and the term proportional to a{sup -4{beta}.} For dark radiation, we find {gamma}=-2/(1+n). Finally, we discuss the issue of conformal frames which naturally arises with scalar-tensor theories. We find that the static internal dimensions in the Jordan frame may become nonstatic in the Einstein frame.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Dark matter relic density in scalar-tensor gravity revisited

We revisit the calculation of dark matter relic abundances in scalar-tensor gravity using a generic form A(φ{sub *}) = e{sup βφ{sub *}{sup 2/2}} for the coupling between the scalar field φ{sub *} and the metric, for which detailed Big Bang Nucleosynthesis constraints are available. We find that BBN constraints restrict the modified expansion rate in these models to be almost degenerate with the standard expansion history at the time of dark matter decoupling. In this case the maximum level of enhancement of the dark matter relic density was found to be a factor of ∼ 3, several orders of magnitude below that found in previous investigations.

79 ASTRONOMY AND ASTROPHYSICS↗