scholarly journals Variations on the Maiani-Testa approach and the inverse problem

2021 ◽  
Vol 2021 (6) ◽  
Author(s):  
M. Bruno ◽  
M. T. Hansen

Abstract We discuss a method to construct hadronic scattering and decay amplitudes from Euclidean correlators, by combining the approach of a regulated inverse Laplace transform with the work of Maiani and Testa [1]. Revisiting the original result of ref. [1], we observe that the key observation, i.e. that only threshold scattering information can be extracted at large separations, can be understood by interpreting the correlator as a spectral function, ρ(ω), convoluted with the Euclidean kernel, e−ωt, which is sharply peaked at threshold. We therefore consider a modification in which a smooth step function, equal to one above a target energy, is inserted in the spectral decomposition. This can be achieved either through Backus-Gilbert-like methods or more directly using the variational approach. The result is a shifted resolution function, such that the large t limit projects onto scattering or decay amplitudes above threshold. The utility of this method is highlighted through large t expansions of both three- and four-point functions that include leading terms proportional to the real and imaginary parts (separately) of the target observable. This work also presents new results relevant for the un-modified correlator at threshold, including expressions for extracting the Nπ scattering length from four-point functions and a new strategy to organize the large t expansion that exhibits better convergence than the expansion in powers of 1/t.

2012 ◽  
Vol 2012 ◽  
pp. 1-12 ◽  
Author(s):  
Chi-Min Liu

The extended Stokes problems, which study the flow suddenly driven by relatively moving half-planes, are reexamined for the Oldroyd-B fluid. This topic has been studied (Liu, 2011) by applying the series expansion to calculate the inverse Laplace transform. The derived solution was correct but tough to perform the calculation due to the series expansion of infinite terms. Herein another approach, the contour integration, is applied to calculate the inversion. Moreover, the Heaviside unit step function is included into the boundary condition to ensure the consistence between boundary and initial conditions. Mathematical methods used herein can be applied to other fluids for the extended Stokes’ problems.


2021 ◽  
Vol 2092 (1) ◽  
pp. 012001
Author(s):  
Yu Jiang ◽  
Gen Nakamura ◽  
Kenji Shirota

Abstract This paper deals with an inverse problem for recovering the viscoelasticity of a living body from MRE (Magnetic Resonance Elastography) data. Based on a viscoelastic partial differential equation whose solution can approximately simulate MRE data, the inverse problem is transformed to a least square variational problem. This is to search for viscoelastic coefficients of this equation such that the solution to a boundary value problem of this equation fits approximately to MRE data with respect to the least square cost function. By computing the Gateaux derivatives of the cost function, we minimize the cost function by the projected gradient method is proposed for recovering the unknown coefficients. The reconstruction results based on simulated data and real experimental data are presented and discussed.


2011 ◽  
Vol 27 (4) ◽  
pp. 045011 ◽  
Author(s):  
B M Brown ◽  
M Jais

2005 ◽  
Vol 21 (6) ◽  
pp. 1953-1973 ◽  
Author(s):  
B M Brown ◽  
M Jais ◽  
I W Knowles

2016 ◽  
Vol 49 (4) ◽  
pp. 1141-1147 ◽  
Author(s):  
Konstantin Tolstikhin ◽  
Berthold Scholtes

The inverse problem of evaluating residual stresses σ(z) in real space using residual stresses σ(τ) in image space is discussed. This problem is ill posed and special solution methods are required in order to obtain a stable solution. Moreover, the real-space solution must be localized in reflecting layers only in multilayer systems. This requirement imposes strong restrictions on the solution methods and does not allow one to use methods based on the inverse Laplace transform employed for compact solid materials. Besides, in the case of solid materials, the use of the inverse Laplace transform often leads to extremely unstable solutions. The stable numerical solution of the discussed inverse problem can be found using a method based on the Tikhonov regularization. Given the measured data and their pointwise error estimation, this method provides stable approximate solutions for both solid materials and thin films in the form of piecewise functions defined solely in diffracting layers. The approximations are shown to converge to the exact function when the noise in the experimental data approaches zero. If the initial data satisfy certain constraints, the method provides a stable exact solution for the inverse problem. A freely available MATLAB package has been developed, and its efficiency was demonstrated in the numerical residual stress calculations carried out for solid materials and thin films.


Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. E209-E223
Author(s):  
Juan Luis Fernández-Martínez ◽  
Zulima Fernández-Muñiz ◽  
Shan Xu ◽  
Ana Cernea ◽  
Colette Sirieix ◽  
...  

We have evaluated the uncertainty analysis of the 3D electrical tomography inverse problem using model reduction via singular-value decomposition and performed sampling of the nonlinear equivalence region via an explorative member of the particle swarm optimization (PSO) family. The procedure begins with the local inversion of the observed data to find a good resistivity model located in the nonlinear equivalence region. Then, the dimensionality is reduced via the spectral decomposition of the 3D geophysical model. Finally, the exploration of the uncertainty space is performed via an exploratory version of PSO (RR-PSO). This sampling methodology does not prejudge where the initial model comes from as long as this model has a geologic meaning. The 3D subsurface conductivity distribution is arranged as a 2D matrix by ordering the conductivity values contained in a given earth section as a column array and stacking parallel sections as columns of the matrix. There are three basic modes of ordering: mode 1 and mode 2, by using vertical sections in two perpendicular directions, and mode 3, by using horizontal sections. The spectral decomposition is then performed using these three 2D modes. Using this approach, it is possible to sample the uncertainty space of the 3D electrical resistivity inverse problem very efficiently. This methodology is intrinsically parallelizable and could be run for different initial models simultaneously. We found the application to a synthetic data set that is well-known in the literature related to this subject, obtaining a set of surviving geophysical models located in the nonlinear equivalence region that can be used to approximate numerically the posterior distribution of the geophysical model parameters (frequentist approach). Based on these models, it is possible to perform the probabilistic segmentation of the inverse solution found, meanwhile answering geophysical questions with its corresponding uncertainty assessment. This methodology has a general character could be applied to any other 3D nonlinear inverse problems by implementing their corresponding forward model.


Author(s):  
Vincent Braibant ◽  
E. Leménager ◽  
Th. Bouet

Abstract This paper discusses the use of Kalman filters as a new approach for building global approximations. Basic ideas and procedures of Kalman filters are first recalled. Next, key elements of how to implement the method for design problems are described. Finally, in order to evaluate the performance of the approach, an inverse problem which consists in optimizing a warhead with respect to constraints on the resulting projectile is solved. It is shown that global approximations are convenient for the solution of complex optimization problems and that Kalman filtering techniques appear a well suited strategy for the construction of global approximations in structural optimization.


Sign in / Sign up

Export Citation Format

Share Document