Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “singular vectors”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

An Efficient and Robust Singular Value Method for Star Pattern Recognition and Attitude Determination

A new star pattern recognition method is developed using singular value decomposition of a measured unit column vector matrix in a measurement frame and the corresponding cataloged vector matrix in a reference frame. It is shown that singular values and right singular vectors are invariant with respect to coordinate transformation and robust under uncertainty. One advantage of singular value comparison is that a pairing process for individual measured and cataloged stars is not necessary, and the attitude estimation and pattern recognition process are not separated. An associated method for mission catalog design is introduced and simulation results are presented.

Juang, Jer-Nan↗
DOE OSTI2024

Efficient generation and amplification of intense vortex and vector laser pulses via strongly-coupled stimulated Brillouin scattering in plasmas

The past decade has seen tremendous progress in the production and utilization of vortex and vector laser pulses. Although both are considered as structured light beams, the vortex lasers have helical phase fronts and phase singularities, while the vector lasers have spatially variable polarization states and polarization singularities. In contrast to the vortex pulses that carry orbital angular momentum (OAM), the vector laser pulses have a complex spin angular momentum (SAM) and OAM coupling. Despite many potential applications enabled by such pulses, the generation of high-power/-intensity vortex and vector beams remains challenging. Here, we demonstrate using theory and three-dimensional simulations that the strongly-coupled stimulated Brillouin scattering (SC-SBS) process in plasmas can be used as a promising amplification technique with up to 65% energy transfer efficiency from the pump beam to the seed beam for both vortex and vector pulses. We also show that SC-SBS is strongly polarization-dependent in plasmas, enabling an all-optical polarization control of the amplified seed beam. Additionally, the interaction of such structured lasers with plasmas leads to various angular momentum couplings and decouplings that produce intense new light structures with controllable OAM and SAM. This scheme paves the way for novel optical devices such as plasma-based amplifiers and light field manipulators.

Physics↗
NASA NTRS2002

Development of an Efficient Binaural Simulation for the Analysis of Structural Acoustic Data

Binaural or "virtual acoustic" representation has been proposed as a method of analyzing acoustic and vibroacoustic data. Unfortunately, this binaural representation can require extensive computer power to apply the Head Related Transfer Functions (HRTFs) to a large number of sources, as with a vibrating structure. This work focuses on reducing the number of real-time computations required in this binaural analysis through the use of Singular Value Decomposition (SVD) and Equivalent Source Reduction (ESR). The SVD method reduces the complexity of the HRTF computations by breaking the HRTFs into dominant singular values (and vectors). The ESR method reduces the number of sources to be analyzed in real-time computation by replacing sources on the scale of a structural wavelength with sources on the scale of an acoustic wavelength. It is shown that the effectiveness of the SVD and ESR methods improves as the complexity of the source increases. In addition, preliminary auralization tests have shown that the results from both the SVD and ESR methods are indistinguishable from the results found with the exhaustive method.

Lalime, Aimee L.↗
NASA NTRS2003

Development of an Efficient Binaural Simulation for the Analysis of Structural Acoustic Data

Applying binaural simulation techniques to structural acoustic data can be very computationally intensive as the number of discrete noise sources can be very large. Typically, Head Related Transfer Functions (HRTFs) are used to individually filter the signals from each of the sources in the acoustic field. Therefore, creating a binaural simulation implies the use of potentially hundreds of real time filters. This paper details two methods of reducing the number of real-time computations required by: (i) using the singular value decomposition (SVD) to reduce the complexity of the HRTFs by breaking them into dominant singular values and vectors and (ii) by using equivalent source reduction (ESR) to reduce the number of sources to be analyzed in real-time by replacing sources on the scale of a structural wavelength with sources on the scale of an acoustic wavelength. The ESR and SVD reduction methods can be combined to provide an estimated computation time reduction of 99.4% for the structural acoustic data tested. In addition, preliminary tests have shown that there is a 97% correlation between the results of the combined reduction methods and the results found with the current binaural simulation techniques

Johnson, Marty E.↗
NASA NTRS2002

Geostrophic balance with a full Coriolis Force: implications for low latitutde studies

In its standard form, geostrophic balance uses a partial representation of the Coriolis force. The resulting formation has a singularity at the equator, and violates mass and momentum conservation. When the horizontal projection of the planetary rotation vector is considered, the singularity at the equator disappears, continuity can be preserved, and quasigeostrophy can be formulated at planetary scale.

balance atmosphere geostrophy remote sensing wind↗
DOE OSTI2023

Mallat Scattering Transformation based surrogate for Magnetohydrodynamics

