SEARCH · Search NASA

Results for “inverse problem solving”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1972

The inverse scattering problem at fixed angular momentum for nonlocal separable interactions

The problem of inverse scattering at fixed angular momentum is considered. The problem is particularized to the case of nonlocal separable interactions. A brief survey of the inverse problem for nonlocal separable interactions is presented. This problem can be solved exactly by integration. It amounts to solving singular integral equations of the Hilbert-Mushkhelishvili type, which have been studied extensively in the past and appear in many areas of physics, including theory of elasticity and dispersions relations in high energy physics.

Chadan, K.↗

NASA NTRS1973

Round-off errors in cutting plane algorithms based on the revised simplex procedure

This report statistically analyzes computational round-off errors associated with the cutting plane approach to solving linear integer programming problems. Cutting plane methods require that the inverse of a sequence of matrices be computed. The problem basically reduces to one of minimizing round-off errors in the sequence of inverses. Two procedures for minimizing this problem are presented, and their influence on error accumulation is statistically analyzed. One procedure employs a very small tolerance factor to round computed values to zero. The other procedure is a numerical analysis technique for reinverting or improving the approximate inverse of a matrix. The results indicated that round-off accumulation can be effectively minimized by employing a tolerance factor which reflects the number of significant digits carried for each calculation and by applying the reinversion procedure once to each computed inverse. If 18 significant digits plus an exponent are carried for each variable during computations, then a tolerance value of 0.1 x 10 to the minus 12th power is reasonable.

Moore, J. E.↗

NASA NTRS1988

Analysis of LIMS data by potential vorticity inversion

This paper analyzes LIMS data for a minor warming of January 1979 by potential vorticity inversion, to investigate how flows in the middle atmosphere were determined. The inversion problem of calculating the flow from the potential vorticity is solved repeatedly, including and excluding different portions of the potential vorticity, thus revealing which bits of potential vorticity were important in determining the flow at a time and location of interest. The results of the analysis of LIMS data indicate that the middle stratospheric flow is dominated by potential vorticity that is local in height and in latitude, while planetary-scale mesospheric flows are primarily induced by potential vorticity in the stratosphere. A possible explanation for this difference is proposed.

Robinson, Walter A.↗

NASA NTRS1990

Doppler radar imaging of spherical planetary surfaces

This paper describes a technique for using echo power-spectra for making radar images of spherical planetary targets. In developing the Doppler-radar imaging system, the target's reflectivity distribution is expanded in a truncated spherical harmonic series; the distribution of echo power in rotational phase and the Doppler frequency are obtained in form of a system of linear equations which express the echo spectra as functions of the series coefficients. To estimate the coefficients from an observed phase-Doppler distribution, the inversion is cast as a least-squares problem, and is solved using singular value decompositions, yielding a linear imaging system which can be fully characterized by its impulse response as a function of latitude. Simulations were designed and used to explore the capabilities of Doppler radar imaging, and the sensitivity of imaging to such factors as subradar altitude coverage and signal-to-noise ratio was investigated.

Hudson, R. Scott↗

NASA NTRS1991

Inverse problems and optimal experiment design in unsteady heat transfer processes identification

Experimental-computational methods for estimating characteristics of unsteady heat transfer processes are analyzed. The methods are based on the principles of distributed parameter system identification. The theoretical basis of such methods is the numerical solution of nonlinear ill-posed inverse heat transfer problems and optimal experiment design problems. Numerical techniques for solving problems are briefly reviewed. The results of the practical application of identification methods are demonstrated when estimating effective thermophysical characteristics of composite materials and thermal contact resistance in two-layer systems.

Artyukhin, Eugene A.↗

NASA NTRS1998

A Framework for a Supervisory Expert System for Robotic Manipulators with Joint-Position Limits and Joint-Rate Limits

