Search NASA⌕ Search

SEARCH · Search NASA

Results for “VECTOR SPACE”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Extracting SIMD Parallelism from Recursive Task-Parallel Programs

The pursuit of computational efficiency has led to the proliferation of throughput-oriented hardware, from GPUs to increasingly wide vector units on commodity processors and accelerators. This hardware is designed to efficiently execute data-parallel computations in a vectorized manner. However, many algorithms are more naturally expressed as divide-and-conquer, recursive, task-parallel computations. In the absence of data parallelism, it seems that such algorithms are not well suited to throughput-oriented architectures. This paper presents a set of novel code transformations that expose the data parallelism latent in recursive, task- parallel programs. These transformations facilitate straightforward vectorization of task-parallel programs on commodity hardware. We also present scheduling policies that maintain high utilization of vector resources while limiting space usage. Across several task-parallel benchmarks, we demonstrate both efficient vector resource utilization and substantial speedup on chips using Intel’s SSE4.2 vector units, as well as accelerators using Intel’s AVX512 units. We then show through rigorous sampling that, in practice, our vectorization techniques are effective for a much larger class of programs.

Recursive Programs, task parallelism, vectorizatio↗
DOE OSTI2021

Binning is Sinning (Supernova Version): The Impact of Self-calibration in Cosmological Analyses with Type Ia Supernovae

Recent cosmological analyses (e.g., JLA, Pantheon) of Type Ia supernovae (SNe Ia) have propagated systematic uncertainties into a covariance matrix and either binned or smoothed the systematic uncertainty vectors in redshift space. We demonstrate that systematic error budgets of these analyses can be improved by a factor of ~ 1.5 × with the use of unbinned and unsmoothed covariance matrices. To understand this, we employ a separate approach that simultaneously fits for cosmological parameters and additional self-calibrating scale parameters that constrain the size of each systematic. We show that the covariance-matrix approach and scale-parameter approach indeed yield equivalent results, implying that in both cases the data can self-calibrate certain systematic uncertainties, but that this ability is hindered when information is binned or smoothed in redshift space. We review the top systematic uncertainties in current analyses and find that the reduction of systematic uncertainties in the unbinned case depends on whether a systematic is solely degenerate with the cosmological model in redshift space or whether it can be described by additional correlations between supernova properties and luminosity. Furthermore, we show that the power of self-calibration increases with the size of the data set, which presents a tremendous opportunity for upcoming analyses of photometrically classified samples, like those of Legacy Survey of Space and Time (LSST) and the Nancy Grace Roman Telescope (NGRST). However, to take advantage of self-calibration in large, photometrically classified samples, we must first address the issue that binning is required in currently used photometric analysis methods.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

ModuleOT: A Hardware Security Module for Operational Technology: Preprint

With increasing penetration levels of distributed energy resources (DERs) on the distribution grid, as well as new technological advancements in the cyber space, new cyberattack vectors are being introduced, and the available attack surface is constantly increasing. Despite this increasing risk, the standard IEEE 1547-2018 does not yet recommend cybersecurity measures for DERs. To address this and to better protect data on the distribution grid - from the standpoints information security as well as operational security - ModuleOT has been developed. The module aims to significantly reduce cyberattack vectors by improving data privacy for user applications. This is accomplished by performing the core functions of encryption, authentication, authorization, certificate management, and user access control. The module integrates a custom security application with hardware cryptographic acceleration. The application secures all communications using Transmission Control Protocol over Internet Protocol (TCP/IP). These include the three most commonly used communications protocols for power systems information exchange: Modbus, Distributed Network Protocol 3 (DNP3), and Smart Energy Profile 2.0 (SEP2.0). These three protocols are also supported by IEEE 1547-2018 for all DER devices. This paper tests the data encryption/decryption feature on a physical networking test bed with emulated Modbus devices reporting grid data and presents the results.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2014

A comparison of four direct geometry time-of-flight spectrometers at the Spallation Neutron Source

The Spallation Neutron Source at Oak Ridge National Laboratory now hosts four direct geometry time-of-flight chopper spectrometers. These instruments cover a range of wave-vector and energy transfer space with varying degrees of neutron flux and resolution. The regions of reciprocal and energy space available to measure at these instruments are not exclusive and overlap significantly. We present a direct comparison of the capabilities of this instrumentation, conducted by data mining the instrument usage histories, and specific scanning regimes. In addition, one of the common science missions for these instruments is the study of magnetic excitations in condensed matter systems. We have measured the powder averaged spin wave spectra in one particular sample using each of these instruments, and use these data in our comparisons.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2017

