Search NASASearch

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 127 records · Page 7

Numerical computations on one-dimensional inverse scattering problems

An approximate method to determine the index of refraction of a dielectric obstacle is presented. For simplicity one dimensional models of electromagnetic scattering are treated. The governing equations yield a second order boundary value problem, in which the index of refraction appears as a functional parameter. The availability of reflection coefficients yield two additional boundary conditions. The index of refraction by a k-th order spline which can be written as a linear combination of B-splines is approximated. For N distinct reflection coefficients, the resulting N boundary value problems yield a system of N nonlinear equations in N unknowns which are the coefficients of the B-splines.

Dunn, M. H.

Structure of magnetic field in the solar wind

This work is concerned with empirical data on magnetic field in the solar wind in frame of a concept of dissipative solar wind, developed in papers (Solar Wind 7 Conf., Pergamon Press, 1992, 165 and 1992 STEP/5th COSPAR Coll. Pergamon Press, 1994, 117; 235; 803). Interplanetary magnetic fields should be classified with respect to their origin. It is very important for all the theoretical problems from the necessity to specify correctly boundary and initial conditions: the magnetic field must be sewed with its source. One should select the field, connected directly with the Sun (stretched out from it), and the field of moving electric currents. It occured central in discussion about the velocity of Alfven waves, probably warming up the solar wind, relative to the Sun, the magnetic field and solar wind plasma. The selection problem corresponds to an inverse problem and obviously has no single solution. The dissipative model of the solar wind introduce the slipping and leakage of plasma relative to magnetic field. There are no 'interplanetary current sheets' in it. But temporal fluctuations from the filamentation of electric currents play the key role. As a whole, the new concept requires the re-interpretation of main objects in the interplanetary magnetic field.

Chertkov, A. D.

The ill-posed nature of the satellite temperature retrieval problem and the limits of retrievability

The inverse problem of satellite temperature profile retrieval is well known to be ill-posed. This means that not only is a vertical temperature profile solution not unique, but that two solutions can be very different from each other. A set of atmosphere-like, and true atmospheric examples of significantly dissimilar inverse solutions, were sought and found, using an 11-channel simulated HIRS sounding radiometer. Using the Riemann-Lebesgue Lemma for guidance, it is shown that simultaneous, numerical solutions of an atmospheric character may differ by as much as 10 K between 10-1000 mb. However, an empirical search for dissimilar solutions in the natural atmosphere reveals an extremely low probability of finding two significantly different RAOBs which produce radiance measurements whose differences cannot be resolved by the satellite radiometer. The empirical results are used to derive a first estimate of the limits of retrievability, analogous to the limits of predictability derivable from the ill-posed nature of the numerical weather prediction problem.

Thompson, Owen E.

Simple formulas for designing an offset multibeam parabolic reflector

Theoretical methods (computer programs) which are available for analyzing reflector performance are of the 'forward' type. Performance parameters can be calculated after the configuration of the reflector is given. In many practical applications, however, the problem is often of the 'inverse' type. The reflector has to be designed on the basis of a given set of performance parameters. Attention is given to a systematic procedure based on simple formulas for solving such an inverse problem. Because of the simplicity of the formulas, the final results obtained from the procedure do not have great accuracy. Thus, they represent only a preliminary design, which can be refined by repeatedly using the forward-type programs.

Lee, S.-W.

Bayesian Approach to the Joint Inversion of Gravity and Magnetic Data, with Application to the Ismenius Area of Mars