This report addresses the problem of path planning and control of robotic manipulators which have joint-position limits and joint-rate limits. The manipulators move autonomously and carry out variable tasks in a dynamic, unstructured and cluttered environment. The issue considered is whether the robotic manipulator can achieve all its tasks, and if it cannot, the objective is to identify the closest achievable goal. This problem is formalized and systematically solved for generic manipulators by using inverse kinematics and forward kinematics. Inverse kinematics are employed to define the subspace, workspace and constrained workspace, which are then used to identify when a task is not achievable. The closest achievable goal is obtained by determining weights for an optimal control redistribution scheme. These weights are quantified by using forward kinematics. Conditions leading to joint rate limits are identified, in particular it is established that all generic manipulators have singularities at the boundary of their workspace, while some have loci of singularities inside their workspace. Once the manipulator singularity is identified the command redistribution scheme is used to compute the closest achievable Cartesian velocities. Two examples are used to illustrate the use of the algorithm: A three link planar manipulator and the Unimation Puma 560. Implementation of the derived algorithm is effected by using a supervisory expert system to check whether the desired goal lies in the constrained workspace and if not, to evoke the redistribution scheme which determines the constraint relaxation between end effector position and orientation, and then computes optimal gains.

Mutambara, Arthur G. O.↗

NASA NTRS1986

Kalman filtering, smoothing and recursive robot arm forward and inverse dynamics

The inverse and forward dynamics problems for multi-link serial manipulators are solved by using recursive techniques from linear filtering and smoothing theory. The pivotal step is to cast the system dynamics and kinematics as a two-point boundary-value problem. Solution of this problem leads to filtering and smoothing techniques identical to the equations of Kalman filtering and Bryson-Frazier fixed time-interval smoothing. The solutions prescribe an inward filtering recursion to compute a sequence of constraint moments and forces followed by an outward recursion to determine a corresponding sequence of angular and linear accelerations. In addition to providing techniques to compute joint accelerations from applied joint moments (and vice versa), the report provides an approach to evaluate recursively the composite multi-link system inertia matrix and its inverse. The report lays the foundation for the potential use of filtering and smoothing techniques in robot inverse and forward dynamics and in robot control design.

Rodriguez, G.↗

NASA NTRS2006

Extension of DAMAS Phased Array Processing for Spatial Coherence Determination (DAMAS-C)

The present study reports a new development of the DAMAS microphone phased array processing methodology that allows the determination and separation of coherent and incoherent noise source distributions. In 2004, a Deconvolution Approach for the Mapping of Acoustic Sources (DAMAS) was developed which decoupled the array design and processing influence from the noise being measured, using a simple and robust algorithm. In 2005, three-dimensional applications of DAMAS were examined. DAMAS has been shown to render an unambiguous quantitative determination of acoustic source position and strength. However, an underlying premise of DAMAS, as well as that of classical array beamforming methodology, is that the noise regions under study are distributions of statistically independent sources. The present development, called DAMAS-C, extends the basic approach to include coherence definition between noise sources. The solutions incorporate cross-beamforming array measurements over the survey region. While the resulting inverse problem can be large and the iteration solution computationally demanding, it solves problems no other technique can approach. DAMAS-C is validated using noise source simulations and is applied to airframe flap noise test results.

Brooks, Thomas F.↗

NASA NTRS1985

Adaptive Inverse Control for Rotorcraft Vibration Reduction

This thesis extends the Least Mean Square (LMS) algorithm to solve the mult!ple-input, multiple-output problem of alleviating N/Rev (revolutions per minute by number of blades) helicopter fuselage vibration by means of adaptive inverse control. A frequency domain locally linear model is used to represent the transfer matrix relating the higher harmonic pitch control inputs to the harmonic vibration outputs to be controlled. By using the inverse matrix as the controller gain matrix, an adaptive inverse regulator is formed to alleviate the N/Rev vibration. The stability and rate of convergence properties of the extended LMS algorithm are discussed. It is shown that the stability ranges for the elements of the stability gain matrix are directly related to the eigenvalues of the vibration signal information matrix for the learning phase, but not for the control phase. The overall conclusion is that the LMS adaptive inverse control method can form a robust vibration control system, but will require some tuning of the input sensor gains, the stability gain matrix, and the amount of control relaxation to be used. The learning curve of the controller during the learning phase is shown to be quantitatively close to that predicted by averaging the learning curves of the normal modes. For higher order transfer matrices, a rough estimate of the inverse is needed to start the algorithm efficiently. The simulation results indicate that the factor which most influences LMS adaptive inverse control is the product of the control relaxation and the the stability gain matrix. A small stability gain matrix makes the controller less sensitive to relaxation selection, and permits faster and more stable vibration reduction, than by choosing the stability gain matrix large and the control relaxation term small. It is shown that the best selections of the stability gain matrix elements and the amount of control relaxation is basically a compromise between slow, stable convergence and fast convergence with increased possibility of unstable identification. In the simulation studies, the LMS adaptive inverse control algorithm is shown to be capable of adapting the inverse (controller) matrix to track changes in the flight conditions. The algorithm converges quickly for moderate disturbances, while taking longer for larger disturbances. Perfect knowledge of the inverse matrix is not required for good control of the N/Rev vibration. However it is shown that measurement noise will prevent the LMS adaptive inverse control technique from controlling the vibration, unless the signal averaging method presented is incorporated into the algorithm.

