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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Cormack Research Project: Glasgow University

The aim of this project was to investigate and improve upon existing methods of analysing data from COMITEL on the Gamma Ray Observatory for neutrons emitted during solar flares. In particular, a strategy for placing confidence intervals on neutron energy distributions, due to uncertainties on the response matrix has been developed. We have also been able to demonstrate the superior performance of one of a range of possible statistical regularization strategies. A method of generating likely models of neutron energy distributions has also been developed as a tool to this end. The project involved solving an inverse problem with noise being added to the data in various ways. To achieve this pre-existing C code was used to run Fortran subroutines which performed statistical regularization on the data.

Skinner, Susan↗
NASA NTRS1977

Influence of temperature and the role of chromium on the kinetics of sulfidation of 310 stainless steel

The sulfidation of 310 stainless steel was studied over the temperature range from 910 K to 1285 K. By adjusting the ratio of hydrogen sulfide, variations in sulfur potential were obtained. The effect of temperature on sulfidation was determined at three different sulfur potentials: 39/sqNm, 0.014/sqNm, and 0.00015/sqNm. All sulfide scales contained one or two surface layers in addition to a subscale. The second outer layer (OL-II), furthest from the alloy, contained primarily Fe-Ni-S. The first outer layer (OL-I), nearest the subscale, contained FE-Cr-S. The subscale consisted of sulfide inclusions in the metal matrix. At a given temperature and sulfur potential, the weight gain data obeyed the parabolic rate law after an initial transient period. The parabolic rate constants obtained at the sulfur potential of 39/sqNm did not show a break when the logarithm of the rate constant was plotted as a function of the inverse of absolute temperature. Sulfidation carried out at sulfur potentials below 0.02/sqNm, however, did show a break at 1145 K, which is termed as the transition temperature. This break was found to be associated with the changes which had occurred in the Fe:Cr ratio of OL-I. Below the transition temperature the activation energy was found to be approximately 125 kj/mole. Above the transition temperature the rate of sulfidation decreased with temperature but dependent on the Fe:Cr ratio in the iron-chromium-sulfide layers of the OL-I. A reaction mechanism consistent with the experimental results has been proposed.

Rao, D. B.↗
NASA NTRS1987

On determining important aspects of mathematical models: Application to problems in physics and chemistry

The use of parametric and functional gradient sensitivity analysis techniques is considered for models described by partial differential equations. By interchanging appropriate dependent and independent variables, questions of inverse sensitivity may be addressed to gain insight into the inversion of observational data for parameter and function identification in mathematical models. It may be argued that the presence of a subset of dominantly strong coupled dependent variables will result in the overall system sensitivity behavior collapsing into a simple set of scaling and self similarity relations amongst elements of the entire matrix of sensitivity coefficients. These general tools are generic in nature, but herein their application to problems arising in selected areas of physics and chemistry is presented.

Rabitz, Herschel↗
NASA NTRS2005

Seismic Study of the Solar Interior: Inferences from SOI/MDI Observations During Solar Activity

Work on the structure, asphericity and dynamics of the solar interior from p-mode frequencies and frequency splittings was carried out primarily in collaboration with Dr. Eff-Darwich (University of La Laguna, Tenerife). This ongoing collaboration produced new results for the inversion of the internal solar rotation rate and further development in inversion methodologies. It also resulted in inferences on the solar stratification. Substantial progress towards the characterization of high-degree p-modes has been achieved. In collaboration with Drs. Rabello-Soares and Schou (Stanford University), we have gained a clear conceptual understanding of the various elements that affect the leakage matrix of the SOI/MDI instrument. This work has precise implications on the properties and the characterization of the HMI instrument being developed for the SDO mission.

Korzennik, Sylvain G.↗
DOE OSTI2025

A novel closed-form inversion of the convection–diffusion equation for rapid convection, diffusion, and source profile estimation

To simplify and routinize particle transport analysis in fusion devices, a novel closed form linear inversion of the 1-D convection diffusion equation to estimate diffusion and convection profiles D(r ⃗ ), v(r ⃗ ) and source distribution s(r ⃗ ), of a single species from measured data is derived and demonstrated on synthetic data. Profile estimates of D(r ⃗ ), v(r ⃗ ), s(r ⃗ ) and their uncertainties are given as a matrix expression constructed directly from the incoming density data of the transported species in space and time, as well as physics assumptions such as particle conservation and experimental geometry. The derived matrix expression can be applied to a pumped or non-pumped recycling species, or a non-recycling species that is effectively “pumped” by plasma-facing surfaces.