This viewgraph presentation reviews a Bayesian approach to the inversion of gravity and magnetic data with specific application to the Ismenius Area of Mars. Many inverse problems encountered in geophysics and planetary science are well known to be non-unique (i.e. inversion of gravity the density structure of a body). In hopes of reducing the non-uniqueness of solutions, there has been interest in the joint analysis of data. An example is the joint inversion of gravity and magnetic data, with the assumption that the same physical anomalies generate both the observed magnetic and gravitational anomalies. In this talk, we formulate the joint analysis of different types of data in a Bayesian framework and apply the formalism to the inference of the density and remanent magnetization structure for a local region in the Ismenius area of Mars. The Bayesian approach allows prior information or constraints in the solutions to be incorporated in the inversion, with the "best" solutions those whose forward predictions most closely match the data while remaining consistent with assumed constraints. The application of this framework to the inversion of gravity and magnetic data on Mars reveals two typical challenges - the forward predictions of the data have a linear dependence on some of the quantities of interest, and non-linear dependence on others (termed the "linear" and "non-linear" variables, respectively). For observations with Gaussian noise, a Bayesian approach to inversion for "linear" variables reduces to a linear filtering problem, with an explicitly computable "error" matrix. However, for models whose forward predictions have non-linear dependencies, inference is no longer given by such a simple linear problem, and moreover, the uncertainty in the solution is no longer completely specified by a computable "error matrix". It is therefore important to develop methods for sampling from the full Bayesian posterior to provide a complete and statistically consistent picture of model uncertainty, and what has been learned from observations. We will discuss advanced numerical techniques, including Monte Carlo Markov

data analysis

New Interstellar Dust Models Consistent with Interstellar Extinction, Emission and Abundances Constraints

We present new interstellar dust models that are consistent with both, the FUV to near-IR extinction and infrared (IR) emission measurements from the diffuse interstellar medium. The models are characterized by different dust compositions and abundances. The problem we solve consists of determining the size distribution of the various dust components of the model. This problem is a typical ill-posed inversion problem which we solve using the regularization approach. We reproduce the Li Draine (2001, ApJ, 554, 778) results, however their model requires an excessive amount of interstellar silicon (48 ppM of hydrogen compared to the 36 ppM available for an ISM of solar composition) to be locked up in dust. We found that dust models consisting of PAHs, amorphous silicate, graphite, and composite grains made up from silicates, organic refractory, and water ice, provide an improved fit to the extinction and IR emission measurements, while still requiring a subsolar amount of silicon to be in the dust. This research was supported by NASA Astrophysical Theory Program NRA 99-OSS-01.

Zubko, V.

Off-Axis Nulling Transfer Function Measurement: A First Assessment

We want to study a polychromatic inverse problem method with nulling interferometers to obtain information on the structures of the exozodiacal light. For this reason, during the first semester of 2013, thanks to the support of the consortium PERSEE, we launched a campaign of laboratory measurements with the nulling interferometric test bench PERSEE, operating with 9 spectral channels between J and K bands. Our objective is to characterise the transfer function, i.e. the map of the null as a function of wavelength for an off-axis source, the null being optimised on the central source or on the source photocenter. We were able to reach on-axis null depths better than 10(exp −4). This work is part of a broader project aiming at creating a simulator of a nulling interferometer in which typical noises of a real instrument are introduced. We present here our first results.

interferometric test bench PERSEE

Wave tilt sounding of multilayered structures

The relationship between the wave tilt and the electrical parameters of a multilayered structure is investigated. Particular emphasis is placed on the inverse problem associated with the sounding planetary surfaces. An inversion technique, based on multifrequency wave tilt, is proposed and demonstrated with several computer models. It is determined that there is close agreement between the electrical parameters used in the models and those in the inversion values.

Warne, L.

Regularization for Atmospheric Temperature Retrieval Problems