Abstract A Machine and Deep Learning (MLDL) methodology is developed and applied to give a high fidelity, fast surrogate for 2D resistive MagnetoHydroDynamic (MHD) simulations of Magnetic Liner Inertial Fusion (MagLIF) implosions. The resistive MHD code is used to generate an ensemble of implosions with different liner aspect ratios, initial gas preheat temperatures (that is, different adiabats), and different liner perturbations. The liner density and magnetic field as functions of x , y , and z were generated. The Mallat Scattering Transformation (MST) is taken of the logarithm of both fields and a Principal Components Analysis (PCA) is done on the logarithm of the MST of both fields. The fields are projected onto the PCA vectors and a small number of these PCA vector components are kept. Singular Value Decompositions of the cross correlation of the input parameters to the output logarithm of the MST of the fields, and of the cross correlation of the SVD vector components to the PCA vector components are done. This allows the identification of the PCA vectors vis-a-vis the input parameters. Finally, a Multi Layer Perceptron (MLP) neural network with ReLU activation and a simple three layer encoder/decoder architecture is trained on this dataset to predict the PCA vector components of the fields as a function of time. Details of the implosion, stagnation, and the disassembly are well captured. Examination of the PCA vectors and a permutation importance analysis of the MLP show definitive evidence of an inverse turbulent cascade into a dipole emergent behavior. The orientation of the dipole is set by the initial liner perturbation. The analysis is repeated with a version of the MST which includes phase, called Wavelet Phase Harmonics (WPH). While WPH do not give the physical insight of the MST, they can and are inverted to give field configurations as a function of time, including field-to-field correlations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The infinite level normal forms for non-resonant double Hopf singularities

Here, in this paper, we explore hypernormal forms of vector fields that have non-resonant double Hopf singularities with a non-zero radial cubic part. Our primary focus is on investigating the infinite-level normal form classification of this type of singularities. We provide a normal form decomposition in terms of planar-rotating and planar-radial vector fields, which greatly facilitate the pattern recognition and analysis of the corresponding generalized homological maps. Notably, our paper represents the first instance of the normal form classification for general non-resonant double Hopf singularities without structural symmetry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Chiral rank-$k$ truncations for the multigrid preconditioner of Wilson fermions in lattice QCD

We present a modification to the setup algorithm for the multigrid preconditioner of Wilson fermions in lattice QCD. A larger number of test vectors than that used in conventional multigrid is generated by the smoother. This set of test vectors is then truncated by a singular value decomposition on the chiral components of the test vectors, which are subsequently used to form the prolongation and restriction matrices of the multigrid hierarchy. This modification is demonstrated to improve the convergence of linear equations on an anisotropic lattice with 𝑚𝜋 ≈ 280 MeV from the Hadron Spectrum Collaboration and an isotropic lattice with 𝑚𝜋 ≈ 220 MeV from the MILC Collaboration. The lattice volume dependence of the method is also examined.

