Search NASA⌕ Search

SEARCH · Search NASA

Results for “Inverse Problems”

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 559 records · Page 31

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.↗

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗

The inference of atmospheric ozone using satellite nadir measurements in the 1042/cm band

A description and detailed analysis of a technique for inferring atmospheric ozone information from satellite nadir measurements in the 1042 cm band are presented. A method is formulated for computing the emission from the lower boundary under the satellite which circumvents the difficult analytical problems caused by the presence of atmospheric clouds and the watervapor continuum absorption. The inversion equations are expanded in terms of the eigenvectors and eigenvalues of a least-squares-solution matrix, and an analysis is performed to determine the information content of the radiance measurements. Under favorable conditions there are only two pieces of independent information available from the measurements: (1) the total ozone and (2) the altitude of the primary maximum in the ozone profile.

Russell, J. M., III↗

Asteroid Lightcurve Inversion

The application of convex profile inversion (CPI) to the interpretation of asteroid lightcurves is discussed. This technique investigates the problem of extracting information about an asteroid's shape from its lightcurve. Whenever four ideal conditions are met, P is an estimator for the asteroids mean cross section C, a convex set defined as the average of all cross sections C(Z) cut by planes a distance z above the asteroids equatorial plane. C is therefore a 2-D average of the asteroids 3-D shape. The method is tested by inverting lightcurves generated analytically for geometrically scattered ellipsoids (GSE's) with semiaxes a or = B or = C. Using a defined 'distance measure' to quantify the difference between any two profiles, the deviation of P from C for GSE's as a function of lightcurve noise level, rotation phase sampling interval delta theta, and departure from ideal conditions is calibrated. The distance between P and a circle provides a gauge of the asteroid's nonsphericity and incorporates all the information contained in the lightcurve.

Ostrow, S. J.↗

On the electromagnetic scattering from infinite rectangular conducting grids

The study and development of two numerical techniques for the analysis of electromagnetic scattering from a rectangular wire mesh are described. Both techniques follow from one basic formulation and they are both solved in the spectral domain. These techniques were developed as a result of an investigation towards more efficient numerical computation for mesh scattering. These techniques are efficient for the following reasons: (a1) make use of the Fast Fourier Transform; (b2) they avoid any convolution problems by converting integrodifferential equations into algebraic equations; and (c3) they do not require inversions of any matrices. The first method, the SIT or Spectral Iteration Technique, is applied for regions where the spacing between wires is not less than two wavelengths. The second method, the SDCG or Spectral Domain Conjugate Gradient approach, can be used for any spacing between adjacent wires. A study of electromagnetic wave properties, such as reflection coefficient, induced currents and aperture fields, as functions of frequency, angle of incidence, polarization and thickness of wires is presented. Examples and comparisons or results with other methods are also included to support the validity of the new algorithms.

Christodoulou, C.↗

Connections between conventional and singular-value-based multi-variable flight control system design techniques

It is proposed that frequency-domain multivariable robustness techniques, when combined with classical multivariable procedures, can offer an additional means of evaluating FCS designs. A lateral-directional FCS for an advanced fighter is used as an example. Robustness to unstructured aircraft-input uncertainties is assessed using purely numerical singular-value procedures. Literal approximations for the singular values of the open-loop plant and controller and for the inverse return difference are shown to provide a means of decomposing and diagnosing robustness problems that are insoluble via purely numerical methods.

Mcruer, D. T.↗

An inverse dynamics approach to trajectory optimization for an aerospace plane

An inverse dynamics approach for trajectory optimization is proposed. This technique can be useful in many difficult trajectory optimization and control problems. The application of the approach is exemplified by ascent trajectory optimization for an aerospace plane. Both minimum-fuel and minimax types of performance indices are considered. When rocket augmentation is available for ascent, it is shown that accurate orbital insertion can be achieved through the inverse control of the rocket in the presence of disturbances.

Lu, Ping↗

Sparse Regression as a Sparse Eigenvalue Problem

