Efficiency of generalized matrix inversion methods
Efficiency of computer programs to calculate generalized inverse of arbitrary real matrix
Delaney, F. C.↗
SEARCH · Search NASA
Results for “inverse calculation”
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.
Quantum mechanical calculations of vibrational population inversion in chemical reactions - Numerically exact L-squared-amplitude-density study of the H2Br reactive system
Numerically exact, fully three-dimensional quantum mechanicl reactive scattering calculations are reported for the H2Br system. Both the exchange (H + H-prime Br to H-prime + HBr) and abstraction (H + HBR to H2 + Br) reaction channels are included in the calculations. The present results are the first completely converged three-dimensional quantum calculations for a system involving a highly exoergic reaction channel (the abstraction process). It is found that the production of vibrationally hot H2 in the abstraction reaction, and hence the extent of population inversion in the products, is a sensitive function of initial HBr rotational state and collision energy.
Zhang, Y. C.↗
HYMPS: Numerical Techniques
The digital calculations that drive the classification portion of the Hybrid Pattern Recognition System (HYMPS) are considered. A revision of these calculations is suggested that involves three items: matrix inversion, det calculations and singularity, and covariance factorization. It is shown that it is more economical to first factor the covariance matrix, and by so doing, delete the matrix inversion routine. In addition, the necessary det calculations can be more easily realized by use of simple theoretical facts about the factorization.
Decell, H. P., Jr.↗
Accelerating and stabilizing the convergence of vibrational self-consistent field calculations via the direct inversion of the iterative subspace (vDIIS) algorithm
The vibrational self-consistent field (VSCF) method yields anharmonic states and spectra for molecular vibrations, and it serves as the starting point for more sophisticated correlated-vibration methods. Convergence of the iterative, non-linear optimization in VSCF calculations can be erratic or altogether unsuccessful, particularly for chemical systems involving low-frequency motions. In this work, a vibrational formulation of the Direct Inversion of the Iterative Subspace method of Pulay is presented and investigated. This formulation accounts for distinct attributes of the vibrational and electronic cases, including the expansion of each single-mode vibrational wavefunction in its own basis set. The resulting Direct Inversion of the Iterative Subspace method is shown to substantially accelerate VSCF convergence in all convergent cases as well as rectify many cases where Roothaan-based methods fail. Performance across systems ranging from small, rigid molecules to weakly bound molecular clusters is investigated in this analysis.
Chemistry↗
Reconciling calculations and measurements of inverse bremsstrahlung absorption
It was recently shown that the use of Coulomb logarithms appropriate for bremsstrahlung radiation (rather than transport processes) along with corrections for the Langdon effect and ion screening reproduced measurements of collisional absorption in well-characterized underdense plasmas [D. Turnbull et al., Phys. Rev. Lett. 130, 145103 (2023)]. However, it was recognized at the time that the use of the standard absorption-reduction factor from Langdon's seminal paper was inconsistent with the use of Coulomb logarithms that are thermally averaged over a Maxwellian electron-velocity distribution function. A more accurate approach would be to average over the expected super-Gaussian distribution function while accounting for the Gaunt factor's velocity dependence, which somewhat mitigates the Langdon effect; however, at that time, this theory matched the data less well. This conflict is now eliminated with the additional insight that the ionization state of our mid-Z ion species (when present) was lower than had been assumed, as evidenced by the Thomson-scattering data and time-dependent Cretin simulations. We are now able to show that an improved treatment of the Langdon effect provides the best match to data. Otherwise, the prior conclusions remain unchanged. We also show an example of the substantial expected impact to the absorption rate in calculations of indirect-drive hohlraums.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Fast Computation of the Inverse CMH Model
A fast computational method based on differential equation approach for inverse DOK model has been extended for the inverse CMH model. Also, a cobweb technique for calculating the inverse CMH model is also presented. The two techniques are differed from the point of view of flexibility and computation time.
Patel, Umesh D.↗
Nonlinear Inverse Scattering Applied to Calculation of the Effective Permittivity of Random Collections of Scatterers
Explore the source record for details and available documents.
two-dimensional scatterers cylinder Monte Carlo Si↗
Reconfigurable Control with Neural Network Augmentation for a Modified F-15 Aircraft
Description of the performance of a simplified dynamic inversion controller with neural network augmentation follows. Simulation studies focus on the results with and without neural network adaptation through the use of an F-15 aircraft simulator that has been modified to include canards. Simulated control law performance with a surface failure, in addition to an aerodynamic failure, is presented. The aircraft, with adaptation, attempts to minimize the inertial cross-coupling effect of the failure (a control derivative anomaly associated with a jammed control surface). The dynamic inversion controller calculates necessary surface commands to achieve desired rates. The dynamic inversion controller uses approximate short period and roll axis dynamics. The yaw axis controller is a sideslip rate command system. Methods are described to reduce the cross-coupling effect and maintain adequate tracking errors for control surface failures. The aerodynamic failure destabilizes the pitching moment due to angle of attack. The results show that control of the aircraft with the neural networks is easier (more damped) than without the neural networks. Simulation results show neural network augmentation of the controller improves performance with aerodynamic and control surface failures in terms of tracking error and cross-coupling reduction.
Burken, John J.↗
Adaptive Control Using Neural Network Augmentation for a Modified F-15 Aircraft
Description of the performance of a simplified dynamic inversion controller with neural network augmentation follows. Simulation studies focus on the results with and without neural network adaptation through the use of an F-15 aircraft simulator that has been modified to include canards. Simulated control law performance with a surface failure, in addition to an aerodynamic failure, is presented. The aircraft, with adaptation, attempts to minimize the inertial cross-coupling effect of the failure (a control derivative anomaly associated with a jammed control surface). The dynamic inversion controller calculates necessary surface commands to achieve desired rates. The dynamic inversion controller uses approximate short period and roll axis dynamics. The yaw axis controller is a sideslip rate command system. Methods are described to reduce the cross-coupling effect and maintain adequate tracking errors for control surface failures. The aerodynamic failure destabilizes the pitching moment due to angle of attack. The results show that control of the aircraft with the neural networks is easier (more damped) than without the neural networks. Simulation results show neural network augmentation of the controller improves performance with aerodynamic and control surface failures in terms of tracking error and cross-coupling reduction.
Burken, John J.↗
Robust inverse kinematics using damped least squares with dynamic weighting
This paper presents a general method for calculating the inverse kinematics with singularity and joint limit robustness for both redundant and non-redundant serial-link manipulators. Damped least squares inverse of the Jacobian is used with dynamic weighting matrices in approximating the solution. This reduces specific joint differential vectors. The algorithm gives an exact solution away from the singularities and joint limits, and an approximate solution at or near the singularities and/or joint limits. The procedure is here implemented for a six d.o.f. teleoperator and a well behaved slave manipulator resulted under teleoperational control.
Schinstock, D. E.↗
Computing Fault Displacements from Surface Deformations
Simplex is a computer program that calculates locations and displacements of subterranean faults from data on Earth-surface deformations. The calculation involves inversion of a forward model (given a point source representing a fault, a forward model calculates the surface deformations) for displacements, and strains caused by a fault located in isotropic, elastic half-space. The inversion involves the use of nonlinear, multiparameter estimation techniques. The input surface-deformation data can be in multiple formats, with absolute or differential positioning. The input data can be derived from multiple sources, including interferometric synthetic-aperture radar, the Global Positioning System, and strain meters. Parameters can be constrained or free. Estimates can be calculated for single or multiple faults. Estimates of parameters are accompanied by reports of their covariances and uncertainties. Simplex has been tested extensively against forward models and against other means of inverting geodetic data and seismic observations. This work
Lyzenga, Gregory↗
Transfer Function Models and Sensitivity Analysis
In many situations, real or induced flaws such as tight cracks with known morphology cannot be manufactured in part configuration specimens or in real parts. Typically, fatigue cracks are manufactured in simple geometry specimens such as flat plates, dog-bone shaped flat or round specimens. If a nondestructive evaluation (NDE) technique is required to provide a reliably detectable target flaw size denoted as α_(90/95) for induced flaws in real part, then the direct method for qualifying the NDE procedure is to use the appropriate induced flaw specimens and perform probability of detection analysis using these flaws. This can be described as direct POD demonstration testing, which may follow guidelines of MIL-HDBK-1823. This paper considers a case, where induced flaws are not available in part configuration specimens. Therefore, a direct POD demonstration study cannot be undertaken. In such situation, general practice for NDE procedure qualification is to use artificial flaws in simple geometry and part configuration specimens, and induced flaws in the chosen simple geometry specimens. Signal response data is taken on all sets of artificial and induced flaws. NDE procedure testing on induced flaw in simple geometry specimen is called NDE demonstration testing here. A transfer function NDE procedure qualification method for forward case calculates predicted induced flaw size for demonstration using a chosen target flaw size. Another transfer function method for inverse case, calculates the target flaw size using a given demonstration flaw size. The transfer function analysis assumes relationship of artificial flaw signal responses in real parts and simple geometry specimens; and induced flaw responses in simple geometry specimens to induced flaws in real parts. The signal response transfer relationships model should be defined before transfer function models can be devised. Assuming that the signal response transfer relationships model is valid, forward and inverse case transfer function methods have been devised. Because of lack of signal response data from induced flaws in real part, 90/95% POD/confidence (P/C) cannot be demonstrated directly. However, the transfer function method may be assessed using simulation to evaluate whether the resulting target flaw size or demonstration flaw size provides adequate confidence to the assumed signal response transfer relationships model. Therefore, the transfer function approach is a risk assessment approach. Both the signal response transfer relationships model and the transfer function model are important in managing risk in results provided by the transfer function NDE technique qualification or assessment. The signal response transfer relationships model needs to be validated with empirical data and then transfer function model needs to be validated for desired P/C on case-by-case basis.
Nondestructive evaluation↗
Radiative Transfer and Retrievals in EOF Domain
The Infrared Atmospheric Sounding Interferometer (IASI) is a hyperspectral sensor with 8461 spectral channels and a nominal spectral resolution of 0.25 cm(sup -1). It is computationally intensive to perform radiative transfer calculations and inversions using all these channels. We will present a Principal Component-based Radiative Transfer Model (PCRTM) and a retrieval algorithm which perform all the necessary calculations in EOF domain. Since the EOFs are orthogonal to each other, only about 100 principal components are needed to represent the information content of the 8461 channels. The PCRTM provides the EOF coefficients and associated derivatives with respect to atmospheric and surface parameters needed by the inversion algorithm. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance and a priori information as constraints. The retrieved parameters include atmospheric temperature, moisture and ozone profiles, cloud parameters, surface skin temperature, and surface emissivities. To make the retrieval system even more compact and stable. The atmospheric vertical profiles are compressed into the EOF space as well. The surface emissivities are also compressed into EOF space.
Liu, Xu↗
Potential flow around two-dimensional airfoils using a singular integral method
The problem of potential flow around two-dimensional airfoils is solved by using a new singular integral method. The potential flow equations for incompressible potential flow are written in a singular integral equation. The equation is solved at N collocation points on the airfoil surface. A unique feature of this method is that the airfoil geometry is specified as an independent variable in the exact integral equation. Compared to other numerical methods, the present calculation procedure is much simpler and gives remarkable accuracy for many body shapes. An advantage of the present method is that it allows the inverse design calculation and the results are extremely accurate.
Nguyen, Yves↗
Enhancing Spin Filters by Use of Bulk Inversion Asymmetry
Theoretical calculations have shown that the degrees of spin polarization in proposed nonmagnetic semiconductor resonant tunneling spin filters could be increased through exploitation of bulk inversion asymmetry (BIA). These enhancements would be effected through suitable orientation of spin collectors (or spin-polarization- inducing lateral electric fields), as described below. Spin filters -- more precisely, sources of spin-polarized electron currents -- have been sought for research on, and development of, the emerging technological discipline of spintronics (spin-transport electronics). The proposed spin filters were to be based on the Rashba effect, which is an energy splitting of what would otherwise be degenerate quantum states, caused by a spinorbit interaction in conjunction with a structural-inversion asymmetry (SIA) in the presence of interfacial electric fields in a semiconductor heterostructure. The magnitude of the energy split is proportional to the electron wave number. In a spin filter, the spin-polarized currents produced by the Rashba effect would be extracted by quantum-mechanical resonant tunneling.
Ting, David↗
Polarimetry with multiple mirror telescopes
The polarizations of multiple mirror telescopes are calculated using Mueller calculus. It is found that the Multiple Mirror Telescope (MMT) produces a constant depolarization that is a function of wavelength and independent of sky position. The efficiency and crosstalk are modeled and experimentally verified. The two- and four-mirror new generation telescopes are found to produce sinusoidal depolarization for which an accurate interpretation of the incident Stokes vector requires inverse matrix calculations. Finally, the depolarization of f/1 paraboloids is calculated and found to be less than 0.1 percent at 3000 A.
West, S. C.↗
Stratospheric Temperatures and Water Loss From Moist Greenhouse Atmospheres of Earth-Like Planets
A radiative-convective climate model is used to calculate stratospheric temperatures and water vapor concentrations for ozone-free atmospheres warmer than that of modern Earth. Cold, dry stratospheres are predicted at low surface temperatures, in agreement with recent 3D calculations. However, at surface temperatures above 350 K, the stratosphere warms and water vapor becomes a major upper atmospheric constituent, allowing water to be lost by photodissociation and hydrogen escape. Hence, a moist greenhouse explanation for loss of water from Venus, or some exoplanet receiving a comparable amount of stellar radiation, remains a viable hypothesis. Temperatures in the upper parts of such atmospheres are well below those estimated for a gray atmosphere, and this factor should be taken into account when performing inverse climate calculations to determine habitable zone boundaries using 1D models.
Astrobiology↗
An Efficient and Accurate Algorithm for Computing Grid-Averaged Solar Fluxes for Horizontally Inhomogeneous Clouds
A computationally efficient method is presented to account for the horizontal cloud inhomogeneity by using a radiatively equivalent plane parallel homogeneous (PPH) cloud. The algorithm can accurately match the calculations of the reference (rPPH) independent column approximation (ICA) results, but use only the same computational time required for a single plane parallel computation. The effective optical depth of this synthetic sPPH cloud is derived by exactly matching the direct transmission to that of the inhomogeneous ICA cloud. The ffective9 scattering asymmetry factor is found from a pre-calculated albedo inverse look-up-table that is allowed to vary over the range from -1.0 to 1.0. In the special cases of conservative scattering and total absorption, the synthetic method is exactly equivalent to the ICA, with only a small bias (about 0.2% in flux) relative to ICA due to imperfect interpolation in using the look-up tables. In principle, the ICA albedo can be approximated accurately regardless of cloud inhomogeneity. For a more complete comparison, the broadband shortwave albedo and transmission calculated from the synthetic sPPH cloud and averaged over all incident directions, have the RMS biases of 0.26% and 0.76%, respectively, for inhomogeneous clouds over a wide variation of particle size. The advantages of the synthetic PPH method are that (1) it is not required that all the cloud subcolumns have uniform microphysical characteristic, (2) it is applicable to any 1D radiative transfer scheme, and (3) it can handle arbitrary cloud optical depth distributions and an arbitrary number of cloud subcolumns with uniform computational efficiency.
cloud inhomogeneity↗