Whyte, Travis [Jülich Supercomputing Center, Jülic↗
NASA NTRS1993

Kalman filter estimation of attitude and gyro bias with the QUEST observation model

The loss function for the Wahba attitude estimation problem employs unit vector observations with scalar weights. It is usually associated with the QUEST observation model, where the actual sensor noise is assumed the same for all components of the observed vector, regardless of the position in the field of view. The QUEST model has the great advantage of being sensor-independent except for the scalar parameter that characterizes the sensor errors. Although efficient algorithms for solving the Wahba problem exist, extending these algorithms to estimate gyro biases or sensor alignments has had mixed success. However, it is straightforward to estimate bias and alignment parameters with a Kalman filter. This paper investigates the use of an extended Kalman filter for the attitude and gyro bias that incorporates the QUEST observation model, to be referred to as the Unit Vector Filter (UVF). The UVF results are compared with those from a more conventional filter, the Real-Time Sequential Filter (RTSF), for which the residual is the two-dimensional projection of the unit vector onto the plane perpendicular to the sensor boresight. The RTSF is similar in design to that used by Multimission Modular Spacecraft for onboard attitude determination. An apparent obstacle to the use of unit vectors as measurements is their singular noise covariance matrix. Shuster has shown that this problem should not affect filter performance. The UVF and RTSF are tested using actual flight data from the Extreme Ultraviolet Explorer (EUVE). It is found that these filters generate nearly identical attitude and gyro bias estimates, thus validating the use of unit vectors and the QUEST noise model.

Sedlak, J.↗
DOE OSTI2026

Accelerating multigrid with streaming chiral SVD for Wilson fermions in lattice QCD

A modification to the setup algorithm for the multigrid preconditioner of Wilson fermions in lattice QCD is presented. A larger basis of test vectors than that used in regular multigrid is calculated by the smoother and truncated by singular value decomposition on the chiral components of the test vectors. The truncated basis is used to form the prolongation and restriction matrices of the multigrid hierarchy. This modification of the setup method is demonstrated to increase the convergence of linear solvers on an anisotropic lattice with m π ≈ 239 MeV from the Hadron Spectrum Collaboration and an isotropic lattice with m π ≈ 220 MeV from the MILC Collaboration. The lattice volume dependence of the method is also examined. Increasing the number of test vectors improves speedup up to a point, but storing these vectors becomes impossible in limited memory resources such as GPUs. To address storage cost, we implement a streaming singular value decomposition of the basis of test vectors on the chiral components and demonstrate a decrease in the number of fine level iterations by a factor of 1.7 for m q ≈ m crit

Iterative methods↗
DOE OSTI2025

Spin density wave and van Hove singularity in the kagome metal CeTi 3 Bi 4

Kagome metals with van Hove singularities near the Fermi level can host intriguing quantum phenomena such as chiral loop currents, electronic nematicity, and unconventional superconductivity. However, to our best knowledge, unconventional magnetic states driven by van Hove singularities–like spin-density waves–have not been observed experimentally in kagome metals. Here, we report the magnetic and electronic structure of the layered kagome metal CeTi 3 Bi 4 , where Ti kagome electronic structure interacts with a magnetic sublattice of Ce 3+ J eff = 1/2 moments. Neutron diffraction reveals an incommensurate spin-density wave ground state of the Ce 3+ moments, coexisting with commensurate antiferromagnetic order across most of the temperature-field phase diagram. The commensurate component is preferentially suppressed by thermal fluctuations and magnetic field, yielding a rich phase diagram involving an intermediate single-Q spin-density wave phase. First-principles calculations and angle-resolved photoemission spectroscopy identify van Hove singularities near the Fermi level, with the observed magnetic propagation vectors connecting their high density of states, strongly suggesting a van Hove singularity-assisted spin-density wave. These findings establish kagome metals LnTi 3 Bi 4 as a model platform where the characteristic electronic structure of the kagome lattice plays a pivotal role in magnetic order.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Collisionless tearing in a field-reversed sheet pinch assuming nonparallel propagation

The problem of collisionless linear tearing is examined assuming a wave vector with a component normal to the equilibrium field. The geometry is defined and the general form of the linear dispersion equation is calculated. The linear theory results when k is parallel to B are reviewed, and Ampere's law is calculated for the external adiabatic region when k times B does not equal zero, using two-fluid theory. A solution is obtained for the approximate form of the perturbed currents and vector potential assuming quasi-parallel k. The resonant current contributions within the singular layer are calculated, obtaining an estimate of the dispersion equation. The form of the adiabatic currents within the singular layer is calculated, showing that an x-z current system persists even in the limit k perpendicular to B goes to zero. Finally, the perturbed vector potential solutions across the singular layer are matched to obtain the shape of the complete eigenfunction.

Quest, K. B.↗
NASA NTRS1975

Minimum-fuel rocket trajectories involving intermediate-thrust arcs

The optimal trajectories in the neighborhood of an optimal intermediate-thrust arc are investigated for the minimum-fuel orbit rendezvous problem with fixed specific impulse. Since such an arc is singular, the thrust acceleration magnitude being the singular control component, a second-variation analysis leads to the identification of a field of neighboring, singular arcs in a state space of dimension four rather than six, provided that a suitable Jacobi condition is met. A given neighboring initial six-dimensional state vector does not generally lie on a neighboring singular arc, and junction onto the appropriate singular arc must be accomplished by a short period of strong variations in the acceleration. The neighboring singular arc meets the final condition in 4 dimensions, rather than 6 dimensions, and rendezvous must be completed by another, terminal short period of strong variations in the acceleration. Implications for midcourse guidance near a singular arc are discussed.

Breakwell, J. V.↗
DOE OSTI2023

Classification of gaseous UF 6 assay by femtosecond LIBS in the 424.4 nm spectral region using numerical HOGSVD-DTW features

This technical note presents experimental results using numerical features of fs-LIBS data to classify the assay value of a gaseous UF 6 material. Here, the data-driven feature vectors are computed by Higher Order Generalized Singular Value Decomposition (HOGSVD) and Dynamic Time Warp (DTW). The method achieves 96.97% accuracy in spectral classification testing with fs-LIBS samples obtained from a UF 6 material with five known assay values ranging from 0.287% to 61.740%, with 100% accuracy for the four largest assay values ranging from 4.615% to 61.740%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Characterization and integration of the singular test integrals in the method-of-moments implementation of the electric-field integral equation

In this paper, we characterize the logarithmic singularities arising in the method of moments from the Green’s function in integrals over the test domain, and we use two approaches for designing geometrically symmetric quadrature rules to integrate these singular integrands. These rules exhibit better convergence properties than quadrature rules for polynomials and, in general, lead to better accuracy with a lower number of quadrature points. In this work, we demonstrate their effectiveness for several examples encountered in both the scalar and vector potentials of the electric-field integral equation (singular, near-singular, and far interactions) as compared to the commonly employed polynomial scheme and the double Ma–Rokhlin–Wandzura (DMRW) rules, whose sample points are located asymmetrically within triangles.

42 ENGINEERING↗
NASA NTRS1979

A vector-continuous loading concept for aerodynamic panel methods

An approach to the reduction of discretization errors in aerodynamic panel methods is presented. The approach is based on preventing the occurence of induced velocity singularities at panel slope discontinuities by maintaining continuity of the velocity jump vector across the panels. The approach was implemented in a two-dimensional incompressible panel method formulation and evaluated by application to several external and internal flow problems. The method is shown to exhibit a second order accuracy trend and to produce smaller errors with velocity component boundary conditions imposed on the real flow than with equipotential boundary conditions imposed on the imaginary flow behind the panels. For flows around airfoil sections with either sharp or blunt trailing edges, the method gives excellent agreement with results from a well developed finite difference method. The method is well behaved and is insensitive to irregularities in panel size distribution.

Kemp, W. B., Jr.↗