We extend the l0-norm "subspectral" algorithms for sparse-LDA [5] and sparse-PCA [6] to general quadratic costs such as MSE in linear (kernel) regression. The resulting "Sparse Least Squares" (SLS) problem is also NP-hard, by way of its equivalence to a rank-1 sparse eigenvalue problem (e.g., binary sparse-LDA [7]). Specifically, for a general quadratic cost we use a highly-efficient technique for direct eigenvalue computation using partitioned matrix inverses which leads to dramatic x103 speed-ups over standard eigenvalue decomposition. This increased efficiency mitigates the O(n4) scaling behaviour that up to now has limited the previous algorithms' utility for high-dimensional learning problems. Moreover, the new computation prioritizes the role of the less-myopic backward elimination stage which becomes more efficient than forward selection. Similarly, branch-and-bound search for Exact Sparse Least Squares (ESLS) also benefits from partitioned matrix inverse techniques. Our Greedy Sparse Least Squares (GSLS) generalizes Natarajan's algorithm [9] also known as Order-Recursive Matching Pursuit (ORMP). Specifically, the forward half of GSLS is exactly equivalent to ORMP but more efficient. By including the backward pass, which only doubles the computation, we can achieve lower MSE than ORMP. Experimental comparisons to the state-of-the-art LARS algorithm [3] show forward-GSLS is faster, more accurate and more flexible in terms of choice of regularization

Exact Sparse Least Squares (ESLS)↗

A Comparison of Control Allocation Methods in the Presence of Parametric Model Uncertainty

When allocating redundant effectors to virtual control commands, linear (generalized inverse) allocators have historically been used on aircraft and spacecraft. While simple to implement, generalized inverses are unable to realize a significant portion of the attainable moments. To address this drawback, the control allocation problem can also be formulated as a linear programming or quadratic programming problem and solved using convex optimization based solvers. These approaches have been shown to access a larger set of attainable moments, however, little work has been done to understand the performance of convex optimization-based control allocation in the presence of parametric model uncertainty. This paper seeks to compare the performance of several control allocation approaches, including two forms of generalized inverse allocators, the pseudo inverse and minimum-variance, and the linear programming and quadratic programming approach in the presence of parametric model uncertainty. The performance of these four allocators were tested on an aircraft model in the presence of realistic parametric model uncertainty and the convex optimization approaches were shown to outperform the generalized inverses.

Luke J Miller↗

Numerical solutions of the triple-deck equations for laminar trailing-edge stall

The problem of determining the effect of laminar boundary layers on the lift of thin wings in subsonic flow at high Reynolds numbers is considered. The boundary value problem is formulated in the framework of the triple-deck theory of Brown and Stewartson. The resulting fourth-order boundary value was solved by an iterative finite-difference technique. An inverse iteration procedure provides proper treatment of the trailing-edge singularity, and asymptotic far-field expansions and coordinate stretchings are used to deal with the problem of the slow algebraic decay of the solution.

Chow, R.↗

Multidisciplinary Design Optimization of Low-Boom Supersonic Aircraft with Mission Constraints

Conceptual design of low-boom supersonic aircraft is heavily dictated by aircraft volume and lift distributions. These unique design characteristics make it a challenge to enforce mission requirements (such as static margins and trim requirements) during design optimization. This low-boom design challenge is resolved by using reversed equivalent area targets for low-fidelity low-boom inverse design and a block coordinate optimization (BCO) method for multidisciplinary design optimization (MDO). The corresponding low-boom MDO problem includes aircraft mission constraints on ranges, cruise speeds, trim for low-boom cruise, static margins for takeoff/cruise/landing, takeoff/landing field lengths, approach velocity, and tail rotation angles for trim at takeoff/landing, as well as fuselage volume constraints for passengers and main gear storage. The BCO method is developed to optimally resolve the conflicts between the low-boom inverse design objective and other design constraints. This method is successfully applied to design a low-boom supersonic configuration that carries 40 passengers, flies a low-boom mission with cruise Mach of 1.6 and range of 2,500 nm, and cruises overwater at Mach 1.8 with range of 3,600 nm. The generated configuration satisfies all specified mission constraints and has the potential to match a reversed equivalent area target with ground noise level below 70 PLdB.

Multidisciplinary design optimization↗

Identification of M2/F3 stability and control derivatives from flight data containing gust effects

