Abstract:This paper is concerned with an inverse wavenumber/frequency-dependent source problem for the Helmholtz equation. In two and three dimensions, the unknown source term is supposed to be compactly supported in spatial variables but independent on one spatial variable. The dependence of the source function on wavenumber/frequency is supposed to be unknown. Based on the Dirichlet-Laplacian and Fourier-Transform methods, we develop two effcient non-iterative numerical algorithms to recover the wavenumber-dependent source. Uniqueness proof and increasing stability analysis are carried out in terms of the boundary measurement data of Dirichlet kind. Numerical experiments are conducted to illustrate the effectiveness and efficiency of the proposed methods.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the inverse problem regarding wavenumber - (or frequency -) dependent sources. Specifically, the researchers focus on the inverse problem of wavenumber - dependent source terms in the Helmholtz equation in two - and three - dimensional spaces. The unknown source term has a compact support in the spatial variable, is independent of one spatial variable, and its dependence on the wavenumber (or frequency) is unknown. Based on the Dirichlet - Laplacian method and the Fourier transform method, the researchers developed two efficient non - iterative numerical algorithms to recover the wavenumber - dependent source term, and carried out uniqueness proofs and stability analyses. In addition, the effectiveness and efficiency of the proposed methods were verified through numerical experiments.
### Core Problems of the Paper
1. **Inverse Problem of Wavenumber - Dependent Source Terms**: How to determine the wavenumber - dependent source function \( f(\mathbf{x}, k) \) from boundary measurement data, where \( \mathbf{x} \) is the spatial variable and \( k \) is the wavenumber.
2. **Uniqueness and Stability**: Prove the uniqueness and stability of the wavenumber - dependent source function \( f(\mathbf{x}, k) \) given multi - frequency near - field measurement data.
3. **Numerical Algorithms**: Develop effective numerical algorithms to recover the wavenumber - dependent source terms and verify the effectiveness and efficiency of these algorithms through numerical experiments.
### Main Contributions
- **Uniqueness Proof**: Through the Fourier transform method and the Dirichlet - Laplacian method, the uniqueness of the wavenumber - dependent source function \( f(\mathbf{x}, k) \) under multi - frequency near - field measurement data was proven.
- **Stability Analysis**: Estimates for increasing stability were derived, and in particular, in the two - dimensional case, an explicit upper - bound estimate was given.
- **Numerical Algorithms**: Two efficient non - iterative numerical algorithms for recovering the wavenumber - dependent source terms were proposed, and the effectiveness of these algorithms was verified through numerical experiments.
### Mathematical Background
- **Helmholtz Equation**: \(\Delta u + k^2 u = f(\mathbf{x}, k)\), where \( u \) is the wave field, \( k \) is the wavenumber, and \( f(\mathbf{x}, k) \) is the wavenumber - dependent source term.
- **Sommerfeld Radiation Condition**: \(\lim_{r \to \infty} r^{(d - 1)/2} (\partial_r u - iku) = 0\), where \( r = |\mathbf{x}| \).
- **Fourier Transform**: Used to handle time - dependent source terms and convert them into problems in the frequency domain.
- **Dirichlet - Laplacian Method**: Utilize the completeness and orthogonality of Dirichlet eigenfunctions to construct test functions for uniqueness proofs and numerical solutions.
### Application Background
- **Acoustic Imaging**: In acoustic imaging, the inverse problem of wavenumber - dependent source terms can help determine the position and shape of the sound source.
- **Medical Imaging**: In medical imaging techniques, such as photoacoustic tomography, the inverse problem of wavenumber - dependent source terms can improve the quality and resolution of images.
Through these methods and analyses, the researchers provided a new theoretical framework and numerical tools for solving the inverse problem of wavenumber - dependent source terms.