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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

The general linear inverse problem - Implication of surface waves and free oscillations for earth structure.

The discrete general linear inverse problem reduces to a set of m equations in n unknowns. There is generally no unique solution, but we can find k linear combinations of parameters for which restraints are determined. The parameter combinations are given by the eigenvectors of the coefficient matrix. The number k is determined by the ratio of the standard deviations of the observations to the allowable standard deviations in the resulting solution. Various linear combinations of the eigenvectors can be used to determine parameter resolution and information distribution among the observations. Thus we can determine where information comes from among the observations and exactly how it constraints the set of possible models. The application of such analyses to surface-wave and free-oscillation observations indicates that (1) phase, group, and amplitude observations for any particular mode provide basically the same type of information about the model; (2) observations of overtones can enhance the resolution considerably; and (3) the degree of resolution has generally been overestimated for many model determinations made from surface waves.

Wiggins, R. A.↗
NASA NTRS1995

Optimization methods, flux conserving methods for steady state Navier-Stokes equation

Navier-Stokes equation as discretized by new flux conserving method proposed by Chang and Scott results in the system: vector F(vector x) = 0, where F is a vector valued function. The Optimization method we use is based on Quasi-Newton methods: given a nonlinear function vector F(vector x) = 0, we solve, Delta(vector x) = -BF(vector x), where Delta(vector x) is the correction term and B is the inverse Jacobian of F(x). Then, iteratively, vector(x(sub (i+1))) = vector(x (sub i)) + alpha.Delta(vector x(sub i)), where alpha is a line search correction term determined by a line search routine. We use the BFCG's update the Jacobian matrix B(sub k) at each iteration. It is well known that B(sub k) approaches B(*) at the solution X(*). This algorithm has several advantages over the Newton-Raphson method. For example, we do not need to calculate the Jacobian matrix at each iteration which is computationally very expensive.

Adeyeye, John↗
NASA NTRS2010

GLGM-3: A Degree-ISO Lunar Gravity Model from the Historical Tracking Data of NASA Moon Orbiters

In preparation for the radio science experiment of the Lunar Reconnaissance Orbiter (LRO) mission, we analyzed the available radio tracking data of previous NASA lunar orbiters. Our goal was to use these historical observations in combination with the new low-altitude data to be obtained by LRO. We performed Precision Orbit Determination on trajectory arcs from Lunar Orbiter 1 in 1966 to Lunar Prospector in 1998, using the GEODYN II program developed at NASA Goddard Space Flight Center. We then created a set of normal equations and solved for the coefficients of a spherical harmonics expansion of the lunar gravity potential up to degree and order 150. The GLGM-3 solution obtained with a global Kaula constraint (2.5 x 10(exp -4)/sq l) shows good agreement with model LP150Q from the Jet Propulsion Laboratory, especially over the nearside. The levels of data fit with both gravity models are very similar (Doppler RMS of approx.0.2 and approx. 1-2 mm/s in the nominal and extended phases, respectiVely). Orbit overlaps and uncertainties estimated from the covariance matrix also agree well. GLGM-3 shows better correlation with lunar topography and admittance over the nearside at high degrees of expansion (l > 100), particularly near the poles. We also present three companion solutions, obtained with the same data set but using alternate inversion strategies that modify the power law constraint and expectation of the individual spherical harmonics coefficients. We give a detailed discussion of the performance of this family of gravity field solutions in terms of observation fit, orbit quality, and geophysical consistency.

Mazarico, E.↗
NASA NTRS1974

Development of technology for fluid-structure interaction modeling of a 1/8-scale dynamic model of shuttle External Tank (ET). Volume i: Technical report

