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 631 records · Page 35

Analytical design of intelligent machines

The problem of designing 'intelligent machines' to operate in uncertain environments with minimum supervision or interaction with a human operator is examined. The structure of an 'intelligent machine' is defined to be the structure of a Hierarchically Intelligent Control System, composed of three levels hierarchically ordered according to the principle of 'increasing precision with decreasing intelligence', namely: the organizational level, performing general information processing tasks in association with a long-term memory; the coordination level, dealing with specific information processing tasks with a short-term memory; and the control level, which performs the execution of various tasks through hardware using feedback control methods. The behavior of such a machine may be managed by controls with special considerations and its 'intelligence' is directly related to the derivation of a compatible measure that associates the intelligence of the higher levels with the concept of entropy, which is a sufficient analytic measure that unifies the treatment of all the levels of an 'intelligent machine' as the mathematical problem of finding the right sequence of internal decisions and controls for a system structured in the order of intelligence and inverse order of precision such that it minimizes its total entropy. A case study on the automatic maintenance of a nuclear plant illustrates the proposed approach.

Saridis, George N.↗

Rendering Three-Dimensional Solar Coronal Structures

An X-ray or EUV image of the corona or chromosphere is a 2D representation of an extended 3D complex for which a general inversion process is impossible. A specific model must be incorporated in order to understand the full 3D structure. We approach this problem by modeling a set of optically-thin 3D plasma flux tubes which we render these as synthetic images. The resulting images allow the interpretation of the X-ray/EUV observations to obtain information on (1) the 3D structure of X-ray images, i.e., the geometric structure of the flux tubes, and on (2) the internal structure using specific plasma characteristics, i.e., the physical structure of the flux tubes. The data-analysis technique uses magnetograms to characterize photospheric magnetic fields and extrapolation techniques to form the field lines. Using a new set of software tools, we have generated 3D flux tube structures around these field lines and integrated the plasma emission along the line of sight to obtain a rendered image. A set of individual flux-tube images is selected by a non-negative least-squares technique to Provide a match with an observed X-ray image. The scheme minimizes the squares of the differences between the synthesized image and the observed image with a non-negative constraint on the coefficients of the brightness of the individual flux-tube loops. The derived images are used to determine the specific photospheric foot points and physical data, i.e., scaling laws for densities and loop lengths. The development has led to Computer efficient integration and display software that is compatible for comparison with observations (e.g., Yohkoh SXT data, NIXT, or EIT). This analysis is important in determining directly the magnetic field configuration, which provides the structure of coronal loops, and indirectly the electric currents or waves, which provide the energy for the heating of the plasma. We have used very simple assumptions (i.e., potential magnetic fields and isothermal corona) to provide an initial test of the techniques before complex models are introduced. We have separated the physical and geometric contributions of the emission for a set of flux tubes and concentrated, in this initial study, on the geometric contributions by making approximations to the physical contributions. The initial results are consistent with the scaling laws derived from Yohkoh SXT data.

Gary, G. Allen↗

Adjoint Sensitivity Analysis of Orbital Mechanics: Application to Computations of Observables' Partials with Respect to Harmonics of the Planetary Gravity Fields

An approach is presented to the inversion of gravity fields based on evaluation of partials of observables with respect to gravity harmonics using the solution of adjoint problem of orbital dynamics of the spacecraft. Corresponding adjoint operator is derived directly from the linear operator of the linearized forward problem of orbital dynamics. The resulting adjoint problem is similar to the forward problem and can be solved by the same methods. For given highest degree N of gravity harmonics desired, this method involves integration of N adjoint solutions as compared to integration of N2 partials of the forward solution with respect to gravity harmonics in the conventional approach. Thus, for higher resolution gravity models, this approach becomes increasingly more effective in terms of computer resources as compared to the approach based on the solution of the forward problem of orbital dynamics.

gravity fields↗

TPSAS-NF1676L-33433-DND

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction, and are characterized by a complex valued frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic scattering of geometric bodies treated with acoustic liners. This work considers using either an impedance boundary condition or admittance boundary condition to allow acoustic scattering problems to be modeled with geometries consisting of both un-lined and lined surfaces; admittance is defined to be the inverse of impedance. Three different acoustic liner models will be discussed: the Extended Helmholtz Resonator Model, the Three-Parameter Impedance Model, and a Broadband Model. In the Extended Helmholtz Resonator Model and Three-Parameter Impedance Model, the liner impedance is specified at a single frequency; the Broadband Model allows for the investigation of multiple frequencies simultaneously. The impedance and admittance boundary conditions for each model will be derived and coupled with a time-domain boundary integral equation. The scattering solution will be obtained iteratively using spatial and temporal basis functions. Stability will be demonstrated through eigenvalue analysis. Moreover, fast algorithms and high performance computing are considered to both assess and reduce the computational cost of the simulation.

Michelle E Rodio↗

An inverse method for computation of structural stiffness distributions of aeroelastically optimized wings

