Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matrix inversion”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Higher Order, Hybrid BEM/FEM Methods Applied to Antenna Modeling

In this presentation, the authors address topics relevant to higher order modeling using hybrid BEM/FEM formulations. The first of these is the limitation on convergence rates imposed by geometric modeling errors in the analysis of scattering by a dielectric sphere. The second topic is the application of an Incomplete LU Threshold (ILUT) preconditioner to solve the linear system resulting from the BEM/FEM formulation. The final tOpic is the application of the higher order BEM/FEM formulation to antenna modeling problems. The authors have previously presented work on the benefits of higher order modeling. To achieve these benefits, special attention is required in the integration of singular and near-singular terms arising in the surface integral equation. Several methods for handling these terms have been presented. It is also well known that achieving ~he high rates of convergence afforded by higher order bases may als'o require the employment of higher order geometry models. A number of publications have described the use of quadratic elements to model curved surfaces. The authors have shown in an EFIE formulation, applied to scattering by a PEC .sphere, that quadratic order elements may be insufficient to prevent the domination of modeling errors. In fact, on a PEC sphere with radius r = 0.58 Lambda(sub 0), a quartic order geometry representation was required to obtain a convergence benefi.t from quadratic bases when compared to the convergence rate achieved with linear bases. Initial trials indicate that, for a dielectric sphere of the same radius, - requirements on the geometry model are not as severe as for the PEC sphere. The authors will present convergence results for higher order bases as a function of the geometry model order in the hybrid BEM/FEM formulation applied to dielectric spheres. It is well known that the system matrix resulting from the hybrid BEM/FEM formulation is ill -conditioned. For many real applications, a good preconditioner is required to obtain usable convergence from an iterative solver. The authors have examined the use of an Incomplete LU Threshold (ILUT) preconditioner . to solver linear systems stemming from higher order BEM/FEM formulations in 2D scattering problems. Although the resulting preconditioner provided aD excellent approximation to the system inverse, its size in terms of non-zero entries represented only a modest improvement when compared with the fill-in associated with a sparse direct solver. Furthermore, the fill-in of the preconditioner could not be substantially reduced without the occurrence of instabilities. In addition to the results for these 2D problems, the authors will present iterative solution data from the application of the ILUT preconditioner to 3D problems.

Fink, P. W.↗
NASA NTRS1982

On the inversion of block tridiagonals without storage constraints

A strategy was developed to permit trade-offs between the number of floating point operations required and the storage requirements for the solution of certain difference problems, such as block tridiagonal systems of equations. This is done by recomputing some intermediate results instead of storing them. Reducing the storage to the square root of the current requirement roughly doubles the number of computations. Reducing the storage more than this tends to make the number of computations prohibitively large. In theory, if m is the order of each sub-matrix in the block tridiagonal matrix, one can solve any linear system with only 5m(2) + 1 temporary storage cells. In many cases m is a constant and quite small. For example, in solving a factored form of the three-dimensional Navier-Stokes equations, the size m of the block tridiagonals is 5. This method lends itself to efficient use on computers with parallel processing or vector processing architectures. On these computers the larger number of floating point operations is more than offset by the decrease in I/O and the increased percentage of vector operations made possible by this algorithm.

Merriam, M. L.↗
DOE OSTI2026

Extracting $B_s\to D_s^*\ellν_\ell$ form factors

Semileptonic $B_{(s)}$ decays are of great phenomenological interest because they allow to determine e.g. CKM matrix elements or test lepton flavor universality. Taking advantage of already existing lattice data, we demonstrate the analysis steps to extract the four form factors describing exclusive semileptonic $B_s\to D_s^*\ellν_\ell$ decays using the narrow width approximation. Our data are based on RBC/UKQCD's set of 2+1 flavor gauge field ensembles with Shamir domain-wall fermion and Iwasaki gauge field action featuring inverse lattice spacings of $a^{-1}=$1.785, 2.383, and 2.785 GeV as well as pion masses between 268 and 433 MeV. Light, strange and charm quarks are simulated using domain-wall fermions, whereas bottom quarks are generated with the relativistic heavy quark (RHQ) action.