A NASTRAN analysis of the external tank (ET) substructure of the 1/8-scale space shuttle structural dynamics model is described. The NASTRAN hydroelastic procedures were used to form a model of the liquid oxygen portion of the ET. Large computer storage requirements and running times were required unless these procedures were modified. Several possibilities were demonstrated including the substitution of the real for the complex eigenvalue routine and the use of the OMIT capability to reduce the number of fluid coordinates. A NASTRAN model of the complete ET was formed and reduced to 252 degrees of freedom. A review of the eigenvector extracted, using the unsymmetrical inverse power method, indicated that the structural OMIT's resulted in unsatisfactory modal deflections. Modifications to the NASTRAN hydroelastic capability appeared necessary to generate a successful ET model. An approach is described which, by assuming incompressibility, reduces the fluid representation to a symmetric mass matrix which can be added to the structural mass. The problem is solved using faster and more efficient eigenvalue routines. Using this approach the ET NASTRAN model was analyzed for three separate weight conditions.

Bernstein, M.↗
NASA NTRS1991

Modeling the spatiotemporal organization of velocity storage in the vestibuloocular reflex by optokinetic studies

1. A generalized three-dimensional state space model of visual vestibular interaction was developed. Matrix and dynamical system operators associated with inputs from the semicircular canals, otolith velocity estimator, and the visual system have been incorporated into the model, which focus on their relationship to the velocity storage integrator. 2. A relationship was postulated between the eigenvalues and the direction of the eigenvectors of the system matrix and the orientation of the spatial vertical. It was assumed that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix, whereas another rotates the eigenvectors of the system matrix. The pitch and roll eigenvectors rotate with the head, whereas the yaw axis eigenvector remains approximately spatially invariant. 3. Based on the three-dimensional model, a computational procedure was formulated to identify the eigenvalues and eigenvectors of the system matrix with the use of a modification of the marquardt algorithm. With the use of data obtained from a monkey, it was shown that the three-dimensional behavior of velocity storage cannot be predicted solely in terms of its time constants, i.e., the inverse of its eigenvalues. With the use of the same eigenvalues the data could either be fit or not fit, depending on the eigenvector directions. Therefore, it is necessary to specify eigenvector directions when characterizing velocity storage in three dimensions. 4. Parameters found with the use of the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced for coupling the visual system to the integrator and to the direct optokinetic pathway. Simulations of optokinetic nystagmus (OKN) and optokinetic after-nystagmus (OKAN) were run. The model predicted the behavior of yaw and pitch OKN and OKAN when the animal is upright. It also predicted the cross-coupling in the side down position. The trajectories in velocity space were also accurately simulated. 5. One of the predictions of the model is that when the stimulus direction is along an eigenvector, the trajectory in velocity space is a straight line. Using the "spectral width" of the residuals from a straight line sequence during OKAN, we developed a methodology to estimate how close the OKAN decay was to an eigenvector trajectory. 6. Thus we have developed a model-based approach for studying and interpreting the response characteristics of velocity storage in three dimensions.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗
DOE OSTI2025

LuGo: An enhanced quantum phase estimation implementation

Quantum Phase Estimation (QPE) is a cardinal algorithm in quantum computing that plays a crucial role in various applications, including cryptography, molecular simulation, and solving systems of linear equations. However, the standard implementation of QPE faces challenges related to time complexity and circuit depth, which limit its practicality for large-scale computations. We introduce LuGo, a novel framework designed to enhance the performance of QPE by reducing circuit duplication, as well as using parallelization techniques to achieve faster generation of the QPE circuit and gate reduction. We validate the effectiveness of our framework by generating quantum linear solver circuits, which require both QPE and inverse QPE, to solve linear systems of equations. LuGo achieves significant improvements in both computational efficiency and hardware requirements without compromising on accuracy. Compared to a standard QPE implementation, LuGo reduces time consumption to generate a circuit that solves a 2 6 × 2 6 system matrix by a factor of 50.68 and over 31× reduction of quantum gates and circuit depth, with no fidelity loss on an ideal quantum simulator. Furthermore, we demonstrated the versatility and scalability of LuGo enabled HHL algorithm by simulating a canonical Hele-Shaw fluid problem using a quantum simulator. With these advantages, LuGo paves the way for more efficient implementations of QPE, enabling broader applications across several quantum computing domains.

Quantum algorithm↗
NASA NTRS1975

The effect of microstructure and strength on the fracture toughness of an 18 Ni, 300 grade maraging steel