An inverse method has been developed to compute the structural stiffness properties of wings given a specified wing loading and aeroelastic twist distribution. The method directly solves for the bending and torsional stiffness distribution of the wing using a modal representation of these properties. An aeroelastic design problem involving the use of a computational aerodynamics method to optimize the aeroelastic twist distribution of a tighter wing operating at maneuver flight conditions is used to demonstrate the application of the method. This exercise verifies the ability of the inverse scheme to accurately compute the structural stiffness distribution required to generate a specific aeroelastic twist under a specified aeroelastic load.

Schuster, David M.↗

The 3D Recognition, Generation, Fusion, Update and Refinement (RG4) Concept

This paper describes an active (real time) recognition strategy whereby information is inferred iteratively across several viewpoints in descent imagery. We will show how we use inverse theory within the context of parametric model generation, namely height and spectral reflection functions, to generate model assertions. Using this strategy in an active context implies that, from every viewpoint, the proposed system must refine its hypotheses taking into account the image and the effect of uncertainties as well. The proposed system employs probabilistic solutions to the problem of iteratively merging information (images) from several viewpoints. This involves feeding the posterior distribution from all previous images as a prior for the next view. Novel approaches will be developed to accelerate the inversion search using novel statistic implementations and reducing the model complexity using foveated vision. Foveated vision refers to imagery where the resolution varies across the image. In this paper, we allow the model to be foveated where the highest resolution region is called the foveation region. Typically, the images will have dynamic control of the location of the foveation region. For descent imagery in the Entry, Descent, and Landing (EDL) process, it is possible to have more than one foveation region. This research initiative is directed towards descent imagery in connection with NASA's EDL applications. Three-Dimensional Model Recognition, Generation, Fusion, Update, and Refinement (RGFUR or RG4) for height and the spectral reflection characteristics are in focus for various reasons, one of which is the prospect that their interpretation will provide for real time active vision for automated EDL.

Maluf, David A.↗

A general explicit solution of the central overlap problem

The matrix of the normal equations of a central overlapping-plate solution with general models for frame parameters and the star parameters is inverted explicitly in terms of the inverses of the matrices in the isolated systems for the frame and star parameters. This makes it possible to calculate from explicit formulas the numerical estimates of the star and frame parameters which would follow from a rigorous overlapping-plate solution. Also, the covariance matrix of the adjustment parameters are calculated.

Eichhorn, Heinrich↗

Interplanetary Trajectory Optimization with Powerlimited Propulsion Systems

A trajectory-optimization process is described in which the optimum­ thrust equations are derived using the calculus of variations. The mag­nitude of the thrust is constrained within an upper and a lower bound, but the thrust direction is arbitrary. This formulation allows both the constant-thrust program and the variable-thrust program to be con­sidered. For the constant-thrust program, certain propulsion-system parameters are optimized for maximum final vehicle mass. This theory has been used to study interplanetary missions to Venus and Mars using a power-limited propulsion system. Both one-way and round­ trip rendezvous trajectories are considered. The analysis employs a two-body inverse-square force-field model of three dimensions. An iterative routine used to solve the two-point boundary-value problem is described in the Appendix.

TRAJECTORY↗

Thermal Analysis Methods and Basic Heat-Transfer Data for a Turbulent Heating Flight Experiment at Mach 20 (Reentry F)

A heat-transfer experiment was flight conducted on a 5 deg half-angle cone, 396.2 cm (13 ft) in length, as it entered the sensible atmosphere under laminar, transitional, and turbulent boundary-layer conditions at a free-stream Mach number of about 20. Accurate turbulent-heat-transfer data with natural transition were obtained for correlation with theories in regions of simultaneous high Mach number, Reynolds number, enthalpy, and total-to-wall temperature ratio. Temperatures were measured at four depths through the 15.24-mm-thick (0.600-in.) beryllium wall. Experimental heating rates at 20 stations on the cone were determined independently from the outermost temperature measurement and from the temperature measurement at the second depth by a single-thermocouple inverse method and also from the temperature histories at all four depths by an integral method. The thermal data analysis procedure, associated problems, and results are presented herein.

Howard, Floyd G.↗

Stellar model chromospheres. IV - The formation of the H-epsilon feature in the sun /G2 V/ and Arcturus /K2 III/

The formation of the Balmer-series member H-epsilon in the near-red wing of the Ca II H line is discussed for two cases: the sun (H-epsilon absorption profile) and Arcturus (H-epsilon emission profile). It is shown that although the H-epsilon source functions in both stars are dominated by the Balmer-continuum radiation field through photoionizations, the line-formation problems in the two stars are quantitatively different, owing to a substantial difference in the relative importance of the stellar chromosphere temperature inversion as compared with the stellar photosphere.

Ayres, T. R.↗

High-Accuracy Compact MacCormack-Type Schemes for Computational Aeroacoustics

Using MacCormack-type methods, a new class of highly accurate compact MacCormack-type schemes is derived which does not require a tridiagonal matrix inversion to obtain the spatial derivatives. Two examples are shown, and results of these schemes for three linear and nonlinear CAA Benchmark Problems are presented.

