INVERSE PROBLEMS NEWSLETTER
DOI: https://doi.org/10.1088/0266-5611/14/3/025
IF: 2.408
1998-06-01
Inverse Problems
Abstract:The Newsletter is a key element in further enhancing the value of the journal to the inverse problems community. So why not be a part of this exciting forum by sending to our Bristol office material suitable for inclusion under any of the categories mentioned above. Your contributions will be very welcome. Forthcoming event Conference announcement 2nd International Conference on Identification in Engineering Systems, University of Wales Swansea, 29 - 31 March 1999 Parameter estimation and system identification are used extensively to obtain dynamic models of engineering systems. A large number of methods have been developed and a huge amount of experience gained in their application, particularly in control engineering and structural dynamics. Although aspects of the methods are different there is a substantial overlap in the methodology and practice used in the identification. Following the success of the first conference held in March 1996 in Swansea, this conference will provide a forum for researchers and practitioners in the art and science of identification from a range of disciplines and provide a further impetus to the cross-fertilization of ideas in this area. Further details, including the procedure for the submission of abstracts, are available at URL: http://www.swan.ac.uk/mecheng/ies99/ or contact Dr M I Friswell (e-mail: m.i.friswell@swansea.ac.uk, fax: +44 (0)1792 295676). Book reviews 75 Years of Radon Transform S Gindikin and P Michor (ed) 1994 Boston, MA: International 303 pp ISBN 1-57146-008-X $42.00 In 1917 Johann Radon published his now celebrated formula for reconstructing a function (in many applications, this means the density of an object) from line integrals through it in different directions [1]. This paper now stands out as one of the most fundamental contributions of all time to the subject of inverse problems. It inspired much fundamental work in mathematics, such as a notable series of papers by Fritz John in the 1930s. However, the inventors of computerized tomography (CT) for medical imaging could not trace any of this when they needed exactly this inversion in the 1960s. The full significance of Radon's pioneering contribution did not become universally recognized until the 1970s. In 1917, early on in his career, Johann Radon was an assistant and a Privatdozent at the University of Vienna. He would later spend 25 years at various German universities before finishing his career as Dean, and later Rector, at his alma mater, the University of Vienna. In 1992, the 75th anniversary of the transform was celebrated at the University of Vienna with a conference. This book forms the proceedings of that conference. The book has three major parts: Biographical contributions. This part contains, among other items, a biography written by Radon's daughter Brigitte Burkovics, reminiscences by one of his students and by Fritz John, and a description by Allen Cormack about how he approached the line-integral inversion problem when developing CT. Scientific contributions. This part constitutes the bulk of the book; it contains selected presentations from the technical sessions of the conference - not easy reading for the mathematically faint-hearted. Reprinted papers. We find here some notes by the Editor of the volume, followed by a now classical (1938) contribution by Fritz John, and the original 1917 paper by Radon. One thing the book does not contain (in line with the mathematical theme of the conference) is any summary of how the transform - in its many forms of implementation - is now used throughout a vast array of applications. For a general scientific audience, the first part contains fascinating historical reading, and the third part provides for pleasant casual browsing (as well as for in-depth study). The second part is strictly for specialists in the mathematical aspects of the Radon transform. Reference[1]Radon J 1917 Über die Bestimmung von Functionen durch ihre integralwerte längs gewisser Mannigfaltigkeiten Ber. Sächsische Akad. der Wissenschaften, Leipzig, Math.-Nat. Kl. 69 262 - 77 B Fornberg University of Colorado, Boulder Inverse Problems of Wave Propagation and Diffraction Proceedings, Aix-les-Bains, France 1996, Lecture Notes in Physics, Volume 486 G Chavent and P C Sabatier (ed) 1997 Berlin: Springer 379 pp ISBN 3-540-62865-7 DM106.00, öS774.00, sFr96.50, £ 41.00, $79.00 This book is a collection of lectures that were presented at the `Conference on Inverse Problems of Wave Propagation and Diffraction'. This meeting was organized by two renowned French experts in inverse problems, Guy Chavent and Pierre C Sabatier, in Aix-les-Bains in September 1996. What impressed me about this book was the clear orientation of the vast majority of the papers towards applications. The range of these applications is quite broad. In terms of fields of study, it varies from `electromagnetic' to `acoustic' inverse problems. Applications naturally dictate a broad variety of mathematical models as well as the appropriate mathematical tools. Regarding the models, the majority of authors consider the frequency domain regime, thus dealing with the Helmholz-like equation. Some authors, however, also work with the inverse problems for the hyperbolic (i.e. time dependent) equations. The goal of almost all the papers is clear and can be appreciated: to develop and test effective numerical methods for the inverse problems under consideration. The problems studied here can mainly be split into two categories, which have rather blurry boundaries: inverse obstacle problems and inverse coefficient problems. In the first class, one is recovering the shape and location of an obstacle given scattered data. Often, this is an impenetrable obstacle. In the second class, the goal is to find an internal structure of a penetrable object, not only its shape. As a result, the second class of problems amounts to the determination of an unknown coefficient(s) of a PDE. An interesting observation is that the majority of authors discussing the inverse obstacle problems are profoundly influenced by, I would say, a classic book by D Colton and R Kress: Inverse Acoustic and Electromagnetic Scattering Theory. A point which is made, implicitly or explicitly, in almost all the papers is the importance of reasonable a priori constraints on the solution of an inverse problem. While a wide variety of numerical approaches is discussed (given the broad scope of inverse problems under consideration), this is perhaps the only point which is common to all the approaches. This falls well into the fundamental Tikhonov's principle of an a priori choice of an appropriate compact set of solutions. Here are some examples of these a priori constraints. The paper by D Colton describes a fresh idea to recover only support of a penetrable anomaly rather than the value of the unknown coefficient within it. In this way the original nonlinear problem is reduced to a linear one, which can be solved in an elegant way using the fact that the fundamental solution of the Helmholz-like PDE has a singularity at the source location. The paper by F Natterer represents another example of a novel and effective numerical approach which was developed using the concept of a priori constraints on solutions. An intruiging core idea of this paper is to solve a classical ill-posed Cauchy problem for the elliptic equation in a well-posed fashion, assuming a priori that the grid size should not be too small. This, in turn, leads to a rapid image reconstruction algorithm, in 3D, though only in the high frequency regime. While the works of Colton, Natterer and some other authors consider the cases of full view illumination and measurements of the target medium, which is quite acceptable in medical imaging, other authors consider the limited view case, which adds up to the ill-posedness of the problems under consideration. This is certainly true for a number of works devoted to inverse seismic problems, which are notoriously challenging. The difficulty of these problems is clearly demonstrated in the paper by L Fatone, P Maponi, C Pignotti and F Zirilli, in which the 1D inverse problem for the 2D (in space) hyperbolic equation is studied. Probably the major challenge, even in the 1D case, consists of figuring out discontinuities of the speed ) in the layered medium, given measurements of the backscattered data on the surface. Although the paper by A Litman, D Lesselier and F Santosa deals with an electromagnetic, rather than with a seismic, inverse problem, it still works with the incomplete data collection. The goal of these authors is to reconstruct the shape and location of a defect given electrical parameters of this object and incomplete view measurements of the scattered electrical field. This problem is solved by introducing a level set function which describes the shapes of a family of defects (including the correct one) and which satisfies a certain Hamilton - Jacobi-type PDE. This PDE, in turn, is connected with a least-squares cost functional. The tutorial paper by M Bertero, P Boccacci and M Piana stays somewhat outside of the common scope of the other papers. Nevertheless, this is a very interesting work. It clarifies a rather ambiguous issue of the achievable resolution R versus the wavelength . The common view is that . The authors show, however, that this is true only in the far-field data. Contrary to this, in the case of the near-field data, resolution can be very much less than the wavelength, because of the evanescent waves, whose intensity decays exponentially with the distance from the source. It is impossible for me to comment on all the papers included in this book. My overall comment, however, is that each of the works presented in this book is interesting in its own right, and I enjoyed reading all of them, thanks to the contributors and the editors. In my opinion, this book represents a state of the art (as of 1996) collection of numerical and some theoretical approaches to inverse problems of wave propagation. M V Klibanov University of North Carolina, Charlotte Inverse Nodal Problems: Finding the Potential from Nodal Lines Memoirs of the American Mathematical Society, Number 572 O H Hald and J R McLaughlin 1996 Providence, RI: American Mathematical Society 148 pp ISBN 0-8218-0486-3 £ 30.00 This book describes new and sometimes unexpected properties of an eigenvalue problem for an elliptic equation, in potential form, on a rectangle with Dirichlet boundary conditions. First the Sturm - Liouville problem for a Laplacian is considered, and eigenvalues (and thus eigenfunctions) can be easily found. The authors studied the behaviour of eigenvalues and proved that for most rectangles almost all eigenvalues are separated from other eigenvalues by a specified gap; they also found other useful features. Next the authors describe properties of the same problem with a sufficiently smooth potential function, using perturbation theory. The main result is the proof that the potential is uniquely determined, up to an additive constant, by a subset of the nodal lines of the eigenfunctions. The authors present a formula that yields an approximation to the potential at a dense set of points. The book is interesting to read; it could be recommended to all specialists in mathematical physics and inverse problems. A Yagola Moscow State University Formulas in Inverse and Ill-Posed Problems Yu E Anikonov 1997 Utrecht: VSP 240 pp ISBN 90-6764-216-9 £ 82.00 The rather unusual title of this book is actually quite apt; it is a collection of representation formulae for a wide variety of inverse problems, but primarily those of evolution type. Lest the reader assume that the holy grail has at last been found, I should point out that not all of the formulae give explicit representations; indeed for most of the inverse problems that are inherently nonlinear the representation is an implicit one. Also, many of the problems are not one of the standards one finds in the current literature or motivated by practical considerations (at least as far as the reviewer is aware). The book's first example is illustative: to recover the coefficient in the parabolic equation , given an initial condition and boundary values on . The overposed data considered are the values of the Fourier transform , for each . While this allows the coefficient to be (formally) recovered, it is not the sort of data that one expects to provide in diffusion problems. Despite this, the book is a mine of information and useful ideas. It would certainly repay anyone interested in a particular inverse problem to examine it. It is not a book one reads but one that is consulted with the understanding that, while it will very likely not contain ready-made answers, there is reasonable hope that it might provide some insight. Proofs of the claims and justification of the formulae are rarely complete, but there is an extensive list of references, although these are certainly not representative of the areas claimed to be covered in the text. To amplify the last remark, more than two thirds of the citations are to the Russian literature; half of these are by the book's author. Of the remaining third, half are to books or articles published before 1975 and only 20% to works of the last ten years. Thus the recent explosion in the inverse problems literature on a worldwide basis over the last decade or so is largely unaddressed. Who will or should buy this book? As the above remarks indicate, it is not really suitable as a textbook for graduate courses, at least not for programmes in western Europe or the USA. I suspect that the market here is the research libraries, as I don't imagine many individuals will find it to be a must-have item for their bookshelves, and the very high cost per page will not help in this regard. W Rundell Texas A&M University, College Station Inverse Problems in Wave Propagation The IMA Volumes in Mathematics and its Applications, Volume 90 G Chavent, G Papanicolaou, P Sacks and W Symes (ed) 1997 Berlin: Springer 498 pp ISBN 0-387-94976-3 DM118.00, öS861.40, sFr104.00, £ 48.50, $96.95 This book is based on the proceedings of a two-week workshop which was an integral part of the 1994 - 1995 IMA programme `Waves and Scattering'. This workshop took place from 6 - 17 March 1995. There are 24 papers. This book represents most applications in inverse wave propagation problems, together with fundamental mathematical investigations of the relation between waves and scatterers. The following contributions are included. R A Albanese discusses wave propagation issues in medicine and environmental health. The uniqueness of one-dimensional reconstruction for orthogonal and oblique incidence is touched upon. J G Berryman shows that the inverse problem with the data serving as constraints is most easily analysed when it is possible to segment the solution space into regions that are feasible (satisfying all the constraints) and infeasible (violating some of the constraints). The variational structure of three inverse problems has been investigated. R W Brookes and K P Bube investigate the numerical convergence of a finite-difference method in space and time for one-dimensional models in reflection seismology, and discuss the order of convergence in relation to the source wavelet. M D Collins discusses topics in ocean acoustic inverse problems, including remote sensing and localization problems. Efficient techniques are discussed for the solution of the forward problem, including poro-elastic media, and techniques for solving global scale acoustic problems have been adapted. D Colton discusses in general the three-dimensional electromagnetic inverse scattering problem for inhomogeneous objects and in particular an example in medical imaging. A two-dimensional problem is solved numerically with the dual space method. E Croc and Y Dermenjian study the generalized modes in an acoustic strip, which is a simplified model of a seismic experiment, where both source and receiver are situated in a well. G Eskin and J Ralston consider the inverse problems for Schrödinger operators with magnetic and electric potentials. A Faridani shows old and new results in computed tomography, which entails the reconstruction of a function f from line integrals of f. Problems of uniqueness, reconstruction formulae for the x-ray transform, filtered backprojection, error estimates and incomplete data problems are treated in detail. D J Foster, R G Keys and D P Schmitt propose a theoretical basis for interpreting amplitude versus offset (AVO) inversion. It defines the background seismic response and characterizes anomalous events by distance from this background. F A Grunbaum and S K Patch investigate how many parameters one can solve for in diffuse tomography. Some bounds of the range of the parameters are given. J G Harris constructs some physically reasonable models of acoustic imaging. In particular, to investigate small surface-breaking cracks, the leaky Rayleigh wave is used in the imaging mechanism. V Isakov analyses the reconstruction of the diffusion and of the principal coefficient of a hyperbolic equation, in particular the diffusion coefficient, by the use of beam solutions and the recovery of discontinuity of the wave speed. V G Khajdukov, V I Kostin and V A Tcheverda consider the r-solution and its applications in a linearized waveform inversion for a layered background. Y Kurylev and A Starkov deal with an inverse boundary problem for a second-order elliptic operator and its nonstationary counterparts. The uniqueness of the reconstruction of a density parameter is proved and a direct procedure of its reconstruction using directional moments is described. Ching-Ju Ashraf Lee and J R McLaughlin solve an inverse nodal problem for a rectangular membrane using the ratio method and the method of parameter identification. Changmei Liu derives a uniqueness theorem for any ball and for two balls, when the scattering amplitudes for some independent incident directions are given. J R McLaughlin, P E Sacks and M Somasundaram discuss the inverse problem of the determination of acoustic parameters using far-field data, by using the properties of the interior eigenvalues, and show how the sound speed depending on the vertical coordinate can be determined. G Hanamura and G Uhlmann describe a layer stripping algorithm in elastic impedance tomography, in which an elastic tensor has to be reconstructed from measurements of the displacement and the traction at the boundary of the domain. G Nolet uses the WKBJ or path integral approximation to the solution of the elastodynamic equations in a slightly heterogeneous earth to partition the inverse problem for a large set of observed seismograms along different wavepaths. R L Nowack discusses an inverse method for the analysis of refraction and wide-angle seismic data and illustrates it through the inversion of a shallow crustal structure. F R Pijpers presents a brief overview of applications of inversions within astronomy, including a recent modification of the Backus and Gilbert method. E L Ritman, J H Dunsmuir, A Faridani, D V Finch, K T Smith and P J Thomas present an example in which local reconstruction extends the capability of a micro-CT scanner beyond the physical limits imposed by global tomographic reconstruction techniques. J Sylvester discusses the layer stripping problem. Some theorems are cited as evidence that it provides a productive theoretical method which yields new insights into an old problem. M E Taylor studies estimates for approximate solutions to acoustic inverse scattering problems, in particular the recovery of the near-field wave from the scattering amplitude and the consequences of linearization of the inverse problems. This book is certainly an important contribution to the theory and application of inverse wave propagation problems and is a good reference work that should be present in the library. P M van den Berg Delft University of Technology
mathematics, applied,physics, mathematical
What problem does this paper attempt to address?
-
Recently Developed Formulations of the Inverse Problem in Acoustics and Electromagnetics
Norman Bleistein,Norbert N. Bojarski
DOI: https://doi.org/10.21236/ada003588
1974-12-01
Abstract:Recently developed formulations of the inverse problems in acoustics and electromagnetics are described. There are two types of formulations, one in the geometrical optics limit and the other, an exact formulation for the inverse source problem. Both basic formulations are extended to include the realistic problem of a 'limited aperture' of observations. It is also shown that the inverse source formulation can be applied to the problem of reconstruction of media inhomogeneities from remotely sensed field data. The basic physical optics result is that the characteristic function of the scattering obstacle and the phase and range normalized scattering amplitude are a Fourier transform pair. All other formulations lead to Fredholm integral equations of the first kind.
-
Some Reconstruction Methods for Inverse Scattering Problems
jijun liu,haibing wang
DOI: https://doi.org/10.1007/978-3-642-13742-6_10
2010-01-01
Abstract:Inverse scattering problems are one of the main research areas in the optimization techniques. The main purpose of inverse scattering problems is to detect the physical properties of an obstacle from some information related to the scattered waves of the obstacle for given incident wave. Generally, if the incident plane waves are given from the finite number of directions, which are indeed the practical situations, there is no uniqueness for reconstructing the obstacle properties such as the boundary shape. In these cases, the optimization techniques can be applied to reconstructing the obstacle approximately. That is, the obstacle shape is approximated by a minimizer of some cost functional which measures the defect between the measurement data of the scattered wave and the computational scattered wave related to the approximate obstacle. Of course, for this optimization problems in infinite dimensional space, some regularizing term should be introduced to the cost functional. Although these general optimization techniques have been applied widely in the last century, they also suffer from many disadvantages theoretically and numerically. From the theoretical point of view, the lack of uniqueness makes the obtained approximate obstacle from the optimization procedure ambiguous, namely, we do not know whether or not the computational result indeed approximates the true one. From the numerical respects, such an optimization procedure needs to solve the direct scattering problem at each iteration step, which entails huge time cost. Moreover, to get the convergence of the iteration, a good initial guess is required. Even if in the case of convergence, the approximate sequence generally approaches to some local minimizer These shortcomings of the direct optimizations used in inverse scattering problems cause some problems for the numerical reconstruction of the obstacles. In recent years, some modified version of optimization algorithms based on the potential methods for inverse scattering problems are proposed, in the sense that the direct problems are not required to solve in the iteration procedures. On the other hand, some new schemes combining the advantages of optimizations with the exact reconstruction formulas have been developed for inverse scattering problems. In this chapter, we give an overview on these two directions. We begin with a brief introduction to the physical background to acoustic scattering problems as well as some well-known inverse scattering models. Then we review some classical and recently developed inversion methods for detecting the information about the unknown scatterer from a knowledge of the far-field pattern u(infinity) for one or several incident plane waves. For each method, the basic idea is described, and some main results are presented. To test their validity, the numerical implementation of all inversion methods needs to be studied. So we finally focus on the numerical realizations of these existing methods, pointing out the main difficulties encountered in numerical realization. In addition, the advantages and disadvantages of these methods are also analyzed.
-
Preface to the Special Issue on Inverse Problems.
Jin Cheng,Benny Y. C. Hon,Masahiro Yamamoto
DOI: https://doi.org/10.1007/s10444-010-9166-z
2010-01-01
Advances in Computational Mathematics
Abstract:The investigation of inverse problems is important due to its significance in practice to the real application fields such as medical diagnosis, materials sciences, geophysics, and industrial engineering.The recent rapid development in numerical methods for solving inverse problems has been attracting much attention from mathematicians, engineers, and scientific researchers.Inverse problems are in nature concerned with the detection of hidden "reasons" from incomplete observations of "outcomes".In practice, most of these detections, testings, and diagnoses possess the character of inverse problems: intrinsic instability.In view of the significance in practical applications, exploitations of feasible and robust numerical methods that can handle well the instability of inverse problem are indispensable.Some of these numerical methods have been exploited due to the pioneering "regularization" work by Tikhonov.Witnessing the recent rapid development in computational technologies and measurement equipments, we are now at the edge towards the next stage when more variety of numerical methods can be developed for real and applied inverse problems.Moreover, with the substantial accumulation of theoretical results for inverse problems, novel numerical methods can be invented.
-
Inverse problems for PDEs: Models, computations and applications
Cheng Jin,Liu Jijun,Zhang Bo
DOI: https://doi.org/10.1360/N012018-00076
2019-01-01
Abstract:Inverse problems for partial differential equations (PDEs) are of great importance in the areas of applied mathematics, whichcover different mathematical branches including PDEs, functional analysis, nonlinear analysis,optimizations, regularization and numerical analysis. These problems have found wide applicationsin many important engineering areas such as media imaging, remote sensing and image processing.Due to the nature of ill-posedness of such kinds of problems, the techniques of regularizationshould be applied for efficiently solving these problems. However, it is very hard to establisha unified framework for inverse problems of PDEs, due to the variety and complexity of the problems.This paper aims to give an overview on several important inverse problems of PDEs models.Based on the systematic recalls on the origins and specialities of inverse problems for PDEs,we focus on three kinds of PDEs models for inverse problems: electrical impedance tomograph, inverse wave scattering, thermal imaging. The crucial problems of fundamental interests, existing resultsand methods as well as further possible research directions are reviewed. Moreover, we also give asystematic analysis of numerical methods for solving inverse problems for PDEs.
-
Bayesian inverse problems for functions and applications to fluid mechanics
S L Cotter,M Dashti,J C Robinson,A M Stuart
DOI: https://doi.org/10.1088/0266-5611/25/11/115008
IF: 2.408
2009-10-29
Inverse Problems
Abstract:In this paper we establish a mathematical framework for a range of inverse problems for functions, given a finite set of noisy observations. The problems are hence underdetermined and are often ill-posed. We study these problems from the viewpoint of Bayesian statistics, with the resulting posterior probability measure being defined on a space of functions. We develop an abstract framework for such problems which facilitates application of an infinite-dimensional version of Bayes theorem, leads to a well-posedness result for the posterior measure (continuity in a suitable probability metric with respect to changes in data), and also leads to a theory for the existence of maximizing the posterior probability (MAP) estimators for such Bayesian inverse problems on function space. A central idea underlying these results is that continuity properties and bounds on the forward model guide the choice of the prior measure for the inverse problem, leading to the desired results on well-posedness and MAP estimators; the PDE analysis and probability theory required are thus clearly dileneated, allowing a straightforward derivation of results. We show that the abstract theory applies to some concrete applications of interest by studying problems arising from data assimilation in fluid mechanics. The objective is to make inference about the underlying velocity field, on the basis of either Eulerian or Lagrangian observations. We study problems without model error, in which case the inference is on the initial condition, and problems with model error in which case the inference is on the initial condition and on the driving noise process or, equivalently, on the entire time-dependent velocity field. In order to undertake a relatively uncluttered mathematical analysis we consider the two-dimensional Navier–Stokes equation on a torus. The case of Eulerian observations—direct observations of the velocity field itself—is then a model for weather forecasting. The case of Lagrangian observations—observations of passive tracers advected by the flow—is then a model for data arising in oceanography. The methodology which we describe herein may be applied to many other inverse problems in which it is of interest to find, given observations, an infinite-dimensional object, such as the initial condition for a PDE. A similar approach might be adopted, for example, to determine an appropriate mathematical setting for the inverse problem of determining an unknown tensor arising in a constitutive law for a PDE, given observations of the solution. The paper is structured so that the abstract theory can be read independently of the particular problems in fluid mechanics which are subsequently studied by application of the theory.
mathematics, applied,physics, mathematical
-
A mathematical framework for inverse wave problems in heterogeneous media
Kirk D. Blazek,Christiaan C. Stolk,William W. Symes
DOI: https://doi.org/10.1088/0266-5611/29/6/065001
2012-11-20
Abstract:This paper provides a theoretical foundation for some common formulations of inverse problems in wave propagation, based on hyperbolic systems of linear integro-differential equations with bounded and measurable coefficients. The coefficients of these time-dependent partial differential equations respresent parametrically the spatially varying mechanical properties of materials. Rocks, manufactured materials, and other wave propagation environments often exhibit spatial heterogeneity in mechanical properties at a wide variety of scales, and coefficient functions representing these properties must mimic this heterogeneity. We show how to choose domains (classes of nonsmooth coefficient functions) and data definitions (traces of weak solutions) so that optimization formulations of inverse wave problems satisfy some of the prerequisites for application of Newton's method and its relatives. These results follow from the properties of a class of abstract first-order evolution systems, of which various physical wave systems appear as concrete instances. Finite speed of propagation for linear waves with bounded, measurable mechanical parameter fields is one of the by-products of this theory.
Mathematical Physics,Analysis of PDEs
-
An introduction to finite element methods for inverse coefficient problems in elliptic PDEs
Bastian Harrach
DOI: https://doi.org/10.1365/s13291-021-00236-2
2021-05-27
Abstract:Several novel imaging and non-destructive testing technologies are based on reconstructing the spatially dependent coefficient in an elliptic partial differential equation from measurements of its solution(s). In practical applications, the unknown coefficient is often assumed to be piecewise constant on a given pixel partition (corresponding to the desired resolution), and only finitely many measurements can be made. This leads to the problem of inverting a finite-dimensional non-linear forward operator $\mathcal F:\ \mathcal D(\mathcal F)\subseteq \mathbb R^n\to \mathbb R^m$, where evaluating $\mathcal F$ requires one or several PDE solutions.
Numerical inversion methods require the implementation of this forward operator and its Jacobian. We show how to efficiently implement both using a standard FEM package and prove convergence of the FEM approximations against their true-solution counterparts. We present simple example codes for Comsol with the Matlab Livelink package, and numerically demonstrate the challenges that arise from non-uniqueness, non-linearity and instability issues. We also discuss monotonicity and convexity properties of the forward operator that arise for symmetric measurement settings.
Numerical Analysis,Analysis of PDEs
-
Hybrid inverse problems and internal functionals
Guillaume Bal
DOI: https://doi.org/10.48550/arXiv.1110.4733
2011-10-21
Analysis of PDEs
Abstract:This paper reviews recent results on hybrid inverse problems, which are also called coupled-physics inverse problems of multi-wave inverse problems. Inverse problems tend to be most useful in, e.g., medical and geophysical imaging, when they combine high contrast with high resolution. In some settings, a single modality displays either high contrast or high resolution but not both. In favorable situations, physical effects couple one modality with high contrast with another modality with high resolution. The mathematical analysis of such couplings forms the class of hybrid inverse problems. Hybrid inverse problems typically involve two steps. In a first step, a well-posed problem involving the high-resolution low-contrast modality is solved from knowledge of boundary measurements. In a second step, a quantitative reconstruction of the parameters of interest is performed from knowledge of the point-wise, internal, functionals of the parameters reconstructed during the first step. This paper reviews mathematical techniques that have been developed in recent years to address the second step. Mathematically, many hybrid inverse problems find interpretations in terms of linear and nonlinear (systems of) equations. In the analysis of such equations, one often needs to verify that qualitative properties of solutions to elliptic linear equations are satisfied, for instance the absence of any critical points. This paper reviews several methods to prove that such qualitative properties hold, including the method based on the construction of complex geometric optics solutions.
-
Parametric Level Set Methods for Inverse Problems
Alireza Aghasi,Misha Kilmer,Eric L. Miller
DOI: https://doi.org/10.1137/100800208
2011-10-06
Abstract:In this paper, a parametric level set method for reconstruction of obstacles in general inverse problems is considered. General evolution equations for the reconstruction of unknown obstacles are derived in terms of the underlying level set parameters. We show that using the appropriate form of parameterizing the level set function results a significantly lower dimensional problem, which bypasses many difficulties with traditional level set methods, such as regularization, re-initialization and use of signed distance function. Moreover, we show that from a computational point of view, low order representation of the problem paves the path for easier use of Newton and quasi-Newton methods. Specifically for the purposes of this paper, we parameterize the level set function in terms of adaptive compactly supported radial basis functions, which used in the proposed manner provides flexibility in presenting a larger class of shapes with fewer terms. Also they provide a "narrow-banding" advantage which can further reduce the number of active unknowns at each step of the evolution. The performance of the proposed approach is examined in three examples of inverse problems, i.e., electrical resistance tomography, X-ray computed tomography and diffuse optical tomography.
Analysis of PDEs,Mathematical Physics,Geophysics
-
Revisit to the Inverse Exponential Radon Transform
Jason You
DOI: https://doi.org/10.48550/arXiv.2002.01622
2020-02-05
Abstract:This revisit gives a survey on the analytical methods for the inverse exponential Radon transform which has been investigated in the past three decades from both mathematical interests and medical applications such as nuclear medicine emission imaging. The derivation of the classical inversion formula is through the recent argument developed for the inverse attenuated Radon transform. That derivation allows the exponential parameter to be a complex constant, which is useful to other applications such as magnetic resonance imaging and tensor field imaging. The survey also includes the new technique of using the finite Hilbert transform to handle the exact reconstruction from 180 degree data. Special treatment has been paid on two practically important subjects. One is the exact reconstruction from partial measurements such as half-scan and truncated-scan data, and the other is the reconstruction from diverging-beam data. The noise propagation in the reconstruction is touched upon with more heuristic discussions than mathematical inference. The numerical realizations of several classical reconstruction algorithms are included. In the conclusion, several topics are discussed for more investigations in the future.
Image and Video Processing,Information Theory,Machine Learning,Numerical Analysis
-
Inverse Scattering Problems with Multi-Frequencies
Gang Bao,Peijun Li,Junshan Lin,Faouzi Triki
DOI: https://doi.org/10.1088/0266-5611/31/9/093001
IF: 2.408
2015-01-01
Inverse Problems
Abstract: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.
-
On reconstruction formulas and algorithms for the thermoacoustic tomography
M. Agranovsky,P. Kuchment,L. Kunyansky
DOI: https://doi.org/10.48550/arXiv.0706.1303
2007-10-22
Abstract:The paper surveys recent progress in establishing uniqueness and developing inversion formulas and algorithms for the thermoacoustic tomography. In mathematical terms, one deals with a rather special inverse problem for the wave equation. In the case of constant sound speed, it can also be interpreted as a problem concerning the spherical mean transform.
Analysis of PDEs,Numerical Analysis
-
Imaging of nonlinear materials via the Monotonicity Principle
Vincenzo Mottola,Antonio Corbo Esposito,Gianpaolo Piscitelli,Antonello Tamburrino
DOI: https://doi.org/10.1088/1361-6420/ad22e9
IF: 2.408
2024-01-27
Inverse Problems
Abstract:Inverse problems, which are related to Maxwell's equations, in the presence of nonlinear materials is a quite new topic in the literature. The lack of contributions in this area can be ascribed to the significant challenges that such problems pose. Retrieving the spatial behaviour of some unknown physical property, from boundary measurements, is a nonlinear and highly ill-posed problem even in the presence of linear materials. Furthermore, this complexity grows exponentially in the presence of nonlinear materials.
In the tomography of linear materials, the Monotonicity Principle (MP) is the foundation of a class of non-iterative algorithms able to guarantee excellent performances and compatibility with real-time applications. Recently, the MP has been extended to nonlinear materials under very general assumptions. Starting from the theoretical background for this extension, we develop a first real-time inversion method for the inverse obstacle problem in the presence of nonlinear materials.
The proposed method is intendend for all problems governed by the quasilinear Laplace equation, i.e. static problems involving nonlinear materials.
In this paper, we provide some preliminary results which give the foundation of our method and some extended numerical examples.
mathematics, applied,physics, mathematical
-
Less is often more: Applied inverse problems using -forward models.
Danny Smyl,Dong Liu
DOI: https://doi.org/10.1016/j.jcp.2019.108949
2019-01-01
Abstract:To solve an applied inverse problem, a numerical forward model for the problem's physics is required. Commonly, the finite element method is employed with discretizations consisting of elements with variable size h and polynomial degree p. Solutions to hp-forward models are known to converge exponentially by simultaneously decreasing h and increasing p. On the other hand, applied inverse problems are often ill-posed and their minimization rate exhibits uncertainty. Presently, the behavior of applied inverse problems incorporating hp elements of differing p, h, and geometry is not fully understood. Nonetheless, recent research suggests that employing increasingly higher-order hp-forward models (increasing mesh density and p) decreases reconstruction errors compared to inverse regimes using lower-order hp-forward models (coarser meshes and lower p). However, an affirmative or negative answer to following question has not been provided, “Does the use of higher order hp-forward models in applied inverse problems always result in lower error reconstructions than approaches using lower order hp-forward models?”In this article we aim to reduce the current knowledge gap and answer the open question by conducting extensive numerical investigations in the context of two contemporary applied inverse problems: elasticity imaging and hydraulic tomography – nonlinear inverse problems with a PDE describing the underlying physics. Our results support a negative answer to the question – i.e. decreasing h (increasing mesh density), increasing p, or simultaneously decreasing h and increasing p does not guarantee lower error reconstructions in applied inverse problems. Rather, there is complex balance between the accuracy of the hp-forward model, noise, prior knowledge (regularization), Jacobian accuracy, and ill-conditioning of the Jacobian matrix which ultimately contribute to reconstruction errors. As demonstrated herein, it is often more advantageous to use lower-order hp-forward models than higher-order hp-forward models in applied inverse problems. These realizations and other counterintuitive behavior of applied inverse problems using hp-forward models are described in detail herein.
-
Three dimensional through-wall imaging: Inverse scattering problems with an inhomogeneous background medium
Rui Chen,Zhun Wei,Xudong Chen
DOI: https://doi.org/10.1109/apcap.2015.7374465
2015-01-01
Abstract:This paper investigates the three dimensional inverse scattering problem with an inhomogeneous background medium where the known inhomogeneities are bounded in a finite domain. The forward scattering problem is solved using the finite element-boundary integral method (FEM-BI) while the inverse problem is oriented to minimize a cost function using nonlinear conjugate gradient optimization algorithm. The background Green's function at each discrete point on the surface of the computational domain can be numerically obtained using commercial software. Such an approach shrinks the domain of inversion to a much smaller region than the standard finite difference approach. The algorithm presented in this paper finds wide applications in nondestructive evaluation, such as through-wall imaging.
-
Stabilizing inverse problems by internal data. II. Non-local internal data and generic linearized uniqueness
Peter Kuchment,Dustin Steinhauer
DOI: https://doi.org/10.48550/arXiv.1407.0763
2014-11-03
Abstract:In the previous paper "Stabilizing Inverse Problems by Internal Data", 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 Frechet 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 "Stabilizing Inverse Problems by Internal Data" can be successfully adopted to this situation as well. We also obtain results on generic uniqueness for the linearized problem in a variety of situations, including those arising in acousto-electric and quantitative photoacoustic tomography.
Analysis of PDEs
-
On mathematical problems of two-coefficient inverse problems of ultrasonic tomography
Alexander V Goncharsky,Sergey Romanov,Sergey Seryozhnikov
DOI: https://doi.org/10.1088/1361-6420/ad2aa9
IF: 2.408
2024-02-20
Inverse Problems
Abstract:This paper proves the theorem of uniqueness for the solution of a coefficient inverse problem for the wave equation in R^2 with two unknown coefficients: speed of sound and absorption. The original nonlinear coefficient inverse problem is reduced to an equivalent system of two uniquely solvable linear integral equations of the first kind with respect to the sound speed and absorption coefficients. Estimates are made, substantiating the multistage method for two unknown coefficients. These estimates show that given sufficiently low frequencies and small inhomogeneities, the residual functional for the nonlinear inverse problem approaches a convex one. This solution method for nonlinear coefficient inverse problems is not linked to the limit approach as frequency tends to zero, but assumes solving the inverse problem using sufficiently low, but not zero, frequencies at the first stage. For small inhomogeneities that are typical, for instance, for medical tasks, carrying out real experiments at such frequencies does not present major difficulties. The capabilities of the method are demonstrated on a model inverse problem with unknown sound speed and absorption coefficients. The method effectively solves the nonlinear problem with parameter values typical for tomographic diagnostics of soft tissues in medicine. A resolution of approximately 2 mm was achieved using an average sounding pulse wavelength of 5 mm.
mathematics, applied,physics, mathematical
-
A Numerical Algorithm for Inverse Problems in Photothermal Detection
JC CHENG,SY ZHANG
DOI: https://doi.org/10.1051/jp4:19947131
1994-01-01
Abstract:An iterative algorithm for solving nonlinear inverse problems in the photothermal detection is presented. The basic idea of the algorithm is that the nonlinear inverse problem in the frequency domain is solved by a hybrid of Newton-like iterative method, pseudoinverse technique, and finite difference method. The numerical simulation shows the algorithm is effective for profile reconstruction of the thermal conductivity in the case of errorless data.
-
Inversion Problem for Radon Transforms Defined on Pseudoconvex Sets
DOI: https://doi.org/10.1134/s1064562424702004
2024-06-11
Doklady Mathematics
Abstract:Some questions concerning the inversion of the classical and generalized integral Radon transforms are discussed. The main issue is to determine information about the integrand if the values of some integrals are known. A feature of this work is that a function is integrated over hyperplanes in a finite-dimensional Euclidean space and the integrands depend not only on the variables of integration, but also on some of the variables characterizing the hyperplanes. The independent variables describing the known integrals are fewer than those in the unknown integrand. We consider discontinuous integrands defined on specifically introduced pseudoconvex sets. A Stefan-type problem of finding discontinuity surfaces of the integrand is posed. Formulas for solving the problem under study are derived by applying special integro-differential operators to known data.
mathematics