Hinson, Edward [ORNL] (ORCID:000000019713140X)↗
NASA NTRS1977

Numerical analysis of free vibrations of damped rotating structures

This paper is concerned with the efficient numerical solution of damped and undamped free vibration problems of rotating structures. While structural discretization is achieved by the finite element method, the associated eigenproblem solution is effected by a combined Sturm sequence and inverse iteration technique that enables the computation of a few required roots only without having to compute any other. For structures of complex configurations, a modal synthesis technique is also presented, which is based on appropriate combinations of eigenproblem solution of various structural components. Such numerical procedures are general in nature, which fully exploit matrix sparsity inherent in finite element discretizations, and prove to be most efficient for the vibration analysis of any damped rotating structure, such as rotating machineries, helicopter and turbine blades, spinning space stations, among others.

Gupta, K. K.↗
NASA NTRS1990

Inverse dynamics of adaptive structures used as space cranes

As a precursor to the real-time control of fast moving adaptive structures used as space cranes, a formulation is given for the flexibility induced motion relative to the nominal motion (i.e., the motion that assumes no flexibility) and for obtaining the open loop time varying driving forces. An algorithm is proposed for the computation of the relative motion and driving forces. The governing equations are given in matrix form with explicit functional dependencies. A simulator is developed to implement the algorithm on a digital computer. In the formulations, the distributed mass of the crane is lumped by two schemes, vz., 'trapezoidal' lumping and 'Simpson's rule' lumping. The effects of the mass lumping schemes are shown by simulator runs.

Das, S. K.↗
NASA NTRS1992

Maximum likelihood techniques applied to quasi-elastic light scattering

There is a necessity of having an automatic procedure for reliable estimation of the quality of the measurement of particle size from QELS (Quasi-Elastic Light Scattering). Getting the measurement itself, before any error estimates can be made, is a problem because it is obtained by a very indirect measurement of a signal derived from the motion of particles in the system and requires the solution of an inverse problem. The eigenvalue structure of the transform that generates the signal is such that an arbitrarily small amount of noise can obliterate parts of any practical inversion spectrum. This project uses the Maximum Likelihood Estimation (MLE) as a framework to generate a theory and a functioning set of software to oversee the measurement process and extract the particle size information, while at the same time providing error estimates for those measurements. The theory involved verifying a correct form of the covariance matrix for the noise on the measurement and then estimating particle size parameters using a modified histogram approach.