Passive remote sensing of the atmosphere is used to determine the atmospheric state. A radiometer measures microwave emissions from earth's atmosphere and surface. The radiance measured by the radiometer is proportional to the brightness temperature. This brightness temperature can be used to estimate atmospheric parameters such as temperature and water vapor content. These quantities are of primary importance for different applications in meteorology, oceanography, and geophysical sciences. Depending on the range in the electromagnetic spectrum being measured by the radiometer and the atmospheric quantities to be estimated, the retrieval or inverse problem of determining atmospheric parameters from brightness temperature might be linear or nonlinear. In most applications, the retrieval problem requires the inversion of a Fredholm integral equation of the first kind making this an ill-posed problem. The numerical solution of the retrieval problem requires the transformation of the continuous problem into a discrete problem. The ill-posedness of the continuous problem translates into ill-conditioning or ill-posedness of the discrete problem. Regularization methods are used to convert the ill-posed problem into a well-posed one. In this paper, we present some results of our work in applying different regularization techniques to atmospheric temperature retrievals using brightness temperatures measured with the SSM/T-1 sensor. Simulation results are presented which show the potential of these techniques to improve temperature retrievals. In particular, no statistical assumptions are needed and the algorithms were capable of correctly estimating the temperature profile corner at the tropopause independent of the initial guess.

Velez-Reyes, Miguel

An approximation theory for the identification of linear thermoelastic systems

An abstract approximation framework and convergence theory for the identification of thermoelastic systems is developed. Starting from an abstract operator formulation consisting of a coupled second order hyperbolic equation of elasticity and first order parabolic equation for heat conduction, well-posedness is established using linear semigroup theory in Hilbert space, and a class of parameter estimation problems is then defined involving mild solutions. The approximation framework is based upon generic Galerkin approximation of the mild solutions, and convergence of solutions of the resulting sequence of approximating finite dimensional parameter identification problems to a solution of the original infinite dimensional inverse problem is established using approximation results for operator semigroups. An example involving the basic equations of one dimensional linear thermoelasticity and a linear spline based scheme are discussed. Numerical results indicate how the approach might be used in a study of damping mechanisms in flexible structures.

Rosen, I. G.

Inverse methods for assessing ship-of-opportunity networks and estimating circulation and winds from tropical expendable bathythermograph data

Inverse methods for estimating the surface ciculation of the equatorial Pacific by combining a linear reduced-gravity shallow-water model with the Tropical Ocean-Global Atmosphere ship-of-opportunity expendable bathythermograph (TOGA SOP XBT) observing program are examined. It is demonstrated that a simple linear model of the upper circulation of the equatorial Pacific can be successfully used as a weak constraint when smoothing the TOGA SOP XBT data. A circulation is sought as the weighted least squares fit to the dynamics and the data. The solution method is an expansion in representer functions, and the generalized inverse problem is thereby reduced from a functional problem to an algebraic problem for the coefficients of the representer. A specific inverse calculation using synthetic forcing and data is presented.

Bennett, Andrew F.

Vertical resolution of middle atmospheric measurements by ground-based microwave radiometry

The vertical resolution obtainable through measurement of trace constituents in the middle atmosphere by ground-based microwave spectroscopy has remained somewhat ambiguous. In order to explore this question, the Backus-Gilbert (1967, 1968, 1970) inversion technique, which automatically yields quantitative estimates of the inversion spatial resolution was applied to this particular inverse problem. This indicated that the optimum resolution of Backus-Gilbert inversions of microwave spectroscopic measurements is about 10 km. A general technique, based on inversion of delta function simulated profiles, was then developed, for evaluating the resolution of any inversion technique and applied to the Chahine (1970) inversion technique. These results indicated that the optimum resolution of the Chahine technique is about 6 or 7 km, or nearly a factor of 2 better than the equivalent Backus-Gilbert results.

Bevilacqua, Richard M.

Estimation on nonlinear damping in second order distributed parameter systems

An approximation and convergence theory for the identification of nonlinear damping in abstract wave equations is developed. It is assumed that the unknown dissipation mechanism to be identified can be described by a maximal monotone operator acting on the generalized velocity. The stiffness is assumed to be linear and symmetric. Functional analytic techniques are used to establish that solutions to a sequence of finite dimensional (Galerkin) approximating identification problems in some sense approximate a solution to the original infinite dimensional inverse problem.

Banks, H. T.