A 300 grade maraging steel was chosen as a vehicle by which to understand the inverse relationship between strength and toughness in high strength alloys such as the 18 Ni maraging steels. The 18 Ni, 300 grade maraging material was a commercial grade consumable-electrode, vacuum arc remelted heat obtained in the form of forged and annealed plate. The matrix contained a population of second-phase impurity inclusions which was a product of the casting and hot working processes. These inclusions did not change with subsequent precipitation hardening. Changes in microstructure resulting in strength increases were brought about by variations in aging temperature and time. Maximum strength was attained in the 300 grade maraging steel by aging at 427 C (800 F) for 100 hours. Tensile, fatigue precracked Charpy impact, and plane-strain fracture toughness tests were performed at room temperature, 20 C (68 F). With increasing strength the fracture toughness decreases as smaller and smaller inclusions act as sites for void initiation.

Psioda, J. A.↗
NASA NTRS1995

A Hybrid Parallel Preconditioning Algorithm For CFD

A new hybrid preconditioning algorithm will be presented which combines the favorable attributes of incomplete lower-upper (ILU) factorization with the favorable attributes of the approximate inverse method recently advocated by numerous researchers. The quality of the preconditioner is adjustable and can be increased at the cost of additional computation while at the same time the storage required is roughly constant and approximately equal to the storage required for the original matrix. In addition, the preconditioning algorithm suggests an efficient and natural parallel implementation with reduced communication. Sample calculations will be presented for the numerical solution of multi-dimensional advection-diffusion equations. The matrix solver has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The full paper will show numerous examples in CFD to demonstrate the efficiency and robustness of the method.

Barth,Timothy J.↗
DOE OSTI2025

Laser powder bed fusion of ODS Fe–Cr–Al (0.3Zr, 0.3Y 2 O 3 ): Unveiling processing-microstructure- mechanical property relationships

Here, this study investigates the fabrication of oxide dispersion strengthened (ODS) Fe-Cr-Al alloys via laser powder bed fusion (LPBF) with strategic additions of 0.3 wt% Zr and 0.3 wt% Y 2 O 3 for enhanced mechanical performance in nuclear applications. Systematic processing parameter optimization yielded three distinct conditions: one low-density product with significant defects and two near-full-density materials with improved consolidation. Comprehensive characterization confirmed single-phase α-ferrite matrix formation with successful incorporation of Y-, Zr-, O-, and C-rich precipitates characteristic of ODS alloys. However, precipitate density remained low (∼10 7 cm −3 ), resulting in sink strength values substantially below optimal levels for radiation resistance. Microhardness values (mid-200s HV) correlated inversely with grain size following the Hall-Petch relationship, indicating grain boundary strengthening as the dominant mechanism rather than precipitation strengthening. The optimized processing conditions achieved excellent mechanical properties with room temperature yield strength of approximately 500 MPa and 30 % elongation, demonstrating superior strength-ductility synergy compared to other additively manufactured ODS materials and performance consistent with literature values for LPBF-processed ODS-FeCrAl alloys. This investigation reveals both the potential and limitations of LPBF processing for ODS Fe-Cr-Al alloys. While successful defect-free fabrication was achieved, results highlight the critical need for systematic optimization of processing parameters and post-processing heat treatments to enhance precipitate density for effective dispersion strengthening and radiation resistance while maintaining additive manufacturing advantages.

Additive manufacturing↗
NASA NTRS1991

Development of advanced control schemes for telerobot manipulators

To study space applications of telerobotics, Goddard Space Flight Center (NASA) has recently built a testbed composed mainly of a pair of redundant slave arms having seven degrees of freedom and a master hand controller system. The mathematical developments required for the computerized simulation study and motion control of the slave arms are presented. The slave arm forward kinematic transformation is presented which is derived using the D-H notation and is then reduced to its most simplified form suitable for real-time control applications. The vector cross product method is then applied to obtain the slave arm Jacobian matrix. Using the developed forward kinematic transformation and quaternions representation of the slave arm end-effector orientation, computer simulation is conducted to evaluate the efficiency of the Jacobian in converting joint velocities into Cartesian velocities and to investigate the accuracy of the Jacobian pseudo-inverse for various sampling times. In addition, the equivalence between Cartesian velocities and quaternion is also verified using computer simulation. The motion control of the slave arm is examined. Three control schemes, the joint-space adaptive control scheme, the Cartesian adaptive control scheme, and the hybrid position/force control scheme are proposed for controlling the motion of the slave arm end-effector. Development of the Cartesian adaptive control scheme is presented and some preliminary results of the remaining control schemes are presented and discussed.

