On quasiclassical approximation in the inverse scattering method

International Nuclear Information System (INIS)

Geogdzhaev, V.V.

1985-01-01

Using as an example quasiclassical limits of the Korteweg-de Vries equation and nonlinear Schroedinger equation, the quasiclassical limiting variant of the inverse scattering problem method is presented. In quasiclassical approximation the inverse scattering problem for the Schroedinger equation is reduced to the classical inverse scattering problem

On the inverse problem of dissipative scattering theory. 3

International Nuclear Information System (INIS)

Neidhardt, H.

1988-01-01

Considering a scattering theory in the class of contractions on Hilbert spaces one solves the inverse problem in an operaor-theoretical manner. The solution is obtained underthe very general assumptions that the free evolutions are different for different time directions that not only the perturbed or full evolutions but also the free evolutions are given by contractions. It is shown that the class of contractive Hankel operators can be viewed as a set of scattering operators. This implies the possibility that the scattering operator can be compact. Moreover, the result is applied to the so-called Lax-Phillips scattering theory with losses restoring a result of B.S. Pavlov on the completion of this theory in a quite different manner. 15 refs

A domain derivative-based method for solving elastodynamic inverse obstacle scattering problems

International Nuclear Information System (INIS)

Le Louër, Frédérique

2015-01-01

The present work is concerned with the shape reconstruction problem of isotropic elastic inclusions from far-field data obtained by the scattering of a finite number of time-harmonic incident plane waves. This paper aims at completing the theoretical framework which is necessary for the application of geometric optimization tools to the inverse transmission problem in elastodynamics. The forward problem is reduced to systems of boundary integral equations following the direct and indirect methods initially developed for solving acoustic transmission problems. We establish the Fréchet differentiability of the boundary to far-field operator and give a characterization of the first Fréchet derivative and its adjoint operator. Using these results we propose an inverse scattering algorithm based on the iteratively regularized Gauß–Newton method and show numerical experiments in the special case of star-shaped obstacles. (paper)

Multiple scattering processes: inverse and direct

International Nuclear Information System (INIS)

Kagiwada, H.H.; Kalaba, R.; Ueno, S.

1975-01-01

The purpose of the work is to formulate inverse problems in radiative transfer, to introduce the functions b and h as parameters of internal intensity in homogeneous slabs, and to derive initial value problems to replace the more traditional boundary value problems and integral equations of multiple scattering with high computational efficiency. The discussion covers multiple scattering processes in a one-dimensional medium; isotropic scattering in homogeneous slabs illuminated by parallel rays of radiation; the theory of functions b and h in homogeneous slabs illuminated by isotropic sources of radiation either at the top or at the bottom; inverse and direct problems of multiple scattering in slabs including internal sources; multiple scattering in inhomogeneous media, with particular reference to inverse problems for estimation of layers and total thickness of inhomogeneous slabs and to multiple scattering problems with Lambert's law and specular reflectors underlying slabs; and anisotropic scattering with reduction of the number of relevant arguments through axially symmetric fields and expansion in Legendre functions. Gaussian quadrature data for a seven point formula, a FORTRAN program for computing the functions b and h, and tables of these functions supplement the text

A direct sampling method to an inverse medium scattering problem

KAUST Repository

Ito, Kazufumi

2012-01-10

In this work we present a novel sampling method for time harmonic inverse medium scattering problems. It provides a simple tool to directly estimate the shape of the unknown scatterers (inhomogeneous media), and it is applicable even when the measured data are only available for one or two incident directions. A mathematical derivation is provided for its validation. Two- and three-dimensional numerical simulations are presented, which show that the method is accurate even with a few sets of scattered field data, computationally efficient, and very robust with respect to noises in the data. © 2012 IOP Publishing Ltd.

Quantum method of the inverse scattering problem. Pt. 1

International Nuclear Information System (INIS)

Sklyamin, E.K.; Takhtadzhyan, L.A.; Faddeev, L.D.

1978-12-01