Extended vector-tensor theories

Recently, several extensions of massive vector theory in curved space-time have been proposed in many literatures. In this paper, we consider the most general vector-tensor theories that contain up to two derivatives with respect to metric and vector field. By imposing a degeneracy condition of the Lagrangian in the context of ADM decomposition of space-time to eliminate an unwanted mode, we construct a new class of massive vector theories where five degrees of freedom can propagate, corresponding to three for massive vector modes and two for massless tensor modes. We find that the generalized Proca and the beyond generalized Proca theories up to the quartic Lagrangian, which should be included in this formulation, are degenerate theories even in curved space-time. Finally, introducing new metric and vector field transformations, we investigate the properties of thus obtained theories under such transformations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Towards space from Hilbert space: finding lattice structure in finite-dimensional quantum systems

We report that field theories place one or more degrees of freedom at every point in space. Hilbert spaces describing quantum field theories, or their finite-dimensional discretizations on lattices, therefore have large amounts of structure: they are isomorphic to the tensor product of a smaller Hilbert space for each lattice site or point in space. Local field theories respecting this structure have interactions which preferentially couple nearby points. The emergence of classicality through decoherence relies on this framework of tensor-product decomposition and local interactions. We explore the emergence of such lattice structure from Hilbert-space considerations alone. We point out that the vast majority of finite-dimensional Hilbert spaces cannot be isomorphic to the tensor product of Hilbert-space subfactors that describes a lattice theory. A generic Hilbert space can only be split into a direct sum corresponding to a basis of state vectors spanning the Hilbert space, we consider setups in which the direct sum is naturally decomposed into two pieces. We define a notion of direct-sum locality which characterizes states and decompositions compatible with Hamiltonian time evolution. Finally, we illustrate these notions for a toy model that is the finite-dimensional discretization of the quantum-mechanical double-well potential. We discuss their relevance in cosmology and field theory, especially for theories which describe a landscape of vacua with different spacetime geometries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Curvature and geometric modules of noncommutative spheres and tori

When considered as submanifolds of Euclidean space, the Riemannian geometry of the round sphere and the Clifford torus may be formulated in terms of Poisson algebraic expressions involving the embedding coordinates, and a central object is the projection operator, projecting tangent vectors in the ambient space onto the tangent space of the submanifold. In this note, we point out that there exist noncommutative analogues of these projection operators, which implies a very natural definition of noncommutative tangent spaces as particular projective modules. These modules carry an induced connection from Euclidean space, and we compute its scalar curvature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Photonic Crystals Based on Media with Arbitrary Anisotropy of Permittivity and Permeability

A general approach to analyzing eigenmodes in anisotropic and gyrotropic 3D photonic crystals based on dielectric and magnetic media is proposed. This approach is based on the representation of stationary macroscopic Maxwell equations in the operator form corresponding to the quantum-mechanical equation for a photon with spin s = 1. In these equations, the strengths of electric and magnetic fields are put into correspondence with state vectors in the complex Hilbert space. The permittivity and permeability serve as operators acting on these vectors. It is shown that the problem of determining the eigenmodes of a photonic crystal is reduced to searching for eigenvectors and eigenvalues of the Hermitian operator characterizing the spin–orbit interaction of a photon in the periodic anisotropic medium under study. It is proposed to use photon states with a certain wave vector (certain momentum) and a certain linear or circular spin polarization as a basis for the representation operator equations. One-dimensional photonic crystals are considered as an example. The influence of anisotropy and gyrotropy on the dispersion of eigenmodes in these crystals is investigated. The group velocity of the eigenmodes, momentum transferred by them, and spin angular momentum are analyzed for the case of gyrotropic media.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Counting topological interface modes using simplicial characteristic classes

A computational approach for predicting the number of topological interface modes (TIMs) in hermitian systems using the spectral flow—monopole correspondence is presented. The number of TIMs is determined by calculating the Chern number of a complex line bundle of local polarisation vectors over a phase space sphere surrounding a Weyl point. The Chern number is computed by constructing the simplicial first Chern class of a discrete vector bundle on a simplicial mesh. This approach is gauge invariant, derivative free, structure preserving, and robust to noise. The algorithm is shown to reproduce the expected number of TIMs for the case of equatorial fluid waves and the topological Langmuir cyclotron wave. The possibility of using this algorithm to analyse experimental measurements of bulk wave polarisations and predict the associated number of TIMs is explored in a synthetic example.