Nguyen, Charles C.↗
NASA NTRS1992

Error analysis applied to several inversion techniques used for the retrieval of middle atmospheric constituents from limb-scanning MM-wave spectroscopic measurements

The formal retrieval error analysis of Rodgers (1990) allows the quantitative determination of such retrieval properties as measurement error sensitivity, resolution, and inversion bias. This technique was applied to five numerical inversion techniques and two nonlinear iterative techniques used for the retrieval of middle atmospheric constituent concentrations from limb-scanning millimeter-wave spectroscopic measurements. It is found that the iterative methods have better vertical resolution, but are slightly more sensitive to measurement error than constrained matrix methods. The iterative methods converge to the exact solution, whereas two of the matrix methods under consideration have an explicit constraint, the sensitivity of the solution to the a priori profile. Tradeoffs of these retrieval characteristics are presented.

Puliafito, E.↗
NASA NTRS1984

Composition and origin of enstatite in E chondrites

A petrologic and electron-probe study of enstatite in six EH3, 4, and 5 chondrites, including the previously undescribed Allan Hills A77156, provides no support for Leitch and Smith's (1982) model in which crystals and liquids from two chemically distinct, partially molten planetesimals are mixed. It is found that Cr and Mn concentrations in enstatite grains in matrix and chondrules are not uniquely defined by their cathodoluminescence color. Wide variations in enstatite compositions among chondrules may partly reflect heterogeneities in the precursor material. However, correlated variations of Cr, Mn, and Fe concentrations in enstatite are probably largely due to fractionation during chondrule crystallization. The inverse correlation of petrologic type with minor element concentrations in enstatite appears to be a metamorphic effect.

Mckinley, S. G.↗
DOE OSTI2024

Advancing attenuation estimation through integration of the Hessian in multiparameter viscoacoustic full-waveform inversion

Accurate seismic attenuation models of subsurface structures not only enhance subsequent migration processes by improving fidelity, resolution, and facilitating amplitude-compliant angle gather generation but also provide valuable constraints on subsurface physical properties. Leveraging full-wavefield information, multiparameter viscoacoustic full-waveform inversion ( Q-FWI) simultaneously estimates seismic velocity and attenuation ( Q) models. However, a major challenge in Q-FWI is the contamination of crosstalk artifacts, where inaccuracies in the velocity model are mistakenly mapped to the inverted attenuation model. While incorporating the Hessian is expected to mitigate these artifacts, the explicit implementation is prohibitively expensive due to its formidable computational cost. In this study, we formulate and develop a Q-FWI algorithm via the Newton-conjugate gradient (CG) framework, where the search direction at each iteration is determined through an internal CG loop. In particular, the Hessian is integrated into each CG step in a matrix-free fashion using the second-order adjoint-state method. We find through synthetic experiments that our Newton-CG Q-FWI significantly mitigates crosstalk artifacts compared with the limited-memory Broyden-Fletcher-Goldfarb-Shanno method and the CG method, albeit with a notable computational cost. In the discussion of several key implementation details, we also determine the significance of the approximate Gauss-Newton Hessian, the second-order adjoint-state method, and the two-stage inversion strategy.

Geochemistry & Geophysics↗
NASA NTRS1979

Off-axis impact of unidirectional composites with cracks: Dynamic stress intensification

The dynamic response of unidirectional composites under off axis (angle loading) impact is analyzed by assuming that the composite contains an initial flaw in the matrix material. The analytical method utilizes Fourier transform for the space variable and Laplace transform for the time variable. The off axis impact is separated into two parts, one being symmetric and the other skew-symmetric with reference to the crack plane. Transient boundary conditions of normal and shear tractions are applied to a crack embedded in the matrix of the unidirectional composite. The two boundary conditions are solved independently and the results superimposed. Mathematically, these conditions reduce the problem to a system of dual integral equations which are solved in the Laplace transform plane for the transformation of the dynamic stress intensity factor. The time inversion is carried out numerically for various combinations of the material properties of the composite and the results are displayed graphically.

