Search NASA⌕ Search

SEARCH · Search NASA

Results for “vector spaces”

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

Medial axis and local thickness computation using the Fast Sweeping Method

This report describes an efficient and robust voxel-based methodology for computing the medial axis, local thickness, and distance-to-skeleton of arbitrary three-dimensional geometries. It is assumed that the object can be represented by an exact or approximate signed distance function on a discrete grid. The gradient of such function is used to formulate a hyperbolic partial differential equation (PDE) that models the collapse of the position vector in space. By exploiting the causality property of the PDE, the Fast Sweeping Method is able to obtain the solution in a finite number of sweeps independent of the mesh resolution. The intersection of characteristic lines leads to the formation of shocks and a discrete bisector function is used to identify the medial axis. The same PDE approach is used to compute the local thickness inside the object and obtain the distance-to-skeleton field. Multiple examples are given in two and three dimensions along with a resolution study. The methodology has optimal complexity and yields subsecond computational times for geometries with over a million zones on a single core. The methodology is also capable of parallelization across shared and distributed memory architectures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Thrifty Array Format (TAF) file specifications

Thrifty Array Foram (TAF) files store numeric data in a binary format, minimizing storage requirements while preserving quick read access. Real data of any size and dimensionality can be stored in this format at varying degrees of numeric precision. Implicit array are associated with each dimension, eliminating the need to explicitly store uniformly-spaced grid vectors. Unlimited text comments may be included with the array for user documentation, and every file begins with a text synopsis of the binary structure. The format is deliberately designed for memory mapping, where portions of the array can be read without loading the entire file at once.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Theory and application of the vector pair correlation function for real-space crystallographic analysis of order/disorder correlations from STEM images

Deviations of local structure and chemistry from the average crystalline unit cell are increasingly recognized to have a significant influence on the properties of many technologically important materials. Here, we present the vector pair correlation function (vPCF) as a new real-space crystallographic analysis method, which can be applied to atomic-resolution scanning transmission electron microscopy (STEM) images to quantify and analyze structural order/disorder correlations. Our STEM-based vPCFs have several advantages over radial PCFs and/or 3D pair distribution functions from x-ray total scattering: vPCFs explicitly retain crystallographic orientation information, are spatially resolved, can be applied directly on a sublattice basis, and are suitable for any material that can be imaged with STEM. To show the utility of our approach, we measure partial vPCFs in Ba 5 SmSn 3 Nb 7 O 30 (BSSN), a tetragonal tungsten bronze (TTB) structured complex oxide. Many TTBs are known to be classical or relaxor ferroelectrics, and these properties have been correlated with the presence of superlattice ordering. BSSN, specifically, exhibits relaxor behavior and an incommensurate structural modulation. From the vPCF data, we observe that, of the cation sites, only the Ba (A2) sublattice is structurally modulated. We then infer the local modulation vector and reveal a marked anisotropy in its correlation length. Finally, short-range correlated polar displacements on the B2 cation sites are observed. This work introduces the vPCF as a powerful real-space crystallography technique, which enables direct, robust quantification of short-to-long range order on a sublattice-specific basis and is applicable to a wide range of complex material types.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

GenMod: A generative modeling approach for spectral representation of PDEs with random inputs

Here, we propose a method for quantifying uncertainty in high-dimensional PDE systems with random parameters, where the number of solution evaluations is small. Parametric PDE solutions are often approximated using a spectral decomposition based on polynomial chaos expansions. For the class of systems we consider (i.e., high dimensional with limited solution evaluations) the coefficients are given by an underdetermined linear system in a regression formulation. This implies additional assumptions, such as sparsity of the coefficient vector, are needed to approximate the solution. Here, we present an approach where we assume the coefficients are close to the range of a generative model that maps from a low to a high dimensional space of coefficients. Our approach is inspired be recent work examining how generative models can be used for compressed sensing in systems with random Gaussian measurement matrices. Using results from PDE theory on coefficient decay rates, we construct an explicit generative model that predicts the polynomial chaos coefficient magnitudes. The algorithm we developed to find the coefficients, which we call GenMod, is composed of two main steps. First, we predict the coefficient signs using Orthogonal Matching Pursuit. Then, we assume the coefficients are within a sparse deviation from the range of a sign-adjusted generative model. This allows us to find the coefficients by solving a nonconvex optimization problem, over the input space of the generative model and the space of sparse vectors. We obtain theoretical recovery results for a Lipschitz continuous generative model and for a more specific generative model, based on coefficient decay rate bounds. We examine three high-dimensional problems and show that, for all three examples, the generative model approach outperforms sparsity promoting methods at small sample sizes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Femtosecond photoemission electron microscopy of surface plasmon polariton beam steering via nanohole arrays