Jacklin, Stephen A.↗

DOE OSTI2024

Paired autoencoders for likelihood-free estimation in inverse problems

Abstract We consider the solution of nonlinear inverse problems where the forward problem is a discretization of a partial differential equation. Such problems are notoriously difficult to solve in practice and require minimizing a combination of a data-fit term and a regularization term. The main computational bottleneck of typical algorithms is the direct estimation of the data misfit. Therefore, likelihood-free approaches have become appealing alternatives. Nonetheless, difficulties in generalization and limitations in accuracy have hindered their broader utility and applicability. In this work, we use a paired autoencoder framework as a likelihood-free estimator (LFE) for inverse problems. We show that the use of such an architecture allows us to construct a solution efficiently and to overcome some known open problems when using LFEs. In particular, our framework can assess the quality of the solution and improve on it if needed. We demonstrate the viability of our approach using examples from full waveform inversion and inverse electromagnetic imaging.

Chung, Matthias (ORCID:0000000178224539)↗

NASA NTRS1974

A study to determine the usefulness of interval analysis in solving problems in celestial mechanics

This investigation was undertaken to determine the usefulness of interval analysis to numerical integration and matrix inversion techniques and to combine these results to determine the value of interval analysis in bounding computational errors in the two-body problem. Conclusions were that interval analysis may be worthwhile in certain small scale isolated problems, but its usefulness in any large scale problem is doubtful.

Walling, D.↗

NASA NTRS1992

A conjugate gradients/trust regions algorithms for training multilayer perceptrons for nonlinear mapping

This paper addresses the issue of applying a globally convergent optimization algorithm to the training of multilayer perceptrons, a class of Artificial Neural Networks. The multilayer perceptrons are trained towards the solution of two highly nonlinear problems: (1) signal detection in a multi-user communication network, and (2) solving the inverse kinematics for a robotic manipulator. The research is motivated by the fact that a multilayer perceptron is theoretically capable of approximating any nonlinear function to within a specified accuracy. The algorithm that has been employed in this study combines the merits of two well known optimization algorithms, the Conjugate Gradients and the Trust Regions Algorithms. The performance is compared to a widely used algorithm, the Backpropagation Algorithm, that is basically a gradient-based algorithm, and hence, slow in converging. The performances of the two algorithms are compared with the convergence rate. Furthermore, in the case of the signal detection problem, performances are also benchmarked by the decision boundaries drawn as well as the probability of error obtained in either case.

Madyastha, Raghavendra K.↗

NASA NTRS2021

Experimental Flight Validation of the Prandtl 1933 Bell Spanload

This report describes the validation of the 1933 Prandtl bell spanload. This spanload is the minimum induced drag of a wing for a given structural weight with properties that eliminate adverse yaw. Aircraft using the Prandtl bell spanload were flown and investigated. The results of this research show that many previously held assumptions should be rethought, and the creation of aircraft using the Prandtl bell spanload will require considerable new techniques. Part of this work centered on the use of inverse methods. The usual first-step approach to a computational fluids problem is to create a geometry of the aircraft. Once this geometry exists, the computational fluids solution has been solved; however, the problem of creating the geometry still exists. A very small segment of the computational fluids world has concentrated on inverse solutions. This design approach begins with an end result - the computational fluids solution; from this end result, the geometry is sought. A more generic inverse tool was desired that would allow for the design of wings - specifically Prandtl bell spanload wings. Such a tool has been developed and allows for variable taper, aspect ratio, sweep, airfoils, and design-lift coefficients; output from the tool results in twist distribution of wings.

Albion H. Bowers↗

NASA NTRS1992

Initial value and two point boundary value solutions to the Clohessy-Wiltshire equations