discrete vector bundles↗
DOE OSTI2015

SOLAR FLARE PREDICTION USING SDO/HMI VECTOR MAGNETIC FIELD DATA WITH A MACHINE-LEARNING ALGORITHM

We attempt to forecast M- and X-class solar flares using a machine-learning algorithm, called support vector machine (SVM), and four years of data from the Solar Dynamics Observatory's Helioseismic and Magnetic Imager, the first instrument to continuously map the full-disk photospheric vector magnetic field from space. Most flare forecasting efforts described in the literature use either line-of-sight magnetograms or a relatively small number of ground-based vector magnetograms. This is the first time a large data set of vector magnetograms has been used to forecast solar flares. We build a catalog of flaring and non-flaring active regions sampled from a database of 2071 active regions, comprised of 1.5 million active region patches of vector magnetic field data, and characterize each active region by 25 parameters. We then train and test the machine-learning algorithm and we estimate its performances using forecast verification metrics with an emphasis on the true skill statistic (TSS). We obtain relatively high TSS scores and overall predictive abilities. We surmise that this is partly due to fine-tuning the SVM for this purpose and also to an advantageous set of features that can only be calculated from vector magnetic field data. We also apply a feature selection algorithm to determine which of our 25 features are useful for discriminating between flaring and non-flaring active regions and conclude that only a handful are needed for good predictive abilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Faxén form of time-domain force on a sphere in unsteady spatially varying viscous compressible flows

An explicit expression for the time-dependent force on a stationary, finite-sized spherical particle located in an unsteady inhomogeneous ambient flow is presented. The force expression accounts for both viscous and compressible effects. Towards this end, a time-harmonic plane travelling wave of a given frequency propagating in a viscous compressible flow over a sphere is considered. Linearized compressible Navier–Stokes equations are solved to obtain an analytical expression for the force exerted on the particle in the frequency domain. The force obtained in the Laplace space due to a travelling wave of a given frequency and wavenumber is then generalized to any arbitrary incoming flow. This is achieved by relating the radial and tangential velocity components in the Laplace space to the surface-averaged radial velocity and volume-averaged velocity vectors respectively in the time space. Moreover an expression relating the surface-averaged radial velocity and volume-averaged velocity vector has been provided. The total force is written as a summation of the undisturbed and disturbed force (quasi-steady, inviscid-unsteady and viscous-unsteady) contributions. The force contributions thus obtained are expressed as comprising of two parts – that arising due to spatial variation in the ambient flow and the other arising due to temporal variation. The current formulation is applicable to inhomogeneous ambient flows, however in the limit of negligible Reynolds and Mach numbers. The results are applicable even for particles of sizes larger than the acoustic wavelength. The accuracy of the explicit time-domain force expression is first tested by computing the force on an 80 mm diameter particle due to a weak planar expansion fan. Extension of this formulation when nonlinear effects become important is also proposed and tested by considering strong expansion fans. The results thus obtained are compared against corresponding axisymmetric numerical simulations.

Mechanics↗
DOE OSTI2020

Subspace Methods in Multi-Parameter Seismic Full Waveform Inversion

In full waveform inversion (FWI) high-resolution subsurface model parameters are sought. FWI is normally treated as a nonlinear least-squares inverse problem, in which the minimum of the corresponding misfit function is found by updating the model parameters. When multiple elastic or acoustic properties are solved for, simple gradient methods tend to confuse parameter classes. This is referred to as parameter cross-talk; it leads to incorrect model solutions, poor convergence and strong dependence on the scaling of the different parameter types. Determining step lengths in a subspace domain, rather than directly in terms of gradients of different parameters, is a potentially valuable approach to address this problem. The particular subspace used can be defined over a span of different sets of data or different parameter classes, provided it involves a small number of vectors compared to those contained in the whole model space. Additionally, in a subspace method, the basis vectors are defined first, and a local minimum is found in the space spanned by these. We examine the application of the subspace method within acoustic FWI in determining simultaneously updates for velocity and density. We first discuss the choice of basis vectors to construct the spanned space, from linear updates by distinguishing only the contributions of different parameter classes towards nonlinear updates by adding the contributions of higher-order perturbations of each parameter class. The numerical character of FWI solutions generated via subspace methods involving different basis vectors is then analyzed and compared with traditional FWI methods. The subspace methods can provide better reconstructions of the model, especially for the velocity, as well as improved convergence rates, while the computational costs are still comparable with the traditional FWI methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Space-Time Block Preconditioning for Incompressible Flow

