Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tensor Approximation”

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 55 records · Page 3

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↗

The Average Spectrum Norm and Near-Optimal Tensor Completion

We propose the average spectrum norm to study the minimum number of measurements required to approximate a multidimensional array (i.e., sample complexity) via low-rank tensor recovery. Our focus is on the tensor completion problem, where the aim is to estimate a multiway array using a subset of tensor entries corrupted by noise. Our average spectrum norm-based analysis provides near-optimal sample complexities, exhibiting dependence on the ambient dimensions and rank that do not suffer from exponential scaling as the order increases.

97 MATHEMATICS AND COMPUTING↗

A quadratic-tensor model algorithm for nonlinear least-squares problems with linear constraints

A new algorithm for solving nonlinear least-squares and nonlinear equation problems is proposed which is based on approximating the nonlinear functions using the quadratic-tensor model by Schnabel and Frank. The algorithm uses a trust region defined by a box containing the current values of the unknowns. The algorithm is found to be effective for problems with linear constraints and dense Jacobian matrices.

Hanson, R. J.↗

High Temperature Elastic Properties of Single Crystal Mullite (Approximately 2.5Al2O3.SiO2) by Brillouin Spectroscopy

The complete elastic tensor of mullite has been determined by brillouin spectroscopy at room temperature and elevated temperatures up to 1200C. Equivalent, isotropic moduli (bulk, shear, and Young's) have been calculated. The room temperature values obtained using Voigt-Reuss-Hill averaging are: K(sub VRH) = 173.5 + 6.9 GPa, G(sub VRH) = 88.0 + 3.5 GPa, E(sub VRH) = 225.9 + 9.0 GPa. All moduli show relatively gradual decreases with temperature. The temperature derivatives obtained for the equivalent, isotropic moduli are: dK(sub VRH)/dT = - 17.5 + 2.5 MPa/deg. C, dG(sub VRH)/dT = -8.8 + 1.4 MPa/deg. C, dE(sub VRH)/dT = -22.6 + 2.8 MPa/deg C. Substantial differences between bulk properties calculated from the single crystal measurements in this study and the properties reported in the literature for polycrystalline sintered mullite are identified, indicating the importance of factors such as microstructure, intergranular phases, and composition to the elasticity of mullite ceramics.

Palko, James W.↗

Geodesic synchrotron radiation in the Kerr geometry by the method of asymptotically factorized Green's functions

The scalar, electromagnetic, and gravitational geodesic-synchrotron-radiation (GSR) spectra are determined for the case of a test particle moving on a highly relativistic circular orbit about a rotating (Kerr) black hole. It is found that the spectral shape depends only weakly on the value of the angular-momentum parameter (a/M) of the black hole, but the total radiated power drops unexpectedly for a value of at least 0.95 and vanishes as the value approaches unity. A spin-dependent factor (involving the inner product of the polarization of a radiated quantum with the source) is isolated to explain the dependence of the spectral shape on the spin of the radiated field. Although the scalar wave equation is solved by separation of variables, this procedure is avoided for the vector and tensor cases by postulating a sum-over-states expansion for the Green's function similar to that found to hold in the scalar case. The terms in this sum, significant for GSR, can then be evaluated in the geometric-optics approximation without requiring the use of vector or tensor spherical harmonics.

Chrzanowski, P. L.↗

Efficient method for approximating nonlinear dynamics: applications to uncertainty propagation and estimation

High-order Taylor series expansions can be used to model nonlinear dynamics at the cost of integrating a large set of variational equations to obtain high-order state-transition tensors (STTs). This paper presents an innovative technique for approximating the high-order STTs that reduces significantly the computational cost by retaining only the dominant secular terms. We propagate the low-order partial derivatives of Kepler’s equation, which only requires the integration of six additional equations to extend an n-th order approximation to order (n + 1). The approximation stems from the Lindstedt-Poincare procedure and exploits the stability properties of orbital motion. Since the method makes no dynamical assumptions, it can accommodate any source of orbital perturbations. We show how the approximation of the second-order STT significantly increases the accuracy of the linear method for uncertainty propagation with only a small computational overhead. Finally, we derive a high-order approximate extended Kalman filter that implements the proposed approximation of the STT and improves the performance of linear filters. Examples of application with different perturbation sources include the heliocentric orbit of an asteroid, an orbiter around Europa, and an Earth-orbiting satellite.

Park, Ryan S↗

Lowering the Scaling of Self-Consistent Field Methods by Combining Tensor Hypercontraction and a Density Difference Ansatz

We present the tensor hypercontraction difference self-consistent field (SCF) method, an approach that reduces the formal computational scaling of traditional naive self-consistent field methods from 𝑂(𝑁 4 ) to 𝑂(𝑁 3 ) with system size 𝑁. The scaling reduction is achieved by developing a new technique for constructing the tensor hypercontraction decomposition based on the fundamental approximation made in density fitting. Combining this scheme with the difference self-consistent field methodology, we achieve a method that enables 𝑂(𝑁 3 ) scaling SCF calculations with only 𝑂(𝑁 2 ) storage requirements. In conclusion, our proof-of-concept numerical tests demonstrate robust performance with errors in total energies below 8 × 10 –4 E h and with sub 1 kcal/mol errors for relative energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Randomized Algorithms for Low-Rank Matrix and Tensor Decompositions

This paper surveys randomized algorithms in numerical linear algebra for low-rank decompositions of matrices and tensors. The survey begins with a review of classical matrix algorithms that can be accelerated by randomized dimensionality reduction, such as the singular value decomposition (SVD) or interpolative (ID) and CUR decompositions. Recent advances in randomized dimensionality reduction are discussed, including new methods of fast matrix sketching and sampling techniques, which are incorporated into classical matrix algorithms for fast low-rank matrix approximations. The extension of randomized matrix algorithms to tensors is then explored for several low-rank tensor decompositions in the CP and Tucker formats, including the higher-order SVD, ID, and CUR decomposition.

Pearce, Katherine J. [The University of Texas at A↗

Ten-moment fluid modeling of the Weibel instability

We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.

Kuldinow, D. A. (ORCID:0000000319730196)↗

Probing the Isospin Composition of Short-Range Correlated Pairs at Jefferson Lab Hall B

Nucleons in short-range correlated (SRC) pairs, due to their close proximity and high relative momentum, can provide insight into the short-range part of the strong nuclear interaction. In particular, the prevalence of np pairs is due to the dominance of a tensor term for correlated nucleons with momenta of approximately 400?600 MeV/c. This dissertation comprises two studies advancing the community?s understanding of the isospin composition of SRC pairs. First, I performed a study of proton and neutron knockout from initially low-momentum and high-momentum states in 3He. Previous work has shown that protons are disproportionately represented in high-momentum states in neutron-rich nuclei. I demonstrate that spectral functions for the proton-rich nucleus 3He predict, in agreement with data, that neutrons are disproportionately represented in high-momentum states, but that 3He does not display the same strong prevalence of np pairs that is observed in larger nuclei. Second, Generalized Contact Formalism (GCF), a well-supported theory for predicting SRC behavior, predicts the transition from an isospin-dependent, tensor-dominant interaction at intermediate distances to a scalar-dominant, isospin-independent interaction at very short distances. This dissertation uses data from the CLAS12 Nuclear Targets Experiments in Hall B at Jefferson Lab to measure the relative abundances of pp and pn pairs for increasing relative momentum and decreasing separation. I provide an independent confirmation of the previously-observed increase in pp pairs at increasing momentum of the struck nucleon. I also contribute to the application of the new CLAS12 Central Neutron Detector by precisely measuring the neutron detection efficiency and developing a machine learning model for rejecting charged particle background.

Seroka, Erin↗

Tensor network representation of non-abelian gauge theory coupled to reduced staggered fermions

We show how to construct a tensor network representation of the path integral for reduced staggered fermions coupled to a non-abelian gauge field in two dimensions. The resulting formulation is both memory and computation efficient because reduced staggered fermions can be represented in terms of a minimal number of tensor indices while the gauge sector can be approximated using Gaussian quadrature with a truncation. Numerical results obtained using the Grassmann TRG algorithm are shown for the case of SU(2) lattice gauge theory and compared to Monte Carlo results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Covariant Quantum Error-Correcting Codes with Metrological Entanglement Advantage

Here, we show that a subset of the basis for the irreducible representations of a tensor-product SU(2) rotation forms a covariant approximate quantum error-correcting code with transversal U(1) logical gates. Generalizing previous work on “thermodynamic codes” to general local spin and different irreducible representations using only properties of the angular momentum algebra, we obtain bounds on the code inaccuracy under generic noise on any known 𝑑 sites, under independent and identically distributed noise, and under heralded 𝑑-local erasures. We demonstrate that this family of codes protects a probe state with quantum Fisher information surpassing the standard quantum limit when the sensing parameter couples to the generator of the U(1) logical gate.

quantum error correction↗

Sound production due to large-scale coherent structures

The acoustic pressure fluctuations due to large-scale finite amplitude disturbances in a free turbulent shear flow are calculated. The flow is decomposed into three component scales; the mean motion, the large-scale wave-like disturbance, and the small-scale random turbulence. The effect of the large-scale structure on the flow is isolated by applying both a spatial and phase average on the governing differential equations and by initially taking the small-scale turbulence to be in energetic equilibrium with the mean flow. The subsequent temporal evolution of the flow is computed from global energetic rate equations for the different component scales. Lighthill's theory is then applied to the region with the flowfield as the source and an observer located outside the flowfield in a region of uniform velocity. Since the time history of all flow variables is known, a minimum of simplifying assumptions for the Lighthill stress tensor is required, including no far-field approximations. A phase average is used to isolate the pressure fluctuations due to the large-scale structure, and also to isolate the dynamic process responsible. Variation of mean square pressure with distance from the source is computed to determine the acoustic far-field location and decay rate, and, in addition, spectra at various acoustic field locations are computed and analyzed. Also included are the effects of varying the growth and decay of the large-scale disturbance on the sound produced.

Gatski, T. B.↗

Development and testing of stable, invariant, isoparametric curvilinear 2- and 3-D hybrid-stress elements

Linear and quadratic Serendipity hybrid-stress elements are examined in respect of stability, coordinate invariance, and optimality. A formulation based upon symmetry group theory successfully addresses these issues in undistorted geometries and is fully detailed for plane elements. The resulting least-order stable invariant stress polynomials can be applied as astute approximations in distorted cases through a variety of tensor components and variational principles. A distortion sensitivity study for two- and three-dimensional elements provides favorable numerical comparisons with the assumed displacement method.

Punch, E. F.↗

Higher-Order Approximations for Stabilizing Zero-Energy Modes in Peridynamics Crystal Plasticity Models with Large Horizon Interactions

The non-ordinary state-based peridynamics theory combines non-local dynamic techniques with a desirable correspondence material principle, allowing for the use of continuum mechanics constitutive models. Such an approach presents a unique capability for solving problems involving discontinuities (e.g., strain localization, fracture, and fragmentation). However, the correspondence-based peridynamics models often suffer from zero-energy mode instabilities in numerical implementation, primarily due to the approximations of the non-local deformation gradient tensor. This paper focuses on a computational scheme for eliminating the zero-energy mode oscillations using a choice of influence functions that improve the truncation error in a higher-order Taylor series expansion of the deformation gradient. The novelty here is a tensor-based derivation of the linear constraint equations, which can be used to systematically identify the particle interaction weight functions for various user-specified horizon radii. In this paper, the proposed higher-order stabilization scheme is demonstrated for multi-dimensional examples involving polycrystalline and composite microstructures, along with comparisons against conventional finite element methods. The proposed stabilization scheme is shown to be highly effective in suppressing the spurious zero-energy mode oscillations in all numerical examples while enabling efficient simulations of strain localizations across material interfaces.

Non-Ordinary State-Based Peridynamics↗

Power Flow Geometry and Approximation

Here, the power flow equations are important in numerous power systems problems of practical interest which consider alternating current power flow (ACPF) physics. Perhaps the most well studied being the alternating current optimal power flow problem (ACOPF), seeking to optimize the operation of an electric power system. Due to their non-linearity, problems which include the power flow equations are typically challenging, particularly in optimization. Interestingly, the set of solutions to the power flow equations forms a smooth manifold. As a result, differential geometry can be used to describe and analyze this set of equations. This approach has proven effective in several engineering applications (e.g., solving ACOPF and analyzing the solution space boundary). Central to the success of this approach is an understanding of the power flow manifold's geometry. In this work, we develop the geometric and topological properties of this manifold using concepts from differential geometry. After demonstrating the convenience of this manifold's representation as a function's graph, computational methods are emphasized: we develop retractions, error bounds for linear approximation, and formulas for evaluating the Riemannian metric (including associated objects such as geodesics and the curvature tensor). Scalar curvature and the second fundamental form play a new role in quantifying the quality of linear approximations, like the popular direct current approximation. All functions are implemented in Julia and available in an online repository. Proofs are included for completeness.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Lectures on the scattering of light

The exact (Mie) theory for the scattering of a plane wave by a dielectric sphere is presented. Since this infinite series solution is computationally impractical for large spheres, another formulation is given in terms of an integral equation valid for a bounded, but otherwise general array of scatterers. This equation is applied to the scattering by a single sphere, and several methods are suggested for approximating the scattering cross section in closed form. A tensor scattering matrix is introduced, in terms of which some general scattering theorems are derived. The application of the formalism to multiple scattering is briefly considered.

Saxon, D. S.↗

Remote sensing of Earth terrain

The layered random medium model is used to investigate the fully polarimetric scattering of electromagnetic waves from vegetation. The vegetation canopy is modeled as an anisotropic random medium containing nonspherical scatterers with preferred alignment. The underlying medium is considered as a homogeneous half space. The scattering effect of the vegetation canopy are characterized by 3-D correlation functions with variances and correlation lengths respectively corresponding to the fluctuation strengths and the physical geometries of the scatterers. The strong fluctuation theory is used to calculate the anisotropic effective permittivity tensor of the random medium and the distorted Born approximation is then applied to obtain the covariance matrix which describes the fully polarimetric scattering properties of the vegetation field. This model accounts for all the interaction processes between the boundaries and the scatterers and includes all the coherent effects due to wave propagation in different directions such as the constructive and destructive interferences. For a vegetation canopy with low attenuation, the boundary between the vegetation and the underlying medium can give rise to significant coherent effects.

Kong, Jin AU↗