Science.gov

Sample records for inverse problem arising

  1. Generalized emissivity inverse problem.

    PubMed

    Ming, DengMing; Wen, Tao; Dai, XianXi; Dai, JiXin; Evenson, William E

    2002-04-01

    Inverse problems have recently drawn considerable attention from the physics community due to of potential widespread applications [K. Chadan and P. C. Sabatier, Inverse Problems in Quantum Scattering Theory, 2nd ed. (Springer Verlag, Berlin, 1989)]. An inverse emissivity problem that determines the emissivity g(nu) from measurements of only the total radiated power J(T) has recently been studied [Tao Wen, DengMing Ming, Xianxi Dai, Jixin Dai, and William E. Evenson, Phys. Rev. E 63, 045601(R) (2001)]. In this paper, a new type of generalized emissivity and transmissivity inverse (GETI) problem is proposed. The present problem differs from our previous work on inverse problems by allowing the unknown (emissivity) function g(nu) to be temperature dependent as well as frequency dependent. Based on published experimental information, we have developed an exact solution formula for this GETI problem. A universal function set suggested for numerical calculation is shown to be robust, making this inversion method practical and convenient for realistic calculations. PMID:12005916

  2. Inverse Problems of Thermoelectricity

    NASA Astrophysics Data System (ADS)

    Anatychuk, L. I.; Luste, O. J.; Kuz, R. V.; Strutinsky, M. N.

    2011-05-01

    Classical thermoelectricity is based on the use of the Seebeck and Thomson effects that occur in the near-contact areas between n- and p-type materials. A conceptually different approach to thermoelectric power converter design that is based on the law of thermoelectric induction of currents is also known. The efficiency of this approach has already been demonstrated by its first applications. More than 10 basically new types of thermoelements were discovered with properties that cannot be achieved by thermocouple power converters. Therefore, further development of this concept is of practical interest. This paper provides a classification and theory for solving the inverse problems of thermoelectricity that form the basis for devising new thermoelement types. Computer methods for their solution for anisotropic and inhomogeneous media are elaborated. Regularities related to thermoelectric current excitation in anisotropic and inhomogeneous media are established. The possibility of obtaining eddy currents of a particular configuration through control of the temperature field and material parameters for the creation of new thermo- element types is demonstrated for three-dimensional (3D) models of anisotropic and inhomogeneous media.

  3. Boundary estimation problems arising in thermal tomography

    NASA Technical Reports Server (NTRS)

    Banks, H. T.; Kojima, Fumio; Winfree, W. P.

    1989-01-01

    Problems on the identification of two-dimensional spatial domains arising in the detection and characterization of structural flaws in materials are considered. For a thermal diffusion system with external boundary input, observations of the temperature on the surface are used in a output least squares approach. Parameter estimation techniques based on the method of mappings are discussed and approximation schemes are developed based on a finite element Galerkin approach. Theoretical convergence results for computational techniques are given and the results are applied to experimental data for the identification of flaws in the thermal testing of materials.

  4. Inverse problem in hydrogeology

    NASA Astrophysics Data System (ADS)

    Carrera, Jesús; Alcolea, Andrés; Medina, Agustín; Hidalgo, Juan; Slooten, Luit J.

    2005-03-01

    The state of the groundwater inverse problem is synthesized. Emphasis is placed on aquifer characterization, where modelers have to deal with conceptual model uncertainty (notably spatial and temporal variability), scale dependence, many types of unknown parameters (transmissivity, recharge, boundary conditions, etc.), nonlinearity, and often low sensitivity of state variables (typically heads and concentrations) to aquifer properties. Because of these difficulties, calibration cannot be separated from the modeling process, as it is sometimes done in other fields. Instead, it should be viewed as one step in the process of understanding aquifer behavior. In fact, it is shown that actual parameter estimation methods do not differ from each other in the essence, though they may differ in the computational details. It is argued that there is ample room for improvement in groundwater inversion: development of user-friendly codes, accommodation of variability through geostatistics, incorporation of geological information and different types of data (temperature, occurrence and concentration of isotopes, age, etc.), proper accounting of uncertainty, etc. Despite this, even with existing codes, automatic calibration facilitates enormously the task of modeling. Therefore, it is contended that its use should become standard practice. L'état du problème inverse des eaux souterraines est synthétisé. L'accent est placé sur la caractérisation de l'aquifère, où les modélisateurs doivent jouer avec l'incertitude des modèles conceptuels (notamment la variabilité spatiale et temporelle), les facteurs d'échelle, plusieurs inconnues sur différents paramètres (transmissivité, recharge, conditions aux limites, etc.), la non linéarité, et souvent la sensibilité de plusieurs variables d'état (charges hydrauliques, concentrations) des propriétés de l'aquifère. A cause de ces difficultés, le calibrage ne peut êtreséparé du processus de modélisation, comme c'est le

  5. Inverse scattering problems with multi-frequencies

    NASA Astrophysics Data System (ADS)

    Bao, Gang; Li, Peijun; Lin, Junshan; Triki, Faouzi

    2015-09-01

    This paper is concerned with computational approaches and mathematical analysis for solving inverse scattering problems in the frequency domain. The problems arise in a diverse set of scientific areas with significant industrial, medical, and military applications. In addition to nonlinearity, there are two common difficulties associated with the inverse problems: ill-posedness and limited resolution (diffraction limit). Due to the diffraction limit, for a given frequency, only a low spatial frequency part of the desired parameter can be observed from measurements in the far field. The main idea developed here is that if the reconstruction is restricted to only the observable part, then the inversion will become stable. The challenging task is how to design stable numerical methods for solving these inverse scattering problems inspired by the diffraction limit. Recently, novel recursive linearization based algorithms have been presented in an attempt to answer the above question. These methods require multi-frequency scattering data and proceed via a continuation procedure with respect to the frequency from low to high. The objective of this paper is to give a brief review of these methods, their error estimates, and the related mathematical analysis. More attention is paid to the inverse medium and inverse source problems. Numerical experiments are included to illustrate the effectiveness of these methods.

  6. Optimization and geophysical inverse problems

    SciTech Connect

    Barhen, J.; Berryman, J.G.; Borcea, L.; Dennis, J.; de Groot-Hedlin, C.; Gilbert, F.; Gill, P.; Heinkenschloss, M.; Johnson, L.; McEvilly, T.; More, J.; Newman, G.; Oldenburg, D.; Parker, P.; Porto, B.; Sen, M.; Torczon, V.; Vasco, D.; Woodward, N.B.

    2000-10-01

    A fundamental part of geophysics is to make inferences about the interior of the earth on the basis of data collected at or near the surface of the earth. In almost all cases these measured data are only indirectly related to the properties of the earth that are of interest, so an inverse problem must be solved in order to obtain estimates of the physical properties within the earth. In February of 1999 the U.S. Department of Energy sponsored a workshop that was intended to examine the methods currently being used to solve geophysical inverse problems and to consider what new approaches should be explored in the future. The interdisciplinary area between inverse problems in geophysics and optimization methods in mathematics was specifically targeted as one where an interchange of ideas was likely to be fruitful. Thus about half of the participants were actively involved in solving geophysical inverse problems and about half were actively involved in research on general optimization methods. This report presents some of the topics that were explored at the workshop and the conclusions that were reached. In general, the objective of a geophysical inverse problem is to find an earth model, described by a set of physical parameters, that is consistent with the observational data. It is usually assumed that the forward problem, that of calculating simulated data for an earth model, is well enough understood so that reasonably accurate synthetic data can be generated for an arbitrary model. The inverse problem is then posed as an optimization problem, where the function to be optimized is variously called the objective function, misfit function, or fitness function. The objective function is typically some measure of the difference between observational data and synthetic data calculated for a trial model. However, because of incomplete and inaccurate data, the objective function often incorporates some additional form of regularization, such as a measure of smoothness

  7. Kapteyn series arising in radiation problems

    NASA Astrophysics Data System (ADS)

    Lerche, I.; Tautz, R. C.

    2008-01-01

    In discussing radiation from multiple point charges or magnetic dipoles, moving in circles or ellipses, a variety of Kapteyn series of the second kind arises. Some of the series have been known in closed form for a hundred years or more, others appear not to be available to analytic persuasion. This paper shows how 12 such generic series can be developed to produce either closed analytic expressions or integrals that are not analytically tractable. In addition, the method presented here may be of benefit when one has other Kapteyn series of the second kind to consider, thereby providing an additional reason to consider such series anew.

  8. Computationally efficient Bayesian inference for inverse problems.

    SciTech Connect

    Marzouk, Youssef M.; Najm, Habib N.; Rahn, Larry A.

    2007-10-01

    Bayesian statistics provides a foundation for inference from noisy and incomplete data, a natural mechanism for regularization in the form of prior information, and a quantitative assessment of uncertainty in the inferred results. Inverse problems - representing indirect estimation of model parameters, inputs, or structural components - can be fruitfully cast in this framework. Complex and computationally intensive forward models arising in physical applications, however, can render a Bayesian approach prohibitive. This difficulty is compounded by high-dimensional model spaces, as when the unknown is a spatiotemporal field. We present new algorithmic developments for Bayesian inference in this context, showing strong connections with the forward propagation of uncertainty. In particular, we introduce a stochastic spectral formulation that dramatically accelerates the Bayesian solution of inverse problems via rapid evaluation of a surrogate posterior. We also explore dimensionality reduction for the inference of spatiotemporal fields, using truncated spectral representations of Gaussian process priors. These new approaches are demonstrated on scalar transport problems arising in contaminant source inversion and in the inference of inhomogeneous material or transport properties. We also present a Bayesian framework for parameter estimation in stochastic models, where intrinsic stochasticity may be intermingled with observational noise. Evaluation of a likelihood function may not be analytically tractable in these cases, and thus several alternative Markov chain Monte Carlo (MCMC) schemes, operating on the product space of the observations and the parameters, are introduced.

  9. Testing Times: Problems Arising from Misdiagnosis.

    ERIC Educational Resources Information Center

    Vialle, Wilma; Konza, Deslea

    1997-01-01

    Three case studies illustrate problems in the identification of gifted students when tests are not used appropriately. The paper concludes that testing must occur within the context of intensive observations of and discussions with the child and family. The importance of all teachers receiving training in gifted education is stressed. (DB)

  10. The continuation inverse problem revisited

    NASA Astrophysics Data System (ADS)

    Huestis, Stephen P.

    1998-06-01

    The non-uniqueness of the continuation of a finite collection of harmonic potential field data to a level surface in the source-free region forces its treatment as an inverse problem. A formalism is proposed for the construction of continuation functions which are extremal by various measures. The problem is cast in such a form that the inverse problem solution is the potential function on the lowest horizontal surface above all sources, serving as the boundary function for the Dirichlet problem in the upper half-plane. The desired continuation, at the higher level of interest, must then be in the range of the upward continuation operator acting on this boundary function, rather than being allowed the full freedom of itself being part of a Dirichlet problem boundary function. Extremal solutions minimize non-linear functionals of the continuation function, which are re-expressed as different functionals of the boundary function. A crux of the method is that there is no essential distinction between the upward and downward continuation inverse problems to levels above or below data locations. Casting the optimization as a Lagrange multiplier problem leads to an integral equation for the boundary function, which is readily solved in the Fourier domain for a certain class of functionals. The desired extremal continuation is then given by upward continuation. It is found that for some functionals, application of the Lagrange multiplier theorem requires a further restriction on the set of allowable boundary functions: bandlimitedness is a natural choice for the continuation problem. With this imposition, the theory is developed in detail for semi-norm functionals penalizing departure from a constant potential, in the 2-norm and Sobelev norm senses, and illustrated by application for a small synthetic Deep Tow magnetic field data set.

  11. Uncertainty quantification for ice sheet inverse problems

    NASA Astrophysics Data System (ADS)

    Petra, N.; Ghattas, O.; Stadler, G.; Zhu, H.

    2011-12-01

    Modeling the dynamics of polar ice sheets is critical for projections of future sea level rise. Yet, there remain large uncertainties in the basal boundary conditions and in the non-Newtonian constitutive relations employed within ice sheet models. In this presentation, we consider the problem of estimating uncertainty in the solution of (large-scale) ice sheet inverse problems within the framework of Bayesian inference. Computing the general solution of the inverse problem-i.e., the posterior probability density-is intractable with current methods on today's computers, due to the expense of solving the forward model (3D full Stokes flow with nonlinear rheology) and the high dimensionality of the uncertain parameters (which are discretizations of the basal slipperiness field and the Glen's law exponent field). However, under the assumption of Gaussian noise and prior probability densities, and after linearizing the parameter-to-observable map, the posterior density becomes Gaussian, and can therefore be characterized by its mean and covariance. The mean is given by the solution of a nonlinear least squares optimization problem, which is equivalent to a deterministic inverse problem with appropriate interpretation and weighting of the data misfit and regularization terms. To obtain this mean, we solve a deterministic ice sheet inverse problem; here, we infer parameters arising from discretizations of basal slipperiness and rheological exponent fields. For this purpose, we minimize a regularized misfit functional between observed and modeled surface flow velocities. The resulting least squares minimization problem is solved using an adjoint-based inexact Newton method, which uses first and second derivative information. The posterior covariance matrix is given (in the linear-Gaussian case) by the inverse of the Hessian of the least squares cost functional of the deterministic inverse problem. Direct computation of the Hessian matrix is prohibitive, since it would

  12. Estimating uncertainties in complex joint inverse problems

    NASA Astrophysics Data System (ADS)

    Afonso, Juan Carlos

    2016-04-01

    Sources of uncertainty affecting geophysical inversions can be classified either as reflective (i.e. the practitioner is aware of her/his ignorance) or non-reflective (i.e. the practitioner does not know that she/he does not know!). Although we should be always conscious of the latter, the former are the ones that, in principle, can be estimated either empirically (by making measurements or collecting data) or subjectively (based on the experience of the researchers). For complex parameter estimation problems in geophysics, subjective estimation of uncertainty is the most common type. In this context, probabilistic (aka Bayesian) methods are commonly claimed to offer a natural and realistic platform from which to estimate model uncertainties. This is because in the Bayesian approach, errors (whatever their nature) can be naturally included as part of the global statistical model, the solution of which represents the actual solution to the inverse problem. However, although we agree that probabilistic inversion methods are the most powerful tool for uncertainty estimation, the common claim that they produce "realistic" or "representative" uncertainties is not always justified. Typically, ALL UNCERTAINTY ESTIMATES ARE MODEL DEPENDENT, and therefore, besides a thorough characterization of experimental uncertainties, particular care must be paid to the uncertainty arising from model errors and input uncertainties. We recall here two quotes by G. Box and M. Gunzburger, respectively, of special significance for inversion practitioners and for this session: "…all models are wrong, but some are useful" and "computational results are believed by no one, except the person who wrote the code". In this presentation I will discuss and present examples of some problems associated with the estimation and quantification of uncertainties in complex multi-observable probabilistic inversions, and how to address them. Although the emphasis will be on sources of uncertainty related

  13. Inverse problem in radionuclide transport

    SciTech Connect

    Yu, C.

    1988-01-01

    The disposal of radioactive waste must comply with the performance objectives set forth in 10 CFR 61 for low-level waste (LLW) and 10 CFR 60 for high-level waste (HLW). To determine probable compliance, the proposed disposal system can be modeled to predict its performance. One of the difficulties encountered in such a study is modeling the migration of radionuclides through a complex geologic medium for the long term. Although many radionuclide transport models exist in the literature, the accuracy of the model prediction is highly dependent on the model parameters used. The problem of using known parameters in a radionuclide transport model to predict radionuclide concentrations is a direct problem (DP); whereas the reverse of DP, i.e., the parameter identification problem of determining model parameters from known radionuclide concentrations, is called the inverse problem (IP). In this study, a procedure to solve IP is tested, using the regression technique. Several nonlinear regression programs are examined, and the best one is recommended. 13 refs., 1 tab.

  14. Minimax approach to inverse problems of geophysics

    NASA Astrophysics Data System (ADS)

    Balk, P. I.; Dolgal, A. S.; Balk, T. V.; Khristenko, L. A.

    2016-03-01

    A new approach is suggested for solving the inverse problems that arise in the different fields of applied geophysics (gravity, magnetic, and electrical prospecting, geothermy) and require assessing the spatial region occupied by the anomaly-generating masses in the presence of different types of a priori information. The interpretation which provides the maximum guaranteed proximity of the model field sources to the real perturbing object is treated as the best interpretation. In some fields of science (game theory, economics, operations research), the decision-making principle that lies in minimizing the probable losses which cannot be prevented if the situation develops by the worst-case scenario is referred to as minimax. The minimax criterion of choice is interesting as, instead of being confined to the indirect (and sometimes doubtful) signs of the "optimal" solution, it relies on the actual properties of the information in the results of a particular interpretation. In the hierarchy of the approaches to the solution of the inverse problems of geophysics ordered by the volume and quality of the retrieved information about the sources of the field, the minimax approach should take special place.

  15. Inverse problem for Bremsstrahlung radiation

    SciTech Connect

    Voss, K.E.; Fisch, N.J.

    1991-10-01

    For certain predominantly one-dimensional distribution functions, an analytic inversion has been found which yields the velocity distribution of superthermal electrons given their Bremsstrahlung radiation. 5 refs.

  16. Parabolic Perturbation of a Nonlinear Hyperbolic Problem Arising in Physiology

    NASA Astrophysics Data System (ADS)

    Colli, P.; Grasselli, M.

    We study a transport-diffusion initial value problem where the diffusion codlicient is "small" and the transport coefficient is a time function depending on the solution in a nonlinear and nonlocal way. We show the existence and the uniqueness of a weak solution of this problem. Moreover we discuss its asymptotic behaviour as the diffusion coefficient goes to zero, obtaining a well-posed first-order nonlinear hyperbolic problem. These problems arise from mathematical models of muscle contraction in the framework of the sliding filament theory.

  17. BOOK REVIEW: Inverse Problems. Activities for Undergraduates

    NASA Astrophysics Data System (ADS)

    Yamamoto, Masahiro

    2003-06-01

    This book is a valuable introduction to inverse problems. In particular, from the educational point of view, the author addresses the questions of what constitutes an inverse problem and how and why we should study them. Such an approach has been eagerly awaited for a long time. Professor Groetsch, of the University of Cincinnati, is a world-renowned specialist in inverse problems, in particular the theory of regularization. Moreover, he has made a remarkable contribution to educational activities in the field of inverse problems, which was the subject of his previous book (Groetsch C W 1993 Inverse Problems in the Mathematical Sciences (Braunschweig: Vieweg)). For this reason, he is one of the most qualified to write an introductory book on inverse problems. Without question, inverse problems are important, necessary and appear in various aspects. So it is crucial to introduce students to exercises in inverse problems. However, there are not many introductory books which are directly accessible by students in the first two undergraduate years. As a consequence, students often encounter diverse concrete inverse problems before becoming aware of their general principles. The main purpose of this book is to present activities to allow first-year undergraduates to learn inverse theory. To my knowledge, this book is a rare attempt to do this and, in my opinion, a great success. The author emphasizes that it is very important to teach inverse theory in the early years. He writes; `If students consider only the direct problem, they are not looking at the problem from all sides .... The habit of always looking at problems from the direct point of view is intellectually limiting ...' (page 21). The book is very carefully organized so that teachers will be able to use it as a textbook. After an introduction in chapter 1, sucessive chapters deal with inverse problems in precalculus, calculus, differential equations and linear algebra. In order to let one gain some insight

  18. An inverse problem in thermal imaging

    NASA Technical Reports Server (NTRS)

    Bryan, Kurt; Caudill, Lester F., Jr.

    1994-01-01

    This paper examines uniqueness and stability results for an inverse problem in thermal imaging. The goal is to identify an unknown boundary of an object by applying a heat flux and measuring the induced temperature on the boundary of the sample. The problem is studied both in the case in which one has data at every point on the boundary of the region and the case in which only finitely many measurements are available. An inversion procedure is developed and used to study the stability of the inverse problem for various experimental configurations.

  19. An inverse acoustic waveguide problem in the time domain

    NASA Astrophysics Data System (ADS)

    Monk, Peter; Selgas, Virginia

    2016-05-01

    We consider the problem of locating an obstacle in a waveguide from time domain measurements of causal waves. More precisely, we assume that we are given the scattered field due to point sources placed on a surface located inside the waveguide away from the obstacle, where the scattered field is measured on the same surface. From this multi-static scattering data we wish to determine the position and shape of an obstacle in the waveguide. To deal with this inverse problem, we adapt and analyze the time domain linear sampling method. This involves proving new time domain estimates for the forward problem, as well as analyzing several time domain operators arising in the inversion scheme. We also implement the inversion algorithm and provide numerical results in two-dimensions using synthetic data.

  20. A scatterometry inverse problem in optical mask metrology

    NASA Astrophysics Data System (ADS)

    Model, R.; Rathsfeld, A.; Gross, H.; Wurm, M.; Bodermann, B.

    2008-11-01

    We discuss the solution of the inverse problem in scatterometry i.e. the determination of periodic surface structures from light diffraction patterns. With decreasing details of lithography masks, increasing demands on metrology techniques arise. By scatterometry as a non-imaging indirect optical method critical dimensions (CD) like side-wall angles, heights, top and bottom widths are determined. The numerical simulation of diffraction is based on the finite element solution of the Helmholtz equation. The inverse problem seeks to reconstruct the grating geometry from measured diffraction patterns. The inverse operator maps efficiencies of diffracted plane wave modes to the grating parameters. We employ a Newton type iterative method to solve the resulting minimum problem. The reconstruction quality surely depends on the angles of incidence, on the wave lengths and/or the number of propagating scattered wave modes and will be discussed by numerical examples.

  1. Comparative evolution of the inverse problems (Introduction to an interdisciplinary study of the inverse problems)

    NASA Technical Reports Server (NTRS)

    Sabatier, P. C.

    1972-01-01

    The progressive realization of the consequences of nonuniqueness imply an evolution of both the methods and the centers of interest in inverse problems. This evolution is schematically described together with the various mathematical methods used. A comparative description is given of inverse methods in scientific research, with examples taken from mathematics, quantum and classical physics, seismology, transport theory, radiative transfer, electromagnetic scattering, electrocardiology, etc. It is hoped that this paper will pave the way for an interdisciplinary study of inverse problems.

  2. Direct and Inverse problems in Electrocardiography

    NASA Astrophysics Data System (ADS)

    Boulakia, M.; Fernández, M. A.; Gerbeau, J. F.; Zemzemi, N.

    2008-09-01

    We present numerical results related to the direct and the inverse problems in electrocardiography. The electrical activity of the heart is described by the bidomain equations. The electrocardiograms (ECGs) recorded in different points on the body surface are obtained by coupling the bidomain equation to a Laplace equation in the torso. The simulated ECGs are quite satisfactory. As regards the inverse problem, our goal is to estimate the parameters of the bidomain-torso model. Here we present some preliminary results of a parameter estimation for the torso model.

  3. Inverse problem of electro-seismic conversion

    NASA Astrophysics Data System (ADS)

    Chen, Jie; Yang, Yang

    2013-11-01

    When a porous rock is saturated with an electrolyte, electrical fields are coupled with seismic waves via the electro-seismic conversion. Pride (1994 Phys. Rev. B 50 15678-96) derived the governing models, in which Maxwell equations are coupled with Biot's equations through the electro-kinetic mobility parameter. The inverse problem of the linearized electro-seismic conversion consists in two steps, namely the inversion of Biot's equations and the inversion of Maxwell equations. We analyze the reconstruction of conductivity and electro-kinetic mobility parameter in Maxwell equations with internal measurements, while the internal measurements are provided by the results of the inversion of Biot's equations. We show that knowledge of two internal data based on well-chosen boundary conditions uniquely determines these two parameters. Moreover, a Lipschitz-type stability is proved based on the same sets of well-chosen boundary conditions.

  4. Solving inversion problems with neural networks

    NASA Technical Reports Server (NTRS)

    Kamgar-Parsi, Behzad; Gualtieri, J. A.

    1990-01-01

    A class of inverse problems in remote sensing can be characterized by Q = F(x), where F is a nonlinear and noninvertible (or hard to invert) operator, and the objective is to infer the unknowns, x, from the observed quantities, Q. Since the number of observations is usually greater than the number of unknowns, these problems are formulated as optimization problems, which can be solved by a variety of techniques. The feasibility of neural networks for solving such problems is presently investigated. As an example, the problem of finding the atmospheric ozone profile from measured ultraviolet radiances is studied.

  5. MAP estimators and their consistency in Bayesian nonparametric inverse problems

    NASA Astrophysics Data System (ADS)

    Dashti, M.; Law, K. J. H.; Stuart, A. M.; Voss, J.

    2013-09-01

    We consider the inverse problem of estimating an unknown function u from noisy measurements y of a known, possibly nonlinear, map {G} applied to u. We adopt a Bayesian approach to the problem and work in a setting where the prior measure is specified as a Gaussian random field μ0. We work under a natural set of conditions on the likelihood which implies the existence of a well-posed posterior measure, μy. Under these conditions, we show that the maximum a posteriori (MAP) estimator is well defined as the minimizer of an Onsager-Machlup functional defined on the Cameron-Martin space of the prior; thus, we link a problem in probability with a problem in the calculus of variations. We then consider the case where the observational noise vanishes and establish a form of Bayesian posterior consistency for the MAP estimator. We also prove a similar result for the case where the observation of {G}(u) can be repeated as many times as desired with independent identically distributed noise. The theory is illustrated with examples from an inverse problem for the Navier-Stokes equation, motivated by problems arising in weather forecasting, and from the theory of conditioned diffusions, motivated by problems arising in molecular dynamics.

  6. A Stochastic Problem Arising in the Storage of Radioactive Waste

    SciTech Connect

    Williams, M.M.R.

    2004-07-15

    Nuclear waste drums can contain a collection of radioactive components of uncertain activity and randomly dispersed in position. This implies that the dose-rate at the surface of different drums in a large assembly of similar drums can have significant variations according to the physical makeup and configuration of the waste components. The present paper addresses this problem by treating the drum, and its waste, as a stochastic medium. It is assumed that the sources in the drum contribute a dose-rate to some external point. The strengths and positions are chosen by random numbers, the dose-rate is calculated and, from several thousand realizations, a probability distribution for the dose-rate is obtained. It is shown that a very close approximation to the dose-rate probability function is the log-normal distribution. This allows some useful statistical indicators, which are of environmental importance, to be calculated with little effort.As an example of a practical situation met in the storage of radioactive waste containers, we study the problem of 'hotspots'. These arise in drums in which most of the activity is concentrated on one radioactive component and hence can lead to the possibility of large surface dose-rates. It is shown how the dose-rate, the variance, and some other statistical indicators depend on the relative activities on the sources. The results highlight the importance of such hotspots and the need to quantify their effect.

  7. Numerical linear algebra for reconstruction inverse problems

    NASA Astrophysics Data System (ADS)

    Nachaoui, Abdeljalil

    2004-01-01

    Our goal in this paper is to discuss various issues we have encountered in trying to find and implement efficient solvers for a boundary integral equation (BIE) formulation of an iterative method for solving a reconstruction problem. We survey some methods from numerical linear algebra, which are relevant for the solution of this class of inverse problems. We motivate the use of our constructing algorithm, discuss its implementation and mention the use of preconditioned Krylov methods.

  8. PUBLISHER'S ANNOUNCEMENT: New developments for Inverse Problems

    NASA Astrophysics Data System (ADS)

    2006-12-01

    2006 has proved to be a very successful year for Inverse Problems. After an increase for the fourth successive year, we achieved our highest impact factor to date, 1.541 (Source: 2005 ISI® Journal Citation Report), and the Editorial Board is keen to build on this success by continuing to improve the service we offer to our readers and authors. The Board has observed that Inverse Problems receives very few Letters to the Editor submissions, and that moreover those that we do receive rarely conform to the requirements for Letters to the Editor set out in the journal's editorial policy. The Board has therefore decided to merge the current Letters to the Editor section into our regular Papers section, which will now accommodate all research articles that meet the journal's high quality standards. Any submissions that would previously have been Letters to the Editor are still very welcome as Papers, and can be submitted by e-mail to ip@iop.org or online using our online submissions form at authors.iop.org/submit. Inverse Problems' processing times are already among the fastest in the field—on average, authors receive our decision on their paper in less than three months. Thanks to our easy-to-use online refereeing system, publishing a Paper is now just as fast as publishing a Letter to the Editor, and we are striving to ensure that the journal's high standards are applied consistently to all our Papers, maintaining Inverse Problems' position as the leading journal in the field. Our highly acclaimed Topical Review section will also continue and grow; providing timely insights into the development of all topical fields within Inverse Problems. We have many exciting Topical Reviews currently in preparation for 2007 and will continue to commission articles at the cutting edge of research. We look forward to receiving your contributions and to continuing to provide the best publication service available.

  9. An efficient method for inverse problems

    NASA Technical Reports Server (NTRS)

    Daripa, Prabir

    1987-01-01

    A new inverse method for aerodynamic design of subcritical airfoils is presented. The pressure distribution in this method can be prescribed in a natural way, i.e. as a function of arclength of the as yet unknown body. This inverse problem is shown to be mathematically equivalent to solving a single nonlinear boundary value problem subject to known Dirichlet data on the boundary. The solution to this problem determines the airfoil, the free stream Mach number M(sub x) and the upstream flow direction theta(sub x). The existence of a solution for any given pressure distribution is discussed. The method is easy to implement and extremely efficient. We present a series of results for which comparisons are made with the known airfoils.

  10. Urban surface water pollution problems arising from misconnections.

    PubMed

    Revitt, D Michael; Ellis, J Bryan

    2016-05-01

    The impacts of misconnections on the organic and nutrient loadings to surface waters are assessed using specific household appliance data for two urban sub-catchments located in the London metropolitan region and the city of Swansea. Potential loadings of biochemical oxygen demand (BOD), soluble reactive phosphorus (PO4-P) and ammoniacal nitrogen (NH4-N) due to misconnections are calculated for three different scenarios based on the measured daily flows from specific appliances and either measured daily pollutant concentrations or average pollutant concentrations for relevant greywater and black water sources obtained from an extensive review of the literature. Downstream receiving water concentrations, together with the associated uncertainties, are predicted from derived misconnection discharge concentrations and compared to existing freshwater standards for comparable river types. Consideration of dilution ratios indicates that these would need to be of the order of 50-100:1 to maintain high water quality with respect to BOD and NH4-N following typical misconnection discharges but only poor quality for PO4-P is likely to be achievable. The main pollutant loading contributions to misconnections arise from toilets (NH4-N and BOD), kitchen sinks (BOD and PO4-P) washing machines (PO4-P and BOD) and, to a lesser extent, dishwashers (PO4-P). By completely eliminating toilet misconnections and ensuring misconnections from all other appliances do not exceed 2%, the potential pollution problems due to BOD and NH4-N discharges would be alleviated but this would not be the case for PO4-P. In the event of a treatment option being preferred to solve the misconnection problem, it is shown that for an area the size of metropolitan Greater London, a sewage treatment plant with a Population Equivalent value approaching 900,000 would be required to efficiently remove BOD and NH4-N to safely dischargeable levels but such a plant is unlikely to have the capacity to deal

  11. A variational Bayesian method to inverse problems with impulsive noise

    NASA Astrophysics Data System (ADS)

    Jin, Bangti

    2012-01-01

    We propose a novel numerical method for solving inverse problems subject to impulsive noises which possibly contain a large number of outliers. The approach is of Bayesian type, and it exploits a heavy-tailed t distribution for data noise to achieve robustness with respect to outliers. A hierarchical model with all hyper-parameters automatically determined from the given data is described. An algorithm of variational type by minimizing the Kullback-Leibler divergence between the true posteriori distribution and a separable approximation is developed. The numerical method is illustrated on several one- and two-dimensional linear and nonlinear inverse problems arising from heat conduction, including estimating boundary temperature, heat flux and heat transfer coefficient. The results show its robustness to outliers and the fast and steady convergence of the algorithm.

  12. Inverse problems biomechanical imaging (Conference Presentation)

    NASA Astrophysics Data System (ADS)

    Oberai, Assad A.

    2016-03-01

    It is now well recognized that a host of imaging modalities (a list that includes Ultrasound, MRI, Optical Coherence Tomography, and optical microscopy) can be used to "watch" tissue as it deforms in response to an internal or external excitation. The result is a detailed map of the deformation field in the interior of the tissue. This deformation field can be used in conjunction with a material mechanical response to determine the spatial distribution of material properties of the tissue by solving an inverse problem. Images of material properties thus obtained can be used to quantify the health of the tissue. Recently, they have been used to detect, diagnose and monitor cancerous lesions, detect vulnerable plaque in arteries, diagnose liver cirrhosis, and possibly detect the onset of Alzheimer's disease. In this talk I will describe the mathematical and computational aspects of solving this class of inverse problems, and their applications in biology and medicine. In particular, I will discuss the well-posedness of these problems and quantify the amount of displacement data necessary to obtain a unique property distribution. I will describe an efficient algorithm for solving the resulting inverse problem. I will also describe some recent developments based on Bayesian inference in estimating the variance in the estimates of material properties. I will conclude with the applications of these techniques in diagnosing breast cancer and in characterizing the mechanical properties of cells with sub-cellular resolution.

  13. Inverse scattering problem for quantum graph vertices

    SciTech Connect

    Cheon, Taksu; Turek, Ondrej; Exner, Pavel

    2011-06-15

    We demonstrate how the inverse scattering problem of a quantum star graph can be solved by means of diagonalization of the Hermitian unitary matrix when the vertex coupling is of the scale-invariant (or Fueloep-Tsutsui) form. This enables the construction of quantum graphs with desired properties in a tailor-made fashion. The procedure is illustrated on the example of quantum vertices with equal transmission probabilities.

  14. 3D Inverse problem: Seawater intrusions

    NASA Astrophysics Data System (ADS)

    Steklova, K.; Haber, E.

    2013-12-01

    Modeling of seawater intrusions (SWI) is challenging as it involves solving the governing equations for variable density flow, multiple time scales and varying boundary conditions. Due to the nonlinearity of the equations and the large aquifer domains, 3D computations are a costly process, particularly when solving the inverse SWI problem. In addition the heads and concentration measurements are difficult to obtain due to mixing, saline wedge location is sensitive to aquifer topography, and there is general uncertainty in initial and boundary conditions and parameters. Some of these complications can be overcome by using indirect geophysical data next to standard groundwater measurements, however, the inverse problem is usually simplified, e.g. by zonation for the parameters based on geological information, steady state substitution of the unknown initial conditions, decoupling the equations or reducing the amount of unknown parameters by covariance analysis. In our work we present a discretization of the flow and solute mass balance equations for variable groundwater (GW) flow. A finite difference scheme is to solve pressure equation and a Semi - Lagrangian method for solute transport equation. In this way we are able to choose an arbitrarily large time step without losing stability up to an accuracy requirement coming from the coupled character of the variable density flow equations. We derive analytical sensitivities of the GW model for parameters related to the porous media properties and also the initial solute distribution. Analytically derived sensitivities reduce the computational cost of inverse problem, but also give insight for maximizing information in collected data. If the geophysical data are available it also enables simultaneous calibration in a coupled hydrogeophysical framework. The 3D inverse problem was tested on artificial time dependent data for pressure and solute content coming from a GW forward model and/or geophysical forward model. Two

  15. Numerical strategies for the solution of inverse problems

    NASA Astrophysics Data System (ADS)

    Hebber (Haber), Eldad

    This thesis deals with the numerical solutions of linear and nonlinear inverse problems. The goal of this thesis is to review and develop new techniques for solving such problems. In so doing, the computations tools for solving inverse problems are comprehensively studied. The thesis can be divided into two parts. In the first part, linear inverse theory is dealt with. Methods to estimate noise and efficiently invert large and full matrixes are reviewed and developed. Emphasis is given to Generalized Cross Validation (GCV) for noise estimation, and to Krylov space methods for efficient methods to invert large systems. This part is summarized by applying and comparing the methods developed on linear inverse problems which arise in gravity and tomography. In the second part of this thesis, extensive use of the linear algebra and the noise estimation methods which were developed in the first part of the thesis is made. A review of the current methods to carry out nonlinear inverse problems is given. A test example is constructed to demonstrate that these methods may fail. Next, a new algorithm for solving nonlinear inverse problems is developed. The algorithm is based on the ability to differentiate between correlated errors which comes from the linearization, and non-correlated noise which comes from the measurement. Based on these two types of noise, a regularization procedure which has two parts is developed. The first part is made of global regulation, to deal with the measurement noise, and the second part is made from a local regularization, to deal with the nonlinearity. The thesis demonstrates that GCV can be used in order to determine the measurement noise, and the Damped Gauss-Newton method can be used in order to deal with the local nonlinear terms. Another aspect of nonlinear inverse theory which is developed in this work concerns approximate sensitivities. A new formulation is suggested for the approximate sensitivities and bounds are calculated using

  16. Discrete Inverse and State Estimation Problems

    NASA Astrophysics Data System (ADS)

    Wunsch, Carl

    2006-06-01

    The problems of making inferences about the natural world from noisy observations and imperfect theories occur in almost all scientific disciplines. This book addresses these problems using examples taken from geophysical fluid dynamics. It focuses on discrete formulations, both static and time-varying, known variously as inverse, state estimation or data assimilation problems. Starting with fundamental algebraic and statistical ideas, the book guides the reader through a range of inference tools including the singular value decomposition, Gauss-Markov and minimum variance estimates, Kalman filters and related smoothers, and adjoint (Lagrange multiplier) methods. The final chapters discuss a variety of practical applications to geophysical flow problems. Discrete Inverse and State Estimation Problems is an ideal introduction to the topic for graduate students and researchers in oceanography, meteorology, climate dynamics, and geophysical fluid dynamics. It is also accessible to a wider scientific audience; the only prerequisite is an understanding of linear algebra. Provides a comprehensive introduction to discrete methods of inference from incomplete information Based upon 25 years of practical experience using real data and models Develops sequential and whole-domain analysis methods from simple least-squares Contains many examples and problems, and web-based support through MIT opencourseware

  17. Structural information in the inverse problem

    NASA Astrophysics Data System (ADS)

    Karaoulis, M.; Larson, T. H.; Ahmed, I.; Revil, A.; Thomason, J.

    2013-12-01

    Data integration in geophysics provides additional information to elucidate subsurface structure and reduce non-uniqueness of inverted models. There are several strategies for incorporating data integration into numerical models. A traditional approach is to use this information as a-priori knowledge, as an initial model or layer boundary specified on the mesh before the inversion. Although in some cases this approach has proven effective, the information is not directly incorporated into the inverse problem and might be lost in the final model. Another strategy is through joint inversion, where data and models are inverted simultaneous. The data integration comes from the joint operator as structural similarity or petrophysical equations. There are some limitations with this approach. In particular, structural similarity doesn't take into account different sensitivity patterns which differ for different geophysical methods, e.g. in a crosswell configuration electrical resistivity tomography is sensitive close to the borehole region while seismic waves are sensitive towards the center part. Moreover different methods have differing resolution. Therefore, a single joint operator might not be effective in all cases. In this work we demonstrate the use of image-guided inversion, where structural information is taken directly from a high resolution geophysical image (e.g. ground penetrating radar or seismic reflection) or from a geological cross-section. This structural information is introduced into the inverse problem through a weighted smoothing matrix, where it correlates and favors formations related to a specific structural feature and not just uniformly across the entire model. Both sharp and smooth features can be imaged and the recovered models can have a more realistic distribution of values. As an example of the method we use migrated seismic reflection images to extract the structural information and resistivity imaging to recover the resistivity

  18. Direct and inverse problems of infrared tomography.

    PubMed

    Sizikov, Valery S; Evseev, Vadim; Fateev, Alexander; Clausen, Sønnik

    2016-01-01

    The problems of infrared tomography-direct (the modeling of measured functions) and inverse (the reconstruction of gaseous medium parameters)-are considered with a laboratory burner flame as an example of an application. The two measurement modes are used: active (ON) with an external IR source and passive (OFF) without one. Received light intensities on detectors are modeled in the direct problem or measured in the experiment whereas integral equations with respect to the absorption coefficient and Planck function (which yields the temperature profile of the medium) are solved in the inverse problem with (1) modeled and (2) measured received intensities as the input data. An axisymmetric flame and parallel scanning scheme of measurements considered in this work yield singular integral equations that are solved numerically using the generalized quadrature method, spline smoothing, and Tikhonov regularization. A software package in MATLAB has been developed. Two numerical examples-with modeled and real input data-were solved. The proposed methodology avoids the necessity of elaborate determination of the absorption coefficient by direct (point) measurements or calculation using spectroscopic databases (e.g., HITRAN/HITEMP). PMID:26835642

  19. Inverse problem of nonlinear dynamical systems: a constructive approach

    SciTech Connect

    Gonzalez-Gascon, F.; Moreno-Insertis, F.; Rodriguez-Camino, E.

    1980-08-01

    A quite simple and practical method is developed for the construction of two dimensional nonlinear dynamical systems (plane vector fields) possessing an arbitrary number of given limit cycles. The method is applied to the construction of n-dimensional dynamical systems (R/sup n/ vector fields) possessing at least one limit cycle and, under certain circumstances, more than one, or even a numerable infinity. Interesting open problems arise when n is greater than two, or where more than one limit cycle appears. Our constructive algorithm for this type of inverse problem is also applied to the construction of second order differential equations (Newtonian differential equations) possessing a finite or infinite number of invariant speeds. This last problem is relevant for certain aspects of the special theory of relativity.

  20. The inverse problem of the optimal regulator.

    NASA Technical Reports Server (NTRS)

    Yokoyama, R.; Kinnen, E.

    1972-01-01

    The inverse problem of the optimal regulator is considered for a general class of multi-input systems with integral-type performance indices. A new phase variable canonical form is shown to be convenient for this analysis. The advantage of the canonical form is to separate the state variables into subvectors of directly controlled, indirectly controlled, and uncontrollable components. Necessary and sufficient conditions for optimized performance indices are given. With the nonlinearities of the system restricted to functions of the directly controlled state variables, additional results are developed about the nonnegative property of optimized loss functions.

  1. EDITORIAL: Inverse Problems' 25th year of publication Inverse Problems' 25th year of publication

    NASA Astrophysics Data System (ADS)

    2008-01-01

    2009 is Inverse Problems' 25th year of publication. In this quarter-century, the journal has established itself as the premier publication venue for inverse problems research. It has matured from its beginnings as a niche journal serving the emerging field of inverse and ill-posed problems to a monthly publication in 2009 covering all aspects of a well-established, vibrant and still-expanding subject. Along with its core readership of pure and applied mathematicians and physicists, Inverse Problems has become widely known across a broad range of researchers in areas such as geophysics, optics, radar, acoustics, communication theory, signal processing and medical imaging, amongst others. The journal's appeal to the inverse problems community and those researchers from the varied fields that encounter such problems can be attributed to our commitment to publishing only the very best papers, and to offering unique services to the community. Besides our regular research papers, which average a remarkably short five months from submission to electronic publication, we regularly publish heavily cited topical review papers and topic-specific special sections, which first appeared in 2004. These highly-downloaded invited articles focus on the latest developments and hot topics in all areas of inverse problems. No other journal in the field offers these features. I am very pleased to take Inverse Problems into its 25th year as Editor-in-Chief. The journal has an impressive tradition of scholarship, established at its inception by the founder and first Editor-in-Chief, Professor Pierre Sabatier. Professor Sabatier envisioned the journal in 1985 as providing a medium for publication of exemplary research in our intrinsically interdisciplinary field. I am glad to say that the support of our authors, readers, referees, Editors-in-Chief, Editorial Boards and Advisory Panels over the years, has resulted in Inverse Problems becoming the top publication in this field, publishing

  2. An inverse problem by boundary element method

    SciTech Connect

    Tran-Cong, T.; Nguyen-Thien, T.; Graham, A.L.

    1996-02-01

    Boundary Element Methods (BEM) have been established as useful and powerful tools in a wide range of engineering applications, e.g. Brebbia et al. In this paper, we report a particular three dimensional implementation of a direct boundary integral equation (BIE) formulation and its application to numerical simulations of practical polymer processing operations. In particular, we will focus on the application of the present boundary element technology to simulate an inverse problem in plastics processing.by extrusion. The task is to design profile extrusion dies for plastics. The problem is highly non-linear due to material viscoelastic behaviours as well as unknown free surface conditions. As an example, the technique is shown to be effective in obtaining the die profiles corresponding to a square viscoelastic extrudate under different processing conditions. To further illustrate the capability of the method, examples of other non-trivial extrudate profiles and processing conditions are also given.

  3. Inverse Variational Problem for Nonstandard Lagrangians

    NASA Astrophysics Data System (ADS)

    Saha, A.; Talukdar, B.

    2014-06-01

    In the mathematical physics literature the nonstandard Lagrangians (NSLs) were introduced in an ad hoc fashion rather than being derived from the solution of the inverse problem of variational calculus. We begin with the first integral of the equation of motion and solve the associated inverse problem to obtain some of the existing results for NSLs. In addition, we provide a number of alternative Lagrangian representations. The case studies envisaged by us include (i) the usual modified Emden-type equation, (ii) Emden-type equation with dissipative term quadratic in velocity, (iii) Lotka-Volterra model and (vi) a number of the generic equations for dissipative-like dynamical systems. Our method works for nonstandard Lagrangians corresponding to the usual action integral of mechanical systems but requires modification for those associated with the modified actions like S =∫abe L(x ,x˙ , t) dt and S =∫abL 1 - γ(x ,x˙ , t) dt because in the latter case one cannot construct expressions for the Jacobi integrals.

  4. A local-order regularization for geophysical inverse problems

    NASA Astrophysics Data System (ADS)

    Gheymasi, H. Mohammadi; Gholami, A.

    2013-11-01

    Different types of regularization have been developed to obtain stable solutions to linear inverse problems. Among these, total variation (TV) is known as an edge preserver method, which leads to piecewise constant solutions and has received much attention for solving inverse problems arising in geophysical studies. However, the method shows staircase effects and is not suitable for the models including smooth regions. To overcome the staircase effect, we present a method, which employs a local-order difference operator in the regularization term. This method is performed in two steps: First, we apply a pre-processing step to find the edge locations in the regularized solution using a properly defined minmod limiter, where the edges are determined by a comparison of the solutions obtained using different order regularizations of the TV types. Then, we construct a local-order difference operator based on the information obtained from the pre-processing step about the edge locations, which is subsequently used as a regularization operator in the final sparsity-promoting regularization. Experimental results from the synthetic and real seismic traveltime tomography show that the proposed inversion method is able to retain the smooth regions of the regularized solution, while preserving sharp transitions presented in it.

  5. PREFACE: International Conference on Inverse Problems 2010

    NASA Astrophysics Data System (ADS)

    Hon, Yiu-Chung; Ling, Leevan

    2011-03-01

    Following the first International Conference on Inverse Problems - Recent Theoretical Development and Numerical Approaches held at the City University of Hong Kong in 2002, the fifth International Conference was held again at the City University during December 13-17, 2010. This fifth conference was jointly organized by Professor Yiu-Chung Hon (Co-Chair, City University of Hong Kong, HKSAR), Dr Leevan Ling (Co-Chair, Hong Kong Baptist University, HKSAR), Professor Jin Cheng (Fudan University, China), Professor June-Yub Lee (Ewha Womans University, South Korea), Professor Gui-Rong Liu (University of Cincinnati, USA), Professor Jenn-Nan Wang (National Taiwan University, Taiwan), and Professor Masahiro Yamamoto (The University of Tokyo, Japan). It was agreed to alternate holding the conference among the above places (China, Japan, Korea, Taiwan, and Hong Kong) once every two years. The next conference has been scheduled to be held at the Southeast University (Nanjing, China) in 2012. The purpose of this series of conferences is to establish a strong collaborative link among the universities of the Asian-Pacific regions and worldwide leading researchers in inverse problems. The conference addressed both theoretical (mathematics), applied (engineering) and developmental aspects of inverse problems. The conference was intended to nurture Asian-American-European collaborations in the evolving interdisciplinary areas and it was envisioned that the conference would lead to long-term commitments and collaborations among the participating countries and researchers. There was a total of more than 100 participants. A call for the submission of papers was sent out after the conference, and a total of 19 papers were finally accepted for publication in this proceedings. The papers included in the proceedings cover a wide scope, which reflects the current flourishing theoretical and numerical research into inverse problems. Finally, as the co-chairs of the Inverse Problems

  6. Bayesian inference tools for inverse problems

    NASA Astrophysics Data System (ADS)

    Mohammad-Djafari, Ali

    2013-08-01

    In this paper, first the basics of Bayesian inference with a parametric model of the data is presented. Then, the needed extensions are given when dealing with inverse problems and in particular the linear models such as Deconvolution or image reconstruction in Computed Tomography (CT). The main point to discuss then is the prior modeling of signals and images. A classification of these priors is presented, first in separable and Markovien models and then in simple or hierarchical with hidden variables. For practical applications, we need also to consider the estimation of the hyper parameters. Finally, we see that we have to infer simultaneously on the unknowns, the hidden variables and the hyper parameters. Very often, the expression of this joint posterior law is too complex to be handled directly. Indeed, rarely we can obtain analytical solutions to any point estimators such the Maximum A posteriori (MAP) or Posterior Mean (PM). Three main tools are then can be used: Laplace approximation (LAP), Markov Chain Monte Carlo (MCMC) and Bayesian Variational Approximations (BVA). To illustrate all these aspects, we will consider a deconvolution problem where we know that the input signal is sparse and propose to use a Student-t prior for that. Then, to handle the Bayesian computations with this model, we use the property of Student-t which is modelling it via an infinite mixture of Gaussians, introducing thus hidden variables which are the variances. Then, the expression of the joint posterior of the input signal samples, the hidden variables (which are here the inverse variances of those samples) and the hyper-parameters of the problem (for example the variance of the noise) is given. From this point, we will present the joint maximization by alternate optimization and the three possible approximation methods. Finally, the proposed methodology is applied in different applications such as mass spectrometry, spectrum estimation of quasi periodic biological signals and

  7. An Inverse Problem of Derivative Security Pricing

    NASA Astrophysics Data System (ADS)

    Zhang, Guanquan; Li, Peijun

    2003-04-01

    Suppose that interest rate is governed by a stochastic differential equation, a partial differential equation for the price of bond can be derived in a similar way to the derivation of the Black-Scholes equation. Valuation of bond with implied function in the equation, which is called the risk market price of interest rate, is known as the model of bond pricing. An inverse problem of bond pricing is to determined the risk market price of interest rate implied by current prices of bonds with different expirations. In this paper, numerical algorithm to solve this system is constructed and some numerical experiments are performed. The numerical results show that the algorithm is quite efficient and robust.

  8. Stochastic inverse problems: Models and metrics

    SciTech Connect

    Sabbagh, Elias H.; Sabbagh, Harold A.; Murphy, R. Kim; Aldrin, John C.; Annis, Charles; Knopp, Jeremy S.

    2015-03-31

    In past work, we introduced model-based inverse methods, and applied them to problems in which the anomaly could be reasonably modeled by simple canonical shapes, such as rectangular solids. In these cases the parameters to be inverted would be length, width and height, as well as the occasional probe lift-off or rotation. We are now developing a formulation that allows more flexibility in modeling complex flaws. The idea consists of expanding the flaw in a sequence of basis functions, and then solving for the expansion coefficients of this sequence, which are modeled as independent random variables, uniformly distributed over their range of values. There are a number of applications of such modeling: 1. Connected cracks and multiple half-moons, which we have noted in a POD set. Ideally we would like to distinguish connected cracks from one long shallow crack. 2. Cracks of irregular profile and shape which have appeared in cold work holes during bolt-hole eddy-current inspection. One side of such cracks is much deeper than other. 3. L or C shaped crack profiles at the surface, examples of which have been seen in bolt-hole cracks. By formulating problems in a stochastic sense, we are able to leverage the stochastic global optimization algorithms in NLSE, which is resident in VIC-3D®, to answer questions of global minimization and to compute confidence bounds using the sensitivity coefficient that we get from NLSE. We will also address the issue of surrogate functions which are used during the inversion process, and how they contribute to the quality of the estimation of the bounds.

  9. Stability analysis of the inverse transmembrane potential problem in electrocardiography

    NASA Astrophysics Data System (ADS)

    Burger, Martin; Mardal, Kent-André; Nielsen, Bjørn Fredrik

    2010-10-01

    In this paper we study some mathematical properties of an inverse problem arising in connection with electrocardiograms (ECGs). More specifically, we analyze the possibility for recovering the transmembrane potential in the heart from ECG recordings, a challenge currently investigated by a growing number of groups. Our approach is based on the bidomain model for the electrical activity in the myocardium, and leads to a parameter identification problem for elliptic partial differential equations (PDEs). It turns out that this challenge can be split into two subproblems: the task of recovering the potential at the heart surface from body surface recordings; the problem of computing the transmembrane potential inside the heart from the potential determined at the heart surface. Problem (1), which can be formulated as the Cauchy problem for an elliptic PDE, has been extensively studied and is well known to be severely ill-posed. The main purpose of this paper is to prove that problem (2) is stable and well posed if a suitable prior is available. Moreover, our theoretical findings are illuminated by a series of numerical experiments. Finally, we discuss some aspects of uniqueness related to the anisotropy in the heart.

  10. The relativistic inverse stellar structure problem

    SciTech Connect

    Lindblom, Lee

    2014-01-14

    The observable macroscopic properties of relativistic stars (whose equations of state are known) can be predicted by solving the stellar structure equations that follow from Einstein’s equation. For neutron stars, however, our knowledge of the equation of state is poor, so the direct stellar structure problem can not be solved without modeling the highest density part of the equation of state in some way. This talk will describe recent work on developing a model independent approach to determining the high-density neutron-star equation of state by solving an inverse stellar structure problem. This method uses the fact that Einstein’s equation provides a deterministic relationship between the equation of state and the macroscopic observables of the stars which are composed of that material. This talk illustrates how this method will be able to determine the high-density part of the neutron-star equation of state with few percent accuracy when high quality measurements of the masses and radii of just two or three neutron stars become available. This talk will also show that this method can be used with measurements of other macroscopic observables, like the masses and tidal deformabilities, which can (in principle) be measured by gravitational wave observations of binary neutron-star mergers.

  11. Direct and Inverse Problems in Statistical Wavefields

    SciTech Connect

    Wolf, Emil

    2002-09-01

    In this report account is presented of research carried out during the period September 1, 1999-August 31, 2002 under the sponsorship of the Department of Energy, grant DE-FG02-90ER14119. The research covered several areas of modern optical physics, particularly propagation of partially coherent light and its interaction with deterministic and with random media, spectroscopy with partially coherent light, polarization properties of statistical wave fields, effects of moving diffusers on coherence and on the spectra of light transmitted and scattered by them, reciprocity inequalities involving spatial and angular correlations of partially coherent beams, spreading of partially coherent beams in-random media, inverse source problems, computed and diffraction tomography and partially coherent solitons. We have discovered a new phenomenon in an emerging field of physical optics, known as singular optics; specifically we found that the spectrum of light changes drastically in the neighborhood of points where the intensity has zero value and where, consequently, the phase becomes singular, We noted some potential applications of this phenomenon. The results of our investigations were reported in 39 publications. They are listed on pages 3 to 5. Summaries of these publications are given on pages 6-13. Scientists who have participated in this research are listed on page 14.

  12. PREFACE: International Conference on Inverse Problems 2010

    NASA Astrophysics Data System (ADS)

    Hon, Yiu-Chung; Ling, Leevan

    2011-03-01

    Following the first International Conference on Inverse Problems - Recent Theoretical Development and Numerical Approaches held at the City University of Hong Kong in 2002, the fifth International Conference was held again at the City University during December 13-17, 2010. This fifth conference was jointly organized by Professor Yiu-Chung Hon (Co-Chair, City University of Hong Kong, HKSAR), Dr Leevan Ling (Co-Chair, Hong Kong Baptist University, HKSAR), Professor Jin Cheng (Fudan University, China), Professor June-Yub Lee (Ewha Womans University, South Korea), Professor Gui-Rong Liu (University of Cincinnati, USA), Professor Jenn-Nan Wang (National Taiwan University, Taiwan), and Professor Masahiro Yamamoto (The University of Tokyo, Japan). It was agreed to alternate holding the conference among the above places (China, Japan, Korea, Taiwan, and Hong Kong) once every two years. The next conference has been scheduled to be held at the Southeast University (Nanjing, China) in 2012. The purpose of this series of conferences is to establish a strong collaborative link among the universities of the Asian-Pacific regions and worldwide leading researchers in inverse problems. The conference addressed both theoretical (mathematics), applied (engineering) and developmental aspects of inverse problems. The conference was intended to nurture Asian-American-European collaborations in the evolving interdisciplinary areas and it was envisioned that the conference would lead to long-term commitments and collaborations among the participating countries and researchers. There was a total of more than 100 participants. A call for the submission of papers was sent out after the conference, and a total of 19 papers were finally accepted for publication in this proceedings. The papers included in the proceedings cover a wide scope, which reflects the current flourishing theoretical and numerical research into inverse problems. Finally, as the co-chairs of the Inverse Problems

  13. A Riemann–Hilbert approach to the inverse problem for the Stark operator on the line

    NASA Astrophysics Data System (ADS)

    Its, A.; Sukhanov, V.

    2016-05-01

    The paper is concerned with the inverse scattering problem for the Stark operator on the line with a potential from the Schwartz class. In our study of the inverse problem, we use the Riemann–Hilbert formalism. This allows us to overcome the principal technical difficulties which arise in the more traditional approaches based on the Gel’fand–Levitan–Marchenko equations, and indeed solve the problem. We also produce a complete description of the relevant scattering data (which have not been obtained in the previous works on the Stark operator) and establish the bijection between the Schwartz class potentials and the scattering data.

  14. Inverse Modelling Problems in Linear Algebra Undergraduate Courses

    ERIC Educational Resources Information Center

    Martinez-Luaces, Victor E.

    2013-01-01

    This paper will offer an analysis from a theoretical point of view of mathematical modelling, applications and inverse problems of both causation and specification types. Inverse modelling problems give the opportunity to establish connections between theory and practice and to show this fact, a simple linear algebra example in two different…

  15. A Forward Glimpse into Inverse Problems through a Geology Example

    ERIC Educational Resources Information Center

    Winkel, Brian J.

    2012-01-01

    This paper describes a forward approach to an inverse problem related to detecting the nature of geological substrata which makes use of optimization techniques in a multivariable calculus setting. The true nature of the related inverse problem is highlighted. (Contains 2 figures.)

  16. Numerical boundary condition procedure for the transonic axisymmetric inverse problem

    NASA Technical Reports Server (NTRS)

    Shankar, V.

    1981-01-01

    Two types of boundary condition procedures for the axisymmetric inverse problem are described. One is a Neumann type boundary condition (analogous to the analysis problem) and the other is a Dirichlet type boundary conditon, both requiring special treatments to make the inverse scheme numerically stable. The dummy point concept is utilized in implementing both. Results indicate the Dirichlet type inverse boundary condition is more robust and conceptually simpler to implement than the Neumann type procedure. A few results demonstrating the powerful capability of the newly developed inverse method that can handle both shocked as well as shockless body design are included.

  17. Solving ill-posed magnetic inverse problem using a Parameterized Trust-Region Sub-problem

    NASA Astrophysics Data System (ADS)

    Abdelazeem, Maha Mohamed

    2013-06-01

    The aim of this paper is to find a plausible and stable solution for the inverse geophysical magnetic problem. Most of the inverse problems in geophysics are considered as ill-posed ones. This is not necessarily due to complex geological situations, but it may arise because of ill-conditioned kernel matrix. To deal with such ill-conditioned matrix, one may truncate the most ill part as in truncated singular value decomposition method (TSVD). In such a method, the question will be where to truncate? In this paper, for comparison, we first try the adaptive pruning algorithm for the discrete L-curve criterion to estimate the regularization parameter for TSVD method. Linear constraints have been added to the ill-conditioned matrix. The same problem is then solved using a global optimizing and regularizing technique based on Parameterized Trust Region Sub-problem (PTRS). The criteria of such technique are to choose a trusted region of the solutions and then to find the satisfying minimum to the objective function. The ambiguity is controlled mainly by proper choosing the trust region. To overcome the natural decay in kernel with depth, a specific depth weighting function is used. A Matlab-based inversion code is implemented and tested on two synthetic total magnetic fields contaminated with different levels of noise to simulate natural fields. The results of PTRS are compared with those of TSVD with adaptive pruning L-curve. Such a comparison proves the high stability of the PTRS method in dealing with potential field problems. The capability of such technique has been further tested by applying it to real data from Saudi Arabia and Italy.

  18. An inverse problem for a class of conditional probability measure-dependent evolution equations

    NASA Astrophysics Data System (ADS)

    Mirzaev, Inom; Byrne, Erin C.; Bortz, David M.

    2016-09-01

    We investigate the inverse problem of identifying a conditional probability measure in measure-dependent evolution equations arising in size-structured population modeling. We formulate the inverse problem as a least squares problem for the probability measure estimation. Using the Prohorov metric framework, we prove existence and consistency of the least squares estimates and outline a discretization scheme for approximating a conditional probability measure. For this scheme, we prove general method stability. The work is motivated by partial differential equation models of flocculation for which the shape of the post-fragmentation conditional probability measure greatly impacts the solution dynamics. To illustrate our methodology, we apply the theory to a particular PDE model that arises in the study of population dynamics for flocculating bacterial aggregates in suspension, and provide numerical evidence for the utility of the approach.

  19. Uniqueness in inverse boundary value problems for fractional diffusion equations

    NASA Astrophysics Data System (ADS)

    Li, Zhiyuan; Imanuvilov, Oleg Yu; Yamamoto, Masahiro

    2016-01-01

    We consider an inverse boundary value problem for diffusion equations with multiple fractional time derivatives. We prove the uniqueness in determining the number of fractional time-derivative terms, the orders of the derivatives and spatially varying coefficients.

  20. Inverse kinematics problem in robotics using neural networks

    NASA Technical Reports Server (NTRS)

    Choi, Benjamin B.; Lawrence, Charles

    1992-01-01

    In this paper, Multilayer Feedforward Networks are applied to the robot inverse kinematic problem. The networks are trained with endeffector position and joint angles. After training, performance is measured by having the network generate joint angles for arbitrary endeffector trajectories. A 3-degree-of-freedom (DOF) spatial manipulator is used for the study. It is found that neural networks provide a simple and effective way to both model the manipulator inverse kinematics and circumvent the problems associated with algorithmic solution methods.

  1. Inverse radiation problem in axisymmetric cylindrical scattering media

    NASA Astrophysics Data System (ADS)

    Menguc, M. P.; Manickavasagam, S.

    1993-09-01

    A semianalytical technique has been developed to solve the inverse radiation problem in absorbing and scattering cylindrical media. The radiative properties in the medium are allowed to vary radially. Isotropic, linearly anisotropic, and Rayleigh scattering phase functions are considered, and both the first- and second-order scattering of radiation are accounted for in the analysis. The angular radiosity distribution obtained from the solution of the forward problem is employed as input to the inverse analysis. A numerical inversion scheme is followed to determine the profiles of extinction coefficient and the single-scattering albedo. For an anisotropically scattering medium, the asymmetry factor is also recovered. It is shown that the method is simple and accurate, even though the inversion is limited to three- or four-layer media. This inversion procedure can easily be used in experiments to determine the effective radiative property distributions in cylindrical systems.

  2. Expectation propagation for nonlinear inverse problems – with an application to electrical impedance tomography

    SciTech Connect

    Gehre, Matthias; Jin, Bangti

    2014-02-15

    In this paper, we study a fast approximate inference method based on expectation propagation for exploring the posterior probability distribution arising from the Bayesian formulation of nonlinear inverse problems. It is capable of efficiently delivering reliable estimates of the posterior mean and covariance, thereby providing an inverse solution together with quantified uncertainties. Some theoretical properties of the iterative algorithm are discussed, and the efficient implementation for an important class of problems of projection type is described. The method is illustrated with one typical nonlinear inverse problem, electrical impedance tomography with complete electrode model, under sparsity constraints. Numerical results for real experimental data are presented, and compared with that by Markov chain Monte Carlo. The results indicate that the method is accurate and computationally very efficient.

  3. The metric-restricted inverse design problem

    NASA Astrophysics Data System (ADS)

    Acharya, Amit; Lewicka, Marta; Pakzad, Mohammad Reza

    2016-06-01

    We study a class of design problems in solid mechanics, leading to a variation on the classical question of equi-dimensional embeddability of Riemannian manifolds. In this general new context, we derive a necessary and sufficient existence condition, given through a system of total differential equations, and discuss its integrability. In the classical context, the same approach yields conditions of immersibility of a given metric in terms of the Riemann curvature tensor. In the present situation, the equations do not close in a straightforward manner, and successive differentiation of the compatibility conditions leads to a new algebraic description of integrability. We also recast the problem in a variational setting and analyze the infimum of the appropriate incompatibility energy, resembling the non-Euclidean elasticity. We then derive a Γ -convergence result for dimension reduction from 3d to 2d in the Kirchhoff energy scaling regime.

  4. Minimax theory for a class of nonlinear statistical inverse problems

    NASA Astrophysics Data System (ADS)

    Ray, Kolyan; Schmidt-Hieber, Johannes

    2016-06-01

    We study a class of statistical inverse problems with nonlinear pointwise operators motivated by concrete statistical applications. A two-step procedure is proposed, where the first step smoothes the data and inverts the nonlinearity. This reduces the initial nonlinear problem to a linear inverse problem with deterministic noise, which is then solved in a second step. The noise reduction step is based on wavelet thresholding and is shown to be minimax optimal (up to logarithmic factors) in a pointwise function-dependent sense. Our analysis is based on a modified notion of Hölder smoothness scales that are natural in this setting.

  5. Solving inverse problems of identification type by optimal control methods

    SciTech Connect

    Lenhart, S.; Protopopescu, V.; Yong, J.

    1997-05-01

    Inverse problems of identification type for nonlinear equations are considered within the framework of optimal control theory. The rigorous solution of any particular problem depends on the functional setting, type of equation, and unknown quantity (or quantities) to be determined. Here we present only the general articulations of the formalism. Compared to classical regularization methods (e.g. Tikhonov coupled with optimization schemes), our approach presents several advantages, namely: (i) a systematic procedure to solve inverse problems of identification type; (ii) an explicit expression for the approximations of the solution; and (iii) a convenient numerical solution of these approximations. {copyright} {ital 1997 American Institute of Physics.}

  6. Solving inverse problems of identification type by optimal control methods

    SciTech Connect

    Lenhart, S.; Protopopescu, V.; Jiongmin Yong

    1997-06-01

    Inverse problems of identification type for nonlinear equations are considered within the framework of optimal control theory. The rigorous solution of any particular problem depends on the functional setting, type of equation, and unknown quantity (or quantities) to be determined. Here the authors present only the general articulations of the formalism. Compared to classical regularization methods (e.g. Tikhonov coupled with optimization schemes), their approach presents several advantages, namely: (i) a systematic procedure to solve inverse problems of identification type; (ii) an explicit expression for the approximations of the solution; and (iii) a convenient numerical solution of these approximations.

  7. Correct averaging in transmission radiography: Analysis of the inverse problem

    NASA Astrophysics Data System (ADS)

    Wagner, Michael; Hampel, Uwe; Bieberle, Martina

    2016-05-01

    Transmission radiometry is frequently used in industrial measurement processes as a means to assess the thickness or composition of a material. A common problem encountered in such applications is the so-called dynamic bias error, which results from averaging beam intensities over time while the material distribution changes. We recently reported on a method to overcome the associated measurement error by solving an inverse problem, which in principle restores the exact average attenuation by considering the Poisson statistics of the underlying particle or photon emission process. In this paper we present a detailed analysis of the inverse problem and its optimal regularized numerical solution. As a result we derive an optimal parameter configuration for the inverse problem.

  8. Finite element neural networks for electromagnetic inverse problems

    NASA Astrophysics Data System (ADS)

    Ramuhalli, P.; Udpa, L.; Udpa, S.

    2002-05-01

    Iterative approaches using numerical forward models are commonly used for solving inverse problems in nondestructive evaluation. The drawbacks of these approaches include their high computational cost and the difficulty in computing gradients for updating defect profiles. This paper proposes a finite element neural network (FENN) that embeds finite element models into a neural network format. This approach enables fast and accurate solution of the forward problem. The FENN can then be used as the forward model in an iterative approach to solve the inverse problem. Gradient-based optimization methods are easily applied since the FENN provides an explicit functional mapping between the defect profile and the measured signal. Results of applying the FENN to several simple electromagnetic forward and inverse problems are presented.

  9. Inverse problem in nondestructive testing using arrayed eddy current sensors.

    PubMed

    Zaoui, Abdelhalim; Menana, Hocine; Feliachi, Mouloud; Berthiau, Gérard

    2010-01-01

    A fast crack profile reconstitution model in nondestructive testing is developed using an arrayed eddy current sensor. The inverse problem is based on an iterative solving of the direct problem using genetic algorithms. In the direct problem, assuming a current excitation, the incident field produced by all the coils of the arrayed sensor is obtained by the translation and superposition of the 2D axisymmetric finite element results obtained for one coil; the impedance variation of each coil, due to the crack, is obtained by the reciprocity principle involving the dyadic Green's function. For the inverse problem, the surface of the crack is subdivided into rectangular cells, and the objective function is expressed only in terms of the depth of each cell. The evaluation of the dyadic Green's function matrix is made independently of the iterative procedure, making the inversion very fast. PMID:22163680

  10. Inverse Problem in Nondestructive Testing Using Arrayed Eddy Current Sensors

    PubMed Central

    Zaoui, Abdelhalim; Menana, Hocine; Feliachi, Mouloud; Berthiau, Gérard

    2010-01-01

    A fast crack profile reconstitution model in nondestructive testing is developed using an arrayed eddy current sensor. The inverse problem is based on an iterative solving of the direct problem using genetic algorithms. In the direct problem, assuming a current excitation, the incident field produced by all the coils of the arrayed sensor is obtained by the translation and superposition of the 2D axisymmetric finite element results obtained for one coil; the impedance variation of each coil, due to the crack, is obtained by the reciprocity principle involving the dyadic Green’s function. For the inverse problem, the surface of the crack is subdivided into rectangular cells, and the objective function is expressed only in terms of the depth of each cell. The evaluation of the dyadic Green’s function matrix is made independently of the iterative procedure, making the inversion very fast. PMID:22163680

  11. On the Inverse Problems of Nonlinear Acoustics and Acoustic Turbulence

    NASA Astrophysics Data System (ADS)

    Gurbatov, S. N.; Rudenko, O. V.

    2015-12-01

    We consider the problem of retrieval of the radiated acoustic signal parameters from the measured wave field in some cross section of the nonlinear medium. The possibilities of solving regular and statistical inverse problems are discussed on the basis of the solution of the Burgers equation for zero and infinitesimal viscosities.

  12. Analytic semigroups: Applications to inverse problems for flexible structures

    NASA Technical Reports Server (NTRS)

    Banks, H. T.; Rebnord, D. A.

    1990-01-01

    Convergence and stability results for least squares inverse problems involving systems described by analytic semigroups are presented. The practical importance of these results is demonstrated by application to several examples from problems of estimation of material parameters in flexible structures using accelerometer data.

  13. Variational structure of inverse problems in wave propagation and vibration

    SciTech Connect

    Berryman, J.G.

    1995-03-01

    Practical algorithms for solving realistic inverse problems may often be viewed as problems in nonlinear programming with the data serving as constraints. Such problems are most easily analyzed when it is possible to segment the solution space into regions that are feasible (satisfying all the known constraints) and infeasible (violating some of the constraints). Then, if the feasible set is convex or at least compact, the solution to the problem will normally lie on the boundary of the feasible set. A nonlinear program may seek the solution by systematically exploring the boundary while satisfying progressively more constraints. Examples of inverse problems in wave propagation (traveltime tomography) and vibration (modal analysis) will be presented to illustrate how the variational structure of these problems may be used to create nonlinear programs using implicit variational constraints.

  14. From inverse problems in mathematical physiology to quantitative differential diagnoses.

    PubMed

    Zenker, Sven; Rubin, Jonathan; Clermont, Gilles

    2007-11-01

    The improved capacity to acquire quantitative data in a clinical setting has generally failed to improve outcomes in acutely ill patients, suggesting a need for advances in computer-supported data interpretation and decision making. In particular, the application of mathematical models of experimentally elucidated physiological mechanisms could augment the interpretation of quantitative, patient-specific information and help to better target therapy. Yet, such models are typically complex and nonlinear, a reality that often precludes the identification of unique parameters and states of the model that best represent available data. Hypothesizing that this non-uniqueness can convey useful information, we implemented a simplified simulation of a common differential diagnostic process (hypotension in an acute care setting), using a combination of a mathematical model of the cardiovascular system, a stochastic measurement model, and Bayesian inference techniques to quantify parameter and state uncertainty. The output of this procedure is a probability density function on the space of model parameters and initial conditions for a particular patient, based on prior population information together with patient-specific clinical observations. We show that multimodal posterior probability density functions arise naturally, even when unimodal and uninformative priors are used. The peaks of these densities correspond to clinically relevant differential diagnoses and can, in the simplified simulation setting, be constrained to a single diagnosis by assimilating additional observations from dynamical interventions (e.g., fluid challenge). We conclude that the ill-posedness of the inverse problem in quantitative physiology is not merely a technical obstacle, but rather reflects clinical reality and, when addressed adequately in the solution process, provides a novel link between mathematically described physiological knowledge and the clinical concept of differential diagnoses

  15. Inverse problems of ultrasound tomography in models with attenuation

    NASA Astrophysics Data System (ADS)

    Goncharsky, Alexander V.; Romanov, Sergey Y.

    2014-04-01

    We develop efficient methods for solving inverse problems of ultrasound tomography in models with attenuation. We treat the inverse problem as a coefficient inverse problem for unknown coordinate-dependent functions that characterize both the speed cross section and the coefficients of the wave equation describing attenuation in the diagnosed region. We derive exact formulas for the gradient of the residual functional in models with attenuation, and develop efficient algorithms for minimizing the gradient of the residual by solving the conjugate problem. These algorithms are easy to parallelize when implemented on supercomputers, allowing the computation time to be reduced by a factor of several hundred compared to a PC. The numerical analysis of model problems shows that it is possible to reconstruct not only the speed cross section, but also the properties of the attenuating medium. We investigate the choice of the initial approximation for iterative algorithms used to solve inverse problems. The algorithms considered are primarily meant for the development of ultrasound tomographs for differential diagnosis of breast cancer.

  16. PREFACE: Inverse Problems in Applied Sciences—towards breakthrough

    NASA Astrophysics Data System (ADS)

    Cheng, Jin; Iso, Yuusuke; Nakamura, Gen; Yamamoto, Masahiro

    2007-06-01

    These are the proceedings of the international conference `Inverse Problems in Applied Sciences—towards breakthrough' which was held at Hokkaido University, Sapporo, Japan on 3-7 July 2006 (http://coe.math.sci.hokudai.ac.jp/sympo/inverse/). There were 88 presentations and more than 100 participants, and we are proud to say that the conference was very successful. Nowadays, many new activities on inverse problems are flourishing at many centers of research around the world, and the conference has successfully gathered a world-wide variety of researchers. We believe that this volume contains not only main papers, but also conveys the general status of current research into inverse problems. This conference was the third biennial international conference on inverse problems, the core of which is the Pan-Pacific Asian area. The purpose of this series of conferences is to establish and develop constant international collaboration, especially among the Pan-Pacific Asian countries, and to lead the organization of activities concerning inverse problems centered in East Asia. The first conference was held at City University of Hong Kong in January 2002 and the second was held at Fudan University in June 2004. Following the preceding two successes, the third conference was organized in order to extend the scope of activities and build useful bridges to the next conference in Seoul in 2008. Therefore this third biennial conference was intended not only to establish collaboration and links between researchers in Asia and leading researchers worldwide in inverse problems but also to nurture interdisciplinary collaboration in theoretical fields such as mathematics, applied fields and evolving aspects of inverse problems. For these purposes, we organized tutorial lectures, serial lectures and a panel discussion as well as conference research presentations. This volume contains three lecture notes from the tutorial and serial lectures, and 22 papers. Especially at this

  17. Improved TV-CS Approaches for Inverse Scattering Problem

    PubMed Central

    Bevacqua, M. T.; Di Donato, L.

    2015-01-01

    Total Variation and Compressive Sensing (TV-CS) techniques represent a very attractive approach to inverse scattering problems. In fact, if the unknown is piecewise constant and so has a sparse gradient, TV-CS approaches allow us to achieve optimal reconstructions, reducing considerably the number of measurements and enforcing the sparsity on the gradient of the sought unknowns. In this paper, we introduce two different techniques based on TV-CS that exploit in a different manner the concept of gradient in order to improve the solution of the inverse scattering problems obtained by TV-CS approach. Numerical examples are addressed to show the effectiveness of the method. PMID:26495420

  18. History and evolution of methods for solving the inverse problem.

    PubMed

    van Oosterom, A

    1991-10-01

    This article serves as an introduction to the other articles in this issue devoted to the problem of the localization of neural generators. Elements of the theory of electric volume conduction are briefly introduced, as far as these apply to the interpretation of observed scalp potentials. First, some basic methods for display of the different aspects of the spatiotemporal information are described. Next, the most prominent source and volume conductor models that have been postulated for the involved forward problem are summarized. The problems of source identification and source localization, known as the inverse problem, are then formulated in terms of a parameter estimation procedure. The importance of introducing a priori information in the inverse problem, aimed at stabilizing (regularizing) the obtained solution, is emphasized. Methods for imposing such constraints are briefly outlined. PMID:1761703

  19. Solving probabilistic inverse problems rapidly with prior samples

    NASA Astrophysics Data System (ADS)

    Käufl, Paul; Valentine, Andrew; De Wit, Ralph; Trampert, Jeannot

    2016-03-01

    Owing to the increasing availability of computational resources, in recent years the probabilistic solution of non-linear, geophysical inverse problems by means of sampling methods has become increasingly feasible. Nevertheless, we still face situations in which a Monte Carlo approach is not practical. This is particularly true in cases where the evaluation of the forward problem is computationally intensive or where inversions have to be carried out repeatedly or in a timely manner, as in natural hazards monitoring tasks such as earthquake early warning. Here, we present an alternative to Monte Carlo sampling, in which inferences are entirely based on a set of prior samples-i.e. samples that have been obtained independent of a particular observed datum. This has the advantage that the computationally expensive sampling stage becomes separated from the inversion stage, and the set of prior samples-once obtained-can be reused for repeated evaluations of the inverse mapping without additional computational effort. This property is useful if the problem is such that repeated inversions of independent data have to be carried out. We formulate the inverse problem in a Bayesian framework and present a practical way to make posterior inferences based on a set of prior samples. We compare the prior sampling based approach to a Markov Chain Monte Carlo approach that samples from the posterior probability distribution. We show results for both a toy example, and a realistic seismological source parameter estimation problem. We find that the posterior uncertainty estimates obtained based on prior sampling can be considered conservative estimates of the uncertainties obtained by directly sampling from the posterior distribution.

  20. Solving probabilistic inverse problems rapidly with prior samples

    NASA Astrophysics Data System (ADS)

    Käufl, Paul; Valentine, Andrew P.; de Wit, Ralph W.; Trampert, Jeannot

    2016-06-01

    Owing to the increasing availability of computational resources, in recent years the probabilistic solution of non-linear, geophysical inverse problems by means of sampling methods has become increasingly feasible. Nevertheless, we still face situations in which a Monte Carlo approach is not practical. This is particularly true in cases where the evaluation of the forward problem is computationally intensive or where inversions have to be carried out repeatedly or in a timely manner, as in natural hazards monitoring tasks such as earthquake early warning. Here, we present an alternative to Monte Carlo sampling, in which inferences are entirely based on a set of prior samples-that is, samples that have been obtained independent of a particular observed datum. This has the advantage that the computationally expensive sampling stage becomes separated from the inversion stage, and the set of prior samples-once obtained-can be reused for repeated evaluations of the inverse mapping without additional computational effort. This property is useful if the problem is such that repeated inversions of independent data have to be carried out. We formulate the inverse problem in a Bayesian framework and present a practical way to make posterior inferences based on a set of prior samples. We compare the prior sampling based approach to a Markov Chain Monte Carlo approach that samples from the posterior probability distribution. We show results for both a toy example, and a realistic seismological source parameter estimation problem. We find that the posterior uncertainty estimates obtained based on prior sampling can be considered conservative estimates of the uncertainties obtained by directly sampling from the posterior distribution.

  1. A tutorial on inverse problems for anomalous diffusion processes

    NASA Astrophysics Data System (ADS)

    Jin, Bangti; Rundell, William

    2015-03-01

    Over the last two decades, anomalous diffusion processes in which the mean squares variance grows slower or faster than that in a Gaussian process have found many applications. At a macroscopic level, these processes are adequately described by fractional differential equations, which involves fractional derivatives in time or/and space. The fractional derivatives describe either history mechanism or long range interactions of particle motions at a microscopic level. The new physics can change dramatically the behavior of the forward problems. For example, the solution operator of the time fractional diffusion diffusion equation has only limited smoothing property, whereas the solution for the space fractional diffusion equation may contain weak singularity. Naturally one expects that the new physics will impact related inverse problems in terms of uniqueness, stability, and degree of ill-posedness. The last aspect is especially important from a practical point of view, i.e., stably reconstructing the quantities of interest. In this paper, we employ a formal analytic and numerical way, especially the two-parameter Mittag-Leffler function and singular value decomposition, to examine the degree of ill-posedness of several ‘classical’ inverse problems for fractional differential equations involving a Djrbashian-Caputo fractional derivative in either time or space, which represent the fractional analogues of that for classical integral order differential equations. We discuss four inverse problems, i.e., backward fractional diffusion, sideways problem, inverse source problem and inverse potential problem for time fractional diffusion, and inverse Sturm-Liouville problem, Cauchy problem, backward fractional diffusion and sideways problem for space fractional diffusion. It is found that contrary to the wide belief, the influence of anomalous diffusion on the degree of ill-posedness is not definitive: it can either significantly improve or worsen the conditioning of

  2. A time domain sampling method for inverse acoustic scattering problems

    NASA Astrophysics Data System (ADS)

    Guo, Yukun; Hömberg, Dietmar; Hu, Guanghui; Li, Jingzhi; Liu, Hongyu

    2016-06-01

    This work concerns the inverse scattering problems of imaging unknown/inaccessible scatterers by transient acoustic near-field measurements. Based on the analysis of the migration method, we propose efficient and effective sampling schemes for imaging small and extended scatterers from knowledge of time-dependent scattered data due to incident impulsive point sources. Though the inverse scattering problems are known to be nonlinear and ill-posed, the proposed imaging algorithms are totally "direct" involving only integral calculations on the measurement surface. Theoretical justifications are presented and numerical experiments are conducted to demonstrate the effectiveness and robustness of our methods. In particular, the proposed static imaging functionals enhance the performance of the total focusing method (TFM) and the dynamic imaging functionals show analogous behavior to the time reversal inversion but without solving time-dependent wave equations.

  3. Solving the Inverse-Square Problem with Complex Variables

    ERIC Educational Resources Information Center

    Gauthier, N.

    2005-01-01

    The equation of motion for a mass that moves under the influence of a central, inverse-square force is formulated and solved as a problem in complex variables. To find the solution, the constancy of angular momentum is first established using complex variables. Next, the complex position coordinate and complex velocity of the particle are assumed…

  4. Toward precise solution of one-dimensional velocity inverse problems

    SciTech Connect

    Gray, S.; Hagin, F.

    1980-01-01

    A family of one-dimensional inverse problems are considered with the goal of reconstructing velocity profiles to reasonably high accuracy. The travel-time variable change is used together with an iteration scheme to produce an effective algorithm for computation. Under modest assumptions the scheme is shown to be convergent.

  5. Computational methods for inverse problems in geophysics: inversion of travel time observations

    USGS Publications Warehouse

    Pereyra, V.; Keller, H.B.; Lee, W.H.K.

    1980-01-01

    General ways of solving various inverse problems are studied for given travel time observations between sources and receivers. These problems are separated into three components: (a) the representation of the unknown quantities appearing in the model; (b) the nonlinear least-squares problem; (c) the direct, two-point ray-tracing problem used to compute travel time once the model parameters are given. Novel software is described for (b) and (c), and some ideas given on (a). Numerical results obtained with artificial data and an implementation of the algorithm are also presented. ?? 1980.

  6. Solution of 3D inverse scattering problems by combined inverse equivalent current and finite element methods

    SciTech Connect

    Kılıç, Emre Eibert, Thomas F.

    2015-05-01

    An approach combining boundary integral and finite element methods is introduced for the solution of three-dimensional inverse electromagnetic medium scattering problems. Based on the equivalence principle, unknown equivalent electric and magnetic surface current densities on a closed surface are utilized to decompose the inverse medium problem into two parts: a linear radiation problem and a nonlinear cavity problem. The first problem is formulated by a boundary integral equation, the computational burden of which is reduced by employing the multilevel fast multipole method (MLFMM). Reconstructed Cauchy data on the surface allows the utilization of the Lorentz reciprocity and the Poynting's theorems. Exploiting these theorems, the noise level and an initial guess are estimated for the cavity problem. Moreover, it is possible to determine whether the material is lossy or not. In the second problem, the estimated surface currents form inhomogeneous boundary conditions of the cavity problem. The cavity problem is formulated by the finite element technique and solved iteratively by the Gauss–Newton method to reconstruct the properties of the object. Regularization for both the first and the second problems is achieved by a Krylov subspace method. The proposed method is tested against both synthetic and experimental data and promising reconstruction results are obtained.

  7. Spectral collocation for multiparameter eigenvalue problems arising from separable boundary value problems

    NASA Astrophysics Data System (ADS)

    Plestenjak, Bor; Gheorghiu, Călin I.; Hochstenbach, Michiel E.

    2015-10-01

    In numerous science and engineering applications a partial differential equation has to be solved on some fairly regular domain that allows the use of the method of separation of variables. In several orthogonal coordinate systems separation of variables applied to the Helmholtz, Laplace, or Schrödinger equation leads to a multiparameter eigenvalue problem (MEP); important cases include Mathieu's system, Lamé's system, and a system of spheroidal wave functions. Although multiparameter approaches are exploited occasionally to solve such equations numerically, MEPs remain less well known, and the variety of available numerical methods is not wide. The classical approach of discretizing the equations using standard finite differences leads to algebraic MEPs with large matrices, which are difficult to solve efficiently. The aim of this paper is to change this perspective. We show that by combining spectral collocation methods and new efficient numerical methods for algebraic MEPs it is possible to solve such problems both very efficiently and accurately. We improve on several previous results available in the literature, and also present a MATLAB toolbox for solving a wide range of problems.

  8. PREFACE: First International Congress of the International Association of Inverse Problems (IPIA): Applied Inverse Problems 2007: Theoretical and Computational Aspects

    NASA Astrophysics Data System (ADS)

    Uhlmann, Gunther

    2008-07-01

    This volume represents the proceedings of the fourth Applied Inverse Problems (AIP) international conference and the first congress of the Inverse Problems International Association (IPIA) which was held in Vancouver, Canada, June 25 29, 2007. The organizing committee was formed by Uri Ascher, University of British Columbia, Richard Froese, University of British Columbia, Gary Margrave, University of Calgary, and Gunther Uhlmann, University of Washington, chair. The conference was part of the activities of the Pacific Institute of Mathematical Sciences (PIMS) Collaborative Research Group on inverse problems (http://www.pims.math.ca/scientific/collaborative-research-groups/past-crgs). This event was also supported by grants from NSF and MITACS. Inverse Problems (IP) are problems where causes for a desired or an observed effect are to be determined. They lie at the heart of scientific inquiry and technological development. The enormous increase in computing power and the development of powerful algorithms have made it possible to apply the techniques of IP to real-world problems of growing complexity. Applications include a number of medical as well as other imaging techniques, location of oil and mineral deposits in the earth's substructure, creation of astrophysical images from telescope data, finding cracks and interfaces within materials, shape optimization, model identification in growth processes and, more recently, modelling in the life sciences. The series of Applied Inverse Problems (AIP) Conferences aims to provide a primary international forum for academic and industrial researchers working on all aspects of inverse problems, such as mathematical modelling, functional analytic methods, computational approaches, numerical algorithms etc. The steering committee of the AIP conferences consists of Heinz Engl (Johannes Kepler Universität, Austria), Joyce McLaughlin (RPI, USA), William Rundell (Texas A&M, USA), Erkki Somersalo (Helsinki University of Technology

  9. PREFACE: First International Congress of the International Association of Inverse Problems (IPIA): Applied Inverse Problems 2007: Theoretical and Computational Aspects

    NASA Astrophysics Data System (ADS)

    Uhlmann, Gunther

    2008-07-01

    This volume represents the proceedings of the fourth Applied Inverse Problems (AIP) international conference and the first congress of the Inverse Problems International Association (IPIA) which was held in Vancouver, Canada, June 25 29, 2007. The organizing committee was formed by Uri Ascher, University of British Columbia, Richard Froese, University of British Columbia, Gary Margrave, University of Calgary, and Gunther Uhlmann, University of Washington, chair. The conference was part of the activities of the Pacific Institute of Mathematical Sciences (PIMS) Collaborative Research Group on inverse problems (http://www.pims.math.ca/scientific/collaborative-research-groups/past-crgs). This event was also supported by grants from NSF and MITACS. Inverse Problems (IP) are problems where causes for a desired or an observed effect are to be determined. They lie at the heart of scientific inquiry and technological development. The enormous increase in computing power and the development of powerful algorithms have made it possible to apply the techniques of IP to real-world problems of growing complexity. Applications include a number of medical as well as other imaging techniques, location of oil and mineral deposits in the earth's substructure, creation of astrophysical images from telescope data, finding cracks and interfaces within materials, shape optimization, model identification in growth processes and, more recently, modelling in the life sciences. The series of Applied Inverse Problems (AIP) Conferences aims to provide a primary international forum for academic and industrial researchers working on all aspects of inverse problems, such as mathematical modelling, functional analytic methods, computational approaches, numerical algorithms etc. The steering committee of the AIP conferences consists of Heinz Engl (Johannes Kepler Universität, Austria), Joyce McLaughlin (RPI, USA), William Rundell (Texas A&M, USA), Erkki Somersalo (Helsinki University of Technology

  10. Application of inverse heat conduction problem on temperature measurement

    NASA Astrophysics Data System (ADS)

    Zhang, X.; Zhou, G.; Dong, B.; Li, Q.; Liu, L. Q.

    2013-09-01

    For regenerative cooling devices, such as G-M refrigerator, pulse tube cooler or thermoacoustic cooler, the gas oscillating bring about temperature fluctuations inevitably, which is harmful in many applications requiring high stable temperatures. To find out the oscillating mechanism of the cooling temperature and improve the temperature stability of cooler, the inner temperature of the cold head has to be measured. However, it is difficult to measure the inner oscillating temperature of the cold head directly because the invasive temperature detectors may disturb the oscillating flow. Fortunately, the outer surface temperature of the cold head can be measured accurately by invasive temperature measurement techniques. In this paper, a mathematical model of inverse heat conduction problem is presented to identify the inner surface oscillating temperature of cold head according to the measured temperature of the outer surface in a GM cryocooler. Inverse heat conduction problem will be solved using control volume approach. Outer surface oscillating temperature could be used as input conditions of inverse problem and the inner surface oscillating temperature of cold head can be inversely obtained. A simple uncertainty analysis of the oscillating temperature measurement also will be provided.

  11. Inverse AVO problem for a stack of layers

    NASA Astrophysics Data System (ADS)

    Malovichko, Liliya

    2015-09-01

    The problem of estimating thin layered model parameters by amplitude variation with offset (AVO) inversion has been studied. The motivation was resolving of the thin layers in inverted prestack seismic data as it contains more information on elastic properties of the subsurface than poststack seismic data. In this paper, an algorithm for solving the prestack inverse AVO problem in the case of multilayered media is derived. This algorithm is based on iterative corrections to the parameters of the initial model which tend to minimise the misfits between observed and synthetic seismograms. The synthetic seismograms are calculated using the reflection-transmission (RT)-matrices method, assuming a plane-wave with respect to the source position. A regularised Gauss-type algorithm for the inversion of prestack seismic data has been used. A differential seismogram computation algorithm to characterise the sensitivity of the seismic signal to the variations of a model parameter was used. The derived solution of the inverse problem is constructed in the time domain. This gives a slight advantage because it allows for visual control of the solution process. One can monitor the amplitude reduction of the data residual (difference between observed and synthetic seismograms) during the iteration process. Numerical examples show the accuracy and efficiency of the method.

  12. SIAM conference on inverse problems: Geophysical applications. Final technical report

    SciTech Connect

    1995-12-31

    This conference was the second in a series devoted to a particular area of inverse problems. The theme of this series is to discuss problems of major scientific importance in a specific area from a mathematical perspective. The theme of this symposium was geophysical applications. In putting together the program we tried to include a wide range of mathematical scientists and to interpret geophysics in as broad a sense as possible. Our speaker came from industry, government laboratories, and diverse departments in academia. We managed to attract a geographically diverse audience with participation from five continents. There were talks devoted to seismology, hydrology, determination of the earth`s interior on a global scale as well as oceanographic and atmospheric inverse problems.

  13. A reduced basis Landweber method for nonlinear inverse problems

    NASA Astrophysics Data System (ADS)

    Garmatter, Dominik; Haasdonk, Bernard; Harrach, Bastian

    2016-03-01

    We consider parameter identification problems in parametrized partial differential equations (PDEs). These lead to nonlinear ill-posed inverse problems. One way of solving them is using iterative regularization methods, which typically require numerous amounts of forward solutions during the solution process. In this article we consider the nonlinear Landweber method and couple it with the reduced basis method as a model order reduction technique in order to reduce the overall computational time. In particular, we consider PDEs with a high-dimensional parameter space, which are known to pose difficulties in the context of reduced basis methods. We present a new method that is able to handle such high-dimensional parameter spaces by combining the nonlinear Landweber method with adaptive online reduced basis updates. It is then applied to the inverse problem of reconstructing the conductivity in the stationary heat equation.

  14. Effects of geometric uncertainty on the inverse EEG problem

    NASA Astrophysics Data System (ADS)

    Weinstein, David M.; Johnson, Christopher R.

    1997-12-01

    A standard method for noninvasively computing neurocortical potentials from potentials measured on the scalp surface is to solve the problem on a generalized geometry and map the results back to the true model. This solution to the inverse EEG problem has been employed using spherical and, more recently, generic cranial models as templates. In the case of the most complex spherical models, the patient's skin, bone, cerebrospinal fluid, gray matter and white matter surfaces are mapped onto concentric spheres. The simplicity of the spherical domain allows for an analytic solution to the surface mapping inverse problem; however, the inaccuracy of such a solution challenges its clinical value. Similarly, solving the problem on a predefined generic model also holds computational allure--the generic model can be hand-picked to reduce the ill-conditioning of the problem. However, we suggest that such results from generic models are still not sufficiently accurate to be of general clinical use. In our paper, we evaluate the impact of varying both model accuracy and model complexity on the inverse cortical mapping. Small modeling perturbations (as might be introduced from noisy or under-sampled data) are shown to have large and detrimental effects on the quality of the solution.

  15. An Inverse Problem Approach for Elasticity Imaging through Vibroacoustics

    PubMed Central

    Aguilo, Miguel A.; Brigham, J. C.; Aquino, W.; Fatemi, M.

    2011-01-01

    A new methodology for estimating the spatial distribution of elastic moduli using the steady-state dynamic response of solids immersed in fluids is presented. The technique relies on the ensuing acoustic pressure field from a remotely excited solid to inversely estimate the spatial distribution of Young’s modulus. This work proposes the use of Gaussian radial basis functions (GRBF) to represent the spatial variation of elastic moduli. GRBF are shown to possess the advantage of representing smooth functions with quasi-compact support, and can efficiently represent elastic moduli distributions such as those that occur in soft biological tissue in the presence of tumors. The direct problem consists of a coupled acoustic-structure interaction boundary value problem solved in the frequency domain using the finite element method. The inverse problem is cast as an optimization problem in which the objective function is defined as a measure of discrepancy between an experimentally measured response and a finite element representation of the system. Non-gradient based optimization algorithms in combination with a divide and conquer strategy are used to solve the resulting optimization problem. The feasibility of the proposed approach is demonstrated through a series of numerical and a physical experiment. For comparison purposes, the surface velocity response was also used for the inverse characterization as the measured response in place of the acoustic pressure. PMID:20335092

  16. A non-local free boundary problem arising in a theory of financial bubbles

    PubMed Central

    Berestycki, Henri; Monneau, Regis; Scheinkman, José A.

    2014-01-01

    We consider an evolution non-local free boundary problem that arises in the modelling of speculative bubbles. The solution of the model is the speculative component in the price of an asset. In the framework of viscosity solutions, we show the existence and uniqueness of the solution. We also show that the solution is convex in space, and establish several monotonicity properties of the solution and of the free boundary with respect to parameters of the problem. To study the free boundary, we use, in particular, the fact that the odd part of the solution solves a more standard obstacle problem. We show that the free boundary is and describe the asymptotics of the free boundary as c, the cost of transacting the asset, goes to zero. PMID:25288815

  17. Forward and inverse problems in fundamental and applied magnetohydrodynamics

    NASA Astrophysics Data System (ADS)

    Giesecke, Andre; Stefani, Frank; Wondrak, Thomas; Xu, Mingtian

    2013-03-01

    This minireview summarizes the recent efforts to solve forward and inverse problems as they occur in different branches of fundamental and applied magnetohydrodynamics. For the forward problem, the main focus is on the numerical treatment of induction processes, including self-excitation of magnetic fields in non-spherical domains and/or under the influence of non-homogeneous material parameters. As an important application of the developed numerical schemes, the functioning of the von-Kármán-sodium (VKS) dynamo experiment is shown to depend crucially on the presence of soft-iron impellers. As for the inverse problem, the main focus is on the mathematical background and some initial practical applications of contactless inductive flow tomography (CIFT), in which flow induced magnetic field perturbations are utilized to reconstruct the velocity field. The promises of CIFT for flow field monitoring in the continuous casting of steel are substantiated by results obtained at a test rig with a low-melting liquid metal. While CIFT is presently restricted to flows with low magnetic Reynolds numbers, some selected problems from non-linear inverse dynamo theory, with possible applications to geo- and astrophysics, are also discussed.

  18. Two hybrid regularization frameworks for solving the electrocardiography inverse problem

    NASA Astrophysics Data System (ADS)

    Jiang, Mingfeng; Xia, Ling; Shou, Guofa; Liu, Feng; Crozier, Stuart

    2008-09-01

    In this paper, two hybrid regularization frameworks, LSQR-Tik and Tik-LSQR, which integrate the properties of the direct regularization method (Tikhonov) and the iterative regularization method (LSQR), have been proposed and investigated for solving ECG inverse problems. The LSQR-Tik method is based on the Lanczos process, which yields a sequence of small bidiagonal systems to approximate the original ill-posed problem and then the Tikhonov regularization method is applied to stabilize the projected problem. The Tik-LSQR method is formulated as an iterative LSQR inverse, augmented with a Tikhonov-like prior information term. The performances of these two hybrid methods are evaluated using a realistic heart-torso model simulation protocol, in which the heart surface source method is employed to calculate the simulated epicardial potentials (EPs) from the action potentials (APs), and then the acquired EPs are used to calculate simulated body surface potentials (BSPs). The results show that the regularized solutions obtained by the LSQR-Tik method are approximate to those of the Tikhonov method, the computational cost of the LSQR-Tik method, however, is much less than that of the Tikhonov method. Moreover, the Tik-LSQR scheme can reconstruct the epcicardial potential distribution more accurately, specifically for the BSPs with large noisy cases. This investigation suggests that hybrid regularization methods may be more effective than separate regularization approaches for ECG inverse problems.

  19. Stabilizing inverse problems by internal data. II: non-local internal data and generic linearized uniqueness

    NASA Astrophysics Data System (ADS)

    Kuchment, Peter; Steinhauer, Dustin

    2015-12-01

    In the previous paper (Kuchment and Steinhauer in Inverse Probl 28(8):084007, 2012), the authors introduced a simple procedure that allows one to detect whether and explain why internal information arising in several novel coupled physics (hybrid) imaging modalities could turn extremely unstable techniques, such as optical tomography or electrical impedance tomography, into stable, good-resolution procedures. It was shown that in all cases of interest, the Fréchet derivative of the forward mapping is a pseudo-differential operator with an explicitly computable principal symbol. If one can set up the imaging procedure in such a way that the symbol is elliptic, this would indicate that the problem was stabilized. In the cases when the symbol is not elliptic, the technique suggests how to change the procedure (e.g., by adding extra measurements) to achieve ellipticity. In this article, we consider the situation arising in acousto-optical tomography (also called ultrasound modulated optical tomography), where the internal data available involves the Green's function, and thus depends globally on the unknown parameter(s) of the equation and its solution. It is shown that the technique of (Kuchment and Steinhauer in Inverse Probl 28(8):084007, 2012) can be successfully adopted to this situation as well. A significant part of the article is devoted to results on generic uniqueness for the linearized problem in a variety of situations, including those arising in acousto-electric and quantitative photoacoustic tomography.

  20. Inverse problems for homogeneous transport equations: II. The multidimensional case

    NASA Astrophysics Data System (ADS)

    Bal, Guillaume

    2000-08-01

    A companion paper by Bal (Bal G 2000 Inverse Problems 16 997) and this paper are parts I and II of a series dealing with the reconstruction from boundary measurements of the scattering operator of homogeneous linear transport equations. This part II deals with the case of convex bounded domains in dimensions higher than one. We distinguish the analysis of smooth boundaries from that of boundaries with discontinuities such as corners. We propose a reconstruction in the case of degenerate symmetric scattering operators and show the well-posedness of the inverse problem. The proof of well-posedness is based on a decomposition of angular moments of the transport solution into unbounded and bounded components. This decomposition allows us to show the linear independence of a sufficiently large number of angular moments of the transport solution that are used to construct an invertible system for the scattering coefficients to be reconstructed.

  1. Inverse problems in the modeling of vibrations of flexible beams

    NASA Technical Reports Server (NTRS)

    Banks, H. T.; Powers, R. K.; Rosen, I. G.

    1987-01-01

    The formulation and solution of inverse problems for the estimation of parameters which describe damping and other dynamic properties in distributed models for the vibration of flexible structures is considered. Motivated by a slewing beam experiment, the identification of a nonlinear velocity dependent term which models air drag damping in the Euler-Bernoulli equation is investigated. Galerkin techniques are used to generate finite dimensional approximations. Convergence estimates and numerical results are given. The modeling of, and related inverse problems for the dynamics of a high pressure hose line feeding a gas thruster actuator at the tip of a cantilevered beam are then considered. Approximation and convergence are discussed and numerical results involving experimental data are presented.

  2. Lycurgus Cup: inverse problem using photographs for characterization of matter.

    PubMed

    Barchiesi, Dominique

    2015-08-01

    Photographs of the Lycurgus Cup with a source light inside and outside exhibit purple and green colors, respectively (dichroism). A model relying on the scattering of light to colors in the photographs is proposed and used within an inverse problem algorithm, to deduce radius and composition of metallic particles, and the refractive index of the surrounding glass medium. The inverse problem algorithm is based on a hybridization of particle swarm optimization and of the simulated annealing methods. The results are compared to experimental measurements on a small sample of glass. The linear laws that are deduced from sets of possible parameters producing the same color in the photographs help simplify the understanding of phenomena. The proportion of silver to gold in nanoparticles is found to be in agreement, but a large proportion of copper is also found. The retrieved refractive index of the surrounding glass is close to 2. PMID:26367298

  3. Solving the inverse problem of noise-driven dynamic networks.

    PubMed

    Zhang, Zhaoyang; Zheng, Zhigang; Niu, Haijing; Mi, Yuanyuan; Wu, Si; Hu, Gang

    2015-01-01

    Nowadays, massive amounts of data are available for analysis in natural and social systems and the tasks to depict system structures from the data, i.e., the inverse problems, become one of the central issues in wide interdisciplinary fields. In this paper, we study the inverse problem of dynamic complex networks driven by white noise. A simple and universal inference formula of double correlation matrices and noise-decorrelation (DCMND) method is derived analytically, and numerical simulations confirm that the DCMND method can accurately depict both network structures and noise correlations by using available output data only. This inference performance has never been regarded possible by theoretical derivation, numerical computation, and experimental design. PMID:25679664

  4. Changes in habitat of fish populations: An inverse problem.

    PubMed

    Levere, Kimberly M

    2016-08-01

    Mathematical modelling applies to a wide variety of application areas, and is an active area of research in many disciplines. It is often the case that accurate depiction of real-world phenomena require increasingly complex models. Unfortunately, this increased complexity in a model causes great difficulty when seeking solutions. What is more, developing a model with known parameters that produces results consistent with observed behaviors may prove to be a difficult or even impossible task. These difficulties have brought about an interest in inverse problems. In this paper we utilize a collage-based approach to solve an inverse problem for a model for the migration of three fish species through floodplain waters. A derivation of the mathematical model is presented and a generalized collage method is discussed and applied to this model to recover diffusion parameters. Theoretical and numerical particulars are discussed and results are presented. PMID:27245383

  5. Solving the inverse problem of noise-driven dynamic networks

    NASA Astrophysics Data System (ADS)

    Zhang, Zhaoyang; Zheng, Zhigang; Niu, Haijing; Mi, Yuanyuan; Wu, Si; Hu, Gang

    2015-01-01

    Nowadays, massive amounts of data are available for analysis in natural and social systems and the tasks to depict system structures from the data, i.e., the inverse problems, become one of the central issues in wide interdisciplinary fields. In this paper, we study the inverse problem of dynamic complex networks driven by white noise. A simple and universal inference formula of double correlation matrices and noise-decorrelation (DCMND) method is derived analytically, and numerical simulations confirm that the DCMND method can accurately depict both network structures and noise correlations by using available output data only. This inference performance has never been regarded possible by theoretical derivation, numerical computation, and experimental design.

  6. Efficient algorithms for linear dynamic inverse problems with known motion

    NASA Astrophysics Data System (ADS)

    Hahn, B. N.

    2014-03-01

    An inverse problem is called dynamic if the object changes during the data acquisition process. This occurs e.g. in medical applications when fast moving organs like the lungs or the heart are imaged. Most regularization methods are based on the assumption that the object is static during the measuring procedure. Hence, their application in the dynamic case often leads to serious motion artefacts in the reconstruction. Therefore, an algorithm has to take into account the temporal changes of the investigated object. In this paper, a reconstruction method that compensates for the motion of the object is derived for dynamic linear inverse problems. The algorithm is validated at numerical examples from computerized tomography.

  7. An inverse finance problem for estimation of the volatility

    NASA Astrophysics Data System (ADS)

    Neisy, A.; Salmani, K.

    2013-01-01

    Black-Scholes model, as a base model for pricing in derivatives markets has some deficiencies, such as ignoring market jumps, and considering market volatility as a constant factor. In this article, we introduce a pricing model for European-Options under jump-diffusion underlying asset. Then, using some appropriate numerical methods we try to solve this model with integral term, and terms including derivative. Finally, considering volatility as an unknown parameter, we try to estimate it by using our proposed model. For the purpose of estimating volatility, in this article, we utilize inverse problem, in which inverse problem model is first defined, and then volatility is estimated using minimization function with Tikhonov regularization.

  8. Diffuse interface methods for inverse problems: case study for an elliptic Cauchy problem

    NASA Astrophysics Data System (ADS)

    Burger, Martin; Løseth Elvetun, Ole; Schlottbom, Matthias

    2015-12-01

    Many inverse problems have to deal with complex, evolving and often not exactly known geometries, e.g. as domains of forward problems modeled by partial differential equations. This makes it desirable to use methods which are robust with respect to perturbed or not well resolved domains, and which allow for efficient discretizations not resolving any fine detail of those geometries. For forward problems in partial differential equations methods based on diffuse interface representations have gained strong attention in the last years, but so far they have not been considered systematically for inverse problems. In this work we introduce a diffuse domain method as a tool for the solution of variational inverse problems. As a particular example we study ECG inversion in further detail. ECG inversion is a linear inverse source problem with boundary measurements governed by an anisotropic diffusion equation, which naturally cries for solutions under changing geometries, namely the beating heart. We formulate a regularization strategy using Tikhonov regularization and, using standard source conditions, we prove convergence rates. A special property of our approach is that not only operator perturbations are introduced by the diffuse domain method, but more important we have to deal with topologies which depend on a parameter \\varepsilon in the diffuse domain method, i.e. we have to deal with \\varepsilon -dependent forward operators and \\varepsilon -dependent norms. In particular the appropriate function spaces for the unknown and the data depend on \\varepsilon . This prevents the application of some standard convergence techniques for inverse problems, in particular interpreting the perturbations as data errors in the original problem does not yield suitable results. We consequently develop a novel approach based on saddle-point problems. The numerical solution of the problem is discussed as well and results for several computational experiments are reported. In

  9. Eddy-current NDE inverse problem with sparse grid algorithm

    NASA Astrophysics Data System (ADS)

    Zhou, Liming; Sabbagh, Harold A.; Sabbagh, Elias H.; Murphy, R. Kim; Bernacchi, William; Aldrin, John C.; Forsyth, David; Lindgren, Eric

    2016-02-01

    In model-based inverse problems, the unknown parameters (such as length, width, depth) need to be estimated. When the unknown parameters are few, the conventional mathematical methods are suitable. But the increasing number of unknown parameters will make the computation become heavy. To reduce the burden of computation, the sparse grid algorithm was used in our work. As a result, we obtain a powerful interpolation method that requires significantly fewer support nodes than conventional interpolation on a full grid.

  10. Some special cases of the electromagnetic inverse problem

    NASA Technical Reports Server (NTRS)

    Weston, V. H.

    1972-01-01

    A review of exact techniques for determining the surface of a three-dimensional perfectly conducting body is given, followed by some new results on the uniqueness question concerning the number of measurements that may be required to explicitly determine the surface of the body. It is then shown that the inhomogeneous but spherically symmetric dielectric electromagnetic case is reducible to a scalar inverse problem that can be treated by known techniques.

  11. Combined approach to the inverse protein folding problem. Final report

    SciTech Connect

    Ruben A. Abagyan

    2000-06-01

    The main scientific contribution of the project ''Combined approach to the inverse protein folding problem'' submitted in 1996 and funded by the Department of Energy in 1997 is the formulation and development of the idea of the multilink recognition method for identification of functional and structural homologues of newly discovered genes. This idea became very popular after they first announced it and used it in prediction of the threading targets for the CASP2 competition (Critical Assessment of Structure Prediction).

  12. Inverse problems of NEO photometry: Imaging the NEO population

    NASA Astrophysics Data System (ADS)

    Kaasalainen1, Mikko; Durech, Josef

    2007-05-01

    Photometry is the main source of information on NEOs (and other asteroids) en masse. Surveys such as Pan-STARRS and LSST will produce colossal photometric databases that can readily be used for mapping the physical characteristics of the whole asteroid population. These datasets are efficiently enriched by any additional dense photometric or other observations. Due to their quickly changing geometries with respect to the Earth, NEOs are the subpopulation that can be mapped the fastest. I review the state of the art in the construction of physical asteroid models from sparse and/or dense photometric data (that can also be combined with other data modes). The models describe the shapes, spin states, scattering properties and surface structure of the targets, and are the solutions of inverse problems necessarily involving comprehensive mathematical analysis. I sum up what we can and cannot get from photometric data, and how all this is done in practice. I also discuss the new freely available software package for solving photometric inverse problems (soon to be released). The analysis of photometric datasets will very soon become an automated industry, resulting in tens of thousands of asteroid models, a large portion of them NEOs. The computational effort in this is considerable in both computer and human time, which means that a large portion of the targets is likely to be analyzed only once. This, again, means that we have to have a good understanding of the reliability of our models, and this is impossible without a thorough understanding of the mathematical nature of the inverse problem(s) involved. Very important concepts are the uniqueness and stability of the solution, the parameter spaces, the so-called inverse crimes in simulations and error prediction, and the domination of systematic errors over random ones.

  13. Asynchronous global optimization techniques for medium and large inversion problems

    SciTech Connect

    Pereyra, V.; Koshy, M.; Meza, J.C.

    1995-04-01

    We discuss global optimization procedures adequate for seismic inversion problems. We explain how to save function evaluations (which may involve large scale ray tracing or other expensive operations) by creating a data base of information on what parts of parameter space have already been inspected. It is also shown how a correct parallel implementation using PVM speeds up the process almost linearly with respect to the number of processors, provided that the function evaluations are expensive enough to offset the communication overhead.

  14. A new efficient recursive technique for solving singular boundary value problems arising in various physical models

    NASA Astrophysics Data System (ADS)

    Roul, Pradip

    2016-04-01

    The paper deals with a numerical technique for solving nonlinear singular boundary value problems arising in various physical models. First, we convert the original problem to an equivalent integral equation to surmount the singularity and employ afterward the boundary condition to compute the undetermined coefficient. Finally, the integral equation without undetermined coefficient is treated using homotopy perturbation method. The present method is implemented on three physical model examples: i) thermal explosions; ii) steady-state oxygen diffusion in a spherical shell; iii) the equilibrium of the isothermal gas sphere. The results obtained by the present method are compared with that obtained using finite-difference method, B-spline method and a numerical technique based on the direct integration method, and comparison reveals that the proposed method with few solution components produces similar results and the method is computationally efficient than others.

  15. Introduction to the 30th volume of Inverse Problems

    NASA Astrophysics Data System (ADS)

    Louis, Alfred K.

    2014-01-01

    The field of inverse problems is a fast-developing domain of research originating from the practical demands of finding the cause when a result is observed. The woodpecker, searching for insects, is probing a tree using sound waves: the information searched for is whether there is an insect or not, hence a 0-1 decision. When the result has to contain more information, ad hoc solutions are not at hand and more sophisticated methods have to be developed. Right from its first appearance, the field of inverse problems has been characterized by an interdisciplinary nature: the interpretation of measured data, reinforced by mathematical models serving the analyzing questions of observability, stability and resolution, developing efficient, stable and accurate algorithms to gain as much information as possible from the input and to feedback to the questions of optimal measurement configuration. As is typical for a new area of research, facets of it are separated and studied independently. Hence, fields such as the theory of inverse scattering, tomography in general and regularization methods have developed. However, all aspects have to be reassembled to arrive at the best possible solution to the problem at hand. This development is reflected by the first and still leading journal in the field, Inverse Problems. Founded by pioneers Roy Pike from London and Pierre Sabatier from Montpellier, who enjoyably describes the journal's nascence in his book Rêves et Combats d'un Enseignant-Chercheur, Retour Inverse [1], the journal has developed successfully over the last few decades. Neither the Editors-in-Chief, formerly called Honorary Editors, nor the board or authors could have set the path to success alone. Their fruitful interplay, complemented by the efficient and highly competent publishing team at IOP Publishing, has been fundamental. As such it is my honor and pleasure to follow my renowned colleagues Pierre Sabatier, Mario Bertero, Frank Natterer, Alberto Grünbaum and

  16. Material derivatives of boundary integral operators in electromagnetism and application to inverse scattering problems

    NASA Astrophysics Data System (ADS)

    Ivanyshyn Yaman, Olha; Le Louër, Frédérique

    2016-09-01

    This paper deals with the material derivative analysis of the boundary integral operators arising from the scattering theory of time-harmonic electromagnetic waves and its application to inverse problems. We present new results using the Piola transform of the boundary parametrisation to transport the integral operators on a fixed reference boundary. The transported integral operators are infinitely differentiable with respect to the parametrisations and simplified expressions of the material derivatives are obtained. Using these results, we extend a nonlinear integral equations approach developed for solving acoustic inverse obstacle scattering problems to electromagnetism. The inverse problem is formulated as a pair of nonlinear and ill-posed integral equations for the unknown boundary representing the boundary condition and the measurements, for which the iteratively regularized Gauss-Newton method can be applied. The algorithm has the interesting feature that it avoids the numerous numerical solution of boundary value problems at each iteration step. Numerical experiments are presented in the special case of star-shaped obstacles.

  17. The inverse problems of wing panel manufacture processes

    NASA Astrophysics Data System (ADS)

    Oleinikov, A. I.; Bormotin, K. S.

    2013-12-01

    It is shown that inverse problems of steady-state creep bending of plates in both the geometrically linear and nonlinear formulations can be represented in a variational formulation. Steady-state values of the obtained functionals corresponding to the solutions of the problems of inelastic deformation and springback are determined by applying a finite element procedure to the functionals. Optimal laws of creep deformation are formulated using the criterion of minimizing damage in the functionals of the inverse problems. The formulated problems are reduced to the problems solved by the finite element method using MSC.Marc software. Currently, forming of light metals poses tremendous challenges due to their low ductility at room temperature and their unusual deformation characteristics at hot-cold work: strong asymmetry between tensile and compressive behavior, and a very pronounced anisotropy. We used the constitutive models of steady-state creep of initially transverse isotropy structural materials the kind of the stress state has influence. The paper gives basics of the developed computer-aided system of design, modeling, and electronic simulation targeting the processes of manufacture of wing integral panels. The modeling results can be used to calculate the die tooling, determine the panel processibility, and control panel rejection in the course of forming.

  18. Stochastic reduced order models for inverse problems under uncertainty

    PubMed Central

    Warner, James E.; Aquino, Wilkins; Grigoriu, Mircea D.

    2014-01-01

    This work presents a novel methodology for solving inverse problems under uncertainty using stochastic reduced order models (SROMs). Given statistical information about an observed state variable in a system, unknown parameters are estimated probabilistically through the solution of a model-constrained, stochastic optimization problem. The point of departure and crux of the proposed framework is the representation of a random quantity using a SROM - a low dimensional, discrete approximation to a continuous random element that permits e cient and non-intrusive stochastic computations. Characterizing the uncertainties with SROMs transforms the stochastic optimization problem into a deterministic one. The non-intrusive nature of SROMs facilitates e cient gradient computations for random vector unknowns and relies entirely on calls to existing deterministic solvers. Furthermore, the method is naturally extended to handle multiple sources of uncertainty in cases where state variable data, system parameters, and boundary conditions are all considered random. The new and widely-applicable SROM framework is formulated for a general stochastic optimization problem in terms of an abstract objective function and constraining model. For demonstration purposes, however, we study its performance in the specific case of inverse identification of random material parameters in elastodynamics. We demonstrate the ability to efficiently recover random shear moduli given material displacement statistics as input data. We also show that the approach remains effective for the case where the loading in the problem is random as well. PMID:25558115

  19. The inverse problems of wing panel manufacture processes

    SciTech Connect

    Oleinikov, A. I.; Bormotin, K. S.

    2013-12-16

    It is shown that inverse problems of steady-state creep bending of plates in both the geometrically linear and nonlinear formulations can be represented in a variational formulation. Steady-state values of the obtained functionals corresponding to the solutions of the problems of inelastic deformation and springback are determined by applying a finite element procedure to the functionals. Optimal laws of creep deformation are formulated using the criterion of minimizing damage in the functionals of the inverse problems. The formulated problems are reduced to the problems solved by the finite element method using MSC.Marc software. Currently, forming of light metals poses tremendous challenges due to their low ductility at room temperature and their unusual deformation characteristics at hot-cold work: strong asymmetry between tensile and compressive behavior, and a very pronounced anisotropy. We used the constitutive models of steady-state creep of initially transverse isotropy structural materials the kind of the stress state has influence. The paper gives basics of the developed computer-aided system of design, modeling, and electronic simulation targeting the processes of manufacture of wing integral panels. The modeling results can be used to calculate the die tooling, determine the panel processibility, and control panel rejection in the course of forming.

  20. Including geological information in the inverse problem of palaeothermal reconstruction

    NASA Astrophysics Data System (ADS)

    Trautner, S.; Nielsen, S. B.

    2003-04-01

    A reliable reconstruction of sediment thermal history is of central importance to the assessment of hydrocarbon potential and the understanding of basin evolution. However, only rarely do sedimentation history and borehole data in the form of present day temperatures and vitrinite reflectance constrain the past thermal evolution to a useful level of accuracy (Gallagher and Sambridge,1992; Nielsen,1998; Trautner and Nielsen,2003). This is reflected in the inverse solutions to the problem of determining heat flow history from borehole data: The recent heat flow is constrained by data while older values are governed by the chosen a prior heat flow. In this paper we reduce this problem by including geological information in the inverse problem. Through a careful analysis of geological and geophysical data the timing of the tectonic processes, which may influence heat flow, can be inferred. The heat flow history is then parameterised to allow for the temporal variations characteristic of the different tectonic events. The inversion scheme applies a Markov chain Monte Carlo (MCMC) approach (Nielsen and Gallagher, 1999; Ferrero and Gallagher,2002), which efficiently explores the model space and futhermore samples the posterior probability distribution of the model. The technique is demonstrated on wells in the northern North Sea with emphasis on the stretching event in Late Jurassic. The wells are characterised by maximum sediment temperature at the present day, which is the worst case for resolution of the past thermal history because vitrinite reflectance is determined mainly by the maximum temperature. Including geological information significantly improves the thermal resolution. Ferrero, C. and Gallagher,K.,2002. Stochastic thermal history modelling.1. Constraining heat flow histories and their uncertainty. Marine and Petroleum Geology, 19, 633-648. Gallagher,K. and Sambridge, M., 1992. The resolution of past heat flow in sedimentary basins from non-linear inversion

  1. Inverse problems with Poisson data: statistical regularization theory, applications and algorithms

    NASA Astrophysics Data System (ADS)

    Hohage, Thorsten; Werner, Frank

    2016-09-01

    Inverse problems with Poisson data arise in many photonic imaging modalities in medicine, engineering and astronomy. The design of regularization methods and estimators for such problems has been studied intensively over the last two decades. In this review we give an overview of statistical regularization theory for such problems, the most important applications, and the most widely used algorithms. The focus is on variational regularization methods in the form of penalized maximum likelihood estimators, which can be analyzed in a general setup. Complementing a number of recent convergence rate results we will establish consistency results. Moreover, we discuss estimators based on a wavelet-vaguelette decomposition of the (necessarily linear) forward operator. As most prominent applications we briefly introduce Positron emission tomography, inverse problems in fluorescence microscopy, and phase retrieval problems. The computation of a penalized maximum likelihood estimator involves the solution of a (typically convex) minimization problem. We also review several efficient algorithms which have been proposed for such problems over the last five years.

  2. Source localization in electromyography using the inverse potential problem

    NASA Astrophysics Data System (ADS)

    van den Doel, Kees; Ascher, Uri M.; Pai, Dinesh K.

    2011-02-01

    We describe an efficient method for reconstructing the activity in human muscles from an array of voltage sensors on the skin surface. MRI is used to obtain morphometric data which are segmented into muscle tissue, fat, bone and skin, from which a finite element model for volume conduction is constructed. The inverse problem of finding the current sources in the muscles is solved using a careful regularization technique which adds a priori information, yielding physically reasonable solutions from among those that satisfy the basic potential problem. Several regularization functionals are considered and numerical experiments on a 2D test model are performed to determine which performs best. The resulting scheme leads to numerical difficulties when applied to large-scale 3D problems. We clarify the nature of these difficulties and provide a method to overcome them, which is shown to perform well in the large-scale problem setting.

  3. Inference in infinite-dimensional inverse problems - Discretization and duality

    NASA Technical Reports Server (NTRS)

    Stark, Philip B.

    1992-01-01

    Many techniques for solving inverse problems involve approximating the unknown model, a function, by a finite-dimensional 'discretization' or parametric representation. The uncertainty in the computed solution is sometimes taken to be the uncertainty within the parametrization; this can result in unwarranted confidence. The theory of conjugate duality can overcome the limitations of discretization within the 'strict bounds' formalism, a technique for constructing confidence intervals for functionals of the unknown model incorporating certain types of prior information. The usual computational approach to strict bounds approximates the 'primal' problem in a way that the resulting confidence intervals are at most long enough to have the nominal coverage probability. There is another approach based on 'dual' optimization problems that gives confidence intervals with at least the nominal coverage probability. The pair of intervals derived by the two approaches bracket a correct confidence interval. The theory is illustrated with gravimetric, seismic, geomagnetic, and helioseismic problems and a numerical example in seismology.

  4. Inverse problem for in vivo NMR spatial localization

    SciTech Connect

    Hasenfeld, A.C.

    1985-11-01

    The basic physical problem of NMR spatial localization is considered. To study diseased sites, one must solve the problem of adequately localizing the NMR signal. We formulate this as an inverse problem. As the NMR Bloch equations determine the motion of nuclear spins in applied magnetic fields, a theoretical study is undertaken to answer the question of how to design magnetic field configurations to achieve these localized excited spin populations. Because of physical constraints in the production of the relevant radiofrequency fields, the problem factors into a temporal one and a spatial one. We formulate the temporal problem as a nonlinear transformation, called the Bloch Transform, from the rf input to the magnetization response. In trying to invert this transformation, both linear (for the Fourier Transform) and nonlinear (for the Bloch Transform) modes of radiofrequency excitation are constructed. The spatial problem is essentially a statics problem for the Maxwell equations of electromagnetism, as the wavelengths of the radiation considered are on the order of ten meters, and so propagation effects are negligible. In the general case, analytic solutions are unavailable, and so the methods of computer simulation are used to map the rf field spatial profiles. Numerical experiments are also performed to verify the theoretical analysis, and experimental confirmation of the theory is carried out on the 0.5 Tesla IBM/Oxford Imaging Spectrometer at the LBL NMR Medical Imaging Facility. While no explicit inverse is constructed to ''solve'' this problem, the combined theoretical/numerical analysis is validated experimentally, justifying the approximations made. 56 refs., 31 figs.

  5. Efficient solution of an inverse problem in cell population dynamics

    NASA Astrophysics Data System (ADS)

    Groh, Andreas; Krebs, Jochen; Wagner, Mathias

    2011-06-01

    In this paper, a size-structured model for cell division is examined and the question of determining the division (birth) rate from a measurable stable size distribution of the population is addressed. This inverse problem can be formulated as a differential-dilation equation. We propose a novel solution scheme based on mollification. The method of approximate inverse allows us to shift the derivative from the data to a precomputable reconstruction kernel. To comprise all available a priori information, a presmoothing step based on regression in reproducing kernel Hilbert spaces is introduced. We establish an error theory for the emerging algorithm, prove convergence and deduce a parameter strategy. The results are substantiated with extensive numerical tests both for artificial and real data based on proliferating tumor cells.

  6. Inverse problems in heterogeneous and fractured media using peridynamics

    SciTech Connect

    Turner, Daniel Z.; van Bloemen Waanders, Bart G.; Parks, Michael L.

    2015-12-10

    The following work presents an adjoint-based methodology for solving inverse problems in heterogeneous and fractured media using state-based peridynamics. We show that the inner product involving the peridynamic operators is self-adjoint. The proposed method is illustrated for several numerical examples with constant and spatially varying material parameters as well as in the context of fractures. We also present a framework for obtaining material parameters by integrating digital image correlation (DIC) with inverse analysis. This framework is demonstrated by evaluating the bulk and shear moduli for a sample of nuclear graphite using digital photographs taken during the experiment. The resulting measured values correspond well with other results reported in the literature. Lastly, we show that this framework can be used to determine the load state given observed measurements of a crack opening. Furthermore, this type of analysis has many applications in characterizing subsurface stress-state conditions given fracture patterns in cores of geologic material.

  7. Principal Component Geostatistical Approach for large-dimensional inverse problems

    PubMed Central

    Kitanidis, P K; Lee, J

    2014-01-01

    The quasi-linear geostatistical approach is for weakly nonlinear underdetermined inverse problems, such as Hydraulic Tomography and Electrical Resistivity Tomography. It provides best estimates as well as measures for uncertainty quantification. However, for its textbook implementation, the approach involves iterations, to reach an optimum, and requires the determination of the Jacobian matrix, i.e., the derivative of the observation function with respect to the unknown. Although there are elegant methods for the determination of the Jacobian, the cost is high when the number of unknowns, m, and the number of observations, n, is high. It is also wasteful to compute the Jacobian for points away from the optimum. Irrespective of the issue of computing derivatives, the computational cost of implementing the method is generally of the order of m2n, though there are methods to reduce the computational cost. In this work, we present an implementation that utilizes a matrix free in terms of the Jacobian matrix Gauss-Newton method and improves the scalability of the geostatistical inverse problem. For each iteration, it is required to perform K runs of the forward problem, where K is not just much smaller than m but can be smaller that n. The computational and storage cost of implementation of the inverse procedure scales roughly linearly with m instead of m2 as in the textbook approach. For problems of very large m, this implementation constitutes a dramatic reduction in computational cost compared to the textbook approach. Results illustrate the validity of the approach and provide insight in the conditions under which this method perform best. PMID:25558113

  8. Principal Component Geostatistical Approach for large-dimensional inverse problems

    NASA Astrophysics Data System (ADS)

    Kitanidis, P. K.; Lee, J.

    2014-07-01

    The quasi-linear geostatistical approach is for weakly nonlinear underdetermined inverse problems, such as Hydraulic Tomography and Electrical Resistivity Tomography. It provides best estimates as well as measures for uncertainty quantification. However, for its textbook implementation, the approach involves iterations, to reach an optimum, and requires the determination of the Jacobian matrix, i.e., the derivative of the observation function with respect to the unknown. Although there are elegant methods for the determination of the Jacobian, the cost is high when the number of unknowns, m, and the number of observations, n, is high. It is also wasteful to compute the Jacobian for points away from the optimum. Irrespective of the issue of computing derivatives, the computational cost of implementing the method is generally of the order of m2n, though there are methods to reduce the computational cost. In this work, we present an implementation that utilizes a matrix free in terms of the Jacobian matrix Gauss-Newton method and improves the scalability of the geostatistical inverse problem. For each iteration, it is required to perform K runs of the forward problem, where K is not just much smaller than m but can be smaller that n. The computational and storage cost of implementation of the inverse procedure scales roughly linearly with m instead of m2 as in the textbook approach. For problems of very large m, this implementation constitutes a dramatic reduction in computational cost compared to the textbook approach. Results illustrate the validity of the approach and provide insight in the conditions under which this method perform best.

  9. A comparative study among different regularization techniques for solving ill-posed magnetic inverse problem

    NASA Astrophysics Data System (ADS)

    Abdelazeem, Maha; Gobashy, Mohamed

    2015-04-01

    The magnetic inverse problem is, intrinsically, non-unique and its numerical solution is unstable. This means that any small perturbation in the data (noise) causes large variation in the solution. This ill-posedness is not only due to complex geological situations, but it may arise because of ill-conditioned kernel matrix. Procedures adopted to stabilize the inversion of ill-posed problem are called regularization, so the selection of regularization parameter is very important to invert the earth model causing the measured magnetic field. Two strategies are commonly used, techniques based on Tikhonov formula and techniques using the trust region sub-problem TRS and the controlling factor will be the radius of such region. In this study, the two categories are compared to examine the stability of solutions with noise. A MATLAB-based inversion code is implemented and tested on some synthetic total magnetic fields with different noise levels added to simulate real fields. The capability of such techniques have been further tested by applying it to real data.

  10. Numerical solution of the imprecisely defined inverse heat conduction problem

    NASA Astrophysics Data System (ADS)

    Smita, Tapaswini; Chakraverty, S.; Diptiranjan, Behera

    2015-05-01

    This paper investigates the numerical solution of the uncertain inverse heat conduction problem. Uncertainties present in the system parameters are modelled through triangular convex normalized fuzzy sets. In the solution process, double parametric forms of fuzzy numbers are used with the variational iteration method (VIM). This problem first computes the uncertain temperature distribution in the domain. Next, when the uncertain temperature measurements in the domain are known, the functions describing the uncertain temperature and heat flux on the boundary are reconstructed. Related example problems are solved using the present procedure. We have also compared the present results with those in [Inf. Sci. (2008) 178 1917] along with homotopy perturbation method (HPM) and [Int. Commun. Heat Mass Transfer (2012) 39 30] in the special cases to demonstrate the validity and applicability.

  11. Inverse atmospheric radiative transfer problems - A nonlinear minimization search method of solution. [aerosol pollution monitoring

    NASA Technical Reports Server (NTRS)

    Fymat, A. L.

    1976-01-01

    The paper studies the inversion of the radiative transfer equation describing the interaction of electromagnetic radiation with atmospheric aerosols. The interaction can be considered as the propagation in the aerosol medium of two light beams: the direct beam in the line-of-sight attenuated by absorption and scattering, and the diffuse beam arising from scattering into the viewing direction, which propagates more or less in random fashion. The latter beam has single scattering and multiple scattering contributions. In the former case and for single scattering, the problem is reducible to first-kind Fredholm equations, while for multiple scattering it is necessary to invert partial integrodifferential equations. A nonlinear minimization search method, applicable to the solution of both types of problems has been developed, and is applied here to the problem of monitoring aerosol pollution, namely the complex refractive index and size distribution of aerosol particles.

  12. Structure-approximating inverse protein folding problem in the 2D HP model.

    PubMed

    Gupta, Arvind; Manuch, Ján; Stacho, Ladislav

    2005-12-01

    The inverse protein folding problem is that of designing an amino acid sequence which has a particular native protein fold. This problem arises in drug design where a particular structure is necessary to ensure proper protein-protein interactions. In this paper, we show that in the 2D HP model of Dill it is possible to solve this problem for a broad class of structures. These structures can be used to closely approximate any given structure. One of the most important properties of a good protein (in drug design) is its stability--the aptitude not to fold simultaneously into other structures. We show that for a number of basic structures, our sequences have a unique fold. PMID:16379538

  13. Introduction to the 30th volume of Inverse Problems

    NASA Astrophysics Data System (ADS)

    Louis, Alfred K.

    2014-01-01

    The field of inverse problems is a fast-developing domain of research originating from the practical demands of finding the cause when a result is observed. The woodpecker, searching for insects, is probing a tree using sound waves: the information searched for is whether there is an insect or not, hence a 0-1 decision. When the result has to contain more information, ad hoc solutions are not at hand and more sophisticated methods have to be developed. Right from its first appearance, the field of inverse problems has been characterized by an interdisciplinary nature: the interpretation of measured data, reinforced by mathematical models serving the analyzing questions of observability, stability and resolution, developing efficient, stable and accurate algorithms to gain as much information as possible from the input and to feedback to the questions of optimal measurement configuration. As is typical for a new area of research, facets of it are separated and studied independently. Hence, fields such as the theory of inverse scattering, tomography in general and regularization methods have developed. However, all aspects have to be reassembled to arrive at the best possible solution to the problem at hand. This development is reflected by the first and still leading journal in the field, Inverse Problems. Founded by pioneers Roy Pike from London and Pierre Sabatier from Montpellier, who enjoyably describes the journal's nascence in his book Rêves et Combats d'un Enseignant-Chercheur, Retour Inverse [1], the journal has developed successfully over the last few decades. Neither the Editors-in-Chief, formerly called Honorary Editors, nor the board or authors could have set the path to success alone. Their fruitful interplay, complemented by the efficient and highly competent publishing team at IOP Publishing, has been fundamental. As such it is my honor and pleasure to follow my renowned colleagues Pierre Sabatier, Mario Bertero, Frank Natterer, Alberto Grünbaum and

  14. Negative Compressibility and Inverse Problem for Spinning Gas

    SciTech Connect

    Vasily Geyko and Nathaniel J. Fisch

    2013-01-11

    A spinning ideal gas in a cylinder with a smooth surface is shown to have unusual properties. First, under compression parallel to the axis of rotation, the spinning gas exhibits negative compressibility because energy can be stored in the rotation. Second, the spinning breaks the symmetry under which partial pressures of a mixture of gases simply add proportional to the constituent number densities. Thus, remarkably, in a mixture of spinning gases, an inverse problem can be formulated such that the gas constituents can be determined through external measurements only.

  15. Solving the inverse problem of magnetisation-stress resolution

    NASA Astrophysics Data System (ADS)

    Staples, S. G. H.; Vo, C.; Cowell, D. M. J.; Freear, S.; Ives, C.; Varcoe, B. T. H.

    2013-04-01

    Magnetostriction in various metals has been known since 1842, recently the focus has shifted away from ferrous metals, towards materials with a straightforward or exaggerated stress magnetostriction relationship. However, there is an increasing interest in understanding ferrous metal relationships, especially steels, because of its widespread use in building structures, transportation infrastructure, and pipelines. The aim of this paper is to solve the inverse problem of determining stress from an observed magnetic field which implies a given magnetic structure and to demonstrate that theoretical calculations using a multi-physics modeling technique agree with this experimental observation.

  16. An analytic method for the inverse problem of MREPT

    NASA Astrophysics Data System (ADS)

    Palamodov, V.

    2016-03-01

    Magnetic resonance electric properties tomography (MREPT) is a medical imaging modality for visualizing the electrical tissue properties of the human body using radio-frequency magnetic fields. This method consists of reconstructing the admittivity distribution from the positive rotating component of the magnetic field. In the newest paper of Ammari et al (2015 Inverse Problems 31 105001) an approximate method of reconstruction of variable admittivity was proposed. In this paper a method for exact reconstruction of the admittivity from data of the positive rotating component of the field is given.

  17. On Lambda and Time Operators: the Inverse Intertwining Problem Revisited

    NASA Astrophysics Data System (ADS)

    Gómez-Cubillo, F.; Suchanecki, Z.; Villullas, S.

    2011-07-01

    An exact theory of irreversibility was proposed by Misra, Prigogine and Courbage, based on non-unitary similarity transformations Λ that intertwine reversible dynamics and irreversible ones. This would advocate the idea that irreversible behavior would originate at the microscopic level. Reversible evolution with an internal time operator have the intertwining property. Recently the inverse intertwining problem has been answered in the negative, that is, not every unitary evolution allowing such Λ-transformation has an internal time. This work contributes new results in this direction.

  18. A spatiotemporal dynamic distributed solution to the MEG inverse problem

    PubMed Central

    Lamus, Camilo; Hämäläinen, Matti S.; Temereanca, Simona; Brown, Emery N.; Purdon, Patrick L.

    2012-01-01

    MEG/EEG are non-invasive imaging techniques that record brain activity with high temporal resolution. However, estimation of brain source currents from surface recordings requires solving an ill-conditioned inverse problem. Converging lines of evidence in neuroscience, from neuronal network models to resting-state imaging and neurophysiology, suggest that cortical activation is a distributed spatiotemporal dynamic process, supported by both local and long-distance neuroanatomic connections. Because spatiotemporal dynamics of this kind are central to brain physiology, inverse solutions could be improved by incorporating models of these dynamics. In this article, we present a model for cortical activity based on nearest-neighbor autoregression that incorporates local spatiotemporal interactions between distributed sources in a manner consistent with neurophysiology and neuroanatomy. We develop a dynamic Maximum a Posteriori Expectation-Maximization (dMAP-EM) source localization algorithm for estimation of cortical sources and model parameters based on the Kalman Filter, the Fixed Interval Smoother, and the EM algorithms. We apply the dMAP-EM algorithm to simulated experiments as well as to human experimental data. Furthermore, we derive expressions to relate our dynamic estimation formulas to those of standard static models, and show how dynamic methods optimally assimilate past and future data. Our results establish the feasibility of spatiotemporal dynamic estimation in large-scale distributed source spaces with several thousand source locations and hundreds of sensors, with resulting inverse solutions that provide substantial performance improvements over static methods. PMID:22155043

  19. Inverse problem of quadratic time-dependent Hamiltonians

    NASA Astrophysics Data System (ADS)

    Guo, Guang-Jie; Meng, Yan; Chang, Hong; Duan, Hui-Zeng; Di, Bing

    2015-08-01

    Using an algebraic approach, it is possible to obtain the temporal evolution wave function for a Gaussian wave-packet obeying the quadratic time-dependent Hamiltonian (QTDH). However, in general, most of the practical cases are not exactly solvable, for we need general solutions of the Riccatti equations which are not generally known. We therefore bypass directly solving for the temporal evolution wave function, and study its inverse problem. We start with a particular evolution of the wave-packet, and get the required Hamiltonian by using the inverse method. The inverse approach opens up a new way to find new exact solutions to the QTDH. Some typical examples are studied in detail. For a specific time-dependent periodic harmonic oscillator, the Berry phase is obtained exactly. Project supported by the National Natural Science Foundation of China (Grant No. 11347171), the Natural Science Foundation of Hebei Province of China (Grant No. A2012108003), and the Key Project of Educational Commission of Hebei Province of China (Grant No. ZD2014052).

  20. A verifiable solution to the MEG inverse problem.

    PubMed

    Barnes, Gareth R; Furlong, Paul L; Singh, Krish D; Hillebrand, Arjan

    2006-06-01

    Magnetoencephalography (MEG) is a non-invasive brain imaging technique with the potential for very high temporal and spatial resolution of neuronal activity. The main stumbling block for the technique has been that the estimation of a neuronal current distribution, based on sensor data outside the head, is an inverse problem with an infinity of possible solutions. Many inversion techniques exist, all using different a-priori assumptions in order to reduce the number of possible solutions. Although all techniques can be thoroughly tested in simulation, implicit in the simulations are the experimenter's own assumptions about realistic brain function. To date, the only way to test the validity of inversions based on real MEG data has been through direct surgical validation, or through comparison with invasive primate data. In this work, we constructed a null hypothesis that the reconstruction of neuronal activity contains no information on the distribution of the cortical grey matter. To test this, we repeatedly compared rotated sections of grey matter with a beamformer estimate of neuronal activity to generate a distribution of mutual information values. The significance of the comparison between the un-rotated anatomical information and the electrical estimate was subsequently assessed against this distribution. We found that there was significant (P < 0.05) anatomical information contained in the beamformer images across a number of frequency bands. Based on the limited data presented here, we can say that the assumptions behind the beamformer algorithm are not unreasonable for the visual-motor task investigated. PMID:16480896

  1. A spatiotemporal dynamic distributed solution to the MEG inverse problem.

    PubMed

    Lamus, Camilo; Hämäläinen, Matti S; Temereanca, Simona; Brown, Emery N; Purdon, Patrick L

    2012-11-01

    MEG/EEG are non-invasive imaging techniques that record brain activity with high temporal resolution. However, estimation of brain source currents from surface recordings requires solving an ill-conditioned inverse problem. Converging lines of evidence in neuroscience, from neuronal network models to resting-state imaging and neurophysiology, suggest that cortical activation is a distributed spatiotemporal dynamic process, supported by both local and long-distance neuroanatomic connections. Because spatiotemporal dynamics of this kind are central to brain physiology, inverse solutions could be improved by incorporating models of these dynamics. In this article, we present a model for cortical activity based on nearest-neighbor autoregression that incorporates local spatiotemporal interactions between distributed sources in a manner consistent with neurophysiology and neuroanatomy. We develop a dynamic maximum a posteriori expectation-maximization (dMAP-EM) source localization algorithm for estimation of cortical sources and model parameters based on the Kalman Filter, the Fixed Interval Smoother, and the EM algorithms. We apply the dMAP-EM algorithm to simulated experiments as well as to human experimental data. Furthermore, we derive expressions to relate our dynamic estimation formulas to those of standard static models, and show how dynamic methods optimally assimilate past and future data. Our results establish the feasibility of spatiotemporal dynamic estimation in large-scale distributed source spaces with several thousand source locations and hundreds of sensors, with resulting inverse solutions that provide substantial performance improvements over static methods. PMID:22155043

  2. Evaluation of simplified evaporation duct refractivity models for inversion problems

    NASA Astrophysics Data System (ADS)

    Saeger, J. T.; Grimes, N. G.; Rickard, H. E.; Hackett, E. E.

    2015-10-01

    To assess a radar system's instantaneous performance on any given day, detailed knowledge of the meteorological conditions is required due to the dependency of atmospheric refractivity on thermodynamic properties such as temperature, water vapor, and pressure. Because of the significant challenges involved in obtaining these data, recent efforts have focused on development of methods to obtain the refractivity structure inversely using radar measurements and radar wave propagation models. Such inversion techniques generally use simplified refractivity models in order to reduce the parameter space of the solution. Here the accuracy of three simple refractivity models is examined for the case of an evaporation duct. The models utilize the basic log linear shape classically associated with evaporation ducts, but each model depends on various parameters that affect different aspects of the profile, such as its shape and duct height. The model parameters are optimized using radiosonde data, and their performance is compared to these atmospheric measurements. The optimized models and data are also used to predict propagation using a parabolic equation code with the refractivity prescribed by the models and measured data, and the resulting propagation patterns are compared. The results of this study suggest that the best log linear model formulation for an inversion problem would be a two-layer model that contains at least three parameters: duct height, duct curvature, and mixed layer slope. This functional form permits a reasonably accurate fit to atmospheric measurements as well as embodies key features of the profile required for correct propagation prediction with as few parameters as possible.

  3. Topological inversion for solution of geodesy-constrained geophysical problems

    NASA Astrophysics Data System (ADS)

    Saltogianni, Vasso; Stiros, Stathis

    2015-04-01

    Geodetic data, mostly GPS observations, permit to measure displacements of selected points around activated faults and volcanoes, and on the basis of geophysical models, to model the underlying physical processes. This requires inversion of redundant systems of highly non-linear equations with >3 unknowns; a situation analogous to the adjustment of geodetic networks. However, in geophysical problems inversion cannot be based on conventional least-squares techniques, and is based on numerical inversion techniques (a priori fixing of some variables, optimization in steps with values of two variables each time to be regarded fixed, random search in the vicinity of approximate solutions). Still these techniques lead to solutions trapped in local minima, to correlated estimates and to solutions with poor error control (usually sampling-based approaches). To overcome these problems, a numerical-topological, grid-search based technique in the RN space is proposed (N the number of unknown variables). This technique is in fact a generalization and refinement of techniques used in lighthouse positioning and in some cases of low-accuracy 2-D positioning using Wi-Fi etc. The basic concept is to assume discrete possible ranges of each variable, and from these ranges to define a grid G in the RN space, with some of the gridpoints to approximate the true solutions of the system. Each point of hyper-grid G is then tested whether it satisfies the observations, given their uncertainty level, and successful grid points define a sub-space of G containing the true solutions. The optimal (minimal) space containing one or more solutions is obtained using a trial-and-error approach, and a single optimization factor. From this essentially deterministic identification of the set of gridpoints satisfying the system of equations, at a following step, a stochastic optimal solution is computed corresponding to the center of gravity of this set of gridpoints. This solution corresponds to a

  4. Review on solving the inverse problem in EEG source analysis

    PubMed Central

    Grech, Roberta; Cassar, Tracey; Muscat, Joseph; Camilleri, Kenneth P; Fabri, Simon G; Zervakis, Michalis; Xanthopoulos, Petros; Sakkalis, Vangelis; Vanrumste, Bart

    2008-01-01

    In this primer, we give a review of the inverse problem for EEG source localization. This is intended for the researchers new in the field to get insight in the state-of-the-art techniques used to find approximate solutions of the brain sources giving rise to a scalp potential recording. Furthermore, a review of the performance results of the different techniques is provided to compare these different inverse solutions. The authors also include the results of a Monte-Carlo analysis which they performed to compare four non parametric algorithms and hence contribute to what is presently recorded in the literature. An extensive list of references to the work of other researchers is also provided. This paper starts off with a mathematical description of the inverse problem and proceeds to discuss the two main categories of methods which were developed to solve the EEG inverse problem, mainly the non parametric and parametric methods. The main difference between the two is to whether a fixed number of dipoles is assumed a priori or not. Various techniques falling within these categories are described including minimum norm estimates and their generalizations, LORETA, sLORETA, VARETA, S-MAP, ST-MAP, Backus-Gilbert, LAURA, Shrinking LORETA FOCUSS (SLF), SSLOFO and ALF for non parametric methods and beamforming techniques, BESA, subspace techniques such as MUSIC and methods derived from it, FINES, simulated annealing and computational intelligence algorithms for parametric methods. From a review of the performance of these techniques as documented in the literature, one could conclude that in most cases the LORETA solution gives satisfactory results. In situations involving clusters of dipoles, higher resolution algorithms such as MUSIC or FINES are however preferred. Imposing reliable biophysical and psychological constraints, as done by LAURA has given superior results. The Monte-Carlo analysis performed, comparing WMN, LORETA, sLORETA and SLF, for different noise levels

  5. Inverse Problems in Complex Models and Applications to Earth Sciences

    NASA Astrophysics Data System (ADS)

    Bosch, M. E.

    2015-12-01

    The inference of the subsurface earth structure and properties requires the integration of different types of data, information and knowledge, by combined processes of analysis and synthesis. To support the process of integrating information, the regular concept of data inversion is evolving to expand its application to models with multiple inner components (properties, scales, structural parameters) that explain multiple data (geophysical survey data, well-logs, core data). The probabilistic inference methods provide the natural framework for the formulation of these problems, considering a posterior probability density function (PDF) that combines the information from a prior information PDF and the new sets of observations. To formulate the posterior PDF in the context of multiple datasets, the data likelihood functions are factorized assuming independence of uncertainties for data originating across different surveys. A realistic description of the earth medium requires modeling several properties and structural parameters, which relate to each other according to dependency and independency notions. Thus, conditional probabilities across model components also factorize. A common setting proceeds by structuring the model parameter space in hierarchical layers. A primary layer (e.g. lithology) conditions a secondary layer (e.g. physical medium properties), which conditions a third layer (e.g. geophysical data). In general, less structured relations within model components and data emerge from the analysis of other inverse problems. They can be described with flexibility via direct acyclic graphs, which are graphs that map dependency relations between the model components. Examples of inverse problems in complex models can be shown at various scales. At local scale, for example, the distribution of gas saturation is inferred from pre-stack seismic data and a calibrated rock-physics model. At regional scale, joint inversion of gravity and magnetic data is applied

  6. Flexible hydraulic DFN for direct and inverse flow problems

    NASA Astrophysics Data System (ADS)

    de Dreuzy, J.; Le Goc, R.; Davy, P.; Bour, O.

    2006-12-01

    Natural fractured media are complex because of the large number of fractures of widely-scattered characteristics (length, transmissivity). Simulating hydraulic processes by accounting for the natural complexity requires flexible software and high performance computing. We develop a parallel fracture network software MPFRAC (for Massive Parallel FRACture network) designed to handle 2D and 3D large networks with fractures of different shapes (ellipses and polygons) delimited by a polyhedron. MPFRAC contains a broad range of fracture generation possibilities to which can be easily added other methods. In 3D, a specific module creates a consistent meshing of the fractures matching at the fracture intersections. The quality of the meshing required for finite element methods is ensured by projecting and the discretizing fracture intersections on a regular grid. Flow is computed on the meshed structure by a mixed finite element scheme and the yielded system can be solved using parallel solvers. Particle transport is subsequently simulated by a particle transport method. All modules of MPFRAC are independent and parallel enhancing the software flexibility and performances. High performance computing on PC-clusters is especially required for 3D flow and transport simulations. We use this software first for direct simulation to determine the influence of the fracture characteristics on the equivalent permeability, flow structure and pump test responses. Based on stochastic fracture distribution characteristics (fracture length, transmissivity, orientation and position distributions) Monte-Carlo simulations are performed in order to extract equivalent simplified flow models conditioned by a very limited number of parameters that could be taken as the a priori reduced search space of the inverse problem. Secondly, we use it for inverse problems currently in 2D for studying the influence of the data density on the identification of a limited number of fracture characteristics

  7. Galerkin approximation for inverse problems for nonautonomous nonlinear distributed systems

    NASA Technical Reports Server (NTRS)

    Banks, H. T.; Reich, Simeon; Rosen, I. G.

    1988-01-01

    An abstract framework and convergence theory is developed for Galerkin approximation for inverse problems involving the identification of nonautonomous nonlinear distributed parameter systems. A set of relatively easily verified conditions is provided which are sufficient to guarantee the existence of optimal solutions and their approximation by a sequence of solutions to a sequence of approximating finite dimensional identification problems. The approach is based on the theory of monotone operators in Banach spaces and is applicable to a reasonably broad class of nonlinear distributed systems. Operator theoretic and variational techniques are used to establish a fundamental convergence result. An example involving evolution systems with dynamics described by nonstationary quasilinear elliptic operators along with some applications are presented and discussed.

  8. Regularized total least squares approach for nonconvolutional linear inverse problems.

    PubMed

    Zhu, W; Wang, Y; Galatsanos, N P; Zhang, J

    1999-01-01

    In this correspondence, a solution is developed for the regularized total least squares (RTLS) estimate in linear inverse problems where the linear operator is nonconvolutional. Our approach is based on a Rayleigh quotient (RQ) formulation of the TLS problem, and we accomplish regularization by modifying the RQ function to enforce a smooth solution. A conjugate gradient algorithm is used to minimize the modified RQ function. As an example, the proposed approach has been applied to the perturbation equation encountered in optical tomography. Simulation results show that this method provides more stable and accurate solutions than the regularized least squares and a previously reported total least squares approach, also based on the RQ formulation. PMID:18267442

  9. Inverse problems and computational cell metabolic models: a statistical approach

    NASA Astrophysics Data System (ADS)

    Calvetti, D.; Somersalo, E.

    2008-07-01

    In this article, we give an overview of the Bayesian modelling of metabolic systems at the cellular and subcellular level. The models are based on detailed description of key biochemical reactions occurring in tissue, which may in turn be compartmentalized into cytosol and mitochondria, and of transports between the compartments. The classical deterministic approach which models metabolic systems as dynamical systems with Michaelis-Menten kinetics, is replaced by a stochastic extension where the model parameters are interpreted as random variables with an appropriate probability density. The inverse problem of cell metabolism in this setting consists of estimating the density of the model parameters. After discussing some possible approaches to solving the problem, we address the issue of how to assess the reliability of the predictions of a stochastic model by proposing an output analysis in terms of model uncertainties. Visualization modalities for organizing the large amount of information provided by the Bayesian dynamic sensitivity analysis are also illustrated.

  10. Inverse spin glass and related maximum entropy problems.

    PubMed

    Castellana, Michele; Bialek, William

    2014-09-12

    If we have a system of binary variables and we measure the pairwise correlations among these variables, then the least structured or maximum entropy model for their joint distribution is an Ising model with pairwise interactions among the spins. Here we consider inhomogeneous systems in which we constrain, for example, not the full matrix of correlations, but only the distribution from which these correlations are drawn. In this sense, what we have constructed is an inverse spin glass: rather than choosing coupling constants at random from a distribution and calculating correlations, we choose the correlations from a distribution and infer the coupling constants. We argue that such models generate a block structure in the space of couplings, which provides an explicit solution of the inverse problem. This allows us to generate a phase diagram in the space of (measurable) moments of the distribution of correlations. We expect that these ideas will be most useful in building models for systems that are nonequilibrium statistical mechanics problems, such as networks of real neurons. PMID:25260004

  11. Fisher information for inverse problems and trace class operators

    NASA Astrophysics Data System (ADS)

    Nordebo, S.; Gustafsson, M.; Khrennikov, A.; Nilsson, B.; Toft, J.

    2012-12-01

    This paper provides a mathematical framework for Fisher information analysis for inverse problems based on Gaussian noise on infinite-dimensional Hilbert space. The covariance operator for the Gaussian noise is assumed to be trace class, and the Jacobian of the forward operator Hilbert-Schmidt. We show that the appropriate space for defining the Fisher information is given by the Cameron-Martin space. This is mainly because the range space of the covariance operator always is strictly smaller than the Hilbert space. For the Fisher information to be well-defined, it is furthermore required that the range space of the Jacobian is contained in the Cameron-Martin space. In order for this condition to hold and for the Fisher information to be trace class, a sufficient condition is formulated based on the singular values of the Jacobian as well as of the eigenvalues of the covariance operator, together with some regularity assumptions regarding their relative rate of convergence. An explicit example is given regarding an electromagnetic inverse source problem with "external" spherically isotropic noise, as well as "internal" additive uncorrelated noise.

  12. Comparison of optimal design methods in inverse problems

    NASA Astrophysics Data System (ADS)

    Banks, H. T.; Holm, K.; Kappel, F.

    2011-07-01

    Typical optimal design methods for inverse or parameter estimation problems are designed to choose optimal sampling distributions through minimization of a specific cost function related to the resulting error in parameter estimates. It is hoped that the inverse problem will produce parameter estimates with increased accuracy using data collected according to the optimal sampling distribution. Here we formulate the classical optimal design problem in the context of general optimization problems over distributions of sampling times. We present a new Prohorov metric-based theoretical framework that permits one to treat succinctly and rigorously any optimal design criteria based on the Fisher information matrix. A fundamental approximation theory is also included in this framework. A new optimal design, SE-optimal design (standard error optimal design), is then introduced in the context of this framework. We compare this new design criterion with the more traditional D-optimal and E-optimal designs. The optimal sampling distributions from each design are used to compute and compare standard errors; the standard errors for parameters are computed using asymptotic theory or bootstrapping and the optimal mesh. We use three examples to illustrate ideas: the Verhulst-Pearl logistic population model (Banks H T and Tran H T 2009 Mathematical and Experimental Modeling of Physical and Biological Processes (Boca Raton, FL: Chapman and Hall/CRC)), the standard harmonic oscillator model (Banks H T and Tran H T 2009) and a popular glucose regulation model (Bergman R N, Ider Y Z, Bowden C R and Cobelli C 1979 Am. J. Physiol. 236 E667-77 De Gaetano A and Arino O 2000 J. Math. Biol. 40 136-68 Toffolo G, Bergman R N, Finegood D T, Bowden C R and Cobelli C 1980 Diabetes 29 979-90).

  13. Inverse zombies, anesthesia awareness, and the hard problem of unconsciousness.

    PubMed

    Mashour, George A; LaRock, Eric

    2008-12-01

    Philosophical (p-) zombies are constructs that possess all of the behavioral features and responses of a sentient human being, yet are not conscious. P-zombies are intimately linked to the hard problem of consciousness and have been invoked as arguments against physicalist approaches. But what if we were to invert the characteristics of p-zombies? Such an inverse (i-) zombie would possess all of the behavioral features and responses of an insensate being, yet would nonetheless be conscious. While p-zombies are logically possible but naturally improbable, an approximation of i-zombies actually exists: individuals experiencing what is referred to as "anesthesia awareness." Patients under general anesthesia may be intubated (preventing speech), paralyzed (preventing movement), and narcotized (minimizing response to nociceptive stimuli). Thus, they appear--and typically are--unconscious. In 1-2 cases/1000, however, patients may be aware of intraoperative events, sometimes without any objective indices. Furthermore, a much higher percentage of patients (22% in a recent study) may have the subjective experience of dreaming during general anesthesia. P-zombies confront us with the hard problem of consciousness--how do we explain the presence of qualia? I-zombies present a more practical problem--how do we detect the presence of qualia? The current investigation compares p-zombies to i-zombies and explores the "hard problem" of unconsciousness with a focus on anesthesia awareness. PMID:18635380

  14. Basis set expansion for inverse problems in plasma diagnostic analysis

    NASA Astrophysics Data System (ADS)

    Jones, B.; Ruiz, C. L.

    2013-07-01

    A basis set expansion method [V. Dribinski, A. Ossadtchi, V. A. Mandelshtam, and H. Reisler, Rev. Sci. Instrum. 73, 2634 (2002)], 10.1063/1.1482156 is applied to recover physical information about plasma radiation sources from instrument data, which has been forward transformed due to the nature of the measurement technique. This method provides a general approach for inverse problems, and we discuss two specific examples relevant to diagnosing fast z pinches on the 20-25 MA Z machine [M. E. Savage, L. F. Bennett, D. E. Bliss, W. T. Clark, R. S. Coats, J. M. Elizondo, K. R. LeChien, H. C. Harjes, J. M. Lehr, J. E. Maenchen, D. H. McDaniel, M. F. Pasik, T. D. Pointon, A. C. Owen, D. B. Seidel, D. L. Smith, B. S. Stoltzfus, K. W. Struve, W. A. Stygar, L. K. Warne, J. R. Woodworth, C. W. Mendel, K. R. Prestwich, R. W. Shoup, D. L. Johnson, J. P. Corley, K. C. Hodge, T. C. Wagoner, and P. E. Wakeland, in Proceedings of the Pulsed Power Plasma Sciences Conference (IEEE, 2007), p. 979]. First, Abel inversion of time-gated, self-emission x-ray images from a wire array implosion is studied. Second, we present an approach for unfolding neutron time-of-flight measurements from a deuterium gas puff z pinch to recover information about emission time history and energy distribution. Through these examples, we discuss how noise in the measured data limits the practical resolution of the inversion, and how the method handles discontinuities in the source function and artifacts in the projected image. We add to the method a propagation of errors calculation for estimating uncertainties in the inverted solution.

  15. Basis set expansion for inverse problems in plasma diagnostic analysis.

    PubMed

    Jones, B; Ruiz, C L

    2013-07-01

    A basis set expansion method [V. Dribinski, A. Ossadtchi, V. A. Mandelshtam, and H. Reisler, Rev. Sci. Instrum. 73, 2634 (2002)] is applied to recover physical information about plasma radiation sources from instrument data, which has been forward transformed due to the nature of the measurement technique. This method provides a general approach for inverse problems, and we discuss two specific examples relevant to diagnosing fast z pinches on the 20-25 MA Z machine [M. E. Savage, L. F. Bennett, D. E. Bliss, W. T. Clark, R. S. Coats, J. M. Elizondo, K. R. LeChien, H. C. Harjes, J. M. Lehr, J. E. Maenchen, D. H. McDaniel, M. F. Pasik, T. D. Pointon, A. C. Owen, D. B. Seidel, D. L. Smith, B. S. Stoltzfus, K. W. Struve, W. A. Stygar, L. K. Warne, J. R. Woodworth, C. W. Mendel, K. R. Prestwich, R. W. Shoup, D. L. Johnson, J. P. Corley, K. C. Hodge, T. C. Wagoner, and P. E. Wakeland, in Proceedings of the Pulsed Power Plasma Sciences Conference (IEEE, 2007), p. 979]. First, Abel inversion of time-gated, self-emission x-ray images from a wire array implosion is studied. Second, we present an approach for unfolding neutron time-of-flight measurements from a deuterium gas puff z pinch to recover information about emission time history and energy distribution. Through these examples, we discuss how noise in the measured data limits the practical resolution of the inversion, and how the method handles discontinuities in the source function and artifacts in the projected image. We add to the method a propagation of errors calculation for estimating uncertainties in the inverted solution. PMID:23902066

  16. Basis set expansion for inverse problems in plasma diagnostic analysis

    SciTech Connect

    Jones, B.; Ruiz, C. L.

    2013-07-15

    A basis set expansion method [V. Dribinski, A. Ossadtchi, V. A. Mandelshtam, and H. Reisler, Rev. Sci. Instrum. 73, 2634 (2002)] is applied to recover physical information about plasma radiation sources from instrument data, which has been forward transformed due to the nature of the measurement technique. This method provides a general approach for inverse problems, and we discuss two specific examples relevant to diagnosing fast z pinches on the 20–25 MA Z machine [M. E. Savage, L. F. Bennett, D. E. Bliss, W. T. Clark, R. S. Coats, J. M. Elizondo, K. R. LeChien, H. C. Harjes, J. M. Lehr, J. E. Maenchen, D. H. McDaniel, M. F. Pasik, T. D. Pointon, A. C. Owen, D. B. Seidel, D. L. Smith, B. S. Stoltzfus, K. W. Struve, W. A. Stygar, L. K. Warne, J. R. Woodworth, C. W. Mendel, K. R. Prestwich, R. W. Shoup, D. L. Johnson, J. P. Corley, K. C. Hodge, T. C. Wagoner, and P. E. Wakeland, in Proceedings of the Pulsed Power Plasma Sciences Conference (IEEE, 2007), p. 979]. First, Abel inversion of time-gated, self-emission x-ray images from a wire array implosion is studied. Second, we present an approach for unfolding neutron time-of-flight measurements from a deuterium gas puff z pinch to recover information about emission time history and energy distribution. Through these examples, we discuss how noise in the measured data limits the practical resolution of the inversion, and how the method handles discontinuities in the source function and artifacts in the projected image. We add to the method a propagation of errors calculation for estimating uncertainties in the inverted solution.

  17. Theoretical study of Laplacian electrocardiography forward and inverse problem

    NASA Astrophysics Data System (ADS)

    Wu, Dongsheng

    The present study is concerned with a fundamental problem of cardiac electrophysiology, that is relating in a quantitative way the electrical activity within the heart to the signals recorded over the body surface. By computer simulation, a rigorous evaluation of the performance of the body surface Laplacian electrocardiographic maps in a physiologically reasonable and well-controlled computational setting is provided in this dissertation. The present forward heart-torso model, a three-dimensional ventricular conduction model embedded in a realistically shaped inhomogeneous torso volume conductor model, represents, up to date, the most advanced computer model, which is available for studying the Laplacian electrocardiographic fields corresponding to normal and abnormal ventricular conduction processes. Theoretical studies show that one can achieve enhanced spatial resolution of mapping cardiac electrical activity by obtaining the Laplacian ECG over the body surface. The present work demonstrates the excellent performance of the body surface Laplacian electrocardiographic maps in resolving and imaging the underlying regional myocardial electrical activity. The biophysics underlying this is that the Laplacian ECG heavily weights the contributions from the myocardial bioelectric sources that are closest to the recording location, whereas the potential ECG sums up the contributions from a large area of activated myocardial tissue. It is this regional nature of the Laplacian ECG that makes it possible to provide a more localized body surface manifestation of the underlying regional myocardial electrical activity. The feasibility of applying the Laplacian ECG to the inverse problems has also been investigated. Theoretical studies of the Laplacian electrocardiogram based inverse problem by using a homogeneous spherical volume conductor and a realistically shaped volume conductor have been conducted. The present work shows encouraging results which suggest the feasibility

  18. Isotropic probability measures in infinite dimensional spaces: Inverse problems/prior information/stochastic inversion

    NASA Technical Reports Server (NTRS)

    Backus, George

    1987-01-01

    Let R be the real numbers, R(n) the linear space of all real n-tuples, and R(infinity) the linear space of all infinite real sequences x = (x sub 1, x sub 2,...). Let P sub n :R(infinity) approaches R(n) be the projection operator with P sub n (x) = (x sub 1,...,x sub n). Let p(infinity) be a probability measure on the smallest sigma-ring of subsets of R(infinity) which includes all of the cylinder sets P sub n(-1) (B sub n), where B sub n is an arbitrary Borel subset of R(n). Let p sub n be the marginal distribution of p(infinity) on R(n), so p sub n(B sub n) = p(infinity)(P sub n to the -1(B sub n)) for each B sub n. A measure on R(n) is isotropic if it is invariant under all orthogonal transformations of R(n). All members of the set of all isotropic probability distributions on R(n) are described. The result calls into question both stochastic inversion and Bayesian inference, as currently used in many geophysical inverse problems.

  19. Solution accelerators for large scale 3D electromagnetic inverse problems

    SciTech Connect

    Newman, Gregory A.; Boggs, Paul T.

    2004-04-05

    We provide a framework for preconditioning nonlinear 3D electromagnetic inverse scattering problems using nonlinear conjugate gradient (NLCG) and limited memory (LM) quasi-Newton methods. Key to our approach is the use of an approximate adjoint method that allows for an economical approximation of the Hessian that is updated at each inversion iteration. Using this approximate Hessian as a preconditoner, we show that the preconditioned NLCG iteration converges significantly faster than the non-preconditioned iteration, as well as converging to a data misfit level below that observed for the non-preconditioned method. Similar conclusions are also observed for the LM iteration; preconditioned with the approximate Hessian, the LM iteration converges faster than the non-preconditioned version. At this time, however, we see little difference between the convergence performance of the preconditioned LM scheme and the preconditioned NLCG scheme. A possible reason for this outcome is the behavior of the line search within the LM iteration. It was anticipated that, near convergence, a step size of one would be approached, but what was observed, instead, were step lengths that were nowhere near one. We provide some insights into the reasons for this behavior and suggest further research that may improve the performance of the LM methods.

  20. Inverse problems in heterogeneous and fractured media using peridynamics

    DOE PAGESBeta

    Turner, Daniel Z.; van Bloemen Waanders, Bart G.; Parks, Michael L.

    2015-12-10

    The following work presents an adjoint-based methodology for solving inverse problems in heterogeneous and fractured media using state-based peridynamics. We show that the inner product involving the peridynamic operators is self-adjoint. The proposed method is illustrated for several numerical examples with constant and spatially varying material parameters as well as in the context of fractures. We also present a framework for obtaining material parameters by integrating digital image correlation (DIC) with inverse analysis. This framework is demonstrated by evaluating the bulk and shear moduli for a sample of nuclear graphite using digital photographs taken during the experiment. The resulting measuredmore » values correspond well with other results reported in the literature. Lastly, we show that this framework can be used to determine the load state given observed measurements of a crack opening. Furthermore, this type of analysis has many applications in characterizing subsurface stress-state conditions given fracture patterns in cores of geologic material.« less

  1. The LHC Inverse Problem, Supersymmetry and the ILC

    SciTech Connect

    Berger, C.F.; Gainer, J.S.; Hewett, J.L.; Lillie, B.; Rizzo, T.G.

    2007-11-12

    We address the question whether the ILC can resolve the LHC Inverse Problem within the framework of the MSSM. We examine 242 points in the MSSM parameter space which were generated at random and were found to give indistinguishable signatures at the LHC. After a realistic simulation including full Standard Model backgrounds and a fast detector simulation, we find that roughly only one third of these scenarios lead to visible signatures of some kind with a significance {ge} 5 at the ILC with {radical}s = 500 GeV. Furthermore, we examine these points in parameter space pairwise and find that only one third of the pairs are distinguishable at the ILC at 5{sigma}.

  2. Multi-scale Plasmonic Nanoparticles and the Inverse Problem

    PubMed Central

    Odom, Teri W.; You, Eun-Ah; Sweeney, Christina M.

    2012-01-01

    This Perspective describes how multi-scale plasmonic structures with two or more length scales (fine, medium, coarse) provide an experimental route for addressing the inverse problem. Specific near-field and far-field optical properties can be targeted and compiled into a plasmon resonance library by taking advantage of length scales spanning three orders of magnitude, from 1 nm to greater than 1000 nm, in a single particle. Examples of multi-scale 1D, 2D, and 3D gold structures created by nanofabrication tools and templates are discussed, and unexpected optical properties compared to those from their smaller counterparts are emphasized. One application of multi-scale particle dimers for surface-enhanced Raman spectroscopy is also described. PMID:23066451

  3. Geomagnetic inverse problem and data assimilation: a progress report

    NASA Astrophysics Data System (ADS)

    Aubert, Julien; Fournier, Alexandre

    2013-04-01

    In this presentation I will present two studies recently undertaken by our group in an effort to bring the benefits of data assimilation to the study of Earth's magnetic field and the dynamics of its liquid iron core, where the geodynamo operates. In a first part I will focus on the geomagnetic inverse problem, which attempts to recover the fluid flow in the core from the temporal variation of the magnetic field (known as the secular variation). Geomagnetic data can be downward continued from the surface of the Earth down to the core-mantle boundary, but not further below, since the core is an electrical conductor. Historically, solutions to the geomagnetic inverse problem in such a sparsely observed system were thus found only for flow immediately below the core mantle boundary. We have recently shown that combining a numerical model of the geodynamo together with magnetic observations, through the use of Kalman filtering, now allows to present solutions for flow throughout the core. In a second part, I will present synthetic tests of sequential geomagnetic data assimilation aiming at evaluating the range at which the future of the geodynamo can be predicted, and our corresponding prospects to refine the current geomagnetic predictions. Fournier, Aubert, Thébault: Inference on core surface flow from observations and 3-D dynamo modelling, Geophys. J. Int. 186, 118-136, 2011, doi: 10.1111/j.1365-246X.2011.05037.x Aubert, Fournier: Inferring internal properties of Earth's core dynamics and their evolution from surface observations and a numerical geodynamo model, Nonlinear Proc. Geoph. 18, 657-674, 2011, doi:10.5194/npg-18-657-2011 Aubert: Flow throughout the Earth's core inverted from geomagnetic observations and numerical dynamo models, Geophys. J. Int., 2012, doi: 10.1093/gji/ggs051

  4. Nonlocal regularization of inverse problems: a unified variational framework

    PubMed Central

    Yang, Zhili; Jacob, Mathews

    2014-01-01

    We introduce a unifying energy minimization framework for nonlocal regularization of inverse problems. In contrast to the weighted sum of square differences between image pixels used by current schemes, the proposed functional is an unweighted sum of inter-patch distances. We use robust distance metrics that promote the averaging of similar patches, while discouraging the averaging of dissimilar patches. We show that the first iteration of a majorize-minimize algorithm to minimize the proposed cost function is similar to current non-local methods. The reformulation thus provides a theoretical justification for the heuristic approach of iterating non-local schemes, which re-estimate the weights from the current image estimate. Thanks to the reformulation, we now understand that the widely reported alias amplification associated with iterative non-local methods are caused by the convergence to local minimum of the nonconvex penalty. We introduce an efficient continuation strategy to overcome this problem. The similarity of the proposed criterion to widely used non-quadratic penalties (eg. total variation and `p semi-norms) opens the door to the adaptation of fast algorithms developed in the context of compressive sensing; we introduce several novel algorithms to solve the proposed non-local optimization problem. Thanks to the unifying framework, these fast algorithms are readily applicable for a large class of distance metrics. PMID:23014745

  5. Nonlocal regularization of inverse problems: a unified variational framework.

    PubMed

    Yang, Zhili; Jacob, Mathews

    2013-08-01

    We introduce a unifying energy minimization framework for nonlocal regularization of inverse problems. In contrast to the weighted sum of square differences between image pixels used by current schemes, the proposed functional is an unweighted sum of inter-patch distances. We use robust distance metrics that promote the averaging of similar patches, while discouraging the averaging of dissimilar patches. We show that the first iteration of a majorize-minimize algorithm to minimize the proposed cost function is similar to current nonlocal methods. The reformulation thus provides a theoretical justification for the heuristic approach of iterating nonlocal schemes, which re-estimate the weights from the current image estimate. Thanks to the reformulation, we now understand that the widely reported alias amplification associated with iterative nonlocal methods are caused by the convergence to local minimum of the nonconvex penalty. We introduce an efficient continuation strategy to overcome this problem. The similarity of the proposed criterion to widely used nonquadratic penalties (e.g., total variation and lp semi-norms) opens the door to the adaptation of fast algorithms developed in the context of compressive sensing; we introduce several novel algorithms to solve the proposed nonlocal optimization problem. Thanks to the unifying framework, these fast algorithms are readily applicable for a large class of distance metrics. PMID:23014745

  6. Solution of the inverse problem of magnetic induction tomography (MIT).

    PubMed

    Merwa, Robert; Hollaus, Karl; Brunner, Patricia; Scharfetter, Hermann

    2005-04-01

    Magnetic induction tomography (MIT) of biological tissue is used to reconstruct the changes in the complex conductivity distribution inside an object under investigation. The measurement principle is based on determining the perturbation DeltaB of a primary alternating magnetic field B0, which is coupled from an array of excitation coils to the object under investigation. The corresponding voltages DeltaV and V0 induced in a receiver coil carry the information about the passive electrical properties (i.e. conductivity, permittivity and permeability). The reconstruction of the conductivity distribution requires the solution of a 3D inverse eddy current problem. As in EIT the inverse problem is ill-posed and on this account some regularization scheme has to be applied. We developed an inverse solver based on the Gauss-Newton-one-step method for differential imaging, and we implemented and tested four different regularization schemes: the first and second approaches employ a classical smoothness criterion using the unit matrix and a differential matrix of first order as the regularization matrix. The third method is based on variance uniformization, and the fourth method is based on the truncated singular value decomposition. Reconstructions were carried out with synthetic measurement data generated with a spherical perturbation at different locations within a conducting cylinder. Data were generated on a different mesh and 1% random noise was added. The model contained 16 excitation coils and 32 receiver coils which could be combined pairwise to give 16 planar gradiometers. With 32 receiver coils all regularization methods yield fairly good 3D-images of the modelled changes of the conductivity distribution, and prove the feasibility of difference imaging with MIT. The reconstructed perturbations appear at the right location, and their size is in the expected range. With 16 planar gradiometers an additional spurious feature appears mirrored with respect to the median

  7. The geometry of discombinations and its applications to semi-inverse problems in anelasticity

    PubMed Central

    Yavari, Arash; Goriely, Alain

    2014-01-01

    The geometrical formulation of continuum mechanics provides us with a powerful approach to understand and solve problems in anelasticity where an elastic deformation is combined with a non-elastic component arising from defects, thermal stresses, growth effects or other effects leading to residual stresses. The central idea is to assume that the material manifold, prescribing the reference configuration for a body, has an intrinsic, non-Euclidean, geometrical structure. Residual stresses then naturally arise when this configuration is mapped into Euclidean space. Here, we consider the problem of discombinations (a new term that we introduce in this paper), that is, a combined distribution of fields of dislocations, disclinations and point defects. Given a discombination, we compute the geometrical characteristics of the material manifold (curvature, torsion, non-metricity), its Cartan's moving frames and structural equations. This identification provides a powerful algorithm to solve semi-inverse problems with non-elastic components. As an example, we calculate the residual stress field of a cylindrically symmetric distribution of discombinations in an infinite circular cylindrical bar made of an incompressible hyperelastic isotropic elastic solid. PMID:25197257

  8. The geometry of discombinations and its applications to semi-inverse problems in anelasticity.

    PubMed

    Yavari, Arash; Goriely, Alain

    2014-09-01

    The geometrical formulation of continuum mechanics provides us with a powerful approach to understand and solve problems in anelasticity where an elastic deformation is combined with a non-elastic component arising from defects, thermal stresses, growth effects or other effects leading to residual stresses. The central idea is to assume that the material manifold, prescribing the reference configuration for a body, has an intrinsic, non-Euclidean, geometrical structure. Residual stresses then naturally arise when this configuration is mapped into Euclidean space. Here, we consider the problem of discombinations (a new term that we introduce in this paper), that is, a combined distribution of fields of dislocations, disclinations and point defects. Given a discombination, we compute the geometrical characteristics of the material manifold (curvature, torsion, non-metricity), its Cartan's moving frames and structural equations. This identification provides a powerful algorithm to solve semi-inverse problems with non-elastic components. As an example, we calculate the residual stress field of a cylindrically symmetric distribution of discombinations in an infinite circular cylindrical bar made of an incompressible hyperelastic isotropic elastic solid. PMID:25197257

  9. A regularizing iterative ensemble Kalman method for PDE-constrained inverse problems

    NASA Astrophysics Data System (ADS)

    Iglesias, Marco A.

    2016-02-01

    . The numerical investigation is carried out with synthetic experiments on two model inverse problems: (i) identification of conductivity on a Darcy flow model and (ii) electrical impedance tomography with the complete electrode model. We further demonstrate the potential application of the method in solving shape identification problems that arises from the aforementioned forward models by means of a level-set approach for the parameterization of unknown geometries.

  10. Methodes entropiques appliquees au probleme inverse en magnetoencephalographie

    NASA Astrophysics Data System (ADS)

    Lapalme, Ervig

    2005-07-01

    This thesis is devoted to biomagnetic source localization using magnetoencephalography. This problem is known to have an infinite number of solutions. So methods are required to take into account anatomical and functional information on the solution. The work presented in this thesis uses the maximum entropy on the mean method to constrain the solution. This method originates from statistical mechanics and information theory. This thesis is divided into two main parts containing three chapters each. The first part reviews the magnetoencephalographic inverse problem: the theory needed to understand its context and the hypotheses for simplifying the problem. In the last chapter of this first part, the maximum entropy on the mean method is presented: its origins are explained and also how it is applied to our problem. The second part is the original work of this thesis presenting three articles; one of them already published and two others submitted for publication. In the first article, a biomagnetic source model is developed and applied in a theoretical con text but still demonstrating the efficiency of the method. In the second article, we go one step further towards a realistic modelization of the cerebral activation. The main priors are estimated using the magnetoencephalographic data. This method proved to be very efficient in realistic simulations. In the third article, the previous method is extended to deal with time signals thus exploiting the excellent time resolution offered by magnetoencephalography. Compared with our previous work, the temporal method is applied to real magnetoencephalographic data coming from a somatotopy experience and results agree with previous physiological knowledge about this kind of cognitive process.

  11. Subspace-based analysis of the ERT inverse problem

    NASA Astrophysics Data System (ADS)

    Ben Hadj Miled, Mohamed Khames; Miller, Eric L.

    2004-05-01

    In a previous work, we proposed a source-type formulation to the electrical resistance tomography (ERT) problem. Specifically, we showed that inhomogeneities in the medium can be viewed as secondary sources embedded in the homogeneous background medium and located at positions associated with variation in electrical conductivity. Assuming a piecewise constant conductivity distribution, the support of equivalent sources is equal to the boundary of the inhomogeneity. The estimation of the anomaly shape takes the form of an inverse source-type problem. In this paper, we explore the use of subspace methods to localize the secondary equivalent sources associated with discontinuities in the conductivity distribution. Our first alternative is the multiple signal classification (MUSIC) algorithm which is commonly used in the localization of multiple sources. The idea is to project a finite collection of plausible pole (or dipole) sources onto an estimated signal subspace and select those with largest correlations. In ERT, secondary sources are excited simultaneously but in different ways, i.e. with distinct amplitude patterns, depending on the locations and amplitudes of primary sources. If the number of receivers is "large enough", different source configurations can lead to a set of observation vectors that span the data subspace. However, since sources that are spatially close to each other have highly correlated signatures, seperation of such signals becomes very difficult in the presence of noise. To overcome this problem we consider iterative MUSIC algorithms like R-MUSIC and RAP-MUSIC. These recursive algorithms pose a computational burden as they require multiple large combinatorial searches. Results obtained with these algorithms using simulated data of different conductivity patterns are presented.

  12. Final Technical Report for "Applied Mathematics Research: Simulation Based Optimization and Application to Electromagnetic Inverse Problems"

    SciTech Connect

    Haber, Eldad

    2014-03-17

    The focus of research was: Developing adaptive mesh for the solution of Maxwell's equations; Developing a parallel framework for time dependent inverse Maxwell's equations; Developing multilevel methods for optimization problems with inequal- ity constraints; A new inversion code for inverse Maxwell's equations in the 0th frequency (DC resistivity); A new inversion code for inverse Maxwell's equations in low frequency regime. Although the research concentrated on electromagnetic forward and in- verse problems the results of the research was applied to the problem of image registration.

  13. Digital image correlation involves an inverse problem: A regularization scheme based on subset size constraint

    NASA Astrophysics Data System (ADS)

    Zhan, Qin; Yuan, Yuan; Fan, Xiangtao; Huang, Jianyong; Xiong, Chunyang; Yuan, Fan

    2016-06-01

    Digital image correlation (DIC) is essentially implicated in a class of inverse problem. Here, a regularization scheme is developed for the subset-based DIC technique to effectively inhibit potential ill-posedness that likely arises in actual deformation calculations and hence enhance numerical stability, accuracy and precision of correlation measurement. With the aid of a parameterized two-dimensional Butterworth window, a regularized subpixel registration strategy is established, in which the amount of speckle information introduced to correlation calculations may be weighted through equivalent subset size constraint. The optimal regularization parameter associated with each individual sampling point is determined in a self-adaptive way by numerically investigating the curve of 2-norm condition number of coefficient matrix versus the corresponding equivalent subset size, based on which the regularized solution can eventually be obtained. Numerical results deriving from both synthetic speckle images and actual experimental images demonstrate the feasibility and effectiveness of the set of newly-proposed regularized DIC algorithms.

  14. Efficient Calculation of Jacobian and Adjoint Vector Products in the Wave Propagational Inverse Problem Using Automatic Differentiation

    NASA Astrophysics Data System (ADS)

    Coleman, Thomas F.; Santosa, Fadil; Verma, Arun

    2000-01-01

    Wave propagational inverse problems arise in several applications including medical imaging and geophysical exploration. In these problems, one is interested in obtaining the parameters describing the medium from its response to excitations. The problems are characterized by their large size, and by the hyperbolic equation which models the physical phenomena. The inverse problems are often posed as a nonlinear data-fitting where the unknown parameters are found by minimizing the misfit between the predicted data and the actual data. In order to solve the problem numerically using a gradient-type approach, one must calculate the action of the Jacobian and its adjoint on a given vector. In this paper, we explore the use of automatic differentiation (AD) to develop codes that perform these calculations. We show that by exploiting structure at 2 scales, we can arrive at a very efficient code whose main components are produced by AD. In the first scale we exploite the time-stepping nature of the hyperbolic solver by using the “Extended Jacobian” framework. In the second (finer) scale, we exploit the finite difference stencil in order to make explicit use of the sparsity in the dependence of the output variables to the input variables. The main ideas in this work are illustrated with a simpler, one-dimensional version of the problem. Numerical results are given for both one- and two- dimensional problems. We present computational templates that can be used in conjunction with optimization packages to solve the inverse problem.

  15. Inverse problem for the current loop model: Possibilities and restrictions

    NASA Astrophysics Data System (ADS)

    Demina, I. M.; Farafonova, Yu. G.

    2016-07-01

    The possibilities of determining arbitrary current loop parameters based on the spatial structures of the magnetic field components generated by this loop on a sphere with a specified radius have been considered with the use of models. The model parameters were selected such that anomalies created by current loops on a sphere with a radius of 6378 km would be comparable in value with the different-scale anomalies of the observed main geomagnetic field (MGF). The least squares method was used to solve the inverse problem. Estimates close to the specified values were obtained for all current loop parameters except the current strength and radius. The radius determination error can reach ±120 km; at the same time, the magnetic moment value is determined with an accuracy of ±1%. The resolvability of the current force and radius can to a certain degree be improved by decreasing the observation sphere radius such that the ratio of the source distance to the current loop radius would be at least smaller than eight, which can be difficult to reach when modeling MGF.

  16. Toward the solution of the inverse problem in neutron reflectometry

    SciTech Connect

    Haan, V.O. de; Well, A.A. van; Sacks, P.E.; Adenwalla, S.; Felcher, G.P.

    1995-08-01

    The authors show that the chemical depth profile of a film of unknown structure can be retrieved unambiguously from neutron reflection data by adding to the system a known magnetic layer. Three independent reflectivities are obtained by taking measurements with the sample magnetized in a magnetic field perpendicular to the surface and subsequently parallel to it, and using in the latter geometry neutrons polarized either in the direction of the field or opposite to it. The procedure consists of two steps. First, from the three reflectivities both the real and imaginary parts of the reflection coefficient of the unknown film are extracted within the framework of the rigorous dynamical theory. Second, the neutron scattering-length density (and consequently the chemical depth profile) is obtained by a suitable numerical technique for the conventional Schroedinger inverse scattering problem. Computer experiments were conducted for selected cases: starting from the profiles the reflectivities were calculated in a limited range of q and then the original profiles were successfully recovered. The influence on the accuracy of the recovered depth profile of the counting statistics and the cutoffs at low and high q are discussed.

  17. Solving Inverse Detection Problems Using Passive Radiation Signatures

    SciTech Connect

    Favorite, Jeffrey A.; Armstrong, Jerawan C.; Vaquer, Pablo A.

    2012-08-15

    The ability to reconstruct an unknown radioactive object based on its passive gamma-ray and neutron signatures is very important in homeland security applications. Often in the analysis of unknown radioactive objects, for simplicity or speed or because there is no other information, they are modeled as spherically symmetric regardless of their actual geometry. In these presentation we discuss the accuracy and implications of this approximation for decay gamma rays and for neutron-induced gamma rays. We discuss an extension of spherical raytracing (for uncollided fluxes) that allows it to be used when the exterior shielding is flat or cylindrical. We revisit some early results in boundary perturbation theory, showing that the Roussopolos estimate is the correct one to use when the quantity of interest is the flux or leakage on the boundary. We apply boundary perturbation theory to problems in which spherically symmetric systems are perturbed in asymmetric nonspherical ways. We apply mesh adaptive direct search (MADS) algorithms to object reconstructions. We present a benchmark test set that may be used to quantitatively evaluate inverse detection methods.

  18. Inverse Problem for Harmonic Oscillator Perturbed by Potential, Characterization

    NASA Astrophysics Data System (ADS)

    Chelkak, Dmitri; Kargaev, Pavel; Korotyaev, Evgeni

    Consider the perturbed harmonic oscillator Ty=-y''+x2y+q(x)y in L2(R), where the real potential q belongs to the Hilbert space H={q', xq∈ L2(R)}. The spectrum of T is an increasing sequence of simple eigenvalues λn(q)=1+2n+μn, n >= 0, such that μn--> 0 as n-->∞. Let ψn(x,q) be the corresponding eigenfunctions. Define the norming constants νn(q)=limx↑∞log |ψn (x,q)/ψn (-x,q)|. We show that for some real Hilbert space and some subspace Furthermore, the mapping ψ:q|-->ψ(q)=({λn(q)}0∞, {νn(q)}0∞) is a real analytic isomorphism between H and is the set of all strictly increasing sequences s={sn}0∞ such that The proof is based on nonlinear functional analysis combined with sharp asymptotics of spectral data in the high energy limit for complex potentials. We use ideas from the analysis of the inverse problem for the operator -y''py, p∈ L2(0,1), with Dirichlet boundary conditions on the unit interval. There is no literature about the spaces We obtain their basic properties, using their representation as spaces of analytic functions in the disk.

  19. Sparse and redundant representations for inverse problems and recognition

    NASA Astrophysics Data System (ADS)

    Patel, Vishal M.

    Sparse and redundant representation of data enables the description of signals as linear combinations of a few atoms from a dictionary. In this dissertation, we study applications of sparse and redundant representations in inverse problems and object recognition. Furthermore, we propose two novel imaging modalities based on the recently introduced theory of Compressed Sensing (CS). This dissertation consists of four major parts. In the first part of the dissertation, we study a new type of deconvolution algorithm that is based on estimating the image from a shearlet decomposition. Shearlets provide a multi-directional and multi-scale decomposition that has been mathematically shown to represent distributed discontinuities such as edges better than traditional wavelets. We develop a deconvolution algorithm that allows for the approximation inversion operator to be controlled on a multi-scale and multi-directional basis. Furthermore, we develop a method for the automatic determination of the threshold values for the noise shrinkage for each scale and direction without explicit knowledge of the noise variance using a generalized cross validation method. In the second part of the dissertation, we study a reconstruction method that recovers highly undersampled images assumed to have a sparse representation in a gradient domain by using partial measurement samples that are collected in the Fourier domain. Our method makes use of a robust generalized Poisson solver that greatly aids in achieving a significantly improved performance over similar proposed methods. We will demonstrate by experiments that this new technique is more flexible to work with either random or restricted sampling scenarios better than its competitors. In the third part of the dissertation, we introduce a novel Synthetic Aperture Radar (SAR) imaging modality which can provide a high resolution map of the spatial distribution of targets and terrain using a significantly reduced number of needed

  20. The Inverse Phase Stability Problem as a Constraint Satisfaction Problem: Application to Materials Design

    NASA Astrophysics Data System (ADS)

    Arróyave, R.; Gibbons, S. L.; Galvan, E.; Malak, R. J.

    2016-05-01

    In general, the forward phase stability problem consists of mapping thermodynamic conditions (e.g., composition, temperature, pressure) to corresponding equilibrium states. In this paper, we instead focus on the generalized inverse phase stability problem (GIPSP) that deals with mapping a set of phase constitutions to a set of corresponding thermodynamic conditions. Specifically, we define the GIPSP as mapping of sets of phase constitution definitions in a multidimensional phase constitution search space to corresponding ranges of thermodynamic conditions. Mathematically, the solution to the GIPSP corresponds to all solutions to a continuous constraint satisfaction problem (CCSP). We present novel algorithms combining computational thermodynamics, evolutionary computation, and machine learning to approximate solution sets to the GIPSP as a CCSP. Some preliminary examples demonstrating the algorithms are presented. Moreover, the implications of the proposed framework for the larger problem of materials design are discussed, and future work is suggested.

  1. Definition and solution of a stochastic inverse problem for the Manning’s n parameter field in hydrodynamic models

    PubMed Central

    Butler, T.; Graham, L.; Estep, D.; Westerink, J.J.

    2015-01-01

    The uncertainty in spatially heterogeneous Manning’s n fields is quantified using a novel formulation and numerical solution of stochastic inverse problems for physics-based models. The uncertainty is quantified in terms of a probability measure and the physics-based model considered here is the state-of-the-art ADCIRC model although the presented methodology applies to other hydrodynamic models. An accessible overview of the formulation and solution of the stochastic inverse problem in a mathematically rigorous framework based on measure theory is presented. Technical details that arise in practice by applying the framework to determine the Manning’s n parameter field in a shallow water equation model used for coastal hydrodynamics are presented and an efficient computational algorithm and open source software package are developed. A new notion of “condition” for the stochastic inverse problem is defined and analyzed as it relates to the computation of probabilities. This notion of condition is investigated to determine effective output quantities of interest of maximum water elevations to use for the inverse problem for the Manning’s n parameter and the effect on model predictions is analyzed. PMID:25937695

  2. Forward- vs. Inverse Problems in Modeling Seismic Attenuation

    NASA Astrophysics Data System (ADS)

    Morozov, I. B.

    2015-12-01

    Seismic attenuation is an important property of wave propagation used in numerous applications. However, the attenuation is also a complex phenomenon, and it is important to differentiate between its two typical uses: 1) in forward problems, to model the amplitudes and spectral contents of waves required for hazard assessment and geotechnical engineering, and 2) in inverse problems, to determine the physical properties of the subsurface. In the forward-problem sense, the attenuation is successfully characterized in terms of empirical parameters of geometric spreading, radiation patterns, scattering amplitudes, t-star, alpha, kappa, or Q. Arguably, the predicted energy losses can be correct even if the underlying attenuation model is phenomenological and not sufficiently based on physics. An example of such phenomenological model is the viscoelasticity based on the correspondence principle and the Q-factor assigned to the material. By contrast, when used to invert for in situ material properties, models addressing the specific physics are required. In many studies (including in this session), a Q-factor is interpreted as a property of a point within the subsurface; however this property is only phenomenological and may be physically insufficient or inconsistent. For example, the bulk or shear Q at the same point can be different when evaluated from different wave modes. The cases of frequency-dependent Q are particularly prone of ambiguities such as trade-off with the assumed background geometric spreading. To rigorously characterize the in situ material properties responsible for seismic-wave attenuation, it is insufficient to only focus on the seismic energy loss. Mechanical models of the material need to be considered. Such models can be constructed by using Lagrangian mechanics. These models should likely contain no Q but will be based on parameters of microstructure such as heterogeneity, fractures, or fluids. I illustrate several such models based on viscosity

  3. Application of evolution strategies for the solution of an inverse problem in near-field optics.

    PubMed

    Macías, Demetrio; Vial, Alexandre; Barchiesi, Dominique

    2004-08-01

    We introduce an inversion procedure for the characterization of a nanostructure from near-field intensity data. The method proposed is based on heuristic arguments and makes use of evolution strategies for the solution of the inverse problem as a nonlinear constrained-optimization problem. By means of some examples we illustrate the performance of our inversion method. We also discuss its possibilities and potential applications. PMID:15330475

  4. Entire nodal solutions to the pure critical exponent problem arising from concentration

    NASA Astrophysics Data System (ADS)

    Clapp, Mónica

    2016-09-01

    We obtain new sign changing solutions to the problem We exhibit solutions up to (℘p) which blow up at a single point as p →2*, developing a peak whose asymptotic profile is a rescaling of a nonradial sign changing solution to problem (℘∞). We also obtain existence and multiplicity of sign changing nonradial solutions to the Bahri-Coron problem (℘2*) in annuli.

  5. Inverse coefficient problems for one-dimensional heat transfer with a preservation of medium temperature condition

    NASA Astrophysics Data System (ADS)

    Oralsyn, Gulaym

    2016-08-01

    We study an inverse coefficient problem for a model equation for one-dimensional heat transfer with a preservation of medium temperature. It is needed (together with finding its solution) to find time dependent unknown coefficient of the equation. So, for this inverse problem, existence of an unique generalized solution is proved. The main difficulty of the considered problems is that the eigenfunction system of the corresponding boundary value problems does not have the basis property.

  6. Model error estimation and correction by solving a inverse problem

    NASA Astrophysics Data System (ADS)

    Xue, Haile

    2016-04-01

    Nowadays, the weather forecasts and climate predictions are increasingly relied on numerical models. Yet, errors inevitably exist in model due to the imperfect numeric and parameterizations. From the practical point of view, model correction is an efficient strategy. Despite of the different complexity of forecast error correction algorithms, the general idea is to estimate the forecast errors by considering the NWP as a direct problem. Chou (1974) suggested an alternative view by considering the NWP as an inverse problem. The model error tendency term (ME) due to the model deficiency is assumed as an unknown term in NWP model, which can be discretized into short intervals (for example 6 hour) and considered as a constant or linear form in each interval. Given the past re-analyses and NWP model, the discretized MEs in the past intervals can be solved iteratively as a constant or linear-increased tendency term in each interval. These MEs can be further used as the online corrections. In this study, an iterative method for obtaining the MEs in past intervals was presented, and its convergence had been confirmed with sets of experiments in the global forecast system of the Global and Regional Assimilation and Prediction System (GRAPES-GFS) for July-August (JA) 2009 and January-February (JF) 2010. Then these MEs were used to get online model corretions based of systematic errors of GRAPES-GFS for July 2009 and January 2010. The data sets associated with initial condition and sea surface temperature (SST) used in this study are both based on NCEP final (FNL) data. According to the iterative numerical experiments, the following key conclusions can be drawn:(1) Batches of iteration test results indicated that the hour 6 forecast errors were reduced to 10% of their original value after 20 steps of iteration.(2) By offlinely comparing the error corrections estimated by MEs to the mean forecast errors, the patterns of estimated errors were considered to agree well with those

  7. Restarted inverse Born series for the Schrödinger problem with discrete internal measurements

    NASA Astrophysics Data System (ADS)

    Bardsley, Patrick; Guevara Vasquez, Fernando

    2014-04-01

    Convergence and stability results for the inverse Born series (Moskow and Schotland 2008 Inverse Problems 24 065005) are generalized to mappings between Banach spaces. We show that by restarting the inverse Born series one obtains a class of iterative methods containing the Gauss-Newton and Chebyshev-Halley methods. We use the generalized inverse Born series results to show convergence of the inverse Born series for the Schrödinger problem with discrete internal measurements. In this problem, the Schrödinger potential is to be recovered from a few measurements of solutions to the Schrödinger equation resulting from a few different source terms. An application of this method to a problem related to transient hydraulic tomography is given, where the source terms model injection and measurement wells.

  8. Particle-Filter-Based Multisensor Fusion For Solving Low-Frequency Electromagnetic NDE Inverse Problems

    SciTech Connect

    Khan, T.; Ramuhalli, Pradeep; Dass, Sarat

    2011-06-30

    Flaw profile characterization from NDE measurements is a typical inverse problem. A novel transformation of this inverse problem into a tracking problem, and subsequent application of a sequential Monte Carlo method called particle filtering, has been proposed by the authors in an earlier publication [1]. In this study, the problem of flaw characterization from multi-sensor data is considered. The NDE inverse problem is posed as a statistical inverse problem and particle filtering is modified to handle data from multiple measurement modes. The measurement modes are assumed to be independent of each other with principal component analysis (PCA) used to legitimize the assumption of independence. The proposed particle filter based data fusion algorithm is applied to experimental NDE data to investigate its feasibility.

  9. Inversions

    ERIC Educational Resources Information Center

    Brown, Malcolm

    2009-01-01

    Inversions are fascinating phenomena. They are reversals of the normal or expected order. They occur across a wide variety of contexts. What do inversions have to do with learning spaces? The author suggests that they are a useful metaphor for the process that is unfolding in higher education with respect to education. On the basis of…

  10. A boundary integral method for an inverse problem in thermal imaging

    NASA Technical Reports Server (NTRS)

    Bryan, Kurt

    1992-01-01

    An inverse problem in thermal imaging involving the recovery of a void in a material from its surface temperature response to external heating is examined. Uniqueness and continuous dependence results for the inverse problem are demonstrated, and a numerical method for its solution is developed. This method is based on an optimization approach, coupled with a boundary integral equation formulation of the forward heat conduction problem. Some convergence results for the method are proved, and several examples are presented using computationally generated data.

  11. Inverse problems in the design, modeling and testing of engineering systems

    NASA Technical Reports Server (NTRS)

    Alifanov, Oleg M.

    1991-01-01

    Formulations, classification, areas of application, and approaches to solving different inverse problems are considered for the design of structures, modeling, and experimental data processing. Problems in the practical implementation of theoretical-experimental methods based on solving inverse problems are analyzed in order to identify mathematical models of physical processes, aid in input data preparation for design parameter optimization, help in design parameter optimization itself, and to model experiments, large-scale tests, and real tests of engineering systems.

  12. Stability estimates for solution of Bitsadze-Samarskii type inverse elliptic problem with Dirichlet conditions

    NASA Astrophysics Data System (ADS)

    Ashyralyyev, Charyyar; Akyüz, Gulzipa

    2016-08-01

    In this study, we discuss well-posedness of Bitsadze-Samarskii type inverse elliptic problem with Dirichlet conditions. We establish abstract results on stability and coercive stability estimates for the solution of this inverse problem. Then, the abstract results are applied to three overdetermined problems for the multi-dimensional elliptic equation with different boundary conditions. Stability inequalities for solutions of these applications are obtained.

  13. Reinforcement learning solution for HJB equation arising in constrained optimal control problem.

    PubMed

    Luo, Biao; Wu, Huai-Ning; Huang, Tingwen; Liu, Derong

    2015-11-01

    The constrained optimal control problem depends on the solution of the complicated Hamilton-Jacobi-Bellman equation (HJBE). In this paper, a data-based off-policy reinforcement learning (RL) method is proposed, which learns the solution of the HJBE and the optimal control policy from real system data. One important feature of the off-policy RL is that its policy evaluation can be realized with data generated by other behavior policies, not necessarily the target policy, which solves the insufficient exploration problem. The convergence of the off-policy RL is proved by demonstrating its equivalence to the successive approximation approach. Its implementation procedure is based on the actor-critic neural networks structure, where the function approximation is conducted with linearly independent basis functions. Subsequently, the convergence of the implementation procedure with function approximation is also proved. Finally, its effectiveness is verified through computer simulations. PMID:26356598

  14. Inversion problem for ion-atom differential elastic scattering.

    NASA Technical Reports Server (NTRS)

    Rich, W. G.; Bobbio, S. M.; Champion, R. L.; Doverspike, L. D.

    1971-01-01

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

  15. The incomplete inverse and its applications to the linear least squares problem

    NASA Technical Reports Server (NTRS)

    Morduch, G. E.

    1977-01-01

    A modified matrix product is explained, and it is shown that this product defiles a group whose inverse is called the incomplete inverse. It was proven that the incomplete inverse of an augmented normal matrix includes all the quantities associated with the least squares solution. An answer is provided to the problem that occurs when the data residuals are too large and when insufficient data to justify augmenting the model are available.

  16. A new computational method for ambiguity assessment of solutions to assignment problems arising in target tracking

    NASA Astrophysics Data System (ADS)

    Mont, Alexander D.; Calderon, Christopher P.; Poore, Aubrey B.

    2014-06-01

    We present a new approach to estimating the probability of each association in a 2D assignment problem defined by likelihood ratios. Our method divides the set of feasible hypotheses into clusters, and converts a collection of hypotheses into a collection of clusters containing them, reducing the variance of the estimate. Simulations show that our method often generates substantially more accurate probability estimates in less time than traditional methods. Our method can obtain reasonably accurate probabilities of association based on only the input cost matrix and single best candidate solution, eliminating the need for a K-best solution method or MCMC sampling.

  17. Discontinuous finite volume element discretization for coupled flow-transport problems arising in models of sedimentation

    NASA Astrophysics Data System (ADS)

    Bürger, Raimund; Kumar, Sarvesh; Ruiz-Baier, Ricardo

    2015-10-01

    The sedimentation-consolidation and flow processes of a mixture of small particles dispersed in a viscous fluid at low Reynolds numbers can be described by a nonlinear transport equation for the solids concentration coupled with the Stokes problem written in terms of the mixture flow velocity and the pressure field. Here both the viscosity and the forcing term depend on the local solids concentration. A semi-discrete discontinuous finite volume element (DFVE) scheme is proposed for this model. The numerical method is constructed on a baseline finite element family of linear discontinuous elements for the approximation of velocity components and concentration field, whereas the pressure is approximated by piecewise constant elements. The unique solvability of both the nonlinear continuous problem and the semi-discrete DFVE scheme is discussed, and optimal convergence estimates in several spatial norms are derived. Properties of the model and the predicted space accuracy of the proposed formulation are illustrated by detailed numerical examples, including flows under gravity with changing direction, a secondary settling tank in an axisymmetric setting, and batch sedimentation in a tilted cylindrical vessel.

  18. Combination of the LSQR method and a genetic algorithm for solving the electrocardiography inverse problem

    NASA Astrophysics Data System (ADS)

    Jiang, Mingfeng; Xia, Ling; Shou, Guofa; Tang, Min

    2007-03-01

    Computing epicardial potentials from body surface potentials constitutes one form of ill-posed inverse problem of electrocardiography (ECG). To solve this ECG inverse problem, the Tikhonov regularization and truncated singular-value decomposition (TSVD) methods have been commonly used to overcome the ill-posed property by imposing constraints on the magnitudes or derivatives of the computed epicardial potentials. Such direct regularization methods, however, are impractical when the transfer matrix is large. The least-squares QR (LSQR) method, one of the iterative regularization methods based on Lanczos bidiagonalization and QR factorization, has been shown to be numerically more reliable in various circumstances than the other methods considered. This LSQR method, however, to our knowledge, has not been introduced and investigated for the ECG inverse problem. In this paper, the regularization properties of the Krylov subspace iterative method of LSQR for solving the ECG inverse problem were investigated. Due to the 'semi-convergence' property of the LSQR method, the L-curve method was used to determine the stopping iteration number. The performance of the LSQR method for solving the ECG inverse problem was also evaluated based on a realistic heart-torso model simulation protocol. The results show that the inverse solutions recovered by the LSQR method were more accurate than those recovered by the Tikhonov and TSVD methods. In addition, by combing the LSQR with genetic algorithms (GA), the performance can be improved further. It suggests that their combination may provide a good scheme for solving the ECG inverse problem.

  19. Identification of weakly coupled multiphysics problems. Application to the inverse problem of electrocardiography

    NASA Astrophysics Data System (ADS)

    Corrado, Cesare; Gerbeau, Jean-Frédéric; Moireau, Philippe

    2015-02-01

    This work addresses the inverse problem of electrocardiography from a new perspective, by combining electrical and mechanical measurements. Our strategy relies on the definition of a model of the electromechanical contraction which is registered on ECG data but also on measured mechanical displacements of the heart tissue typically extracted from medical images. In this respect, we establish in this work the convergence of a sequential estimator which combines for such coupled problems various state of the art sequential data assimilation methods in a unified consistent and efficient framework. Indeed, we aggregate a Luenberger observer for the mechanical state and a Reduced-Order Unscented Kalman Filter applied on the parameters to be identified and a POD projection of the electrical state. Then using synthetic data we show the benefits of our approach for the estimation of the electrical state of the ventricles along the heart beat compared with more classical strategies which only consider an electrophysiological model with ECG measurements. Our numerical results actually show that the mechanical measurements improve the identifiability of the electrical problem allowing to reconstruct the electrical state of the coupled system more precisely. Therefore, this work is intended to be a first proof of concept, with theoretical justifications and numerical investigations, of the advantage of using available multi-modal observations for the estimation and identification of an electromechanical model of the heart.

  20. Remark on boundary data for inverse boundary value problems for the Navier-Stokes equations

    NASA Astrophysics Data System (ADS)

    Imanuvilov, O. Yu; Yamamoto, M.

    2015-10-01

    In this note, we prove that for the Navier-Stokes equations, a pair of Dirichlet and Neumann data and pressure uniquely correspond to a pair of Dirichlet data and surface stress on the boundary. Hence the two inverse boundary value problems in Imanuvilov and Yamamoto (2015 Inverse Probl. 31 035004) and Lai et al (Arch. Rational Mech. Anal.) are the same.

  1. New problems arising from old drugs: second-generation effects of acetaminophen.

    PubMed

    Tiegs, Gisa; Karimi, Khalil; Brune, Kay; Arck, Petra

    2014-09-01

    Acetaminophen (APAP)/paracetamol is one of the most commonly used over-the-counter drugs taken worldwide for treatment of pain and fever. Although considered as safe when taken in recommended doses not higher than 4 g/day, APAP overdose is currently the most important cause of acute liver failure (ALF). ALF may require liver transplantation and can be fatal. The reasons for APAP-related ALF are mostly intentional (suicidal) or unintentional overdose. However, results from large scale epidemiological studies provide increasing evidence for second generation effects of APAP, even when taken in pharmacological doses. Most strikingly, APAP medication during pregnancy has been associated with health problems including neurodevelopmental and behavioral disorders such as attention deficit hyperactivity disorder and increase in the risk of wheezing and incidence of asthma among offspring. This article reviews the epidemiological findings and aims to shed light into the molecular and cellular mechanisms responsible for APAP-mediated prenatal risk for asthma. PMID:25075430

  2. Application of spectral Lanczos decomposition method to large scale problems arising geophysics

    SciTech Connect

    Tamarchenko, T.

    1996-12-31

    This paper presents an application of Spectral Lanczos Decomposition Method (SLDM) to numerical modeling of electromagnetic diffusion and elastic waves propagation in inhomogeneous media. SLDM approximates an action of a matrix function as a linear combination of basis vectors in Krylov subspace. I applied the method to model electromagnetic fields in three-dimensions and elastic waves in two dimensions. The finite-difference approximation of the spatial part of differential operator reduces the initial boundary-value problem to a system of ordinary differential equations with respect to time. The solution to this system requires calculating exponential and sine/cosine functions of the stiffness matrices. Large scale numerical examples are in a good agreement with the theoretical error bounds and stability estimates given by Druskin, Knizhnerman, 1987.

  3. Graph coarsening methods for Karush-Kuhn-Tucker matrices arising in orthogonal collocation of optimal control problems

    NASA Astrophysics Data System (ADS)

    Senses, Begum

    A state-defect constraint pairing graph coarsening method is described for improving computational efficiency during the numerical factorization of large sparse Karush-Kuhn-Tucker matrices that arise from the discretization of optimal control problems via a Legendre-Gauss-Radau orthogonal collocation method. The method takes advantage of the particular sparse structure of the Karush-Kuhn-Tucker matrix that arises from the orthogonal collocation method. The state-defect constraint pairing graph coarsening method pairs each component of the state with its corresponding defect constraint and forces paired rows to be adjacent in the reordered Karush-Kuhn-Tucker matrix. Aggregate state-defect constraint pairing results are presented using a wide variety of benchmark optimal control problems where it is found that the proposed state-defect constraint pairing graph coarsening method significantly reduces both the number of delayed pivots and the number of floating point operations and increases the computational efficiency by performing more floating point operations per unit time. It is then shown that the state-defect constraint pairing graph coarsening method is less effective on Karush-Kuhn-Tucker matrices arising from Legendre-Gauss-Radau collocation when the optimal control problem contains state and control equality path constraints because such matrices may have delayed pivots that correspond to both defect and path constraints. An unweighted alternate graph coarsening method that employs maximal matching and a weighted alternate graph coarsening method that employs Hungarian algorithm on a weighting matrix are then used to attempt to further reduce the number of delayed pivots. It is found, however, that these alternate graph coarsening methods provide no further advantage over the state-defect constraint pairing graph coarsening method.

  4. A parameter identification problem arising from a two-dimensional airfoil section model

    SciTech Connect

    Cerezo, G.M.

    1994-12-31

    The development of state space models for aeroelastic systems, including unsteady aerodynamics, is particularly important for the design of highly maneuverable aircraft. In this work we present a state space formulation for a special class of singular neutral functional differential equations (SNFDE) with initial data in C(-1, 0). This work is motivated by the two-dimensional airfoil model presented by Burns, Cliff and Herdman in. In the same authors discuss the validity of the assumptions under which the model was formulated. They pay special attention to the derivation of the evolution equation for the circulation on the airfoil. This equation was coupled to the rigid-body dynamics of the airfoil in order to obtain a complete set of functional differential equations that describes the composite system. The resulting mathematical model for the aeroelastic system has a weakly singular component. In this work we consider a finite delay approximation to the model presented in. We work with a scalar model in which we consider the weak singularity appearing in the original problem. The main goal of this work is to develop numerical techniques for the identification of the parameters appearing in the kernel of the associated scalar integral equation. Clearly this is the first step in the study of parameter identification for the original model and the corresponding validation of this model for the aeroelastic system.

  5. PREFACE: The Second International Conference on Inverse Problems: Recent Theoretical Developments and Numerical Approaches

    NASA Astrophysics Data System (ADS)

    Cheng, Jin; Hon, Yiu-Chung; Seo, Jin Keun; Yamamoto, Masahiro

    2005-01-01

    The Second International Conference on Inverse Problems: Recent Theoretical Developments and Numerical Approaches was held at Fudan University, Shanghai from 16-21 June 2004. The first conference in this series was held at the City University of Hong Kong in January 2002 and it was agreed to hold the conference once every two years in a Pan-Pacific Asian country. The next conference is scheduled to be held at Hokkaido University, Sapporo, Japan in July 2006. The purpose of this series of biennial conferences is to establish and develop constant international collaboration, especially among the Pan-Pacific Asian countries. In recent decades, interest in inverse problems has been flourishing all over the globe because of both the theoretical interest and practical requirements. In particular, in Asian countries, one is witnessing remarkable new trends of research in inverse problems as well as the participation of many young talents. Considering these trends, the second conference was organized with the chairperson Professor Li Tat-tsien (Fudan University), in order to provide forums for developing research cooperation and to promote activities in the field of inverse problems. Because solutions to inverse problems are needed in various applied fields, we entertained a total of 92 participants at the second conference and arranged various talks which ranged from mathematical analyses to solutions of concrete inverse problems in the real world. This volume contains 18 selected papers, all of which have undergone peer review. The 18 papers are classified as follows: Surveys: four papers give reviews of specific inverse problems. Theoretical aspects: six papers investigate the uniqueness, stability, and reconstruction schemes. Numerical methods: four papers devise new numerical methods and their applications to inverse problems. Solutions to applied inverse problems: four papers discuss concrete inverse problems such as scattering problems and inverse problems in

  6. Comparing hard and soft prior bounds in geophysical inverse problems

    NASA Technical Reports Server (NTRS)

    Backus, George E.

    1988-01-01

    In linear inversion of a finite-dimensional data vector y to estimate a finite-dimensional prediction vector z, prior information about X sub E is essential if y is to supply useful limits for z. The one exception occurs when all the prediction functionals are linear combinations of the data functionals. Two forms of prior information are compared: a soft bound on X sub E is a probability distribution p sub x on X which describes the observer's opinion about where X sub E is likely to be in X; a hard bound on X sub E is an inequality Q sub x(X sub E, X sub E) is equal to or less than 1, where Q sub x is a positive definite quadratic form on X. A hard bound Q sub x can be softened to many different probability distributions p sub x, but all these p sub x's carry much new information about X sub E which is absent from Q sub x, and some information which contradicts Q sub x. Both stochastic inversion (SI) and Bayesian inference (BI) estimate z from y and a soft prior bound p sub x. If that probability distribution was obtained by softening a hard prior bound Q sub x, rather than by objective statistical inference independent of y, then p sub x contains so much unsupported new information absent from Q sub x that conclusions about z obtained with SI or BI would seen to be suspect.

  7. Comparing hard and soft prior bounds in geophysical inverse problems

    NASA Technical Reports Server (NTRS)

    Backus, George E.

    1987-01-01

    In linear inversion of a finite-dimensional data vector y to estimate a finite-dimensional prediction vector z, prior information about X sub E is essential if y is to supply useful limits for z. The one exception occurs when all the prediction functionals are linear combinations of the data functionals. Two forms of prior information are compared: a soft bound on X sub E is a probability distribution p sub x on X which describeds the observer's opinion about where X sub E is likely to be in X; a hard bound on X sub E is an inequality Q sub x(X sub E, X sub E) is equal to or less than 1, where Q sub x is a positive definite quadratic form on X. A hard bound Q sub x can be softened to many different probability distributions p sub x, but all these p sub x's carry much new information about X sub E which is absent from Q sub x, and some information which contradicts Q sub x. Both stochastic inversion (SI) and Bayesian inference (BI) estimate z from y and a soft prior bound p sub x. If that probability distribution was obtained by softening a hard prior bound Q sub x, rather than by objective statistical inference independent of y, then p sub x contains so much unsupported new information absent from Q sub x that conclusions about z obtained with SI or BI would seen to be suspect.

  8. Global solution to a hyperbolic problem arising in the modeling of blood flow in circulatory systems

    NASA Astrophysics Data System (ADS)

    Ruan, Weihua; Clark, M. E.; Zhao, Meide; Curcio, Anthony

    2007-07-01

    This paper considers a system of first-order, hyperbolic, partial differential equations in the domain of a one-dimensional network. The system models the blood flow in human circulatory systems as an initial-boundary-value problem with boundary conditions of either algebraic or differential type. The differential equations are nonhomogeneous with frictional damping terms and the state variables are coupled at internal junctions. The existence and uniqueness of the local classical solution have been established in our earlier work [W. Ruan, M.E. Clark, M. Zhao, A. Curcio, A hyperbolic system of equations of blood flow in an arterial network, J. Appl. Math. 64 (2) (2003) 637-667; W. Ruan, M.E. Clark, M. Zhao, A. Curcio, Blood flow in a network, Nonlinear Anal. Real World Appl. 5 (2004) 463-485; W. Ruan, M.E. Clark, M. Zhao, A. Curcio, A quasilinear hyperbolic system that models blood flow in a network, in: Charles V. Benton (Ed.), Focus on Mathematical Physics Research, Nova Science Publishers, Inc., New York, 2004, pp. 203-230]. This paper continues the analysis and gives sufficient conditions for the global existence of the classical solution. We prove that the solution exists globally if the boundary data satisfy the dissipative condition (2.3) or (3.2), and the norms of the initial and forcing functions in a certain Sobolev space are sufficiently small. This is only the first step toward establishing the global existence of the solution to physiologically realistic models, because, in general, the chosen dissipative conditions (2.3) and (3.2) do not appear to hold for the originally proposed boundary conditions (1.3)-(1.12).

  9. Methods of noninvasive electrophysiological heart examination basing on solution of inverse problem of electrocardiography

    NASA Astrophysics Data System (ADS)

    Grigoriev, M.; Babich, L.

    2015-09-01

    The article represents the main noninvasive methods of heart electrical activity examination, theoretical bases of solution of electrocardiography inverse problem, application of different methods of heart examination in clinical practice, and generalized achievements in this sphere in global experience.

  10. Numerical solution of an inverse electrocardiography problem for a medium with piecewise constant electrical conductivity

    NASA Astrophysics Data System (ADS)

    Denisov, A. M.; Zakharov, E. V.; Kalinin, A. V.; Kalinin, V. V.

    2010-07-01

    A numerical method is proposed for solving an inverse electrocardiography problem for a medium with a piecewise constant electrical conductivity. The method is based on the method of boundary integral equations and Tikhonov regularization.

  11. Average synaptic activity and neural networks topology: a global inverse problem

    NASA Astrophysics Data System (ADS)

    Burioni, Raffaella; Casartelli, Mario; di Volo, Matteo; Livi, Roberto; Vezzani, Alessandro

    2014-03-01

    The dynamics of neural networks is often characterized by collective behavior and quasi-synchronous events, where a large fraction of neurons fire in short time intervals, separated by uncorrelated firing activity. These global temporal signals are crucial for brain functioning. They strongly depend on the topology of the network and on the fluctuations of the connectivity. We propose a heterogeneous mean-field approach to neural dynamics on random networks, that explicitly preserves the disorder in the topology at growing network sizes, and leads to a set of self-consistent equations. Within this approach, we provide an effective description of microscopic and large scale temporal signals in a leaky integrate-and-fire model with short term plasticity, where quasi-synchronous events arise. Our equations provide a clear analytical picture of the dynamics, evidencing the contributions of both periodic (locked) and aperiodic (unlocked) neurons to the measurable average signal. In particular, we formulate and solve a global inverse problem of reconstructing the in-degree distribution from the knowledge of the average activity field. Our method is very general and applies to a large class of dynamical models on dense random networks.

  12. Some Inverse Problems in Periodic Homogenization of Hamilton-Jacobi Equations

    NASA Astrophysics Data System (ADS)

    Luo, Songting; Tran, Hung V.; Yu, Yifeng

    2016-03-01

    We look at the effective Hamiltonian {overline{H}} associated with the Hamiltonian {H(p,x)=H(p)+V(x)} in the periodic homogenization theory. Our central goal is to understand the relation between {V} and {overline{H}} . We formulate some inverse problems concerning this relation. Such types of inverse problems are, in general, very challenging. In this paper, we discuss several special cases in both convex and nonconvex settings.

  13. Review of the inverse scattering problem at fixed energy in quantum mechanics

    NASA Technical Reports Server (NTRS)

    Sabatier, P. C.

    1972-01-01

    Methods of solution of the inverse scattering problem at fixed energy in quantum mechanics are presented. Scattering experiments of a beam of particles at a nonrelativisitic energy by a target made up of particles are analyzed. The Schroedinger equation is used to develop the quantum mechanical description of the system and one of several functions depending on the relative distance of the particles. The inverse problem is the construction of the potentials from experimental measurements.

  14. Some Inverse Problems in Periodic Homogenization of Hamilton-Jacobi Equations

    NASA Astrophysics Data System (ADS)

    Luo, Songting; Tran, Hung V.; Yu, Yifeng

    2016-09-01

    We look at the effective Hamiltonian {overline{H}} associated with the Hamiltonian {H(p,x)=H(p)+V(x)} in the periodic homogenization theory. Our central goal is to understand the relation between {V} and {overline{H}}. We formulate some inverse problems concerning this relation. Such types of inverse problems are, in general, very challenging. In this paper, we discuss several special cases in both convex and nonconvex settings.

  15. Sedimentary Facies Analysis Using AVIRIS Data: A Geophysical Inverse Problem

    NASA Technical Reports Server (NTRS)

    Boardmann, Joe W.; Goetz, Alexander F. H.

    1990-01-01

    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.

  16. On solvability of inverse coefficient problems for nonlinear convection—diffusion—reaction equation

    NASA Astrophysics Data System (ADS)

    Brizitskii, R. V.; Saritskaya, Zh Yu

    2016-06-01

    An inverse coefficient problem is considered for a stationary nonlinear convection- diffusion-reaction equation, in which reaction coefficient has a rather common dependence on substance concentration and on spacial variable. The solvability of the considered nonlinear boundary value problem is proved in a general case. The existence of solutions of the inverse problem is proved for the reaction coefficients, which are defined by the product of two functions. The first function depends on a spatial variable, the second one depends nonlinearly on the solution of the boundary value problem. The mentioned inverse problem consists in reconstructing the first function with the help of additional information provided by the solution of the boundary value problem.

  17. FOREWORD: 5th International Workshop on New Computational Methods for Inverse Problems

    NASA Astrophysics Data System (ADS)

    Vourc'h, Eric; Rodet, Thomas

    2015-11-01

    This volume of Journal of Physics: Conference Series is dedicated to the scientific research presented during the 5th International Workshop on New Computational Methods for Inverse Problems, NCMIP 2015 (http://complement.farman.ens-cachan.fr/NCMIP_2015.html). This workshop took place at Ecole Normale Supérieure de Cachan, on May 29, 2015. The prior editions of NCMIP also took place in Cachan, France, firstly within the scope of ValueTools Conference, in May 2011, and secondly at the initiative of Institut Farman, in May 2012, May 2013 and May 2014. The New Computational Methods for Inverse Problems (NCMIP) workshop focused on recent advances in the resolution of inverse problems. Indeed, inverse problems appear in numerous scientific areas such as geophysics, biological and medical imaging, material and structure characterization, electrical, mechanical and civil engineering, and finances. The resolution of inverse problems consists of estimating the parameters of the observed system or structure from data collected by an instrumental sensing or imaging device. Its success firstly requires the collection of relevant observation data. It also requires accurate models describing the physical interactions between the instrumental device and the observed system, as well as the intrinsic properties of the solution itself. Finally, it requires the design of robust, accurate and efficient inversion algorithms. Advanced sensor arrays and imaging devices provide high rate and high volume data; in this context, the efficient resolution of the inverse problem requires the joint development of new models and inversion methods, taking computational and implementation aspects into account. During this one-day workshop, researchers had the opportunity to bring to light and share new techniques and results in the field of inverse problems. The topics of the workshop were: algorithms and computational aspects of inversion, Bayesian estimation, Kernel methods, learning methods

  18. FOREWORD: 4th International Workshop on New Computational Methods for Inverse Problems (NCMIP2014)

    NASA Astrophysics Data System (ADS)

    2014-10-01

    This volume of Journal of Physics: Conference Series is dedicated to the scientific contributions presented during the 4th International Workshop on New Computational Methods for Inverse Problems, NCMIP 2014 (http://www.farman.ens-cachan.fr/NCMIP_2014.html). This workshop took place at Ecole Normale Supérieure de Cachan, on May 23, 2014. The prior editions of NCMIP also took place in Cachan, France, firstly within the scope of ValueTools Conference, in May 2011 (http://www.ncmip.org/2011/), and secondly at the initiative of Institut Farman, in May 2012 and May 2013, (http://www.farman.ens-cachan.fr/NCMIP_2012.html), (http://www.farman.ens-cachan.fr/NCMIP_2013.html). The New Computational Methods for Inverse Problems (NCMIP) Workshop focused on recent advances in the resolution of inverse problems. Indeed, inverse problems appear in numerous scientific areas such as geophysics, biological and medical imaging, material and structure characterization, electrical, mechanical and civil engineering, and finances. The resolution of inverse problems consists of estimating the parameters of the observed system or structure from data collected by an instrumental sensing or imaging device. Its success firstly requires the collection of relevant observation data. It also requires accurate models describing the physical interactions between the instrumental device and the observed system, as well as the intrinsic properties of the solution itself. Finally, it requires the design of robust, accurate and efficient inversion algorithms. Advanced sensor arrays and imaging devices provide high rate and high volume data; in this context, the efficient resolution of the inverse problem requires the joint development of new models and inversion methods, taking computational and implementation aspects into account. During this one-day workshop, researchers had the opportunity to bring to light and share new techniques and results in the field of inverse problems. The topics of the

  19. The inverse problem of refraction travel times, part I: Types of Geophysical Nonuniqueness through Minimization

    USGS Publications Warehouse

    Ivanov, J.; Miller, R.D.; Xia, J.; Steeples, D.; Park, C.B.

    2005-01-01

    In a set of two papers we study the inverse problem of refraction travel times. The purpose of this work is to use the study as a basis for development of more sophisticated methods for finding more reliable solutions to the inverse problem of refraction travel times, which is known to be nonunique. The first paper, "Types of Geophysical Nonuniqueness through Minimization," emphasizes the existence of different forms of nonuniqueness in the realm of inverse geophysical problems. Each type of nonuniqueness requires a different type and amount of a priori information to acquire a reliable solution. Based on such coupling, a nonuniqueness classification is designed. Therefore, since most inverse geophysical problems are nonunique, each inverse problem must be studied to define what type of nonuniqueness it belongs to and thus determine what type of a priori information is necessary to find a realistic solution. The second paper, "Quantifying Refraction Nonuniqueness Using a Three-layer Model," serves as an example of such an approach. However, its main purpose is to provide a better understanding of the inverse refraction problem by studying the type of nonuniqueness it possesses. An approach for obtaining a realistic solution to the inverse refraction problem is planned to be offered in a third paper that is in preparation. The main goal of this paper is to redefine the existing generalized notion of nonuniqueness and a priori information by offering a classified, discriminate structure. Nonuniqueness is often encountered when trying to solve inverse problems. However, possible nonuniqueness diversity is typically neglected and nonuniqueness is regarded as a whole, as an unpleasant "black box" and is approached in the same manner by applying smoothing constraints, damping constraints with respect to the solution increment and, rarely, damping constraints with respect to some sparse reference information about the true parameters. In practice, when solving geophysical

  20. Guest Editors' introduction to the special section on statistical and computational issues in inverse problems

    NASA Astrophysics Data System (ADS)

    Tenorio, L.; Haber, E.; Symes, W. W.; Stark, P. B.; Cox, D.; Ghattas, O.

    2008-06-01

    In the words of D D Jackson, the data of real-world inverse problems tend to be inaccurate, insufficient and inconsistent (1972 Geophys. J. R. Astron. Soc. 28 97-110). In view of these features, the characterization of solution uncertainty is an essential aspect of the study of inverse problems. The development of computational technology, in particular of multiscale and adaptive methods and robust optimization algorithms, has combined with advances in statistical methods in recent years to create unprecedented opportunities to understand and explore the role of uncertainty in inversion. Following this introductory article, the special section contains 16 papers describing recent statistical and computational advances in a variety of inverse problem settings.

  1. Using a derivative-free optimization method for multiple solutions of inverse transport problems

    DOE PAGESBeta

    Armstrong, Jerawan C.; Favorite, Jeffrey A.

    2016-01-14

    Identifying unknown components of an object that emits radiation is an important problem for national and global security. Radiation signatures measured from an object of interest can be used to infer object parameter values that are not known. This problem is called an inverse transport problem. An inverse transport problem may have multiple solutions and the most widely used approach for its solution is an iterative optimization method. This paper proposes a stochastic derivative-free global optimization algorithm to find multiple solutions of inverse transport problems. The algorithm is an extension of a multilevel single linkage (MLSL) method where a meshmore » adaptive direct search (MADS) algorithm is incorporated into the local phase. Furthermore, numerical test cases using uncollided fluxes of discrete gamma-ray lines are presented to show the performance of this new algorithm.« less

  2. An adaptive large neighborhood search heuristic for Two-Echelon Vehicle Routing Problems arising in city logistics

    PubMed Central

    Hemmelmayr, Vera C.; Cordeau, Jean-François; Crainic, Teodor Gabriel

    2012-01-01

    In this paper, we propose an adaptive large neighborhood search heuristic for the Two-Echelon Vehicle Routing Problem (2E-VRP) and the Location Routing Problem (LRP). The 2E-VRP arises in two-level transportation systems such as those encountered in the context of city logistics. In such systems, freight arrives at a major terminal and is shipped through intermediate satellite facilities to the final customers. The LRP can be seen as a special case of the 2E-VRP in which vehicle routing is performed only at the second level. We have developed new neighborhood search operators by exploiting the structure of the two problem classes considered and have also adapted existing operators from the literature. The operators are used in a hierarchical scheme reflecting the multi-level nature of the problem. Computational experiments conducted on several sets of instances from the literature show that our algorithm outperforms existing solution methods for the 2E-VRP and achieves excellent results on the LRP. PMID:23483764

  3. Extremal norms of the potentials recovered from inverse Dirichlet problems

    NASA Astrophysics Data System (ADS)

    Qi, Jiangang; Chen, Shaozhu

    2016-03-01

    Consider the Sturm-Liouville eigenvalue problem -y\\prime\\prime (x)+q(x)y(x)=λ y(x),x\\in [0,1],y(0)=y(1)=0, where q\\in {L}1[0,1], and its spectrum is denoted by σ (q). For a real number λ, define {{Ω }}(λ )=\\{q\\in {L}1[0,1] :λ \\in σ (q)\\} and E(λ )={inf}\\{\\parallel q\\parallel :q\\in {{Ω }}(λ )\\}. We will set up a formula for E(λ ) explicitly in terms of λ and specify where the infimum can be attained. As an application, we will give the extremal values of the nth eigenvalue of the Dirichlet problem for potentials on a sphere {L}1[0,1], n≥slant 1. The proofs are based on a new Lyapunov-type inequality for Sturm-Liouville equations with potentials.

  4. Statistical method for resolving the photon-photoelectron-counting inversion problem

    SciTech Connect

    Wu Jinlong; Li Tiejun; Peng, Xiang; Guo Hong

    2011-02-01

    A statistical inversion method is proposed for the photon-photoelectron-counting statistics in quantum key distribution experiment. With the statistical viewpoint, this problem is equivalent to the parameter estimation for an infinite binomial mixture model. The coarse-graining idea and Bayesian methods are applied to deal with this ill-posed problem, which is a good simple example to show the successful application of the statistical methods to the inverse problem. Numerical results show the applicability of the proposed strategy. The coarse-graining idea for the infinite mixture models should be general to be used in the future.

  5. Unrealistic parameter estimates in inverse modelling: A problem or a benefit for model calibration?

    USGS Publications Warehouse

    Poeter, E.P.; Hill, M.C.

    1996-01-01

    Estimation of unrealistic parameter values by inverse modelling is useful for constructed model discrimination. This utility is demonstrated using the three-dimensional, groundwater flow inverse model MODFLOWP to estimate parameters in a simple synthetic model where the true conditions and character of the errors are completely known. When a poorly constructed model is used, unreasonable parameter values are obtained even when using error free observations and true initial parameter values. This apparent problem is actually a benefit because it differentiates accurately and inaccurately constructed models. The problems seem obvious for a synthetic problem in which the truth is known, but are obscure when working with field data. Situations in which unrealistic parameter estimates indicate constructed model problems are illustrated in applications of inverse modelling to three field sites and to complex synthetic test cases in which it is shown that prediction accuracy also suffers when constructed models are inaccurate.

  6. Geophysical approaches to inverse problems: A methodological comparison. Part 1: A Posteriori approach

    NASA Technical Reports Server (NTRS)

    Seidman, T. I.; Munteanu, M. J.

    1979-01-01

    The relationships of a variety of general computational methods (and variances) for treating illposed problems such as geophysical inverse problems are considered. Differences in approach and interpretation based on varying assumptions as to, e.g., the nature of measurement uncertainties are discussed along with the factors to be considered in selecting an approach. The reliability of the results of such computation is addressed.

  7. Moment inversion problem for piecewise D-finite functions

    NASA Astrophysics Data System (ADS)

    Batenkov, Dmitry

    2009-10-01

    We consider the problem of exact reconstruction of univariate functions with jump discontinuities at unknown positions from their moments. These functions are assumed to satisfy an a priori unknown linear homogeneous differential equation with polynomial coefficients on each continuity interval. Therefore, they may be specified by a finite amount of information. This reconstruction problem has practical importance in signal processing and other applications. It is somewhat of a 'folklore' that the sequence of the moments of such 'piecewise D-finite' functions satisfies a linear recurrence relation of bounded order and degree. We derive this recurrence relation explicitly. It turns out that the coefficients of the differential operator which annihilates every piece of the function, as well as the locations of the discontinuities, appear in this recurrence in a precisely controlled manner. This leads to the formulation of a generic algorithm for reconstructing a piecewise D-finite function from its moments. We investigate the conditions for solvability of resulting linear systems in the general case, as well as analyse a few particular examples. We provide results of numerical simulations for several types of signals, which test the sensitivity of the proposed algorithm to noise.

  8. Some speed-up strategies for solving inverse radiative transfer problems

    NASA Astrophysics Data System (ADS)

    Favennec, Y.; Le Hardy, D.; Dubot, F.; Rousseau, B.; Rousse, D. R.

    2016-01-01

    Inversion based on the radiative transfer equation (RTE) is generally highly CPU time consuming because the forward model itself is complicated to solve when the space dimension is greater than one, and because the inversion is based on a large number of forward model runs until convergence is reached. The goal of this paper is to set up some speed-up strategies specific of inversion when radiative transfer problems are dealt with. More specifically, the accurate identification of the volumetric radiative properties i.e. both the absorption and scattering coefficients is the objective of this study.

  9. From Bayes to Tarantola: New insights to understand uncertainty in inverse problems

    NASA Astrophysics Data System (ADS)

    Fernández-Martínez, J. L.; Fernández-Muñiz, Z.; Pallero, J. L. G.; Pedruelo-González, L. M.

    2013-11-01

    Anyone working on inverse problems is aware of their ill-posed character. In the case of inverse problems, this concept (ill-posed) proposed by J. Hadamard in 1902, admits revision since it is somehow related to their ill-conditioning and the use of local optimization methods to find their solution. A more general and interesting approach regarding risk analysis and epistemological decision making would consist in analyzing the existence of families of equivalent model parameters that are compatible with the prior information and predict the observed data within the same error bounds. Otherwise said, the ill-posed character of discrete inverse problems (ill-conditioning) originates that their solution is uncertain. Traditionally nonlinear inverse problems in discrete form have been solved via local optimization methods with regularization, but linear analysis techniques failed to account for the uncertainty in the solution that it is adopted. As a result of this fact uncertainty analysis in nonlinear inverse problems has been approached in a probabilistic framework (Bayesian approach), but these methods are hindered by the curse of dimensionality and by the high computational cost needed to solve the corresponding forward problems. Global optimization techniques are very attractive, but most of the times are heuristic and have the same limitations than Monte Carlo methods. New research is needed to provide uncertainty estimates, especially in the case of high dimensional nonlinear inverse problems with very costly forward problems. After the discredit of deterministic methods and some initial years of Bayesian fever, now the pendulum seems to return back, because practitioners are aware that the uncertainty analysis in high dimensional nonlinear inverse problems cannot (and should not be) solved via random sampling methodologies. The main reason is that the uncertainty “space” of nonlinear inverse problems has a mathematical structure that is embedded in the

  10. Fast and accurate analytical model to solve inverse problem in SHM using Lamb wave propagation

    NASA Astrophysics Data System (ADS)

    Poddar, Banibrata; Giurgiutiu, Victor

    2016-04-01

    Lamb wave propagation is at the center of attention of researchers for structural health monitoring of thin walled structures. This is due to the fact that Lamb wave modes are natural modes of wave propagation in these structures with long travel distances and without much attenuation. This brings the prospect of monitoring large structure with few sensors/actuators. However the problem of damage detection and identification is an "inverse problem" where we do not have the luxury to know the exact mathematical model of the system. On top of that the problem is more challenging due to the confounding factors of statistical variation of the material and geometric properties. Typically this problem may also be ill posed. Due to all these complexities the direct solution of the problem of damage detection and identification in SHM is impossible. Therefore an indirect method using the solution of the "forward problem" is popular for solving the "inverse problem". This requires a fast forward problem solver. Due to the complexities involved with the forward problem of scattering of Lamb waves from damages researchers rely primarily on numerical techniques such as FEM, BEM, etc. But these methods are slow and practically impossible to be used in structural health monitoring. We have developed a fast and accurate analytical forward problem solver for this purpose. This solver, CMEP (complex modes expansion and vector projection), can simulate scattering of Lamb waves from all types of damages in thin walled structures fast and accurately to assist the inverse problem solver.

  11. Observation and inverse problems in coupled cell networks

    NASA Astrophysics Data System (ADS)

    Joly, Romain

    2012-03-01

    A coupled cell network is a model for many situations such as food webs in ecosystems, cellular metabolism and economic networks. It consists in a directed graph G, each node (or cell) representing an agent of the network and each directed arrow representing which agent acts on which. It yields a system of differential equations \\dot x(t)=f(x(t)) , where the component i of f depends only on the cells xj(t) for which the arrow j → i exists in G. In this paper, we investigate the observation problems in coupled cell networks: can one deduce the behaviour of the whole network (oscillations, stabilization, etc) by observing only one of the cells? We show that the natural observation properties hold for almost all the interactions f.

  12. Analytical and numerical analysis of inverse optimization problems: conditions of uniqueness and computational methods

    PubMed Central

    Zatsiorsky, Vladimir M.

    2011-01-01

    One of the key problems of motor control is the redundancy problem, in particular how the central nervous system (CNS) chooses an action out of infinitely many possible. A promising way to address this question is to assume that the choice is made based on optimization of a certain cost function. A number of cost functions have been proposed in the literature to explain performance in different motor tasks: from force sharing in grasping to path planning in walking. However, the problem of uniqueness of the cost function(s) was not addressed until recently. In this article, we analyze two methods of finding additive cost functions in inverse optimization problems with linear constraints, so-called linear-additive inverse optimization problems. These methods are based on the Uniqueness Theorem for inverse optimization problems that we proved recently (Terekhov et al., J Math Biol 61(3):423–453, 2010). Using synthetic data, we show that both methods allow for determining the cost function. We analyze the influence of noise on the both methods. Finally, we show how a violation of the conditions of the Uniqueness Theorem may lead to incorrect solutions of the inverse optimization problem. PMID:21311907

  13. Analytical and numerical analysis of inverse optimization problems: conditions of uniqueness and computational methods.

    PubMed

    Terekhov, Alexander V; Zatsiorsky, Vladimir M

    2011-02-01

    One of the key problems of motor control is the redundancy problem, in particular how the central nervous system (CNS) chooses an action out of infinitely many possible. A promising way to address this question is to assume that the choice is made based on optimization of a certain cost function. A number of cost functions have been proposed in the literature to explain performance in different motor tasks: from force sharing in grasping to path planning in walking. However, the problem of uniqueness of the cost function(s) was not addressed until recently. In this article, we analyze two methods of finding additive cost functions in inverse optimization problems with linear constraints, so-called linear-additive inverse optimization problems. These methods are based on the Uniqueness Theorem for inverse optimization problems that we proved recently (Terekhov et al., J Math Biol 61(3):423-453, 2010). Using synthetic data, we show that both methods allow for determining the cost function. We analyze the influence of noise on the both methods. Finally, we show how a violation of the conditions of the Uniqueness Theorem may lead to incorrect solutions of the inverse optimization problem. PMID:21311907

  14. Reification of galaxies: cognitive astrophysics and the multiwavelength inverse problem

    NASA Astrophysics Data System (ADS)

    Madore, Barry F.

    2012-08-01

    Lessons learned in the history and philosophy of science have generally had little immediate impact on how we as individual astronomers conduct our research. And yet we do share many common views on how we undertake basic research, and how we translate observations and theory into communicable knowledge. In this introductory talk I will illustrate how we as extragalactic astronomers have already violated some of the basic tenets of what constitutes ``science'' as seen from a philosophical point of view, and I will predict what the future of astronomy as a science may soon look like. Simple examples of how we are already ``cognitively closed'' to many immediate and tangible aspects of the Universe will be given and some solutions to this dilemma will be proposed. We may be at a point in time where more data is not necessarily the best solution to our problems. Discovering that familiar concepts and even certain objects may not exist in the traditional sense of the word could provide a motivation for broadening our way of conceptualizing the extragalactic Universe, more as a continuum of processes and phase transitions rather than an assembly of discrete objects. Once again the Universe may be ``forcing us to think''.

  15. Removal of numerical instability in the solution of an inverse heat conduction problem

    NASA Astrophysics Data System (ADS)

    Pourgholi, R.; Azizi, N.; Gasimov, Y. S.; Aliev, F.; Khalafi, H. K.

    2009-06-01

    In this paper, we consider an inverse heat conduction problem (IHCP). A set of temperature measurements at a single sensor location inside the heat conduction body is required. Using a transformation, the ill-posed IHCP becomes a Cauchy problem. Since the solution of Cauchy problem, exists and is unique but not always stable, the ill-posed problem is closely approximated by a well-posed problem. For this new well-posed problem, the existence, uniqueness, and stability of the solution are proved.

  16. Analysis and algorithms for a regularized Cauchy problem arising from a non-linear elliptic PDE for seismic velocity estimation

    SciTech Connect

    Cameron, M.K.; Fomel, S.B.; Sethian, J.A.

    2009-01-01

    In the present work we derive and study a nonlinear elliptic PDE coming from the problem of estimation of sound speed inside the Earth. The physical setting of the PDE allows us to pose only a Cauchy problem, and hence is ill-posed. However we are still able to solve it numerically on a long enough time interval to be of practical use. We used two approaches. The first approach is a finite difference time-marching numerical scheme inspired by the Lax-Friedrichs method. The key features of this scheme is the Lax-Friedrichs averaging and the wide stencil in space. The second approach is a spectral Chebyshev method with truncated series. We show that our schemes work because of (1) the special input corresponding to a positive finite seismic velocity, (2) special initial conditions corresponding to the image rays, (3) the fact that our finite-difference scheme contains small error terms which damp the high harmonics; truncation of the Chebyshev series, and (4) the need to compute the solution only for a short interval of time. We test our numerical scheme on a collection of analytic examples and demonstrate a dramatic improvement in accuracy in the estimation of the sound speed inside the Earth in comparison with the conventional Dix inversion. Our test on the Marmousi example confirms the effectiveness of the proposed approach.

  17. Interlocked optimization and fast gradient algorithm for a seismic inverse problem

    SciTech Connect

    Metivier, Ludovic

    2011-08-10

    Highlights: {yields} A 2D extension of the 1D nonlinear inversion of well-seismic data is given. {yields} Appropriate regularization yields a well-determined large scale inverse problem. {yields} An interlocked optimization loop acts as an efficient preconditioner. {yields} The adjoint state method is used to compute the misfit function gradient. {yields} Domain decomposition method yields an efficient parallel implementation. - Abstract: We give a nonlinear inverse method for seismic data recorded in a well from sources at several offsets from the borehole in a 2D acoustic framework. Given the velocity field, approximate values of the impedance are recovered. This is a 2D extension of the 1D inversion of vertical seismic profiles . The inverse problem generates a large scale undetermined ill-conditioned problem. Appropriate regularization terms render the problem well-determined. An interlocked optimization algorithm yields an efficient preconditioning. A gradient algorithm based on the adjoint state method and domain decomposition gives a fast parallel numerical method. For a realistic test case, convergence is attained in an acceptable time with 128 processors.

  18. Sensitivity analysis of solutions of the harmonic inversion problem: are all data points created equal?

    PubMed

    Celik, Hasan; Shaka, A J; Mandelshtam, V A

    2010-09-01

    We consider the harmonic inversion problem, and the associated spectral estimation problem, both of which are key numerical problems in NMR data analysis. Under certain conditions (in particular, in exact arithmetic) these problems have unique solutions. Therefore, these solutions must not depend on the inversion algorithm, as long as it is exact in principle. In this paper, we are not concerned with the algorithmic aspects of harmonic inversion, but rather with the sensitivity of the solutions of the problem to perturbations of the time-domain data. A sensitivity analysis was performed and the counterintuitive results call into question the common assumption made in "super-resolution" methods using non-uniform data sampling, namely, that the data should be sampled more often where the time signal has the largest signal-to-noise ratio. The numerical analysis herein demonstrates that the spectral parameters (such as the peak positions and amplitudes) resulting from the solution of the harmonic inversion problem are least susceptible to the perturbations in the values of data points at the edges of the time interval and most susceptible to the perturbations in the values at intermediate times. PMID:20663693

  19. The inverse problem of refraction travel times, part II: Quantifying refraction nonuniqueness using a three-layer model

    USGS Publications Warehouse

    Ivanov, J.; Miller, R.D.; Xia, J.; Steeples, D.

    2005-01-01

    This paper is the second of a set of two papers in which we study the inverse refraction problem. The first paper, "Types of Geophysical Nonuniqueness through Minimization," studies and classifies the types of nonuniqueness that exist when solving inverse problems depending on the participation of a priori information required to obtain reliable solutions of inverse geophysical problems. In view of the classification developed, in this paper we study the type of nonuniqueness associated with the inverse refraction problem. An approach for obtaining a realistic solution to the inverse refraction problem is offered in a third paper that is in preparation. The nonuniqueness of the inverse refraction problem is examined by using a simple three-layer model. Like many other inverse geophysical problems, the inverse refraction problem does not have a unique solution. Conventionally, nonuniqueness is considered to be a result of insufficient data and/or error in the data, for any fixed number of model parameters. This study illustrates that even for overdetermined and error free data, nonlinear inverse refraction problems exhibit exact-data nonuniqueness, which further complicates the problem of nonuniqueness. By evaluating the nonuniqueness of the inverse refraction problem, this paper targets the improvement of refraction inversion algorithms, and as a result, the achievement of more realistic solutions. The nonuniqueness of the inverse refraction problem is examined initially by using a simple three-layer model. The observations and conclusions of the three-layer model nonuniqueness study are used to evaluate the nonuniqueness of more complicated n-layer models and multi-parameter cell models such as in refraction tomography. For any fixed number of model parameters, the inverse refraction problem exhibits continuous ranges of exact-data nonuniqueness. Such an unfavorable type of nonuniqueness can be uniquely solved only by providing abundant a priori information

  20. Taming the non-linearity problem in GPR full-waveform inversion for high contrast media

    NASA Astrophysics Data System (ADS)

    Meles, Giovanni; Greenhalgh, Stewart; van der Kruk, Jan; Green, Alan; Maurer, Hansruedi

    2012-03-01

    We present a new algorithm for the inversion of full-waveform ground-penetrating radar (GPR) data. It is designed to tame the non-linearity issue that afflicts inverse scattering problems, especially in high contrast media. We first investigate the limitations of current full-waveform time-domain inversion schemes for GPR data and then introduce a much-improved approach based on a combined frequency-time-domain analysis. We show by means of several synthetic tests and theoretical considerations that local minima trapping (common in full bandwidth time-domain inversion) can be avoided by starting the inversion with only the low frequency content of the data. Resolution associated with the high frequencies can then be achieved by progressively expanding to wider bandwidths as the iterations proceed. Although based on a frequency analysis of the data, the new method is entirely implemented by means of a time-domain forward solver, thus combining the benefits of both frequency-domain (low frequency inversion conveys stability and avoids convergence to a local minimum; whereas high frequency inversion conveys resolution) and time-domain methods (simplicity of interpretation and recognition of events; ready availability of FDTD simulation tools).

  1. Taming the non-linearity problem in GPR full-waveform inversion for high contrast media

    NASA Astrophysics Data System (ADS)

    Meles, Giovanni; Greenhalgh, Stewart; van der Kruk, Jan; Green, Alan; Maurer, Hansruedi

    2011-02-01

    We present a new algorithm for the inversion of full-waveform ground-penetrating radar (GPR) data. It is designed to tame the non-linearity issue that afflicts inverse scattering problems, especially in high contrast media. We first investigate the limitations of current full-waveform time-domain inversion schemes for GPR data and then introduce a much-improved approach based on a combined frequency-time-domain analysis. We show by means of several synthetic tests and theoretical considerations that local minima trapping (common in full bandwidth time-domain inversion) can be avoided by starting the inversion with only the low frequency content of the data. Resolution associated with the high frequencies can then be achieved by progressively expanding to wider bandwidths as the iterations proceed. Although based on a frequency analysis of the data, the new method is entirely implemented by means of a time-domain forward solver, thus combining the benefits of both frequency-domain (low frequency inversion conveys stability and avoids convergence to a local minimum; whereas high frequency inversion conveys resolution) and time-domain methods (simplicity of interpretation and recognition of events; ready availability of FDTD simulation tools).

  2. Joint analysis of refractions with surface waves: An inverse solution to the refraction-traveltime problem

    USGS Publications Warehouse

    Ivanov, J.; Miller, R.D.; Xia, J.; Steeples, D.; Park, C.B.

    2006-01-01

    We describe a possible solution to the inverse refraction-traveltime problem (IRTP) that reduces the range of possible solutions (nonuniqueness). This approach uses a reference model, derived from surface-wave shear-wave velocity estimates, as a constraint. The application of the joint analysis of refractions with surface waves (JARS) method provided a more realistic solution than the conventional refraction/tomography methods, which did not benefit from a reference model derived from real data. This confirmed our conclusion that the proposed method is an advancement in the IRTP analysis. The unique basic principles of the JARS method might be applicable to other inverse geophysical problems. ?? 2006 Society of Exploration Geophysicists.

  3. Decomposing Large Inverse Problems with an Augmented Lagrangian Approach: Application to Joint Inversion of Body-Wave Travel Times and Surface-Wave Dispersion Measurements

    NASA Astrophysics Data System (ADS)

    Reiter, D. T.; Rodi, W. L.

    2015-12-01

    Constructing 3D Earth models through the joint inversion of large geophysical data sets presents numerous theoretical and practical challenges, especially when diverse types of data and model parameters are involved. Among the challenges are the computational complexity associated with large data and model vectors and the need to unify differing model parameterizations, forward modeling methods and regularization schemes within a common inversion framework. The challenges can be addressed in part by decomposing the inverse problem into smaller, simpler inverse problems that can be solved separately, providing one knows how to merge the separate inversion results into an optimal solution of the full problem. We have formulated an approach to the decomposition of large inverse problems based on the augmented Lagrangian technique from optimization theory. As commonly done, we define a solution to the full inverse problem as the Earth model minimizing an objective function motivated, for example, by a Bayesian inference formulation. Our decomposition approach recasts the minimization problem equivalently as the minimization of component objective functions, corresponding to specified data subsets, subject to the constraints that the minimizing models be equal. A standard optimization algorithm solves the resulting constrained minimization problems by alternating between the separate solution of the component problems and the updating of Lagrange multipliers that serve to steer the individual solution models toward a common model solving the full problem. We are applying our inversion method to the reconstruction of the·crust and upper-mantle seismic velocity structure across Eurasia.· Data for the inversion comprise a large set of P and S body-wave travel times·and fundamental and first-higher mode Rayleigh-wave group velocities.

  4. FOREWORD: Tackling inverse problems in a Banach space environment: from theory to applications Tackling inverse problems in a Banach space environment: from theory to applications

    NASA Astrophysics Data System (ADS)

    Schuster, Thomas; Hofmann, Bernd; Kaltenbacher, Barbara

    2012-10-01

    Inverse problems can usually be modelled as operator equations in infinite-dimensional spaces with a forward operator acting between Hilbert or Banach spaces—a formulation which quite often also serves as the basis for defining and analyzing solution methods. The additional amount of structure and geometric interpretability provided by the concept of an inner product has rendered these methods amenable to a convergence analysis, a fact which has led to a rigorous and comprehensive study of regularization methods in Hilbert spaces over the last three decades. However, for numerous problems such as x-ray diffractometry, certain inverse scattering problems and a number of parameter identification problems in PDEs, the reasons for using a Hilbert space setting seem to be based on conventions rather than an appropriate and realistic model choice, so often a Banach space setting would be closer to reality. Furthermore, non-Hilbertian regularization and data fidelity terms incorporating a priori information on solution and noise, such as general Lp-norms, TV-type norms, or the Kullback-Leibler divergence, have recently become very popular. These facts have motivated intensive investigations on regularization methods in Banach spaces, a topic which has emerged as a highly active research field within the area of inverse problems. Meanwhile some of the most well-known regularization approaches, such as Tikhonov-type methods requiring the solution of extremal problems, and iterative ones like the Landweber method, the Gauss-Newton method, as well as the approximate inverse method, have been investigated for linear and nonlinear operator equations in Banach spaces. Convergence with rates has been proven and conditions on the solution smoothness and on the structure of nonlinearity have been formulated. Still, beyond the existing results a large number of challenging open questions have arisen, due to the more involved handling of general Banach spaces and the larger variety

  5. ARISE antenna

    NASA Astrophysics Data System (ADS)

    Chmielewski, Arthur B.; Noca, Muriel; Ulvestad, James

    2000-03-01

    Supermassive black holes are among the most spectacular objects in the Universe, and are laboratories for physics in extreme conditions. Understanding the physics of massive black holes and related phenomena is a primary goal of the ARISE mission. The scientific goals of the mission are described in detail on the ARISE web site http://arise.ipl.nasa.gov and in the ARISE Science Goals document. The following paper, as the title suggests, is not intended to be a comprehensive description of ARISE, but deals only with one aspect of the ARISE mission-the inflatable antenna which is the key element of the ARISE spacecraft. This spacecraft,due to the extensive reliance on inflatables, may be considered as the first generation Gossamer spacecraft

  6. A numerical study of the inverse problem of breast infrared thermography modeling

    NASA Astrophysics Data System (ADS)

    Jiang, Li; Zhan, Wang; Loew, Murray H.

    2010-03-01

    Infrared thermography has been shown to be a useful adjunctive tool for breast cancer detection. Previous thermography modeling techniques generally dealt with the "forward problem", i.e., to estimate the breast thermogram from known properties of breast tissues. The present study aims to deal with the so-called "inverse problem", namely to estimate the thermal properties of the breast tissues from the observed surface temperature distribution. By comparison, the inverse problem is a more direct way of interpreting a breast thermogram for specific physiological and/or pathological information. In tumor detection, for example, it is particularly important to estimate the tumor-induced thermal contrast, even though the corresponding non-tumor thermal background usually is unknown due to the difficulty of measuring the individual thermal properties. Inverse problem solving is technically challenging due to its ill-posed nature, which is evident primarily by its sensitivity to imaging noise. Taking advantage of our previously developed forward-problemsolving techniques with comprehensive thermal-elastic modeling, we examine here the feasibility of solving the inverse problem of the breast thermography. The approach is based on a presumed spatial constraint applied to three major thermal properties, i.e., thermal conductivity, blood perfusion, and metabolic heat generation, for each breast tissue type. Our results indicate that the proposed inverse-problem-solving scheme can be numerically stable under imaging noise of SNR ranging 32 ~ 40 dB, and that the proposed techniques can be effectively used to improve the estimation to the tumor-induced thermal contrast, especially for smaller and deeper tumors.

  7. Numerical solution of 2D-vector tomography problem using the method of approximate inverse

    NASA Astrophysics Data System (ADS)

    Svetov, Ivan; Maltseva, Svetlana; Polyakova, Anna

    2016-08-01

    We propose a numerical solution of reconstruction problem of a two-dimensional vector field in a unit disk from the known values of the longitudinal and transverse ray transforms. The algorithm is based on the method of approximate inverse. Numerical simulations confirm that the proposed method yields good results of reconstruction of vector fields.

  8. Integro-differential method of solving the inverse coefficient heat conduction problem

    NASA Astrophysics Data System (ADS)

    Baranov, V. L.; Zasyad'Ko, A. A.; Frolov, G. A.

    2010-03-01

    On the basis of differential transformations, a stable integro-differential method of solving the inverse heat conduction problem is suggested. The method has been tested on the example of determining the thermal diffusivity on quasi-stationary fusion and heating of a quartz glazed ceramics specimen.

  9. The application of subspace preconditioned LSQR algorithm for solving the electrocardiography inverse problem.

    PubMed

    Jiang, Mingfeng; Xia, Ling; Huang, Wenqing; Shou, Guofa; Liu, Feng; Crozier, Stuart

    2009-10-01

    Regularization is an effective method for the solution of ill-posed ECG inverse problems, such as computing epicardial potentials from body surface potentials. The aim of this work was to explore more robust regularization-based solutions through the application of subspace preconditioned LSQR (SP-LSQR) to the study of model-based ECG inverse problems. Here, we presented three different subspace splitting methods, i.e., SVD, wavelet transform and cosine transform schemes, to the design of the preconditioners for ill-posed problems, and to evaluate the performance of algorithms using a realistic heart-torso model simulation protocol. The results demonstrated that when compared with the LSQR, LSQR-Tik and Tik-LSQR method, the SP-LSQR produced higher efficiency and reconstructed more accurate epcicardial potential distributions. Amongst the three applied subspace splitting schemes, the SVD-based preconditioner yielded the best convergence rate and outperformed the other two in seeking the inverse solutions. Moreover, when optimized by the genetic algorithms (GA), the performances of SP-LSQR method were enhanced. The results from this investigation suggested that the SP-LSQR was a useful regularization technique for cardiac inverse problems. PMID:19564127

  10. ON THE GEOSTATISTICAL APPROACH TO THE INVERSE PROBLEM. (R825689C037)

    EPA Science Inventory

    Abstract

    The geostatistical approach to the inverse problem is discussed with emphasis on the importance of structural analysis. Although the geostatistical approach is occasionally misconstrued as mere cokriging, in fact it consists of two steps: estimation of statist...