Abstract:In the present paper we describe a method for solving inverse problems for the Helmholtz equation in radially-symmetric domains given multi-frequency data. Our approach is based on the construction of suitable trace formulas which relate the impedance of the total field at multiple frequencies to derivatives of the potential. Using this trace formula we obtain a system of coupled differential equations which can be solved to obtain the potential in a stable manner. Finally, the performance of the reconstruction algorithm is illustrated with several numerical examples.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to solve the acoustic inverse scattering problem with radial symmetry in the two - dimensional plane given multi - frequency data. Specifically, the goal of the paper is to recover the scattering potential from the measured data of the scattering field. This problem is very important in multiple application fields, such as geophysics, medical imaging, non - destructive testing and sonar, etc.
### Problem Background
The inverse scattering problem is highly nonlinear, which makes its solution very challenging. In some cases, the problem can be simplified by transforming the wave equation into the Schrödinger equation, but this transformation may be numerically unstable and difficult to achieve in two - or three - dimensional cases. Another method is to deal with the nonlinear problem through linearization techniques (such as the Born approximation), but this usually fails in high - contrast environments. Therefore, much research focuses on methods for directly solving the fully nonlinear problem, including the use of nonlinear optimization techniques and signal and image processing techniques.
### Main Contributions of the Paper
1. **New Trace Formula**: The paper proposes a new trace formula that relates the impedance of the field to the scattering potential. This formula is applicable to radially symmetric scattering potentials.
2. **Numerical Solution Method**: Based on the new trace formula, the paper develops a stable and high - order numerical method to solve the inverse scattering problem. This method can effectively recover the scattering potential from multi - frequency data.
3. **Numerical Experiments**: The paper verifies the effectiveness and accuracy of the proposed method through multiple numerical experiments.
### Mathematical Tools
- **Radial Symmetry**: Assume that the scattering potential is radially symmetric, that is, \( Q(r)=q(\|r\|) \), where \( q \) is a continuous function.
- **Wave Equation**: The time - harmonic acoustic wave equation is:
\[
\Delta u(x)+k^{2}(1 + Q(x))u(x)=f(x)
\]
where \( k \) is the frequency, \( f(x) \) is the source term, satisfying the Sommerfeld radiation condition.
- **Fourier Coefficients**: Transform the two - dimensional problem into a series of one - dimensional problems through the Fourier transform.
- **Bessel and Hankel Functions**: These special functions play a key role in the analysis and numerical solution process.
- **Impedance**: Define the impedance \( \phi_{n}(r,k) \) and derive the Riccati equation it satisfies.
### Conclusion
The method proposed in the paper shows high accuracy and efficiency in dealing with two - dimensional acoustic inverse scattering problems with radial symmetry. Future work will be extended to three - dimensional problems and waveguide problems with constant cross - section parameters, as well as non - radially symmetric scattering potential problems.
Hope this summary is helpful for you to understand the content of the paper! If you have any further questions or need more detailed explanations, please feel free to let me know.