The application of the maximum likelihood identification technique to M2/F3 lifting body flight data containing wind gust effects is discussed. With the objective of this effort being the identification of the stability and control derivatives, it is shown first that the output error technique (or modified Newton-Raphson) fails to fit the recorded data accurately. The means of applying the maximum likelihood technique to this problem are then discussed and the results given which indicate an accurate fit to the data. The question of derivative signs opposite to the wind tunnel values is then addressed and the results of three techniques for dealing with this problem are presented. These techniques are a priori weighting, fixing parameter values, and rank deficient inverses.

Stepner, D. E.↗

Recursive partitioned inversion of large (1500 x 1500) symmetric matrices

A recursive algorithm was designed to invert large, dense, symmetric, positive definite matrices using small amounts of computer core, i.e., a small fraction of the core needed to store the complete matrix. The described algorithm is a generalized Gaussian elimination technique. Other algorithms are also discussed for the Cholesky decomposition and step inversion techniques. The purpose of the inversion algorithm is to solve large linear systems of normal equations generated by working geodetic problems. The algorithm was incorporated into a computer program called SOLVE. In the past the SOLVE program has been used in obtaining solutions published as the Goddard earth models.

Putney, B. H.↗

Inversion layer MOS solar cells

Inversion layer (IL) Metal Oxide Semiconductor (MOS) solar cells were fabricated. The fabrication technique and problems are discussed. A plan for modeling IL cells is presented. Future work in this area is addressed.

Ho, Fat Duen↗

The Hubble Space Telescope optical systems failure report

The findings of the Hubble Space Telescope Optical Systems Board of Investigation are reported. The Board was formed to determine the cause of the flaw in the telescope, how it occurred, and why it was not detected before launch. The Board conducted its investigation to include interviews with personnel involved in the fabrication and test of the telescope, review of documentation, and analysis and test of the equipment used in the fabrication of the telescope's mirrors. The investigation proved that the primary mirror was made in the wrong shape (a 0.4-wave rms wavefront error at 632.8 nm). The primary mirror was manufactured by the Perkin-Elmer Corporation (Hughes Danbury Optical Systems, Inc.). The critical optics used as a template in shaping the mirror, the reflective null corrector (RNC), consisted of two small mirrors and a lens. This unit had been preserved by the manufacturer exactly as it was during the manufacture of the mirror. When the Board measured the RNC, the lens was incorrectly spaced from the mirrors. Calculations of the effect of such displacement on the primary mirror show that the measured amount, 1.3 mm, accounts in detail for the amount and character of the observed image blurring. No verification of the reflective null corrector's dimensions was carried out by Perkin-Elmer after the original assembly. There were, however, clear indications of the problem from auxiliary optical tests made at the time. A special optical unit called an inverse null corrector, designed to mimic the reflection from a perfect primary mirror, was built and used to align the apparatus; when so used, it clearly showed the error in the reflective null corrector. A second null corrector was used to measure the vertex radius of the finished primary mirror. It, too, clearly showed the error in the primary mirror. Both indicators of error were discounted at the time as being themselves flawed. The Perkin-Elmer plan for fabricating the primary mirror placed complete reliance on the reflective null corrector as the only test to be used in both manufacturing and verifying the mirror's surface with the required precision. This methodology should have alerted NASA management to the fragility of the process and the possibility of gross error. Such errors had been seen in other telescope programs, yet no independent tests were planned, although some simple tests to protect against major error were considered and rejected. During the critical time period, there was great concern about cost and schedule, which further inhibited consideration of independent tests.

Source record↗

Contact control for advanced applications of light weight arms

Many applications of robotic and teleoperated manipulator arms require operation in contact and non-contact regimes. This paper deals with both regimes and the transition between them with special attention given to problems of flexibility in the links and drives. This is referred to as contact control. Inverse dynamics is used to plan the tip motion of the flexible link so that the free motion can stop very near the contact surface without collision due to overshoot. Contact must occur at a very low speed since the high frequency impact forces are too sudden to be affected by any feedback generated torques applied to a joint at the other end of the link. The effect of approach velocity and surface properties are discussed. Force tracking is implemented by commands to the deflection states of the link and the contact force. This enables a natural transition between tip position and tip force control that is not possible when the arm is treated as rigid. The effect of feedback gain, force trajectory, and desired final force are of particular interest and are studied. Experimental results are presented on a one link arm and the system performance in the overall contact task is analyzed. Extension to multi-link cases with potential applications are discussed.

Book, Wayne J.↗

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↗