Iterative Regularization Methods for Inverse Problems in Acoustics

Author(s):  
Thomas K. DeLillo ◽  
Tomasz Hrycak ◽  
Nicolas Valdivia

We consider the use of conjugate-gradient-like iterative methods for the solution of integral equations arising from an inverse problem in acoustics in a bounded three dimensional region. The inverse problem is the computation of the normal velocities on the boundary of a region from pressure measurements on an interior surface. The pressure satisfies the Helmholtz equation in the region. Two formulations are considered: one based on the representation of pressures by a single layer potential and the other based on the Helmholtz-Kirchhoff integral equation. Both formulations can be used to approximate the Neumann Green’s function as an alternative. The integral equations are all ill-posed and are discretized by a boundary element method. The resulting liner systems are ill-conditoned and a (smooth) regularized solutions must be sought. Two regularization rules, including a new one, for conjugate-gradient-like methods are applied and found to have advantages over a standard method based on the truncated singular value decomposition using generalized cross validation. Due to the occurence of multiple singular values for our integral operators, conjugate gradient methods compute the optimal solution in the first few iterations and prove to be particularly fast for these large scale acoustics problems.

2016 ◽  
Vol 26 (3) ◽  
pp. 623-640 ◽  
Author(s):  
Sara Beddiaf ◽  
Laurent Autrique ◽  
Laetitia Perez ◽  
Jean-Claude Jolly

Abstract Inverse three-dimensional heat conduction problems devoted to heating source localization are ill posed. Identification can be performed using an iterative regularization method based on the conjugate gradient algorithm. Such a method is usually implemented off-line, taking into account observations (temperature measurements, for example). However, in a practical context, if the source has to be located as fast as possible (e.g., for diagnosis), the observation horizon has to be reduced. To this end, several configurations are detailed and effects of noisy observations are investigated.


2014 ◽  
Vol 2014 ◽  
pp. 1-14
Author(s):  
San-Yang Liu ◽  
Yuan-Yuan Huang

This paper investigates a general form of guaranteed descent conjugate gradient methods which satisfies the descent conditiongkTdk≤-1-1/4θkgk2  θk>1/4and which is strongly convergent whenever the weak Wolfe line search is fulfilled. Moreover, we present several specific guaranteed descent conjugate gradient methods and give their numerical results for large-scale unconstrained optimization.


1994 ◽  
Vol 261 ◽  
pp. 199-222 ◽  
Author(s):  
C. Pozrikidis

The force and torque on a particle that translates, rotates, or is held stationary in an incident flow within a channel with parallel-sided walls, are considered in the limit of Stokes flow. Assuming that the particle has an axisymmetric shape with axis perpendicular to the channel walls, the problem is formulated in terms of a boundary integral equation that is capable of describing arbitrary three-dimensional Stokes flow in an axisymmetric domain. The method involves: (a) representing the flow in terms of a single-layer potential that is defined over the physical boundaries of the flow as well as other external surfaces, (b) decomposing the polar cylindrical components of the velocity, boundary surface force, and single-layer potential in complex Fourier series, and (c) collecting same-order Fourier coefficients to obtain a system of one-dimensional Fredholm integral equations of the first kind for the coefficients of the surface force over the traces of the natural boundaries of the flow in an azimuthal plane. In the particular case where the polar cylindrical components of the boundary velocity exhibit a first harmonic dependence on the azimuthal angle, we obtain a reduced system of three real integral equations. A numerical method of solution that is based on a standard boundary element-collocation procedure is developed and tested. For channel flow, the effect of domain truncation on the nature of the far flow is investigated with reference to plane Hagen–Poiseuille flow past a cylindrical post. Numerical results are presented for the force and torque exerted on a family of oblate spheroids located above a single plane wall or within a parallel-sided channel. The effect of particle shape on the structure of the flow is illustrated, and some novel features of the motion are discussed. The numerical computations reveal the range of accuracy of previous asymptotic solutions for small or tightly fitting spherical particles.


1995 ◽  
Vol 62 (2) ◽  
pp. 273-281 ◽  
Author(s):  
C. H. Kuo ◽  
L. M. Keer

The three-dimensional problem of a multilayered composite containing an arbitrarily oriented crack is considered in this paper. The crack problem is analyzed by the equivalent body force method, which reduces the problem to a set of singular integral equations. To compute the kernels of the integral equations, the stiffness matrix for the layered medium is formulated in the Hankel transformed domain. The transformed components of the Green’s functions and derivatives are determined by solving the stiffness matrix equations, and the kernels are evaluated by performing the inverse Hankel transform. The crack-opening displacements and the three modes of the stress intensity factor at the crack front are obtained by numerically solving the integral equations. Examples are given for a penny-shaped crack in a bimaterial and a three-material system, and for a semicircular crack in a single layer adhered to an elastic half-space.


Sign in / Sign up

Export Citation Format

Share Document