Hixon, R.↗

Optimization of Time-Dependent Particle Tracing Using Tetrahedral Decomposition

An efficient algorithm is presented for computing particle paths, streak lines and time lines in time-dependent flows with moving curvilinear grids. The integration, velocity interpolation and step-size control are all performed in physical space which avoids the need to transform the velocity field into computational space. This leads to higher accuracy because there are no Jacobian matrix approximations or expensive matrix inversions. Integration accuracy is maintained using an adaptive step-size control scheme which is regulated by the path line curvature. The problem of cell-searching, point location and interpolation in physical space is simplified by decomposing hexahedral cells into tetrahedral cells. This enables the point location to be done analytically and substantially faster than with a Newton-Raphson iterative method. Results presented show this algorithm is up to six times faster than particle tracers which operate on hexahedral cells yet produces almost identical particle trajectories.

Kenwright, David↗

Maximally Informative Statistics for Localization and Mapping

This paper presents an algorithm for localization and mapping for a mobile robot using monocular vision and odometry as its means of sensing. The approach uses the Variable State Dimension filtering (VSDF) framework to combine aspects of Extended Kalman filtering and nonlinear batch optimization. This paper describes two primary improvements to the VSDF. The first is to use an interpolation scheme based on Gaussian quadrature to linearize measurements rather than relying on analytic Jacobians. The second is to replace the inverse covariance matrix in the VSDF with its Cholesky factor to improve the computational complexity. Results of applying the filter to the problem of localization and mapping with omnidirectional vision are presented.

Deans, Matthew C.↗

Elastic lithosphere thickness on the moon from mare tectonic features - A formal inversion

The thickness (T) of the lunar elastic lithosphere at the time (3.6 to 3.8 billion years ago) of the earliest preserved basalt flows in circular mare basins can be estimated by inverting the observed locations of extensional tectonic features in and surrounding the maria. In performing the inversion, the lithosphere is modeled as an elastic shell with a liquid interior, and the basalt load for each mare is approximated by a set of concentric cylinders. To permit solving the forward problem of placing radial limits on the positions of the rilles around a given mare, an additional parameter F, the ratio of the radial stress at the radial limits to the maximum radial stress, is introduced. T and F are chosen to give the best weighted-squares fit of the radial limits to the observations, and are used as the initial values in a linearized matrix inversion to check the resolution and estimate errors. The application of the procedure to three maria with prominent extensional features, Humorum, Orientale, and Serenitatis, gives values of T from about 40 + or - 10 to 50 + or - 10 km, and in each case the linearized matrix equation has an exact inverse.

Comer, R. P.↗

Structural control design based on reduced-order observer

An observer-based structural control design method is proposed in this paper. The method is a semi-inverse design procedure in that the control law is not designed before the observer system, but is a result that comes from the observer design. However, the observer design is not completely independent of the control design either, but seeks to yield a control law that is close to a prescribed control law. First, the observer design problem is considered as the reconstruction of a linear function of the state vector. The linear state function to be reconstructed is the given control law. Then, based on the derivation for linear state function observers, the observer design is formulated as a parameter optimization problem. The optimization objective is to generate a matrix that is close to the optimal feedback gain matrix. Based on that matrix, the form of the observer and a new control law can be determined. The semi-inverse design procedure can yield a reduced-order observer with dimension considerably smaller than that of the system. Two examples are used to demonstrate the proposed design procedure.

Su, Tzu-Jeng↗

Observability for two dimensional systems

Sufficient conditions that a two-dimensional system with output is locally observable are presented. Known results depend on time derivatives of the output and the inverse function theorem. In some cases, no informaton is provided by these theories, and one must study observability by other methods. The observability problem is dualized to the controllability problem, and the deep results of Hermes on local controllability are applied to prove a theorem concerning local observability.

Hunt, L. R.↗

Observability for two-dimensional systems

Sufficient conditions that a two-dimensinal system with output is locally observable are presented. Known results depend on time derivatives of the output and the inverse function theorem. In some cases, no information is provided by these theories, and one must study observability by other methods. The observability problem is dualized to the controllability problems, and the deep results of Hermes on local controllability are applied to prove a theorem concerning local observability.

Hunt, L. R.↗

Profiling Saturn's rings by radio occultation

The development of reconstruction algorithms that correct for diffraction effects in radio occultation measurements is described. The reciprocal Fresnel transform relationship between the complex amplitude of the observed coherent signal and the complex microwave transmittance of the rings is derived using the Huygens-Fresnel formulation of the diffraction problem. The effects of the finite data segment width, the uncertainties in the Fresnel scale, systematic phase errors in the kernel of the inverse transform, reference oscillator instabilities, and random noise measurements on the resolution of the reconstructed transmittance are analyzed. Examples of reconstructed opacity profiles for some regions of Saturn's rings derived by applying the reconstruction theory to Voyager 1 at Saturn data are presented.

Marouf, Essam A.↗