Directional control over surface plasmon polariton (SPP) waves is a prerequisite for the development of miniaturized optical circuitry. Here, the efficacy of single and dual component SPP steering elements is explored through photoemission electron microscopy. Our imaging scheme relies on two-color photoemission and counter-propagating SPP generation, which collectively allow SPPs to be visualized in real space. Wave-vector difference mixing between two-dimensional arrays of nanoholes and photon momenta enables SPP steering with directionality governed by the array lattice constant. The magnitude of grating momentum required to achieve counter-propagation guarantees efficient steering over the full range of the negative half-space. In our dual component configuration, separate SPP generation and Bragg diffraction based steering optics are employed. We find that array Bragg planes principally influence the SPP angles through the array band structure, which allows us to visualize both positive and negative refractory waves.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Accelerating eigenvalue computation for nuclear structure calculations via perturbative corrections

Subspace projection methods utilizing perturbative corrections have been proposed for computing the lowest few eigenvalues and corresponding eigenvectors of large Hamiltonian matrices. In this paper, we build upon these methods and introduce the term Subspace Projection with Perturbative Corrections (SPPC) method to refer to this approach. We tailor the SPPC for nuclear many-body Hamiltonians represented in a truncated configuration interaction subspace, i.e., the no-core shell model (NCSM). We use the hierarchical structure of the NCSM Hamiltonian to partition the Hamiltonian as the sum of two matrices. The first matrix corresponds to the Hamiltonian represented in a small configuration space, whereas the second is viewed as the perturbation to the first matrix. Eigenvalues and eigenvectors of the first matrix can be computed efficiently. Because of the split, perturbative corrections to the eigenvectors of the first matrix can be obtained efficiently from the solutions of a sequence of linear systems of equations defined in the small configuration space. These correction vectors can be combined with the approximate eigenvectors of the first matrix to construct a subspace from which more accurate approximations of the desired eigenpairs can be obtained. We show by numerical examples that the SPPC method can be more efficient than conventional iterative methods for solving large-scale eigenvalue problems such as the Lanczos, block Lanczos and the locally optimal block preconditioned conjugate gradient (LOBPCG) method. The method can also be combined with other methods to avoid convergence stagnation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 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 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 OSTI2024

The Cauchy Problem for Boltzmann Bi-linear Systems: The Mixing of Monatomic and Polyatomic Gases

Abstract From a unified vision of vector valued solutions in weighted Banach spaces, this paper establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such models can be viewed as extensions of scalar Boltzmann binary elastic flows, as much as monatomic gas mixtures with disparate masses and single polyatomic gases, providing a unified approach for vector valued solutions in weighted Banach spaces. Novel aspects of this work include developing the extension of a general ODE theory in vector valued weighted Banach spaces, precise lower bounds for the collision frequency in terms of the weighted Banach norm, energy identities, angular or compact manifold averaging lemmas which provide coerciveness resulting into global in time stability, a new combinatorics estimate for p -binomial forms producing sharper estimates for the k -moments of bi-linear collisional forms. These techniques enable the Cauchy problem improvement that resolves the model with initial data corresponding to strictly positive and bounded initial vector valued mass and total energy, in addition to only a $$2^+$$ 2 + moment determined by the hard potential rates discrepancy, a result comparable in generality to the classical Cauchy theory of the scalar homogeneous Boltzmann equation.

Physics↗
DOE OSTI2022

