Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 109 records · Page 6

Improved Bayesian regularization of inverse problems in vibrations and acoustics using noise-only measurements

Here, this paper studies Tikhonov regularization (ridge regression) parameter selection for problems in vibrations and acoustics. The selection method is based on a popular Bayesian method, but it incorporates measurements of sensor noise. The regularization parameter is closely related to the ratio of system input energy to noise energy, so noise measurements inform the inference procedure and improve parameter identification. In cases where standard Bayesian regularization identifies zero as the optimal regularization parameter, noise measurements guarantee a unique nonzero optimum. Sufficient theoretical criteria are developed for this guarantee. The method is verified in even-determined and under-determined configurations in an acoustic source localization simulation and a vibration load identification experiment. It is shown to yield significant improvements over existing empirical Bayesian regularization. Improvements are larger in the even-determined case and smaller in the under-determined case, wherein the inverse solution is less sensitive to the regularization parameter.

42 ENGINEERING↗

Sequential ensemble transform for Bayesian inverse problems

In this work, we present the Sequential Ensemble Transform (SET) method, an approach for generating approximate samples from a Bayesian posterior distribution. The method explores the posterior distribution by solving a sequence of discrete optimal transport problems to produce a series of transport plans which map prior samples to posterior samples. We prove that the sequence of Dirac mixture distributions produced by the SET method converges weakly to the true posterior as the sample size approaches infinity. Furthermore, our numerical results indicate that, when compared to standard Sequential Monte Carlo (SMC) methods, the SET approach is more robust to the choice of Markov mutation kernels and requires less computational efforts to reach a similar accuracy when used to explore complex posterior distributions. Finally, we describe adaptive schemes that allow to completely automate the use of the SET method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The P prime(f) to N(h) inversion problem in ionospheric soundings

A general review is given of the inversion techniques used to derive the ionospheric electron density N as a function of altitude h from group path P versus frequency f measurements obtained by vertical incidence ionospheric sounders. The paper discusses the medium under investigation, the experimental techniques used to obtain the P(f) data, the theoretical considerations leading to the integral equation relating P(f) to N(h), and the assumptions made in the inversion process. The lamination inversion technique is then presented, with special attention given to mathematical difficulties arising from discontinuities in the P(f) function, infinities in the integrand, and in some cases unknown integration limits. Methods outlined for minimizing the uncertainties due to discontinuities include the use of redundant information - that is the two distinct P(f) functions available for a given N(h) profile - and the use of models based on statistical data. Mathematical procedures are discussed that increase significantly the efficiency and accuracy of the required numerical integrations. The accuracy of the inversion technique is deduced by comparing the resulting N(h) profile with N(h) data obtained by simultaneous but independent observations.

Jackson, J. E.↗

Alpert multi-wavelets for functional inverse problems: direct optimization and deep learning

Computational engineering models often contain unknown entities (e.g. parameters, initial and boundary conditions) that require estimation from other measured observable data. Estimating such unknown entities is challenging when they involve spatio-temporal fields because such functional variables often require an infinite-dimensional representation. Here, we address this problem by transforming an unknown functional field using Alpert wavelet bases and truncating the resulting spectrum. Hence the problem reduces to the estimation of few coefficients that can be performed using common optimization methods. We apply this method on a one-dimensional heat transfer problem where we estimate the heat source field varying in both time and space. The observable data is comprised of temperature measured at several thermocouples in the domain. This latter is composed of either copper or stainless steel. The optimization using our method based on wavelets is able to estimate the heat source with an error between 5% and 7%. We analyze the effect of the domain material and number of thermocouples as well as the sensitivity to the initial guess of the heat source. Finally, we estimate the unknown heat source using a different approach based on deep learning techniques where we consider the input and output of a multi-layer perceptron in wavelet form. We find that this deep learning approach is more accurate than the optimization approach with errors below 4%.

97 MATHEMATICS AND COMPUTING↗

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

The inverse problem: Ocean tides derived from earth tide observations

Indirect mapping ocean tides by means of land and island-based tidal gravity measurements is presented. The inverse scheme of linear programming is used for indirect mapping of ocean tides. Open ocean tides were measured by the numerical integration of Laplace's tidal equations.

Kuo, J. T.↗

The inverse problem of constructing a gravimetric geoid

Computation of a single geoidal height from gravity acceleration data formally requires that the latter be known everywhere on the earth. A computational procedure based on linear inverse theory for estimating geoidal heights from incomplete sets of data is presented. The same scheme can be used to estimate gravity accelerations from altimetry-derived geoids. The systematic error owing to lack of data and the choice of a particular inverse operator is described by using resolution functions and their spherical harmonic expansions. An rms value of this error is also estimated by assuming a spectrum for the unknown geoid. The influence of the size of the data region, the spacing between data, the filtering applied to the data, and the model weighting function chosen are all quantified in a spherical geometry. The examples presented show that when low degree spherical harmonic coefficients are available - from satellite orbit analysis - a band-passed version of the geoid can be constructed from local gravity data, even with a relatively restricted data set.

Zlotnicki, V.↗

Estimating reaction parameters in mechanism-enabled population balance models of nanoparticle size distributions: A Bayesian inverse problem approach

In order to quantitatively predict nano- as well as other particle-size distributions, one needs to have both a mathematical model and estimates of the parameters that appear in these models. Here, we show how one can use Bayesian inversion to obtain statistical estimates for the parameters that appear in recently derived mechanism-enabled population balance models (ME-PBM) of nanoparticle growth. The Bayesian approach addresses the question of “how well do we know our parameters, along with their uncertainties?.” The results reveal that Bayesian inversion statistical analysis on an example, prototype $\mathrm{lr(0)_n}$ nanoparticle formation system allows one to estimate not just the most likely rate constants and other parameter values, but also their SDs, confidence intervals, and other statistical information. Moreover, knowing the reliability of the mechanistic model's parameters in turn helps inform one about the reliability of the proposed mechanism, as well as the reliability of its predictions. Importantly, the paper can also be seen as a tutorial with the additional goal of achieving a “Gold Standard” Bayesian inversion ME-PBM benchmark that others can use as a control to check their own use of this methodology for other systems of interest throughout nature. Overall, the results provide strong support for the hypothesis that there is substantial value in using a Bayesian inversion methodology for parameter estimation in particle formation systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inversion problem for ion-atom differential elastic scattering.

The paper describes a practical application of Remler's (1971) method by which one constructs a set of phase shifts from high resolution measurements of the differential elastic scattering of protons by rare-gas atoms. These JWKB phase shifts are then formally inverted to determine the corresponding intermolecular potentials. The validity of the method is demonstrated by comparing an intermolecular potential obtained by direct inversion of experimental data with a fairly accurate calculation by Wolniewicz (1965).

Rich, W. G.↗

Sedimentary Facies Analysis Using AVIRIS Data: A Geophysical Inverse Problem

AVIRIS data can be used to quantitatively analyze and map sedimentary lithofacies. The observed radiance spectra can be reduced to 'apparent reflectance' spectra by topographic and reflectance characterization of several field sites within the image. These apparent reflectance spectra correspond to the true reflectance at each pixel, multiplied by an unknown illumination factor (ranging in value from zero to one). The spatial abundance patterns of spectrally defined lithofacies and the unknown illumination factors can be simultaneously derived using constrained linear spectral unmixing methods. Estimates of the minimum uncertainty in the final results (due to noise, instrument resolutions, degree of illumination and mixing systematics) can be made by forward and inverse modeling. Specific facies studies in the Rattlesnake Hills region of Wyoming illustrate the successful application of these methods.

Boardmann, Joe W.↗