Sih, G. C.↗
NASA NTRS1986

Microwave remote sensing: Active and passive. Volume 3 - From theory to applications

Aspects of volume scattering and emission theory are discussed, taking into account a weakly scattering medium, the Born approximation, first-order renormalization, the radiative transfer method, and the matrix-doubling method. Other topics explored are related to scatterometers and probing systems, the passive microwave sensing of the atmosphere, the passive microwave sensing of the ocean, the passive microwave sensing of land, the active microwave sensing of land, and radar remote sensing applications. Attention is given to inversion techniques, atmospheric attenuation and emission, a temperature profile retrieval from ground-based observations, mapping rainfall rates, the apparent temperature of the sea, the emission behavior of bare soil surfaces, the emission behavior of vegetation canopies, the emission behavior of snow, wind-vector radar scatterometry, radar measurements of sea ice, and the back-scattering behavior of cultural vegetation canopies.

Ulaby, F. T.↗
NASA NTRS1988

Random field estimation approach to multibody dynamics

The use of recursive techniques similar to random field models to factor and invert the inertia matrix of a multibody system is discussed. An equivalence is established between the composite multibody system inertia matrix and the covariance of the output of a described linear system model. Conditional mean estimation and sequential estimation problems are solved along with problems of filtering and smoothing. Formulas are developed to compute the covariance of several relevant quantities. The foregoing results are used to obtain the inverse of the composite multibody system inertia in closed form.

Rodriguez, G.↗
DOE OSTI2024

Brittle failure analysis and modeling of high-burnup PWR fuel cladding alloys

The aim of this research is the development of methods for predicting mechanical behavior and identification of limiting conditions to prevent brittle failure of high-burnup (HBU) pressure water reactor (PWR) fuel cladding alloys. A finite element (FE) model of the ring compression test (RCT) was created to analyze the failure behavior of zirconium-based alloys with radial hydrides during the RCT. An elastic-plastic material model describes the zirconium alloy. The stress-strain curve needed for the elastic-plastic material model was derived by inverse finite element analyses. Cohesive zone modeling is used to reproduce sudden load drops during RCT loading. Based on the failure mechanism in non-irradiated ZIRLO (R) claddings, a micro-mechanical model was developed that distinguishes between brittle failure along hydrides and ductile failure of the zirconium matrix. Two different cohesive laws representing these types of failure are present in the same cohesive interface. The key differences between these constitutive laws are the cohesive strength, the stress at which damage initiates, and the cohesive energy, which is the damage energy dissipated by the cohesive zone. Statistically generated matrix-hydride distributions were mapped onto the cohesive elements and simulations with focus on the first load drop were performed. Computational results are in good agreement with the RCT results conducted on high-burnup M5 (R) samples. It could be shown that crack initiation and propagation strongly depend on the specific configuration of hydrides and matrix material in the fracture area.

Simbruner, Kai↗
NASA NTRS1982

Measurement of the rugged invariants of magnetohydrodynamic turbulence in the solar wind

Measurements of the total energy, cross helicity, and magnetic helicity of the solar wind at 1, 2.8, and 5 AU are presented. These quantities are the three rugged invariants of three-dimensional ideal incompressible MHD turbulence theory. The theoretical technique for measuring the magnetic helicity from the matrix of two-point correlations is shown. The length scales characterizing the magnetic helicity are found to be equal to or greater than those which characterize the magnetic energy. The magnetic helicity typically lies at scales larger than the magnetic correlation length, consistent with the expectations of the inverse cascade and selective decay hypotheses of three-dimensional MHD turbulence. At smaller scales, the magnetic helicity oscillates in sign. Measurements of the cross helicity are not fully consistent with the usual interpretation in terms of outward propagating Alfvenic functions. Especially during the interval at 5 AU the cross helicity is found to oscillate in sign indicating fluctuations propagating both outward and inward.

Matthaeus, W. H.↗