A detailed analysis of the dynamics of fast neutrino flavor conversions with scattering effects

We calculate the dynamics of fast neutrino flavor conversions with Boltzmann collisions of neutrino scatterings in a homogeneous system. We find the enhancement or suppression of the flavor conversions in various setups of the collision terms. We analyze the mechanism of fast flavor conversions based on the motion of polarization vectors in the cylindrical coordinate analogous to a pendulum motion. The phases of the all the polarization vectors synchronize in the linear evolution phase, and the phase deviation from the Hamiltonian governs the neutrino flavor conversions. In a non-linear regime of flavor conversions, the collision terms induce a spiral motion of the polarization vector and gradually make the phase space smaller. The collision terms align all of the polarization vectors, and the flavor conversions eventually settle into equilibrium when the distributions of neutrinos become isotropic. Though our current analysis does not fully clarify the non-linear phenomena of fast flavor conversions, the framework of the pendulum motion gives a new insight into this complicated phenomenon that will be helpful in further studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

System and method for compact adaptive free space optical communications

The present disclosure relates to a free-space optical communications system that may make use of a controller and a digital micro mirror (DMM) assembly in communication with the controller. The DMM may have a plurality of independently controllable micromirror elements forming both an emitter and a receiver. A laser is included which is configured to generate an optical beam which may be applied to select ones of the first subpluralities of micromirror elements to generate a transmitted free space optical signal along a selected vector. A detector is included for receiving an incoming free space optical signal imaged by at least one of the micromirror elements.

42 ENGINEERING↗
DOE OSTI2024

Charm physics with overlap fermions on 2+1-flavor domain wall fermion configurations*

Decay constants of pseudoscalar mesons D, D s , η c , and vector mesons D*, D$^{*}_{s}$, J/ψ are determined from the N f = 2 + 1 lattice QCD at a lattice spacing a ~ 0.08 fm. For vector mesons, the decay constants defined by tensor currents are given in the $\overline{MS}$ scheme at 2 GeV. The calculation is performed on domain wall fermion configurations generated by the RBC-UKQCD collaborations and the overlap fermion action is used for the valence quarks. Comparing the current results with our previous results at a coarser lattice spacing a ~ 0.11 fm provides a better understanding of the discretization error. We obtain $f$$^{T}_{D^*_s}$($\overline{MS}$, 2 GeV)/f D$^{*}_{s}$ = 0.909(18) with a better precision than our previous result. Combining our f D$^{*}_{s}$ = 277(11) MeV with the total width of D$^{*}_{s}$ determined in a recent study gives a branching fraction 4.26(52) x 10 –5 for D$^{*}_{s}$ leptonic decay.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Photon topology

The topology of photons in vacuum is interesting because there are no photons with k = 0, creating a hole in momentum space. We show that while the set of all photons forms a trivial vector bundle $γ$ over this momentum space, the R and L photons form topologically nontrivial subbundles $γ±$ with first Chern numbers ∓2. In contrast, $γ$ has no linearly polarized subbundles, and there is no Chern number associated with linear polarizations. It is a known difficulty that the standard version of Wigner’s little group method produces singular representations of the Poincaré group for massless particles. By considering representations of the Poincaré group on vector bundles we obtain a version of Wigner’s little group method for massless particles which avoids these singularities. Here we show that any massless bundle representation of the Poincaré group can be canonically decomposed into irreducible bundle representations labeled by helicity, which in turn can be associated to smooth irreducible Hilbert space representations. This proves that the R and L photons are globally well defined as particles and that the photon wave function can be uniquely split into R and L components. This formalism offers a method of quantizing the electromagnetic field without invoking discontinuous polarization vectors as in the traditional scheme. We also demonstrate that the spin-Chern number of photons is not a purely topological quantity. Lastly, there has been an extended debate on whether photon angular momentum can be split into spin and orbital parts. Our work explains the precise issues that prevent this splitting. Photons do not admit a spin operator; instead, the angular momentum associated with photons’ internal degree of freedom is described by a helicity-induced subalgebra corresponding to the translational symmetry of $γ$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