Boushmelev, Anastasia [Siegen U.] (ORCID:000000033↗
NASA NTRS1983

The threshold laws for electron-atom and positron-atom impact ionization

The Coulomb-dipole theory is employed to derive a threshold law for the lowest energy needed for the separation of three particles from one another. The study focuses on an electron impinging on a neutral atom, and the dipole is formed between an inner electron and the nucleus. The analytical dependence of the transition matrix element on energy is reduced to lowest order to obtain the threshold law, with the inner electron providing a shield for the nucleus. Experimental results using the LAMPF accelerator to produce a high energy beam of H- ions, which are then exposed to an optical laser beam to detach the negative H- ion, are discussed. The threshold level is found to be confined to the region defined by the upper bound of the inverse square of the Coulomb-dipole region. Difficulties in exact experimental confirmation of the threshold are considered.

Temkin, A.↗
NASA NTRS1978

Inversion of data from diffraction-limited multiwavelength remote sensors. I - Linear case

The remote sensing inverse problem is considered in which the sensor does not necessarily view the same area on the earth at each wavelength, and spatial correlations of the geophysical parameters may be present. Under the conditions of linearity and stationary statistics, the minimum mean-square error solution to the problem of inverting such data is a spatial filter of the Wiener-Kolmogorov class. The resulting remote-sensing system can be characterized by an impulse response matrix in ordinary space or by a transfer matrix in frequency space. A signal-to-noise matrix for the geophysical parameters to be sensed is also defined; this matrix depends on the postulated a priori statistics and on the characteristics of the remote-sensing instrument. The system transfer matrix and the signal-to-noise matrix are simultaneously diagonalizable. The optimum transfer matrix filters out of the estimate vector those eigenparameters for which eigenvalues of the signal-to-noise matrix are less than unity.

Rosenkranz, P. W.↗
DOE OSTI2025

Inversion of dislocation-impurity interactions in α -Fe under magnetic state changes

Here, in this work, we investigate the dislocation-impurity interaction energies and their profiles for various 3d elements —V, Cr, Mn, Cu, Ni, and Co —in and around 1/2 $\langle111\rangle$ screw dislocations in α -Fe using ab initio methods. We consider the ferromagnetic and paramagnetic states, with the latter being modeled through both the disordered local moment model and a spin-wave approach. Our findings reveal that (1) magnetic effects are large compared to size misfit effects of substitutional impurities, and (2) dislocation-impurity interactions are dependent on the magnetic state of the matrix and thermal lattice expansion. In particular, Cu changes from core-attractive in the ferromagnetic state to repulsive in the paramagnetic state.

36 MATERIALS SCIENCE↗
NASA NTRS1976

On a numerical solution of the supersonic panel flutter eigenproblem

An automated digital computer procedure is presented in this paper which enables efficient solution of the eigenvalue problem associated with the supersonic panel flutter phenomena. The step-by-step incremental solution procedure is based on an inverse iteration technique which effectively utilizes solution results from the previous step in determining such results during the current solution step. Also, the computations are limited to the determination of a few specific roots only, which are expected to contain the flutter mode, and this is achieved at each step without having to compute any other root. The structural discretization achieved by the finite-element method yields highly banded stiffness, mass, and aerodynamic matrices; the aerodynamic matrix evaluated by the linearized piston theory is real but unsymmetric in nature. The solution algorithm presented in this paper fully exploits the banded form of the associated matrices, and the resulting computer program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be most efficient and economical when compared to existing procedures of such analysis. Numerical results are presented for a two-dimensional panel flutter problem.

Gupta, K. K.↗
DOE OSTI2025

Imaging of dynamic processes in materials with a laser-wakefield accelerator

Betatron x rays generated from laser-wakefield accelerators are a promising source for imaging dynamic processes in materials. Here, we present proof-of-concept imaging of microstructural evolution in a hypermonotectic Al-Bi alloy, which consists of liquid Bi particles in a solid Al matrix. The images capture the elongation and fragmentation of Bi particles upon isothermal annealing. Because of the femtosecond time scale of the betatron source, the images are not subject to motion blur, whereas the accessibility of the source allows for studies of long-term processes such as annealing. The high-resolution data reveal that the evolution of Bi particles is mediated by an interplay of grain-boundary wetting and morphological instability, in stark contrast to an earlier proposal for (inverse) coarsening.

Alloys↗
NASA NTRS1979

Modeling of terrain gradient for stochastically spaced rows of a measurement matrix

Terrain gradients are employed to evaluate passable regions for unmanned martian roving vehicle. Range data matrix is displaced randomly row wise at the shallow elevation angles near the skyline. The magnitude of the measurement noise in the elevation angles can approach that of the spacing of the same angle. By using a variable incremental data spacing scanning scheme, one can estimate this signal noise ratio. It is found that the error in slope estimate at far distance becomes large for a given elevation angle error. Evaluation of the in-path slopes can be expressed in terms of the inverse of the range slopes. This is because of the fact that the elevation angle is considered as a random variable while the range data are relatively less noisy. An error analysis is performed and it is found that the change of slope is a nonlinear function of the error in elevation angle.

Mediavilla, R.↗
NASA NTRS1987

Spatially random models, estimation theory, and robot arm dynamics

Spatially random models provide an alternative to the more traditional deterministic models used to describe robot arm dynamics. These alternative models can be used to establish a relationship between the methodologies of estimation theory and robot dynamics. A new class of algorithms for many of the fundamental robotics problems of inverse and forward dynamics, inverse kinematics, etc. can be developed that use computations typical in estimation theory. The algorithms make extensive use of the difference equations of Kalman filtering and Bryson-Frazier smoothing to conduct spatial recursions. The spatially random models are very easy to describe and are based on the assumption that all of the inertial (D'Alembert) forces in the system are represented by a spatially distributed white-noise model. The models can also be used to generate numerically the composite multibody system inertia matrix. This is done without resorting to the more common methods of deterministic modeling involving Lagrangian dynamics, Newton-Euler equations, etc. These methods make substantial use of human knowledge in derivation and minipulation of equations of motion for complex mechanical systems.

Rodriguez, G.↗
NASA NTRS2001

Particle Engulfment and Pushing (PEP): Past Micro-Gravity Experiments and Future Experimental Plan on the International Space Station (ISS)

The interaction of an insoluble particle with a growing solid-liquid interface (SLI) has been a subject of investigation for the four decades. For a metallurgist or a material scientist understanding the fundamental physics of such an interaction is relevant for applications that include distribution of reinforcement particles in metal matrix composites, inclusion management in castings, and distribution of Y2Ba1Cu1O5 (211) precipitates (flux pinning sites) in Y1Ba2Cu3O7 (123) superconducting crystals. The same physics is also applicable to other areas including geological applications (frost heaving in soils) and preservation of biological cells. Experimentally this interaction can be quantified in terms of a critical growth velocity, Vcr, of the SLI below which particles are pushed ahead of the advancing interface, and above which the particles are engulfed. Past experimental evidence suggests that this Vcr is an inverse function of the particle radius, R. In order to isolate the fundamental physics that governs such a relationship it is necessary to minimize natural convection at the SLI that is inherent in ground based experiments. Hence for the purpose of producing benchmark data (Vcr vs. R) PEP is a natural candidate for micro-gravity experimentation. Accordingly, experiments with pure Al containing a dispersion of ZrO2 particles and an organic analogue, succinonitrile (SCN) containing polystyrene particles have been performed on the LMS and USMP-4 mission respectively. In this paper we will summarize the experimental data that was obtained during these two micro-gravity missions and show that the results differ compared to terrestrial experiments. We will also discuss the basic elements of our analytical and numerical model and present a comparison of the predictions of these models against micro-gravity experimental data. Finally. we will discuss our future experimental plan that includes the ISS glovebox and MSRRl.

Sen, Subhayu↗
DOE OSTI2024

Thermal quasiparticle theory

The widely used thermal Hartree–Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is postulated in consultation with the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasiparticle) energies in such a way that all fundamental thermodynamic relations are obeyed. The associated density matrix is of a one-electron type, whose diagonal elements take the form of the Fermi–Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi–Dirac or thermal HF theory. The theory thus stipulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green’s function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperatures, the theory approaches finite-temperature MBPT of the same order, but it may outperform the latter at intermediate temperatures by including additional electron-correlation effects through orbital energies. Here, a physical meaning of these thermal orbital energies is proposed (encompassing that of thermal HF orbital energies, which has been elusive) as a finite-temperature version of Janak’s theorem.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1973

A geometric theory for the QR, LU and power iterations.

Consideration of the task of computing the invariant subspaces of a given matrix. For this purpose the LU, QR, treppen and bi-iterations have been presented, used, and studied more or less independently of the old-fashioned power method. Each of these methods generates implicitly a sequence of subspaces which determines the convergence properties of the method. The iterations differ in the way in which a basis is constructed to represent each subspace. This aspect largely determines the usefulness of the method. It is shown that the first four iterations produce exactly the same sequence of subspaces as do direct and inverse iteration started from appropriate subspaces. Their convergence properties are therefore the same, and a complete geometric convergence theory is presented in terms of the power method. It is shown that Hessenberg matrices are associated with ideal starting spaces.

Parlett, B. N.↗
DOE OSTI2024

S -matrix positivity without Lorentz invariance: a case study

We investigate the analytic structure of scattering amplitudes in theories in which Lorentz invariance is spontaneously broken. We do so by computing and studying the S-matrix for a simple example: a superfluid described by a complex scalar with quartic interactions. The computation is confined to tree-level, for there are no absolutely stable single-particle states, though the lifetime can be made long by lowering the chemical potential. For the 2 → 2 amplitude in center-of-mass configurations, not only is crossing symmetry violated, there appears a tree level branch cut for unphysical kinematics. Its appearance is a consequence of non-analyticity in the dispersion relation. The branch point defines a new scale in the problem, which scales inversely with the chemical potential. In this example, even derivatives of the forward amplitude are positive while odd derivatives are negative. This pattern can be understood in a general way in the limit of a small chemical potential, or weak Lorentz breaking.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Identification Uncertainty in Inverse Material Model Parameter Determination: A Sensitivity‐Based Decision Process for Load Path Selection

This research proposes a sensitivity-based framework for selecting the optimal prescribed loading path for a biaxial cruciform specimen. Optimality here is determined by the direction and magnitude of the prescribed displacement that minimizes the influence of random noise on the material model parameter identification. Using simulated experimental data based on finite element simulation, in this work, we identify the material model parameters of a Ludwik hardening model and plane stress implementation of the Hill-48 yield criterion using finite element model updating (FEMU). Our analysis reveals that the identification (or estimator) uncertainty of model parameters depends on the displacement boundary conditions (i.e., loading sequence) and the ground-truth value of the individual parameters. Optimal experimental design (OED) criteria based on the Fisher information matrix were investigated to mitigate indecision in the choice of optimal load path when the identification uncertainty of different material model parameters optimized at different load paths. The determinant of the Fisher information matrix was chosen here as the more useful metric due to its ability to capture uncertainty of the most influential material model parameters. The proposed framework demonstrates potential for real-time automated load step selection using scalar criteria derived prior to mechanical loading. The framework can be generalized to other geometries, boundary conditions and material models, allowing this procedure to be utilized for different experimental configurations and materials.

Fayad, Samuel S. [University of Illinois at Urbana↗
NASA NTRS1999

How to Estimate Attitude from Vector Observations

In many spacecraft attitude systems, the attitude observations are naturally represented as unit vectors. Typical examples are the unit vectors giving the direction to the sun or a star and the unit vector in the direction of the Earth's magnetic field. In 1965, Grace Wahba, proposed the following problem: Find the orthogonal matrix A with determinant +1 that minimizes the loss function L(A) is identity with 1/2(Sum from i a(sub i) (absolute value of b(sub i - A(r(sub i))(exp 2))) where the set of b(sub i) is a set of unit vectors measured in a spacecraft's body frame, the set of r(sub i) are the corresponding unit vectors in a reference frame, and the set of a(sub i) are non-negative weights. Wahba's problem can be related to Maximum Likelihood Estimation if the weights are chosen to be inverse variances, a(sub i) = sigma((sub -2). Wahba didn't assume this, but it will be convenient to assume it in this paper. Wahba'soptimality condition has provided the basis for many attitude determination algorithms. The purpose of this paper is to give an overview of the most popular and most promising algorithm and to provide accuracy and speed comparisons.

Markley,F. Landis↗
NASA NTRS2018

Active and Passive Radiative Transfer Modeling of the Olympic Mountains Experiment

Sensor forward models are an important tool for interpreting remote sensing observations of geophysical phenomena. By implementing a three-dimensional framework, we can simulate and analyze observations from various sensors on disparate platforms. To demonstrate our model framework, we simulate observations from the Olympic Mountains Experiment (OLYMPEX). The use of cloud model simulations allows us to understand sensor response to cloud ice, falling snow, and other processes and features, and the application of model tools to observations allows us to quantify precipitation.MIIST 3D Forward ModelThe Multi-Instrument Inverse Solver Testbed(MIIST) uses the Atmospheric Radiative TransferSimulator (ARTS) for solving the vector radiativetransfer (RT) equation in up to three spatialdimensions within a spherical geometry• Gas absorptiono Line-by-line calculationso Fast transmittance tables• Hydrometeor scattering solverso Discrete ordinateo RT4 (Evans, 1D)o Radar Single Scattering (1D or 3D)o Monte Carlo (3D)Scattering TablesHigh-fidelity hydrometeor scatteringtables are necessary for accurateand consistent forward modeling ofmulti-frequency observations• Requires full Stokes matriceso And absorption vector• Randomly oriented particleso Discrete Dipole Approximationo Characteristic Basis Function Method(coming soon)• Horizontally-oriented plateso Invariant Imbedding T-matrix MethodCloud Resolving SimulationsCloud resolving simulations (e.g.,NU-WRF) supply output consistentwith ARTS needs• Atmospheric Informationo Temperatureo Pressure / heighto Water vapor• Hydrometeor Profileso ARTS architecture ripe for explicit binmicrophysics• Examples use Morrison 2M schemeThe Olympic Mountains Experiment (OLYMPEX)Validation for GPM of mid-latitudefrontal systems approaching nearcoastalmountains from the ocean• Large collection of ground-based andairborne sensorso Radarso Radiometerso In situ• Contemporaneous with RADEXo Two sets of radar at same frequenciesRadiometer Simulation (3 km NUWRF, 20151203, 15:00)2018.12.14 7Simulate 166 GHz polarizationdifference• Corresponds to the presence of aligned icecrystals• Look at trends for both simulations andobservations• Simulations can tolerate lower resolutiono Larger domainSimulations from Observations: OLYMPEXSimulate sensor response usinggeophysical retrievals as input• Single frequency radar retrievals• Multiple scattering enhancementapparent at W band• Spatially dependent phenomenonModeling Application: 1D Retrievals03 December 2015• DC-8 and ER-2 flightso Focus on APR-3 (DC-8)• Citationo Stacked microphysics legso Qualitative comparisonso Range of frozen habitso Presence of supercooledliquid cloudsResults• Retrievals match probeso Good qualitative match• Bands of increasedreflectivity correspond tolarge Dm and highaggregate fraction• Significant amounts ofsupercooled liquid water

Adams, Ian S.↗
DOE OSTI2024

Inducing a tunable skyrmion-antiskyrmion system through ion beam modification of FeGe films

Abstract Skyrmions and antiskyrmions are nanoscale swirling textures of magnetic moments formed by chiral interactions between atomic spins in magnetic noncentrosymmetric materials and multilayer films with broken inversion symmetry. These quasiparticles are of interest for use as information carriers in next-generation, low-energy spintronic applications. To develop skyrmion-based memory and logic, we must understand skyrmion-defect interactions with two main goals—determining how skyrmions navigate intrinsic material defects and determining how to engineer disorder for optimal device operation. Here, we introduce a tunable means of creating a skyrmion-antiskyrmion system by engineering the disorder landscape in FeGe using ion irradiation. Specifically, we irradiate epitaxial B20-phase FeGe films with 2.8 MeV Au 4+ ions at varying fluences, inducing amorphous regions within the crystalline matrix. Using low-temperature electrical transport and magnetization measurements, we observe a strong topological Hall effect with a double-peak feature that serves as a signature of skyrmions and antiskyrmions. These results are a step towards the development of information storage devices that use skyrmions and antiskyrmions as storage bits, and our system may serve as a testbed for theoretically predicted phenomena in skyrmion-antiskyrmion crystals.

74 ATOMIC AND MOLECULAR PHYSICS↗