In this work the authors use a formulation for the method of the inverse scattering problem for quantum-mechanical models of the field theory, that can be found in a quantization of these fully integrable systems. As the most important example serves the system (sinγ) 2 with the movement equation: γtt -γxx + m 2 /β sinβγ = 0 that is known under the specification Sine-Gordon-equation. (orig.) [de

On the quantum inverse scattering problem

International Nuclear Information System (INIS)

Maillet, J.M.; Terras, V.

2000-01-01

A general method for solving the so-called quantum inverse scattering problem (namely the reconstruction of local quantum (field) operators in term of the quantum monodromy matrix satisfying a Yang-Baxter quadratic algebra governed by an R-matrix) for a large class of lattice quantum integrable models is given. The principal requirement being the initial condition (R(0)=P, the permutation operator) for the quantum R-matrix solving the Yang-Baxter equation, it applies not only to most known integrable fundamental lattice models (such as Heisenberg spin chains) but also to lattice models with arbitrary number of impurities and to the so-called fused lattice models (including integrable higher spin generalizations of Heisenberg chains). Our method is then applied to several important examples like the sl n XXZ model, the XYZ spin-((1)/(2)) chain and also to the spin-s Heisenberg chains

One-dimensional scattering problem for inverse square potential

International Nuclear Information System (INIS)

Mineev, V.S.

1990-01-01

Analytical continuation of the solution for the Schroedinger equation of inverse square potential, together with the modified method for variation of constants makes it possible to construct admittable self-adjoint extensions and to completely analyze the respective scattering problem along the entire line. In this case, the current density conservation and the wave function continuity when passing through the singular point x=0 require, that a 8-shaped induced potential should be introduced in the Schroedinger equation. The relevant calculations have shown that the potential x -2 can be either absolutely penetrable or absolutely impenetrable. 16 refs

Multi-frequency direct sampling method in inverse scattering problem

Science.gov (United States)

Kang, Sangwoo; Lambert, Marc; Park, Won-Kwang

2017-10-01

We consider the direct sampling method (DSM) for the two-dimensional inverse scattering problem. Although DSM is fast, stable, and effective, some phenomena remain unexplained by the existing results. We show that the imaging function of the direct sampling method can be expressed by a Bessel function of order zero. We also clarify the previously unexplained imaging phenomena and suggest multi-frequency DSM to overcome traditional DSM. Our method is evaluated in simulation studies using both single and multiple frequencies.

Angle-domain inverse scattering migration/inversion in isotropic media

Science.gov (United States)

Li, Wuqun; Mao, Weijian; Li, Xuelei; Ouyang, Wei; Liang, Quan

2018-07-01

The classical seismic asymptotic inversion can be transformed into a problem of inversion of generalized Radon transform (GRT). In such methods, the combined parameters are linearly attached to the scattered wave-field by Born approximation and recovered by applying an inverse GRT operator to the scattered wave-field data. Typical GRT-style true-amplitude inversion procedure contains an amplitude compensation process after the weighted migration via dividing an illumination associated matrix whose elements are integrals of scattering angles. It is intuitional to some extent that performs the generalized linear inversion and the inversion of GRT together by this process for direct inversion. However, it is imprecise to carry out such operation when the illumination at the image point is limited, which easily leads to the inaccuracy and instability of the matrix. This paper formulates the GRT true-amplitude inversion framework in an angle-domain version, which naturally degrades the external integral term related to the illumination in the conventional case. We solve the linearized integral equation for combined parameters of different fixed scattering angle values. With this step, we obtain high-quality angle-domain common-image gathers (CIGs) in the migration loop which provide correct amplitude-versus-angle (AVA) behavior and reasonable illumination range for subsurface image points. Then we deal with the over-determined problem to solve each parameter in the combination by a standard optimization operation. The angle-domain GRT inversion method keeps away from calculating the inaccurate and unstable illumination matrix. Compared with the conventional method, the angle-domain method can obtain more accurate amplitude information and wider amplitude-preserved range. Several model tests demonstrate the effectiveness and practicability.

Uniqueness and numerical methods in inverse obstacle scattering

International Nuclear Information System (INIS)

Kress, Rainer

2007-01-01

The inverse problem we consider in this tutorial is to determine the shape of an obstacle from the knowledge of the far field pattern for scattering of time-harmonic plane waves. In the first part we will concentrate on the issue of uniqueness, i.e., we will investigate under what conditions an obstacle and its boundary condition can be identified from a knowledge of its far field pattern for incident plane waves. We will review some classical and some recent results and draw attention to open problems. In the second part we will survey on numerical methods for solving inverse obstacle scattering problems. Roughly speaking, these methods can be classified into three groups. Iterative methods interpret the inverse obstacle scattering problem as a nonlinear ill-posed operator equation and apply iterative schemes such as regularized Newton methods, Landweber iterations or conjugate gradient methods for its solution. Decomposition methods, in principle, separate the inverse scattering problem into an ill-posed linear problem to reconstruct the scattered wave from its far field and the subsequent determination of the boundary of the scatterer from the boundary condition. Finally, the third group consists of the more recently developed sampling methods. These are based on the numerical evaluation of criteria in terms of indicator functions that decide whether a point lies inside or outside the scatterer. The tutorial will give a survey by describing one or two representatives of each group including a discussion on the various advantages and disadvantages

Shrinkage-thresholding enhanced born iterative method for solving 2D inverse electromagnetic scattering problem

KAUST Repository

Desmal, Abdulla; Bagci, Hakan

2014-01-01

A numerical framework that incorporates recently developed iterative shrinkage thresholding (IST) algorithms within the Born iterative method (BIM) is proposed for solving the two-dimensional inverse electromagnetic scattering problem. IST

Inverse scattering with supersymmetric quantum mechanics

International Nuclear Information System (INIS)

Baye, Daniel; Sparenberg, Jean-Marc

2004-01-01

The application of supersymmetric quantum mechanics to the inverse scattering problem is reviewed. The main difference with standard treatments of the inverse problem lies in the simple and natural extension to potentials with singularities at the origin and with a Coulomb behaviour at infinity. The most general form of potentials which are phase-equivalent to a given potential is discussed. The use of singular potentials allows adding or removing states from the bound spectrum without contradicting the Levinson theorem. Physical applications of phase-equivalent potentials in nuclear reactions and in three-body systems are described. Derivation of a potential from the phase shift at fixed orbital momentum can also be performed with the supersymmetric inversion by using a Bargmann-type approximation of the scattering matrix or phase shift. A unique singular potential without bound states can be obtained from any phase shift. A limited number of bound states depending on the singularity can then be added. This inversion procedure is illustrated with nucleon-nucleon scattering

Direct sampling methods for inverse elastic scattering problems

Science.gov (United States)

Ji, Xia; Liu, Xiaodong; Xi, Yingxia

2018-03-01

We consider the inverse elastic scattering of incident plane compressional and shear waves from the knowledge of the far field patterns. Specifically, three direct sampling methods for location and shape reconstruction are proposed using the different component of the far field patterns. Only inner products are involved in the computation, thus the novel sampling methods are very simple and fast to be implemented. With the help of the factorization of the far field operator, we give a lower bound of the proposed indicator functionals for sampling points inside the scatterers. While for the sampling points outside the scatterers, we show that the indicator functionals decay like the Bessel functions as the sampling point goes away from the boundary of the scatterers. We also show that the proposed indicator functionals continuously dependent on the far field patterns, which further implies that the novel sampling methods are extremely stable with respect to data error. For the case when the observation directions are restricted into the limited aperture, we firstly introduce some data retrieval techniques to obtain those data that can not be measured directly and then use the proposed direct sampling methods for location and shape reconstructions. Finally, some numerical simulations in two dimensions are conducted with noisy data, and the results further verify the effectiveness and robustness of the proposed sampling methods, even for multiple multiscale cases and limited-aperture problems.

Using mixed data in the inverse scattering problem

International Nuclear Information System (INIS)

Lassaut, M.; Larsen, S.Y.; Sofianos, S.A.; Wallet, J.C.

2008-01-01

Consider the fixed-l inverse scattering problem. We show that the zeros of the regular solution of the Schroedinger equation, τ n (E), which are monotonic functions of the energy, determine a unique potential when the domain of the energy is such that the τ n (E) range from zero to infinity. This suggest that the use of the mixed data of phase-shifts (δ(l 0 , k),k ≥ k 0 ) set-theoretic union (δ(l,k 0 ),l ≥ l 0 ), for which the zeros of the regular solution are monotonic in both domains, and range from zero to infinity, offers the possibility of determining the potential in a unique way. (author)

On the solution of the inverse scattering problem for the quadratic bundle of the one-dimensional Schroedinger operators of the whole axis

International Nuclear Information System (INIS)

Maksudov, F.G.; Gusejnov, G.Sh.

1986-01-01

Inverse scattering problem for the quadratic bundle of the Schroedinger one-dimensional operators in the whole axis is solved. The problem solution is given on the assumption of the discrete spectrum absence. In the discrete spectrum presence the inverse scattering problem solution is known for the Shroedinger differential equation considered

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

International Nuclear Information System (INIS)

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

2015-01-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

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

Energy Technology Data Exchange (ETDEWEB)

Kılıç, Emre, E-mail: emre.kilic@tum.de; 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.

On the method of inverse scattering problem and Baecklund transformations for supersymmetric equations

International Nuclear Information System (INIS)

Chaichian, M.; Kulish, P. P.

1978-04-01

Supersymmetric Liouville and sine-Gordon equations are studied. We write down for these models the system of linear equations for which the method of inverse scattering problem should be applicable. Expressions for an infinite set of conserved currents are explicitly given. Supersymmetric Baecklund transformations and generalized conservation laws are constructed. (author)

The Mathematical Basis of the Inverse Scattering Problem for Cracks from Near-Field Data

Directory of Open Access Journals (Sweden)

Yao Mao

2015-01-01

Full Text Available We consider the acoustic scattering problem from a crack which has Dirichlet boundary condition on one side and impedance boundary condition on the other side. The inverse scattering problem in this paper tries to determine the shape of the crack and the surface impedance coefficient from the near-field measurements of the scattered waves, while the source point is placed on a closed curve. We firstly establish a near-field operator and focus on the operator’s mathematical analysis. Secondly, we obtain a uniqueness theorem for the shape and surface impedance. Finally, by using the operator’s properties and modified linear sampling method, we reconstruct the shape and surface impedance.

Direct and inverse scattering for viscoelastic media

International Nuclear Information System (INIS)

Ammicht, E.; Corones, J.P.; Krueger, R.J.

1987-01-01

A time domain approach to direct and inverse scattering problems for one-dimensional viscoelastic media is presented. Such media can be characterized as having a constitutive relation between stress and strain which involves the past history of the strain through a memory function, the relaxation modulus. In the approach in this article, the relaxation modulus of a material is shown to be related to the reflection properties of the material. This relation provides a constructive algorithm for direct and inverse scattering problems. A numerical implementation of this algorithm is tested on several problems involving realistic relaxation moduli

Inverse Scattering Problem For The Schrödinger Equation With An Additional Quadratic Potential On The Entire Axis

Science.gov (United States)

Guseinov, I. M.; Khanmamedov, A. Kh.; Mamedova, A. F.

2018-04-01

We consider the Schrödinger equation with an additional quadratic potential on the entire axis and use the transformation operator method to study the direct and inverse problems of the scattering theory. We obtain the main integral equations of the inverse problem and prove that the basic equations are uniquely solvable.

Point sources and multipoles in inverse scattering theory

CERN Document Server

Potthast, Roland

2001-01-01

Over the last twenty years, the growing availability of computing power has had an enormous impact on the classical fields of direct and inverse scattering. The study of inverse scattering, in particular, has developed rapidly with the ability to perform computational simulations of scattering processes and led to remarkable advances in a range of applications, from medical imaging and radar to remote sensing and seismic exploration. Point Sources and Multipoles in Inverse Scattering Theory provides a survey of recent developments in inverse acoustic and electromagnetic scattering theory. Focusing on methods developed over the last six years by Colton, Kirsch, and the author, this treatment uses point sources combined with several far-reaching techniques to obtain qualitative reconstruction methods. The author addresses questions of uniqueness, stability, and reconstructions for both two-and three-dimensional problems.With interest in extracting information about an object through scattered waves at an all-ti...

A two-stage method for inverse medium scattering

KAUST Repository

Ito, Kazufumi; Jin, Bangti; Zou, Jun

2013-01-01

We present a novel numerical method to the time-harmonic inverse medium scattering problem of recovering the refractive index from noisy near-field scattered data. The approach consists of two stages, one pruning step of detecting the scatterer

Indoor detection of passive targets recast as an inverse scattering problem

Science.gov (United States)

Gottardi, G.; Moriyama, T.

2017-10-01

The wireless local area networks represent an alternative to custom sensors and dedicated surveillance systems for target indoor detection. The availability of the channel state information has opened the exploitation of the spatial and frequency diversity given by the orthogonal frequency division multiplexing. Such a fine-grained information can be used to solve the detection problem as an inverse scattering problem. The goal of the detection is to reconstruct the properties of the investigation domain, namely to estimate if the domain is empty or occupied by targets, starting from the measurement of the electromagnetic perturbation of the wireless channel. An innovative inversion strategy exploiting both the frequency and the spatial diversity of the channel state information is proposed. The target-dependent features are identified combining the Kruskal-Wallis test and the principal component analysis. The experimental validation points out the detection performance of the proposed method when applied to an existing wireless link of a WiFi architecture deployed in a real indoor scenario. False detection rates lower than 2 [%] have been obtained.

Inverse and Ill-posed Problems Theory and Applications

CERN Document Server

Kabanikhin, S I

2011-01-01

The text demonstrates the methods for proving the existence (if et all) and finding of inverse and ill-posed problems solutions in linear algebra, integral and operator equations, integral geometry, spectral inverse problems, and inverse scattering problems. It is given comprehensive background material for linear ill-posed problems and for coefficient inverse problems for hyperbolic, parabolic, and elliptic equations. A lot of examples for inverse problems from physics, geophysics, biology, medicine, and other areas of application of mathematics are included.

A direct sampling method for inverse electromagnetic medium scattering

KAUST Repository

Ito, Kazufumi; Jin, Bangti; Zou, Jun

2013-01-01

In this paper, we study the inverse electromagnetic medium scattering problem of estimating the support and shape of medium scatterers from scattered electric/magnetic near-field data. We shall develop a novel direct sampling method based

Inverse radiative transfer problems in two-dimensional heterogeneous media

International Nuclear Information System (INIS)

Tito, Mariella Janette Berrocal

2001-01-01

The analysis of inverse problems in participating media where emission, absorption and scattering take place has several relevant applications in engineering and medicine. Some of the techniques developed for the solution of inverse problems have as a first step the solution of the direct problem. In this work the discrete ordinates method has been used for the solution of the linearized Boltzmann equation in two dimensional cartesian geometry. The Levenberg - Marquardt method has been used for the solution of the inverse problem of internal source and absorption and scattering coefficient estimation. (author)

Music algorithm for imaging of a sound-hard arc in limited-view inverse scattering problem

Science.gov (United States)

Park, Won-Kwang

2017-07-01

MUltiple SIgnal Classification (MUSIC) algorithm for a non-iterative imaging of sound-hard arc in limited-view inverse scattering problem is considered. In order to discover mathematical structure of MUSIC, we derive a relationship between MUSIC and an infinite series of Bessel functions of integer order. This structure enables us to examine some properties of MUSIC in limited-view problem. Numerical simulations are performed to support the identified structure of MUSIC.

Inverse problems and inverse scattering of plane waves

CERN Document Server

Ghosh Roy, Dilip N

2001-01-01

The purpose of this text is to present the theory and mathematics of inverse scattering, in a simple way, to the many researchers and professionals who use it in their everyday research. While applications range across a broad spectrum of disciplines, examples in this text will focus primarly, but not exclusively, on acoustics. The text will be especially valuable for those applied workers who would like to delve more deeply into the fundamentally mathematical character of the subject matter.Practitioners in this field comprise applied physicists, engineers, and technologists, whereas the theory is almost entirely in the domain of abstract mathematics. This gulf between the two, if bridged, can only lead to improvement in the level of scholarship in this highly important discipline. This is the book''s primary focus.

Discrete inverse scattering theory and the continuum limit

International Nuclear Information System (INIS)

Berryman, J.G.; Greene, R.R.

1978-01-01

The class of satisfactory difference approximations for the Schroedinger equation in discrete inverse scattering theory is shown smaller than previously supposed. A fast algorithm (analogous to the Levinson algorithm for Toeplitz matrices) is found for solving the discrete inverse problem. (Auth.)

Solution of the inverse scattering problem at fixed energy with non-physical S matrix elements

International Nuclear Information System (INIS)

Eberspaecher, M.; Amos, K.; Apagyi, B.

1999-12-01

The quantum mechanical inverse elastic scattering problem is solved with the modified Newton-Sabatier method. A set of S matrix elements calculated from a realistic analytic optical model potential serves as input data. It is demonstrated that the quality of the inversion potential can be improved by including non-physical S matrix elements to half, quarter and eighth valued partial waves if the original set does not contain enough information to determine the interaction potential. We demonstrate that results can be very sensitive to the choice of those non-physical S matrix values both with the analytic potential model and in a real application in which the experimental cross section for the symmetrical scattering system of 12 C+ 12 C at E=7.998 MeV is analyzed

Solution of the Cox-Thompson inverse scattering problem using finite set of phase shifts

CERN Document Server

Apagyi, B; Scheid, W

2003-01-01

A system of nonlinear equations is presented for the solution of the Cox-Thompson inverse scattering problem (1970 J. Math. Phys. 11 805) at fixed energy. From a given finite set of phase shifts for physical angular momenta, the nonlinear equations determine related sets of asymptotic normalization constants and nonphysical (shifted) angular momenta from which all quantities of interest, including the inversion potential itself, can be calculated. As a first application of the method we use input data consisting of a finite set of phase shifts calculated from Woods-Saxon and box potentials representing interactions with diffuse or sharp surfaces, respectively. The results for the inversion potentials, their first moments and asymptotic properties are compared with those provided by the Newton-Sabatier quantum inversion procedure. It is found that in order to achieve inversion potentials of similar quality, the Cox-Thompson method requires a smaller set of phase shifts than the Newton-Sabatier procedure.

Solution of the Cox-Thompson inverse scattering problem using finite set of phase shifts

International Nuclear Information System (INIS)

Apagyi, Barnabas; Harman, Zoltan; Scheid, Werner

2003-01-01

A system of nonlinear equations is presented for the solution of the Cox-Thompson inverse scattering problem (1970 J. Math. Phys. 11 805) at fixed energy. From a given finite set of phase shifts for physical angular momenta, the nonlinear equations determine related sets of asymptotic normalization constants and nonphysical (shifted) angular momenta from which all quantities of interest, including the inversion potential itself, can be calculated. As a first application of the method we use input data consisting of a finite set of phase shifts calculated from Woods-Saxon and box potentials representing interactions with diffuse or sharp surfaces, respectively. The results for the inversion potentials, their first moments and asymptotic properties are compared with those provided by the Newton-Sabatier quantum inversion procedure. It is found that in order to achieve inversion potentials of similar quality, the Cox-Thompson method requires a smaller set of phase shifts than the Newton-Sabatier procedure

Inverse Scattering in a Multipath Environment

Directory of Open Access Journals (Sweden)

A. Cuccaro

2016-09-01

Full Text Available In this contribution an inverse scattering problem is ad- dressed in a multipath environment. In particular, multipath is created by known ”extra” point-like scatterers (passive elements expressely deployed between the scene under in- vestigation and the source/measurement domains. Through a back-projection imaging scheme, the role of the passive elements on the achievable performance is shown and com- pared to the free-space case.

Prize for Industrial Applications of Physics Talk: The Inverse Scattering Problem and the role of measurements in its solution

Science.gov (United States)

Wyatt, Philip

2009-03-01

The electromagnetic inverse scattering problem suggests that if a homogeneous and non-absorbing object be illuminated with a monochromatic light source and if the far field scattered light intensity is known at sufficient scattering angles, then, in principle, one could derive the dielectric structure of the scattering object. In general, this is an ill-posed problem and methods must be developed to regularize the search for unique solutions. An iterative procedure often begins with a model of the scattering object, solves the forward scattering problem using this model, and then compares these calculated results with the measured values. Key to any such solution is instrumentation capable of providing adequate data. To this end, the development of the first laser based absolute light scattering photometers is described together with their continuing evolution and some of the remarkable discoveries made with them. For particles much smaller than the wavelength of the incident light (e.g. macromolecules), the inverse scattering problems are easily solved. Among the many solutions derived with this instrumentation are the in situ structure of bacterial cells, new drug delivery mechanisms, the development of new vaccines and other biologicals, characterization of wines, the possibility of custom chemotherapy, development of new polymeric materials, identification of protein crystallization conditions, and a variety discoveries concerning protein interactions. A new form of the problem is described to address bioterrorist threats. Over the many years of development and refinement, one element stands out as essential for the successes that followed: the R and D teams were always directed and executed by physics trained theorists and experimentalists. 14 Ph. D. physicists each made his/her unique contribution to the development of these evolving instruments and the interpretation of their results.

A necessary condition for applying MUSIC algorithm in limited-view inverse scattering problem

Science.gov (United States)

Park, Taehoon; Park, Won-Kwang

2015-09-01

Throughout various results of numerical simulations, it is well-known that MUltiple SIgnal Classification (MUSIC) algorithm can be applied in the limited-view inverse scattering problems. However, the application is somehow heuristic. In this contribution, we identify a necessary condition of MUSIC for imaging of collection of small, perfectly conducting cracks. This is based on the fact that MUSIC imaging functional can be represented as an infinite series of Bessel function of integer order of the first kind. Numerical experiments from noisy synthetic data supports our investigation.

A necessary condition for applying MUSIC algorithm in limited-view inverse scattering problem

International Nuclear Information System (INIS)

Park, Taehoon; Park, Won-Kwang

2015-01-01

Throughout various results of numerical simulations, it is well-known that MUltiple SIgnal Classification (MUSIC) algorithm can be applied in the limited-view inverse scattering problems. However, the application is somehow heuristic. In this contribution, we identify a necessary condition of MUSIC for imaging of collection of small, perfectly conducting cracks. This is based on the fact that MUSIC imaging functional can be represented as an infinite series of Bessel function of integer order of the first kind. Numerical experiments from noisy synthetic data supports our investigation. (paper)

Inverse problems in linear transport theory

International Nuclear Information System (INIS)

Dressler, K.

1988-01-01

Inverse problems for a class of linear kinetic equations are investigated. The aim is to identify the scattering kernel of a transport equation (corresponding to the structure of a background medium) by observing the 'albedo' part of the solution operator for the corresponding direct initial boundary value problem. This means to get information on some integral operator in an integrodifferential equation through on overdetermined boundary value problem. We first derive a constructive method for solving direct halfspace problems and prove a new factorization theorem for the solutions. Using this result we investigate stationary inverse problems with respect to well posedness (e.g. reduce them to classical ill-posed problems, such as integral equations of first kind). In the time-dependent case we show that a quite general inverse problem is well posed and solve it constructively. (orig.)

Inverse radiation problem of temperature distribution in one-dimensional isotropically scattering participating slab with variable refractive index

International Nuclear Information System (INIS)

Namjoo, A.; Sarvari, S.M. Hosseini; Behzadmehr, A.; Mansouri, S.H.

2009-01-01

In this paper, an inverse analysis is performed for estimation of source term distribution from the measured exit radiation intensities at the boundary surfaces in a one-dimensional absorbing, emitting and isotropically scattering medium between two parallel plates with variable refractive index. The variation of refractive index is assumed to be linear. The radiative transfer equation is solved by the constant quadrature discrete ordinate method. The inverse problem is formulated as an optimization problem for minimizing an objective function which is expressed as the sum of square deviations between measured and estimated exit radiation intensities at boundary surfaces. The conjugate gradient method is used to solve the inverse problem through an iterative procedure. The effects of various variables on source estimation are investigated such as type of source function, errors in the measured data and system parameters, gradient of refractive index across the medium, optical thickness, single scattering albedo and boundary emissivities. The results show that in the case of noisy input data, variation of system parameters may affect the inverse solution, especially at high error values in the measured data. The error in measured data plays more important role than the error in radiative system parameters except the refractive index distribution; however the accuracy of source estimation is very sensitive toward error in refractive index distribution. Therefore, refractive index distribution and measured exit intensities should be measured accurately with a limited error bound, in order to have an accurate estimation of source term in a graded index medium.

Solution of the nonlinear inverse scattering problem by T-matrix completion. I. Theory.

Science.gov (United States)

Levinson, Howard W; Markel, Vadim A

2016-10-01

We propose a conceptually different method for solving nonlinear inverse scattering problems (ISPs) such as are commonly encountered in tomographic ultrasound imaging, seismology, and other applications. The method is inspired by the theory of nonlocality of physical interactions and utilizes the relevant formalism. We formulate the ISP as a problem whose goal is to determine an unknown interaction potential V from external scattering data. Although we seek a local (diagonally dominated) V as the solution to the posed problem, we allow V to be nonlocal at the intermediate stages of iterations. This allows us to utilize the one-to-one correspondence between V and the T matrix of the problem. Here it is important to realize that not every T corresponds to a diagonal V and we, therefore, relax the usual condition of strict diagonality (locality) of V. An iterative algorithm is proposed in which we seek T that is (i) compatible with the measured scattering data and (ii) corresponds to an interaction potential V that is as diagonally dominated as possible. We refer to this algorithm as to the data-compatible T-matrix completion. This paper is Part I in a two-part series and contains theory only. Numerical examples of image reconstruction in a strongly nonlinear regime are given in Part II [H. W. Levinson and V. A. Markel, Phys. Rev. E 94, 043318 (2016)10.1103/PhysRevE.94.043318]. The method described in this paper is particularly well suited for very large data sets that become increasingly available with the use of modern measurement techniques and instrumentation.

An analytical approach to estimate the number of small scatterers in 2D inverse scattering problems

International Nuclear Information System (INIS)

Fazli, Roohallah; Nakhkash, Mansor

2012-01-01

This paper presents an analytical method to estimate the location and number of actual small targets in 2D inverse scattering problems. This method is motivated from the exact maximum likelihood estimation of signal parameters in white Gaussian noise for the linear data model. In the first stage, the method uses the MUSIC algorithm to acquire all possible target locations and in the next stage, it employs an analytical formula that works as a spatial filter to determine which target locations are associated to the actual ones. The ability of the method is examined for both the Born and multiple scattering cases and for the cases of well-resolved and non-resolved targets. Many numerical simulations using both the coincident and non-coincident arrays demonstrate that the proposed method can detect the number of actual targets even in the case of very noisy data and when the targets are closely located. Using the experimental microwave data sets, we further show that this method is successful in specifying the number of small inclusions. (paper)

Applications of inverse and algebraic scattering theories

Energy Technology Data Exchange (ETDEWEB)

Amos, K. [Qinghua Univ., Beijing, BJ (China). Dept. of Physics

1997-06-01

Inverse scattering theories, algebraic scattering theory and exactly solvable scattering potentials are diverse ways by which scattering potentials can be defined from S-functions specified by fits to fixed energy, quantal scattering data. Applications have been made in nuclear (heavy ion and nucleon-nucleus scattering), atomic and molecular (electron scattering from simple molecules) systems. Three inverse scattering approaches are considered in detail; the semiclassical WKB and fully quantal Lipperheide-Fiedeldey method, than algebraic scattering theory is applied to heavy ion scattering and finally the exactly solvable Ginocchio potentials. Some nuclear results are ambiguous but the atomic and molecular inversion potentials are in good agreement with postulated forms. 21 refs., 12 figs.

Lectures on the inverse scattering method

International Nuclear Information System (INIS)

Zakharov, V.E.

1983-06-01

In a series of six lectures an elementary introduction to the theory of inverse scattering is given. The first four lectures contain a detailed theory of solitons in the framework of the KdV equation, together with the inverse scattering theory of the one-dimensional Schroedinger equation. In the fifth lecture the dressing method is described, while the sixth lecture gives a brief review of the equations soluble by the inverse scattering method. (author)

Estimation of biological parameters of marine organisms using linear and nonlinear acoustic scattering model-based inversion methods.

Science.gov (United States)

Chu, Dezhang; Lawson, Gareth L; Wiebe, Peter H

2016-05-01

The linear inversion commonly used in fisheries and zooplankton acoustics assumes a constant inversion kernel and ignores the uncertainties associated with the shape and behavior of the scattering targets, as well as other relevant animal parameters. Here, errors of the linear inversion due to uncertainty associated with the inversion kernel are quantified. A scattering model-based nonlinear inversion method is presented that takes into account the nonlinearity of the inverse problem and is able to estimate simultaneously animal abundance and the parameters associated with the scattering model inherent to the kernel. It uses sophisticated scattering models to estimate first, the abundance, and second, the relevant shape and behavioral parameters of the target organisms. Numerical simulations demonstrate that the abundance, size, and behavior (tilt angle) parameters of marine animals (fish or zooplankton) can be accurately inferred from the inversion by using multi-frequency acoustic data. The influence of the singularity and uncertainty in the inversion kernel on the inversion results can be mitigated by examining the singular values for linear inverse problems and employing a non-linear inversion involving a scattering model-based kernel.

On increasing stability in the two dimensional inverse source scattering problem with many frequencies

Science.gov (United States)

Entekhabi, Mozhgan Nora; Isakov, Victor

2018-05-01

In this paper, we will study the increasing stability in the inverse source problem for the Helmholtz equation in the plane when the source term is assumed to be compactly supported in a bounded domain Ω with a sufficiently smooth boundary. Using the Fourier transform in the frequency domain, bounds for the Hankel functions and for scattering solutions in the complex plane, improving bounds for the analytic continuation, and the exact observability for the wave equation led us to our goals which are a sharp uniqueness and increasing stability estimate when the wave number interval is growing.

Inverse scattering transform for the vector nonlinear Schroedinger equation with nonvanishing boundary conditions

International Nuclear Information System (INIS)

Prinari, Barbara; Ablowitz, Mark J.; Biondini, Gino

2006-01-01

The inverse scattering transform for the vector defocusing nonlinear Schroedinger (NLS) equation with nonvanishing boundary values at infinity is constructed. The direct scattering problem is formulated on a two-sheeted covering of the complex plane. Two out of the six Jost eigenfunctions, however, do not admit an analytic extension on either sheet of the Riemann surface. Therefore, a suitable modification of both the direct and the inverse problem formulations is necessary. On the direct side, this is accomplished by constructing two additional analytic eigenfunctions which are expressed in terms of the adjoint eigenfunctions. The discrete spectrum, bound states and symmetries of the direct problem are then discussed. In the most general situation, a discrete eigenvalue corresponds to a quartet of zeros (poles) of certain scattering data. The inverse scattering problem is formulated in terms of a generalized Riemann-Hilbert (RH) problem in the upper/lower half planes of a suitable uniformization variable. Special soliton solutions are constructed from the poles in the RH problem, and include dark-dark soliton solutions, which have dark solitonic behavior in both components, as well as dark-bright soliton solutions, which have one dark and one bright component. The linear limit is obtained from the RH problem and is shown to correspond to the Fourier transform solution obtained from the linearized vector NLS system

Mittag–Leffler's function, Vekua transform and an inverse obstacle scattering problem

International Nuclear Information System (INIS)

Ikehata, Masaru

2010-01-01

This paper studies a prototype of inverse obstacle scattering problems whose governing equation is the Helmholtz equation in two dimensions. An explicit method to extract information about the location and shape of unknown obstacles from the far-field operator with a fixed wave number is given. The method is based on an explicit construction of a modification of Mittag–Leffler's function via the Vekua transform and the study of the asymptotic behaviour; an explicit density in the Herglotz wavefunction that approximates the modification of Mittag–Leffler's function in the bounded domain surrounding unknown obstacles; a system of inequalities derived from Kirsch's factorization formula of the far-field operator. Then an indicator function which can be calculated from the far-field operator acting on the density is introduced. It is shown that the asymptotic behaviour of the indicator function yields information about the visible part of the exterior of the obstacles

Scattering angle base filtering of the inversion gradients

KAUST Repository

Alkhalifah, Tariq Ali

2014-01-01

Full waveform inversion (FWI) requires a hierarchical approach based on the availability of low frequencies to maneuver the complex nonlinearity associated with the problem of velocity inversion. I develop a model gradient filter to help us access the parts of the gradient more suitable to combat this potential nonlinearity. The filter is based on representing the gradient in the time-lag normalized domain, in which low scattering angles of the gradient update are initially muted. The result are long-wavelength updates controlled by the ray component of the wavefield. In this case, even 10 Hz data can produce near zero wavelength updates suitable for a background correction of the model. Allowing smaller scattering angle to contribute provides higher resolution information to the model.

Multi-scattering inversion for low model wavenumbers

KAUST Repository

Alkhalifah, Tariq Ali

2015-08-19

A successful full wavenumber inversion (FWI) implementation updates the low wavenumber model components first for proper wavefield propagation description, and slowly adds the high-wavenumber potentially scattering parts of the model. The low-wavenumber components can be extracted from the transmission parts of the recorded data given by direct arrivals or the transmission parts of the single and double-scattering wave-fields developed from a predicted scatter field. We develop a combined inversion of data modeled from the source and those corresponding to single and double scattering to update both the velocity model and the component of the velocity (perturbation) responsible for the single and double scattering. The combined inversion helps us access most of the potential model wavenumber information that may be embedded in the data. A scattering angle filter is used to divide the gradient of the combined inversion so initially the high wavenumber (low scattering angle) components of the gradient is directed to the perturbation model and the low wavenumber (high scattering angle) components to the velocity model. As our background velocity matures, the scattering angle divide is slowly lowered to allow for more of the higher wavenumbers to contribute the velocity model.

Single- and coupled-channel radial inverse scattering with supersymmetric transformations

International Nuclear Information System (INIS)

Baye, Daniel; Sparenberg, Jean-Marc; Pupasov-Maksimov, Andrey M; Samsonov, Boris F

2014-01-01

The present status of the three-dimensional inverse-scattering method with supersymmetric transformations is reviewed for the coupled-channel case. We first revisit in a pedagogical way the single-channel case, where the supersymmetric approach is shown to provide a complete, efficient and elegant solution to the inverse-scattering problem for the radial Schrödinger equation with short-range interactions. A special emphasis is put on the differences between conservative and non-conservative transformations, i.e. transformations that do or do not conserve the behaviour of solutions of the radial Schrödinger equation at the origin. In particular, we show that for the zero initial potential, a non-conservative transformation is always equivalent to a pair of conservative transformations. These single-channel results are illustrated on the inversion of the neutron–proton triplet eigenphase shifts for the S- and D-waves. We then summarize and extend our previous works on the coupled-channel case, i.e. on systems of coupled radial Schrödinger equations, and stress remaining difficulties and open questions of this problem by putting it in perspective with the single-channel case. We mostly concentrate on two-channel examples to illustrate general principles while keeping mathematics as simple as possible. In particular, we discuss the important difference between the equal-threshold and different-threshold problems. For equal thresholds, conservative transformations can provide non-diagonal Jost and scattering matrices. Iterations of such transformations in the two-channel case are studied and shown to lead to practical algorithms for inversion. A convenient particular technique where the mixing parameter can be fitted without modifying the eigenphases is developed with iterations of pairs of conjugate transformations. This technique is applied to the neutron–proton triplet S–D scattering matrix, for which exactly-solvable matrix potential models are constructed

Uniqueness in inverse elastic scattering with finitely many incident waves

International Nuclear Information System (INIS)

Elschner, Johannes; Yamamoto, Masahiro

2009-01-01

We consider the third and fourth exterior boundary value problems of linear isotropic elasticity and present uniqueness results for the corresponding inverse scattering problems with polyhedral-type obstacles and a finite number of incident plane elastic waves. Our approach is based on a reflection principle for the Navier equation. (orig.)

Semi-analytic equations to the Cox-Thompson inverse scattering method at fixed energy for special cases

International Nuclear Information System (INIS)

Palmai, T.; Apagyi, B.; Horvath, M.

2008-01-01

Solution of the Cox-Thompson inverse scattering problem at fixed energy 1-3 is reformulated resulting in semi-analytic equations. The new set of equations for the normalization constants and the nonphysical (shifted) angular momenta are free of matrix inversion operations. This simplification is a result of treating only the input phase shifts of partial waves of a given parity. Therefore, the proposed method can be applied for identical particle scattering of the bosonic type (or for certain cases of identical fermionic scattering). The new formulae are expected to be numerically more efficient than the previous ones. Based on the semi-analytic equations an approximate method is proposed for the generic inverse scattering problem, when partial waves of arbitrary parity are considered. (author)

Subspace-based optimization method for inverse scattering problems with an inhomogeneous background medium

International Nuclear Information System (INIS)

Chen, Xudong

2010-01-01

This paper proposes a version of the subspace-based optimization method to solve the inverse scattering problem with an inhomogeneous background medium where the known inhomogeneities are bounded in a finite domain. Although the background Green's function at each discrete point in the computational domain is not directly available in an inhomogeneous background scenario, the paper uses the finite element method to simultaneously obtain the Green's function at all discrete points. The essence of the subspace-based optimization method is that part of the contrast source is determined from the spectrum analysis without using any optimization, whereas the orthogonally complementary part is determined by solving a lower dimension optimization problem. This feature significantly speeds up the convergence of the algorithm and at the same time makes it robust against noise. Numerical simulations illustrate the efficacy of the proposed algorithm. The algorithm presented in this paper finds wide applications in nondestructive evaluation, such as through-wall imaging

A two-stage method for inverse medium scattering

KAUST Repository

Ito, Kazufumi

2013-03-01

We present a novel numerical method to the time-harmonic inverse medium scattering problem of recovering the refractive index from noisy near-field scattered data. The approach consists of two stages, one pruning step of detecting the scatterer support, and one resolution enhancing step with nonsmooth mixed regularization. The first step is strictly direct and of sampling type, and it faithfully detects the scatterer support. The second step is an innovative application of nonsmooth mixed regularization, and it accurately resolves the scatterer size as well as intensities. The nonsmooth model can be efficiently solved by a semi-smooth Newton-type method. Numerical results for two- and three-dimensional examples indicate that the new approach is accurate, computationally efficient, and robust with respect to data noise. © 2012 Elsevier Inc.

An inverse boundary value problem for the Schroedinger operator with vector potentials in two dimensions

International Nuclear Information System (INIS)

Ziqi Sun

1993-01-01

During the past few years a considerable interest has been focused on the inverse boundary value problem for the Schroedinger operator with a scalar (electric) potential. The popularity gained by this subject seems to be due to its connection with the inverse scattering problem at fixed energy, the inverse conductivity problem and other important inverse problems. This paper deals with an inverse boundary value problem for the Schroedinger operator with vector (electric and magnetic) potentials. As in the case of the scalar potential, results of this study would have immediate consequences in the inverse scattering problem for magnetic field at fixed energy. On the other hand, inverse boundary value problems for elliptic operators are of independent interest. The study is partly devoted to the understanding of the inverse boundary value problem for a class of general elliptic operator of second order. Note that a self-adjoint elliptic operator of second order with Δ as its principal symbol can always be written as a Schroedinger operator with vector potentials

Relativistic Inverse Scattering Problem for a Superposition of a Nonlocal Separable and a Local Quasipotential

International Nuclear Information System (INIS)

Chernichenko, Yu.D.

2005-01-01

Within the relativistic quasipotential approach to quantum field theory, the relativistic inverse scattering problem is solved for the case where the total quasipotential describing the interaction of two relativistic spinless particles having different masses is a superposition of a nonlocal separable and a local quasipotential. It is assumed that the local component of the total quasipotential is known and that there exist bound states in this local component. It is shown that the nonlocal separable component of the total interaction can be reconstructed provided that the local component, an increment of the phase shift, and the energies of bound states are known

A direct sampling method for inverse electromagnetic medium scattering

KAUST Repository

Ito, Kazufumi

2013-09-01

In this paper, we study the inverse electromagnetic medium scattering problem of estimating the support and shape of medium scatterers from scattered electric/magnetic near-field data. We shall develop a novel direct sampling method based on an analysis of electromagnetic scattering and the behavior of the fundamental solution. It is applicable to a few incident fields and needs only to compute inner products of the measured scattered field with the fundamental solutions located at sampling points. Hence, it is strictly direct, computationally very efficient and highly robust to the presence of data noise. Two- and three-dimensional numerical experiments indicate that it can provide reliable support estimates for multiple scatterers in the case of both exact and highly noisy data. © 2013 IOP Publishing Ltd.

Support minimized inversion of acoustic and elastic wave scattering

International Nuclear Information System (INIS)

Safaeinili, A.

1994-01-01

This report discusses the following topics on support minimized inversion of acoustic and elastic wave scattering: Minimum support inversion; forward modelling of elastodynamic wave scattering; minimum support linearized acoustic inversion; support minimized nonlinear acoustic inversion without absolute phase; and support minimized nonlinear elastic inversion

Inverse scattering scheme for the Dirac equation at fixed energy

International Nuclear Information System (INIS)

Leeb, H.; Lehninger, H.; Schilder, C.

2001-01-01

Full text: Based on the concept of generalized transformation operators a new hierarchy of Dirac equations with spherical symmetric scalar and fourth component vector potentials is presented. Within this hierarchy closed form expressions for the solutions, the potentials and the S-matrix can be given in terms of solutions of the original Dirac equation. Using these transformations an inverse scattering scheme has been constructed for the Dirac equation which is the analog to the rational scheme in the non-relativistic case. The given method provides for the first time an inversion scheme with closed form expressions for the S-matrix for non-relativistic scattering problems with central and spin-orbit potentials. (author)

Direct and inverse problems of studying the properties of multilayer nanostructures based on a two-dimensional model of X-ray reflection and scattering

Science.gov (United States)

Khachaturov, R. V.

2014-06-01

A mathematical model of X-ray reflection and scattering by multilayered nanostructures in the quasi-optical approximation is proposed. X-ray propagation and the electric field distribution inside the multilayered structure are considered with allowance for refraction, which is taken into account via the second derivative with respect to the depth of the structure. This model is used to demonstrate the possibility of solving inverse problems in order to determine the characteristics of irregularities not only over the depth (as in the one-dimensional problem) but also over the length of the structure. An approximate combinatorial method for system decomposition and composition is proposed for solving the inverse problems.

Inverse scattering transform for the time dependent Schroedinger equation with applications to the KPI equation

Energy Technology Data Exchange (ETDEWEB)

Xin, Zhou [Wisconsin Univ., Madison (USA). Dept. of Mathematics

1990-03-01

For the direct-inverse scattering transform of the time dependent Schroedinger equation, rigorous results are obtained based on an operator-triangular-factorization approach. By viewing the equation as a first order operator equation, similar results as for the first order n x n matrix system are obtained. The nonlocal Riemann-Hilbert problem for inverse scattering is shown to have solution. (orig.).

Inverse scattering transform for the time dependent Schroedinger equation with applications to the KPI equation

International Nuclear Information System (INIS)

Zhou Xin

1990-01-01

For the direct-inverse scattering transform of the time dependent Schroedinger equation, rigorous results are obtained based on an operator-triangular-factorization approach. By viewing the equation as a first order operator equation, similar results as for the first order n x n matrix system are obtained. The nonlocal Riemann-Hilbert problem for inverse scattering is shown to have solution. (orig.)

Inverse problem in nuclear physics

International Nuclear Information System (INIS)

Zakhariev, B.N.

1976-01-01

The method of reconstruction of interaction from the scattering data is formulated in the frame of the R-matrix theory in which the potential is determined by position of resonance Esub(lambda) and their reduced widths γ 2 lambda. In finite difference approximation for the Schroedinger equation this new approach allows to make the logics of the inverse problem IP more clear. A possibility of applications of IP formalism to various nuclear systems is discussed. (author)

Inverse random source scattering for the Helmholtz equation in inhomogeneous media

Science.gov (United States)

Li, Ming; Chen, Chuchu; Li, Peijun

2018-01-01

This paper is concerned with an inverse random source scattering problem in an inhomogeneous background medium. The wave propagation is modeled by the stochastic Helmholtz equation with the source driven by additive white noise. The goal is to reconstruct the statistical properties of the random source such as the mean and variance from the boundary measurement of the radiated random wave field at multiple frequencies. Both the direct and inverse problems are considered. We show that the direct problem has a unique mild solution by a constructive proof. For the inverse problem, we derive Fredholm integral equations, which connect the boundary measurement of the radiated wave field with the unknown source function. A regularized block Kaczmarz method is developed to solve the ill-posed integral equations. Numerical experiments are included to demonstrate the effectiveness of the proposed method.

Analysis of MUSIC-type imaging functional for single, thin electromagnetic inhomogeneity in limited-view inverse scattering problem

Science.gov (United States)

Ahn, Chi Young; Jeon, Kiwan; Park, Won-Kwang

2015-06-01

This study analyzes the well-known MUltiple SIgnal Classification (MUSIC) algorithm to identify unknown support of thin penetrable electromagnetic inhomogeneity from scattered field data collected within the so-called multi-static response matrix in limited-view inverse scattering problems. The mathematical theories of MUSIC are partially discovered, e.g., in the full-view problem, for an unknown target of dielectric contrast or a perfectly conducting crack with the Dirichlet boundary condition (Transverse Magnetic-TM polarization) and so on. Hence, we perform further research to analyze the MUSIC-type imaging functional and to certify some well-known but theoretically unexplained phenomena. For this purpose, we establish a relationship between the MUSIC imaging functional and an infinite series of Bessel functions of integer order of the first kind. This relationship is based on the rigorous asymptotic expansion formula in the existence of a thin inhomogeneity with a smooth supporting curve. Various results of numerical simulation are presented in order to support the identified structure of MUSIC. Although a priori information of the target is needed, we suggest a least condition of range of incident and observation directions to apply MUSIC in the limited-view problem.

A Riemann-Hilbert approach to the inverse problem for the Stark operator on the line

Science.gov (United States)

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.

Inversion of electron-water elastic scattering data

International Nuclear Information System (INIS)

Lun, A.; Chen, X.J.; Allen, L.J.; Amos, K.

1994-01-01

Fixed energy inverse scattering theory has been used to analyse the differential cross-sections for the elastic scattering of electrons from water molecules. Both semiclassical (WKB) and fully quantal inversion methods have been used with data taken in the energy range 100 to 1000 eV. Constrained to be real, the local inversion potentials are found to be energy dependent; a dependence that can be interpreted as the local equivalence of true nonlocality in the actual interaction. 14 refs., 4 tabs., 8 figs

Shrinkage-thresholding enhanced born iterative method for solving 2D inverse electromagnetic scattering problem

KAUST Repository

Desmal, Abdulla

2014-07-01

A numerical framework that incorporates recently developed iterative shrinkage thresholding (IST) algorithms within the Born iterative method (BIM) is proposed for solving the two-dimensional inverse electromagnetic scattering problem. IST algorithms minimize a cost function weighted between measurement-data misfit and a zeroth/first-norm penalty term and therefore promote "sharpness" in the solution. Consequently, when applied to domains with sharp variations, discontinuities, or sparse content, the proposed framework is more efficient and accurate than the "classical" BIM that minimizes a cost function with a second-norm penalty term. Indeed, numerical results demonstrate the superiority of the IST-BIM over the classical BIM when they are applied to sparse domains: Permittivity and conductivity profiles recovered using the IST-BIM are sharper and more accurate and converge faster. © 1963-2012 IEEE.

Inverse radiative transfer problems in two-dimensional heterogeneous media; Problemas inversos em transferencia radiativa em meios heterogeneos bidimensionais

Energy Technology Data Exchange (ETDEWEB)

Tito, Mariella Janette Berrocal

2001-01-01

The analysis of inverse problems in participating media where emission, absorption and scattering take place has several relevant applications in engineering and medicine. Some of the techniques developed for the solution of inverse problems have as a first step the solution of the direct problem. In this work the discrete ordinates method has been used for the solution of the linearized Boltzmann equation in two dimensional cartesian geometry. The Levenberg - Marquardt method has been used for the solution of the inverse problem of internal source and absorption and scattering coefficient estimation. (author)

Uniqueness of inverse scattering problem in local quantum physics

Energy Technology Data Exchange (ETDEWEB)

Schroer, Bert [Centro Brasileiro de Pesquisas Fisicas (CBPF), Rio de Janeiro, RJ (Brazil)]. E-mail: schroer@cbpf.br

2001-06-01

It is shown that the a Bisognano-Wichmann-Unruh inspired formulation of local quantum physics which starts from wedge-localized algebras, leads to a uniqueness proof for the scattering problem. The important mathematical tool is the thermal KMS aspect of localization and its strengthening by the requirement of crossing symmetry for generalized formfactors. (author)

Inverse feasibility problems of the inverse maximum flow problems

Indian Academy of Sciences (India)

199–209. c Indian Academy of Sciences. Inverse feasibility problems of the inverse maximum flow problems. ADRIAN DEACONU. ∗ and ELEONOR CIUREA. Department of Mathematics and Computer Science, Faculty of Mathematics and Informatics, Transilvania University of Brasov, Brasov, Iuliu Maniu st. 50,. Romania.

Inversion assuming weak scattering

DEFF Research Database (Denmark)

Xenaki, Angeliki; Gerstoft, Peter; Mosegaard, Klaus

2013-01-01

due to the complex nature of the field. A method based on linear inversion is employed to infer information about the statistical properties of the scattering field from the obtained cross-spectral matrix. A synthetic example based on an active high-frequency sonar demonstrates that the proposed...

Integrating the Toda Lattice with Self-Consistent Source via Inverse Scattering Method

International Nuclear Information System (INIS)

Urazboev, Gayrat

2012-01-01

In this work, there is shown that the solutions of Toda lattice with self-consistent source can be found by the inverse scattering method for the discrete Sturm-Liuville operator. For the considered problem the one-soliton solution is obtained.

Inverse scattering transform for the time dependent Schrödinger equation with applications to the KPI equation

Science.gov (United States)

Zhou, Xin

1990-03-01

For the direct-inverse scattering transform of the time dependent Schrödinger equation, rigorous results are obtained based on an opertor-triangular-factorization approach. By viewing the equation as a first order operator equation, similar results as for the first order n x n matrix system are obtained. The nonlocal Riemann-Hilbert problem for inverse scattering is shown to have solution.

Inverse scattering solution of the Chew-Low equation

International Nuclear Information System (INIS)

Nakano, K.

1985-01-01

Techniques for solving the inverse scattering problem are applied to the Chew-Low equation to obtain the nucleon form factor directly from the experimental phase shifts. A new dispersion relation is derived for the P 11 wave because of its sign-changing phase shift. A self-consistent solution for each channel is obtained, but the universality of form factor is not confirmed. Also, an iterative procedure based on Omnes' method is developed in order to solve coupled-channel, singular integral equations. (orig.)

Inverse problems of geophysics

International Nuclear Information System (INIS)

Yanovskaya, T.B.

2003-07-01

This report gives an overview and the mathematical formulation of geophysical inverse problems. General principles of statistical estimation are explained. The maximum likelihood and least square fit methods, the Backus-Gilbert method and general approaches for solving inverse problems are discussed. General formulations of linearized inverse problems, singular value decomposition and properties of pseudo-inverse solutions are given

Introduction to ground penetrating radar inverse scattering and data processing

CERN Document Server

Persico, Raffaele

2014-01-01

This book presents a comprehensive treatment of ground penetrating radar using both forward and inverse scattering mathematical techniques. Use of field data instead of laboratory data enables readers to envision real-life underground imaging; a full color insert further clarifies understanding. Along with considering the practical problem of achieving interpretable underground images, this book also features significant coverage of the problem's mathematical background. This twofold approach provides a resource that will appeal both to application oriented geologists and testing specialists,

Scattering of particles with inclusions. Modeling and inverse problem solution in the Rayleigh-Gans approximation

International Nuclear Information System (INIS)

Otero, F A; Frontini, G L; Elicabe, G E

2011-01-01

An analytic model for the scattering of a spherical particle with spherical inclusions has been proposed under the RG approximation. The model can be used without limitations to describe an X-ray scattering experiment. However, for light scattering several conditions must be fulfilled. Based on this model an inverse methodology is proposed to estimate the radii of host particle and inclusions, the number of inclusions and the Distance Distribution Functions (DDF's) of the distances between inclusions and the distances between inclusions and the origin of coordinates. The methodology is numerically tested in a light scattering example in which the host particle is eliminated by matching the refractive indices of host particle and medium. The results obtained for this cluster particle are very satisfactory.

The possibilities of linearized inversion of internally scattered seismic data

KAUST Repository

Aldawood, Ali

2014-08-05

Least-square migration is an iterative linearized inversion scheme that tends to suppress the migration artifacts and enhance the spatial resolution of the migrated image. However, standard least-square migration, based on imaging single scattering energy, may not be able to enhance events that are mainly illuminated by internal multiples such as vertical and nearly vertical faults. To alleviate this problem, we propose a linearized inversion framework to migrate internally multiply scattered energy. We applied this least-square migration of internal multiples to image a vertical fault. Tests on synthetic data demonstrate the ability of the proposed method to resolve a vertical fault plane that is poorly resolved by least-square imaging using primaries only. We, also, demonstrate the robustness of the proposed scheme in the presence of white Gaussian random observational noise and in the case of imaging the fault plane using inaccurate migration velocities.

The possibilities of linearized inversion of internally scattered seismic data

KAUST Repository

Aldawood, Ali; Alkhalifah, Tariq Ali; Hoteit, Ibrahim; Zuberi, Mohammad; Turkiyyah, George

2014-01-01

Least-square migration is an iterative linearized inversion scheme that tends to suppress the migration artifacts and enhance the spatial resolution of the migrated image. However, standard least-square migration, based on imaging single scattering energy, may not be able to enhance events that are mainly illuminated by internal multiples such as vertical and nearly vertical faults. To alleviate this problem, we propose a linearized inversion framework to migrate internally multiply scattered energy. We applied this least-square migration of internal multiples to image a vertical fault. Tests on synthetic data demonstrate the ability of the proposed method to resolve a vertical fault plane that is poorly resolved by least-square imaging using primaries only. We, also, demonstrate the robustness of the proposed scheme in the presence of white Gaussian random observational noise and in the case of imaging the fault plane using inaccurate migration velocities.

Multi-scattering inversion for low model wavenumbers

KAUST Repository

Alkhalifah, Tariq Ali; Wu, Zedong

2015-01-01

modeled from the source and those corresponding to single and double scattering to update both the velocity model and the component of the velocity (perturbation) responsible for the single and double scattering. The combined inversion helps us access most

Support Minimized Inversion of Acoustic and Elastic Wave Scattering

Science.gov (United States)

Safaeinili, Ali

Inversion of limited data is common in many areas of NDE such as X-ray Computed Tomography (CT), Ultrasonic and eddy current flaw characterization and imaging. In many applications, it is common to have a bias toward a solution with minimum (L^2)^2 norm without any physical justification. When it is a priori known that objects are compact as, say, with cracks and voids, by choosing "Minimum Support" functional instead of the minimum (L^2)^2 norm, an image can be obtained that is equally in agreement with the available data, while it is more consistent with what is most probably seen in the real world. We have utilized a minimum support functional to find a solution with the smallest volume. This inversion algorithm is most successful in reconstructing objects that are compact like voids and cracks. To verify this idea, we first performed a variational nonlinear inversion of acoustic backscatter data using minimum support objective function. A full nonlinear forward model was used to accurately study the effectiveness of the minimized support inversion without error due to the linear (Born) approximation. After successful inversions using a full nonlinear forward model, a linearized acoustic inversion was developed to increase speed and efficiency in imaging process. The results indicate that by using minimum support functional, we can accurately size and characterize voids and/or cracks which otherwise might be uncharacterizable. An extremely important feature of support minimized inversion is its ability to compensate for unknown absolute phase (zero-of-time). Zero-of-time ambiguity is a serious problem in the inversion of the pulse-echo data. The minimum support inversion was successfully used for the inversion of acoustic backscatter data due to compact scatterers without the knowledge of the zero-of-time. The main drawback to this type of inversion is its computer intensiveness. In order to make this type of constrained inversion available for common use, work

Effective theories of scattering with an attractive inverse-square potential and the three-body problem

International Nuclear Information System (INIS)

Barford, Thomas; Birse, Michael C

2005-01-01

A distorted-wave version of the renormalization group is applied to scattering by an inverse-square potential and to three-body systems. In attractive three-body systems, the short-distance wavefunction satisfies a Schroedinger equation with an attractive inverse-square potential, as shown by Efimov. The resulting oscillatory behaviour controls the renormalization of the three-body interactions, with the renormalization-group flow tending to a limit cycle as the cut-off is lowered. The approach used here leads to single-valued potentials with discontinuities as the bound states are cut off. The perturbations around the cycle start with a marginal term whose effect is simply to change the phase of the short-distance oscillations, or the self-adjoint extension of the singular Hamiltonian. The full power counting in terms of the energy and two-body scattering length is constructed for short-range three-body forces

The seismic reflection inverse problem

International Nuclear Information System (INIS)

Symes, W W

2009-01-01

The seismic reflection method seeks to extract maps of the Earth's sedimentary crust from transient near-surface recording of echoes, stimulated by explosions or other controlled sound sources positioned near the surface. Reasonably accurate models of seismic energy propagation take the form of hyperbolic systems of partial differential equations, in which the coefficients represent the spatial distribution of various mechanical characteristics of rock (density, stiffness, etc). Thus the fundamental problem of reflection seismology is an inverse problem in partial differential equations: to find the coefficients (or at least some of their properties) of a linear hyperbolic system, given the values of a family of solutions in some part of their domains. The exploration geophysics community has developed various methods for estimating the Earth's structure from seismic data and is also well aware of the inverse point of view. This article reviews mathematical developments in this subject over the last 25 years, to show how the mathematics has both illuminated innovations of practitioners and led to new directions in practice. Two themes naturally emerge: the importance of single scattering dominance and compensation for spectral incompleteness by spatial redundancy. (topical review)

Inverse scattering theory foundations of tomography with diffracting wavefields

International Nuclear Information System (INIS)

Devaney, A.J.

1987-01-01

The underlying mathematical models employed in reflection and transmission computed tomography using diffracting wavefields (called diffraction tomography) are reviewed and shown to have a rigorous basis in inverse scattering theory. In transmission diffraction tomography the underlying wave model is shown to be the Rytov approximation to the complex phase of the wavefield transmitted by the object being probed while in reflection diffraction tomography the underlying wave model is shown to be the Born approximation to the backscattered wavefield from the object. In both cases the goal of the reconstruction process is the determination of the objects's complex index of refraction as a function of position r/sup →/ and, possibly, the frequency ω of the probing wavefield. By use of these approximations the reconstruction problem for both transmission and reflection diffraction tomography can be cast into the simple and elegant form of linearized inverse scattering theory. Linearized inverse scattering theory is shown to lead directly to generalized projection-slice theorems for both reflection and transmission diffraction tomography that provide a simple mathematical relationship between the object's complex index of refraction (the unknown) and the data (the complex phase of the transmitted wave or the complex amplitude of the reflected wave). The conventional projection-slice theorem of X-ray CT is shown to result from the generalized projection-slice theorem for transmission diffraction tomography in the limit of vanishing wavelength (in the absence of wave effects). Fourier based and back-projection type reconstruction algorithms are shown to be directly derivable from the generalized projection-slice theorems

Continuity of the direct and inverse problems in one-dimensional scattering theory and numerical solution of the inverse problem

International Nuclear Information System (INIS)

Moura, C.A. de.

1976-09-01

We propose an algorithm for computing the potential V(x) associated to the one-dimensional Schroedinger operator E identical to - d 2 /dx 2 + V(x) -infinite < x< infinite from knowledge of the S.matrix, more exactly, of one of the reelection coefficients. The convergence of the algorithm is guaranteed by the stability results obtained for both the direct and inverse problems

Inverse Scattering Method and Soliton Solution Family for String Effective Action

International Nuclear Information System (INIS)

Ya-Jun, Gao

2009-01-01

A modified Hauser–Ernst-type linear system is established and used to develop an inverse scattering method for solving the motion equations of the string effective action describing the coupled gravity, dilaton and Kalb–Ramond fields. The reduction procedures in this inverse scattering method are found to be fairly simple, which makes the proposed inverse scattering method applied fine and effective. As an application, a concrete family of soliton solutions for the considered theory is obtained

A Literature Survey on Inverse Scattering for Electron Density Profile Determination. Volume II.

Science.gov (United States)

1981-09-24

THE INVERSE SCATTERING PROBLEM4 FOR THE EQUAT ION Of ACOUSTIC$ AVILA, G.S.S. DEPT. DE MATEMATICA . INST. DE CIENCIAS EXATAS. UNIV. Of BRASILIA...of Colict support Portinari. Joao C. Departamento do Matematica . Pontificia Universidade Catolica do Rio de Janeiro, Rio do Janeiro. Brasil J. Math

Extended resolvent and inverse scattering with an application to KPI

International Nuclear Information System (INIS)

Boiti, M.; Pempinelli, F.; Pogrebkov, A.K.; Prinari, B.

2003-01-01

We present in detail an extended resolvent approach for investigating linear problems associated to 2+1 dimensional integrable equations. Our presentation is based as an example on the nonstationary Schroedinger equation with potential being a perturbation of the one-soliton potential by means of a decaying two-dimensional function. Modification of the inverse scattering theory as well as properties of the Jost solutions and spectral data as follows from the resolvent approach are given

Full waveform inversion based on scattering angle enrichment with application to real dataset

KAUST Repository

Wu, Zedong

2015-08-19

Reflected waveform inversion (RWI) provides a method to reduce the nonlinearity of the standard full waveform inversion (FWI). However, the drawback of the existing RWI methods is inability to utilize diving waves and the extra sensitivity to the migrated image. We propose a combined FWI and RWI optimization problem through dividing the velocity into the background and perturbed components. We optimize both the background and perturbed components, as independent parameters. The new objective function is quadratic with respect to the perturbed component, which will reduce the nonlinearity of the optimization problem. Solving this optimization provides a true amplitude image and utilizes the diving waves to update the velocity of the shallow parts. To insure a proper wavenumber continuation, we use an efficient scattering angle filter to direct the inversion at the early stages to direct energy corresponding to large (smooth velocity) scattering angles to the background velocity update and the small (high wavenumber) scattering angles to the perturbed velocity update. This efficient implementation of the filter is fast and requires less memory than the conventional approach based on extended images. Thus, the new FWI procedure updates the background velocity mainly along the wavepath for both diving and reflected waves in the initial stages. At the same time, it updates the perturbation with mainly reflections (filtering out the diving waves). To demonstrate the capability of this method, we apply it to a real 2D marine dataset.

The shifting zoom: new possibilities for inverse scattering on electrically large domains

Science.gov (United States)

Persico, Raffaele; Ludeno, Giovanni; Soldovieri, Francesco; De Coster, Alberic; Lambot, Sebastien

2017-04-01

Inverse scattering is a subject of great interest in diagnostic problems, which are in their turn of interest for many applicative problems as investigation of cultural heritage, characterization of foundations or subservices, identification of unexploded ordnances and so on [1-4]. In particular, GPR data are usually focused by means of migration algorithms, essentially based on a linear approximation of the scattering phenomenon. Migration algorithms are popular because they are computationally efficient and do not require the inversion of a matrix, neither the calculation of the elements of a matrix. In fact, they are essentially based on the adjoint of the linearised scattering operator, which allows in the end to write the inversion formula as a suitably weighted integral of the data [5]. In particular, this makes a migration algorithm more suitable than a linear microwave tomography inversion algorithm for the reconstruction of an electrically large investigation domain. However, this computational challenge can be overcome by making use of investigation domains joined side by side, as proposed e.g. in ref. [3]. This allows to apply a microwave tomography algorithm even to large investigation domains. However, the joining side by side of sequential investigation domains introduces a problem of limited (and asymmetric) maximum view angle with regard to the targets occurring close to the edges between two adjacent domains, or possibly crossing these edges. The shifting zoom is a method that allows to overcome this difficulty by means of overlapped investigation and observation domains [6-7]. It requires more sequential inversion with respect to adjacent investigation domains, but the really required extra-time is minimal because the matrix to be inverted is calculated ones and for all, as well as its singular value decomposition: what is repeated more time is only a fast matrix-vector multiplication. References [1] M. Pieraccini, L. Noferini, D. Mecatti, C

Solution of inverse localization problem associated to multistatic radar system

Directory of Open Access Journals (Sweden)

Boutkhil M.

2016-01-01

Full Text Available This work deals with the problem of inverse localization by a target with the aim to retrieve the position of the target, given the intensity and phase of the electromagnetic waves scattered by this object. Assuming the surface cross section to be known as well as the intensity and phase of the scattered waves, the target position was reconstructed through the echo signals scattered of each bistatic. We develop in the same time a multistatic ambiguity function trough bistatic ambiguity function to investigate several fundamental aspects that determine multistatic radar performance. We used a multistatic radar constructed of two bistatic radars, two transmitters and one receiver.

«Paralipomena» on uniqueness in inverse scattering from a finite number of data

Directory of Open Access Journals (Sweden)

R. Persico

2007-06-01

Full Text Available This paper shows new proof of non-uniqueness of the solution for the retrieving of a compact-supported function starting from a finite number of samples of its spectrum. As will be shown, this is relevant for linear inverse scattering problems, that in many cases can be recast as the reconstruction of a compact supported function from a finite set of samples of its spectrum. Since this reconstruction is not unique, from a practical point of view, any linear inverse scattering algorithm that can be recast in terms of a Fourier relationship between unknowns and data necessarily «trusts» on the absence of invisible objects in the particular situation at hand.

A systematic approach to robust preconditioning for gradient-based inverse scattering algorithms

International Nuclear Information System (INIS)

Nordebo, Sven; Fhager, Andreas; Persson, Mikael; Gustafsson, Mats

2008-01-01

This paper presents a systematic approach to robust preconditioning for gradient-based nonlinear inverse scattering algorithms. In particular, one- and two-dimensional inverse problems are considered where the permittivity and conductivity profiles are unknown and the input data consist of the scattered field over a certain bandwidth. A time-domain least-squares formulation is employed and the inversion algorithm is based on a conjugate gradient or quasi-Newton algorithm together with an FDTD-electromagnetic solver. A Fisher information analysis is used to estimate the Hessian of the error functional. A robust preconditioner is then obtained by incorporating a parameter scaling such that the scaled Fisher information has a unit diagonal. By improving the conditioning of the Hessian, the convergence rate of the conjugate gradient or quasi-Newton methods are improved. The preconditioner is robust in the sense that the scaling, i.e. the diagonal Fisher information, is virtually invariant to the numerical resolution and the discretization model that is employed. Numerical examples of image reconstruction are included to illustrate the efficiency of the proposed technique

Off-shell T-matrices from inverse scattering

International Nuclear Information System (INIS)

Von Geramb, H.V.; Amos, K.A.

1989-01-01

Inverse scattering theory is used to determine local, energy independent, coordinate space nucleon-nucleon potentials. Inversions are made of phase shifts obtained by analyzes of data and from meson exchange theory, in particular the Paris and the Bonn parametrizations. Half off-shell T-matrices are generated to compare the exact meson theoretical results with those of inversion and it is found that phase equivalent interactions have essentially the same off-shell behaviour for any physically significant range of momenta. 8 refs., 8 figs

Extended resolvent and inverse scattering with an application to KPI

Science.gov (United States)

Boiti, M.; Pempinelli, F.; Pogrebkov, A. K.; Prinari, B.

2003-08-01

We present in detail an extended resolvent approach for investigating linear problems associated to 2+1 dimensional integrable equations. Our presentation is based as an example on the nonstationary Schrödinger equation with potential being a perturbation of the one-soliton potential by means of a decaying two-dimensional function. Modification of the inverse scattering theory as well as properties of the Jost solutions and spectral data as follows from the resolvent approach are given.

Constraint on Parameters of Inverse Compton Scattering Model for ...

Indian Academy of Sciences (India)

B2319+60, two parameters of inverse Compton scattering model, the initial Lorentz factor and the factor of energy loss of relativistic particles are constrained. Key words. Pulsar—inverse Compton scattering—emission mechanism. 1. Introduction. Among various kinds of models for pulsar radio emission, the inverse ...

Inverse scattering transform and soliton solutions for square matrix nonlinear Schrödinger equations with non-zero boundary conditions

Science.gov (United States)

Prinari, Barbara; Demontis, Francesco; Li, Sitai; Horikis, Theodoros P.

2018-04-01

The inverse scattering transform (IST) with non-zero boundary conditions at infinity is developed for an m × m matrix nonlinear Schrödinger-type equation which, in the case m = 2, has been proposed as a model to describe hyperfine spin F = 1 spinor Bose-Einstein condensates with either repulsive interatomic interactions and anti-ferromagnetic spin-exchange interactions (self-defocusing case), or attractive interatomic interactions and ferromagnetic spin-exchange interactions (self-focusing case). The IST for this system was first presented by Ieda et al. (2007) , using a different approach. In our formulation, both the direct and the inverse problems are posed in terms of a suitable uniformization variable which allows to develop the IST on the standard complex plane, instead of a two-sheeted Riemann surface or the cut plane with discontinuities along the cuts. Analyticity of the scattering eigenfunctions and scattering data, symmetries, properties of the discrete spectrum, and asymptotics are derived. The inverse problem is posed as a Riemann-Hilbert problem for the eigenfunctions, and the reconstruction formula of the potential in terms of eigenfunctions and scattering data is provided. In addition, the general behavior of the soliton solutions is analyzed in detail in the 2 × 2 self-focusing case, including some special solutions not previously discussed in the literature.

Phaseless tomographic inverse scattering in Banach spaces

International Nuclear Information System (INIS)

Estatico, C.; Fedeli, A.; Pastorino, M.; Randazzo, A.; Tavanti, E.

2016-01-01

In conventional microwave imaging, a hidden dielectric object under test is illuminated by microwave incident waves and the field it scatters is measured in magnitude and phase in order to retrieve the dielectric properties by solving the related non-homogenous Helmholtz equation or its Lippmann-Schwinger integral formulation. Since the measurement of the phase of electromagnetic waves can be still considered expensive in real applications, in this paper only the magnitude of the scattering wave fields is measured in order to allow a reduction of the cost of the measurement apparatus. In this respect, we firstly analyse the properties of the phaseless scattering nonlinear forward modelling operator in its integral form and we provide an analytical expression for computing its Fréchet derivative. Then, we propose an inexact Newton method to solve the associated nonlinear inverse problems, where any linearized step is solved by a L p Banach space iterative regularization method which acts on the dual space L p* . Indeed, it is well known that regularization in special Banach spaces, such us L p with 1 < p < 2, allows to promote sparsity and to reduce Gibbs phenomena and over-smoothness. Preliminary results concerning numerically computed field data are shown. (paper)

Study on the generalized WKB approximation for the inverse scattering problem at fixed energy for complex potentials

International Nuclear Information System (INIS)

Pozdnyakov, Yu.A.; Terenetskij, K.O.

1981-01-01

The approximate method for solution of the inverse scattering problem (ISP) at fixed energy for complex spherically symmetric potentials decreasing faster 1/r is considered. The method is based on using a generalized WKB approximation. For the designed potential V(r) a sufficiently ''close'' reference potential V(r) has been chosen. For both potentials S-matrix elements (ME) have been calculated and inversion procedure has been carried out. S-ME have been calculated for integral-valued and intermediate angular moment values. S-ME are presented in a graphical form for being restored reference, and restored potentials for proton scattering with Esub(p)=49.48 MeV energy on 12 C nuclei. The restoration is the better the ''closer'' the sought-for potential to the reference one. This allows to specify the potential by means of iterations: the restored potential can be used as a reference one, etc. The operation of a restored potential smoothing before the following iteration is introduced. Drawbacks and advantages of the ISP solution method under consideration are pointed out. The method application is strongly limited by the requirement that the energy should be higher than a certain ''critical'' one. The method is applicable in a wider region of particle energies (in the low-energies direction) than the ordinary WKB method. The method is more simple in realization conformably to complex potentials. The investigations carried out of the proposed ISP solution method at fixed energy for complex spherically-symmetric potentials allow to conclude that the method can be successFully applied to specify the central part of interaction of nucleons, α-particles and heavy ions of average and high energies with atomic nuclei [ru

Solving the relativistic inverse scattering problem on the basis of n/d equations and application of the resulting solution to analysis of pion-nucleon interaction at low and intermediate energies

International Nuclear Information System (INIS)

Safronov, A.N.

2007-01-01

Full text: The pion-nucleon dynamics is one of the most fundamental problems in nuclear and particle physics. It is now widely believed that QCD is fundamental theory of strong interactions. On this basis all hadron-hadron interactions are completely determined by the underlying quark-gluon dynamics. However, due to the formidable mathematical problems raised by the non-perturbative character of QCD at low and intermediate energies, we are still far from a quantitative understanding hadron-hadron interactions from this point of view. Recently the relativistic approaches to constructing effective interaction operators between strongly interacting composite particles has been proposed on the basis of analytic S-matrix theory and methods for solving the inverse quantum scattering problem. The kernel of Marchenko equation in theory of inverse scattering problem can be expressed in terms of the discontinuity of the partial wave amplitude on dynamic cut in the complex s=k 2 plane, k being the relative momentum of colliding particles. The discontinuities of partial-wave amplitudes are determined by model-independent quantities (renormalized vertex constants and amplitudes of sub-processes involving on-mass-shell particles off physical region) and can be calculated by methods of relativistic quantum field theory within various dynamical approaches. In particular, effective field theory can be used to calculate the discontinuities across dynamical cuts closest to physical region. In present work a new manifestly Poincare-invariant approach to solving the inverse scattering problem is developed with allowance for inelasticity effects. The equations of the N/D method are used as dynamical equations in this approach. With the help of N/D-equations it was earlier shown that solution of a scattering problem in case of nonzero angular momentum does not exist for arbitrary discontinuity of partial-wave amplitude. The method is elaborated allowing to determine contributions of

Invisibility problem in acoustics, electromagnetism and heat transfer. Inverse design method

Science.gov (United States)

Alekseev, G.; Tokhtina, A.; Soboleva, O.

2017-10-01

Two approaches (direct design and inverse design methods) for solving problems of designing devices providing invisibility of material bodies of detection using different physical fields - electromagnetic, acoustic and static are discussed. The second method is applied for solving problems of designing cloaking devices for the 3D stationary thermal scattering model. Based on this method the design problems under study are reduced to respective control problems. The material parameters (radial and tangential heat conductivities) of the inhomogeneous anisotropic medium filling the thermal cloak and the density of auxiliary heat sources play the role of controls. A unique solvability of direct thermal scattering problem in the Sobolev space is proved and the new estimates of solutions are established. Using these results, the solvability of control problem is proved and the optimality system is derived. Based on analysis of optimality system, the stability estimates of optimal solutions are established and numerical algorithms for solving particular thermal cloaking problem are proposed.

Rayleigh scattering and nonlinear inversion of elastic waves

Energy Technology Data Exchange (ETDEWEB)

Gritto, Roland [Univ. of California, Berkeley, CA (United States)

1995-12-01

Rayleigh scattering of elastic waves by an inclusion is investigated and the limitations determined. In the near field of the inhomogeneity, the scattered waves are up to a factor of 300 stronger than in the far field, excluding the application of the far field Rayleigh approximation for this range. The investigation of the relative error as a function of parameter perturbation shows a range of applicability broader than previously assumed, with errors of 37% and 17% for perturbations of -100% and +100%, respectively. The validity range for the Rayleigh limit is controlled by large inequalities, and therefore, the exact limit is determined as a function of various parameter configurations, resulting in surprisingly high values of up to k_pR = 0.9. The nonlinear scattering problem can be solved by inverting for equivalent source terms (moments) of the scatterer, before the elastic parameters are determined. The nonlinear dependence between the moments and the elastic parameters reveals a strong asymmetry around the origin, which will produce different results for weak scattering approximations depending on the sign of the anomaly. Numerical modeling of cross hole situations shows that near field terms are important to yield correct estimates of the inhomogeneities in the vicinity of the receivers, while a few well positioned sources and receivers considerably increase the angular coverage, and thus the model resolution of the inversion parameters. The pattern of scattered energy by an inhomogeneity is complicated and varies depending on the object, the wavelength of the incident wave, and the elastic parameters involved. Therefore, it is necessary to investigate the direction of scattered amplitudes to determine the best survey geometry.

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

Science.gov (United States)

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.

Statistical perspectives on inverse problems

DEFF Research Database (Denmark)

Andersen, Kim Emil

of the interior of an object from electrical boundary measurements. One part of this thesis concerns statistical approaches for solving, possibly non-linear, inverse problems. Thus inverse problems are recasted in a form suitable for statistical inference. In particular, a Bayesian approach for regularisation...... problem is given in terms of probability distributions. Posterior inference is obtained by Markov chain Monte Carlo methods and new, powerful simulation techniques based on e.g. coupled Markov chains and simulated tempering is developed to improve the computational efficiency of the overall simulation......Inverse problems arise in many scientific disciplines and pertain to situations where inference is to be made about a particular phenomenon from indirect measurements. A typical example, arising in diffusion tomography, is the inverse boundary value problem for non-invasive reconstruction...

Inverse logarithmic potential problem

CERN Document Server

Cherednichenko, V G

1996-01-01

The Inverse and Ill-Posed Problems Series is a series of monographs publishing postgraduate level information on inverse and ill-posed problems for an international readership of professional scientists and researchers. The series aims to publish works which involve both theory and applications in, e.g., physics, medicine, geophysics, acoustics, electrodynamics, tomography, and ecology.

Scattering-angle based filtering of the waveform inversion gradients

KAUST Repository

Alkhalifah, Tariq Ali

2014-01-01

Full waveform inversion (FWI) requires a hierarchical approach to maneuver the complex non-linearity associated with the problem of velocity update. In anisotropic media, the non-linearity becomes far more complex with the potential trade-off between the multiparameter description of the model. A gradient filter helps us in accessing the parts of the gradient that are suitable to combat the potential non-linearity and parameter trade-off. The filter is based on representing the gradient in the time-lag normalized domain, in which the low scattering angle of the gradient update is initially muted out in the FWI implementation, in what we may refer to as a scattering angle continuation process. The result is a low wavelength update dominated by the transmission part of the update gradient. In this case, even 10 Hz data can produce vertically near-zero wavenumber updates suitable for a background correction of the model. Relaxing the filtering at a later stage in the FWI implementation allows for smaller scattering angles to contribute higher-resolution information to the model. The benefits of the extended domain based filtering of the gradient is not only it's ability in providing low wavenumber gradients guided by the scattering angle, but also in its potential to provide gradients free of unphysical energy that may correspond to unrealistic scattering angles.

Scattering-angle based filtering of the waveform inversion gradients

KAUST Repository

Alkhalifah, Tariq Ali

2014-11-22

Full waveform inversion (FWI) requires a hierarchical approach to maneuver the complex non-linearity associated with the problem of velocity update. In anisotropic media, the non-linearity becomes far more complex with the potential trade-off between the multiparameter description of the model. A gradient filter helps us in accessing the parts of the gradient that are suitable to combat the potential non-linearity and parameter trade-off. The filter is based on representing the gradient in the time-lag normalized domain, in which the low scattering angle of the gradient update is initially muted out in the FWI implementation, in what we may refer to as a scattering angle continuation process. The result is a low wavelength update dominated by the transmission part of the update gradient. In this case, even 10 Hz data can produce vertically near-zero wavenumber updates suitable for a background correction of the model. Relaxing the filtering at a later stage in the FWI implementation allows for smaller scattering angles to contribute higher-resolution information to the model. The benefits of the extended domain based filtering of the gradient is not only it\\'s ability in providing low wavenumber gradients guided by the scattering angle, but also in its potential to provide gradients free of unphysical energy that may correspond to unrealistic scattering angles.

Multi-parameter Analysis and Inversion for Anisotropic Media Using the Scattering Integral Method

KAUST Repository

Djebbi, Ramzi

2017-10-24

scattering integral method. Application to various synthetic and real data examples show accurate inversion results. I show that a good background ƞ model is required to accurately recover vh. For 3-D problems, I promote a hybrid approach, where efficient ray tracing is used to compute the sensitivity kernels. The proposed method highly reduces the computational cost.

Efficient scattering angle filtering for Full waveform inversion

KAUST Repository

Alkhalifah, Tariq Ali

2015-01-01

Controlling the scattering angles between the state and the adjoint variables for the energy admitted into an inversion gradient or an image can help improve these functions for objectives in full waveform inversion (FWI) or seismic imaging. However, the access of the scattering angle information usually requires an axis extension that could be costly, especially in 3D. For the purpose of a scattering angle filter, I develop techniques that utilize the mapping nature (no domain extension) of the filter for constant-velocity background models to interpolate between such filtered gradients using the actual velocity. The concept has well known roots in the application of phase-shift-plus-interpolation utilized commonly in the downward continuation process. If the difference between the minimum and maximum velocity of the background medium is large, we obtain filtered gradients corresponding to more constant velocity backgrounds and use linear interpolation between such velocities. The accuracy of this approximation for the Marmousi model gradient demonstrates the e ectiveness of the approach.

Efficient scattering angle filtering for Full waveform inversion

KAUST Repository

Alkhalifah, Tariq Ali

2015-08-19

Controlling the scattering angles between the state and the adjoint variables for the energy admitted into an inversion gradient or an image can help improve these functions for objectives in full waveform inversion (FWI) or seismic imaging. However, the access of the scattering angle information usually requires an axis extension that could be costly, especially in 3D. For the purpose of a scattering angle filter, I develop techniques that utilize the mapping nature (no domain extension) of the filter for constant-velocity background models to interpolate between such filtered gradients using the actual velocity. The concept has well known roots in the application of phase-shift-plus-interpolation utilized commonly in the downward continuation process. If the difference between the minimum and maximum velocity of the background medium is large, we obtain filtered gradients corresponding to more constant velocity backgrounds and use linear interpolation between such velocities. The accuracy of this approximation for the Marmousi model gradient demonstrates the e ectiveness of the approach.

On the Quantum Inverse problem for the continuous Heisenberg spin chain with axial anisotropy

International Nuclear Information System (INIS)

Roy Chowdhury, A.; Chanda, P.K.

1986-06-01

We have considered the Quantum Inverse problem for the continuous form of Heisenberg spin chain with anisotropy. The form of quantum R-matrix, the commutation rules for the scattering data, and the explicit structure of the excitation spectrum are obtained. (author)

Collision-induced Raman scattering and the peculiar case of neon: Anisotropic spectrum, anisotropy, and the inverse scattering problem

International Nuclear Information System (INIS)

Dixneuf, Sophie; Rachet, Florent; Chrysos, Michael

2015-01-01

Owing in part to the p orbitals of its filled L shell, neon has repeatedly come on stage for its peculiar properties. In the context of collision-induced Raman spectroscopy, in particular, we have shown, in a brief report published a few years ago [M. Chrysos et al., Phys. Rev. A 80, 054701 (2009)], that the room-temperature anisotropic Raman lineshape of Ne–Ne exhibits, in the far wing of the spectrum, a peculiar structure with an aspect other than a smooth wing (on a logarithmic plot) which contrasts with any of the existing studies, and whose explanation lies in the distinct way in which overlap and exchange interactions interfere with the classical electrostatic ones in making the polarizability anisotropy, α ∥ − α ⊥ . Here, we delve deeper into that study by reporting data for that spectrum up to 450 cm −1 and for even- and odd-order spectral moments up to M 6 , as well as quantum lineshapes, generated from SCF, CCSD, and CCSD(T) models for α ∥ − α ⊥ , which are critically compared with the experiment. On account of the knowledge of the spectrum over the augmented frequency domain, we show how the inverse scattering problem can be tackled both effectively and economically, and we report an analytic function for the anisotropy whose quantum lineshape faithfully reproduces our observations

Inverse problems in the Bayesian framework

International Nuclear Information System (INIS)

Calvetti, Daniela; Somersalo, Erkki; Kaipio, Jari P

2014-01-01

The history of Bayesian methods dates back to the original works of Reverend Thomas Bayes and Pierre-Simon Laplace: the former laid down some of the basic principles on inverse probability in his classic article ‘An essay towards solving a problem in the doctrine of chances’ that was read posthumously in the Royal Society in 1763. Laplace, on the other hand, in his ‘Memoirs on inverse probability’ of 1774 developed the idea of updating beliefs and wrote down the celebrated Bayes’ formula in the form we know today. Although not identified yet as a framework for investigating inverse problems, Laplace used the formalism very much in the spirit it is used today in the context of inverse problems, e.g., in his study of the distribution of comets. With the evolution of computational tools, Bayesian methods have become increasingly popular in all fields of human knowledge in which conclusions need to be drawn based on incomplete and noisy data. Needless to say, inverse problems, almost by definition, fall into this category. Systematic work for developing a Bayesian inverse problem framework can arguably be traced back to the 1980s, (the original first edition being published by Elsevier in 1987), although articles on Bayesian methodology applied to inverse problems, in particular in geophysics, had appeared much earlier. Today, as testified by the articles in this special issue, the Bayesian methodology as a framework for considering inverse problems has gained a lot of popularity, and it has integrated very successfully with many traditional inverse problems ideas and techniques, providing novel ways to interpret and implement traditional procedures in numerical analysis, computational statistics, signal analysis and data assimilation. The range of applications where the Bayesian framework has been fundamental goes from geophysics, engineering and imaging to astronomy, life sciences and economy, and continues to grow. There is no question that Bayesian

A multifrequency MUSIC algorithm for locating small inhomogeneities in inverse scattering

International Nuclear Information System (INIS)

Griesmaier, Roland; Schmiedecke, Christian

2017-01-01

We consider an inverse scattering problem for time-harmonic acoustic or electromagnetic waves with sparse multifrequency far field data-sets. The goal is to localize several small penetrable objects embedded inside an otherwise homogeneous background medium from observations of far fields of scattered waves corresponding to incident plane waves with one fixed incident direction but several different frequencies. We assume that the far field is measured at a few observation directions only. Taking advantage of the smallness of the scatterers with respect to wavelength we utilize an asymptotic representation formula for the far field to design and analyze a MUSIC-type reconstruction method for this setup. We establish lower bounds on the number of frequencies and receiver directions that are required to recover the number and the positions of an ensemble of scatterers from the given measurements. Furthermore we briefly sketch a possible application of the reconstruction method to the practically relevant case of multifrequency backscattering data. Numerical examples are presented to document the potentials and limitations of this approach. (paper)

An inverse-scattering approach to the physics of transition metals ...

African Journals Online (AJOL)

A method is developed for the deduction of a transition metal ion potential from a knowledge of the phase-shift. The method used is based the distorted plane – wave scattering approximation for the deduction of non singular potentials from scattering phase shifts in an inverse scattering approach. The resulting electron ...

Imaging electrical conductivity, permeability, and/or permittivity contrasts using the Born Scattering Inversion (BSI)

Science.gov (United States)

Darrh, A.; Downs, C. M.; Poppeliers, C.

2017-12-01

Born Scattering Inversion (BSI) of electromagnetic (EM) data is a geophysical imaging methodology for mapping weak conductivity, permeability, and/or permittivity contrasts in the subsurface. The high computational cost of full waveform inversion is reduced by adopting the First Born Approximation for scattered EM fields. This linearizes the inverse problem in terms of Born scattering amplitudes for a set of effective EM body sources within a 3D imaging volume. Estimation of scatterer amplitudes is subsequently achieved by solving the normal equations. Our present BSI numerical experiments entail Fourier transforming real-valued synthetic EM data to the frequency-domain, and minimizing the L2 residual between complex-valued observed and predicted data. We are testing the ability of BSI to resolve simple scattering models. For our initial experiments, synthetic data are acquired by three-component (3C) electric field receivers distributed on a plane above a single point electric dipole within a homogeneous and isotropic wholespace. To suppress artifacts, candidate Born scatterer locations are confined to a volume beneath the receiver array. Also, we explore two different numerical linear algebra algorithms for solving the normal equations: Damped Least Squares (DLS), and Non-Negative Least Squares (NNLS). Results from NNLS accurately recover the source location only for a large dense 3C receiver array, but fail when the array is decimated, or is restricted to horizontal component data. Using all receiver stations and all components per station, NNLS results are relatively insensitive to a sub-sampled frequency spectrum, suggesting that coarse frequency-domain sampling may be adequate for good target resolution. Results from DLS are insensitive to diminishing array density, but contain spatially oscillatory structure. DLS-generated images are consistently centered at the known point source location, despite an abundance of surrounding structure.

Rational reflection coefficient and inverse scattering on the line

International Nuclear Information System (INIS)

Sabatier, P.C.

1983-01-01

Inverse scattering for the Schroedinger equation on the line is studied for reflection and transmission coefficients that satisfy the usual regularity conditions and are rational functions of k. The origin is still a particular point, but the potentials do not need to be cut at this point like in previous studies. Giving up this restriction corresponds to the existence of poles for both reflection coefficients in both upper and lower half k-planes. It is shown that the problem reduces to solving a linear algebraic system. A different algorithm, made of a sequence of Darboux-Backlund transforms, gives also the solution in closed form and enables to study separately modifications of both sides of the potential due to the introduction of poles. Thus it paves the way for approximation studies. Generalizations and particular problems will be studied in forthcoming papers

Direct Problems and Inverse Problems in Biometric Systems

OpenAIRE

Mihailescu Marius Iulian

2013-01-01

The article purpose is to describe the two sides of biometrics technologies, direct problems and inverse problems. The advance that we face today in field of Information Technology makes Information Security an inseparable part. The authentication has a huge role when we deal about security. The problems that can appear in implementing and developing biometrics systems is raising many problems, and one of the goal of this article is to focus on direct and inverse problems which is a new and c...

Remarks on the inverse scattering transform associated with toda equations

Science.gov (United States)

Ablowitz, Mark J.; Villorroel, J.

The Inverse Scattering Transforms used to solve both the 2+1 Toda equation and a novel reduction, the Toda differential-delay equations are outlined. There are a number of interesting features associated with these systems and the related scattering theory.

Parameter estimation and inverse problems

CERN Document Server

Aster, Richard C; Thurber, Clifford H

2005-01-01

Parameter Estimation and Inverse Problems primarily serves as a textbook for advanced undergraduate and introductory graduate courses. Class notes have been developed and reside on the World Wide Web for faciliting use and feedback by teaching colleagues. The authors'' treatment promotes an understanding of fundamental and practical issus associated with parameter fitting and inverse problems including basic theory of inverse problems, statistical issues, computational issues, and an understanding of how to analyze the success and limitations of solutions to these probles. The text is also a practical resource for general students and professional researchers, where techniques and concepts can be readily picked up on a chapter-by-chapter basis.Parameter Estimation and Inverse Problems is structured around a course at New Mexico Tech and is designed to be accessible to typical graduate students in the physical sciences who may not have an extensive mathematical background. It is accompanied by a Web site that...

Heavy ion scattering; a fixed energy inverse problem

International Nuclear Information System (INIS)

Amos, K.

1993-01-01

Heavy ion scattering has been studied quite intensively in the last decade and central in most analyses of data from such experiments be they on fusion, particle transfer or internal state excitations of the colliding pair, is the inter-ion interaction affecting their relative motion. It is customary to use the elastic scattering data to constrain solutions of the (nonrelativistic) Schroedinger equation to ascertain the character of that (central and complex) heavy ion potential. These matters for projectiles ranging from the lightest 'heavy' ion, a proton, to Oxygen nuclei are considered in brief herein. The targets range from 12 C to 208 Pb. The central entity in the analyses to be discussed will be the S-function, and so for completeness, the simple potential scattering theory details are presented that specify the S-function and relate it to measured cross-sections. 20 refs., 18 figs

The factorization method for inverse acoustic scattering in a layered medium

International Nuclear Information System (INIS)

Bondarenko, Oleksandr; Kirsch, Andreas; Liu, Xiaodong

2013-01-01

In this paper, we consider a problem of inverse acoustic scattering by an impenetrable obstacle embedded in a layered medium. We will show that the factorization method can be applied to recover the embedded obstacle; that is, the equation F-tilde g =φ z is solvable if and only if the sampling point z is in the interior of the unknown obstacle. Here, F-tilde is a self-adjoint operator related to the far field operator and ϕ z is the far field pattern of the Green function with respect to the problem of scattering by the background medium for point z. The validity of the factorization method is proven with the help of a mixed reciprocity principle and an application of the scattering operator. Due to the established mixed reciprocity principle, knowledge of the Green function for the background medium is no longer required, which makes the method attractive from the computational point of view. The paper is only concerned with sound-soft obstacles, but the analysis can be easily extended for sound-hard obstacles, or obstacles with separated sound-soft and sound-hard parts. Finally, we provide an explicit example for a radially symmetric case and present some numerical examples. (paper)

Some aspects of the inverse problem of scattering at fixed energy

International Nuclear Information System (INIS)

Coudray, C.

1979-01-01

The first two chapters deal with the Newton-Sabatier method. Numerical tests are performed for real and complex potentials. They allow the study of the respective influences of energy, and of the internal parameters of the potential: its shape, depth and range. Within certain limits, good agreements are obtained. In particular, it is shown that they always require energies larger than a 'critical' energy, the dependence of which in function of the internal parameters of the potential being analyzed. Then the third chapter is devoted to transparent and quasi-transparent potentials in Born approximation. A class of such potentials is exhibited and studied. All of them oscillate, and their decrease at infinity may be chosen according to any arbitrary power of the variable. One of them is the Born approximation of the transparent potential of the Newton-Sabatier method. The last chapter concerns finite range complex potentials belonging to a well-defined class. For such potentials, a set of coherent inverse scattering date is given. The corresponding fundamental equation is written and shown to possess an unique solution [fr

Optimization and inverse problems in electromagnetism

CERN Document Server

Wiak, Sławomir

2003-01-01

From 12 to 14 September 2002, the Academy of Humanities and Economics (AHE) hosted the workshop "Optimization and Inverse Problems in Electromagnetism". After this bi-annual event, a large number of papers were assembled and combined in this book. During the workshop recent developments and applications in optimization and inverse methodologies for electromagnetic fields were discussed. The contributions selected for the present volume cover a wide spectrum of inverse and optimal electromagnetic methodologies, ranging from theoretical to practical applications. A number of new optimal and inverse methodologies were proposed. There are contributions related to dedicated software. Optimization and Inverse Problems in Electromagnetism consists of three thematic chapters, covering: -General papers (survey of specific aspects of optimization and inverse problems in electromagnetism), -Methodologies, -Industrial Applications. The book can be useful to students of electrical and electronics engineering, computer sci...

Connection between Dirac and matrix Schroedinger inverse-scattering transforms

International Nuclear Information System (INIS)

Jaulent, M.; Leon, J.J.P.

1978-01-01

The connection between two applications of the inverse scattering method for solving nonlinear equations is established. The inverse method associated with the massive Dirac system (D) : (iσ 3 d/dx - i q 3 σ 1 - q 1 σ 2 + mσ 2 )Y = epsilonY is rediscovered from the inverse method associated with the 2 x 2 matrix Schroedinger equation (S) : Ysub(xx) + (k 2 -Q)Y = 0. Here Q obeys a nonlinear constraint equivalent to a linear constraint on the reflection coefficient for (S). (author)

Inversion of the total cross sections for electron-molecule and electron-atom scattering

International Nuclear Information System (INIS)

Lun, D.R.; Amos, K.; Allen, L.J.

1994-01-01

Inverse scattering theory has been applied to construct the interaction potentials from total cross sections as a function of energy for electrons scattered off of atoms and molecules. The underlying potentials are assumed to be real and energy independent and are evaluated using the Eikonal approximation and with real phase shifts determined from the total cross sections. The inversion potentials have been determined using either a high energy limit approximation or by using a fixed energy inversion method at select energies. These procedures have been used to analyse e - - CH 4 , e - - SiH 4 , e - -Kr and e - -Xe scattering data in particular. 14 refs., 1 tabs., 3 figs

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

Science.gov (United States)

Vourc'h, Eric; Rodet, Thomas

2015-11-01

, convex optimization, free discontinuity problems, metamodels, proper orthogonal decomposition, reduced models for the inversion, non-linear inverse scattering, image reconstruction and restoration, and applications (bio-medical imaging, non-destructive evaluation...). NCMIP 2015 was a one-day workshop held in May 2015 which attracted around 70 attendees. Each of the submitted papers has been reviewed by two reviewers. There have been 15 accepted papers. In addition, three international speakers were invited to present a longer talk. The workshop was supported by Institut Farman (ENS Cachan, CNRS) and endorsed by the following French research networks: GDR ISIS, GDR MIA, GDR MOA and GDR Ondes. The program committee acknowledges the following research laboratories: CMLA, LMT, LURPA and SATIE.

Inverse problems for Maxwell's equations

CERN Document Server

Romanov, V G

1994-01-01

The Inverse and Ill-Posed Problems Series is a series of monographs publishing postgraduate level information on inverse and ill-posed problems for an international readership of professional scientists and researchers. The series aims to publish works which involve both theory and applications in, e.g., physics, medicine, geophysics, acoustics, electrodynamics, tomography, and ecology.

EDITORIAL: Inverse Problems in Engineering

Science.gov (United States)

West, Robert M.; Lesnic, Daniel

2007-01-01

Presented here are 11 noteworthy papers selected from the Fifth International Conference on Inverse Problems in Engineering: Theory and Practice held in Cambridge, UK during 11-15 July 2005. The papers have been peer-reviewed to the usual high standards of this journal and the contributions of reviewers are much appreciated. The conference featured a good balance of the fundamental mathematical concepts of inverse problems with a diverse range of important and interesting applications, which are represented here by the selected papers. Aspects of finite-element modelling and the performance of inverse algorithms are investigated by Autrique et al and Leduc et al. Statistical aspects are considered by Emery et al and Watzenig et al with regard to Bayesian parameter estimation and inversion using particle filters. Electrostatic applications are demonstrated by van Berkel and Lionheart and also Nakatani et al. Contributions to the applications of electrical techniques and specifically electrical tomographies are provided by Wakatsuki and Kagawa, Kim et al and Kortschak et al. Aspects of inversion in optical tomography are investigated by Wright et al and Douiri et al. The authors are representative of the worldwide interest in inverse problems relating to engineering applications and their efforts in producing these excellent papers will be appreciated by many readers of this journal.

Quantitative photoplethysmography: Lambert-Beer law or inverse function incorporating light scatter.

Science.gov (United States)

Cejnar, M; Kobler, H; Hunyor, S N

1993-03-01

Finger blood volume is commonly determined from measurement of infra-red (IR) light transmittance using the Lambert-Beer law of light absorption derived for use in non-scattering media, even when such transmission involves light scatter around the phalangeal bone. Simultaneous IR transmittance and finger volume were measured over the full dynamic range of vascular volumes in seven subjects and outcomes compared with data fitted according to the Lambert-Beer exponential function and an inverse function derived for light attenuation by scattering materials. Curves were fitted by the least-squares method and goodness of fit was compared using standard errors of estimate (SEE). The inverse function gave a better data fit in six of the subjects: mean SEE 1.9 (SD 0.7, range 0.7-2.8) and 4.6 (2.2, 2.0-8.0) respectively (p < 0.02, paired t-test). Thus, when relating IR transmittance to blood volume, as occurs in the finger during measurements of arterial compliance, an inverse function derived from a model of light attenuation by scattering media gives more accurate results than the traditional exponential fit.

BOOK REVIEW: Inverse Problems. Activities for Undergraduates

Science.gov (United States)

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

Soliton solutions of the mixed discrete modified Korteweg-de Vries hierarchy via the inverse scattering transform

International Nuclear Information System (INIS)

Li Qi; Duan Qiuyuan; Zhang Jianbing

2012-01-01

The mixed discrete modified Korteweg-de Vries (mKdV) hierarchy and the Lax pair are derived. The hierarchy related to the Ablowitz-Ladik spectral problem is reduced to the isospectral discrete mKdV hierarchy and to the non-isospectral discrete mKdV hierarchy. N-soliton solutions of the hierarchies are obtained through inverse scattering transform.

Inverse scattering and GPR data processing: an Introduction

Science.gov (United States)

Persico, Raffaele

2014-05-01

relative answers have been inserted at the end of many chapters. On the other hand, it seemed also well advised to stress the fact that, within a GPR prospecting, the main involved parameters (especially the propagation velocity of the electromagnetic waves in the soil) and the same useful datum, i.e. "the scattered field", in most cases have to be worked out from the same GPR data. In particular, usually we don't have the possibility to measure apart the "incident field" and then retrieve the scattered field by means of an immediate subtraction operation. Indeed, these aspects are an intrinsic part of the GPR data processing, and should not left out from a text on this topic. In the end, GPR data processing has its own specificities within the larger framework of the inverse scattering problems, and the book has tried to put into evidence this fact too. Finally, even if this text is not focused on electronics, it was important to account for the fact that there are two categories of GPR systems, namely those working in time domain and those working in frequency domain. This implies some consequences in terms of the parametric choices in order to gather and process correctly the data, which has been devoted some attention too. The main aim of the book is to resume and gather together things mostly already known, but usually spread within different texts and contexts, often dealt with different approaches and expressed, let say, with different languages. The book is mainly thought of for Ph.D. students, students of master courses and university students at their last year in geophysics end engineering, but it is accessible to any GPR user with some minimal basis (i.e. at university level) on electromagnetism. Some small research work has been performed too, as e.g. with regard to the calculation in closed form of the Hermitian images for stepped frequency systems, or with respect to the introduction of the effective maximum view angle, or in order to propose a new plane

New rigorous asymptotic theorems for inverse scattering amplitudes

International Nuclear Information System (INIS)

Lomsadze, Sh.Yu.; Lomsadze, Yu.M.

1984-01-01

The rigorous asymptotic theorems both of integral and local types obtained earlier and establishing logarithmic and in some cases even power correlations aetdeen the real and imaginary parts of scattering amplitudes Fsub(+-) are extended to the inverse amplitudes 1/Fsub(+-). One also succeeds in establishing power correlations of a new type between the real and imaginary parts, both for the amplitudes themselves and for the inverse ones. All the obtained assertions are convenient to be tested in high energy experiments when the amplitudes show asymptotic behaviour

Seismic inverse scattering in the downward continuation approach

NARCIS (Netherlands)

Stolk, C.C.; de Hoop, M.V.

Seismic data are commonly modeled by a linearization around a smooth background medium in combination with a high frequency approximation. The perturbation of the medium coefficient is assumed to contain the discontinuities. This leads to two inverse problems, first the linearized inverse problem

Inverse problems for the Boussinesq system

International Nuclear Information System (INIS)

Fan, Jishan; Jiang, Yu; Nakamura, Gen

2009-01-01

We obtain two results on inverse problems for a 2D Boussinesq system. One is that we prove the Lipschitz stability for the inverse source problem of identifying a time-independent external force in the system with observation data in an arbitrary sub-domain over a time interval of the velocity and the data of velocity and temperature at a fixed positive time t 0 > 0 over the whole spatial domain. The other one is that we prove a conditional stability estimate for an inverse problem of identifying the two initial conditions with a single observation on a sub-domain

Resolvent approach for two-dimensional scattering problems. Application to the nonstationary Schroedinger problem and the KPI equation

International Nuclear Information System (INIS)

Boiti, M.; Pempinelli, F.; Pogrebkov, A.K.; Polivanov, M.C.

1993-01-01

The resolvent operator of the linear problem is determined as the full Green function continued in the complex domain in two variables. An analog of the known Hilbert identity is derived. The authors demonstrate the role of this identity in the study of two-dimensional scattering. Considering the nonstationary Schroedinger equation as an example, it is shown that all types of solutions of the linear problem, as well as spectral data known in the literature, are given as specific values of this unique function - the resolvent function. A new form of the inverse problem is formulated. 7 refs

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

Science.gov (United States)

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

Simplified solutions of the Cox-Thompson inverse scattering method at fixed energy

International Nuclear Information System (INIS)

Palmai, Tamas; Apagyi, Barnabas; Horvath, Miklos

2008-01-01

Simplified solutions of the Cox-Thompson inverse quantum scattering method at fixed energy are derived if a finite number of partial waves with only even or odd angular momenta contribute to the scattering process. Based on new formulae various approximate methods are introduced which also prove applicable to the generic scattering events

The inverse problem for the one-dimensional Schroedinger equation with an energy-dependent potential. II

International Nuclear Information System (INIS)

Jaulent, M.; Jean, C.

1976-01-01

The one-dimensional Schroedinger equation y + ''+ ) 7k 2 -V + (k,x){y + =0, x belonging to R, was previously considered when the potential V + (k,x) depends on the energy k 2 in the following way: V + (k,x)=U(x)+2kQ(x), (U(x), Q(x)) belonging to a large class of pairs of real potentials admitting no bound state). The two systems of differential and integral equations then introduced are solved. Then, investigating the inverse scattering problem it is found that a necessary and sufficient condition for one of the functions S + (k) and Ssub(-1)sup(+)(k) to be the scattering matrix associated with a pair (U(x), Q(x)) is that S + (k) (or equivalently Ssub(-1)sup(+)(k) belongs to the class S introduced. This pair is the only one admitting this function as its scattering matrix. Investigating the inverse reflection problem, it is found that a necessary and sufficient condition for a function S 21 + (k) to be the reflection coefficient to the right associated with a pair (U(x), Q(x)) is that S 21 + (k) belongs to the class R introduced. This pair is the only one admitting this function as

Elastic scattering of protons on 8Li nucleus in inverse kinematics

International Nuclear Information System (INIS)

Zhusupov, M.A.; Ibraeva, E.T.; Sanfirova, A.B.; Imambekov, O.

2002-01-01

In the present paper the proton elastic scattering on 8 Li in inverse kinematics is studies. The inverse kinematics means that a beam of radioactive nuclei is scattered on a stable hydrogen target. Proton as a target has an advantage during the interaction since it is stable and mechanism of proton-nucleus scattering is quite simple. 8 Li nucleus is considered in the three-body αtn-model with realistic potential of inter-cluster interactions. The wave function of this nucleus is calculated in the work where it was shown that such model well describes the main spectroscopic characteristics of the nucleus, root-mean square radius, binding energy, location of low laying energy levels, magnetic momentum and also total cross section and 7 Li(n, γ) 8 Li reaction rate at a wide energy region. Within Glauber-Sitenko multiply scattering theory, the differential cross section of elastic p 8 Li-scattering has been calculated. The first and the second multiplicities of scattering on nucleons and clusters of the nucleus were taken into account in Ω multiply scattering operator. There were considered several cases when as the initial parameters both amplitudes of nucleon-nucleon and nucleon-cluster scattering were taken. Sensitivity of the differential cross section both to the different wave functions of the target-nucleus and to the parameters of the elementary amplitudes and sensitivity to the scattering multiplicities at several beam energies has been investigated. Comparison with differential cross sections of elastic p 6 Li- and p 7 Li scattering has been carried out

Numerical solution of a coefficient inverse problem with multi-frequency experimental raw data by a globally convergent algorithm

Science.gov (United States)

Nguyen, Dinh-Liem; Klibanov, Michael V.; Nguyen, Loc H.; Kolesov, Aleksandr E.; Fiddy, Michael A.; Liu, Hui

2017-09-01

We analyze in this paper the performance of a newly developed globally convergent numerical method for a coefficient inverse problem for the case of multi-frequency experimental backscatter data associated to a single incident wave. These data were collected using a microwave scattering facility at the University of North Carolina at Charlotte. The challenges for the inverse problem under the consideration are not only from its high nonlinearity and severe ill-posedness but also from the facts that the amount of the measured data is minimal and that these raw data are contaminated by a significant amount of noise, due to a non-ideal experimental setup. This setup is motivated by our target application in detecting and identifying explosives. We show in this paper how the raw data can be preprocessed and successfully inverted using our inversion method. More precisely, we are able to reconstruct the dielectric constants and the locations of the scattering objects with a good accuracy, without using any advanced a priori knowledge of their physical and geometrical properties.

Microlocal analysis of a seismic linearized inverse problem

NARCIS (Netherlands)

Stolk, C.C.

1999-01-01

The seismic inverse problem is to determine the wavespeed c x in the interior of a medium from measurements at the boundary In this paper we analyze the linearized inverse problem in general acoustic media The problem is to nd a left inverse of the linearized forward map F or equivalently to nd the

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

Science.gov (United States)

2014-10-01

workshop were: algorithms and computational aspects of inversion, Bayesian estimation, Kernel methods, learning methods, convex optimization, free discontinuity problems, metamodels, proper orthogonal decomposition, reduced models for the inversion, non-linear inverse scattering, image reconstruction and restoration, and applications (bio-medical imaging, non-destructive evaluation...). NCMIP 2014 was a one-day workshop held in May 2014 which attracted around sixty attendees. Each of the submitted papers has been reviewed by two reviewers. There have been nine accepted papers. In addition, three international speakers were invited to present a longer talk. The workshop was supported by Institut Farman (ENS Cachan, CNRS) and endorsed by the following French research networks (GDR ISIS, GDR MIA, GDR MOA, GDR Ondes). The program committee acknowledges the following research laboratories: CMLA, LMT, LURPA, SATIE. Eric Vourc'h and Thomas Rodet

Introduction to inverse problems for differential equations

CERN Document Server

Hasanov Hasanoğlu, Alemdar

2017-01-01

This book presents a systematic exposition of the main ideas and methods in treating inverse problems for PDEs arising in basic mathematical models, though it makes no claim to being exhaustive. Mathematical models of most physical phenomena are governed by initial and boundary value problems for PDEs, and inverse problems governed by these equations arise naturally in nearly all branches of science and engineering. The book’s content, especially in the Introduction and Part I, is self-contained and is intended to also be accessible for beginning graduate students, whose mathematical background includes only basic courses in advanced calculus, PDEs and functional analysis. Further, the book can be used as the backbone for a lecture course on inverse and ill-posed problems for partial differential equations. In turn, the second part of the book consists of six nearly-independent chapters. The choice of these chapters was motivated by the fact that the inverse coefficient and source problems considered here a...

Obtaining sparse distributions in 2D inverse problems

OpenAIRE

Reci, A; Sederman, Andrew John; Gladden, Lynn Faith

2017-01-01

The mathematics of inverse problems has relevance across numerous estimation problems in science and engineering. L1 regularization has attracted recent attention in reconstructing the system properties in the case of sparse inverse problems; i.e., when the true property sought is not adequately described by a continuous distribution, in particular in Compressed Sensing image reconstruction. In this work, we focus on the application of L1 regularization to a class of inverse problems; relaxat...

LAI inversion from optical reflectance using a neural network trained with a multiple scattering model

Science.gov (United States)

Smith, James A.

1992-01-01

The inversion of the leaf area index (LAI) canopy parameter from optical spectral reflectance measurements is obtained using a backpropagation artificial neural network trained using input-output pairs generated by a multiple scattering reflectance model. The problem of LAI estimation over sparse canopies (LAI 1000 percent for low LAI. Minimization methods applied to merit functions constructed from differences between measured reflectances and predicted reflectances using multiple-scattering models are unacceptably sensitive to a good initial guess for the desired parameter. In contrast, the neural network reported generally yields absolute percentage errors of <30 percent when weighting coefficients trained on one soil type were applied to predicted canopy reflectance at a different soil background.

Efficient scattering-angle enrichment for a nonlinear inversion of the background and perturbations components of a velocity model

KAUST Repository

Wu, Zedong

2017-07-04

Reflection-waveform inversion (RWI) can help us reduce the nonlinearity of the standard full-waveform inversion (FWI) by inverting for the background velocity model using the wave-path of a single scattered wavefield to an image. However, current RWI implementations usually neglect the multi-scattered energy, which will cause some artifacts in the image and the update of the background. To improve existing RWI implementations in taking multi-scattered energy into consideration, we split the velocity model into background and perturbation components, integrate them directly in the wave equation, and formulate a new optimization problem for both components. In this case, the perturbed model is no longer a single-scattering model, but includes all scattering. Through introducing a new cheap implementation of scattering angle enrichment, the separation of the background and perturbation components can be implemented efficiently. We optimize both components simultaneously to produce updates to the velocity model that is nonlinear with respect to both the background and the perturbation. The newly introduced perturbation model can absorb the non-smooth update of the background in a more consistent way. We apply the proposed approach on the Marmousi model with data that contain frequencies starting from 5 Hz to show that this method can converge to an accurate velocity starting from a linearly increasing initial velocity. Also, our proposed method works well when applied to a field data set.

Size Estimates in Inverse Problems

KAUST Repository

Di Cristo, Michele

2014-01-06

Detection of inclusions or obstacles inside a body by boundary measurements is an inverse problems very useful in practical applications. When only finite numbers of measurements are available, we try to detect some information on the embedded object such as its size. In this talk we review some recent results on several inverse problems. The idea is to provide constructive upper and lower estimates of the area/volume of the unknown defect in terms of a quantity related to the work that can be expressed with the available boundary data.

Multiparameter Optimization for Electromagnetic Inversion Problem

Directory of Open Access Journals (Sweden)

M. Elkattan

2017-10-01

Full Text Available Electromagnetic (EM methods have been extensively used in geophysical investigations such as mineral and hydrocarbon exploration as well as in geological mapping and structural studies. In this paper, we developed an inversion methodology for Electromagnetic data to determine physical parameters of a set of horizontal layers. We conducted Forward model using transmission line method. In the inversion part, we solved multi parameter optimization problem where, the parameters are conductivity, dielectric constant, and permeability of each layer. The optimization problem was solved by simulated annealing approach. The inversion methodology was tested using a set of models representing common geological formations.

On two-spectra inverse problems

OpenAIRE

Guliyev, Namig J.

2018-01-01

We consider a two-spectra inverse problem for the one-dimensional Schr\\"{o}dinger equation with boundary conditions containing rational Herglotz--Nevanlinna functions of the eigenvalue parameter and provide a complete solution of this problem.

Fixed energy inversion of 5 eV e-Xe atom scattering

International Nuclear Information System (INIS)

Lovell, A.; Amos, K.

2000-01-01

Fixed energy inverse scattering theory has been used to define central and spin-orbit Schroedinger potentials for the scattering of 5 eV polarized electrons from Xe atoms. The results are typical for a range of such data; including energies above threshold when the potentials become complex. The phase shifts obtained from an analysis of the measured differential cross section and analyzing power has been used as input data. Both semi-classical (WKB) and fully quantal inversion methods have been used to extract central and spin-orbit interactions. The analysis shows that information additional to the set of input phase shifts extracted from this (and similar) data may be needed to ascertain physical potentials

Physical optics far field inverse scattering in the time domain

International Nuclear Information System (INIS)

Bleistein, N.

1976-01-01

The physical optics far field inverse scattering (POFFIS) identity relates the phase and range normalized far field back scattering amplitude to the spatial Fourier transform of the characteristic function of the scattering obstacle. The characteristic function is equal to unity in the region occupied by the obstacle and zero elsewhere. The original identity was derived by Bojarski for impulsive point sources. The result is extended to sources of arbitrary time dependence. One obtains an alternative form of Bojarski's POFFIS identity. One also derives a POFFIS identity in the time domain. Numerically synthesized checks on the method are provided

3rd Annual Workshop on Inverse Problem

CERN Document Server

2015-01-01

This proceeding volume is based on papers presented on the Third Annual Workshop on Inverse Problems which was organized by the Department of Mathematical Sciences, Chalmers University of Technology and University of Gothenburg, and took place in May 2013 in Stockholm. The purpose of this workshop was to present new analytical developments and numerical techniques for solution of inverse problems for a wide range of applications in acoustics, electromagnetics, optical fibers, medical imaging, geophysics, etc. The contributions in this volume reflect these themes and will be beneficial to researchers who are working in the area of applied inverse problems.

Ensemble Kalman methods for inverse problems

International Nuclear Information System (INIS)

Iglesias, Marco A; Law, Kody J H; Stuart, Andrew M

2013-01-01

The ensemble Kalman filter (EnKF) was introduced by Evensen in 1994 (Evensen 1994 J. Geophys. Res. 99 10143–62) as a novel method for data assimilation: state estimation for noisily observed time-dependent problems. Since that time it has had enormous impact in many application domains because of its robustness and ease of implementation, and numerical evidence of its accuracy. In this paper we propose the application of an iterative ensemble Kalman method for the solution of a wide class of inverse problems. In this context we show that the estimate of the unknown function that we obtain with the ensemble Kalman method lies in a subspace A spanned by the initial ensemble. Hence the resulting error may be bounded above by the error found from the best approximation in this subspace. We provide numerical experiments which compare the error incurred by the ensemble Kalman method for inverse problems with the error of the best approximation in A, and with variants on traditional least-squares approaches, restricted to the subspace A. In so doing we demonstrate that the ensemble Kalman method for inverse problems provides a derivative-free optimization method with comparable accuracy to that achieved by traditional least-squares approaches. Furthermore, we also demonstrate that the accuracy is of the same order of magnitude as that achieved by the best approximation. Three examples are used to demonstrate these assertions: inversion of a compact linear operator; inversion of piezometric head to determine hydraulic conductivity in a Darcy model of groundwater flow; and inversion of Eulerian velocity measurements at positive times to determine the initial condition in an incompressible fluid. (paper)

Resolvent approach for two-dimensional scattering problems. Application to the nonstationary Schrödinger problem and the KPI equation

Science.gov (United States)

Boiti, M.; Pempinelli, F.; Pogrebkov, A. K.; Polivanov, M. C.

1992-11-01

The resolvent operator of the linear problem is determined as the full Green function continued in the complex domain in two variables. An analog of the known Hilbert identity is derived. We demonstrate the role of this identity in the study of two-dimensional scattering. Considering the nonstationary Schrödinger equation as an example, we show that all types of solutions of the linear problems, as well as spectral data known in the literature, are given as specific values of this unique function — the resolvent function. A new form of the inverse problem is formulated.

Simulation of inverse Compton scattering and its implications on the scattered linewidth

Science.gov (United States)

Ranjan, N.; Terzić, B.; Krafft, G. A.; Petrillo, V.; Drebot, I.; Serafini, L.

2018-03-01

Rising interest in inverse Compton sources has increased the need for efficient models that properly quantify the behavior of scattered radiation given a set of interaction parameters. The current state-of-the-art simulations rely on Monte Carlo-based methods, which, while properly expressing scattering behavior in high-probability regions of the produced spectra, may not correctly simulate such behavior in low-probability regions (e.g. tails of spectra). Moreover, sampling may take an inordinate amount of time for the desired accuracy to be achieved. In this paper, we present an analytic derivation of the expression describing the scattered radiation linewidth and propose a model to describe the effects of horizontal and vertical emittance on the properties of the scattered radiation. We also present an improved version of the code initially reported in Krafft et al. [Phys. Rev. Accel. Beams 19, 121302 (2016), 10.1103/PhysRevAccelBeams.19.121302], that can perform the same simulations as those present in cain and give accurate results in low-probability regions by integrating over the emissions of the electrons. Finally, we use these codes to carry out simulations that closely verify the behavior predicted by the analytically derived scaling law.

The inverse conductivity problem with limited data and applications

International Nuclear Information System (INIS)

Isakov, Victor

2007-01-01

This paper describes recent uniqueness results in inverse problems for semiconductor devices and in the inverse conductivity problem. We remind basic inverse probelsm in semiconductor theory and outline use of an adjoint equation and a proof of uniqueness of piecewise constant doping profile. For the inverse conductivity problem we give a first uniqueness proof when the Dirichlet-to-Neumann map is given at an arbitrarily small part of the boundary of a three-dimensional domain

Approximation of Bayesian Inverse Problems for PDEs

OpenAIRE

Cotter, S. L.; Dashti, M.; Stuart, A. M.

2010-01-01

Inverse problems are often ill posed, with solutions that depend sensitively on data.n any numerical approach to the solution of such problems, regularization of some form is needed to counteract the resulting instability. This paper is based on an approach to regularization, employing a Bayesian formulation of the problem, which leads to a notion of well posedness for inverse problems, at the level of probability measures. The stability which results from this well posedness may be used as t...

Robust inverse scattering full waveform seismic tomography for imaging complex structure

International Nuclear Information System (INIS)

Nurhandoko, Bagus Endar B.; Sukmana, Indriani; Wibowo, Satryo; Deny, Agus; Kurniadi, Rizal; Widowati, Sri; Mubarok, Syahrul; Susilowati; Kaswandhi

2012-01-01

Seismic tomography becomes important tool recently for imaging complex subsurface. It is well known that imaging complex rich fault zone is difficult. In this paper, The application of time domain inverse scattering wave tomography to image the complex fault zone would be shown on this paper, especially an efficient time domain inverse scattering tomography and their run in cluster parallel computer which has been developed. This algorithm is purely based on scattering theory through solving Lippmann Schwienger integral by using Born's approximation. In this paper, it is shown the robustness of this algorithm especially in avoiding the inversion trapped in local minimum to reach global minimum. A large data are solved by windowing and blocking technique of memory as well as computation. Parameter of windowing computation is based on shot gather's aperture. This windowing technique reduces memory as well as computation significantly. This parallel algorithm is done by means cluster system of 120 processors from 20 nodes of AMD Phenom II. Benchmarking of this algorithm is done by means Marmoussi model which can be representative of complex rich fault area. It is shown that the proposed method can image clearly the rich fault and complex zone in Marmoussi model even though the initial model is quite far from the true model. Therefore, this method can be as one of solution to image the very complex mode.

Robust inverse scattering full waveform seismic tomography for imaging complex structure

Energy Technology Data Exchange (ETDEWEB)

Nurhandoko, Bagus Endar B.; Sukmana, Indriani; Wibowo, Satryo; Deny, Agus; Kurniadi, Rizal; Widowati, Sri; Mubarok, Syahrul; Susilowati; Kaswandhi [Wave Inversion and Subsurface Fluid Imaging Research (WISFIR) Lab., Complex System Research Division, Physics Department, Faculty of Mathematics and Natural Sciences, Institut Teknologi Bandung. and Rock Fluid Imaging Lab., Rock Physics and Cluster C (Indonesia); Rock Fluid Imaging Lab., Rock Physics and Cluster Computing Center, Bandung (Indonesia); Physics Department of Institut Teknologi Bandung (Indonesia); Rock Fluid Imaging Lab., Rock Physics and Cluster Computing Center, Bandung, Indonesia and Institut Teknologi Telkom, Bandung (Indonesia); Rock Fluid Imaging Lab., Rock Physics and Cluster Computing Center, Bandung (Indonesia)

2012-06-20

Seismic tomography becomes important tool recently for imaging complex subsurface. It is well known that imaging complex rich fault zone is difficult. In this paper, The application of time domain inverse scattering wave tomography to image the complex fault zone would be shown on this paper, especially an efficient time domain inverse scattering tomography and their run in cluster parallel computer which has been developed. This algorithm is purely based on scattering theory through solving Lippmann Schwienger integral by using Born's approximation. In this paper, it is shown the robustness of this algorithm especially in avoiding the inversion trapped in local minimum to reach global minimum. A large data are solved by windowing and blocking technique of memory as well as computation. Parameter of windowing computation is based on shot gather's aperture. This windowing technique reduces memory as well as computation significantly. This parallel algorithm is done by means cluster system of 120 processors from 20 nodes of AMD Phenom II. Benchmarking of this algorithm is done by means Marmoussi model which can be representative of complex rich fault area. It is shown that the proposed method can image clearly the rich fault and complex zone in Marmoussi model even though the initial model is quite far from the true model. Therefore, this method can be as one of solution to image the very complex mode.

Formulas in inverse and ill-posed problems

CERN Document Server

Anikonov, Yu E

1997-01-01

The Inverse and Ill-Posed Problems Series is a series of monographs publishing postgraduate level information on inverse and ill-posed problems for an international readership of professional scientists and researchers. The series aims to publish works which involve both theory and applications in, e.g., physics, medicine, geophysics, acoustics, electrodynamics, tomography, and ecology.

On the inverse problem in high-energy elastic hadron scattering and the applicability of a representation for the real part of the amplitude

International Nuclear Information System (INIS)

Fagundes, Daniel Almeida

2010-01-01

The theoretical description of high-energy elastic hadron scattering constitutes an open problem in both, the underlying quantum field theory of strong interactions (QCD) and the phenomenological context. In this work the inverse problem in elastic hadron scattering is discussed in the impact parameter and eikonal frameworks, specifically a study on the empirical extraction of the profile, the inelastic overlap and the eikonal functions, from the experimental data and some principles and high-energy theorems (model independent). The analysis is limited to elastic proton-proton scattering in the center of momentum energy interval 19.4 - 62.5 GeV. In particular, a novel representation for the Martin's Real Part Formula is introduced but without the scaling property and suitable for empirical analysis. By means of this representation, and two other parametrizations previously introduced (constrained and unconstrained), several properties of the inelastic overlap function and the imaginary part of the eikonal (opacity) in the momentum transfer space are determined, in special: (1) evidence of a peripheral effect (tail) in the inelastic overlap function in the parameter impact space above 2 fm; (2) development of analytical parametrizations for this function leading to three gaussian components with centers at 0.0, ∼0.7 and ∼1.3 fm; (3) evidence of a finite zero (change of sign) in the opacity function in the momentum transfer space; (4) development of empirical parametrization for this function consistent with form factors as a product of two monopoles with constrained masses (not a dipole type) and a term with zero; (5) detailed discussion on the determination of the opacity function in the momentum transfer space through the semi-analytical approach. The applicability of these empirical results in the development of eikonal models (mainly those inspired in QCD) is also discussed. (author)

A Survey on Inverse Problems for Applied Sciences

Directory of Open Access Journals (Sweden)

Fatih Yaman

2013-01-01

Full Text Available The aim of this paper is to introduce inversion-based engineering applications and to investigate some of the important ones from mathematical point of view. To do this we employ acoustic, electromagnetic, and elastic waves for presenting different types of inverse problems. More specifically, we first study location, shape, and boundary parameter reconstruction algorithms for the inaccessible targets in acoustics. The inverse problems for the time-dependent differential equations of isotropic and anisotropic elasticity are reviewed in the following section of the paper. These problems were the objects of the study by many authors in the last several decades. The physical interpretations for almost all of these problems are given, and the geophysical applications for some of them are described. In our last section, an introduction with many links into the literature is given for modern algorithms which combine techniques from classical inverse problems with stochastic tools into ensemble methods both for data assimilation as well as for forecasting.

Inverse Raman scattering in silicon: A free-carrier enhanced effect

International Nuclear Information System (INIS)

Solli, D. R.; Koonath, P.; Jalali, B.

2009-01-01

Stimulated Raman scattering has been harnessed to produce the first silicon lasers and amplifiers. The Raman effect can also produce intensity-dependent nonlinear loss through a corollary process, inverse Raman scattering (IRS). This process has never been observed in a semiconductor. We demonstrate IRS in silicon--a process that is substantially modified by optically generated free carriers--achieving attenuation levels >15 dB with a pump intensity of 4 GW/cm 2 . Surprisingly, free-carrier absorption, the detrimental effect that generally suppresses nonlinear effects in silicon, actually facilitates IRS by delaying the onset of contamination from coherent anti-Stokes Raman scattering. Silicon-based IRS could be a valuable tool for chip-scale signal processing.

Reconstruction Methods for Inverse Problems with Partial Data

DEFF Research Database (Denmark)

Hoffmann, Kristoffer

This thesis presents a theoretical and numerical analysis of a general mathematical formulation of hybrid inverse problems in impedance tomography. This includes problems from several existing hybrid imaging modalities such as Current Density Impedance Imaging, Magnetic Resonance Electrical...... Impedance Tomography, and Ultrasound Modulated Electrical Impedance Tomography. After giving an introduction to hybrid inverse problems in impedance tomography and the mathematical tools that facilitate the related analysis, we explain in detail the stability properties associated with the classification...... of a linearised hybrid inverse problem. This is done using pseudo-differential calculus and theory for overdetermined boundary value problem. Using microlocal analysis we then present novel results on the propagation of singularities, which give a precise description of the distinct features of solutions...

Gradient-type methods in inverse parabolic problems

International Nuclear Information System (INIS)

Kabanikhin, Sergey; Penenko, Aleksey

2008-01-01

This article is devoted to gradient-based methods for inverse parabolic problems. In the first part, we present a priori convergence theorems based on the conditional stability estimates for linear inverse problems. These theorems are applied to backwards parabolic problem and sideways parabolic problem. The convergence conditions obtained coincide with sourcewise representability in the self-adjoint backwards parabolic case but they differ in the sideways case. In the second part, a variational approach is formulated for a coefficient identification problem. Using adjoint equations, a formal gradient of an objective functional is constructed. A numerical test illustrates the performance of conjugate gradient algorithm with the formal gradient.

PREFACE: Inverse Problems in Applied Sciences—towards breakthrough

Science.gov (United States)

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

Introduction to the 30th volume of Inverse Problems

Science.gov (United States)

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

New method for solving multidimensional scattering problem

International Nuclear Information System (INIS)

Melezhik, V.S.

1991-01-01

A new method is developed for solving the quantum mechanical problem of scattering of a particle with internal structure. The multichannel scattering problem is formulated as a system of nonlinear functional equations for the wave function and reaction matrix. The method is successfully tested for the scattering from a nonspherical potential well and a long-range nonspherical scatterer. The method is also applicable to solving the multidimensional Schroedinger equation with a discrete spectrum. As an example the known problem of a hydrogen atom in a homogeneous magnetic field is analyzed

The inside–outside duality for inverse scattering problems with near field data

International Nuclear Information System (INIS)

Lechleiter, Armin; Peters, Stefan

2015-01-01

We derive an inside–outside duality for near field scattering data generated by time-harmonic scattering of acoustic point sources from a sound-soft scatterer. This duality in particular rigorously characterizes interior Dirichlet eigenvalues of the scattering object by near field operators for an interval of wave numbers. As a crucial new concept to prove this duality we exploit the numerical ranges of certain modifications of these near field operators. We also show that our theoretical results can be numerically used to approximate interior Dirichlet eigenvalues from multi-frequency near field measurements. (paper)

Inverse Modelling Problems in Linear Algebra Undergraduate Courses

Science.gov (United States)

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…

Tropospheric nitrogen dioxide inversions based on spectral measurements of scattered sunlight

NARCIS (Netherlands)

Vlemmix, T.

2011-01-01

This thesis describes the development of inversion methods for tropospheric nitrogen dioxide (NO2), based on ground based observations of scattered sunlight with themulti-axis differential optical absorption spectroscopy (MAX-DOAS) technique. NO2 is an atmospheric trace gas which, when present near

Particle Swarm Optimization and Uncertainty Assessment in Inverse Problems

Directory of Open Access Journals (Sweden)

José L. G. Pallero

2018-01-01

Full Text Available Most inverse problems in the industry (and particularly in geophysical exploration are highly underdetermined because the number of model parameters too high to achieve accurate data predictions and because the sampling of the data space is scarce and incomplete; it is always affected by different kinds of noise. Additionally, the physics of the forward problem is a simplification of the reality. All these facts result in that the inverse problem solution is not unique; that is, there are different inverse solutions (called equivalent, compatible with the prior information that fits the observed data within similar error bounds. In the case of nonlinear inverse problems, these equivalent models are located in disconnected flat curvilinear valleys of the cost-function topography. The uncertainty analysis consists of obtaining a representation of this complex topography via different sampling methodologies. In this paper, we focus on the use of a particle swarm optimization (PSO algorithm to sample the region of equivalence in nonlinear inverse problems. Although this methodology has a general purpose, we show its application for the uncertainty assessment of the solution of a geophysical problem concerning gravity inversion in sedimentary basins, showing that it is possible to efficiently perform this task in a sampling-while-optimizing mode. Particularly, we explain how to use and analyze the geophysical models sampled by exploratory PSO family members to infer different descriptors of nonlinear uncertainty.

Inverse planning and optimization: a comparison of solutions

Energy Technology Data Exchange (ETDEWEB)

Ringor, Michael [School of Health Sciences, Purdue University, West Lafayette, IN (United States); Papiez, Lech [Department of Radiation Oncology, Indiana University, Indianapolis, IN (United States)

1998-09-01

The basic problem in radiation therapy treatment planning is to determine an appropriate set of treatment parameters that would induce an effective dose distribution inside a patient. One can approach this task as an inverse problem, or as an optimization problem. In this presentation, we compare both approaches. The inverse problem is presented as a dose reconstruction problem similar to tomography reconstruction. We formulate the optimization problem as linear and quadratic programs. Explicit comparisons are made between the solutions obtained by inversion and those obtained by optimization for the case in which scatter and attenuation are ignored (the NS-NA approximation)

6th International Workshop on New Computational Methods for Inverse Problems

International Nuclear Information System (INIS)

2016-01-01

methods, learning methods, convex optimization, free discontinuity problems, metamodels, proper orthogonal decomposition, reduced models for the inversion, non-linear inverse scattering, image reconstruction and restoration, and applications (bio-medical imaging, nondestructive evaluation...). NCMIP 2016 was a one-day workshop held in May 2016 which attracted around seventy attendees. Each of the submitted papers has been reviewed by two reviewers. There have been eleven accepted papers. In addition, three international speakers were invited to present a longer talk. The workshop was supported by Institut Farman (ENS Cachan, CNRS) and endorsed by the following French research networks: GDR ISIS, GDR MIA, GDR MOA, GDR Ondes. The program committee acknowledges the following research laboratories: CMLA, LMT, LURPA and SATIE. Eric Vourc'h and Thomas Rodet Workshop co-chairs: Eric Vourc'h, SATIE laboratory, Ecole Normale Supérieure de Cachan, CNRS, France Thomas Rodet, SATIE laboratory, Ecole Normale Supérieure de Cachan, CNRS, France Technical program committee: Alexandre Baussard, ENSTA Bretagne, Lab-STICC, France Marc Bonnet, ENSTA, ParisTech, France Laure Blanc-Féraud, I3S laboratory and INRIA Nice Sophia-Antipolis, France Antonin Chambolle, CMAP, Ecole Polytechnique, CNRS, France Oliver Dorn, School of Mathematics, University of Manchester, UK Cécile Durieu, SATIE, ENS Cachan, CNRS, France Laurent Fribourg, LSV, ENS Cachan, CNRS, France Jerôme Idier, IRCCyN Laboratory, Ecole Centrale de Nantes, France Pierre-Yves Joubert, IEF, Paris-Sud University, CNRS, France Marc Lambert, Geeps Laboratory, CNRS, CentraleSupElec, Paris-Sud University, France Giacomo Oliveri, eledia research center/eledia@L2S group, University of Trento, Italy Dominique Lesselier, L2S Laboratory, CNRS, CentraleSupElec, Paris-Sud University, France Matteo Pastorino, DIBE, University of Genoa, Italy Gabriel Peyré, Ceremade laboratory, University of Paris Dauphine, France Anthony Quinn

One-dimensional inverse problems of mathematical physics

CERN Document Server

Lavrent'ev, M M; Yakhno, V G; Schulenberger, J R

1986-01-01

This monograph deals with the inverse problems of determining a variable coefficient and right side for hyperbolic and parabolic equations on the basis of known solutions at fixed points of space for all times. The problems are one-dimensional in nature since the desired coefficient of the equation is a function of only one coordinate, while the desired right side is a function only of time. The authors use methods based on the spectral theory of ordinary differential operators of second order and also methods which make it possible to reduce the investigation of the inverse problems to the in

Inverse scattering solution of non-linear evolution equations in one space dimension: an introduction

International Nuclear Information System (INIS)

Alvarez-Estrada, R.F.

1979-01-01

A comprehensive review of the inverse scattering solution of certain non-linear evolution equations of physical interest in one space dimension is presented. We explain in some detail the interrelated techniques which allow to linearize exactly the following equations: (1) the Korteweg and de Vries equation; (2) the non-linear Schrodinger equation; (3) the modified Korteweg and de Vries equation; (4) the Sine-Gordon equation. We concentrate in discussing the pairs of linear operators which accomplish such an exact linearization and the solution of the associated initial value problem. The application of the method to other non-linear evolution equations is reviewed very briefly

Non-local currents in 2D QFT: an alternative To - the quantum inverse scattering method

International Nuclear Information System (INIS)

Bernard, D.; Leclair, A.; Cornell Univ., Ithaca, NY

1990-01-01

The formalism based on non-local charges that we propose provides an alternative to the quantum inverse scattering method for solving integrable quantum field theories in 2D. The content of the paper is: 1. Introduction: historical background. 2. The NLC approach to 2D QFT: a summary. 3 Exchange algebras and on-shell conservation laws: why non-local charges are useful. 4. The lattice construction: the geometrical origin of non-local conserved currents. 5. The continuum construction: how to deal with non-local conserved currents. 6. Examples: Yangian and quantum group currents. 7 Conclusions: open problems. 22 refs., 4 figs

Size Estimates in Inverse Problems

KAUST Repository

Di Cristo, Michele

2014-01-01

Detection of inclusions or obstacles inside a body by boundary measurements is an inverse problems very useful in practical applications. When only finite numbers of measurements are available, we try to detect some information on the embedded

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

Science.gov (United States)

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

On reciprocal Baecklund transformations of inverse scattering schemes

International Nuclear Information System (INIS)

Rogers, C.; Wong, P.

1984-01-01

The notion of reciprocally related inverse scattering schemes is introduced and is shown to be a key component in the link between the AKNS and WKI schemes. Reciprocal auto-Baecklund transformations are represented both for a generalised Harry-Dym equation and an equation descriptive of nonlinear oscillation of elastic beams. Further, the N-loop soliton solution of the KIW equation is generated in a convenient parametric form via reciprocal Baecklund transformations. Finally, an important reduction to canonical spectral form is shown to be a reciprocal transformation. (Auth.)

Inverse problems in vision and 3D tomography

CERN Document Server

Mohamad-Djafari, Ali

2013-01-01

The concept of an inverse problem is a familiar one to most scientists and engineers, particularly in the field of signal and image processing, imaging systems (medical, geophysical, industrial non-destructive testing, etc.) and computer vision. In imaging systems, the aim is not just to estimate unobserved images, but also their geometric characteristics from observed quantities that are linked to these unobserved quantities through the forward problem. This book focuses on imagery and vision problems that can be clearly written in terms of an inverse problem where an estimate for the image a

Inverse problems basics, theory and applications in geophysics

CERN Document Server

Richter, Mathias

2016-01-01

The overall goal of the book is to provide access to the regularized solution of inverse problems relevant in geophysics without requiring more mathematical knowledge than is taught in undergraduate math courses for scientists and engineers. From abstract analysis only the concept of functions as vectors is needed. Function spaces are introduced informally in the course of the text, when needed. Additionally, a more detailed, but still condensed introduction is given in Appendix B. A second goal is to elaborate the single steps to be taken when solving an inverse problem: discretization, regularization and practical solution of the regularized optimization problem. These steps are shown in detail for model problems from the fields of inverse gravimetry and seismic tomography. The intended audience is mathematicians, physicists and engineers having a good working knowledge of linear algebra and analysis at the upper undergraduate level.

REGULARIZED D-BAR METHOD FOR THE INVERSE CONDUCTIVITY PROBLEM

DEFF Research Database (Denmark)

Knudsen, Kim; Lassas, Matti; Mueller, Jennifer

2009-01-01

A strategy for regularizing the inversion procedure for the two-dimensional D-bar reconstruction algorithm based on the global uniqueness proof of Nachman [Ann. Math. 143 (1996)] for the ill-posed inverse conductivity problem is presented. The strategy utilizes truncation of the boundary integral...... the convergence of the reconstructed conductivity to the true conductivity as the noise level tends to zero. The results provide a link between two traditions of inverse problems research: theory of regularization and inversion methods based on complex geometrical optics. Also, the procedure is a novel...

Solving inverse two-point boundary value problems using collage coding

Science.gov (United States)

Kunze, H.; Murdock, S.

2006-08-01

The method of collage coding, with its roots in fractal imaging, is the central tool in a recently established rigorous framework for solving inverse initial value problems for ordinary differential equations (Kunze and Vrscay 1999 Inverse Problems 15 745-70). We extend these ideas to solve the following inverse problem: given a function u(x) on [A, B] (which may be the interpolation of data points), determine a two-point boundary value problem on [A, B] which admits u(x) as a solution as closely as desired. The solution of such inverse problems may be useful in parameter estimation or determination of potential functional forms of the underlying differential equation. We discuss ways to improve results, including the development of a partitioning scheme. Several examples are considered.

Modeling of uncertainties in statistical inverse problems

International Nuclear Information System (INIS)

Kaipio, Jari

2008-01-01

In all real world problems, the models that tie the measurements to the unknowns of interest, are at best only approximations for reality. While moderate modeling and approximation errors can be tolerated with stable problems, inverse problems are a notorious exception. Typical modeling errors include inaccurate geometry, unknown boundary and initial data, properties of noise and other disturbances, and simply the numerical approximations of the physical models. In principle, the Bayesian approach to inverse problems, in which all uncertainties are modeled as random variables, is capable of handling these uncertainties. Depending on the type of uncertainties, however, different strategies may be adopted. In this paper we give an overview of typical modeling errors and related strategies within the Bayesian framework.

Recently Developed Formulations of the Inverse Problem in Acoustics and Electromagnetics

Science.gov (United States)

1974-12-01

solution for scattering by a sphere. The inverse transform of irs?(K) is calculated, this function yielding --y (x). Figure 4.2 is a graph of this...time or decays "sufficiently rapidly", then T+- o. In this case, we may let T -1 in (8.9) and obtain the inverse transform (k = w/c) of (5.6) as the

Inverse scattering and solitons in An-1 affine Toda field theories

International Nuclear Information System (INIS)

Beggs, E.J.; Johnson, P.R.

1997-01-01

We implement the inverse scattering method in the case of the A n affine Toda field theories, by studying the space-time evolution of simple poles in the underlying loop group. We find the known single-soliton solutions, as well as additional solutions with non-linear modes of oscillation around the standard solution, by studying the particularly simple case where the residue at the pole is a rank-one projection. We show that these solutions with extra modes have the same mass and topological charges as the standard solutions, so we do not shed any light on the missing topological charge problem in these models. From the monodromy matrix it is shown that these solutions have the same higher conserved charges as the standard solutions. We also show that the integrated energy-momentum density can be calculated from the central extension of the loop group. (orig.)

Inverse problems in ordinary differential equations and applications

CERN Document Server

Llibre, Jaume

2016-01-01

This book is dedicated to study the inverse problem of ordinary differential equations, that is it focuses in finding all ordinary differential equations that satisfy a given set of properties. The Nambu bracket is the central tool in developing this approach. The authors start characterizing the ordinary differential equations in R^N which have a given set of partial integrals or first integrals. The results obtained are applied first to planar polynomial differential systems with a given set of such integrals, second to solve the 16th Hilbert problem restricted to generic algebraic limit cycles, third for solving the inverse problem for constrained Lagrangian and Hamiltonian mechanical systems, fourth for studying the integrability of a constrained rigid body. Finally the authors conclude with an analysis on nonholonomic mechanics, a generalization of the Hamiltonian principle, and the statement an solution of the inverse problem in vakonomic mechanics.

A simple method for solving the inverse scattering problem

International Nuclear Information System (INIS)

Melnikov, V.N.; Rudyak, B.V.; Zakhariev, V.N.

1977-01-01

A new method is proposed for approximate reconstruction of a potential as a step function from scattering data using the completeness relation of solutions of the Schroedinger equation. The suggested method allows one to take into account exactly the additional centrifugal barrier for partial waves with angular momentum l>0, and also the Coulomb potential. The method admits different generalizations. Numerical calculations for checking the method have been performed

An inverse problem approach to pattern recognition in industry

Directory of Open Access Journals (Sweden)

Ali Sever

2015-01-01

Full Text Available Many works have shown strong connections between learning and regularization techniques for ill-posed inverse problems. A careful analysis shows that a rigorous connection between learning and regularization for inverse problem is not straightforward. In this study, pattern recognition will be viewed as an ill-posed inverse problem and applications of methods from the theory of inverse problems to pattern recognition are studied. A new learning algorithm derived from a well-known regularization model is generated and applied to the task of reconstruction of an inhomogeneous object as pattern recognition. Particularly, it is demonstrated that pattern recognition can be reformulated in terms of inverse problems defined by a Riesz-type kernel. This reformulation can be employed to design a learning algorithm based on a numerical solution of a system of linear equations. Finally, numerical experiments have been carried out with synthetic experimental data considering a reasonable level of noise. Good recoveries have been achieved with this methodology, and the results of these simulations are compatible with the existing methods. The comparison results show that the Regularization-based learning algorithm (RBA obtains a promising performance on the majority of the test problems. In prospects, this method can be used for the creation of automated systems for diagnostics, testing, and control in various fields of scientific and applied research, as well as in industry.

A tutorial on inverse problems for anomalous diffusion processes

International Nuclear Information System (INIS)

Jin, Bangti; Rundell, William

2015-01-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

The inverse spectral problem for pencils of differential operators

International Nuclear Information System (INIS)

Guseinov, I M; Nabiev, I M

2007-01-01

The inverse problem of spectral analysis for a quadratic pencil of Sturm-Liouville operators on a finite interval is considered. A uniqueness theorem is proved, a solution algorithm is presented, and sufficient conditions for the solubility of the inverse problem are obtained. Bibliography: 31 titles.

Applications of elliptic Carleman inequalities to Cauchy and inverse problems

CERN Document Server

Choulli, Mourad

2016-01-01

This book presents a unified approach to studying the stability of both elliptic Cauchy problems and selected inverse problems. Based on elementary Carleman inequalities, it establishes three-ball inequalities, which are the key to deriving logarithmic stability estimates for elliptic Cauchy problems and are also useful in proving stability estimates for certain elliptic inverse problems. The book presents three inverse problems, the first of which consists in determining the surface impedance of an obstacle from the far field pattern. The second problem investigates the detection of corrosion by electric measurement, while the third concerns the determination of an attenuation coefficient from internal data, which is motivated by a problem encountered in biomedical imaging.

Spectral Green’s function nodal method for multigroup SN problems with anisotropic scattering in slab-geometry non-multiplying media

International Nuclear Information System (INIS)

Menezes, Welton A.; Filho, Hermes Alves; Barros, Ricardo C.

2014-01-01

Highlights: • Fixed-source S N transport problems. • Energy multigroup model. • Anisotropic scattering. • Slab-geometry spectral nodal method. - Abstract: A generalization of the spectral Green’s function (SGF) method is developed for multigroup, fixed-source, slab-geometry discrete ordinates (S N ) problems with anisotropic scattering. The offered SGF method with the one-node block inversion (NBI) iterative scheme converges numerical solutions that are completely free from spatial truncation errors for multigroup, slab-geometry S N problems with scattering anisotropy of order L, provided L < N. As a coarse-mesh numerical method, the SGF method generates numerical solutions that generally do not give detailed information on the problem solution profile, as the grid points can be located considerably away from each other. Therefore, we describe in this paper a technique for the spatial reconstruction of the coarse-mesh solution generated by the multigroup SGF method. Numerical results are given to illustrate the method’s accuracy

A flexible, extendable, modular and computationally efficient approach to scattering-integral-based seismic full waveform inversion

Science.gov (United States)

Schumacher, F.; Friederich, W.; Lamara, S.

2016-02-01

We present a new conceptual approach to scattering-integral-based seismic full waveform inversion (FWI) that allows a flexible, extendable, modular and both computationally and storage-efficient numerical implementation. To achieve maximum modularity and extendability, interactions between the three fundamental steps carried out sequentially in each iteration of the inversion procedure, namely, solving the forward problem, computing waveform sensitivity kernels and deriving a model update, are kept at an absolute minimum and are implemented by dedicated interfaces. To realize storage efficiency and maximum flexibility, the spatial discretization of the inverted earth model is allowed to be completely independent of the spatial discretization employed by the forward solver. For computational efficiency reasons, the inversion is done in the frequency domain. The benefits of our approach are as follows: (1) Each of the three stages of an iteration is realized by a stand-alone software program. In this way, we avoid the monolithic, unflexible and hard-to-modify codes that have often been written for solving inverse problems. (2) The solution of the forward problem, required for kernel computation, can be obtained by any wave propagation modelling code giving users maximum flexibility in choosing the forward modelling method. Both time-domain and frequency-domain approaches can be used. (3) Forward solvers typically demand spatial discretizations that are significantly denser than actually desired for the inverted model. Exploiting this fact by pre-integrating the kernels allows a dramatic reduction of disk space and makes kernel storage feasible. No assumptions are made on the spatial discretization scheme employed by the forward solver. (4) In addition, working in the frequency domain effectively reduces the amount of data, the number of kernels to be computed and the number of equations to be solved. (5) Updating the model by solving a large equation system can be

Inverse source problems in elastodynamics

Science.gov (United States)

Bao, Gang; Hu, Guanghui; Kian, Yavar; Yin, Tao

2018-04-01

We are concerned with time-dependent inverse source problems in elastodynamics. The source term is supposed to be the product of a spatial function and a temporal function with compact support. We present frequency-domain and time-domain approaches to show uniqueness in determining the spatial function from wave fields on a large sphere over a finite time interval. The stability estimate of the temporal function from the data of one receiver and the uniqueness result using partial boundary data are proved. Our arguments rely heavily on the use of the Fourier transform, which motivates inversion schemes that can be easily implemented. A Landweber iterative algorithm for recovering the spatial function and a non-iterative inversion scheme based on the uniqueness proof for recovering the temporal function are proposed. Numerical examples are demonstrated in both two and three dimensions.

Inverse planning for x-ray rotation therapy: a general solution of the inverse problem

International Nuclear Information System (INIS)

Oelfke, U.; Bortfeld, T.

1999-01-01

Rotation therapy with photons is currently under investigation for the delivery of intensity modulated radiotherapy (IMRT). An analytical approach for inverse treatment planning of this radiotherapy technique is described. The inverse problem for the delivery of arbitrary 2D dose profiles is first formulated and then solved analytically. In contrast to previously applied strategies for solving the inverse problem, it is shown that the most general solution for the fluence profiles consists of two independent solutions of different parity. A first analytical expression for both fluence profiles is derived. The mathematical derivation includes two different strategies, an elementary expansion of fluence and dose into polynomials and a more practical approach in terms of Fourier transforms. The obtained results are discussed in the context of previous work on this problem. (author)

Acoustic inverse scattering using topological derivative of far-field measurements-based L2 cost functionals

International Nuclear Information System (INIS)

Bellis, Cédric; Bonnet, Marc; Cakoni, Fioralba

2013-01-01

Originally formulated in the context of topology optimization, the concept of topological derivative has also proved effective as a qualitative inversion tool for a wave-based identification of finite-sized objects. This approach remains, however, largely based on a heuristic interpretation of the topological derivative, whereas most other qualitative approaches to inverse scattering are backed by a mathematical justification. As an effort toward bridging this gap, this study focuses on a topological derivative approach applied to the L 2 -norm of the misfit between far-field measurements. Either an inhomogeneous medium or a finite number of point-like scatterers are considered, using either the Born approximation or a full-scattering model. Topological derivative-based imaging functionals are analyzed using a suitable factorization of the far-field operator, for each of the considered cases, in order to characterize their behavior and assess their ability to reconstruct the unknown scatterer(s). Results include the justification of the usual sign heuristic underpinning the method for (i) the Born approximation and (ii) full-scattering models limited to moderately strong scatterers. Semi-analytical and numerical examples are presented. Within the chosen framework, the topological derivative approach is finally discussed and compared to other well-known qualitative methods. (paper)

Integral geometry and inverse problems for hyperbolic equations

CERN Document Server

Romanov, V G

1974-01-01

There are currently many practical situations in which one wishes to determine the coefficients in an ordinary or partial differential equation from known functionals of its solution. These are often called "inverse problems of mathematical physics" and may be contrasted with problems in which an equation is given and one looks for its solution under initial and boundary conditions. Although inverse problems are often ill-posed in the classical sense, their practical importance is such that they may be considered among the pressing problems of current mathematical re­ search. A. N. Tihonov showed [82], [83] that there is a broad class of inverse problems for which a particular non-classical definition of well-posed ness is appropriate. This new definition requires that a solution be unique in a class of solutions belonging to a given subset M of a function space. The existence of a solution in this set is assumed a priori for some set of data. The classical requirement of continuous dependence of the solutio...

The inverse problem of the magnetostatic nondestructive testing

International Nuclear Information System (INIS)

Pechenkov, A.N.; Shcherbinin, V.E.

2006-01-01

The inverse problem of magnetostatic nondestructive testing consists in the calculation of the shape and magnetic characteristics of a flaw in a uniform magnetized body with measurement of static magnetic field beyond the body. If the flaw does not contain any magnetic material, the inverse problem is reduced to identification of the shape and magnetic susceptibility of the substance. This case has been considered in the study [ru

Inverse problems in systems biology

International Nuclear Information System (INIS)

Engl, Heinz W; Lu, James; Müller, Stefan; Flamm, Christoph; Schuster, Peter; Kügler, Philipp

2009-01-01

Systems biology is a new discipline built upon the premise that an understanding of how cells and organisms carry out their functions cannot be gained by looking at cellular components in isolation. Instead, consideration of the interplay between the parts of systems is indispensable for analyzing, modeling, and predicting systems' behavior. Studying biological processes under this premise, systems biology combines experimental techniques and computational methods in order to construct predictive models. Both in building and utilizing models of biological systems, inverse problems arise at several occasions, for example, (i) when experimental time series and steady state data are used to construct biochemical reaction networks, (ii) when model parameters are identified that capture underlying mechanisms or (iii) when desired qualitative behavior such as bistability or limit cycle oscillations is engineered by proper choices of parameter combinations. In this paper we review principles of the modeling process in systems biology and illustrate the ill-posedness and regularization of parameter identification problems in that context. Furthermore, we discuss the methodology of qualitative inverse problems and demonstrate how sparsity enforcing regularization allows the determination of key reaction mechanisms underlying the qualitative behavior. (topical review)

The philosophical aspect of learning inverse problems of mathematical physics

Directory of Open Access Journals (Sweden)

Виктор Семенович Корнилов

2018-12-01

Full Text Available The article describes specific questions student learning inverse problems of mathematical physics. When teaching inverse problems of mathematical physics to the understanding of the students brought the information that the inverse problems of mathematical physics with a philosophical point of view are the problems of determining the unknown causes of known consequences, and the search for their solutions have great scientific and educational potential. The reasons are specified in the form of unknown coefficients, right side, initial conditions of the mathematical model of inverse problems, and as a consequence are functionals of the solution of this mathematical model. In the process of learning the inverse problems of mathematical physics focuses on the philosophical aspects of the phenomenon of information and identify cause-effect relations. It is emphasized that in the process of logical analysis applied and humanitarian character, students realize that information is always related to the fundamental philosophical questions that the analysis applied and the humanitarian aspects of the obtained results the inverse problem of mathematical physics allows students to make appropriate inferences about the studied process and to, ultimately, new information, to study its properties and understand its value. Philosophical understanding of the notion of information opens up to students a new methodological opportunities to comprehend the world and helps us to reinterpret existing science and philosophy of the theory related to the disclosure of the interrelationship of all phenomena of reality.

On the use of the Reciprocity Gap Functional in inverse scattering with near-field data: An application to mammography

International Nuclear Information System (INIS)

Delbary, Fabrice; Piana, Michele; Aramini, Riccardo; Brignone, Massimo; Bozza, Giovanni

2008-01-01

Microwave tomography is a non-invasive approach to the early diagnosis of breast cancer. However the problem of visualizing tumors from diffracted microwaves is a difficult nonlinear ill-posed inverse scattering problem. We propose a qualitative approach to the solution of such a problem, whereby the shape and location of cancerous tissues can be detected by means of a combination of the Reciprocity Gap Functional method and the Linear Sampling method. We validate this approach to synthetic near-fields produced by a finite element method for boundary integral equations, where the breast is mimicked by the axial view of two nested cylinders, the external one representing the skin and the internal one representing the fat tissue.

On the use of the Reciprocity Gap Functional in inverse scattering with near-field data: An application to mammography

Science.gov (United States)

Delbary, Fabrice; Aramini, Riccardo; Bozza, Giovanni; Brignone, Massimo; Piana, Michele

2008-11-01

Microwave tomography is a non-invasive approach to the early diagnosis of breast cancer. However the problem of visualizing tumors from diffracted microwaves is a difficult nonlinear ill-posed inverse scattering problem. We propose a qualitative approach to the solution of such a problem, whereby the shape and location of cancerous tissues can be detected by means of a combination of the Reciprocity Gap Functional method and the Linear Sampling method. We validate this approach to synthetic near-fields produced by a finite element method for boundary integral equations, where the breast is mimicked by the axial view of two nested cylinders, the external one representing the skin and the internal one representing the fat tissue.

Inverse scattering transform method and soliton solutions for Davey-Stewartson II equation

International Nuclear Information System (INIS)

Arkadiev, V.A.; Pogrebkov, A.K.; Polivanov, M.C.

1989-01-01

The inverse scattering method for Davey-Stewartson II (DS-II) equation including both soliton and continuous spectrum solutions is developed. The explicit formulae for N-soliton solutions are given. Note that our solitons decrease as |z| -2 with z tending to infinity. (author). 8 refs

An inverse problem for evolution inclusions

OpenAIRE

Ton, Bui An

2002-01-01

An inverse problem, the determination of the shape and a convective coefficient on a part of the boundary from partial measurements of the solution, is studied using 2-person optimal control techniques.

Level set methods for inverse scattering—some recent developments

International Nuclear Information System (INIS)

Dorn, Oliver; Lesselier, Dominique

2009-01-01

We give an update on recent techniques which use a level set representation of shapes for solving inverse scattering problems, completing in that matter the exposition made in (Dorn and Lesselier 2006 Inverse Problems 22 R67) and (Dorn and Lesselier 2007 Deformable Models (New York: Springer) pp 61–90), and bringing it closer to the current state of the art

Solving inversion problems with neural networks

Science.gov (United States)

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.

Multi-parameter Analysis and Inversion for Anisotropic Media Using the Scattering Integral Method

KAUST Repository

Djebbi, Ramzi

2017-01-01

the model. I study the prospect of applying a scattering integral approach for multi-parameter inversion for a transversely isotropic model with a vertical axis of symmetry. I mainly analyze the sensitivity kernels to understand the sensitivity of seismic

Inverse problems in classical and quantum physics

International Nuclear Information System (INIS)

Almasy, A.A.

2007-01-01

The subject of this thesis is in the area of Applied Mathematics known as Inverse Problems. Inverse problems are those where a set of measured data is analysed in order to get as much information as possible on a model which is assumed to represent a system in the real world. We study two inverse problems in the fields of classical and quantum physics: QCD condensates from tau-decay data and the inverse conductivity problem. Despite a concentrated effort by physicists extending over many years, an understanding of QCD from first principles continues to be elusive. Fortunately, data continues to appear which provide a rather direct probe of the inner workings of the strong interactions. We use a functional method which allows us to extract within rather general assumptions phenomenological parameters of QCD (the condensates) from a comparison of the time-like experimental data with asymptotic space-like results from theory. The price to be paid for the generality of assumptions is relatively large errors in the values of the extracted parameters. Although we do not claim that our method is superior to other approaches, we hope that our results lend additional confidence to the numerical results obtained with the help of methods based on QCD sum rules. EIT is a technology developed to image the electrical conductivity distribution of a conductive medium. The technique works by performing simultaneous measurements of direct or alternating electric currents and voltages on the boundary of an object. These are the data used by an image reconstruction algorithm to determine the electrical conductivity distribution within the object. In this thesis, two approaches of EIT image reconstruction are proposed. The first is based on reformulating the inverse problem in terms of integral equations. This method uses only a single set of measurements for the reconstruction. The second approach is an algorithm based on linearisation which uses more then one set of measurements. A

Inverse problems in classical and quantum physics

Energy Technology Data Exchange (ETDEWEB)

Almasy, A.A.

2007-06-29

The subject of this thesis is in the area of Applied Mathematics known as Inverse Problems. Inverse problems are those where a set of measured data is analysed in order to get as much information as possible on a model which is assumed to represent a system in the real world. We study two inverse problems in the fields of classical and quantum physics: QCD condensates from tau-decay data and the inverse conductivity problem. Despite a concentrated effort by physicists extending over many years, an understanding of QCD from first principles continues to be elusive. Fortunately, data continues to appear which provide a rather direct probe of the inner workings of the strong interactions. We use a functional method which allows us to extract within rather general assumptions phenomenological parameters of QCD (the condensates) from a comparison of the time-like experimental data with asymptotic space-like results from theory. The price to be paid for the generality of assumptions is relatively large errors in the values of the extracted parameters. Although we do not claim that our method is superior to other approaches, we hope that our results lend additional confidence to the numerical results obtained with the help of methods based on QCD sum rules. EIT is a technology developed to image the electrical conductivity distribution of a conductive medium. The technique works by performing simultaneous measurements of direct or alternating electric currents and voltages on the boundary of an object. These are the data used by an image reconstruction algorithm to determine the electrical conductivity distribution within the object. In this thesis, two approaches of EIT image reconstruction are proposed. The first is based on reformulating the inverse problem in terms of integral equations. This method uses only a single set of measurements for the reconstruction. The second approach is an algorithm based on linearisation which uses more then one set of measurements. A

Quantal inversion of cross-section for the elastic scattering of 200 MeV protons from 12C

International Nuclear Information System (INIS)

Allen, L.J.; Amos, K.; Dortmans, P.J.

1994-01-01

Fixed energy quantal inverse scattering theory has been used to analyse the differential cross-section from the elastic scattering of 200 MeV protons from 12 C. Ambiguities in obtaining the scattering function from the differential cross-section are discussed and by means of example it is illustrated that not all scattering functions lead to physically reasonable potentials. 8 refs., 2 tabs., 4 figs

Inverse kinematics problem in robotics using neural networks

Science.gov (United States)

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.

Solution accelerators for large scale 3D electromagnetic inverse problems

International Nuclear Information System (INIS)

Newman, Gregory A.; Boggs, Paul T.

2004-01-01

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

Inverse Scattering, the Coupling Constant Spectrum, and the Riemann Hypothesis

International Nuclear Information System (INIS)

Khuri, N. N.

2002-01-01

It is well known that the s-wave Jost function for a potential, λV, is an entire function of λ with an infinite number of zeros extending to infinity. For a repulsive V, and at zero energy, these zeros of the 'coupling constant', λ, will all be real and negative, λ n (0) n n =1/2+iγ n . Thus, finding a repulsive V whose coupling constant spectrum coincides with the Riemann zeros will establish the Riemann hypothesis, but this will be a very difficult and unguided search.In this paper we make a significant enlargement of the class of potentials needed for a generalization of the above idea. We also make this new class amenable to construction via inverse scattering methods. We show that all one needs is a one parameter class of potentials, U(s;x), which are analytic in the strip, 0≤Res≤1, Ims>T 0 , and in addition have an asymptotic expansion in powers of [s(s-1)] -1 , i.e. U(s;x)=V 0 (x)+gV 1 (x)+g 2 V 2 (x)+...+O(g N ), with g=[s(s-1)] -1 . The potentials V n (x) are real and summable. Under suitable conditions on the V n 's and the O(g N ) term we show that the condition, ∫ 0 ∞ vertical bar f 0 (x) vertical bar 2 V 1 (x) dx≠0, where f 0 is the zero energy and g=0 Jost function for U, is sufficient to guarantee that the zeros g n are real and, hence, s n =1/2+iγ n , for γ n ≥T 0 .Starting with a judiciously chosen Jost function, M(s,k), which is constructed such that M(s,0) is Riemann's ξ(s) function, we have used inverse scattering methods to actually construct a U(s;x) with the above properties. By necessity, we had to generalize inverse methods to deal with complex potentials and a nonunitary S-matrix. This we have done at least for the special cases under consideration.For our specific example, ∫ 0 ∞ vertical bar f 0 (x) vertical bar 2 V 1 (x) dx=0 and, hence, we get no restriction on Img n or Res n . The reasons for the vanishing of the above integral are given, and they give us hints on what one needs to proceed further. The problem

Physics-based Inverse Problem to Deduce Marine Atmospheric Boundary Layer Parameters

Science.gov (United States)

2017-03-07

knowledge and capabilities in the use and development of inverse problem techniques to deduce atmospheric parameters. WORK COMPLETED The research completed...please find the Final Technical Report with SF 298 for Dr. Erin E. Hackett’s ONR grant entitled Physics -based Inverse Problem to Deduce Marine...From- To) 07/03/2017 Final Technica l Dec 2012- Dec 2016 4. TITLE AND SUBTITLE 5a. CONTRACT NUMBER Physics -based Inverse Problem to Deduce Marine

Energy-dependent inversion of p+16O scattering data

International Nuclear Information System (INIS)

Cooper, S.G.

1997-01-01

A fast iterative procedure is developed to determine potentials by inversion from elastic cross section, analysing powers and reaction cross-section measurements covering a wide energy range. The procedure incorporates both energy and parity dependence. The method is applied to extensive p+ 16 O scattering data for an energy range from 27.3 to 46.1 MeV, giving a solution which simultaneously reproduces the data at all energies. The wide angle data is well reproduced by including parity dependence and a linear energy dependence is established for the real potential, including the parity-dependent component. The real terms agree qualitatively with potentials derived from the single channel RGM, but the central and spin-orbit imaginary components have distinct features strongly suggestive of further non-local contributions, possibly arising from channel coupling. The large data set is found essential to reduce the potential ambiguities present when fitting scattering data. (orig.)

Pulsar high energy emission due to inverse Compton scattering

Energy Technology Data Exchange (ETDEWEB)

Lyutikov, Maxim

2013-06-15

We discuss growing evidence that pulsar high energy is emission is generated via Inverse Compton mechanism. We reproduce the broadband spectrum of Crab pulsar, from UV to very high energy gamma-rays - nearly ten decades in energy, within the framework of the cyclotron-self-Compton model. Emission is produced by two counter-streaming beams within the outer gaps, at distances above ∼ 20 NS radii. The outward moving beam produces UV-X-ray photons via Doppler-booster cyclotron emission, and GeV photons by Compton scattering the cyclotron photons produced by the inward going beam. The scattering occurs in the deep Klein-Nishina regime, whereby the IC component provides a direct measurement of particle distribution within the magnetosphere. The required plasma multiplicity is high, ∼10{sup 6} – 10{sup 7}, but is consistent with the average particle flux injected into the pulsar wind nebula.

An Entropic Estimator for Linear Inverse Problems

Directory of Open Access Journals (Sweden)

Amos Golan

2012-05-01

Full Text Available In this paper we examine an Information-Theoretic method for solving noisy linear inverse estimation problems which encompasses under a single framework a whole class of estimation methods. Under this framework, the prior information about the unknown parameters (when such information exists, and constraints on the parameters can be incorporated in the statement of the problem. The method builds on the basics of the maximum entropy principle and consists of transforming the original problem into an estimation of a probability density on an appropriate space naturally associated with the statement of the problem. This estimation method is generic in the sense that it provides a framework for analyzing non-normal models, it is easy to implement and is suitable for all types of inverse problems such as small and or ill-conditioned, noisy data. First order approximation, large sample properties and convergence in distribution are developed as well. Analytical examples, statistics for model comparisons and evaluations, that are inherent to this method, are discussed and complemented with explicit examples.

Geostatistical regularization operators for geophysical inverse problems on irregular meshes

Science.gov (United States)

Jordi, C.; Doetsch, J.; Günther, T.; Schmelzbach, C.; Robertsson, J. OA

2018-05-01

Irregular meshes allow to include complicated subsurface structures into geophysical modelling and inverse problems. The non-uniqueness of these inverse problems requires appropriate regularization that can incorporate a priori information. However, defining regularization operators for irregular discretizations is not trivial. Different schemes for calculating smoothness operators on irregular meshes have been proposed. In contrast to classical regularization constraints that are only defined using the nearest neighbours of a cell, geostatistical operators include a larger neighbourhood around a particular cell. A correlation model defines the extent of the neighbourhood and allows to incorporate information about geological structures. We propose an approach to calculate geostatistical operators for inverse problems on irregular meshes by eigendecomposition of a covariance matrix that contains the a priori geological information. Using our approach, the calculation of the operator matrix becomes tractable for 3-D inverse problems on irregular meshes. We tested the performance of the geostatistical regularization operators and compared them against the results of anisotropic smoothing in inversions of 2-D surface synthetic electrical resistivity tomography (ERT) data as well as in the inversion of a realistic 3-D cross-well synthetic ERT scenario. The inversions of 2-D ERT and seismic traveltime field data with geostatistical regularization provide results that are in good accordance with the expected geology and thus facilitate their interpretation. In particular, for layered structures the geostatistical regularization provides geologically more plausible results compared to the anisotropic smoothness constraints.

Carleman estimates and applications to inverse problems for hyperbolic systems

CERN Document Server

Bellassoued, Mourad

2017-01-01

This book is a self-contained account of the method based on Carleman estimates for inverse problems of determining spatially varying functions of differential equations of the hyperbolic type by non-overdetermining data of solutions. The formulation is different from that of Dirichlet-to-Neumann maps and can often prove the global uniqueness and Lipschitz stability even with a single measurement. These types of inverse problems include coefficient inverse problems of determining physical parameters in inhomogeneous media that appear in many applications related to electromagnetism, elasticity, and related phenomena. Although the methodology was created in 1981 by Bukhgeim and Klibanov, its comprehensive development has been accomplished only recently. In spite of the wide applicability of the method, there are few monographs focusing on combined accounts of Carleman estimates and applications to inverse problems. The aim in this book is to fill that gap. The basic tool is Carleman estimates, the theory of wh...

Inverse periodic problem for the discrete approximation of the Schroedinger nonlinear equation

International Nuclear Information System (INIS)

Bogolyubov, N.N.; Prikarpatskij, A.K.; AN Ukrainskoj SSR, Lvov. Inst. Prikladnykh Problem Mekhaniki i Matematiki)

1982-01-01

The problem of numerical solution of the Schroedinger nonlinear equation (1) iPSIsub(t) = PSIsub(xx)+-2(PSI)sup(2)PSI. The numerical solution of nonlinear differential equation supposes its discrete approximation is required for the realization of the computer calculation process. Tor the equation (1) there exists the following discrete approximation by variable x(2) iPSIsub(n, t) = (PSIsub(n+1)-2PSIsub(n)+PSIsub(n-1))/(Δx)sup(2)+-(PSIsub(n))sup(2)(PSIsub(n+1)+PSIsub(n-1)), n=0, +-1, +-2... where PSIsub(n)(+) is the corresponding value of PSI(x, t) function in the node and divisions with the equilibrium step Δx. The main problem is obtaining analytically exact solutions of the equations (2). The analysis of the equation system (2) is performed on the base of the discrete analogue of the periodic variant of the inverse scattering problem method developed with the aid of nonlinear equations of the Korteweg-de Vries type. Obtained in explicit form are analytical solutions of the equations system (2). The solutions are expressed through the Riemann THETA-function [ru

Transient electromagnetic scattering on anisotropic media

International Nuclear Information System (INIS)

Stewart, R.D.

1990-01-01

This dissertation treats the problem of transient scattering of obliquely incident electromagnetic plane waves on a stratified anisotropic dielectric slab. Scattering operators are derived for the reflective response of the medium. The internal fields are calculated. Wave splitting and invariant imbedding techniques are used. These techniques are first presented for fields normally incident on a stratified, isotropic dielectric medium. The techniques of wave splitting and invariant imbedding are applied to normally incident plane waves on an anisotropic medium. An integro-differential equation is derived for the reflective response and the direct and inverse scattering problems are discussed. These techniques are applied to the case of obliquely incident plane waves. The reflective response is derived and the direct and inverse problems discussed and compared to those for the normal incidence case. The internal fields are investigated for the oblique incidence via a Green's function approach. A numerical scheme is presented to calculate the Green's function. Finally, symmetry relations of the reflective response are discussed

Direct and inverse source problems for a space fractional advection dispersion equation

KAUST Repository

Aldoghaither, Abeer

2016-05-15

In this paper, direct and inverse problems for a space fractional advection dispersion equation on a finite domain are studied. The inverse problem consists in determining the source term from final observations. We first derive the analytic solution to the direct problem which we use to prove the uniqueness and the unstability of the inverse source problem using final measurements. Finally, we illustrate the results with a numerical example.

Bayesian probability theory and inverse problems

International Nuclear Information System (INIS)

Kopec, S.

1994-01-01

Bayesian probability theory is applied to approximate solving of the inverse problems. In order to solve the moment problem with the noisy data, the entropic prior is used. The expressions for the solution and its error bounds are presented. When the noise level tends to zero, the Bayesian solution tends to the classic maximum entropy solution in the L 2 norm. The way of using spline prior is also shown. (author)

Inverse source problems for eddy current equations

International Nuclear Information System (INIS)

Rodríguez, Ana Alonso; Valli, Alberto; Camaño, Jessika

2012-01-01

We study the inverse source problem for the eddy current approximation of Maxwell equations. As for the full system of Maxwell equations, we show that a volume current source cannot be uniquely identified by knowledge of the tangential components of the electromagnetic fields on the boundary, and we characterize the space of non-radiating sources. On the other hand, we prove that the inverse source problem has a unique solution if the source is supported on the boundary of a subdomain or if it is the sum of a finite number of dipoles. We address the applicability of this result for the localization of brain activity from electroencephalography and magnetoencephalography measurements. (paper)

Moebius inverse problem for distorted black holes

International Nuclear Information System (INIS)

Rosu, H.

1993-01-01

Hawking ''thermal'' radiation could be a means to detect black holes of micron sizes, which may be hovering through the universe. We consider these micro-black holes to be distorted by the presence of some distribution of matter representing a convolution factor for their Hawking radiation. One may hope to determine from their Hawking signals the temperature distribution of their material shells by the inverse black body problem. In 1990, Nan-xian Chen has used a so-called modified Moebius transform to solve the inverse black body problem. We discuss and apply this technique to Hawking radiation. Some comments on supersymmetric applications of Moebius function and transform are also added. (author). 22 refs

Understanding of phase modulation in two-level systems through inverse scattering

International Nuclear Information System (INIS)

Hasenfeld, A.; Hammes, S.L.; Warren, W.S.

1988-01-01

Analytical and numerical calculations describe the effects of shaped radiation pulses on two-level systems in terms of quantum-mechanical scattering. Previous results obtained in the reduced case of amplitude modulation are extended to the general case of simultaneous amplitude and phase modulation. We show that an infinite family of phase- and amplitude-modulated pulses all generate rectangular inversion profiles. Experimental measurements also verify the theoretical analysis

Lax-pair operators for squared eigenfunctions in the inverse scattering transformations

International Nuclear Information System (INIS)

Iino, Kazuhiro; Ichikawa, Yoshihiko.

1982-05-01

Modification of the algorithm of Chen, Lee and Liu enables us to construct alternative Lax-pair operators for the Korteweg-de Vries equation and the modified Korteweg-de Vries equation. These Lax-pair operators stand as the Lax-pair operators for the squared eigenfunction and the sum of squared eigenfunctions of the Ablowitz-Kaup-Newell-Segur inverse scattering transformation for these celebrated nonlinear evolution equations. (author)

Collage-based approaches for elliptic partial differential equations inverse problems

Science.gov (United States)

Yodzis, Michael; Kunze, Herb

2017-01-01

The collage method for inverse problems has become well-established in the literature in recent years. Initial work developed a collage theorem, based upon Banach's fixed point theorem, for treating inverse problems for ordinary differential equations (ODEs). Amongst the subsequent work was a generalized collage theorem, based upon the Lax-Milgram representation theorem, useful for treating inverse problems for elliptic partial differential equations (PDEs). Each of these two different approaches can be applied to elliptic PDEs in one space dimension. In this paper, we explore and compare how the two different approaches perform for the estimation of the diffusivity for a steady-state heat equation.

General inverse problems for regular variation

DEFF Research Database (Denmark)

Damek, Ewa; Mikosch, Thomas Valentin; Rosinski, Jan

2014-01-01

Regular variation of distributional tails is known to be preserved by various linear transformations of some random structures. An inverse problem for regular variation aims at understanding whether the regular variation of a transformed random object is caused by regular variation of components ...

On the quantum inverse problem for a new type of nonlinear Schroedinger equation for Alfven waves in plasma

International Nuclear Information System (INIS)

Sen, S.; Roy Chowdhury, A.

1989-06-01

The nonlinear Alfven waves are governed by the Vector Derivative nonlinear Schroedinger (VDNLS) equation, which for parallel or quasi parallel propagation reduces to the Derivative Nonlinear Schroedinger (DNLS) equation for the circularly polarized waves. We have formulated the Quantum Inverse problem for a new type of Nonlinear Schroedinger Equation which has many properties similar to the usual NLS problem but the structure of classical and quantum R matrix are distinctly different. The commutation rules of the scattering data are obtained and the Algebraic Bethe Ansatz is formulated to derive the eigenvalue equation for the energy of the excited states. 10 refs

Inverse problem of solar oscillations

International Nuclear Information System (INIS)

Sekii, T.; Shibahashi, H.

1987-01-01

The authors present some preliminary results of numerical simulation to infer the sound velocity distribution in the solar interior from the oscillation data of the Sun as the inverse problem. They analyze the acoustic potential itself by taking account of some factors other than the sound velocity, and infer the sound velocity distribution in the deep interior of the Sun

Inverse scattering with mixed spectrum from δ-potentials

International Nuclear Information System (INIS)

Lin Jiancheng.

1987-03-01

The inverse problem is studied in a system with mixed spectrum, i.e. the continuous part of the spectrum coincides with that of a repulsive δ-potential and the discrete part coincides with that of an attractive δ-potential. (author). 2 refs, 5 figs

Wave scattering theory and the absorption problem for a black hole

International Nuclear Information System (INIS)

Sanchez, N.

1977-01-01

The general problem of scattering and absorption of waves from a Schwarzschild black hole is investigated. A scattering absorption amplitude is introduced. The unitarity theorem for this problem is derived from the wave equation and its boundary conditions. The formulation of the problem, within the formal scattering theory approach, is also given. The existence of a singularity in space-time is related explicitly to the presence of a nonzero absorption cross section. Another derivation of the unitarity theorem for our problem is given by operator methods. The reciprocity relation is also proved; that is, for the scattering of waves the black hole is a reciprocal system. Finally, the elastic scattering problem is considered, and the elastic scattering amplitude is calculated for high frequencies and small scattering angles

Efficient scattering-angle enrichment for a nonlinear inversion of the background and perturbations components of a velocity model

KAUST Repository

Wu, Zedong; Alkhalifah, Tariq Ali

2017-01-01

Reflection-waveform inversion (RWI) can help us reduce the nonlinearity of the standard full-waveform inversion (FWI) by inverting for the background velocity model using the wave-path of a single scattered wavefield to an image. However, current

Direct and Inverse Problems in Statistical Wavefields

International Nuclear Information System (INIS)

Wolf, Emil

2002-01-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

Efficient generalized Golub-Kahan based methods for dynamic inverse problems

Science.gov (United States)

Chung, Julianne; Saibaba, Arvind K.; Brown, Matthew; Westman, Erik

2018-02-01

We consider efficient methods for computing solutions to and estimating uncertainties in dynamic inverse problems, where the parameters of interest may change during the measurement procedure. Compared to static inverse problems, incorporating prior information in both space and time in a Bayesian framework can become computationally intensive, in part, due to the large number of unknown parameters. In these problems, explicit computation of the square root and/or inverse of the prior covariance matrix is not possible, so we consider efficient, iterative, matrix-free methods based on the generalized Golub-Kahan bidiagonalization that allow automatic regularization parameter and variance estimation. We demonstrate that these methods for dynamic inversion can be more flexible than standard methods and develop efficient implementations that can exploit structure in the prior, as well as possible structure in the forward model. Numerical examples from photoacoustic tomography, space-time deblurring, and passive seismic tomography demonstrate the range of applicability and effectiveness of the described approaches. Specifically, in passive seismic tomography, we demonstrate our approach on both synthetic and real data. To demonstrate the scalability of our algorithm, we solve a dynamic inverse problem with approximately 43 000 measurements and 7.8 million unknowns in under 40 s on a standard desktop.

Compact FEL-driven inverse compton scattering gamma-ray source

Energy Technology Data Exchange (ETDEWEB)

Placidi, M. [Lawrence Berkeley National Laboratory, Berkeley, CA 94720 (United States); Di Mitri, S., E-mail: simone.dimitri@elettra.eu [Elettra - Sincrotrone Trieste S.C.p.A., 34149 Basovizza, Trieste (Italy); Pellegrini, C. [SLAC National Accelerator Laboratory, Menlo Park, CA 94025 (United States); University of California, Los Angeles, CA 90095 (United States); Penn, G. [Lawrence Berkeley National Laboratory, Berkeley, CA 94720 (United States)

2017-05-21

Many research and applications areas require photon sources capable of producing gamma-ray beams in the multi-MeV energy range with reasonably high fluxes and compact footprints. Besides industrial, nuclear physics and security applications, a considerable interest comes from the possibility to assess the state of conservation of cultural assets like statues, columns etc., via visualization and analysis techniques using high energy photon beams. Computed Tomography scans, widely adopted in medicine at lower photon energies, presently provide high quality three-dimensional imaging in industry and museums. We explore the feasibility of a compact source of quasi-monochromatic, multi-MeV gamma-rays based on Inverse Compton Scattering (ICS) from a high intensity ultra-violet (UV) beam generated in a free-electron laser by the electron beam itself. This scheme introduces a stronger relationship between the energy of the scattered photons and that of the electron beam, resulting in a device much more compact than a classic ICS for a given scattered energy. The same electron beam is used to produce gamma-rays in the 10–20 MeV range and UV radiation in the 10–15 eV range, in a ~4×22 m{sup 2} footprint system.

Maximum a posteriori probability estimates in infinite-dimensional Bayesian inverse problems

International Nuclear Information System (INIS)

Helin, T; Burger, M

2015-01-01

A demanding challenge in Bayesian inversion is to efficiently characterize the posterior distribution. This task is problematic especially in high-dimensional non-Gaussian problems, where the structure of the posterior can be very chaotic and difficult to analyse. Current inverse problem literature often approaches the problem by considering suitable point estimators for the task. Typically the choice is made between the maximum a posteriori (MAP) or the conditional mean (CM) estimate. The benefits of either choice are not well-understood from the perspective of infinite-dimensional theory. Most importantly, there exists no general scheme regarding how to connect the topological description of a MAP estimate to a variational problem. The recent results by Dashti and others (Dashti et al 2013 Inverse Problems 29 095017) resolve this issue for nonlinear inverse problems in Gaussian framework. In this work we improve the current understanding by introducing a novel concept called the weak MAP (wMAP) estimate. We show that any MAP estimate in the sense of Dashti et al (2013 Inverse Problems 29 095017) is a wMAP estimate and, moreover, how the wMAP estimate connects to a variational formulation in general infinite-dimensional non-Gaussian problems. The variational formulation enables to study many properties of the infinite-dimensional MAP estimate that were earlier impossible to study. In a recent work by the authors (Burger and Lucka 2014 Maximum a posteriori estimates in linear inverse problems with logconcave priors are proper bayes estimators preprint) the MAP estimator was studied in the context of the Bayes cost method. Using Bregman distances, proper convex Bayes cost functions were introduced for which the MAP estimator is the Bayes estimator. Here, we generalize these results to the infinite-dimensional setting. Moreover, we discuss the implications of our results for some examples of prior models such as the Besov prior and hierarchical prior. (paper)

Obtaining the crystal potential by inversion from electron scattering intensities

International Nuclear Information System (INIS)

Allen, L.T.; Josefsson, T.W.; Leeb, H.

1998-01-01

A method to obtain the crystal potential from the intensities of the diffracted beams in high energy electron diffraction is proposed. It is based on a series of measurements for specific well determined orientations of the incident beam which determine the moduli of all elements of the scattering matrix. Using unitarity and the specific form of the scattering matrix (including symmetries) an overdetermined set of non-linear equations is obtained from these data. Solution of these equations yields the required phase information and allows the determination of a (projected) crystal potential by inversion which is unique up to an arbitrary shift of the origin. The reconstruction of potentials from intensities is illustrated for two realistic examples, a [111] systematic row case in ZnS and a [110] zone axis orientation in GaAs (both noncentrosymmetric crystals)

Wigner representation in scattering problems

International Nuclear Information System (INIS)

Remler, E.A.

1975-01-01

The basic equations of quantum scattering are translated into the Wigner representation. This puts quantum mechanics in the form of a stochastic process in phase space. Instead of complex valued wavefunctions and transition matrices, one now works with real-valued probability distributions and source functions, objects more responsive to physical intuition. Aside from writing out certain necessary basic expressions, the main purpose is to develop and stress the interpretive picture associated with this representation and to derive results used in applications published elsewhere. The quasiclassical guise assumed by the formalism lends itself particularly to approximations of complex multiparticle scattering problems is laid. The foundation for a systematic application of statistical approximations to such problems. The form of the integral equation for scattering as well as its mulitple scattering expansion in this representation are derived. Since this formalism remains unchanged upon taking the classical limit, these results also constitute a general treatment of classical multiparticle collision theory. Quantum corrections to classical propogators are discussed briefly. The basic approximation used in the Monte Carlo method is derived in a fashion that allows for future refinement and includes bound state production. The close connection that must exist between inclusive production of a bound state and of its constituents is brought out in an especially graphic way by this formalism. In particular one can see how comparisons between such cross sections yield direct physical insight into relevant production mechanisms. A simple illustration of scattering by a bound two-body system is treated. Simple expressions for single- and double-scattering contributions to total and differential cross sections, as well as for all necessary shadow corrections thereto, are obtained and compared to previous results of Glauber and Goldberger

Method and software to solution of inverse and inverse design fluid flow and heat transfer problems is compatible with CFD-software

Energy Technology Data Exchange (ETDEWEB)

Krukovsky, P G [Institute of Engineering Thermophysics, National Academy of Sciences of Ukraine, Kiev (Ukraine)

1998-12-31

The description of method and software FRIEND which provide a possibility of solution of inverse and inverse design problems on the basis of existing (base) CFD-software for solution of direct problems (in particular, heat-transfer and fluid-flow problems using software PHOENICS) are presented. FRIEND is an independent additional module that widens the operational capacities of the base software unified with this module. This unifying does not require any change or addition to the base software. Interfacing of FRIEND and the base software takes place through input and output files of the base software. A brief description of the computational technique applied for the inverse problem solution, same detailed information on the interfacing of FRIEND and CFD-software and solution results for testing inverse and inverse design problems, obtained using the tandem CFD-software PHOENICS and FRIEND, are presented. (author) 9 refs.

Method and software to solution of inverse and inverse design fluid flow and heat transfer problems is compatible with CFD-software

Energy Technology Data Exchange (ETDEWEB)

Krukovsky, P.G. [Institute of Engineering Thermophysics, National Academy of Sciences of Ukraine, Kiev (Ukraine)

1997-12-31

The description of method and software FRIEND which provide a possibility of solution of inverse and inverse design problems on the basis of existing (base) CFD-software for solution of direct problems (in particular, heat-transfer and fluid-flow problems using software PHOENICS) are presented. FRIEND is an independent additional module that widens the operational capacities of the base software unified with this module. This unifying does not require any change or addition to the base software. Interfacing of FRIEND and the base software takes place through input and output files of the base software. A brief description of the computational technique applied for the inverse problem solution, same detailed information on the interfacing of FRIEND and CFD-software and solution results for testing inverse and inverse design problems, obtained using the tandem CFD-software PHOENICS and FRIEND, are presented. (author) 9 refs.

Analog fault diagnosis by inverse problem technique

KAUST Repository

Ahmed, Rania F.

2011-12-01

A novel algorithm for detecting soft faults in linear analog circuits based on the inverse problem concept is proposed. The proposed approach utilizes optimization techniques with the aid of sensitivity analysis. The main contribution of this work is to apply the inverse problem technique to estimate the actual parameter values of the tested circuit and so, to detect and diagnose single fault in analog circuits. The validation of the algorithm is illustrated through applying it to Sallen-Key second order band pass filter and the results show that the detecting percentage efficiency was 100% and also, the maximum error percentage of estimating the parameter values is 0.7%. This technique can be applied to any other linear circuit and it also can be extended to be applied to non-linear circuits. © 2011 IEEE.

Time-reversal of electromagnetic scattering for small scatterer classification

International Nuclear Information System (INIS)

Smith, J Torquil; Berryman, James G

2012-01-01

Time-reversal operators, or the alternatively labelled, but equivalent, multistatic response matrix methods, are used to show how to determine the number of scatterers present in an electromagnetic scattering scenario that might be typical of UneXploded Ordinance (UXO) detection, classification and removal applications. Because the nature of the target UXO application differs from that of many other common inversion problems, emphasis is placed here on classification and enumeration rather than on detailed imaging. The main technical issues necessarily revolve around showing that it is possible to find a sufficient number of constraints via multiple measurements (i.e. using several distinct views at the target site) to solve the enumeration problem. The main results show that five measurements with antenna pairs are generally adequate to solve the classification and enumeration problems. However, these results also demonstrate a need for decreasing noise levels in the multistatic matrix as the number n of scatterers increases for the intended practical applications of the method. (paper)

Comparing Two Approaches for Point-Like Scatterer Detection

Directory of Open Access Journals (Sweden)

Angela Dell’Aversano

2015-01-01

Full Text Available Many inverse scattering problems concern the detection and localisation of point-like scatterers which are sparsely enclosed within a prescribed investigation domain. Therefore, it looks like a good option to tackle the problem by applying reconstruction methods that are properly tailored for such a type of scatterers or that naturally enforce sparsity in the reconstructions. Accordingly, in this paper we compare the time reversal-MUSIC and the compressed sensing. The study develops through numerical examples and focuses on the role of noise in data and mutual coupling between the scatterers.

Eigenvectors phase correction in inverse modal problem

Science.gov (United States)

Qiao, Guandong; Rahmatalla, Salam

2017-12-01

The solution of the inverse modal problem for the spatial parameters of mechanical and structural systems is heavily dependent on the quality of the modal parameters obtained from the experiments. While experimental and environmental noises will always exist during modal testing, the resulting modal parameters are expected to be corrupted with different levels of noise. A novel methodology is presented in this work to mitigate the errors in the eigenvectors when solving the inverse modal problem for the spatial parameters. The phases of the eigenvector component were utilized as design variables within an optimization problem that minimizes the difference between the calculated and experimental transfer functions. The equation of motion in terms of the modal and spatial parameters was used as a constraint in the optimization problem. Constraints that reserve the positive and semi-positive definiteness and the inter-connectivity of the spatial matrices were implemented using semi-definite programming. Numerical examples utilizing noisy eigenvectors with augmented Gaussian white noise of 1%, 5%, and 10% were used to demonstrate the efficacy of the proposed method. The results showed that the proposed method is superior when compared with a known method in the literature.

Inverse electronic scattering by Green's functions and singular values decomposition

International Nuclear Information System (INIS)

Mayer, A.; Vigneron, J.-P.

2000-01-01

An inverse scattering technique is developed to enable a sample reconstruction from the diffraction figures obtained by electronic projection microscopy. In its Green's functions formulation, this technique takes account of all orders of diffraction by performing an iterative reconstruction of the wave function on the observation screen. This scattered wave function is then backpropagated to the sample to determine the potential-energy distribution, which is assumed real valued. The method relies on the use of singular values decomposition techniques, thus providing the best least-squares solutions and enabling a reduction of noise. The technique is applied to the analysis of a two-dimensional nanometric sample that is observed in Fresnel conditions with an electronic energy of 25 eV. The algorithm turns out to provide results with a mean relative error of the order of 5% and to be very stable against random noise

Bilinear Inverse Problems: Theory, Algorithms, and Applications

Science.gov (United States)

Ling, Shuyang

We will discuss how several important real-world signal processing problems, such as self-calibration and blind deconvolution, can be modeled as bilinear inverse problems and solved by convex and nonconvex optimization approaches. In Chapter 2, we bring together three seemingly unrelated concepts, self-calibration, compressive sensing and biconvex optimization. We show how several self-calibration problems can be treated efficiently within the framework of biconvex compressive sensing via a new method called SparseLift. More specifically, we consider a linear system of equations y = DAx, where the diagonal matrix D (which models the calibration error) is unknown and x is an unknown sparse signal. By "lifting" this biconvex inverse problem and exploiting sparsity in this model, we derive explicit theoretical guarantees under which both x and D can be recovered exactly, robustly, and numerically efficiently. In Chapter 3, we study the question of the joint blind deconvolution and blind demixing, i.e., extracting a sequence of functions [special characters omitted] from observing only the sum of their convolutions [special characters omitted]. In particular, for the special case s = 1, it becomes the well-known blind deconvolution problem. We present a non-convex algorithm which guarantees exact recovery under conditions that are competitive with convex optimization methods, with the additional advantage of being computationally much more efficient. We discuss several applications of the proposed framework in image processing and wireless communications in connection with the Internet-of-Things. In Chapter 4, we consider three different self-calibration models of practical relevance. We show how their corresponding bilinear inverse problems can be solved by both the simple linear least squares approach and the SVD-based approach. As a consequence, the proposed algorithms are numerically extremely efficient, thus allowing for real-time deployment. Explicit theoretical

Inverse Compton gamma-rays from pulsars

International Nuclear Information System (INIS)

Morini, M.

1983-01-01

A model is proposed for pulsar optical and gamma-ray emission where relativistic electrons beams: (i) scatter the blackbody photons from the polar cap surface giving inverse Compton gamma-rays and (ii) produce synchrotron optical photons in the light cylinder region which are then inverse Compton scattered giving other gamma-rays. The model is applied to the Vela pulsar, explaining the first gamma-ray pulse by inverse Compton scattering of synchrotron photons near the light cylinder and the second gamma-ray pulse partly by inverse Compton scattering of synchrotron photons and partly by inverse Compton scattering of the thermal blackbody photons near the star surface. (author)

Forward and inverse viscoelastic wave scattering by irregular inclusions for shear wave elastography.

Science.gov (United States)

Bernard, Simon; Cloutier, Guy

2017-10-01

Inversion methods in shear wave elastography use simplifying assumptions to recover the mechanical properties of soft tissues. Consequently, these methods suffer from artifacts when applied to media containing strong stiffness contrasts, and do not provide a map of the viscosity. In this work, the shear wave field recorded inside and around an inclusion was used to estimate the viscoelastic properties of the inclusion and surrounding medium, based on an inverse problem approach assuming local homogeneity of both media. An efficient semi-analytical method was developed to model the scattering of an elastic wave by an irregular inclusion, based on a decomposition of the field by Bessel functions and on a decomposition of the boundaries as Fourier series. This model was validated against finite element modeling. Shear waves were experimentally induced by acoustic radiation force in soft tissue phantoms containing stiff and soft inclusions, and the displacement field was imaged at a high frame rate using plane wave imaging. A nonlinear least-squares algorithm compared the model to the experimental data and adjusted the geometrical and mechanical parameters. The estimated shear storage and loss moduli were in good agreement with reference measurements, as well as the estimated inclusion shape. This approach provides an accurate estimation of geometry and viscoelastic properties for a single inclusion in a homogeneous background in the context of radiation force elastography.

Inverse Problems in Systems Biology: A Critical Review.

Science.gov (United States)

Guzzi, Rodolfo; Colombo, Teresa; Paci, Paola

2018-01-01

Systems Biology may be assimilated to a symbiotic cyclic interplaying between the forward and inverse problems. Computational models need to be continuously refined through experiments and in turn they help us to make limited experimental resources more efficient. Every time one does an experiment we know that there will be some noise that can disrupt our measurements. Despite the noise certainly is a problem, the inverse problems already involve the inference of missing information, even if the data is entirely reliable. So the addition of a certain limited noise does not fundamentally change the situation but can be used to solve the so-called ill-posed problem, as defined by Hadamard. It can be seen as an extra source of information. Recent studies have shown that complex systems, among others the systems biology, are poorly constrained and ill-conditioned because it is difficult to use experimental data to fully estimate their parameters. For these reasons was born the concept of sloppy models, a sequence of models of increasing complexity that become sloppy in the limit of microscopic accuracy. Furthermore the concept of sloppy models contains also the concept of un-identifiability, because the models are characterized by many parameters that are poorly constrained by experimental data. Then a strategy needs to be designed to infer, analyze, and understand biological systems. The aim of this work is to provide a critical review to the inverse problems in systems biology defining a strategy to determine the minimal set of information needed to overcome the problems arising from dynamic biological models that generally may have many unknown, non-measurable parameters.

Direct and inverse source problems for a space fractional advection dispersion equation

KAUST Repository

Aldoghaither, Abeer; Laleg-Kirati, Taous-Meriem; Liu, Da Yan

2016-01-01

In this paper, direct and inverse problems for a space fractional advection dispersion equation on a finite domain are studied. The inverse problem consists in determining the source term from final observations. We first derive the analytic

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

Science.gov (United States)

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.

An inverse Sturm–Liouville problem with a fractional derivative

KAUST Repository

Jin, Bangti; Rundell, William

2012-01-01

In this paper, we numerically investigate an inverse problem of recovering the potential term in a fractional Sturm-Liouville problem from one spectrum. The qualitative behaviors of the eigenvalues and eigenfunctions are discussed, and numerical

ISOGEOMETRIC SHAPE OPTIMIZATION FOR ELECTROMAGNETIC SCATTERING PROBLEMS

DEFF Research Database (Denmark)

Nguyen, D. M.; Evgrafov, Anton; Gravesen, Jens

2012-01-01

We consider the benchmark problem of magnetic energy density enhancement in a small spatial region by varying the shape of two symmetric conducting scatterers. We view this problem as a prototype for a wide variety of geometric design problems in electromagnetic applications. Our approach...

An inverse problem for a one-dimensional time-fractional diffusion problem

KAUST Repository

Jin, Bangti; Rundell, William

2012-01-01

We study an inverse problem of recovering a spatially varying potential term in a one-dimensional time-fractional diffusion equation from the flux measurements taken at a single fixed time corresponding to a given set of input sources. The unique

New approaches to old scattering problems

Science.gov (United States)

Goldberg, Joshua David

This thesis is broken into two parts corresponding to the research done before and after the death of Roger Dashen. The first part addresses the problem of high frequency scattering from flat objects. A new formalism, developed by Dashen and Wurmser, is applied to the two dimensional problem of scattering from a soft infinite strip. It is seen that the cross section can be related to a quantity, termed the divergence coefficient, which describes the behavior of the field near the edges. A simple expression for the divergence coefficient and the scattering cross section is derived which, in contrast to traditional results, is uniformly valid in the high frequency limit. The basic procedure is to first express the divergence coefficient in a series form involving Mathieu functions, approximating the terms in this series by their uniform WKB representation, and using the Poisson sum formula to convert the WKB-based series to a more rapidly converging series of integrals which can be evaluated asymptotically. The result is a new expression for the scattering cross section which is compared with previously obtained results. The second part addresses a specific problem in the field of wave propagation in random media: computing the average field for the case of a plane wave incident on a region with a weakly fluctuating sound speed. A review of the existing mathematical methods for treating this problem in both the small and large-scale fluctuation cases is given. In the small-scale regime, previously unrecognized problems with the closure theory are discussed and numerical results are given which illustrate the role played by backscattering in this type of propagation. In the large-scale regime, a new mathematical approach, analogous to the renormalization technique, is described and used to derive a new expression for the mean field valid in this limit. This result is compared with the traditional expressions for this quantity.

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

Science.gov (United States)

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

Separable expansion for realistic multichannel scattering problems

International Nuclear Information System (INIS)

Canton, L.; Cattapan, G.; Pisent, G.

1987-01-01

A new approach to the multichannel scattering problem with realistic local or nonlocal interactions is developed. By employing the negative-energy solutions of uncoupled Sturmian eigenvalue problems referring to simple auxiliary potentials, the coupling interactions appearing to the original multichannel problem are approximated by finite-rank potentials. By resorting to integral-equation tecniques the coupled-channel equations are then reduced to linear algebraic equations which can be straightforwardly solved. Compact algebraic expressions for the relevant scattering matrix elements are thus obtained. The convergence of the method is tasted in the single-channel case with realistic optical potentials. Excellent agreement is obtained with a few terms in the separable expansion for both real and absorptive interactions

Numerical investigation of the inverse blackbody radiation problem

International Nuclear Information System (INIS)

Xin Tan, Guo-zhen Yang, Ben-yuan Gu

1994-01-01

A numerical algorithm for the inverse blackbody radiation problem, which is the determination of the temperature distribution of a thermal radiator (TDTR) from its total radiated power spectrum (TRPS), is presented, based on the general theory of amplitude-phase retrieval. With application of this new algorithm, the ill-posed nature of the Fredholm equation of the first kind can be largely overcome and a convergent solution to high accuracy can be obtained. By incorporation of the hybrid input-output algorithm into our algorithm, the convergent process can be substantially expedited and the stagnation problem of the solution can be averted. From model calculations it is found that the new algorithm can also provide a robust reconstruction of the TDTR from the noise-corrupted data of the TRPS. Therefore the new algorithm may offer a useful approach to solving the ill-posed inverse problem. 18 refs., 9 figs

Variational structure of inverse problems in wave propagation and vibration

Energy Technology Data Exchange (ETDEWEB)

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.

Inversion of real and complex phase shifts to potentials by the generalized Cox-Thompson inverse scattering method at fixed energy

International Nuclear Information System (INIS)

Melchert, O; Scheid, W; Apagyi, B

2006-01-01

The Cox-Thompson inverse scattering method at fixed energy has been generalized to treat complex phase shifts derived from experiments. New formulae for relating phase shifts to shifted angular momenta are derived. The method is applied to phase shifts of known potentials in order to test its quality and stability and, further, it is used to invert experimental n-α and n- 12 C phase shifts

FOREWORD: 2nd International Workshop on New Computational Methods for Inverse Problems (NCMIP 2012)

Science.gov (United States)

Blanc-Féraud, Laure; Joubert, Pierre-Yves

2012-09-01

Conference logo This volume of Journal of Physics: Conference Series is dedicated to the scientific contributions presented during the 2nd International Workshop on New Computational Methods for Inverse Problems, (NCMIP 2012). This workshop took place at Ecole Normale Supérieure de Cachan, in Cachan, France, on 15 May 2012, at the initiative of Institut Farman. The first edition of NCMIP also took place in Cachan, France, within the scope of the ValueTools Conference, in May 2011 (http://www.ncmip.org/2011/). The 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 finance. 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, convex optimization, free discontinuity problems, metamodels, proper orthogonal decomposition

FOREWORD: 3rd International Workshop on New Computational Methods for Inverse Problems (NCMIP 2013)

Science.gov (United States)

Blanc-Féraud, Laure; Joubert, Pierre-Yves

2013-10-01

aspects of inversion, Bayesian estimation, kernel methods, learning methods, convex optimization, free discontinuity problems, metamodels, proper orthogonal decomposition, reduced models for the inversion, non-linear inverse scattering, image reconstruction and restoration, and applications (bio-medical imaging, non-destructive evaluation...). NCMIP 2013 was a one-day workshop held in May 2013 which attracted around 60 attendees. Each of the submitted papers has been reviewed by three reviewers. Among the accepted papers, there are seven oral presentations, five posters and one invited poster (On a deconvolution challenge presented by C Vonesch from EPFL, Switzerland). In addition, three international speakers were invited to present a longer talk. The workshop was supported by Institut Farman (ENS Cachan, CNRS) and endorsed by the following French research networks (GDR ISIS, GDR Ondes, GDR MOA, GDR MSPC). The program committee acknowledges the following research laboratories CMLA, LMT, LSV, LURPA, SATIE. Laure Blanc-Féraud and Pierre-Yves Joubert Workshop co-chair Laure Blanc-Féraud, I3S laboratory and INRIA Nice Sophia-Antipolis, France Pierre-Yves Joubert, IEF, Paris-Sud University, CNRS, France Technical program committee Gilles Aubert, J-A Dieudonné Laboratory, CNRS and University of Nice-Sophia Antipolis, France Nabil Anwer, LURPA, ENS Cachan, France Alexandre Baussard, ENSTA Bretagne, Lab-STICC, France Marc Bonnet, ENSTA, ParisTech, France Antonin Chambolle, CMAP, Ecole Polytechnique, CNRS, France Oliver Dorn, School of Mathematics, University of Manchester, UK Cécile Durieu, SATIE, ENS Cachan, CNRS, France Gérard Favier, I3S Laboratory, University of Nice Sophia-Antipolis, France Mário Figueiredo, Instituto Superior Técnico, Lisbon, Portugal Laurent Fribourg, LSV, ENS Cachan, CNRS, France Marc Lambert, L2S Laboratory, CNRS, SupElec, Paris-Sud University, France Dominique Lesselier, L2S Laboratory, CNRS, SupElec, Paris-Sud University, France Matteo

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

International Nuclear Information System (INIS)

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

2011-01-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.

An Inverse Eigenvalue Problem for a Vibrating String with Two Dirichlet Spectra

KAUST Repository

Rundell, William; Sacks, Paul

2013-01-01

A classical inverse problem is "can you hear the density of a string clamped at both ends?" The mathematical model gives rise to an inverse Sturm-Liouville problem for the unknown density ñ, and it is well known that the answer is negative

Direct and inverse problems of infrared tomography

DEFF Research Database (Denmark)

Sizikov, Valery S.; Evseev, Vadim; Fateev, Alexander

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

Identification of the Thermophysical Properties of the Soil by Inverse Problem

OpenAIRE

Mansour , Salwa; Canot , Édouard; Muhieddine , Mohamad

2016-01-01

International audience; This paper introduces a numerical strategy to estimate the thermophysical properties of a saturated porous medium (volumetric heat capacity (ρC)s , thermal conductivity λs and porosity φ) where a phase change problem (liquid/vapor) appears due strong heating. The estimation of these properties is done by inverse problem knowing the heating curves at selected points of the medium. To solve the inverse problem, we use both the Damped Gauss Newton and the Levenberg Marqua...

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

Science.gov (United States)

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.

A hybrid algorithm for solving inverse problems in elasticity

Directory of Open Access Journals (Sweden)

Barabasz Barbara

2014-12-01

Full Text Available The paper offers a new approach to handling difficult parametric inverse problems in elasticity and thermo-elasticity, formulated as global optimization ones. The proposed strategy is composed of two phases. In the first, global phase, the stochastic hp-HGS algorithm recognizes the basins of attraction of various objective minima. In the second phase, the local objective minimizers are closer approached by steepest descent processes executed singly in each basin of attraction. The proposed complex strategy is especially dedicated to ill-posed problems with multimodal objective functionals. The strategy offers comparatively low computational and memory costs resulting from a double-adaptive technique in both forward and inverse problem domains. We provide a result on the Lipschitz continuity of the objective functional composed of the elastic energy and the boundary displacement misfits with respect to the unknown constitutive parameters. It allows common scaling of the accuracy of solving forward and inverse problems, which is the core of the introduced double-adaptive technique. The capability of the proposed method of finding multiple solutions is illustrated by a computational example which consists in restoring all feasible Young modulus distributions minimizing an objective functional in a 3D domain of a photo polymer template obtained during step and flash imprint lithography.

MAP estimators and their consistency in Bayesian nonparametric inverse problems

KAUST Repository

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

2013-01-01

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. © 2013 IOP Publishing Ltd.

Migration of scattered teleseismic body waves

Science.gov (United States)

Bostock, M. G.; Rondenay, S.

1999-06-01

The retrieval of near-receiver mantle structure from scattered waves associated with teleseismic P and S and recorded on three-component, linear seismic arrays is considered in the context of inverse scattering theory. A Ray + Born formulation is proposed which admits linearization of the forward problem and economy in the computation of the elastic wave Green's function. The high-frequency approximation further simplifies the problem by enabling (1) the use of an earth-flattened, 1-D reference model, (2) a reduction in computations to 2-D through the assumption of 2.5-D experimental geometry, and (3) band-diagonalization of the Hessian matrix in the inverse formulation. The final expressions are in a form reminiscent of the classical diffraction stack of seismic migration. Implementation of this procedure demands an accurate estimate of the scattered wave contribution to the impulse response, and thus requires the removal of both the reference wavefield and the source time signature from the raw record sections. An approximate separation of direct and scattered waves is achieved through application of the inverse free-surface transfer operator to individual station records and a Karhunen-Loeve transform to the resulting record sections. This procedure takes the full displacement field to a wave vector space wherein the first principal component of the incident wave-type section is identified with the direct wave and is used as an estimate of the source time function. The scattered displacement field is reconstituted from the remaining principal components using the forward free-surface transfer operator, and may be reduced to a scattering impulse response upon deconvolution of the source estimate. An example employing pseudo-spectral synthetic seismograms demonstrates an application of the methodology.

Computational electromagnetics and model-based inversion a modern paradigm for eddy-current nondestructive evaluation

CERN Document Server

Sabbagh, Harold A; Sabbagh, Elias H; Aldrin, John C; Knopp, Jeremy S

2013-01-01

Computational Electromagnetics and Model-Based Inversion: A Modern Paradigm for Eddy Current Nondestructive Evaluation describes the natural marriage of the computer to eddy-current NDE. Three distinct topics are emphasized in the book: (a) fundamental mathematical principles of volume-integral equations as a subset of computational electromagnetics, (b) mathematical algorithms applied to signal-processing and inverse scattering problems, and (c) applications of these two topics to problems in which real and model data are used. By showing how mathematics and the computer can solve problems more effectively than current analog practices, this book defines the modern technology of eddy-current NDE. This book will be useful to advanced students and practitioners in the fields of computational electromagnetics, electromagnetic inverse-scattering theory, nondestructive evaluation, materials evaluation and biomedical imaging. Users of eddy-current NDE technology in industries as varied as nuclear power, aerospace,...

Inverse problem in radionuclide transport

International Nuclear Information System (INIS)

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

An accurate solver for forward and inverse transport

International Nuclear Information System (INIS)

Monard, Francois; Bal, Guillaume

2010-01-01

This paper presents a robust and accurate way to solve steady-state linear transport (radiative transfer) equations numerically. Our main objective is to address the inverse transport problem, in which the optical parameters of a domain of interest are reconstructed from measurements performed at the domain's boundary. This inverse problem has important applications in medical and geophysical imaging, and more generally in any field involving high frequency waves or particles propagating in scattering environments. Stable solutions of the inverse transport problem require that the singularities of the measurement operator, which maps the optical parameters to the available measurements, be captured with sufficient accuracy. This in turn requires that the free propagation of particles be calculated with care, which is a difficult problem on a Cartesian grid. A standard discrete ordinates method is used for the direction of propagation of the particles. Our methodology to address spatial discretization is based on rotating the computational domain so that each direction of propagation is always aligned with one of the grid axes. Rotations are performed in the Fourier domain to achieve spectral accuracy. The numerical dispersion of the propagating particles is therefore minimal. As a result, the ballistic and single scattering components of the transport solution are calculated robustly and accurately. Physical blurring effects, such as small angular diffusion, are also incorporated into the numerical tool. Forward and inverse calculations performed in a two-dimensional setting exemplify the capabilities of the method. Although the methodology might not be the fastest way to solve transport equations, its physical accuracy provides us with a numerical tool to assess what can and cannot be reconstructed in inverse transport theory.

Scattering theory some old and new problems

CERN Document Server

Yafaev, Dmitri R

2000-01-01

Scattering theory is, roughly speaking, perturbation theory of self-adjoint operators on the (absolutely) continuous spectrum. It has its origin in mathematical problems of quantum mechanics and is intimately related to the theory of partial differential equations. Some recently solved problems, such as asymptotic completeness for the Schrödinger operator with long-range and multiparticle potentials, as well as open problems, are discussed. Potentials for which asymptotic completeness is violated are also constructed. This corresponds to a new class of asymptotic solutions of the time-dependent Schrödinger equation. Special attention is paid to the properties of the scattering matrix, which is the main observable of the theory. The book is addressed to readers interested in a deeper study of the subject.

Parameterization of light scattering for solving the inverse problem of determining the concentrations of the principal light scattering and absorbing admixtures in shelf waters

Directory of Open Access Journals (Sweden)

Vadim N. Pelevin

2001-12-01

Full Text Available A method for estimating the water backscattering coefficient was put forward on the basis of experimental data of diffuse attenuation coefficient for downwelling irradiance and irradiance reflectance. Calculations were carried out for open sea waters of different types and the spectral dependencies were found ("anomalous" spectra and explained. On this basis, a new model of light backscattering on particles in the sea is proposed. This model may be useful for modelling remote sensing reflectance spectra in order to solve the inverse problems of estimating the concentration of natural admixtures in shelf waters.

Efficient Monte Carlo sampling of inverse problems using a neural network-based forward—applied to GPR crosshole traveltime inversion

Science.gov (United States)

Hansen, T. M.; Cordua, K. S.

2017-12-01

Probabilistically formulated inverse problems can be solved using Monte Carlo-based sampling methods. In principle, both advanced prior information, based on for example, complex geostatistical models and non-linear forward models can be considered using such methods. However, Monte Carlo methods may be associated with huge computational costs that, in practice, limit their application. This is not least due to the computational requirements related to solving the forward problem, where the physical forward response of some earth model has to be evaluated. Here, it is suggested to replace a numerical complex evaluation of the forward problem, with a trained neural network that can be evaluated very fast. This will introduce a modeling error that is quantified probabilistically such that it can be accounted for during inversion. This allows a very fast and efficient Monte Carlo sampling of the solution to an inverse problem. We demonstrate the methodology for first arrival traveltime inversion of crosshole ground penetrating radar data. An accurate forward model, based on 2-D full-waveform modeling followed by automatic traveltime picking, is replaced by a fast neural network. This provides a sampling algorithm three orders of magnitude faster than using the accurate and computationally expensive forward model, and also considerably faster and more accurate (i.e. with better resolution), than commonly used approximate forward models. The methodology has the potential to dramatically change the complexity of non-linear and non-Gaussian inverse problems that have to be solved using Monte Carlo sampling techniques.

Bayesian inverse problems for functions and applications to fluid mechanics

International Nuclear Information System (INIS)

Cotter, S L; Dashti, M; Robinson, J C; Stuart, A M

2009-01-01

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

Moving Least Squares Method for a One-Dimensional Parabolic Inverse Problem

Directory of Open Access Journals (Sweden)

Baiyu Wang

2014-01-01

Full Text Available This paper investigates the numerical solution of a class of one-dimensional inverse parabolic problems using the moving least squares approximation; the inverse problem is the determination of an unknown source term depending on time. The collocation method is used for solving the equation; some numerical experiments are presented and discussed to illustrate the stability and high efficiency of the method.

Inverse problems using ANN in long range atmospheric dispersion with signature analysis picked scattered numerical sensors from CFD

International Nuclear Information System (INIS)

Sharma, Pavan K.; Gera, B.; Ghosh, A.K.; Kushwaha, H.S.

2010-01-01

Scalar dispersion in the atmosphere is an important area wherein different approaches are followed in development of good analytical model. The analyses based on Computational Fluid Dynamics (CFD) codes offer an opportunity of model development based on first principles of physics and hence such models have an edge over the existing models. Both forward and backward calculation methods are being developed for atmospheric dispersion around NPPs at BARC Forward modeling methods, which describe the atmospheric transport from sources to receptors, use forward-running transport and dispersion models or computational fluid dynamics models which are run many times, and the resulting dispersion field is compared to observations from multiple sensors. Backward or inverse modeling methods use only one model run in the reverse direction from the receptors to estimate the upwind sources. Inverse modeling methods include adjoint and tangent linear models, Kalman filters, and variational data assimilation, and neural network. The present paper is aimed at developing a new approach where the identified specific signatures at receptor points form the basis for source estimation or inversions. This approach is expected to reduce the large transient data sets to reduced and meaningful data sets. In fact this reduces the inherently transient data set into a time independent mean data set. Forward computation were carried out with CFD code for various case to generate a large set of data to train the ANN. Specific signature analysis was carried out to find the parameters of interest for ANN training like peak concentration, time to reach peak concentration and time to fall, the ANN was trained with data and source strength and location were predicted from ANN. Inverse problem was performed using ANN approach in long range atmospheric dispersion. An illustration of application of CFD code for atmospheric dispersion studies for a hypothetical case is also included in the paper. (author)

Coefficient Inverse Problem for Poisson's Equation in a Cylinder

NARCIS (Netherlands)

Solov'ev, V. V.

2011-01-01

The inverse problem of determining the coefficient on the right-hand side of Poisson's equation in a cylindrical domain is considered. The Dirichlet boundary value problem is studied. Two types of additional information (overdetermination) can be specified: (i) the trace of the solution to the

NUMERICAL SOLUTION OF SINGULAR INVERSE NODAL PROBLEM BY USING CHEBYSHEV POLYNOMIALS

OpenAIRE

NEAMATY, ABDOLALI; YILMAZ, EMRAH; AKBARPOOR, SHAHRBANOO; DABBAGHIAN, ABDOLHADI

2017-01-01

In this study, we consider Sturm-Liouville problem in two cases: the first case having no singularity and the second case having a singularity at zero. Then, we calculate the eigenvalues and the nodal points and present the uniqueness theorem for the solution of the inverse problem by using a dense subset of the nodal points in two given cases. Also, we use Chebyshev polynomials of the first kind for calculating the approximate solution of the inverse nodal problem in these cases. Finally, we...

An inverse heat transfer problem for optimization of the thermal ...

Indian Academy of Sciences (India)

This paper takes a different approach towards identification of the thermal process in machining, using inverse heat transfer problem. Inverse heat transfer method allows the closest possible experimental and analytical approximation of thermal state for a machining process. Based on a temperature measured at any point ...

Invariant imbedding equations for linear scattering problems

International Nuclear Information System (INIS)

Apresyan, L.

1988-01-01

A general form of the invariant imbedding equations is investigated for the linear problem of scattering by a bounded scattering volume. The conditions for the derivability of such equations are described. It is noted that the possibility of the explicit representation of these equations for a sphere and for a layer involves the separation of variables in the unperturbed wave equation

Connection of Scattering Principles: A Visual and Mathematical Tour

Science.gov (United States)

Broggini, Filippo; Snieder, Roel

2012-01-01

Inverse scattering, Green's function reconstruction, focusing, imaging and the optical theorem are subjects usually studied as separate problems in different research areas. We show a physical connection between the principles because the equations that rule these "scattering principles" have a similar functional form. We first lead the reader…

Differential equations inverse and direct problems

CERN Document Server

Favini, Angelo

2006-01-01

DEGENERATE FIRST ORDER IDENTIFICATION PROBLEMS IN BANACH SPACES A NONISOTHERMAL DYNAMICAL GINZBURG-LANDAU MODEL OF SUPERCONDUCTIVITY. EXISTENCE AND UNIQUENESS THEOREMSSOME GLOBAL IN TIME RESULTS FOR INTEGRODIFFERENTIAL PARABOLIC INVERSE PROBLEMSFOURTH ORDER ORDINARY DIFFERENTIAL OPERATORS WITH GENERAL WENTZELL BOUNDARY CONDITIONSTUDY OF ELLIPTIC DIFFERENTIAL EQUATIONS IN UMD SPACESDEGENERATE INTEGRODIFFERENTIAL EQUATIONS OF PARABOLIC TYPE EXPONENTIAL ATTRACTORS FOR SEMICONDUCTOR EQUATIONSCONVERGENCE TO STATIONARY STATES OF SOLUTIONS TO THE SEMILINEAR EQUATION OF VISCOELASTICITY ASYMPTOTIC BEHA

Numerical solution of the multichannel scattering problem

International Nuclear Information System (INIS)

Korobov, V.I.

1992-01-01

A numerical algorithm for solving the multichannel elastic and inelastic scattering problem is proposed. The starting point is the system of radial Schroedinger equations with linear boundary conditions imposed at some point R=R m placed somewhere in asymptotic region. It is discussed how the obtained linear equation can be splitted into a zero-order operator and its pertturbative part. It is shown that Lentini - Pereyra variable order finite-difference method appears to be very suitable for solving that kind of problems. The derived procedure is applied to dμ+t→tμ+d inelastic scattering in the framework of the adiabatic multichannel approach. 19 refs.; 1 fig.; 1 tab

Solving of L0 norm constrained EEG inverse problem.

Science.gov (United States)

Xu, Peng; Lei, Xu; Hu, Xiao; Yao, Dezhong

2009-01-01

l(0) norm is an effective constraint used to solve EEG inverse problem for a sparse solution. However, due to the discontinuous and un-differentiable properties, it is an open issue to solve the l(0) norm constrained problem, which is usually instead solved by using some alternative functions like l(1) norm to approximate l(0) norm. In this paper, a continuous and differentiable function having the same form as the transfer function of Butterworth low-pass filter is introduced to approximate l(0) norm constraint involved in EEG inverse problem. The new approximation based approach was compared with l(1) norm and LORETA solutions on a realistic head model using simulated sources. The preliminary results show that this alternative approximation to l(0) norm is promising for the estimation of EEG sources with sparse distribution.

On multiple level-set regularization methods for inverse problems

International Nuclear Information System (INIS)

DeCezaro, A; Leitão, A; Tai, X-C

2009-01-01

We analyze a multiple level-set method for solving inverse problems with piecewise constant solutions. This method corresponds to an iterated Tikhonov method for a particular Tikhonov functional G α based on TV–H 1 penalization. We define generalized minimizers for our Tikhonov functional and establish an existence result. Moreover, we prove convergence and stability results of the proposed Tikhonov method. A multiple level-set algorithm is derived from the first-order optimality conditions for the Tikhonov functional G α , similarly as the iterated Tikhonov method. The proposed multiple level-set method is tested on an inverse potential problem. Numerical experiments show that the method is able to recover multiple objects as well as multiple contrast levels

An inverse analysis of a transient 2-D conduction-radiation problem using the lattice Boltzmann method and the finite volume method coupled with the genetic algorithm

International Nuclear Information System (INIS)

Das, Ranjan; Mishra, Subhash C.; Ajith, M.; Uppaluri, R.

2008-01-01

This article deals with the simultaneous estimation of parameters in a 2-D transient conduction-radiation heat transfer problem. The homogeneous medium is assumed to be absorbing, emitting and scattering. The boundaries of the enclosure are diffuse gray. Three parameters, viz. the scattering albedo, the conduction-radiation parameter and the boundary emissivity, are simultaneously estimated by the inverse method involving the lattice Boltzmann method (LBM) and the finite volume method (FVM) in conjunction with the genetic algorithm (GA). In the direct method, the FVM is used for computing the radiative information while the LBM is used to solve the energy equation. The temperature field obtained in the direct method is used in the inverse method for simultaneous estimation of unknown parameters using the LBM-FVM and the GA. The LBM-FVM-GA combination has been found to accurately predict the unknown parameters

Inverse problem for in vivo NMR spatial localization

International Nuclear Information System (INIS)

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

Inverse problem for in vivo NMR spatial localization

Energy Technology Data Exchange (ETDEWEB)

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.

Inverse Problem in Self-assembly

Science.gov (United States)

Tkachenko, Alexei

2012-02-01

By decorating colloids and nanoparticles with DNA, one can introduce highly selective key-lock interactions between them. This leads to a new class of systems and problems in soft condensed matter physics. In particular, this opens a possibility to solve inverse problem in self-assembly: how to build an arbitrary desired structure with the bottom-up approach? I will present a theoretical and computational analysis of the hierarchical strategy in attacking this problem. It involves self-assembly of particular building blocks (``octopus particles''), that in turn would assemble into the target structure. On a conceptual level, our approach combines elements of three different brands of programmable self assembly: DNA nanotechnology, nanoparticle-DNA assemblies and patchy colloids. I will discuss the general design principles, theoretical and practical limitations of this approach, and illustrate them with our simulation results. Our crucial result is that not only it is possible to design a system that has a given nanostructure as a ground state, but one can also program and optimize the kinetic pathway for its self-assembly.

Inverse problem theory methods for data fitting and model parameter estimation

CERN Document Server

Tarantola, A

2002-01-01

Inverse Problem Theory is written for physicists, geophysicists and all scientists facing the problem of quantitative interpretation of experimental data. Although it contains a lot of mathematics, it is not intended as a mathematical book, but rather tries to explain how a method of acquisition of information can be applied to the actual world.The book provides a comprehensive, up-to-date description of the methods to be used for fitting experimental data, or to estimate model parameters, and to unify these methods into the Inverse Problem Theory. The first part of the book deals wi

Numerical modelling of multiple scattering between two elastical particles

DEFF Research Database (Denmark)

Bjørnø, Irina; Jensen, Leif Bjørnø

1998-01-01

in suspension have been studied extensively since Foldy's formulation of his theory for isotropic scattering by randomly distributed scatterers. However, a number of important problems related to multiple scattering are still far from finding their solutions. A particular, but still unsolved, problem......Multiple acoustical signal interactions with sediment particles in the vicinity of the seabed may significantly change the course of sediment concentration profiles determined by inversion from acoustical backscattering measurements. The scattering properties of high concentrations of sediments...... is the question of proximity thresholds for influence of multiple scattering in terms of particle properties like volume fraction, average distance between particles or other related parameters. A few available experimental data indicate a significance of multiple scattering in suspensions where the concentration...

THE DIDACTIC ANALYSIS OF STUDIES ON THE INVERSE PROBLEMS FOR THE DIFFERENTIAL EQUATIONS

Directory of Open Access Journals (Sweden)

В С Корнилов

2017-12-01

Full Text Available In article results of the didactic analysis of the organization and carrying out seminar classes in the inverse problems for the differential equations for students of higher educational institutions of the physical and mathematical directions of preparation are discussed. Such analysis includes a general characteristic of mathematical content of seminar occupations, the analysis of structure of seminar occupation, the analysis of realization of the developing and educational purposes, allocation of didactic units and informative means which have to be acquired by students when training each section of content of training in the inverse problems and other important psychology and pedagogical aspects. The attention to establishment of compliance to those of seminar occupations to lecture material and identification of functions in teaching and educational process which are carried out at the solution of the inverse problems, and also is paid to need to show various mathematical receptions and methods of their decision. Such didactic analysis helps not only to reveal such inverse problems at which solution students can collectively join in creative process of search of their decision, but also effectively organize control of assimilation of knowledge and abilities of students on the inverse problems for the differential equations.

Solving probabilistic inverse problems rapidly with prior samples

NARCIS (Netherlands)

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

2016-01-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

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

Energy Technology Data Exchange (ETDEWEB)

NONE

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.

Reconstruction of absorption and scattering coefficients in two dimensional heterogeneous participating media

International Nuclear Information System (INIS)

Montero, Raul F. Carita; Roberty, Nilson C.; Silva Neto, Antonio J.; Universidade Federal, Rio de Janeiro, RJ

2002-01-01

In the present work it is presented the solution of the two dimensional inverse radiative transfer problem of scattering and absorption coefficients estimation, in heterogeneous media, using the source-detector methodology and a discrete ordinates method consistent with the source-detector system. The mathematical formulation of the direct and inverse problems is presented as well as test case results. (author)

An optical potential from inversion of the 350 MeV 16 O - 16 O scattering data

International Nuclear Information System (INIS)

Allen, L.J.; Berge, L.; Steward, C.; Lipperheide, R.; Froebrich, P.

1992-01-01

A quantal inversion of the 16 O- 16 O scattering data at 350 MeV yields an optical potential which gives an excellent fit (χ 2 /F = 1.65) to the measured cross-section. The real part of this potential is shallower than any potential used by others for distances between 2 and 6 fm. The imaginary potential is also relatively weak. This potential does not favour a rainbow interpretation of the structure in data observed at large scattering angles. 12 refs., 1 tab., 4 figs

A Frequency Matching Method: Solving Inverse Problems by Use of Geologically Realistic Prior Information

DEFF Research Database (Denmark)

Lange, Katrine; Frydendall, Jan; Cordua, Knud Skou

2012-01-01

The frequency matching method defines a closed form expression for a complex prior that quantifies the higher order statistics of a proposed solution model to an inverse problem. While existing solution methods to inverse problems are capable of sampling the solution space while taking into account...... arbitrarily complex a priori information defined by sample algorithms, it is not possible to directly compute the maximum a posteriori model, as the prior probability of a solution model cannot be expressed. We demonstrate how the frequency matching method enables us to compute the maximum a posteriori...... solution model to an inverse problem by using a priori information based on multiple point statistics learned from training images. We demonstrate the applicability of the suggested method on a synthetic tomographic crosshole inverse problem....

Synthetic acceleration methods for linear transport problems with highly anisotropic scattering

International Nuclear Information System (INIS)

Khattab, K.M.; Larsen, E.W.

1992-01-01

The diffusion synthetic acceleration (DSA) algorithm effectively accelerates the iterative solution of transport problems with isotropic or mildly anisotropic scattering. However, DSA loses its effectiveness for transport problems that have strongly anisotropic scattering. Two generalizations of DSA are proposed, which, for highly anisotropic scattering problems, converge at least an order of magnitude (clock time) faster than the DSA method. These two methods are developed, the results of Fourier analysis that theoretically predict their efficiency are described, and numerical results that verify the theoretical predictions are presented. (author). 10 refs., 7 figs., 5 tabs

Synthetic acceleration methods for linear transport problems with highly anisotropic scattering

International Nuclear Information System (INIS)

Khattab, K.M.; Larsen, E.W.

1991-01-01

This paper reports on the diffusion synthetic acceleration (DSA) algorithm that effectively accelerates the iterative solution of transport problems with isotropic or mildly anisotropic scattering. However, DSA loses its effectiveness for transport problems that have strongly anisotropic scattering. Two generalizations of DSA are proposed, which, for highly anisotropic scattering problems, converge at least an order of magnitude (clock time) faster than the DSA method. These two methods are developed, the results of Fourier analyses that theoretically predict their efficiency are described, and numerical results that verify the theoretical predictions are presented

An inverse optimal control problem in the electrical discharge ...

Indian Academy of Sciences (India)

Marin Gostimirovic

2018-05-10

May 10, 2018 ... Keywords. EDM process; discharge energy; heat source parameters; inverse problem; optimization. 1. Introduction .... ation, thermal modeling of the EDM process would become ..... simulation of die-sinking EDM. CIRP Ann.

Hybrid inverse problems for a system of Maxwell’s equations

International Nuclear Information System (INIS)

Bal, Guillaume; Zhou, Ting

2014-01-01

This paper concerns the quantitative step of the medical imaging modality thermo-acoustic tomography (TAT). We model the radiation propagation by a system of Maxwell’s equations. We show that the index of refraction of light and the absorption coefficient (conductivity) can be uniquely and stably reconstructed from a sufficiently large number of TAT measurements. Our method is based on verifying that the linearization of the inverse problem forms a redundant elliptic system of equations. We also observe that the reconstructions are qualitatively quite different from the setting where radiation is modeled by a scalar Helmholtz equation as in Bal G et al (2011 Inverse Problems 27 055007). (paper)

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

International Nuclear Information System (INIS)

Hohage, Thorsten; Werner, Frank

2016-01-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. (topical review)

A study of light scattering of mononuclear blood cells with scanning flow cytometry

International Nuclear Information System (INIS)

Zharinov, Alexey; Tarasov, Peter; Shvalov, Alexander; Semyanov, Konstantin; Bockstaele, Dirk R. van; Maltsev, Valeri

2006-01-01

This study describes the measurement of light scattering of human mononuclear blood cells, the development of an appropriate optical model for those cells, and solution of the inverse light-scattering problem. The angular dependency of light-scattering intensity of mononuclear blood cells was experimentally measured by means of scanning flow cytometry. A sphere consisting of several concentric homogeneous layers with different refractive indices was tested as an optical model for mononuclear blood cells. A five-layer model has given the best agreement between experimental and theoretical light-scattering profiles. The inverse light-scattering problem was solved for a five-layer model with an optimization procedure that allows one to retrieve cell parameters: cell size relates to the outer diameter of the fifth layer; size of the nucleus relates to the outer diameter of the third layer. Mean values of cell size, nuclear size, refractive indices of nucleus and cellular cytoplasm were determined for blood monocytes and lymphocytes

An inverse problem strategy based on forward model evaluations: Gradient-based optimization without adjoint solves

Energy Technology Data Exchange (ETDEWEB)

Aguilo Valentin, Miguel Alejandro [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2016-07-01

This study presents a new nonlinear programming formulation for the solution of inverse problems. First, a general inverse problem formulation based on the compliance error functional is presented. The proposed error functional enables the computation of the Lagrange multipliers, and thus the first order derivative information, at the expense of just one model evaluation. Therefore, the calculation of the Lagrange multipliers does not require the solution of the computationally intensive adjoint problem. This leads to significant speedups for large-scale, gradient-based inverse problems.

Collage-type approach to inverse problems for elliptic PDEs on perforated domains

Directory of Open Access Journals (Sweden)

Herb E. Kunze

2015-02-01

Full Text Available We present a collage-based method for solving inverse problems for elliptic partial differential equations on a perforated domain. The main results of this paper establish a link between the solution of an inverse problem on a perforated domain and the solution of the same model on a domain with no holes. The numerical examples at the end of the paper show the goodness of this approach.

A gradient based algorithm to solve inverse plane bimodular problems of identification

Science.gov (United States)

Ran, Chunjiang; Yang, Haitian; Zhang, Guoqing

2018-02-01

This paper presents a gradient based algorithm to solve inverse plane bimodular problems of identifying constitutive parameters, including tensile/compressive moduli and tensile/compressive Poisson's ratios. For the forward bimodular problem, a FE tangent stiffness matrix is derived facilitating the implementation of gradient based algorithms, for the inverse bimodular problem of identification, a two-level sensitivity analysis based strategy is proposed. Numerical verification in term of accuracy and efficiency is provided, and the impacts of initial guess, number of measurement points, regional inhomogeneity, and noisy data on the identification are taken into accounts.

Physics-based models for measurement correlations: application to an inverse Sturm–Liouville problem

International Nuclear Information System (INIS)

Bal, Guillaume; Ren Kui

2009-01-01

In many inverse problems, the measurement operator, which maps objects of interest to available measurements, is a smoothing (regularizing) operator. Its inverse is therefore unbounded and as a consequence, only the low-frequency component of the object of interest is accessible from inevitably noisy measurements. In many inverse problems however, the neglected high-frequency component may significantly affect the measured data. Using simple scaling arguments, we characterize the influence of the high-frequency component. We then consider situations where the correlation function of such an influence may be estimated by asymptotic expansions, for instance as a random corrector in homogenization theory. This allows us to consistently eliminate the high-frequency component and derive a closed form, more accurate, inverse problem for the low-frequency component of the object of interest. We present the asymptotic expression of the correlation matrix of the eigenvalues in a Sturm–Liouville problem with unknown potential. We propose an iterative algorithm for the reconstruction of the potential from knowledge of the eigenvalues and show that using the approximate correlation matrix significantly improves the reconstructions

Inelastic Proton Scattering on 21Na in Inverse Kinematics

Science.gov (United States)

Austin, Roby

2009-10-01

R.A.E. Austin, R. Kanungo, S. Reeve, Saint Mary's University; D.G. Jenkins, C.Aa.Diget, A. Robinson, A.G. Tuff, O. Roberts, University of York, UK; P.J. Woods, T. Davinson, G. J. Lotay, University of Edinburgh; C.-Y. Wu, Lawrence Livermore National Laboratory; H. Al Falou, G.C. Ball, M. Djongolov, A. Garnsworthy, G. Hackman, J.N. Orce, C.J. Pearson, S. Triambak, S.J. Williams, TRIUMF; C. Andreiou, D.S. Cross, N. Galinski, R. Kshetri, Simon Fraser University; C. Sumithrarachchi, M.A. Schumaker, University of Guelph; M.P. Jones, S.V. Rigby, University of Liverpool; D. Cline, A. Hayes, University of Rochester; T.E. Drake, University of Toronto; We describe an experiment and associated technique [1] to measure resonances of interest in astrophysical reactions. At the TRIUMF ISAC-II radioactive beam accelerator facility in Canada, particles inelastically scattered in inverse kinematics are detected with Bambino, a δE-E silicon telescope spanning 15-40 degrees in the lab. We use the TIGRESS to detect gamma rays in coincidence with the charged particles to cleanly select inelastic scattering events. We measured resonances above the alpha threshold in ^22Mg of relevance to the rate of break-out from the hot-CNO cycle via the reaction ^ 18Ne(α,p)^21Na. [1] PJ Woods et al. Rex-ISOLDE proposal 424 Cern (2003).

Detection of inverse Compton scattering in plasma wakefield experiments

Energy Technology Data Exchange (ETDEWEB)

Bohlen, Simon

2016-12-15

Inverse Compton scattering (ICS) is the process of scattering of photons and electrons, where the photons gain a part of the electrons energy. In combination with plasma wakefield acceleration (PWA), ICS offers a compact MeV γ-ray source. A numerical study of ICS radiation produced in PWA experiments at FLASHForward was performed, using an ICS simulation code and the results from particle-in-cell modelling. The possibility of determining electron beam properties from measurements of the γ-ray source was explored for a wide range of experimental conditions. It was found that information about the electron divergence, the electron spectrum and longitudinal information can be obtained from measurements of the ICS beams for some cases. For the measurement of the ICS profile at FLASHForward, a CsI(Tl) scintillator array was chosen, similar to scintillators used in other ICS experiments. To find a suitable detector for spectrum measurements, an experimental test of a Compton spectrometer at the RAL was conducted. This test showed that a similar spectrometer could also be used at FLASHForward. However, changes to the spectrometer could be needed in order to use the pair production effect. In addition, further studies using Geant4 could lead to a better reconstruction of the obtained data. The studies presented here show that ICS is a promising method to analyse electron parameters from PWA experiments in further detail.

Spectral solution of the inverse Mie problem

Science.gov (United States)

Romanov, Andrey V.; Konokhova, Anastasiya I.; Yastrebova, Ekaterina S.; Gilev, Konstantin V.; Strokotov, Dmitry I.; Chernyshev, Andrei V.; Maltsev, Valeri P.; Yurkin, Maxim A.

2017-10-01

We developed a fast method to determine size and refractive index of homogeneous spheres from the power Fourier spectrum of their light-scattering patterns (LSPs), measured with the scanning flow cytometer. Specifically, we used two spectral parameters: the location of the non-zero peak and zero-frequency amplitude, and numerically inverted the map from the space of particle characteristics (size and refractive index) to the space of spectral parameters. The latter parameters can be reliably resolved only for particle size parameter greater than 11, and the inversion is unique only in the limited range of refractive index with upper limit between 1.1 and 1.25 (relative to the medium) depending on the size parameter and particular definition of uniqueness. The developed method was tested on two experimental samples, milk fat globules and spherized red blood cells, and resulted in accuracy not worse than the reference method based on the least-square fit of the LSP with the Mie theory. Moreover, for particles with significant deviation from the spherical shape the spectral method was much closer to the Mie-fit result than the estimated uncertainty of the latter. The spectral method also showed adequate results for synthetic LSPs of spheroids with aspect ratios up to 1.4. Overall, we present a general framework, which can be used to construct an inverse algorithm for any other experimental signals.

Data-Driven Model Order Reduction for Bayesian Inverse Problems

KAUST Repository

Cui, Tiangang; Youssef, Marzouk; Willcox, Karen

2014-01-01

One of the major challenges in using MCMC for the solution of inverse problems is the repeated evaluation of computationally expensive numerical models. We develop a data-driven projection- based model order reduction technique to reduce

Electromagnetic imaging of multiple-scattering small objects: non-iterative analytical approach

International Nuclear Information System (INIS)

Chen, X; Zhong, Y

2008-01-01

Multiple signal classification (MUSIC) imaging method and the least squares method are applied to solve the electromagnetic inverse scattering problem of determining the locations and polarization tensors of a collection of small objects embedded in a known background medium. Based on the analysis of induced electric and magnetic dipoles, the proposed MUSIC method is able to deal with some special scenarios, due to the shapes and materials of objects, to which the standard MUSIC doesn't apply. After the locations of objects are obtained, the nonlinear inverse problem of determining the polarization tensors of objects accounting for multiple scattering between objects is solved by a non-iterative analytical approach based on the least squares method

Solving inverse problems through a smooth formulation of multiple-point geostatistics

DEFF Research Database (Denmark)

Melnikova, Yulia

be inferred, for instance, from a conceptual geological model termed a training image.The main motivation for this study was the challenge posed by history matching, an inverse problem aimed at estimating rock properties from production data. We addressed two main difficulties of the history matching problem...... corresponding inverse problems. However, noise in data, non-linear relationships and sparse observations impede creation of realistic reservoir models. Including complex a priori information on reservoir parameters facilitates the process of obtaining acceptable solutions. Such a priori knowledge may...... strategies including both theoretical motivation and practical aspects of implementation. Finally, it is complemented by six research papers submitted, reviewed and/or published in the period 2010 - 2013....

Well-posedness of inverse problems for systems with time dependent parameters

DEFF Research Database (Denmark)

Banks, H. T.; Pedersen, Michael

2009-01-01

on the data of the problem. We also consider well-posedness as well as finite element type approximations in associated inverse problems. The problem above is a weak formulation that includes models in abstract differential operator form that include plate, beam and shell equations with several important...

A penalty method for PDE-constrained optimization in inverse problems

International Nuclear Information System (INIS)

Leeuwen, T van; Herrmann, F J

2016-01-01

Many inverse and parameter estimation problems can be written as PDE-constrained optimization problems. The goal is to infer the parameters, typically coefficients of the PDE, from partial measurements of the solutions of the PDE for several right-hand sides. Such PDE-constrained problems can be solved by finding a stationary point of the Lagrangian, which entails simultaneously updating the parameters and the (adjoint) state variables. For large-scale problems, such an all-at-once approach is not feasible as it requires storing all the state variables. In this case one usually resorts to a reduced approach where the constraints are explicitly eliminated (at each iteration) by solving the PDEs. These two approaches, and variations thereof, are the main workhorses for solving PDE-constrained optimization problems arising from inverse problems. In this paper, we present an alternative method that aims to combine the advantages of both approaches. Our method is based on a quadratic penalty formulation of the constrained optimization problem. By eliminating the state variable, we develop an efficient algorithm that has roughly the same computational complexity as the conventional reduced approach while exploiting a larger search space. Numerical results show that this method indeed reduces some of the nonlinearity of the problem and is less sensitive to the initial iterate. (paper)

Comparison of matrix methods for elastic wave scattering problems

International Nuclear Information System (INIS)

Tsao, S.J.; Varadan, V.K.; Varadan, V.V.

1983-01-01

This article briefly describes the T-matrix method and the MOOT (method of optimal truncation) of elastic wave scattering as they apply to A-D, SH- wave problems as well as 3-D elastic wave problems. Two methods are compared for scattering by elliptical cylinders as well as oblate spheroids of various eccentricity as a function of frequency. Convergence, and symmetry of the scattering cross section are also compared for ellipses and spheroidal cavities of different aspect ratios. Both the T-matrix approach and the MOOT were programmed on an AMDHL 470 computer using double precision arithmetic. Although the T-matrix method and MOOT are not always in agreement, it is in no way implied that any of the published results using MOOT are in error

Automatic differentiation in geophysical inverse problems

Science.gov (United States)

Sambridge, M.; Rickwood, P.; Rawlinson, N.; Sommacal, S.

2007-07-01

Automatic differentiation (AD) is the technique whereby output variables of a computer code evaluating any complicated function (e.g. the solution to a differential equation) can be differentiated with respect to the input variables. Often AD tools take the form of source to source translators and produce computer code without the need for deriving and hand coding of explicit mathematical formulae by the user. The power of AD lies in the fact that it combines the generality of finite difference techniques and the accuracy and efficiency of analytical derivatives, while at the same time eliminating `human' coding errors. It also provides the possibility of accurate, efficient derivative calculation from complex `forward' codes where no analytical derivatives are possible and finite difference techniques are too cumbersome. AD is already having a major impact in areas such as optimization, meteorology and oceanography. Similarly it has considerable potential for use in non-linear inverse problems in geophysics where linearization is desirable, or for sensitivity analysis of large numerical simulation codes, for example, wave propagation and geodynamic modelling. At present, however, AD tools appear to be little used in the geosciences. Here we report on experiments using a state of the art AD tool to perform source to source code translation in a range of geoscience problems. These include calculating derivatives for Gibbs free energy minimization, seismic receiver function inversion, and seismic ray tracing. Issues of accuracy and efficiency are discussed.

Resonant Inverse Compton Scattering Spectra from Highly Magnetized Neutron Stars

Science.gov (United States)

Wadiasingh, Zorawar; Baring, Matthew G.; Gonthier, Peter L.; Harding, Alice K.

2018-02-01

Hard, nonthermal, persistent pulsed X-ray emission extending between 10 and ∼150 keV has been observed in nearly 10 magnetars. For inner-magnetospheric models of such emission, resonant inverse Compton scattering of soft thermal photons by ultrarelativistic charges is the most efficient production mechanism. We present angle-dependent upscattering spectra and pulsed intensity maps for uncooled, relativistic electrons injected in inner regions of magnetar magnetospheres, calculated using collisional integrals over field loops. Our computations employ a new formulation of the QED Compton scattering cross section in strong magnetic fields that is physically correct for treating important spin-dependent effects in the cyclotron resonance, thereby producing correct photon spectra. The spectral cutoff energies are sensitive to the choices of observer viewing geometry, electron Lorentz factor, and scattering kinematics. We find that electrons with energies ≲15 MeV will emit most of their radiation below 250 keV, consistent with inferred turnovers for magnetar hard X-ray tails. More energetic electrons still emit mostly below 1 MeV, except for viewing perspectives sampling field-line tangents. Pulse profiles may be singly or doubly peaked dependent on viewing geometry, emission locale, and observed energy band. Magnetic pair production and photon splitting will attenuate spectra to hard X-ray energies, suppressing signals in the Fermi-LAT band. The resonant Compton spectra are strongly polarized, suggesting that hard X-ray polarimetry instruments such as X-Calibur, or a future Compton telescope, can prove central to constraining model geometry and physics.

Quantum-mechanical scattering in one dimension

International Nuclear Information System (INIS)

Boya, Luis J.

2008-01-01

The purpose of this mainly pedagogical review is to fill a lacuna in the usual treatment of scattering in quantum mechanics, by showing the essential of it in the simplest, one-dimensional setting. We define in this situation amplitudes and scattering coefficients and deal with optical and Levinson' theorems as consequences of unitarity in coordinate or momentum space. Parity waves en lieu of partial waves, integral equations and Born series, etc., are defined naturally in this frame. Several solvable examples are shown. Two topics best studied in 1d are transparent potentials and supersymmetric quantum mechanics. Elementary analytical properties and general behaviour of amplitudes give rise to study inverse problems, that is, recovering the potential from scattering data. Isospectral deformations of the wave equation give relations with some nonlinear evolution equations (Lax), solvable by the inverse scattering method (Kruskal), and we consider the KdV equation as an example. We also refer briefly to some singular potentials, where, e.g., the essence of renormalization can be read off again in the simplest setting. The whole paper emphasizes the tutorial and introductory aspects

Analytical description of photon beam phase spaces in inverse Compton scattering sources

Directory of Open Access Journals (Sweden)

C. Curatolo

2017-08-01

Full Text Available We revisit the description of inverse Compton scattering sources and the photon beams generated therein, emphasizing the behavior of their phase space density distributions and how they depend upon those of the two colliding beams of electrons and photons. The main objective is to provide practical formulas for bandwidth, spectral density, brilliance, which are valid in general for any value of the recoil factor, i.e. both in the Thomson regime of negligible electron recoil, and in the deep Compton recoil dominated region, which is of interest for gamma-gamma colliders and Compton sources for the production of multi-GeV photon beams. We adopt a description based on the center of mass reference system of the electron-photon collision, in order to underline the role of the electron recoil and how it controls the relativistic Doppler/boost effect in various regimes. Using the center of mass reference frame greatly simplifies the treatment, allowing us to derive simple formulas expressed in terms of rms momenta of the two colliding beams (emittance, energy spread, etc. and the collimation angle in the laboratory system. Comparisons with Monte Carlo simulations of inverse Compton scattering in various scenarios are presented, showing very good agreement with the analytical formulas: in particular we find that the bandwidth dependence on the electron beam emittance, of paramount importance in Thomson regime, as it limits the amount of focusing imparted to the electron beam, becomes much less sensitive in deep Compton regime, allowing a stronger focusing of the electron beam to enhance luminosity without loss of mono-chromaticity. A similar effect occurs concerning the bandwidth dependence on the frequency spread of the incident photons: in deep recoil regime the bandwidth comes out to be much less dependent on the frequency spread. The set of formulas here derived are very helpful in designing inverse Compton sources in diverse regimes, giving a

Accounting for imperfect forward modeling in geophysical inverse problems — Exemplified for crosshole tomography

DEFF Research Database (Denmark)

Hansen, Thomas Mejer; Cordua, Knud Skou; Holm Jacobsen, Bo

2014-01-01

forward models, can be more than an order of magnitude larger than the measurement uncertainty. We also found that the modeling error is strongly linked to the spatial variability of the assumed velocity field, i.e., the a priori velocity model.We discovered some general tools by which the modeling error...... synthetic ground-penetrating radar crosshole tomographic inverse problems. Ignoring the modeling error can lead to severe artifacts, which erroneously appear to be well resolved in the solution of the inverse problem. Accounting for the modeling error leads to a solution of the inverse problem consistent...

The Cauchy problem for the Pavlov equation

International Nuclear Information System (INIS)

Grinevich, P G; Santini, P M; Wu, D

2015-01-01

Commutation of multidimensional vector fields leads to integrable nonlinear dispersionless PDEs that arise in various problems of mathematical physics and have been intensively studied in recent literature. This report aims to solve the scattering and inverse scattering problem for integrable dispersionless PDEs, recently introduced just at a formal level, concentrating on the prototypical example of the Pavlov equation, and to justify an existence theorem for global bounded solutions of the associated Cauchy problem with small data. (paper)

LinvPy : a Python package for linear inverse problems

OpenAIRE

Beaud, Guillaume François Paul

2016-01-01

The goal of this project is to make a Python package including the tau-estimator algorithm to solve linear inverse problems. The package must be distributed, well documented, easy to use and easy to extend for future developers.

Multiscattering inversion for low-model wavenumbers

KAUST Repository

Alkhalifah, Tariq Ali

2016-09-21

A successful full-waveform inversion implementation updates the low-wavenumber model components first for a proper description of the wavefield propagation and slowly adds the high wavenumber potentially scattering parts of the model. The low-wavenumber components can be extracted from the transmission parts of the recorded wavefield emanating directly from the source or the transmission parts from the single- or double-scattered wavefield computed from a predicted scatter field acting as secondary sources.We use a combined inversion of data modeled from the source and those corresponding to single and double scattering to update the velocity model and the component of the velocity (perturbation) responsible for the single and double scattering. The combined inversion helps us access most of the potential model wavenumber information that may be embedded in the data. A scattering-angle filter is used to divide the gradient of the combined inversion, so initially the high-wavenumber (low-scattering-angle) components of the gradient are directed to the perturbation model and the low-wavenumber (highscattering- angle) components are directed to the velocity model. As our background velocity matures, the scatteringangle divide is slowly lowered to allow for more of the higher wavenumbers to contribute the velocity model. Synthetic examples including the Marmousi model are used to demonstrate the additional illumination and improved velocity inversion obtained when including multiscattered energy. © 2016 Society of Exploration Geophysicists.

Geometric approach to inverse scattering for the Schroedinger equation with magnetic and electric potentials

International Nuclear Information System (INIS)

Arians, S.

1997-01-01

We consider the Hamiltonian H=(p-A(x)) 2 /(2m)+V(x) of a quantum particle in a magnetic field B=rotA and a potential V in space dimensions ν≥2. If V is of short range, then the high-velocity limit of the scattering operator uniquely determines the magnetic field B and the potential V. If, in addition, long-range potentials V l are present, some knowledge of (the far out tail of) V l is needed to define a modified Dollard wave operator and a scattering operator S D . Again its high- velocity limit uniquely determines B and V=V s +V l . Moreover, we give explicit error bounds which are inverse proportional to the velocity. copyright 1997 American Institute of Physics

Point source reconstruction principle of linear inverse problems

International Nuclear Information System (INIS)

Terazono, Yasushi; Matani, Ayumu; Fujimaki, Norio; Murata, Tsutomu

2010-01-01

Exact point source reconstruction for underdetermined linear inverse problems with a block-wise structure was studied. In a block-wise problem, elements of a source vector are partitioned into blocks. Accordingly, a leadfield matrix, which represents the forward observation process, is also partitioned into blocks. A point source is a source having only one nonzero block. An example of such a problem is current distribution estimation in electroencephalography and magnetoencephalography, where a source vector represents a vector field and a point source represents a single current dipole. In this study, the block-wise norm, a block-wise extension of the l p -norm, was defined as the family of cost functions of the inverse method. The main result is that a set of three conditions was found to be necessary and sufficient for block-wise norm minimization to ensure exact point source reconstruction for any leadfield matrix that admit such reconstruction. The block-wise norm that satisfies the conditions is the sum of the cost of all the observations of source blocks, or in other words, the block-wisely extended leadfield-weighted l 1 -norm. Additional results are that minimization of such a norm always provides block-wisely sparse solutions and that its solutions form cones in source space

Solving inverse problems of mathematical physics by means of the PHOENICS software package

Energy Technology Data Exchange (ETDEWEB)

Matsevity, Y; Lushpenko, S [Institute for Problems in Machinery, National Academy of Sciences of Ukraine Pozharskogo, Kharkov (Ukraine)

1998-12-31

Several approaches on organizing solution of inverse problems by means of PHOENICS on the basis of the technique of automated fitting are proposing. A version of a `nondestructive` method of using PHOENICS in the inverse problem solution regime and the ways of altering the program in the case of introducing optimization facilities in it are under consideration. (author) 12 refs.

Solving inverse problems of mathematical physics by means of the PHOENICS software package

Energy Technology Data Exchange (ETDEWEB)

Matsevity, Y.; Lushpenko, S. [Institute for Problems in Machinery, National Academy of Sciences of Ukraine Pozharskogo, Kharkov (Ukraine)

1997-12-31

Several approaches on organizing solution of inverse problems by means of PHOENICS on the basis of the technique of automated fitting are proposing. A version of a `nondestructive` method of using PHOENICS in the inverse problem solution regime and the ways of altering the program in the case of introducing optimization facilities in it are under consideration. (author) 12 refs.

Observation of isoscalar multipole strengths in exotic doubly-magic 56Ni in inelastic α scattering in inverse kinematics

Directory of Open Access Journals (Sweden)

S. Bagchi

2015-12-01

Full Text Available The Isoscalar Giant Monopole Resonance (ISGMR and the Isoscalar Giant Dipole Resonance (ISGDR compression modes have been studied in the doubly-magic unstable nucleus 56Ni. They were measured by inelastic α-particle scattering in inverse kinematics at 50 MeV/u with the MAYA active target at the GANIL facility. The centroid of the ISGMR has been obtained at Ex=19.1±0.5 MeV. Evidence for the low-lying part of the ISGDR has been found at Ex=17.4±0.7 MeV. The strength distribution for the dipole mode shows similarity with the prediction from the Hartree–Fock (HF based random-phase approximation (RPA [1]. These measurements confirm inelastic α-particle scattering as a suitable probe for exciting the ISGMR and the ISGDR modes in radioactive isotopes in inverse kinematics.

Posterior consistency for Bayesian inverse problems through stability and regression results

International Nuclear Information System (INIS)

Vollmer, Sebastian J

2013-01-01

In the Bayesian approach, the a priori knowledge about the input of a mathematical model is described via a probability measure. The joint distribution of the unknown input and the data is then conditioned, using Bayes’ formula, giving rise to the posterior distribution on the unknown input. In this setting we prove posterior consistency for nonlinear inverse problems: a sequence of data is considered, with diminishing fluctuations around a single truth and it is then of interest to show that the resulting sequence of posterior measures arising from this sequence of data concentrates around the truth used to generate the data. Posterior consistency justifies the use of the Bayesian approach very much in the same way as error bounds and convergence results for regularization techniques do. As a guiding example, we consider the inverse problem of reconstructing the diffusion coefficient from noisy observations of the solution to an elliptic PDE in divergence form. This problem is approached by splitting the forward operator into the underlying continuum model and a simpler observation operator based on the output of the model. In general, these splittings allow us to conclude posterior consistency provided a deterministic stability result for the underlying inverse problem and a posterior consistency result for the Bayesian regression problem with the push-forward prior. Moreover, we prove posterior consistency for the Bayesian regression problem based on the regularity, the tail behaviour and the small ball probabilities of the prior. (paper)

MAP estimators and their consistency in Bayesian nonparametric inverse problems

KAUST Repository

Dashti, M.

2013-09-01

We consider the inverse problem of estimating an unknown function u from noisy measurements y of a known, possibly nonlinear, map 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 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. © 2013 IOP Publishing Ltd.

MAP estimators and their consistency in Bayesian nonparametric inverse problems

International Nuclear Information System (INIS)

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

2013-01-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. (paper)

Unfolding in particle physics: A window on solving inverse problems

International Nuclear Information System (INIS)

Spano, F.

2013-01-01

Unfolding is the ensemble of techniques aimed at resolving inverse, ill-posed problems. A pedagogical introduction to the origin and main problems related to unfolding is presented and used as the the stepping stone towards the illustration of some of the most common techniques that are currently used in particle physics experiments. (authors)

Non-regularized inversion method from light scattering applied to ferrofluid magnetization curves for magnetic size distribution analysis

International Nuclear Information System (INIS)

Rijssel, Jos van; Kuipers, Bonny W.M.; Erné, Ben H.

2014-01-01

A numerical inversion method known from the analysis of light scattering by colloidal dispersions is now applied to magnetization curves of ferrofluids. The distribution of magnetic particle sizes or dipole moments is determined without assuming that the distribution is unimodal or of a particular shape. The inversion method enforces positive number densities via a non-negative least squares procedure. It is tested successfully on experimental and simulated data for ferrofluid samples with known multimodal size distributions. The created computer program MINORIM is made available on the web. - Highlights: • A method from light scattering is applied to analyze ferrofluid magnetization curves. • A magnetic size distribution is obtained without prior assumption of its shape. • The method is tested successfully on ferrofluids with a known size distribution. • The practical limits of the method are explored with simulated data including noise. • This method is implemented in the program MINORIM, freely available online

Canonical quantization of gravity and a problem of scattering

International Nuclear Information System (INIS)

Rubakov, V.A.

1980-01-01

Linearized theory of gravity is quantized both in a naive way and as a proper limit of the Dirac-Wheeler-De Witt approach to the quantization of the full theory. The equivalence between the two approaches is established. The problem of scattering in the canonically quantized theory of gravitation is investigated. The concept of the background metric naturally appears in the canonical formalism for this case. The equivalence between canonical and path-integral approaches is established for the problem of scattering. Some kinetical properties of functionals in Wheeler superspace are studied in an appendix. (author)

An inverse hyperbolic heat conduction problem in estimating surface heat flux by the conjugate gradient method

International Nuclear Information System (INIS)

Huang, C.-H.; Wu, H.-H.

2006-01-01

In the present study an inverse hyperbolic heat conduction problem is solved by the conjugate gradient method (CGM) in estimating the unknown boundary heat flux based on the boundary temperature measurements. Results obtained in this inverse problem will be justified based on the numerical experiments where three different heat flux distributions are to be determined. Results show that the inverse solutions can always be obtained with any arbitrary initial guesses of the boundary heat flux. Moreover, the drawbacks of the previous study for this similar inverse problem, such as (1) the inverse solution has phase error and (2) the inverse solution is sensitive to measurement error, can be avoided in the present algorithm. Finally, it is concluded that accurate boundary heat flux can be estimated in this study

Inverse problem of radiofrequency sounding of ionosphere

Science.gov (United States)

Velichko, E. N.; Yu. Grishentsev, A.; Korobeynikov, A. G.

2016-01-01

An algorithm for the solution of the inverse problem of vertical ionosphere sounding and a mathematical model of noise filtering are presented. An automated system for processing and analysis of spectrograms of vertical ionosphere sounding based on our algorithm is described. It is shown that the algorithm we suggest has a rather high efficiency. This is supported by the data obtained at the ionospheric stations of the so-called “AIS-M” type.

An inverse problem in a parabolic equation

Directory of Open Access Journals (Sweden)

Zhilin Li

1998-11-01

Full Text Available In this paper, an inverse problem in a parabolic equation is studied. An unknown function in the equation is related to two integral equations in terms of heat kernel. One of the integral equations is well-posed while another is ill-posed. A regularization approach for constructing an approximate solution to the ill-posed integral equation is proposed. Theoretical analysis and numerical experiment are provided to support the method.

An inverse method for radiation transport

Energy Technology Data Exchange (ETDEWEB)

Favorite, J. A. (Jeffrey A.); Sanchez, R. (Richard)

2004-01-01

Adjoint functions have been used with forward functions to compute gradients in implicit (iterative) solution methods for inverse problems in optical tomography, geoscience, thermal science, and other fields, but only once has this approach been used for inverse solutions to the Boltzmann transport equation. In this paper, this approach is used to develop an inverse method that requires only angle-independent flux measurements, rather than angle-dependent measurements as was done previously. The method is applied to a simplified form of the transport equation that does not include scattering. The resulting procedure uses measured values of gamma-ray fluxes of discrete, characteristic energies to determine interface locations in a multilayer shield. The method was implemented with a Newton-Raphson optimization algorithm, and it worked very well in numerical one-dimensional spherical test cases. A more sophisticated optimization method would better exploit the potential of the inverse method.

Inverse Problems and Uncertainty Quantification

KAUST Repository

Litvinenko, Alexander

2014-01-06

In a Bayesian setting, inverse problems and uncertainty quantification (UQ) - the propagation of uncertainty through a computational (forward) modelare strongly connected. In the form of conditional expectation the Bayesian update becomes computationally attractive. This is especially the case as together with a functional or spectral approach for the forward UQ there is no need for time- consuming and slowly convergent Monte Carlo sampling. The developed sampling- free non-linear Bayesian update is derived from the variational problem associated with conditional expectation. This formulation in general calls for further discretisa- tion to make the computation possible, and we choose a polynomial approximation. After giving details on the actual computation in the framework of functional or spectral approximations, we demonstrate the workings of the algorithm on a number of examples of increasing complexity. At last, we compare the linear and quadratic Bayesian update on the small but taxing example of the chaotic Lorenz 84 model, where we experiment with the influence of different observation or measurement operators on the update.

Inverse Problems and Uncertainty Quantification

KAUST Repository

Litvinenko, Alexander; Matthies, Hermann G.

2014-01-01

In a Bayesian setting, inverse problems and uncertainty quantification (UQ) - the propagation of uncertainty through a computational (forward) modelare strongly connected. In the form of conditional expectation the Bayesian update becomes computationally attractive. This is especially the case as together with a functional or spectral approach for the forward UQ there is no need for time- consuming and slowly convergent Monte Carlo sampling. The developed sampling- free non-linear Bayesian update is derived from the variational problem associated with conditional expectation. This formulation in general calls for further discretisa- tion to make the computation possible, and we choose a polynomial approximation. After giving details on the actual computation in the framework of functional or spectral approximations, we demonstrate the workings of the algorithm on a number of examples of increasing complexity. At last, we compare the linear and quadratic Bayesian update on the small but taxing example of the chaotic Lorenz 84 model, where we experiment with the influence of different observation or measurement operators on the update.

Inverse problems and uncertainty quantification

KAUST Repository

Litvinenko, Alexander

2013-12-18

In a Bayesian setting, inverse problems and uncertainty quantification (UQ)— the propagation of uncertainty through a computational (forward) model—are strongly connected. In the form of conditional expectation the Bayesian update becomes computationally attractive. This is especially the case as together with a functional or spectral approach for the forward UQ there is no need for time- consuming and slowly convergent Monte Carlo sampling. The developed sampling- free non-linear Bayesian update is derived from the variational problem associated with conditional expectation. This formulation in general calls for further discretisa- tion to make the computation possible, and we choose a polynomial approximation. After giving details on the actual computation in the framework of functional or spectral approximations, we demonstrate the workings of the algorithm on a number of examples of increasing complexity. At last, we compare the linear and quadratic Bayesian update on the small but taxing example of the chaotic Lorenz 84 model, where we experiment with the influence of different observation or measurement operators on the update.

Software tool for resolution of inverse problems using artificial intelligence techniques: an application in neutron spectrometry

International Nuclear Information System (INIS)

Castaneda M, V. H.; Martinez B, M. R.; Solis S, L. O.; Castaneda M, R.; Leon P, A. A.; Hernandez P, C. F.; Espinoza G, J. G.; Ortiz R, J. M.; Vega C, H. R.; Mendez, R.; Gallego, E.; Sousa L, M. A.

2016-10-01

The Taguchi methodology has proved to be highly efficient to solve inverse problems, in which the values of some parameters of the model must be obtained from the observed data. There are intrinsic mathematical characteristics that make a problem known as inverse. Inverse problems appear in many branches of science, engineering and mathematics. To solve this type of problem, researches have used different techniques. Recently, the use of techniques based on Artificial Intelligence technology is being explored by researches. This paper presents the use of a software tool based on artificial neural networks of generalized regression in the solution of inverse problems with application in high energy physics, specifically in the solution of the problem of neutron spectrometry. To solve this problem we use a software tool developed in the Mat Lab programming environment, which employs a friendly user interface, intuitive and easy to use for the user. This computational tool solves the inverse problem involved in the reconstruction of the neutron spectrum based on measurements made with a Bonner spheres spectrometric system. Introducing this information, the neural network is able to reconstruct the neutron spectrum with high performance and generalization capability. The tool allows that the end user does not require great training or technical knowledge in development and/or use of software, so it facilitates the use of the program for the resolution of inverse problems that are in several areas of knowledge. The techniques of Artificial Intelligence present singular veracity to solve inverse problems, given the characteristics of artificial neural networks and their network topology, therefore, the tool developed has been very useful, since the results generated by the Artificial Neural Network require few time in comparison to other techniques and are correct results comparing them with the actual data of the experiment. (Author)

Dependence of the forward light scattering on the refractive index of particles

Science.gov (United States)

Guo, Lufang; Shen, Jianqi

2018-05-01

In particle sizing technique based on forward light scattering, the scattered light signal (SLS) is closely related to the relative refractive index (RRI) of the particles to the surrounding, especially when the particles are transparent (or weakly absorbent) and the particles are small in size. The interference between the diffraction (Diff) and the multiple internal reflections (MIR) of scattered light can lead to the oscillation of the SLS on RRI and the abnormal intervals, especially for narrowly-distributed small particle systems. This makes the inverse problem more difficult. In order to improve the inverse results, Tikhonov regularization algorithm with B-spline functions is proposed, in which the matrix element is calculated for a range of particle sizes instead using the mean particle diameter of size fractions. In this way, the influence of abnormal intervals on the inverse results can be eliminated. In addition, for measurements on narrowly distributed small particles, it is suggested to detect the SLS in a wider scattering angle to include more information.

Renormalized nonlinear sensitivity kernel and inverse thin-slab propagator in T-matrix formalism for wave-equation tomography

International Nuclear Information System (INIS)

Wu, Ru-Shan; Wang, Benfeng; Hu, Chunhua

2015-01-01

We derived the renormalized nonlinear sensitivity operator and the related inverse thin-slab propagator (ITSP) for nonlinear tomographic waveform inversion based on the theory of nonlinear partial derivative operator and its De Wolf approximation. The inverse propagator is based on a renormalization procedure to the forward and inverse transition matrix scattering series. The ITSP eliminates the divergence of the inverse Born series for strong perturbations by stepwise partial summation (renormalization). Numerical tests showed that the inverse Born T-series starts to diverge at moderate perturbation (20% for the given model of Gaussian ball with a radius of 5 wavelength), while the ITSP has no divergence problem for any strong perturbations (up to 100% perturbation for test model). In addition, the ITSP is a non-iterative, marching algorithm with only one sweep, and therefore very efficient in comparison with the iterative inversion based on the inverse-Born scattering series. This convergence and efficiency improvement has potential applications to the iterative procedure of waveform inversion. (paper)

Inverse problem of the vibrational band gap of periodically supported beam

Science.gov (United States)

Shi, Xiaona; Shu, Haisheng; Dong, Fuzhen; Zhao, Lei

2017-04-01

The researches of periodic structures have a long history with the main contents confined in the field of forward problem. In this paper, the inverse problem is considered and an overall frame is proposed which includes two main stages, i.e., the band gap criterion and its optimization. As a preliminary investigation, the inverse problem of the flexural vibrational band gap of a periodically supported beam is analyzed. According to existing knowledge of its forward problem, the band gap criterion is given in implicit form. Then, two cases with three independent parameters, namely the double supported case and the triple one, are studied in detail and the explicit expressions of the feasible domain are constructed by numerical fitting. Finally, the parameter optimization of the double supported case with three variables is conducted using genetic algorithm aiming for the best mean attenuation within specified frequency band.

Methods and Algorithms for Solving Inverse Problems for Fractional Advection-Dispersion Equations

KAUST Repository

Aldoghaither, Abeer

2015-11-12

Fractional calculus has been introduced as an e cient tool for modeling physical phenomena, thanks to its memory and hereditary properties. For example, fractional models have been successfully used to describe anomalous di↵usion processes such as contaminant transport in soil, oil flow in porous media, and groundwater flow. These models capture important features of particle transport such as particles with velocity variations and long-rest periods. Mathematical modeling of physical phenomena requires the identification of pa- rameters and variables from available measurements. This is referred to as an inverse problem. In this work, we are interested in studying theoretically and numerically inverse problems for space Fractional Advection-Dispersion Equation (FADE), which is used to model solute transport in porous media. Identifying parameters for such an equa- tion is important to understand how chemical or biological contaminants are trans- ported throughout surface aquifer systems. For instance, an estimate of the di↵eren- tiation order in groundwater contaminant transport model can provide information about soil properties, such as the heterogeneity of the medium. Our main contribution is to propose a novel e cient algorithm based on modulat-ing functions to estimate the coe cients and the di↵erentiation order for space FADE, which can be extended to general fractional Partial Di↵erential Equation (PDE). We also show how the method can be applied to the source inverse problem. This work is divided into two parts: In part I, the proposed method is described and studied through an extensive numerical analysis. The local convergence of the proposed two-stage algorithm is proven for 1D space FADE. The properties of this method are studied along with its limitations. Then, the algorithm is generalized to the 2D FADE. In part II, we analyze direct and inverse source problems for a space FADE. The problem consists of recovering the source term using final

Application of the kernel method to the inverse geosounding problem.

Science.gov (United States)

Hidalgo, Hugo; Sosa León, Sonia; Gómez-Treviño, Enrique

2003-01-01

Determining the layered structure of the earth demands the solution of a variety of inverse problems; in the case of electromagnetic soundings at low induction numbers, the problem is linear, for the measurements may be represented as a linear functional of the electrical conductivity distribution. In this paper, an application of the support vector (SV) regression technique to the inversion of electromagnetic data is presented. We take advantage of the regularizing properties of the SV learning algorithm and use it as a modeling technique with synthetic and field data. The SV method presents better recovery of synthetic models than Tikhonov's regularization. As the SV formulation is solved in the space of the data, which has a small dimension in this application, a smaller problem than that considered with Tikhonov's regularization is produced. For field data, the SV formulation develops models similar to those obtained via linear programming techniques, but with the added characteristic of robustness.

Two-Dimensional Linear Inversion of GPR Data with a Shifting Zoom along the Observation Line

Directory of Open Access Journals (Sweden)

Raffaele Persico

2017-09-01

Full Text Available Linear inverse scattering problems can be solved by regularized inversion of a matrix, whose calculation and inversion may require significant computing resources, in particular, a significant amount of RAM memory. This effort is dependent on the extent of the investigation domain, which drives a large amount of data to be gathered and a large number of unknowns to be looked for, when this domain becomes electrically large. This leads, in turn, to the problem of inversion of excessively large matrices. Here, we consider the problem of a ground-penetrating radar (GPR survey in two-dimensional (2D geometry, with antennas at an electrically short distance from the soil. In particular, we present a strategy to afford inversion of large investigation domains, based on a shifting zoom procedure. The proposed strategy was successfully validated using experimental radar data.

Integral equations of the first kind, inverse problems and regularization: a crash course

International Nuclear Information System (INIS)

Groetsch, C W

2007-01-01

This paper is an expository survey of the basic theory of regularization for Fredholm integral equations of the first kind and related background material on inverse problems. We begin with an historical introduction to the field of integral equations of the first kind, with special emphasis on model inverse problems that lead to such equations. The basic theory of linear Fredholm equations of the first kind, paying particular attention to E. Schmidt's singular function analysis, Picard's existence criterion, and the Moore-Penrose theory of generalized inverses is outlined. The fundamentals of the theory of Tikhonov regularization are then treated and a collection of exercises and a bibliography are provided

Parameterization analysis and inversion for orthorhombic media

KAUST Repository

Masmoudi, Nabil

2018-05-01

Accounting for azimuthal anisotropy is necessary for the processing and inversion of wide-azimuth and wide-aperture seismic data because wave speeds naturally depend on the wave propagation direction. Orthorhombic anisotropy is considered the most effective anisotropic model that approximates the azimuthal anisotropy we observe in seismic data. In the framework of full wave form inversion (FWI), the large number of parameters describing orthorhombic media exerts a considerable trade-off and increases the non-linearity of the inversion problem. Choosing a suitable parameterization for the model, and identifying which parameters in that parameterization could be well resolved, are essential to a successful inversion. In this thesis, I derive the radiation patterns for different acoustic orthorhombic parameterization. Analyzing the angular dependence of the scattering of the parameters of different parameterizations starting with the conventionally used notation, I assess the potential trade-off between the parameters and the resolution in describing the data and inverting for the parameters. In order to build practical inversion strategies, I suggest new parameters (called deviation parameters) for a new parameterization style in orthorhombic media. The novel parameters denoted ∈d, ƞd and δd are dimensionless and represent a measure of deviation between the vertical planes in orthorhombic anisotropy. The main feature of the deviation parameters consists of keeping the scattering of the vertical transversely isotropic (VTI) parameters stationary with azimuth. Using these scattering features, we can condition FWI to invert for the parameters which the data are sensitive to, at different stages, scales, and locations in the model. With this parameterization, the data are mainly sensitive to the scattering of 3 parameters (out of six that describe an acoustic orthorhombic medium): the horizontal velocity in the x1 direction, ∈1 which provides scattering mainly near

Quantum trajectories in complex space: One-dimensional stationary scattering problems

International Nuclear Information System (INIS)

Chou, C.-C.; Wyatt, Robert E.

2008-01-01

One-dimensional time-independent scattering problems are investigated in the framework of the quantum Hamilton-Jacobi formalism. The equation for the local approximate quantum trajectories near the stagnation point of the quantum momentum function is derived, and the first derivative of the quantum momentum function is related to the local structure of quantum trajectories. Exact complex quantum trajectories are determined for two examples by numerically integrating the equations of motion. For the soft potential step, some particles penetrate into the nonclassical region, and then turn back to the reflection region. For the barrier scattering problem, quantum trajectories may spiral into the attractors or from the repellers in the barrier region. Although the classical potentials extended to complex space show different pole structures for each problem, the quantum potentials present the same second-order pole structure in the reflection region. This paper not only analyzes complex quantum trajectories and the total potentials for these examples but also demonstrates general properties and similar structures of the complex quantum trajectories and the quantum potentials for one-dimensional time-independent scattering problems

Data-Driven Model Order Reduction for Bayesian Inverse Problems

KAUST Repository

Cui, Tiangang

2014-01-06

One of the major challenges in using MCMC for the solution of inverse problems is the repeated evaluation of computationally expensive numerical models. We develop a data-driven projection- based model order reduction technique to reduce the computational cost of numerical PDE evaluations in this context.

Variational approach to direct and inverse problems of atmospheric pollution studies

Science.gov (United States)

Penenko, Vladimir; Tsvetova, Elena; Penenko, Alexey

2016-04-01

We present the development of a variational approach for solving interrelated problems of atmospheric hydrodynamics and chemistry concerning air pollution transport and transformations. The proposed approach allows us to carry out complex studies of different-scale physical and chemical processes using the methods of direct and inverse modeling [1-3]. We formulate the problems of risk/vulnerability and uncertainty assessment, sensitivity studies, variational data assimilation procedures [4], etc. A computational technology of constructing consistent mathematical models and methods of their numerical implementation is based on the variational principle in the weak constraint formulation specifically designed to account for uncertainties in models and observations. Algorithms for direct and inverse modeling are designed with the use of global and local adjoint problems. Implementing the idea of adjoint integrating factors provides unconditionally monotone and stable discrete-analytic approximations for convection-diffusion-reaction problems [5,6]. The general framework is applied to the direct and inverse problems for the models of transport and transformation of pollutants in Siberian and Arctic regions. The work has been partially supported by the RFBR grant 14-01-00125 and RAS Presidium Program I.33P. References: 1. V. Penenko, A.Baklanov, E. Tsvetova and A. Mahura . Direct and inverse problems in a variational concept of environmental modeling //Pure and Applied Geoph.(2012) v.169: 447-465. 2. V. V. Penenko, E. A. Tsvetova, and A. V. Penenko Development of variational approach for direct and inverse problems of atmospheric hydrodynamics and chemistry, Izvestiya, Atmospheric and Oceanic Physics, 2015, Vol. 51, No. 3, p. 311-319, DOI: 10.1134/S0001433815030093. 3. V.V. Penenko, E.A. Tsvetova, A.V. Penenko. Methods based on the joint use of models and observational data in the framework of variational approach to forecasting weather and atmospheric composition

Inverse problems for partial differential equations

CERN Document Server

Isakov, Victor

2017-01-01

This third edition expands upon the earlier edition by adding nearly 40 pages of new material reflecting the analytical and numerical progress in inverse problems in last 10 years. As in the second edition, the emphasis is on new ideas and methods rather than technical improvements. These new ideas include use of the stationary phase method in the two-dimensional elliptic problems and of multi frequencies\\temporal data to improve stability and numerical resolution. There are also numerous corrections and improvements of the exposition throughout. This book is intended for mathematicians working with partial differential equations and their applications, physicists, geophysicists, and financial, electrical, and mechanical engineers involved with nondestructive evaluation, seismic exploration, remote sensing, and various kinds of tomography. Review of the second edition: "The first edition of this excellent book appeared in 1998 and became a standard reference for everyone interested in analysis and numerics of...

Inverse Problems in Geosciences: Modelling the Rock Properties of an Oil Reservoir

DEFF Research Database (Denmark)

Lange, Katrine

. We have developed and implemented the Frequency Matching method that uses the closed form expression of the a priori probability density function to formulate an inverse problem and compute the maximum a posteriori solution to it. Other methods for computing models that simultaneously fit data...... of the subsurface of the reservoirs. Hence the focus of this work has been on acquiring models of spatial parameters describing rock properties of the subsurface using geostatistical a priori knowledge and available geophysical data. Such models are solutions to often severely under-determined, inverse problems...

A variational Bayesian method to inverse problems with impulsive noise

KAUST Repository

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. © 2011 Elsevier Inc.

a Proposed Benchmark Problem for Scatter Calculations in Radiographic Modelling

Science.gov (United States)

Jaenisch, G.-R.; Bellon, C.; Schumm, A.; Tabary, J.; Duvauchelle, Ph.

2009-03-01

Code Validation is a permanent concern in computer modelling, and has been addressed repeatedly in eddy current and ultrasonic modeling. A good benchmark problem is sufficiently simple to be taken into account by various codes without strong requirements on geometry representation capabilities, focuses on few or even a single aspect of the problem at hand to facilitate interpretation and to avoid that compound errors compensate themselves, yields a quantitative result and is experimentally accessible. In this paper we attempt to address code validation for one aspect of radiographic modeling, the scattered radiation prediction. Many NDT applications can not neglect scattered radiation, and the scatter calculation thus is important to faithfully simulate the inspection situation. Our benchmark problem covers the wall thickness range of 10 to 50 mm for single wall inspections, with energies ranging from 100 to 500 keV in the first stage, and up to 1 MeV with wall thicknesses up to 70 mm in the extended stage. A simple plate geometry is sufficient for this purpose, and the scatter data is compared on a photon level, without a film model, which allows for comparisons with reference codes like MCNP. We compare results of three Monte Carlo codes (McRay, Sindbad and Moderato) as well as an analytical first order scattering code (VXI), and confront them to results obtained with MCNP. The comparison with an analytical scatter model provides insights into the application domain where this kind of approach can successfully replace Monte-Carlo calculations.

Solving inverse problems for biological models using the collage method for differential equations.

Science.gov (United States)

Capasso, V; Kunze, H E; La Torre, D; Vrscay, E R

2013-07-01

In the first part of this paper we show how inverse problems for differential equations can be solved using the so-called collage method. Inverse problems can be solved by minimizing the collage distance in an appropriate metric space. We then provide several numerical examples in mathematical biology. We consider applications of this approach to the following areas: population dynamics, mRNA and protein concentration, bacteria and amoeba cells interaction, tumor growth.

SCIENTIFIC AND METHODICAL ASPECTS OF FORMATION OF SUBJECT CONTENT OF TRAINING COURSESFOR INVERSE PROBLEMS FOR DIFFERENTIAL EQUATIONS

Directory of Open Access Journals (Sweden)

В С Корнилов

2016-12-01

Full Text Available The article presents scientific and methodical aspects of forming the content of education inverse problems for differential equations for students of higher educational institutions of physical, mathematical and natural science training areas. The goals are formulated and the principles of training are the content of learning inverse problems for differential equations. Attention is drawn to the particular issues of teaching courses inverse problems. Describes the classification criteria and target modules that play the role of tools to create and analyze the model and curriculum, forming learning content inverse problems for differential equations. The content classification features and target modules. Formulate conclusions that learning the inverse problems for differential equations has scientific, educational and humanitarian potential of students and as a result of this training they gain the fundamental knowledge in the applied and computational mathematics, and also develop scientific worldview, applied, environmental, information thinking.

Iterative Reconstruction Methods for Hybrid Inverse Problems in Impedance Tomography

DEFF Research Database (Denmark)

Hoffmann, Kristoffer; Knudsen, Kim

2014-01-01

For a general formulation of hybrid inverse problems in impedance tomography the Picard and Newton iterative schemes are adapted and four iterative reconstruction algorithms are developed. The general problem formulation includes several existing hybrid imaging modalities such as current density...... impedance imaging, magnetic resonance electrical impedance tomography, and ultrasound modulated electrical impedance tomography, and the unified approach to the reconstruction problem encompasses several algorithms suggested in the literature. The four proposed algorithms are implemented numerically in two...

Toward precise solution of one-dimensional velocity inverse problems

International Nuclear Information System (INIS)

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

Numerical approach to the inverse convection-diffusion problem

International Nuclear Information System (INIS)

Yang, X-H; She, D-X; Li, J-Q

2008-01-01

In this paper, the inverse problem on source term identification in convection-diffusion equation is transformed into an optimization problem. To reduce the computational cost and improve computational accuracy for the optimization problem, a new algorithm, chaos real-coded hybrid-accelerating evolution algorithm (CRHAEA), is proposed, in which an initial population is generated by chaos mapping, and new chaos mutation and simplex evolution operation are used. With the shrinking of searching range, CRHAEA gradually directs to an optimal result with the excellent individuals obtained by real-coded evolution algorithm. Its convergence is analyzed. Its efficiency is demonstrated by 15 test functions. Numerical simulation shows that CRHAEA has some advantages over the real-coded accelerated evolution algorithm, the chaos algorithm and the pure random search algorithm

Stochastic inverse problems: Models and metrics

International Nuclear Information System (INIS)

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

2015-01-01

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

Stochastic inverse problems: Models and metrics

Science.gov (United States)

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

2015-03-01

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.

From inverse problems to learning: a Statistical Mechanics approach

Science.gov (United States)

Baldassi, Carlo; Gerace, Federica; Saglietti, Luca; Zecchina, Riccardo

2018-01-01

We present a brief introduction to the statistical mechanics approaches for the study of inverse problems in data science. We then provide concrete new results on inferring couplings from sampled configurations in systems characterized by an extensive number of stable attractors in the low temperature regime. We also show how these result are connected to the problem of learning with realistic weak signals in computational neuroscience. Our techniques and algorithms rely on advanced mean-field methods developed in the context of disordered systems.

An inverse Sturm–Liouville problem with a fractional derivative

KAUST Repository

Jin, Bangti

2012-05-01

In this paper, we numerically investigate an inverse problem of recovering the potential term in a fractional Sturm-Liouville problem from one spectrum. The qualitative behaviors of the eigenvalues and eigenfunctions are discussed, and numerical reconstructions of the potential with a Newton method from finite spectral data are presented. Surprisingly, it allows very satisfactory reconstructions for both smooth and discontinuous potentials, provided that the order . α∈. (1,. 2) of fractional derivative is sufficiently away from 2. © 2012 Elsevier Inc.

A nonlinear inversion for the velocity background and perturbation models

KAUST Repository

Wu, Zedong

2015-08-19

Reflected waveform inversion (RWI) provides a method to reduce the nonlinearity of the standard full waveform inversion (FWI) by inverting for the single scattered wavefield obtained using an image. However, current RWI methods usually neglect diving waves, which is an important source of information for extracting the long wavelength components of the velocity model. Thus, we propose a new optimization problem through breaking the velocity model into the background and the perturbation in the wave equation directly. In this case, the perturbed model is no longer the single scattering model, but includes all scattering. We optimize both components simultaneously, and thus, the objective function is nonlinear with respect to both the background and perturbation. The new introduced w can absorb the non-smooth update of background naturally. Application to the Marmousi model with frequencies that start at 5 Hz shows that this method can converge to the accurate velocity starting from a linearly increasing initial velocity. Application to the SEG2014 demonstrates the versatility of the approach.

Comparison of four software packages applied to a scattering problem

DEFF Research Database (Denmark)

Albertsen, Niels Christian; Chesneaux, Jean-Marie; Christiansen, Søren

1999-01-01

We investigate characteristic features of four different software packages by applying them to the numerical solution of a non-trivial physical problem in computer simulation, viz., scattering of waves from a sinusoidal boundary. The numerical method used is based on boundary collocation. This le......We investigate characteristic features of four different software packages by applying them to the numerical solution of a non-trivial physical problem in computer simulation, viz., scattering of waves from a sinusoidal boundary. The numerical method used is based on boundary collocation...

ITOUGH2: Solving TOUGH inverse problems

Energy Technology Data Exchange (ETDEWEB)

Finsterle, S.; Pruess, K. [Lawrence Berkeley Laboratory, CA (United States)

1995-03-01

ITOUGH2 is a program that provides inverse modeling capabilities for the TOUGH2 code. While the main purpose of ITOUGH2 is to estimate two-phase hydraulic properties of calibrating a TOUGH2 model to laboratory or field data, the information obtained by evaluating parameter sensitivities can also be used to optimize the design of an experiment, and to analyze the uncertainty of model predictions. ITOUGH2 has been applied to a number of laboratory and field experiments on different scales. Three examples are discussed in this paper, demonstrating the code`s capability to support test design, data analysis, and model predictions for a variety of TOUGH problems.

Inverse Problems in a Bayesian Setting

KAUST Repository

Matthies, Hermann G.

2016-02-13

In a Bayesian setting, inverse problems and uncertainty quantification (UQ)—the propagation of uncertainty through a computational (forward) model—are strongly connected. In the form of conditional expectation the Bayesian update becomes computationally attractive. We give a detailed account of this approach via conditional approximation, various approximations, and the construction of filters. Together with a functional or spectral approach for the forward UQ there is no need for time-consuming and slowly convergent Monte Carlo sampling. The developed sampling-free non-linear Bayesian update in form of a filter is derived from the variational problem associated with conditional expectation. This formulation in general calls for further discretisation to make the computation possible, and we choose a polynomial approximation. After giving details on the actual computation in the framework of functional or spectral approximations, we demonstrate the workings of the algorithm on a number of examples of increasing complexity. At last, we compare the linear and nonlinear Bayesian update in form of a filter on some examples.

Inverse Problems in a Bayesian Setting

KAUST Repository

Matthies, Hermann G.; Zander, Elmar; Rosić, Bojana V.; Litvinenko, Alexander; Pajonk, Oliver

2016-01-01

In a Bayesian setting, inverse problems and uncertainty quantification (UQ)—the propagation of uncertainty through a computational (forward) model—are strongly connected. In the form of conditional expectation the Bayesian update becomes computationally attractive. We give a detailed account of this approach via conditional approximation, various approximations, and the construction of filters. Together with a functional or spectral approach for the forward UQ there is no need for time-consuming and slowly convergent Monte Carlo sampling. The developed sampling-free non-linear Bayesian update in form of a filter is derived from the variational problem associated with conditional expectation. This formulation in general calls for further discretisation to make the computation possible, and we choose a polynomial approximation. After giving details on the actual computation in the framework of functional or spectral approximations, we demonstrate the workings of the algorithm on a number of examples of increasing complexity. At last, we compare the linear and nonlinear Bayesian update in form of a filter on some examples.

A mathematical framework for inverse wave problems in heterogeneous media

NARCIS (Netherlands)

Blazek, K.D.; Stolk, C.; Symes, W.W.

2013-01-01

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

Methane combustion kinetic rate constants determination: an ill-posed inverse problem analysis

Directory of Open Access Journals (Sweden)

Bárbara D. L. Ferreira

2013-01-01

Full Text Available Methane combustion was studied by the Westbrook and Dryer model. This well-established simplified mechanism is very useful in combustion science, for computational effort can be notably reduced. In the inversion procedure to be studied, rate constants are obtained from [CO] concentration data. However, when inherent experimental errors in chemical concentrations are considered, an ill-conditioned inverse problem must be solved for which appropriate mathematical algorithms are needed. A recurrent neural network was chosen due to its numerical stability and robustness. The proposed methodology was compared against Simplex and Levenberg-Marquardt, the most used methods for optimization problems.

Some surprises and paradoxes revealed by inverse problem approach and notion about qualitative solutions of Schroedinger equations 'in mind'

International Nuclear Information System (INIS)

Zakhariev, B.N.; Chabanov, V.M.

2008-01-01

It was an important examination to give a review talk at the previous Conference on Inverse Quantum Scattering (1996, Lake Balaton) about computer visualization of this science in front of its fathers - creators, B.M. Levitan and V.A. Marchenko. We have achieved a new understanding that the discovered main rules of transformations of a single wave function bump, e.g., for the ground bound states of one dimensional quantum systems are applicable to any state of any potential with arbitrary number of bumps from finite to unlimited ones as scattering states and bound states embedded into continuum. It appeared that we need only to repeat the rule mentally the necessary number of times. That uttermost simplification and unification of physical notion of spectral, scattering and decay control for any potential have got an obligatory praise from B.M. Levitan at the conference and was a mighty stimulus for our further research. After that we have written both Russian (2002) and improved English editions of 'Submissive Quantum Mechanics. New Status of the Theory in Inverse Problem Approach' (appeared at the very end of 2007). This book was written for correction of the present defect in quantum education throughout the world. Recently the quantum IP intuition helped us to discover a new concept of permanent wave resonance with potential spatial oscillations. This means the constant wave swinging frequency on the whole energy intervals of spectral forbidden zones destroying physical solutions and deepening the theory of waves in periodic potentials. It also shows the other side of strengthening the fundamentally important magic structures. A 'new language' of wave bending will be presented to enrich our quantum intuition, e.g., the paradoxical effective attraction of barriers and repulsion of wells in multichannel systems, etc. (author)

Multiple scattering problems in heavy ion elastic recoil detection analysis

International Nuclear Information System (INIS)

Johnston, P.N.; El Bouanani, M.; Stannard, W.B.; Bubb, I.F.; Cohen, D.D.; Dytlewski, N.; Siegele, R.

1998-01-01

A number of groups use Heavy Ion Elastic Recoil Detection Analysis (HIERDA) to study materials science problems. Nevertheless, there is no standard methodology for the analysis of HIERDA spectra. To overcome this deficiency we have been establishing codes for 2-dimensional data analysis. A major problem involves the effects of multiple and plural scattering which are very significant, even for quite thin (∼100 nm) layers of the very heavy elements. To examine the effects of multiple scattering we have made comparisons between the small-angle model of Sigmund et al. and TRIM calculations. (authors)

A variational Bayesian method to inverse problems with impulsive noise

KAUST Repository

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

Inverse problem for the mean-field monomer-dimer model with attractive interaction

International Nuclear Information System (INIS)

Contucci, Pierluigi; Luzi, Rachele; Vernia, Cecilia

2017-01-01

The inverse problem method is tested for a class of monomer-dimer statistical mechanics models that contain also an attractive potential and display a mean-field critical point at a boundary of a coexistence line. The inversion is obtained by analytically identifying the parameters in terms of the correlation functions and via the maximum-likelihood method. The precision is tested in the whole phase space and, when close to the coexistence line, the algorithm is used together with a clustering method to take care of the underlying possible ambiguity of the inversion. (paper)

NON-INVASIVE INVERSE PROBLEM IN CIVIL ENGINEERING

Directory of Open Access Journals (Sweden)

Jan Havelka

2017-11-01

Full Text Available In this contribution we focus on recovery of spatial distribution of material parameters utilizing only non-invasive boundary measurements. Such methods has gained its importance as imaging techniques in medicine, geophysics or archaeology. We apply similar principles for non-stationary heat transfer in civil engineering. In oppose to standard technique which rely on external loading devices, we assume the natural fluctuation of temperature throughout day and night can provide sufficient information to recover the underlying material parameters. The inverse problem was solved by a modified regularised Gauss-Newton iterative scheme and the underlying forward problem is solved with a finite element space-time discretisation. We show a successful reconstruction of material parameters on a synthetic example with real measurements. The virtual experiment also reveals the insensitivity to practical precision of sensor measurements.

Oblique projections and standard-form transformations for discrete inverse problems

DEFF Research Database (Denmark)

Hansen, Per Christian

2013-01-01

This tutorial paper considers a specific computational tool for the numerical solution of discrete inverse problems, known as the standard-form transformation, by which we can treat general Tikhonov regularization problems efficiently. In the tradition of B. N. Datta's expositions of numerical li...... linear algebra, we use the close relationship between oblique projections, pseudoinverses, and matrix computations to derive a simple geometric motivation and algebraic formulation of the standard-form transformation....

Two numerical methods for an inverse problem for the 2-D Helmholtz equation

CERN Document Server

Gryazin, Y A; Lucas, T R

2003-01-01

Two solution methods for the inverse problem for the 2-D Helmholtz equation are developed, tested, and compared. The proposed approaches are based on a marching finite-difference scheme which requires the solution of an overdetermined system at each step. The preconditioned conjugate gradient method is used for rapid solutions of these systems and an efficient preconditioner has been developed for this class of problems. Underlying target applications include the imaging of land mines, unexploded ordinance, and pollutant plumes in environmental cleanup sites, each formulated as an inverse problem for a 2-D Helmholtz equation. The images represent the electromagnetic properties of the respective underground regions. Extensive numerical results are presented.

Use of Genetic Algorithms to solve Inverse Problems in Relativistic Hydrodynamics

Science.gov (United States)

Guzmán, F. S.; González, J. A.

2018-04-01

We present the use of Genetic Algorithms (GAs) as a strategy to solve inverse problems associated with models of relativistic hydrodynamics. The signal we consider to emulate an observation is the density of a relativistic gas, measured at a point where a shock is traveling. This shock is generated numerically out of a Riemann problem with mildly relativistic conditions. The inverse problem we propose is the prediction of the initial conditions of density, velocity and pressure of the Riemann problem that gave origin to that signal. For this we use the density, velocity and pressure of the gas at both sides of the discontinuity, as the six genes of an organism, initially with random values within a tolerance. We then prepare an initial population of N of these organisms and evolve them using methods based on GAs. In the end, the organism with the best fitness of each generation is compared to the signal and the process ends when the set of initial conditions of the organisms of a later generation fit the Signal within a tolerance.

On the connection between the inverse transform method and the exact quantum eigenstates

International Nuclear Information System (INIS)

Honerkamp, J.; Weber, P.; Wiesler, A.

1979-01-01

The 'inverse scattering transformation', which has been used to solve certain nonlinear field theories classically, is discussed in the context of the quantized version of these theories. In particular the non-linear Schroedinger equation and the massive Thirring model are considered. It is found that certain Jost functions of the associated scattering problem lead already, in quantizing the theory, to creation operators for the exact eigenstates of the corresponding Hamiltonians. (Auth.)

Solving Non-Isospectral mKdV Equation and Sine-Gordon Equation Hierarchies with Self-Consistent Sources via Inverse Scattering Transform

International Nuclear Information System (INIS)