The nonhomogeneous Clohessy-Wiltshire (C-W) equations are formulated and solved as an initial value problem in the form structure of linear systems theory. The state transition matrix (STM) and its inverse are obtained explicitly in both Newtonian and Hamiltonian form. It is shown that the STM for the C-2 equations possesses a special property making its inverse easily obtainable. Since solutions to the C-W equations are needed in two-point boundary value form to construct a good mission design tool for orbit transfer, the Lambert problem is solved in the context of the C-W equations.

Mullins, Larry D.↗

NASA NTRS1987

A direct algorithm for solution of incompressible three-dimensional unsteady Navier-Stokes equations

A direct, implicit, numerical solution algorithm, with second-order accuracy in space and time, is constructed for the three-dimensional unsteady incompressible Navier-Stokes equations formulated in terms of velocity and vorticity, using generalized orthogonal coordinates to achieve the accurate solution of complex viscous flow configurations. A numerically stable, efficient, direct inversion procedure is developed for the computationally intensive divergence-curl elliptic velocity problem. This overdetermined partial differential operator is first formulated as a uniquely determined, nonsingular matrix-vector problem; this aspect of the procedure is a unique feature of the present analysis. The three-dimensional vorticity-transport equation is solved by a modified factorization technique which completely eliminates the need for any block-matrix inversions and only scalar tridiagonal matrices need to be inverted. The method is applied to the test problem of the three-dimensional flow within a shear-driven cubical box. Coherent streamwise vortex structures are observed within the steady-state flow at Re = 100.

Osswald, G. A.↗

DOE OSTI2024

Overview of Large Helical Device experiments of basic plasma physics for solving crucial issues in reaching burning plasma conditions

Recently, experiments on basic plasma physics issues for solving future problems in fusion energy have been performed on a Large Helical Device. There are several problems to be solved in future devices for fusion energy. Emerging issues in burning plasma are: alpha-channeling (ion heating by alpha particles), turbulence and transport in electron dominant heating helium ash exhaust, reduction of the divertor heat load. To solve these problems, understanding the basic plasma physics of (1) wave–particle interaction through (inverse) Landau damping, (2) characteristics of electron-scale (high-k) turbulence, (3) ion mixing and the isotope effect, and (4) turbulence spreading and detachment, is necessary. This overview discusses the experimental studies on these issues and turbulent transport in multi-ion plasma and other issues in the appendix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

Total ozone and aerosol optical depths inferred from radiometric measurements in the Chappuis absorption band

A second-derivative smoothing technique, commonly used in inversion work, is applied to the problem of inferring total columnar ozone amounts and aerosol optical depths. The application is unique in that the unknowns may be solved for directly without employing standard inversion methods. It is shown, however, that by employing inversion constraints, better solutions are normally obtained. The method is quite versatile and able to deal with varying total ozone and various aerosol size distributions. The technique is applied first in simulation, then to 119 days of measurements taken in Tucson, Arizona, that are compared to TOMS values for the same dates. The technique is also applied to two measurements taken at Mauna Loa, Hawaii, for which Dobson ozone values are available in addition to the TOMS values, and the results agree to within 15 percent. It is also shown through simulations that additional information can be obtained from measurements outside the Chappuis band. This approach reduces the bias and spread of the estimated total ozone and is unique in that it uses measurements from both the Chappuis and Huggins absorption bands.

Flittner, D. E.↗

NASA NTRS1993

A novel pulse technique for thermal diffusivity determination of high-temperature levitated materials - Inverse analysis and experimental set-up

In order to determine the thermal diffusivity of materials at high temperatures, a two-step extended flash technique, which is applicable to levitated spherical samples, is proposed. The containerless flash method is modeled as an axisymmetric transient conduction heat transfer problem within the sphere subjected to a radiative boundary condition on its surface. The problem of nonlinearity which arises from the radiative heat transfer boundary condition is solved by replacing it with the measured time-dependent surface temperature data, thus giving rise to an 'Inverse Analysis'. Upon obtaining the analytic solution for the temperature field, the determination of the thermal diffusivity turns into a minimization problem, whereby the difference between the actual surface boundary condition and the analytic solution is minimized. In performing the proposed experiments, there is a need to undertake a cool-down experiment which should be governed by the lumped parameter analysis.

Shen, F.↗