Edwards, Robert V.↗
DOE OSTI2025

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential for specialized applications in segmented inverse beta decay (IBD) neutrino detectors, astronomy, machine learning, and more. Although this method may easily extend to 3D scalar fields, our focus here is on 2D real-valued fields as it directly applies to directionality. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Data Analysis, Statistics and Probability (physics↗
NASA NTRS1983

Lithification opf gas-rich chondrite regolith breccias by grain boundary and localized shock melting

The fine-grained matrices (less than 150 microns) of 14 gas-rich ordinary chondrile regolith breccias were studied in an attempt to decipher the nature of the lithification process that converted loose regolith material into consolidated breccias. It is found that there is a continuouos gradation in matrix textures from nearly completely clastic (class A) to highly cemented (class C) breccias in which the remining clasts are completely surrounded by interstitial, shock-melted material. It is concluded that this interstitial material is formed by shock melting in the porous regolith. In general, the abundances of solar-wind-implanted He-4 and Ne-20 are inversely correlated with the abundance of intenstitial, shock-melted, feldspathic material. Chondrites with the highest abundance of interstitial, melted material (class C) experienced the highest shock pressures and temperatures and suffered the most extensive degassing. It is this interstitial, feldspathic melt that lithifies and cements the breccias together; those breccias with very little interstitial melt (class A) are the most porous and least consolidated.

Bischoff, A.↗
NASA NTRS1999

Efficient Implementations of the Quadrature-Free Discontinuous Galerkin Method

The efficiency of the quadrature-free form of the dis- continuous Galerkin method in two dimensions, and briefly in three dimensions, is examined. Most of the work for constant-coefficient, linear problems involves the volume and edge integrations, and the transformation of information from the volume to the edges. These operations can be viewed as matrix-vector multiplications. Many of the matrices are sparse as a result of symmetry, and blocking and specialized multiplication routines are used to account for the sparsity. By optimizing these operations, a 35% reduction in total CPU time is achieved. For nonlinear problems, the calculation of the flux becomes dominant because of the cost associated with polynomial products and inversion. This component of the work can be reduced by up to 75% when the products are approximated by truncating terms. Because the cost is high for nonlinear problems on general elements, it is suggested that simplified physics and the most efficient element types be used over most of the domain.

Lockard, David P.↗
NASA NTRS2019

Bayesian Inference and the Effects of Varying Uncertainty Models in Charring Ablator Calibration and Uncertainty Quantification Problems

The Mars Science Laboratory (MSL) vehicle utilized a heat shield constructed from NASA’s Phenolic-Impregnated Carbon Ablator (PICA) material to protect the main structure from the high enthalpy environment encountered during hypersonic atmospheric entry. During the vehicle’s descent through Martian atmosphere, multiple thermocouples embedded within the heat shield captured in-depth material temperature data that allow for studies to be conducted on current material response reconstruction tools. In the present work, material temperature data obtained from thermocouples within the MISP-4 plug (MEDLI (Mars Science Laboratory Entry, Descent, and Landing Instrument) Integrated Sensor Plug) are utilized in the calibration of Theoretical Ablative Composite for Open Testing (TACOT) model parameters in conjunction with NASA’s Porous material Analysis Toolbox (PATO) through Bayesian inference where uncertainty due to parametric, modeling, and experimental sources is simultaneously quantified. Prior to the study, a sensitivity analysis is performed through computation of the robust Sobol indices in an effort to study the relationship between input space and model response and to reduce the dimensionality of the statistical inverse problem. The Bayesian inference methodology necessitates an a-priori choice to be made for the uncertainty model for which numerous possibilities are available. Across most works, however, only basic additive or multiplicative models are utilized with pre-defined magnitudes of uncertainty based on a-priori knowledge or to-be-calibrated multipliers of static covariance matrix structures. The present effort explores the effects of informed uncertainty models, ones with temporal dependence that are simultaneously calibrated through Bayesian inference, on calibrated results for parameters that make up the uncertain input space.

Sensitivity Analysis↗
DOE OSTI2024

Comparison of machine learning and electrical resistivity arrays to inverse modeling for locating and characterizing subsurface targets

Here, this study evaluates the performance of multiple machine learning (ML) algorithms and electrical resistivity (ER) arrays for inversion with comparison to a conventional Gauss-Newton numerical inversion method. Four different ML models and four arrays were used for the estimation of only six variables for locating and characterizing hypothetical subsurface targets. The combination of dipole-dipole with Multilayer Perceptron Neural Network (MLP-NN) had the highest accuracy. Evaluation showed that both MLP-NN and Gauss-Newton methods performed well for estimating the matrix resistivity while target resistivity accuracy was lower, and MLP-NN produced sharper contrast at target boundaries for the field and hypothetical data. Both methods exhibited comparable target characterization performance, whereas MLP-NN had increased accuracy compared to Gauss-Newton in prediction of target width and height, which was attributed to numerical smoothing present in the Gauss-Newton approach. MLP-NN was also applied to a field dataset acquired at U.S. DOE Hanford site.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1973

Mixed finite-difference scheme for free vibration analysis of noncircular cylinders

A mixed finite-difference scheme is presented for the free-vibration analysis of simply supported closed noncircular cylindrical shells. The problem is formulated in terms of eight first-order differential equations in the circumferential coordinate which possess a symmetric coefficient matrix and are free of the derivatives of the elastic and geometric characteristics of the shell. In the finite-difference discretization, two interlacing grids are used for the different fundamental unknowns in such a way as to avoid averaging in the difference-quotient expressions used for the first derivative. The resulting finite-difference equations are symmetric. The inverse-power method is used for obtaining the eigenvalues and eigenvectors.

Noor, A. K.↗
DOE OSTI2026

Dual Representations and H ∞ -Optimal Control of Partial Differential Equations

We consider H ∞ -optimal state-feedback control of the class of linear Partial Differential Equations (PDEs) which admit a Partial Integral Equation (PIE) representation. While linear matrix inequalities are commonly used for optimal control of Ordinary Differential Equations (ODEs), the absence of a universal state-space representation and suitable dual form prevents such methods from being applied to optimal control of PDEs. Specifically, for ODEs, the controller synthesis problem is defined in state-space, and duality is used to resolve the bilinearity of that synthesis problem. Recently, the PIE representation was proposed as a universal state-space representation for linear PDE systems. In this paper, we show that any PDE system represented by a PIE admits a dual PIE with identical stability and I/O properties. This result allows us to reformulate the stabilizing and optimal state-feedback control problems as convex optimization over the cone of positive Partial Integral (PI) operators. Operator inversion formulae then allow us to construct feedback gains for the original PDE system. The results are verified through application to several canonical problems in optimal control of PDEs and indicate the resulting bounds on H ∞ norm are not conservative.

42 ENGINEERING↗
DOE OSTI2024

Aligned yet large dipoles: a SMEFT study

We study a non-universal flavor scenario at the level of the Standard Model Effective Field Theory, according to which the matrix of Wilson coefficients cuW of an up-type electroweak quark dipole operator is aligned with the up-type Yukawa coupling. Such an alignment usually follows from the assumption of Minimal Flavor Violation (MFV), away from which we step by allowing the entries of cuW to be sizable along the first quark generations. A particular example, which we refer to as “inverse hierarchy MFV”, features Wilson coefficients inversely proportional to quark masses, and arises from BSM models respecting MFV and containing heavy fields that replicate the mass hierarchy of SM quarks. We then analyze the phenomenology driven by cuW at colliders and at lower-energy flavor experiments. We show that precision measurements of the process pp → Wh → γγℓν at FCC-hh could set an upper bound on |cuW| ≲ O (10 -2 )(Λ/TeV) 2 , with Λ the cutoff of the effective field theory. This bound is an order of magnitude stronger than the existing LHC bounds. Moreover, we estimate that Wh → b$\overline{b}$ℓν at HL-LHC could also give competitive bounds. In the low-energy regime, we consider bounds arising from rare kaon decays, which turn out to be loose, |c$_{uW}^{11}$ | < O (1)(Λ/TeV) 2 . We finally demonstrate that our flavor and operator assumptions can be derived from a weakly-coupled UV model, which we choose to simultaneously illustrate the UV origin of inverse hierarchy MFV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1980

Errors incurred in a plane-wave-type expansion of a Gaussian beam

The multipole expansion obtained by Morita et al. (1968) of the Gaussian laser beam used to levitate an aerosol particle in order that its complete phase matrix may be measured is compared with that of Tsai and Pogorzelski (1975) in order to demonstrate the effect of the incorrect expansion used by Morita. Errors incurred by the use of an equation in which one side satisfies the scalar wave equation while the other side does not and can be reduced to a plane wave amplitude are calculated as functions of the inverse of the wave number times the beam waist, the wave number times the radial spherical coordinate and the angular spherical coordinate. Errors on the order of a few percent, considered undetectable are obtained in the squared-field amplitudes due to the expansion, however, they are found to become significant (several tens of percent) when the angle is zero. It is concluded that the expansion of Morita should only be used in the regions where the spherical angle is less than 0.01 and its product with the wave number and the radial spherical coordinate is less than unity.

Kattawar, G. W.↗
NASA NTRS1987

Steering Law Design for Redundant Single Gimbal Control Moment Gyro Systems

The correspondence between robotic manipulators and single gimbal Control Moment Gyro (CMG) systems was exploited to aid in the understanding and design of single gimbal CMG Steering laws. A test for null motion near a singular CMG configuration was derived which is able to distinguish between escapable and unescapable singular states. Detailed analysis of the Jacobian matrix null-space was performed and results were used to develop and test a variety of single gimbal CMG steering laws. Computer simulations showed that all existing singularity avoidance methods are unable to avoid Elliptic internal singularities. A new null motion algorithm using the Moore-Penrose pseudoinverse, however, was shown by simulation to avoid Elliptic type singularities under certain conditions. The SR-inverse, with appropriate null motion was proposed as a general approach to singularity avoidance, because of its ability to avoid singularities through limited introduction of torque error. Simulation results confirmed the superior performance of this method compared to the other available and proposed pseudoinverse-based Steering laws.

Bedrossian, Nazareth Sarkis↗