Parallel-in-time methods have become increasingly popular in the simulation of time-dependent numerical PDEs, allowing for the efficient use of additional message passing interface processes when spatial parallelism saturates. Most methods treat the solution and parallelism in space and time separately. In contrast, all-at-once methods solve the full space-time system directly, largely treating time as simply another spatial dimension. All-at-once methods offer a number of benefits over separate treatment of space and time, most notably significantly increased parallelism and faster time to solution (when applicable). However, the development of fast, scalable all-at-once methods has largely been limited to time-dependent (advection-)diffusion problems. This paper introduces the concept of space-time block preconditioning for the all-at-once solution of incompressible flow. By extending well-known concepts of spatial block preconditioning to the space-time setting, we develop a block preconditioner whose application requires the solution of a space-time (advection-)diffusion equation in the velocity block, coupled with a pressure Schur complement approximation consisting of independent spatial solves at each time-step, and a space-time matrix-vector multiplication. The new method is tested on four classical models in incompressible flow. Finally, the results indicate perfect scalability in refinement of spatial and temporal mesh spacing, perfect scalability in nonlinear Picard iteration count when applied to a nonlinear Navier--Stokes problem, and minimal overhead in terms of number of preconditioner applications compared with sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The lattice extraction of the TMD soft function using the auxiliary field representation of the Wilson line

The transverse momentum dependent (TMD) soft function can be obtained by formulating the Wilson line in terms of auxiliary 1-dimensional fermion fields on the lattice. In this formulation, the directional vector of the auxiliary field in Euclidean space has the form ˜ 𝑛 = ¹𝑖𝑛0– ®0?– 𝑛3º, where the time component is purely imaginary. The components of these complex directional vectors in the Euclidean space can be mapped directly to the rapidities of the Minkowski space soft function. We present the results of the one-loop calculation of the Euclidean space analog to the soft function using these complex directional vectors. As a result, we show that the calculation is valid only when the directional vectors obey the relation: j𝑟 j = j𝑛3𝑛0j ¡ 1, and that this result corresponds to a computation in Minkowski space with space-like directed Wilson lines. Finally, we show that a lattice calculable object can be constructed that has the desired properties of the soft function.

Francis, Anthony↗
DOE OSTI2024

Kernelized approaches to streaming compression of scientific data

In this paper three algorithms are developed for the streaming compression of scientific data. The algorithms presented are reliant on the theory of vector-valued reproducing kernel Hilbert spaces and operator valued kernel. Further, the scientific data is modeled as a snapshot of time dependent vector field F(x, t) over a manifold M and the recovery of the data is framed as a learning problem. These processes are then appropriately modified and ana lyzed for the streaming scenario in which data is generated without the ability to revisit past entries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Geometric Analysis Based Double Closed-loop Iterative Learning Control of Output PDF Shaping of Fiber Length Distribution in Refining Process

In order to improve the pulp quality and to reduce the energy consumption, the fiber length distribution (FLD) is generally employed as one of the important technological indexes in the refining process. Considering that the traditional mean and variance of fiber length are unable to adequately characterize the non-Gaussian distribution properties, this paper proposes a novel geometric analysis based double closed-loop iterative learning control (ILC) method for probability density function (PDF) shaping of output FLD in the refining process. Primarily, a RBF neural network (NN) with Gaussian-type is utilized to approximate the square root PDF in the inner loop, where the RBF basis function parameters (center and width) are tuned between any two adjacent batches by using an ILC law, and the subspace identification method can be applied to establish the state-space model of weight vector. Then, for the sake of accelerating the convergence rate of the closed-loop sysyem, a geometric analysis based ILC method is adopted in the outer loop. Lastly, both simulation and experiments demonstrate the effectiveness and practicability of the proposed approach.

42 ENGINEERING↗
DOE OSTI2015

Inflation with massive vector fields

We investigate the coupling between the inflaton and massive vector fields. All renormalizable couplings with shift symmetry of the inflaton are considered. The massive vector can be decomposed into a scalar mode and a divergence-free vector mode. We show that the former naturally interacts with the inflaton and the latter decouples at tree level. The model in general predicts f{sub NL}equil = O(1), while in some regions of the parameter space large non-Gaussianity can arise.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